导航卫星的地面站定轨算法

导航星座轨道的长期保持是星座导航系统运营管理的重要组成部分,而现有的导航卫星地面定轨算法又存在精度不高或计算量大不适合工程应用的问题。为此,研究了单向、被动测量模式的导航卫星地面定轨算法。基于单向伪距观测,将导航卫星钟差参数作为状态量,推导了滤波算法的状态方程、测量方程,并最终建立了滤波器模型。以不同轨道面的4颗GPS导航卫星为例进行了2天的仿真试验,考虑卫星的可见性仿真中加入了测量中断,并设计在测量恢复后重启滤波算法。仿真结果表明,4颗卫星的轨道位置估计精度可以达到米级,钟差随机偏差的估计精度可以达到纳秒级,并且在滤波中断后重启滤波器,仍然可以达到此估计精度,表明此定轨算法具有收敛性和稳定性。     

基于Gooding算法的天基光学目标跟踪定轨

开展了基于Gooding算法的400km天基平台光学目标监测的轨道确定研究,当测量误差为3”和6”时,分别对800,1500及36000km轨道高度目标进行初始轨道确定及轨道改进分析.仿真结果表明,利用400km轨道高度平台对800~1500km轨道高度目标进行初定轨,测量数据误差为3”~6”时,4~15min弧段的初定轨精度约在10km量级,1~2min弧段的初定轨精度约在100km量级;15min初定轨弧段轨道改进后误差在100m量级,弧段小于10min时轨道改进误差精度在km量级.利用400km轨道高度平台对36000km轨道高度目标进行初定轨,测量数据误差为3”时,15~20min弧段的初定轨精度约在数十km量级,8~10min弧段的初定轨精度在100km量级;轨道改进后误差在km量级.测量数据误差为6”时,20min弧段初定轨精度在10km量级,8~15min弧段初定轨精度在100km量级,轨道改进后误差精度在10km量级.      

Comparison analyses on the 150 × 150 lunar gravity field models by gravity/topography admittance,correlation and precision orbit determination

We analyzed the 150 × 150 lunar gravity field models, LP150Q, GLGM-3 and SGM150, using the power spectrum on the lunar nearside and farside, the lunar global and localized gravity/topography admittance and correlation, and Chang’E-2 precision orbit determination to investigate which model is a more effective tool to estimate geophysical parameters and determine the lunar satellite precision orbit. Results indicate that all gravity field models can be used to estimate the lunar geophysical parameters of the nearside of the Moon. However, SGM150 is better in such computation of the farside. Additionally, SGM150 is shown to be the most useful model for determining the lunar satellite orbit.     

编队星座相对位置精确测定与自主定轨方法

针对编队飞行星座,提出了一种全自主的高精度定轨与相对位置精确测定方法,对其主要关键技术之一的空间绝对定向的基本方案及其可行性和指标进行了研究,并进行了仿真验证.结果表明,通过编队星座星间高精度的距离测量和绝对定向观测,可以实现无需地面测控站和卫星导航系统支持下的编队星座全自主导航方案;在一定测距和测向误差条件下,绝对定轨精度可以优于20m,相对位置确定精度可优于10cm.相对误差与采样间隔有较明显关系,未来可考虑采样间隔控制在10s以内即可;绝对位置误差大小与采样间隔无明显关系,其中的主要误差是星座的整体平移误差.仿真结果验证了所提方案的正确性.      

GPS-only gravity field recovery with GOCE,CHAMP, and GRACE

Gravity missions such as the Gravity field and steady-state Ocean Circulation Explorer (GOCE) are equipped with onboard Global Positioning System (GPS) receivers for precise orbit determination (POD), instrument time-tagging, and the extraction of the long wavelength part of the Earth’s gravity field. The very low orbital altitude of the GOCE satellite and the availability of dense 1 s GPS tracking data are ideal characteristics to exploit the contribution of GPS high-low Satellite-to-Satellite Tracking (hl-SST) to gravity field determination. We present gravity field solutions based on about 8 months of GOCE GPS hl-SST data from 2009 and compare the results with those obtained from the CHAllenging Minisatellite Payload (CHAMP) and Gravity Recovery And Climate Experiment (GRACE) missions. The very low orbital altitude of GOCE significantly improves gravity field recovery from GPS hl-SST data above degree 20, but not for the degrees below 20, where the quality of the spherical harmonic coefficients remains essentially unchanged. Despite the limited time span of GOCE data used, the gravity field of the Earth can be resolved up to about degree 115 using GPS data only. Empirically determined phase center variations (PCVs) of the GOCE onboard GPS helix antenna are, however, mandatory to achieve this performance.     

基于光学测角数据的风云四号同步轨道卫星精密定轨

针对风云四号同步卫星的精密定轨和精度评估需求,首先利用地面光学测角数据对FY-4A卫星进行精密定轨,定轨后方位角和高度角的残差rms分别为0.25＂和0.45＂。与基于测距数据的轨道相比,位置精度在有测角数据的弧段内小于50m。进一步联合测角数据和测距数据对FY-4A卫星进行联合定轨,定轨后轨道重叠精度优于15m。利用联合定轨结果评估了基于测距数据的实时轨道产品精度,可以明显发现轨道精度随着测距数据的积累而逐步提高。     

利用星地差分GPS的地基测控系统实时标校方法

针对地基测控系统传统标校方法和基于差分GPS事后标校方法的不足,提出了一种基于低轨卫星与地面测控站之间星地差分GPS的地基测控系统测量误差实时标校方法。与基于差分GPS的事后标校相比,实时标校能使地基测控系统及时获取标校后的测量数据,从而实时进行轨道解算和预报,并及时上注以提升卫星运行性能。针对星地长基线、高动态和实时标校场景,系统地分析了影响星地基线估计性能的各项误差及修正效果,并提出相对位置精度因子的概念,由此得到星地基线估计精度预算。采用基于抗差自适应卡尔曼滤波的实时星地基线估计算法,并利用加权最小二乘法求解测控系统测量误差,从而获得校准结果。利用星载双频GPS接收机和导航信号模拟器构建半实物仿真平台,仿真结果表明,实时标校后测距系统误差残差降低到40cm左右,测速系统误差残差降低到1cm/s以下,与理论分析结果一致,可以较好地满足未来航天任务的测控需求。     

Formation control using μ-synthesis for Inner-Formation Gravity Measurement Satellite System

Inner-Formation Gravity Measurement Satellite System (IFGMSS) is used to map the gravity field of Earth. The IFGMSS consists of two satellites in which one is called “inner satellite” and the other one is named as “outer satellite”. To measure the pure Earth gravity, the inner satellite is located in the cavity of the outer satellite. Because of the shield effect of the cavity, the inner satellite is affected only by the gravitational force, so it can sense Earth gravity precisely. To avoid the collision between the inner satellite and the outer satellite, it is best to perform a real-time control on the outer satellite. In orbit, the mass of the outer satellite decreases with the consumption of its propellant. The orbit angular rate of the inner satellite varies with time due to various disturbing forces. These two parameters’ uncertainties make the C–W function be not so accurate to describe the formation behavior of these two satellites. Furthermore, the thrusters also have some uncertainties due to the unmodelled dynamics. To cancel the effects caused by the above uncertainties, we have studied the robust control method based on the μ-synthesis. This μ-synthesis eliminates the conservativeness and improves the control efficiency comparing with the H_∞ method. Finally, to test the control method, we simulate an IFGMSS mission in which the satellite runs in a sun synchronous circular orbit with an altitude of 300 km. The simulation results show the effectiveness of the robust control method. The performances of the closed-loop system with the μ-controller are tested by the μ-analysis. It has found that the nominal performance, the robust stability and the robust performance are all achieved. The transient simulation results further prove the control response is fast and the accuracy of the relative position meets the demand of the gravity measurement.     

A simulation of Martian gravity field recovery by using a near equatorial orbiter

An improvement to the Martian gravity field may be achieved by means of future orbiting spacecraft with small eccentricity and low altitude exemplified through a newly proposed mission design that may be tested in upcoming reconnaissance of Mars. Here, the near equatorial orbital character (with an inclination approximating 10°, eccentricity as 0.01 and semi-major axis as 4000 km) is considered, and its contribution to Martian gravity field solution is analyzed by comparing it with a hypothetical polar circular orbiter. The solution models are evaluated in terms of the following viewpoints: power spectra of gravity field coefficients, correlations of low degree zonal coefficients, precise orbit determination, and error distribution of both Mars free air gravity anomaly and areoid. At the same time, the contributions of the near equatorial orbiters in low degree zonal coefficients time variations are also considered. The present results show that the near equatorial orbiter allows us to improve the accuracy of the Martian gravity field solution, decrease correlation of low degree zonal coefficients, retrieve much better time variable information of low degree zonal coefficients, improve precise orbit determination, and provide more accurate Mars free air gravity anomaly and areoid around the equatorial region.     

The use of Gaussian equations of motions of a satellite for local gravity anomaly recovery

The orbital elements of a low Earth orbiting satellite and their velocities can be used for local determination of gravity anomaly. The important issue is to find direct relations among the anomalies and these parameters. Here, a primary theoretical study is presented for this purpose. The Gaussian equations of motion of a satellite are used to develop integral formulas for recovering the gravity anomalies. The behaviour of kernels of the integrals are investigated for a two-month simulated orbit similar to that of the Gravity field and steady-state ocean circulation explorer (GOCE) mission over Fennoscandia. Numerical investigations show that the integral formulas have neither isotropic nor well-behaved kernels. In such a case, gravity anomaly recovery is not successful due to large spatial truncation error of the integral formulas. Reformulation of the problem by combining the orbital elements and their velocities leads to an integral with a well-behaved kernel which is suitable for our purpose. Also based on these combinations some general relations among the orbital elements and their velocities are obtained which can be used for validation of orbital parameters and their velocities.     

星载GNSS确定GEO卫星轨道的积分滤波方法

采用星载全球导航卫星系统（GNSS）确定地球静止轨道（GEO）,以解决目前应用星载全球定位系统(GPS)时导航卫星可见性差的问题。以风云卫星为例,分析了未来的GNSS相对于GEO卫星的可见性,针对GEO轨道上导航接收机采样间隔较长的问题,综合轨道积分和卡尔曼滤波方法的优点,提出了确定GEO卫星轨道的积分滤波方法。并利用STK软件仿真产生所需数据,用MATLAB对提出的算法编程并进行仿真验证,结果表明,提出的方法性能优越,定轨精度较高。     

一种卫星天文自主定轨定姿方法研究

利用安装在卫星上的太阳敏感器和紫外敏感器测量出的卫星—太阳、卫星—地球和卫星一月球方向矢量，并利用雷达测高仪测出的地心距作为观测量，提出采用广义卡尔曼滤波方法实时地确定卫星绕地球飞行的轨道，同时确定出卫星的对地姿态．对自主定轨进行了数学仿真，分析比较了采样周期、轨道倾角、轨道偏心率和轨道高度等因素对定轨精度的影响．总结了其变化规律，该方法可用于提高卫星自主定轨精度．     

Impact of the orbital eccentricity on the attitude performance before and after the deorbiting phase for Alsat-1

Due to the presence of periodic forcing terms in the gravity gradient torque, orbit eccentricity may produce large response for the roll, yaw and pitch angles. This paper investigates the influence of the orbit eccentricity on the performance of the attitude determination and control subsystem (ADCS) pointing of passive Low Earth Orbit (LEO) satellites stabilized by a gravity gradient boom or having long appendages before and after the deorbiting operation. The contribution of this work is twofold. First, the satellite attitude dynamics and kinematics are modeled by introducing the orbit eccentricity in the equations of motion of a LEO satellite in order to provide the best scenario in which satellite operators can keep the nominal functionality of LEO satellites with a gravity gradient boom after the deorbiting operation. Second, a Quaternion-based Extended Kalman Filter (EKF) is analyzed when the orbit eccentricity is considered in order to determine the influence of this disturbance on the convergence and stability of the filter. The simulations in this work are based on the true parameters of Alsat-1 which is a typical LEO satellite stabilized by a gravity gradient boom. The results show that the orbit eccentricity has a big influence on the pointing system accuracy causing micro-vibrations that affect the geocentric pointing particularly after the deorbiting phase. In this case, satellites have no orbital correction option. The Quaternion-based Extended Kalman Filter analyzed in this paper, achieved satisfactory results for eccentricity values less than 0.4 with respect to pointing system accuracy. However, singularities were observed for eccentricity values greater than 0.4.     

基于会切场推力器无拖曳系统建模及参数优化

随着卫星重力测量技术的突破性进展,对航天器试验环境要求也在不断提高,航天器受到的残余扰动必须尽可能减小。作为中国将来重力场测量卫星备选主推力器的会切场推力器,其推力器的控制精度直接决定了测量的准确性。文章首先通过PID方法设计了位移模式下的无拖曳控制器,该控制器在预估阻力系数、参考质量与卫星本体的位移差、速度差等性能方面有良好的表现,在应对卫星运行时的突发情况时表现出很强的稳定性。但PID参数没有达到最优解,在此基础上对于该模型的控制精度进行优化,用遗传算法对PID控制的参数进行筛选。结果分析表明,会切场推力器的控制精度有所改善,NTW方向上的速度和位移误差均减小;推力阻力和显著减少;控制精度提高,更好地满足使用需求。     

深空卫星重力测量计划研究综述

地球卫星重力测量计划CHAMP(CHAllenging Minisatellite Payload)、GRACE(Gravity Recovery and Climate Experiment)、GOCE(Gravity field and steady-state Ocean Circulation Explorer)和月球卫星重力测量计划(Gravity Recovery and Interior Laboratory,GRAIL)的成功实施,以及下一代地球重力卫星(GRACE Follow-On)的即将发射昭示着我们将迎来一个前所未有的高精度和高空间分辨的深空卫星重力探测时代。围绕深空卫星重力测量的研究背景、必要性、可行性、卫星重力反演软件平台构建、轨道摄动和未来研究方向开展了研究论证。研究表明:深空卫星重力测量作为新世纪重力探测技术,在精化量体重力场、提高惯性导航精度、天体动力学、天体物理学和军事技术的研究,以及促进国民经济发展和提高社会效益等方面具有广泛的应用前景。     

Regional gravity field model in Pakistan area from the combination of CHAMP,GRACE and ground data using least squares collocation: A case study

This study describes a methodology of recovery of the Earth’s gravity field from CHAMP and GRACE satellites data in Pakistan using least squares collocation (LSC) based downward continuation technique. The CHAMP height anomalies and GRACE gravity disturbances derived from the observed satellite data have been used in combination solution using LSC with observed gravity values at the Earth surface. The combined covariance functions of height anomalies and/or gravity disturbances at satellite altitudes and observed gravity anomalies at Earth surface have been used as the basis for combination and downward continuation solution. The variance of predicted gravity anomalies from GRACE gravity disturbances is relatively lower than the corresponding results of gravity anomalies from CHAMP height anomalies. This fact may be attributed partly to the amplification of noise and partly to the unstable inverse transformation process of height anomalies to gravity anomalies. The impact of data error variance has been studied in the context of smoothing and noise reduction in the final solution of downward continuation using least squares collocation. The raising of data error suppresses the noise and as a result a smooth final solution is obtained. The prediction results appear to be dependent on the quality of data and goodness of combined covariance function, which are fairly comparable for the CHAMP and GRACE data. The recovered gravity field from satellite data appears to contribute mainly to medium and long wavelength parts of total gravity field spectrum. Due to flexibility of data handling in least squares collocation, this procedure is applicable to any observable of gravity field being at different altitudes and with different data spacing.     

AIUB-CHAMP02S: The influence of GNSS model changes on gravity field recovery using spaceborne GPS

The gravity field model AIUB-CHAMP02S, which is based on six years of CHAMP GPS data, is presented here. The gravity field parameters were derived using a two step procedure: In a first step a kinematic trajectory of a low Earth orbiting (LEO) satellite is computed using the GPS data from the on-board receiver. In this step the orbits and clock corrections of the GPS satellites as well as the Earth rotation parameters (ERPs) are introduced as known. In the second step this kinematic orbit is represented by a gravitational force model and orbit parameters.     

Characteristics and accuracies of the GRACE inter-satellite pointing

For almost 10 years, the Gravity Recovery and Climate Experiment (GRACE) has provided information about the Earth gravity field with unprecedented accuracy. Efforts are ongoing to approach the GRACE baseline accuracy as there still remains an order of magnitude between the present error level of the gravity field solutions and the GRACE baseline. At the current level of accuracy, thorough investigation of sensor related effects is necessary as they are one of the potential contributors to the error budget. In the science mode operations, the twin satellites are kept precisely pointed with their KBR antennas towards each other. It is the task of the onboard attitude and orbit control system (AOCS) to keep the satellites in the required formation. We analyzed long time series of the inter-satellite pointing variations as they reflect the AOCS performance and characteristics. We present significant systematic effects in the inter-satellite pointing and discuss their possible sources. Prominent features are especially related to the magnetic torquer characteristics, star cameras’ performance and KBR antenna calibration parameters. The relation between the magnetic torquer attitude control and the Earth magnetic field, impact of the different performance of the two star camera heads on the attitude control and the features due to uncertainties in the calibration parameters relating the star camera frame to K-frame are discussed in detail. Proper understanding of these effects will help to reduce their impact on the science data and subsequently increase the accuracy of the gravity field solutions. Moreover, understanding the complexity of the onboard system is essential not only for increasing the accuracy of the GRACE data but also for the development of the future gravity field satellite missions.     

Enhancing the kinematic precise orbit determination of low earth orbiters using GPS receiver clock modelling

Clock error estimation has been the focus of a great deal of research because of the extensive usage of clocks in GPS positioning applications. The receiver clock error in the spacecraft orbit determination is commonly estimated on an epoch-by-epoch basis, along with the spacecraft’s position. However, due to the high correlation between the spacecraft orbit altitude and the receiver clock parameters, estimates of the radial component are degraded in the kinematic approach. Using clocks with high stability, the predictable behaviour of the receiver oscillator can be exploited to improve the positioning accuracy, especially for the radial component. This paper introduces two GPS receiver clock models to describe the deterministic and stochastic property of the receiver clock, both of which can improve the accuracy of kinematic orbit determination for spacecraft in low earth orbit. In particular, the clock parameters are estimated as time offset and frequency offset in the two-state model. The frequency drift is also estimated as an unknown parameter in the three-state model. Additionally, residual non-deterministic random errors such as frequency white noise, frequency random walk noise and frequency random run noise are modelled. Test results indicate that the positioning accuracy could be improved significantly using one day of GRACE flight data. In particular, the error of the radial component was reduced by over 40.0% in the real-time scenario.     

NRLMSISE-00大气模型与GRACE和CHAMP卫星大气密度数据的对比分析

利用GRACE(Gravity Recovery And Climate Experiment)和CHAMP(Challenging Mini-Satellite Payload)卫星2002-2008年的大气密度数据与NRLMSISE-00大气模型密度结果进行比较,分析了模型密度误差及其特点.结果显示,NRLMSISE-00大气模型计算的密度值普遍偏大,其相对误差随经纬度变化,在高纬度相对较小;相对误差随地方时变化,在02:00LT和15:00LT左右较大,10:00LT和20:00LT左右较小.通过模型密度相对误差与太阳F_10.7指数的对比分析发现,在太阳活动低年模型相对误差最大,而在太阳活动高年相对误差较小;将模型结果分别与GRACEA/B双星和CHAMP卫星的密度数据进行比较,发现对于轨道高度更高的GRACE卫星轨道,模型相对误差更大;在地磁平静期,相对误差与地磁ap指数(当前3h)相关性不强,但是在大磁暴发生时,误差急剧增大.