加速仪数据在CHAMP卫星精密定轨中的贡献

通过用几个实例计算，评估了CHAMP卫星加速仪数据对SLR数据定轨精度的贡献．结果显示用加速仪数据代替非保守力对CHAMP卫星进行精密定轨时，其SLR的残差从原来的16．5cm减少到2．7cm，与精密轨道相比，卫星位置误差由18．9cm减少到6．2cm(一天资料)．     

Preliminary Exploration of Orbit Determination Method by Means of Space-based Angle Measurement Data

The determination of the orbit of a target satellite by means of observed angle measurement data has a very important applicative value. In order to solve this problem, an optimization mathematic model labelled STSM is first constructed. Next is given the necessary numerical procedure of finding the initial values by polynomial fitting of the angle measurement data. Thirdly, the properties of STSM are analyzed in detail, from which the basic method for solving the optimization problem is found. Finally, the above-mentioned method is applied to a large number of simulation calculations and a preliminary analysis of the results is made. The calculated results show that for a given accuracy and density of data, it is feasible to determine the orbital elements of the target satellite by using space-based angle measurement data according to the procedure given.     

The Data Depth-weight-kernel Estimation of the Model Error of Satellite Orbit Determination

It is an objective fact that there exists error in the satellite dynamic model and it will be transferred to satellite orbit determination algorithm, forming a part of the connotative model error. Mixed with the systematic error and random error of the measurements, they form the unitive model error and badly restrict the precision of the orbit determination. We deduce in detail the equations of orbit improvement for a system with dynamic model error, construct the parametric model for the explicit part of the model and nonparametric model for the error that can not be explicitly described. We also construct the partially linear orbit determination model, estimate and fit the model error using a two-stage estimation and a kernel function estimation, and finally make the corresponding compensation in the orbit determination. Beginning from the data depth theory, a data depth weight kernel estimator for model error is proposed for the sake of promoting the steadiness of model error estimation. Simulation experiments of SBSS are performed. The results show clearly that the model error is one of the most important effects that will influence the precision of the orbit determination. The kernel function method can effectively estimate the model error, with the window width as a major restrict parameter. A data depth-weight-kernel estimation, however, can improve largely the robustness of the kernel function and therefore improve the precision of orbit determination.     

A Study on Orbit Determination by Using the Data of Crossovers from Satellite-borne Altimeter

Starting with the status of the developments of oceanic altimetry satellites, the significance of orbit determination by using altimeter data is introduced. Then the error correction model of the altimeter data and the calculational method of the data of crossovers are analyzed. The modification quantities of the errors and the adopted model concerned in the files of the altimeter data from the aircraft JASON-1 are also introduced in detail. Finally, through the calculations of the simulated data and the data of actual measurements, the highest accuracies of the orbit determination reached by solely using the altimeter data and the data of crossovers are analyzed, respectively. And this work provides a valuable reference to the practical applications in future.     

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.     

Analysis of Land Water Storage in Southwest China Based on GRACE Data

The inversion of the variation in the land water storage in Southwest China is carried out by taking advantage of the data obtained by the earth gravity satellite GRACE (Gravity Recovery and Climate Experiment) for 64 months from January 2005 to April 2010. The result shows that by selecting an appropriate Gauss radius (R = 600 km) and taking the average gravitational field of the adopted data as the back-ground gravitational field, the land water storage in Southwest China inverted on the basis of the GRACE data reflects the drought in Southwest China at the beginning of 2010 very well.     

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     

基于智能优化算法的小天体初轨确定

经典的初轨确定方法包括Laplace方法和Gauss方法以及它们的各种变化形式. 除这些经典方法之外, 基于当今光学观测数据的特点, 学者们也陆续提出了一些其他的初轨确定方法, 包括双r (目标距离观测者的距离)方法和可行域方法. 双r方法的一种实现方式是通过猜测某两个时刻(通常是定轨弧段的首、末时刻)目标离观测者的距离, 结合观测者在空间中的位置矢量, 即可求解相应的Lambert弧段作为目标轨道的初始猜测. 进一步, 以其他观测时刻的RMS (Root Mean Square)为优化变量可以改进初始猜测从而确定初轨. 可行域方法则是针对一组初始观测参数(包括赤经、赤纬及其变率), 根据一些初始假设将目标(离观测者的)距离及其变率约束在可行域内, 并通过三角划分逐步逼近的方式寻找到使观测RMS最小的猜测解. 针对一系列模拟观测数据以及实测数据, 将智能优化算法(粒子群算法)应用于这两种初轨方法, 并将结果与改进的Laplace算法的结果进行比较. 由于双r方法不仅可以用于短弧定轨还可用于长弧关联, 所以进一步给出了针对长弧段数据的关联结果.     

基于双频星载GPS数据的LEO卫星运动学定轨研究

运动学定轨是星载GPS特有的定轨方法,该方法不依赖于任何力学模型(地球重力场、大气阻力及太阳辐射压等),尤其适用于受大气阻力影响严重的低轨卫星定轨.基于双频星载GPS数据,研究了运动学定轨原理,讨论了数据预处理方法,建立了一套非差运动学定轨算法.并以GRACE (Gravity Recovery And Climate Experiment)-A、B卫星2008年2月实测数据作为试算验证了本研究方法的有效性和可靠性.GRACE 卫星实测数据计算结果表明:运动学定轨能达到5 cm精度(相对于SLR (Satellite Laser Ranging)),与动力学和简化动力学定轨精度相当.     

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.     

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

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

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%.     

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

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

同步卫星短弧定轨

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

分布逼近的卡尔曼滤波及其在星载GPS卫星定轨中的应用

探讨了利用推广卡尔曼滤波估计非线性系统状态时存在的问题，进而简要介绍了目前自动控制领域广泛使用的分布逼近的卡尔曼滤波UKF(Unscented Kalman Filter)，考虑到卫星定轨中系统的动态方程和量测方程存在严重非线性的情形，提出了将UKF用于星载GPS卫星定轨．实例计算结果表明，UKF的性能要优于推广卡尔曼滤波．     

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.     

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

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

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.     

基于GRACE数据的西南陆地水储量分析

利用地球重力卫星GRACE(Gravity Recovey and Climate Experiment)2005年1月～2010年4月期间64个月的数据,对我国西南地区陆地水储量变化进行了反演.结果表明,选取适当的高斯半径(R=600 km)和所采用数据的平均引力场作为背景引力场,则基于GRACE数据反演的西南陆地水...     

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.