粗差检测与抗差估计相结合的方法在动态相对定位中的应用

分析了粗差对GPS动态相对定位的影响,给出了相应的影响函数和抗差估计方法。建议采用粗差探测与抗差估计相结合的方法来处理GPS动态定位中的粗差问题。抗差估计中的等价权分别采用了IGG3方案和双因子等价权函数,对两种方案的效果进行比较和分析,并利用实测数据对所提方法的实际效果进行了检验。结果表明,基于粗差探测和抗差估计相结合的粗差处理方法可有效控制和抵御粗差的影响;分别采用IGG3方案和双因子等价权函数的定位结果差异较小。     

A parametric study on the impact of satellite attitude errors on GOCE gravity field recovery

In the recent design of the Gravity field and steady-state Ocean Circulation Explorer (GOCE) satellite mission, the gravity gradients are defined in the gradiometer reference frame (GRF), which deviates from the actual flight direction (local orbit reference frame, LORF) by up to 3–4°. The main objective of this paper is to investigate the effect of uncertainties in the knowledge of the gradiometer orientation due to attitude reconstitution errors on the gravity field solution. In the framework of several numerical simulations, which are based on a realistic mission configuration, different scenarios are investigated, to provide the accuracy requirements of the orientation information. It turns out that orientation errors have to be seriously considered, because they may represent a significant error component of the gravity field solution. While in a realistic mission scenario (colored gradiometer noise) the gravity field solutions are quite insensitive to small orientation biases, random noise applied to the attitude information can have a considerable impact on the accuracy of the resolved gravity field models.     

一种基于加表零偏稳定性的GNSS/SINS组合导航异常探测方法

在地面车载组合导航中,全球导航卫星系统（global navigation satellite system,GNSS）的观测值容易受地面复杂环境的干扰,导致其定位结果出现异常,严重影响GNSS/捷联惯性导航系统（strap-down inertial navigation system,SINS）组合的滤波解算。从惯导系统误差特性的角度,研究了一种基于加表零偏稳定性的组合导航异常探测新方法。该方法从加表零偏解算的异常来发现GNSS位置、速度等观测值中的粗差,并采取剔除和降权的抗差方法抵御粗差影响。通过一组车载数据的分析表明,观测粗差对加表零偏解算的影响十分显著,以此为判别条件能够准确地发现观测粗差。采用该方法后,位置误差、速度误差和姿态误差的均方根分别减小了70.8%、87.9%和77.7%,显著提高了组合导航的解算精度和鲁棒性,为组合导航数据的抗差处理提供了一种新思路。     

Earth’s gravity field modelling based on satellite accelerations derived from onboard GPS phase measurements

GPS data collected by satellite gravity missions can be used for extracting the long-wavelength part of the Earth’s gravity field. We propose a new data processing method which makes use of the ‘average acceleration’ approach to gravity field modelling. In this method, satellite accelerations are directly derived from GPS carrier phase measurements with an epoch-differenced scheme. As a result, no ambiguity solutions are needed and the systematic errors that do not change much from epoch to epoch are largely eliminated. The GPS data collected by the Gravity Field and Steady-State Ocean Circulation Explorer (GOCE) satellite mission are used to demonstrate the added value of the proposed method. An analysis of the residual accelerations shows that accelerations derived in this way are more precise, with noise being reduced by about 20 and 5% at the cross-track component and the other two components, respectively, as compared to those based on kinematic orbits. The accelerations obtained in this way allow the recovery of the gravity field to a slightly higher maximum degree compared to the solution based on kinematic orbits. Furthermore, the gravity field solution has an overall better performance. Errors in spherical harmonic coefficients are smaller, especially at low degrees. The cumulative geoid height error is reduced by about 15 and 5% up to degree 50 and 150, respectively. An analysis in the spatial domain shows that large errors along the geomagnetic equator, which are caused by a high electron density coupled with large short-term variations, are substantially reduced. Finally, the new method allows for a better observation of mass transport signals. In particular, sufficiently realistic signatures of regional mass anomalies in North America and south-west Africa are obtained.     

抗差估计在多站多星卡尔曼滤波精密定轨中的应用

针对在多站多星卡尔曼滤波定轨中,粗差探测判断标准选择不合适时会影响卡尔曼滤波解精度的问题,采用放宽粗差探测阈值并在滤波阶段采用抗差估计的方法来控制滤波的精度。利用IGS站和BD站的观测数据,分别计算了GPS、BD卫星单天弧度的定轨结果,并与IGS精密星历、武大精密星历作比较。结果表明,利用抗差卡尔曼滤波方法轨道精度得到了提高。     

多种跟踪组合卫星重力场恢复方法初探

将卫星地面跟踪、高低卫星跟踪和低低卫星跟踪恢复重力场的原理进行统一 ,导出了两种能充分综合多种跟踪观测的卫星重力场恢复基本观测方程。在建立基本观测量与重力场参数的函数传输关系的基础上 ,对目前存在的重力场恢复算法给出了直观的解释 ,并以距离二次变率为重力场恢复的基本观测量 ,提出了一种能组合各种卫星跟踪观测量的动态 /动力组合定轨和重力场恢复方法。     

双因子抗差贝叶斯估计及其在GIS数据采集质量控制中的应用

在对GIS扫描数字化过程中,所采集的数据里经常含有系统误差和异常误差。在纠正这些误差的过程中,应同时考虑控制点坐标和数字化坐标的先验协因数。文中以数字化的系统误差拟合模型为例,提出双因子抗差贝叶斯估计误差拟合方法,并分别就控制点含有粗差、数字化点含有粗差和两者均含有粗差3种情况建立相应的模型,实例验证了其抗差效果。     

利用基于抗差垂直向方差分量估计的GPS-InSAR数据融合方法反演三维形变场

GPS-InSAR数据融合解算三维形变场模型易受观测值粗差影响,且基于方差分量估计的定权方法不具备抵御粗差能力,计算效率低下。鉴于此,本文提出了一种基于抗差垂直向方差分量估计的GPS-InSAR融合解算模型,利用方差分量估计方法及抗差估计理论,通过对观测值最优化分类并进行选权迭代,精确分配权重,进而有效计算三维形变场。试验结果表明,该方法能有效抵御观测值粗差不利影响,提高三维形变场反演精度,提升逐点式计算的三维形变场效率。     

一种三步抗差方案的设计

抗差估计的关键是选择恰当的极值函数ρ（．）或ψ（．）函数。本文针对抗差估计涉及的几个重要问题，提出一种新的抗差方案的构想。新的抗差方案试图实现如下三阶段目标；１）选取抗差性能很强的初值，包括定位参数和单位权中误差的初值。２）综合抗差，即要求能实现系数阵空间和观测空间同时抗差，能解决相关观测抗差估计问题，实现定位参数与单位权中误差的联合抗差。３）剔除粗差，给出有关估值的准确统计特性。     

基于卫星加速度恢复地球重力场的去相关滤波法

基于加速度法恢复地球重力场时,卫星加速度是由卫星轨道数值微分得到,而数值微分会放大高频误差,进而降低了重力场解算结果的精度.针对数值微分导出的加速度误差具有有色噪声的特性,提出利用去相关算法构造白化滤波器对加速度有色噪声进行滤波处理,并根据去相关的基本原理分别构造了基于三点差分和ARMA模型的白化滤波器.采用不同噪声背景的CHAMP卫星模拟轨道数据进行解算,结果表明:基于去相关滤波解算的重力场模型精度均要比等权解算的重力场模型精度高,初步验证了去相关滤波方法的有效性.     

Improved GRACE kinematic orbit determination using GPS receiver clock modeling

Kinematic positions of Low Earth Orbiters based on GPS tracking are frequently used as pseudo-observations for single satellite gravity field determination. Unfortunately, the accuracy of the satellite trajectory is partly limited because the receiver synchronization error has to be estimated along with the kinematic coordinates at every observation epoch. We review the requirements for GPS receiver clock modeling in Precise Point Positioning (PPP) and analyze its impact on kinematic orbit determination for the two satellites of the Gravity Recovery and Climate Experiment (GRACE) mission using both simulated and real data. We demonstrate that a piecewise linear parameterization can be used to model the ultra-stable oscillators that drive the GPS receivers on board of the GRACE satellites. Using such a continuous clock model allows position estimation even if the number of usable GPS satellites drops to three and improves the robustness of the solution with respect to outliers. Furthermore, simulations indicate a potential accuracy improvement of the satellite trajectory of at least 40 % in the radial direction and up to 7 % in the along-track and cross-track directions when a 60-s piecewise linear clock model is estimated instead of epoch-wise independent receiver clock offsets. For PPP with real GRACE data, the accuracy evaluation is hampered by the lack of a reference orbit of significantly higher accuracy. However, comparisons with a smooth reduced-dynamic orbit indicate a significant reduction of the high-frequency noise in the radial component of the kinematic orbit.     

天问一号拓展任务对火星低阶重力场解算的潜在贡献分析

天问一号是中国首次独立开展的行星际探测任务,实现了对火星的环绕、着陆、巡视探测。天问一号正常科学任务阶段环绕以极轨设计为主,与历史火星任务类似,对当前火星重力场模型精度改进有限,因此其拓展任务轨道选取至关重要。通过对极轨圆轨道和近赤道大偏心率轨道进行仿真模拟,分析两种典型轨道构型对现有火星重力场模型改进的可能性,基于不同误差考量仿真解算了对应6个重力场模型。借助重力场功率谱分析,发现在测量噪声为0.1 mm/s的情况下,不论采用极轨还是近赤道轨道,一个月的跟踪数据均可较好地反演出42阶次的火星重力场模型;考虑综合误差影响之后,发现两种轨道对于重力场解算精度类似,其中实施近赤道大偏心率轨道对35阶次以上约束略强。     

低轨卫星初始状态误差对反演地球重力场的影响

基于动力学法反演地球重力场的基本理论,研究了卫星初始状态向量误差对应用低轨卫星精密轨道数据反演地球重力场的影响。在仅考虑低轨卫星初始状态误差的情况下进行了模拟计算,结果表明:在利用低轨卫星精密轨道数据反演地球重力场时,卫星初始状态向量误差需要重新进行估计;在目前的轨道精度水平下,若不顾及误差方程二次项的影响,反演弧长不宜过长;卫星初始状态速度误差(约1.5mm/s)的影响要大于位置误差(约10 cm)的影响。     

利用序列平均的高精度GPS基线解算模型

系统分析了改进随机模型和改进函数模型两类GPS基线解算模型的优缺点,在此基础上,提出了一种基于序列平均的高精度GPS基线解算模型,即采用动态单历元技术进行静态基线解算,充分利用多路径效应的低频特性,采用小波变化理论,对坐标序列进行多路径效应的去除,提取低频残差项进行序列平均,得到基线向量解。同时,以动态坐标序列为依据,对出现粗差历元或者卫星进行处理,有效弥补了仅采用残差序列进行粗差判断的不足,提高了基线解算的精度和可靠性。实验表明,新模型可以更为有效地削弱多路径效应的影响,而且对于较短的观测时间尤为突出;结合坐标序列和残差序列,能更为有效地进行粗差的探测和去除,提高基线解的精度和可靠性。     

Outlier detection algorithms and their performance in GOCE gravity field processing

The satellite missions CHAMP, GRACE, and GOCE mark the beginning of a new era in gravity field determination and modeling. They provide unique models of the global stationary gravity field and its variation in time. Due to inevitable measurement errors, sophisticated pre-processing steps have to be applied before further use of the satellite measurements. In the framework of the GOCE mission, this includes outlier detection, absolute calibration and validation of the SGG (satellite gravity gradiometry) measurements, and removal of temporal effects. In general, outliers are defined as observations that appear to be inconsistent with the remainder of the data set. One goal is to evaluate the effect of additive, innovative and bulk outliers on the estimates of the spherical harmonic coefficients. It can be shown that even a small number of undetected outliers (<0.2 of all data points) can have an adverse effect on the coefficient estimates. Consequently, concepts for the identification and removal of outliers have to be developed. Novel outlier detection algorithms are derived and statistical methods are presented that may be used for this purpose. The methods aim at high outlier identification rates as well as small failure rates. A combined algorithm, based on wavelets and a statistical method, shows best performance with an identification rate of about 99%. To further reduce the influence of undetected outliers, an outlier detection algorithm is implemented inside the gravity field solver (the Quick-Look Gravity Field Analysis tool was used). This results in spherical harmonic coefficient estimates that are of similar quality to those obtained without outliers in the input data.     

Exploring gravity field determination from orbit perturbations of the European Gravity Mission GOCE

 A comparison was made between two methods for gravity field recovery from orbit perturbations that can be derived from global positioning system satellite-to-satellite tracking observations of the future European gravity field mission GOCE (Gravity Field and Steady-State Ocean Circulation Explorer). The first method is based on the analytical linear orbit perturbation theory that leads under certain conditions to a block-diagonal normal matrix for the gravity unknowns, significantly reducing the required computation time. The second method makes use of numerical integration to derive the observation equations, leading to a full set of normal equations requiring powerful computer facilities. Simulations were carried out for gravity field recovery experiments up to spherical harmonic degree and order 80 from 10 days of observation. It was found that the first method leads to large approximation errors as soon as the maximum degree surpasses the first resonance orders and great care has to be taken with modeling resonance orbit perturbations, thereby loosing the block-diagonal structure. The second method proved to be successful, provided a proper division of the data period into orbital arcs that are not too long. Received: 28 April 2000 / Accepted: 6 November 2000     

两种渐消滤波与自适应抗差滤波的综合比较分析

给出了两种渐消滤波解及其相应的原则，介绍了自适应抗差滤波原理和相应的解，分别从原理和解的表达式上分析了两种渐消滤波和自适应抗差滤波解算的基本性能。利用一个实际算例比较了渐消滤波与自适应抗差滤波解在控制状态异常影响方面的能力。     

Comparison of Two Fading Filters anc Adaptively Robust Filter

Two kinds of fading filters and their principles are introduced. An adaptive robust filter is given with corresponding principle. The basic abilities of the fading filters and adaptively robust filter in controlling the influences of the kinematic model errors are analyzed. A practical example is given. The results of the fading filter and adaptively robust filter are compared and analyzed.     

两种渐消滤波与自适应抗差滤波的综合比较分析(英文)

Two kinds of fading filters and their principles are introduced. An adaptive robust filter is given with corresponding principle. The basic abilities of the fading filters and adaptively robust filter in controlling the in- fluences of the kinematic model errors are analyzed. A practical example is given. The results of the fading filter and adaptively robust filter are compared and analyzed.     

Comparison of two fading filters and adaptively robust filter

Two kinds of fading filters and their principles are introduced. An adaptive robust filter is given with corresponding principle. The basic abilities of the fading filters and adaptively robust filter in controlling the influences of the kinematic model errors are analyzed. A practical example is given. The results of the fading filter and adaptively robust filter are compared and analyzed.