GIM在LEO卫星单频GPS定轨中的应用

电离层延迟误差是单频GPS(Global Positioning System)数据最主要的误差源,为提高基于单频GPS数据的LEO(Low Earth Orbiting)卫星定轨精度,必须消除/减弱GPS观测数据中电离层延迟影响.研究了全球电离层模型GIM(Global IonosphericMaps)在基于单频GPS伪距数据的低轨卫星运动学和动力学定轨中的应用,并通过估算电离层尺度因子的方法消除C/A码伪距观测量中电离层延迟影响.由于LEO卫星星载GPS信号受电离层延迟影响与卫星轨道高度相关,选取了轨道高度在300～800 km的CHAMP(CHAllenging Mini-satellite Payload)、GRACE(Gravity Recovery AndClimate Experiment)、TerraSAR-X及SAC-C等LEO卫星C/A码伪距观测量作为试算数据.CHAMP等卫星实测数据计算结果表明:以JPL(Jet Propulsion Laboratory)发布的GIM模型作为背景模型,通过电离层比例因子法能很好地消除C/A码伪距观测量中电离层延迟影响,提高LEO卫星运动学和动力学定轨精度,其中,CHAMP卫星轨道最低,受电离层延迟影响最严重,定轨精度提高最显著,分别为55.6%和47.6%;SAC-C卫星轨道高度最高,受电离层延迟影响最小,相应的定轨精度提高幅度也最低,分别为47.8%和38.2%.     

GPS星历对LEO星载GPS精密定轨精度的影响

目前,越来越多的低轨卫星上都搭载了用于精密定轨的星载GPs接收机,星载GPS已成为低轨卫星精密定轨的主要手段之一.星载GPS精密定轨精度依赖于GPS星历及钟差精度.基于SHORDE-Ⅲ非差动力学定轨功能,以2005年8月1日至8月7日一周的GRACE卫星实测数据为例,采用事后精密轨道(igs)、快速轨道(igr)和超快速轨道(igu)三种GPS星历在同等条件下定轨,估计GPS星历精度对低轨卫星定轨精度的影响,实际计算结果表明igs和igr两类GPS星历定轨精度相当,约为9.5 cm,igu星历定轨精度略低于igs和igr星历,约为10.5cm:高频GPS卫星钟差数据对定轨精度会产生1-6cm影响.     

多星混合资料的定轨

在空间目标光学观测资料定轨中常常会遇到多个目标的观测被标记为同一个目标的情况,由于包含了多个目标的数据,定轨过程无法收敛或者完全错误.从极大似然估计角度,采用EM(Expectation Maximum)方法提出一种将轨道改进和识别过程相互融合的处理方法,并在具体实现过程中给出一种稳健估计方法.数值模拟表明方法简便、有效、可行.采用站间差历元差的手段,对当前电离层变化值进行求解与预测,可以将电离层延时误差的变化控制在一定范围,满足频率传递的要求.     

小倾角LEO多星天基平台光学跟踪GEO精密定轨

针对地基卫星测控系统(Tracking Telemetry and Command, TT&C)系统对地球静止轨道(Geostation-\lk ary Earth Orbit, GEO)卫星在空间和时间覆盖上的局限性, 提出小倾角低地球轨道(Low Earth Orbit, LEO)多星组网天基平台对GEO卫星进行跟踪定轨的方法. 根据空间环境和光学可视条件对仿真数据进行筛选以模拟真实的观测场景, 利用光学测角数据, 使用数值方法对GEO卫星的轨道进行确定. 结果与参考轨道进行重叠对比, 在平台轨道精度5 m、测量精度5rq\rq、 定轨弧长12 h的情况下, 两颗LEO卫星对GEO卫星进行跟踪定轨的精度可达到千米量级, 4颗LEO卫星对GEO目标进行跟踪定轨的精度可达到百米量级. 随着LEO组网卫星数量的增加, 定轨精度得到了较大的提高.     

星载GPS低轨卫星运动学定轨及研究进展

重点论述了星载GPS低轨卫星运动学定轨的基本原理和方法,指出了各种运动学定轨法的优点和不足,分析了各种运动学定轨法解算的特点,并给出了相应的处理策略,同时着重介绍了星载GPS低轨卫星运动学定轨及相关研究的进展。     

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

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

关于北斗G2卫星的溯往追终

2022年1月,失效的北斗G2卫星被实践21号卫星从地球静止轨道拖入了坟墓轨道.为了这项捕获任务的安全实施,需要预先确定北斗G2的旋转状态.基于过去10 yr的测光观测数据展示了北斗G2卫星自转的演化过程.根据北斗G2的自转速度和轨道的演化,确认了在过去的10 yr里发生的6次异常事件.据推测, 2012年的小碎片碰撞事件,是随后几年燃油泄漏的导火索. 2017年之后剩余燃油完全释放,再也没有出现转速异常.将2014年太阳能帆板损坏和2016年的解体事件后建立的旋转动力学模型外推1 yr,转轴的标准偏差小于3°,转速标准偏差为0.11°·s^-1,能够有效地满足捕获任务时刻旋转状态的精度要求.     

关于北斗G2卫星的溯往追终

2022年1月, 失效的北斗G2卫星被实践21号卫星从地球静止轨道拖入了坟墓轨道. 为了这项捕获任务的安全实施, 需要预先确定北斗G2的旋转状态. 基于过去10 yr的测光观测数据展示了北斗G2卫星自转的演化过程. 根据北斗G2的自转速度和轨道的演化, 确认了在过去的10 yr里发生的6次异常事件. 据推测, 2012年的小碎片碰撞事件, 是随后几年燃油泄漏的导火索. 2017年之后剩余燃油完全释放, 再也没有出现转速异常. 将2014年太阳能帆板损坏和2016年的解体事件后建立的旋转动力学模型外推1 yr, 转轴的标准偏差小于3°, 转速标准偏差为0.11°· s^-1, 能够有效地满足捕获任务时刻旋转状态的精度要求.     

星载测高仪交叉点数据定轨研究

从海洋测高卫星的发展状况入手,介绍了用测高仪数据进行定轨的意义.然后,分析了测高仪数据的误差修正模型及交叉点数据计算方法,详细介绍了JASON-1测高数据文件中所涉及的误差修正量和采用的模型.最后,通过仿真数据和实测数据的计算,分析了单用测高仪数据和交叉点数据进行轨道确定所能达到的精度,为未来的实际应用提供了有价值的参考.     

一种初轨计算的稳健方法

光学测角资料的初轨计算在空间目标搜索发现中具有重要作用,当观测资料存在野值时,基于最小二乘的经典初轨计算方法不够稳健.采用最小一乘方法建立了一种初轨计算的稳健方法,方法将初轨计算问题转换为线性规划问题求解,并通过bootstrap方法给出估计精度.数值计算结果表明方法稳健有效,并具有较高的崩溃点.     

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     

The Application of GIM in Precise Orbit Determination for LEO Satellites with Single-Frequency GPS Measurements

With the precise GPS ephemeris and clock error available, the iono- spheric delay is left as the dominant error source in the single-frequency GPS data. Thus, the removal of ionospheric effects is a ma jor prerequisite for an improved orbit reconstruction of LEO satellites based on the single-frequency GPS data. In this paper, the use of Global Ionospheric Maps (GIM) in kine- matic and dynamic orbit determinations for LEO satellites with single-frequency GPS pseudorange measurements is discussed first, and then, estimating the iono- spheric scale factor to remove the ionospheric effects from the C/A-code pseu- dorange measurements for both kinematic and dynamic orbit determinations is addressed. As it is known that the ionospheric delay of space-borne GPS sig- nals is strongly dependent on the orbit altitudes of LEO satellites, we select the real C/A-code pseudorange measurement data of the CHAMP, GRACE, TerraSAR-X and SAC-C satellites with altitudes between 300 km and 800 km as sample data in this paper. It is demonstrated that the approach to eliminating ionospheric effects in C/A-code pseudorange measurements by estimating the ionospheric scale factor is highly effective. Employing this approach, the accu- racy of both kinematic and dynamic orbits can be improved notably. Among those five LEO satellites, CHAMP with the lowest orbit altitude has the most remarkable improvements in orbit accuracy, which are 55.6% and 47.6% for kine- matic and dynamic orbits, respectively. SAC-C with the highest orbit altitude has the least improvements in orbit accuracy accordingly, which are 47.8% and 38.2%, respectively.     

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

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

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

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

TLE-Aided Orbit Determination Using Single-Station SLR Data

It is difficult to use the single-station satellite laser ranging (SLR) data for orbit determination, due to the singular geometrical distribution of the observations. The single-station data produced by performing the diffuse- reflection SLR on the earth-orbiting space debris are therefore ineffective for orbit improvement. To solve this problem, we propose an orbit determination method by using single-station SLR data in aid of the two-line element set (TLE). For verifying its feasibility, this method is implemented and applied to the orbit determination of the satellite Ajisai, using the single-station SLR data of five passes in one day and the corresponding TLE. And on this basis, the five-day orbit prediction is generated, the result indicates that the errors of predicted positions are less than 40 m. In addition, the potential application of this method in the orbit improvement of space debris is discussed.     

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.     

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.