Research on the Method of Orbit Determination Based on the Self-adaptive Stable Least p-norm Estimate

The traditional least square estimation (LSE) method for orbit determination will not be optimal if the error of observational data does not obey the Gaussian distribution. In order to solve this problem, the least p-norm (Lp) estimation method is presented in this paper to deal with the non-Gaussian distribution cases. We show that a suitable selection of parameter p may guarantee a reasonable orbit determination result. The character of Lp estimation is analyzed. It is shown that the traditional Lp estimation method is not a robust method. And a stable Lp estimating based on data depth weighting is put forward to deal with the model error and outlier. In the orbit determination process, the outlier of observational data and coarse model error can be quantitatively described by their weights. The farther is the data from the data center, the smaller is the value of data depth and the smaller is the weighted value accordingly. The result of the new Lp method is stabler than that of the traditional Lp estimation and the breakdown point could be up to 1/2. In addition, the orbit parameter is adaptively estimated by residual analysis and matrix estimation method, and the estimation efficiency is enhanced. Finally, by taking the Space-based Space Surveillance System as an example and performing simulation experiments, we show that if there are system error or abnormal value in the observational data or system error in satellite dynamical model and space-based observation platform, LSE will not be optimal even though the observational data obeys the Gaussian distribution, and the orbit determination precision by the self-adaptive robust Lp estimation method is much better than that by the traditional LSE method.     

Homogeneous and Heterogeneous Observational Data Weighting Method of Combined Orbit Determination of Near-earth Satellites of Bi-satellite Positioning System and its Applications

In the light of the problem of amalgamation and processing of multisource observational data in the combined orbit determination of near-earth satellites of the bi-satellite positioning system, the optimal weighting method of the improved variance component estimation of the two-step systematic error correction of homogeneous observational data is proposed. Analyses show that the multi-source amalgamation measurement model of the heterogeneous observational data essentially is a multi-structure, multi-parameter non-linear regression model, and the optimal weighting method of the combination of model structure characteristic analysis and variance component estimation of the heterogeneous observational data is established. The realization algorithms of the optimal weighting and the combined orbit determination parameter estimation of the two sorts of observational data are designed, and the simulation experiments of the combined orbit determination are carried out by taking the distances among the two satellites and the backup satellite and the homogeneous observational data and the distance between the two satellites and the heterogeneous observational data of satellite sensor angle measurements as the examples. The results of theoretical analysis and simulation calculation show that for the combined orbit determination of homogeneous observational data, the accuracy of orbit determination obtained by adopting the variance component estimation method of the two-step systematic error correction can be more superior than that obtained by means of the traditional empirical weighting method. For the combined orbit determination of heterogeneous observational data, through the introduction of the weighting factor by which the model structure is characterized the accuracies of the combined orbit determination of the near-earth satellite and geostationary satellite are both improved to a certain extent in comparison with the mean weighting mode.     

星载GPS自主定轨的长期稳定性

通过模拟长达100天的星载GPS伪距观测资料,进行卡尔曼滤波定轨仿真计算．重点研究：1)采用简化的动力学模型与简化状态转移矩阵,是否保证滤波的长期稳定性；2)模型误差矩阵Q的选取对滤波定轨精度的影响；3)与事后最小二乘批处理相比较,在简化模型下自主定轨的精度．同时给出了相应的结论．     

同步卫星短弧定轨

同步卫星受到摄动力的影响,它的实际轨道有一点漂移．卫星需要不断的调轨调姿,以保证其正常运行．为了研究卫星在几小时,甚至更短的时间内的轨迹情况,采用短弧段定轨法．用动力学方法进行短弧定轨,分别研究1小时和15分钟定轨并进行比较,目的是为了在同步轨道卫星变轨后,能尽快地为卫星提供精密的预报轨道．此外,在系列短弧定轨后,得到精密轨道系列,为研究轨道变化的力学因素及研究短弧中卫星转发器时延变化规律等提供依据．     

The Effect of GPS Ephemeris on the Accuracy of Precise Orbit Determination of LEO Satellite-borne GPS

The satellite-borne GPS receivers dedicated to precise orbit determination are now being carried by more and more low earth orbit (LEO) satellites and the satellite-borne GPS has become one of the main means for the precise orbit determination of low earth orbit satellites. The accuracy of satellite-borne GPS precise orbit determination depends on the accuracies of the GPS ephemeris and the clock error. Based on the orbit determination function of SHORDEIII zero-difference dynamics and using the observational data obtained by the GRACE satellites for the week from 2005 August 1 to 7 as an example, three versions of GPS ephemerides (igs, igr and igu) are used to carry out orbit determination under the same conditions and to estimate the effect of the GPS ephemeris accuracy on the accuracy of orbit determination of low earth orbit satellites. Our calculated results show that the two ephemerides, igs and igr, are equivalent to each other in orbit determination accuracy (about 9.5 cm), while igu is slightly less accurate, at about 10.5 cm. The effect produced by the data of the high frequency GPS satellite clock error on the accuracy of orbit determination is 1–6 cm.     

LP估计在星载GPS运动学定轨中的应用及精度分析

为得到高精度的星载GPS运动学定轨,必须利用观测精度高的相位观测值,但是相位观测值预处理后,仍然存在残余小周跳．在残差服从正态分布情况下LS法是最佳参数解算方法,但该方法不能解决资料的系统误差消除问题,LP估计是处理资料残差分布含有系统误差的有效方法之一．基于LS、LP方法的有效条件和GPS数据预处理的特性,将LP估计方法引入星载GPS运动学定轨数据处理中,以CHAMP卫星资料为例,研究了LP估计在星载GPS运动学定轨中的应用及其精度分析．实践表明:在处理含有残余小周跳的相位观测值时,LP估计比LS更有效,提高了星载GPS运动学定轨精度,但随着残余周跳的进一步修复,LP估计相对于LS估计的优越性越来越弱,在资料完全没有系统误差,残差服从正态分布的情况下,LP估计不能很好地体现其优越性,精度反而低于LS估计．     

基于GPS与三轴磁强计的联合导航算法

为了修正近地轨道(小于1000 km)地磁场模型,提高导航的精度,在地磁导航系统中引入GPS作为一个新的测量设备,提出了一种基于三轴磁强计与GPS的联合导航算法．该算法取卫星的位置和速度向量作为状态向量,建立状态方程;取卫星周围的磁场强度和GPS接收的信号作为观测量,建立观测方程;并以GPS确定的轨道状态量为标准量,去估计磁场模型参数的修正量,构成冷备冗余导航算法．仿真结果表明,提出的导航算法对轨道位置的估计误差小于50 m,速度的估计误差小于0．1 m／s,导航算法的精度和收敛性都优于使用单一地磁导航的系统．     

星载GPS精密测轨研究及应用

星载ＧＰＳ定轨系统由于其全天侯、价格低、不受卫星高度的影响可达到米级、分米级乃至厘米级的测轨精度等特点已成为低轨道卫星精密定轨的要求,概述了星载ＧＰＳ系统的组成和定轨原理,给出了星载ＧＰＳ测轨的几种主要方法和数学模型,同时根据ＴＯＰＥＸ卫星星载ＧＰＳ实测数据分析了各种定轨方法的测轨精度以及影响定轨精度的各种因素。     

Kinematic Precise Orbit Determination for LEO Satellites Using Space-borne Dual-frequency GPS Measurements

As a special approach to orbit determination for satellites with spaceborne GPS receivers, the kinematic Precise Orbit Determination (POD) is independent of any mechanical model (e.g., the Earth gravity ?eld, atmospheric drag, solar radiation pressure, etc.), and thus especially suitable for the orbit determination of Low Earth Orbiting (LEO)satellites perturbed strongly bythe atmosphere. In this paper, based on the space-borne dual-frequency GPS data, we study the kinematic POD, discuss the pre-processing of the data, and construct an algorithm of zero-difference kinematic POD. Using the observational data from GRACE (Gravity Recovery And Climate Experiment) satellites covering the whole month of February 2008, we verify the effectiveness and reliability of this algorithm. The results show that the kinematic POD may attain an accuracy of about 5 cm (with respect to satellite laser ranging data), which is at the same level as the dynamic and reduced-dynamic PODs     

Analysis of the wide area differential correction for BeiDou global satellite navigation system

The regional BeiDou Satellite System, or BDS2, broadcasts a differential correction as Equivalent Satellite Clock Correction to correct both orbit and satellite clock errors. For the global BDS, or BDS3, satellite orbit and clock corrections conforming with RTCA standards will be broadcast to authorized users. The hybrid constellation and regional monitoring network pose challenges for the high precision separation of orbit and satellite clock corrections. Three correction models of kinematic,dynamic and Two-way Satellite Time Frequency Transfer(TWSTFT)-based dynamic were studied to estimate the satellite orbit and clock corrections. The correction accuracy of the three models is compared and analyzed based on the BDS observation data. Results show that the accuracies(root mean square, RMS) of dual-frequency real-time positioning for the three models are about 1.76 m, 1.78 m and 2.08 m respectively, which are comparable with the performance of WAAS and EGNOS. With dynamic corrections, the precision of Precise Point Positioning(PPP) experiments may reach about 23 cm after convergence.     

两行根数辅助的SLR单站定轨

单站测距资料定轨的困难限制了漫反射SLR(Satellite Laser Ranging)测距资料的应用.为此,提出利用两行根数模拟多站SLR测距资料作为辅助,实现单站SLR测距资料定轨的方法.该方法对卫星Ajisai单站SLR测距资料定轨并生成5 d预报轨道,误差小于40 m,实现利用单站测距资料的轨道改进,验证了方法的可行性.     

月球探测器过渡轨道的短弧定轨方法

论述的短弧定轨,是指在无先验信息情况下又避开多变元迭代的初轨计算方法,它需要相应的动力学问题有一能反映短弧内达到一定精度的近似分析解．探测器进入月球引力作用范围后接近月球时可以处理成相对月球的受摄二体问题,而在地球附近,则可处理成相对地球的受摄二体问题,但在整个过渡段的力模型只能处理成一个受摄的限制性三体问题．而限制性三体问题无分析解,即使在月球引力作用范围外,对于大推力脉冲式的过渡方式,相对地球的变化椭圆轨道的偏心率很大(超过Laplace极限),在考虑月球引力摄动时亦无法构造摄动分析解．就此问题,考虑在地球非球形引力(只包含J2项)和月球引力共同作用下,构造了探测器飞抵月球过渡轨道段的时间幂级数解,在此基础上给出一种受摄二体问题意义下的初轨计算方法,经数值验证,定轨方法有效,可供地面测控系统参考．     

基于SLR精密轨道的天文定位精度分析

利用全球卫星激光测距服务系统(ILRS,International Laser Ranging Service)标准点资料对Ajisai卫星进行精密定轨,残差均方根(RMS)优于3 cm,得到该星的精密轨道.进而对长春站40 cm空间碎片光电望远镜获得的Ajisai卫星的天文定位资料进行精度分析,外符合精度约3″左右.单独利用天文定位数据进行轨道改进,内符合精度优于3″.改进轨道的x、y、z坐标3分量在观测数据覆盖范围内的精度在100 m之内.同样地对Jason-1卫星作数据分析,结果和Ajisai卫星精度相当.分析各个弧段的精度变化,发现定标星个数减少,会导致天文定位精度下降.据此提出可以把最少定标星比例作为评定数据质量的参考指标之一.     

基于GEO卫星分段测距资料的定轨研究

为了满足工程需求,充分利用观测资料信息,快速提供较好的卫星位置预报精度是必要的。在精密定轨的基础上,结合某GEO(地球同步轨道)卫星连续3天的测距观测资料,研究了使用GEO卫星的不同弧长的观测资料进行定轨和轨道预报的精度。结果表明,每天间隔11小时选取该段观测开始30分钟后的5分钟资料进行定轨,可以得到一个为期2天、精度好于15 m的预报轨道。     

Estimation of Distribution Algorithm for Initial Orbit Determination of Too-short-arc Based on Kernel Density Estimation

A new approach of the initial orbit determination for too-short-arc with angular measurements is implemented by building the probabilistic model in the solution space with the estimation of distribution algorithm. Without any assumption about distribution, the non-parametric kernel density estimation is employed in the model building. The method, unlike other evolutionary algorithms, such as the genetic algorithm and particle swarm optimization, considers the fitness as well as the characteristics of the solution space. Numerical experiments with real observations indicate that without any constraints, the proposed technique has a good performance for the observations of usual accuracy.     

北斗三号卫星精密定轨中的光压模型研究

对于在轨运行的BDS (BeiDou Navigation Satellite System)卫星, 太阳光压是作用在卫星上主要的非引力摄动. 受多种因素的影响, 太阳光压摄动力难以精确建模, 是BDS卫星精密定轨和轨道预报过程中重要的误差来源. 由于ECOMC (Empirical CODE Orbit Model 1 and 2 Combined)模型兼顾了ECOM1 (Empirical CODE Orbit Model 1)和ECOM2 (Empirical CODE Orbit Model 2)模型的特点, 在模型中引入了较多的待估参数, 使得参数之间存在强相关性. 针对ECOMC模型的这一缺陷, 文中收集了2019年1月至2022年4月武汉大学分析中心提供的BDS-3卫星精密星历, 采用动力学轨道拟合方法得到了ECOMC模型的13个光压参数. 通过对该模型的光压参数进行时间序列分析, 分别给出了BDS-3 IGSO (Inclined Geosynchronous Orbit)和MEO (Medium Earth Orbit)卫星光压模型的参数选择策略. 并利用轨道拟合和轨道预报试验, 验证了光压模型参数选择策略的合理性. 结果表明, 采用改进型ECOMC模型进行BDS-3 IGSO和MEO卫星轨道拟合的效果最佳, 同时, 也能够提升BDS-3 IGSO和MEO卫星中长期轨道预报的精度.     

Orbit determination of satellite “Tance 1” with VLBI data

The satellite “Tance 1” of the “Double-Star Program” is the first truly scientific experimentation satellite of China. Its orbit is the farthest so far launched in China, with a geocentric apogee reaching 78 thousand kilometers. The tracking of “Tance 1” and of more distant space targets, such as the lunar exploration craft, can be realized with the VLBI technique of radio astronomy. In order to test and verify the role which the VLBI technique plays in the lunar exploration program of China, Shanghai Astronomical Observatory organized the only 3 tracking stations in China (located at Shanghai, Urumqi and Kunming), to carry out test tracking of “Tance 1,” and used the time delay data obtained to determine the orbit of “Tance 1” over a two-day period, so providing a preliminary assessment of the possibility of VLBI orbit determination. The fitting error of the orbit so obtained is about 5.5 m in the time delay and about 2 cm/s in the delay rate (this for checking only), much better than is provided by the preliminary orbit (used merely for ensuring tracking) in which the corresponding figures are around 2 km and 15 cm/s. Further, if the orbit is determined by using both the time delay and time delay rate data (with weights according to their internal accuracies), then the residuals are 5.5 m in the time delay and 2 cm/s in the delay rate. For an appreciation of the true accuracy of the VLBI orbit determination, we used simulation data (of the observed two-day VLBI data) and found the results depended greatly on the error in the dynamic model of the satellite which, however, is difficult to assess, while the formal residuals are of the order of 1 kin in the delay and of cm/s in the delay rate. The simulation computation also indicates that a joint determination using both VLBI and USB data will have an improved accuracy.     

方差分量估计方法在ERS—2卫星精密定轨中的应用

将方差分量估计（VCE）方法应用于ERS－2卫星的精密定轨,用SLR和PRARE资料计算了1998年前3个月的23个长度为5天的弧段（除了调整轨道的时段外,相邻弧段有两天的重叠）,从观测值残差分析和重叠弧段比较两个方面,考察VCE方法对定轨计算的影响,并给出了各组观测值的平均验后均方差,对观测值残差的分析表明,使用VCE方法明显地改善了观测值的拟合程度,但从阿卑（Abbey)标准对观测值残差的检验结果来看,使用VCE方法不能消除轨道中由力学模型和几何模型误差引入的系统差,重叠弧段比较的结果表明：（1）使用VCE方法缩小了重叠弧段的平均距离差,并改善了一部分权段明显不合理的偏离,使最后得到的轨道具有更均匀的精度,（2）相比较而言,VCE方法使相邻弧段靠拢的趋势在轨道切向体现得较为明显,由各组观测值的平均验后方差可见,说单个标准点观测值而言,部分SLR台站的观测资料在定轨计算中占有比其他观测资料更重的地位,纵观全文,使用VCE后得到的观测值的平衡验后均方差来确定资料的双重将比使用均方差更为合理。     

环月飞行器精密定轨的模拟仿真

以中国正在实施的探月计划“嫦娥1号”工程为背景，分析了在中国联合S波段(USB)测控网和甚长基线射电干涉(VLBI)跟踪网的现有空间分布、观测精度水平下的环月飞行器精密定轨．采用的方法是模拟仿真计算，即首先模拟观测数据，然后在计入各误差源的影响后进行求解，并对解算结果进行比较．模拟仿真的工具是美国宇航局哥达德飞行中心的空间数据分析软件系统GEODYN．环月飞行的主要误差源是月球重力场，为此首先讨论了目前精度最高的月球重力场模型JGL165P1的(形式)误差．在模拟了测距、测速以及VLBI的时延、时延率数据后，计入月球重力场的误差进行精密轨道确定．定轨时采用了减缩动力学(reduced dynamic)方法，即选用合适的经验加速度参数吸收重力场误差对定轨的影响．结果表明对于一个不将月球重力场作为主要科学目标的探月计划(如“嫦娥1号”)，减缩动力学方法是一个简单、有效地提高环月飞行器定轨精度的方法．     

Application of variance component estimation to precise orbit determination for ERS-2

Variance component estimation (VCE) is applied to precise orbit determination (POD) of the ERS-2 satellite. Twenty 5-day long arcs in the early three months in 1998 were calculated using the SLR and PRARE data. In the data the adjacent arcs overlap for two days except the intervals for orbit maneuver. The effect of VCE orbit determination on the calculation is investigated by an analysis of residuals and comparison of overlapping arcs, and the mean a posteriori standard deviation of each group of measured residuals is given. It is shown by the residuals analysis that the fitting of the measurements is significantly improved by VCD. However, according to Abbey criterion, VCD is not able to eliminate the systematic errors due to errors in the dynamic and geometric models. The results of the comparison of the overlapping arcs show that (1) VCE reduces the mean range deviation of overlapping arcs, especially where there are obviously unreasonable deviations, so that the orbit obtained has a more uniform precision; (2) By using VCE, adjacent arcs tend to close up and this is more apparent in the transverse direction. From the mean a posteriori standard error of each group of measurements, it can be seen that as far as the single normal point measurement is concerned, the data of some SLR stations are more important than other measurements in POD calculation. Generally speaking, determination of weighting by using VCE is more reasonable than by using initial standard deviation.