两种GPS广播星历参数算法的比较

介绍了新GPS广播星历参数,分析了新旧GPS广播星历参数及其用户算法的区别,导出了卫星位置矢量对新GPS广播星历参数的偏导数,并与旧GPS广播星历参数偏导数进行了简要的比较分析.计算结果表明,两种GPS广播星历参数算法都能很好地描述MEO卫星的轨道特性,且广播星历参数拟合算法的精度损失仅为数厘米.      

GEO卫星广播星历的拟合参数设计

基于第一类无奇点要素的拟合参数和5°倾角旋转法,可以实现GEO发播电文与MEO和IGSO的电文一致性,但是没有完全避免数学奇异性.文章从无奇点要素的选取和主要摄动力的影响分析两方面出发,提出了一套16参数的GEO星历表示法.首先,利用同步要素描述卫星的星下点运动轨迹;其次,摄动参数中顾及了地球引力位田谐项引起的定点漂移.拟合试验表明,2h星历拟合的用户距离误差的均方根优于0.01m.     

Earth–Moon low energy trajectory optimization in the real system

The problem of the Earth–Moon low energy trajectory optimization in the real system (the model defined by the JPL ephemeris DE405) is considered in this paper. First, this problem is investigated in the model of circular restricted three-body problem, since the fuel consumption is closely related to the Jacobi integral of the transfer trajectory, a method based on Jacobi integral is proposed and eight optimal trajectories are obtained. These optimal trajectories provide initial information (the flight time and the braking velocity impulse) to search the optimal low energy trajectories in the real system through optimization techniques. Considering the merit and drawback of particle swarm optimization and differential evolution algorithm in solving the space trajectory problem, an improved cooperative evolutionary algorithm is put forward. Result shows that the low energy trajectories in the real system are more fuel-efficient than the corresponding ones under the circular restricted three-body problem.     

基于自适应滤波的编队卫星实时相对定轨

提出基于自适应滤波的编队卫星实时相对定轨算法,利用2005-12-09—10两颗GRACE（Gravity Recovery and Climate Experiment）卫星的GPS（Global Positioning System）实测数据进行实时相对定轨试验计算,采用JPL（Jet Propulsion Laboratory）轨道对试验结果外部检核,结果表明:①自适应滤波相对定轨通过自适应因子,可以较好地平衡编队卫星的观测信息和相对动力学信息,其相对定轨结果精度优于Kalman滤波相对定轨结果;②自适应滤波相对定轨结果随着星间基线缩短而精度提高;③两颗GRACE卫星采用单频伪距和广播星历进行自适应滤波相对定轨,可以得到精度优于6cm的星间基线。     

一种24参数GPS广播星历方案研究

针对进一步降低广播星历拟合误差,文章提出了一种改进的GPS 24参数广播星历,并给出了广播星历的改进算法。实验结果表明:对于4h拟合长度,GPS全星座所有卫星、所有弧段的三维位置拟合性能均有显著提升,三维位置拟合误差的均方根值和最大值分别降低87%和47%,以三维位置拟合误差的均方根值作为评价指标,改进的GPS 24参数广播星历拟合精度达到毫米级。     

一类适用于各种轨道类型的导航卫星广播星历研究

导航卫星广播星历是卫星导航定位的基础,要求有较高的精度和实时性。GPS和GLONASS广播星历形式适合于某一类卫星轨道,但存在一定局限性。利用切比雪夫多项式拟合导航卫星轨道,拟合系数作为广播星历的参数发布,这是一种新的广播星历形式,可作为各种轨道类型的导航卫星的广播星历参数,而且具有精度高、计算速度快等特点。     

一种大椭圆轨道卫星星历预报方法

针对大偏心率轨道长时间星历预报,稠密输出星历的低效率问题,提出了一种新的星历预报方法。该方法通过建立卫星摄动力模型,对卫星运动方程进行数值积分来获取卫星星历。该方法的轨道积分器以标准的8阶Adams-Cowell多步法为基础,针对大偏心率轨道的特性,在一个轨道周期的不同时段内采用不同的积分步长,同时引入了用于生成小间隔等间距卫星星历的插值公式。该方法降低了卫星运动方程右函数的计算次数,尤其适用于需要稠密输出卫星星历的情况。以STK(Satellite Tool Kit,卫星工具包)的HPOP(High Precision Orbit Propagator,高精度轨道预报)模型为验证基准,通过数学仿真校验了该方法的有效性和精度。算例表明,该方法在预报精度损失很小的情况下,计算时间远小于标准的多步法和HPOP模型。     

一种高精度实时GPS卫星钟差预测算法

为获得高精度实时GPS卫星钟差,文章提出一种基于多项式和最小二乘支持向量机(Least Squares Support Vector Machines,LS-SVM)相结合的钟差预报方法.该方法采用国际GNSS服务发布的超快速观测星历建模进行短期预报,首先根据卫星钟的物理特性用附有周期项的多项式模型进行拟合以提取趋势项和周期项,然后用LS-SVM对多项式拟合残差进行建模预报,最后将预报结果加上趋势项和周期项,得到最终的钟差预报值.试验结果表明,所提算法能够实时有效地对GPS卫星钟差进行预报,且精度优于超快速预报星历.     

Refinement of global ionospheric coefficients for GNSS applications: Methodology and results

We developed the methodology for the optimal estimation of global ionospheric coefficients of the current Global Navigation Satellite Systems (GNSSs), including the eight- and ten-parameter Klobuchar-like as well as NeQuick models. The ionospheric coefficients of those correction models are calculated from two sets of globally distributed tracking stations of the International GNSS Services (IGS). Performance of the re-estimated Klobuchar-like and NeQuick coefficients are validated during 2002–2014 over the continental and oceanic areas, respectively. Over the continental areas, GPS TECs derived from 40 ground GPS receivers are selected as reference. The eight-, ten-parameter Klobuchar-like and NeQuick models can mitigate the ionospheric delay by 65.8, 67.3 and 75.0%, respectively. Over the global oceans, the independent TECs derived from Jason-1&2 altimeters are used as reference. The re-estimated ionospheric correction models can mitigate 56.1–66.7% of the delay errors. Compared to the original GPS Ionospheric Correction Algorithm (ICA), performance of those eight-, ten-parameter Klobuchar-like and NeQuick models has improved 3.4, 5.9 and 13.4% during the whole test period, respectively. The methodology developed here takes the advantage of high-quality ionospheric TECs derived from the global network of GNSS receivers. The re-estimated ionospheric coefficients can be used as precise ionospheric products to monitor and assess GNSS broadcast ionospheric parameters and to improve the performance of various single-frequency GNSS applications.     

Application of an evolution strategy in planetary ephemeris modeling

Classical planetary ephemeris construction comprises three major steps which are to be performed iteratively: numerical integration of coupled equations of motion of a multi-body system (propagator step), reduction of observations (reduction step), and optimization of model parameters (adjustment step). In future, this approach may become challenged by further refinements in force modeling (e.g. inclusion of much more significant minor bodies than in the past), an ever-growing number of planetary observations (e.g. the vast amount of spacecraft tracking data), and big data issues in general. In order to circumvent the need for both the inversion of normal equation matrices and the determination of partial derivatives, and to prepare the ephemeris for applications apart from stand-alone solar-system planetary orbit calculations, here we propose an alternative ephemeris construction method. The main idea is to solve it as an optimization problem by straightforward direct evaluation of the whole set of mathematical formulas, rather than to solve it as an inverse problem with all its tacit mathematical assumptions and potential numerical difficulties. The usual gradient search is replaced by a stochastic search, namely an evolution strategy, the latter of which is perfect for the exploitation of parallel computing capabilities. Furthermore, this new approach allows for multi-criteria optimization and time-varying optima. These issues will become important in future once ephemeris construction is just one part of even larger optimization problems, e.g. the combined and consistent determination of a generalized physical state (orbit, size, shape, rotation, gravity, $\dots$) of celestial bodies (planets, satellites, asteroids, or comets), and/or if one seeks near real-time solutions. Here, we outline the general idea and exemplarily optimize high-correlated asteroidal ring model parameters (total mass and heliocentric radius), and individual asteroid masses, based on simulated observations.