基于预测滤波的捷联惯导任意双位置对准方法

针对捷联惯性导航系统(SINS,Strapdown Inertial Navigation System)在大失准角情况下的初始对准问题,建立基于加性四元数误差模型的非线性滤波方程,并提出一种基于模型预测滤波(MPF, Model Predictive Filter)与扩展卡尔曼滤波(EKF, Extended Kalman Filter)相结合的地面任意双位置初始对准方法.该方法将部分惯性器件误差作为模型误差,在线实时估计并修正系统模型,提高了状态估计的精度,并克服了将模型误差假设为高斯白噪声的局限性.半物理仿真结果表明,该方法有效提高了SINS姿态误差角的估计精度,而且也降低了系统状态变量的维数,提高了对准解算的实时性.     

无轨道动力学模型的增强型GPS实时星间相对定位方法

虽然基于卫星轨道动力学模型的增强型GPS实时星间相对定位可获得厘米级精度,但通常在一般条件下难以获取高精度轨道动力学模型.本文提出了一种无轨道动力学模型条件下的高精度实时相对定位方法,即利用GPS三差观测方程获得状态矢量的时间转移矩阵,采用自适应UKF(Unsented Kalman Filter,无迹卡尔曼滤波器)实时估计星间相对位置.仿真结果表明,相比传统的序贯点估计方法分米级的相对定位精度,新方法在星间距离小于5 km时,其相对定位精度可达厘米量级.     

基于混沌粒子群优化的北斗/GPS组合导航选星算法

全球卫星导航系统（GNSS）接收机在接收信号的过程中会受到诸如建筑物遮挡、信号干扰等因素的影响,无法得到全部可见星。为减轻多星座组合接收机的处理负担,研究利用部分可见卫星进行定位的快速选星算法,提出了一种基于混沌粒子群优化（CPSO）的北斗/GPS组合导航选星算法。首先,对当前历元时刻可见卫星进行连续编码,按照选星数目分组,每个分组视为一个粒子。然后,通过混沌映射初始化粒子种群,选取几何精度因子（GDOP）作为评价粒子优劣的适应度函数;粒子通过粒子群优化算法的速度-位移模型更新自身位置,逐渐趋近空间卫星几何分布较好的卫星组合全局最优解。最后,采集北斗/GPS实际数据对选星算法进行仿真验证和性能比较,结果表明,所提算法在选星颗数多于5颗时,单次选星耗时为遍历法选星的37.5%,选星结果的几何精度因子计算误差在0~0.6之间。该算法可适用于北斗/GPS组合导航定位不同选星颗数的情况。      

Iridium/ORBCOMM机会信号融合定位技术

在全球导航卫星系统(GNSS)不可用情况下，低地球轨道(LEO)卫星机会信号(SOP)定位技术是一种有效的导航定位解决方案。单LEO星座机会信号定位技术面临星座构型不足或可见卫星偏少等问题，多LEO星座机会信号融合定位技术可有效解决该问题。通过分析瞬时多普勒定位原理，建立了Iridium/ORBCOMM机会信号融合定位模型，引入基于Helmert方差估计的加权最小二乘算法进一步提高定位精度。实测数据表明:基于Helmert方差估计的Iridium/ORBCOMM机会信号融合定位精度优于70 m，验证了多LEO星座机会信号融合定位的可行性和有效性。      

GPS 时间同步技术及其在航天领域的应用前景

导航星全球定位系统(简称 GPS)是一项全新的卫星导航定位的高技术工程,是全球可共享的技术资源。它具有全球覆盖、全天候工作、全天24小时连续而实时地为无限多个用户提供高精度七维(三维位置、三维速度和时间)信息的能力。实时定位精度可达10m 以内,事后处理的定位精度可达1m 左右。测速精度可达0.1—0.01m/s。时间传递和时间同步精度可达10—1ns(10~(-9)s).本文简要介绍了 GPS 的概况和国内外开发研究的动态。着重就 GPS 时间传递的原理、力法和误差做了较深入的分析和论证。报道了国内首次 GPS 共视法时间同步大型试验的情况和取得的结果。最后,作者专门论述了 GPS 技术在航天领域的应用前景。     

地标自主导航的高精度地标库建立方法

利用地标信息的卫星自主导航方法是通过星载光学摄像机实时拍摄地标景象来计算卫星的位置、速度等参数,具有完全自主的特点,并且适用范围广泛,可用于能够周期性获得地面图像的航天器,但受到地标选取方法以及地标库建立策略的制约,影响到导航精度,阻碍了推广使用。针对上述问题,提出了建立全局及高性能的地标库,即采用全球地标控制点选取的原则选取性能优良的全球地标样本,并自动生成全局地标库。通过计算机仿真验证了该方法的有效性,结果表明,使用全球地标控制点选取地标样本,并自动生成全局地标库,可以快速而准确地为航天器进行高精度的自主导航,导航位置误差约为99 m,速度误差约为0.08 m/s。      

一种用于GPS/DR组合定位的非线性滤波算法

建立了适用于车辆导航系统的基于UKF(Unscented Kalman Filter)的GPS/DR(Global Positioning System/Dead Reckoning)组合定位滤波模型及算法.针对系统状态方程为线性、观测方程为非线性的特点,提出了一种将UKF和EKF(Extended Kalman Filter)相结合的非线性滤波算法.结合后的算法和原有UKF算法相比减少了在时间更新阶段的运算量,并且由于采用基于Unscented变换的思想来处理系统观测方程的非线性问题,避免了EKF引入的线性化误差,提高了滤波精度.仿真结果证明:算法在减少运算量的同时,仍具有较高的滤波精度,且明显优于EKF,因而能够满足车辆导航系统占用资源少、滤波精度高的要求.     

大方位失准角MIMU空中对准建模及半物理仿真

由于快速性的要求,微小型无人机不经过地面精确初始对准就升空作业,因此MIMU(Micro Inertial Measurement Unit)空中对准在大失准角下进行. 为了提高微小型无人机空中的反应速度和作业精度,把非线性误差部分作为状态变量,建立MIMU在大方位失准角下无需小角度近似的空中对准的线性模型,同时为解决噪声不确定导致滤波器发散的问题,提出将AKF (Adaptive Kalman Filter)应用在GPS(Global Positioning System)辅助MIMU的空中对准中,半物理仿真结果证实其取得了比基于非线性误差模型的EKF(Extended Kalman Filter)精度高且速度快的结果,不仅使MIMU的方位失准角由60° 快速下降到2° 左右,且所需时间仅为EKF的67%.     

一种模型误差鲁棒的卫星姿态确定方法

提出一种结合非线性预测滤波和二阶插值滤波实现基于星光/陀螺的高精度姿态确定的新算法.该算法用非线性预测滤波估计模型误差,再对补偿后的模型用高精度的二阶插值滤波来估计姿态参数.解决了在卫星实际运行中难以获得姿态确定系统的精确动力学模型,采用传统EKF(Extended Kalman Filter)将模型误差作为零均值白噪声处理,导致滤波精度降低甚至发散的问题.同时,二阶插值滤波将非线性模型按照二阶近似,无需计算函数偏导数,得到高精度的卫星姿态估计.仿真验证了该方法能有效地实时估计并补偿模型误差,提高了姿态估计的精度,且估计精度受滤波周期的影响不大,从而验证了算法的鲁棒性和有效性.      

基于ORBCOMM卫星机会信号的定位技术

为摆脱对全球导航卫星系统（GNSS）的依赖，克服其有意或无意干扰情况下无法工作等问题，可采用机会信号（SOP）实现定位，低轨卫星机会信号具备信号功率高、覆盖性广及无需增建基础设施等优点。提出了利用轨道通信卫星（ORBCOMM）系统实现天基机会信号定位。通过对ORBCOMM卫星机会信号的通信体制进行深入研究，实现了利用ORBCOMM卫星机会信号获取多普勒测量信息，建立了瞬时多普勒定位及其几何精度因子的数学模型，并采用卫星TLE数据结合轨道预测模型获得的卫星轨道信息实现ORBCOMM卫星机会信号定位。实测结果表明:利用ORBCOMM卫星机会信号可实现精度优于140 m的定位。研究成果对基于天基机会信号定位技术的理论研究及应用具有重要意义。      

The impact of orbital errors on the estimation of satellite clock errors and PPP

Precise satellite orbit and clocks are essential for providing high accuracy real-time PPP (Precise Point Positioning) service. However, by treating the predicted orbits as fixed, the orbital errors may be partially assimilated by the estimated satellite clock and hence impact the positioning solutions. This paper presents the impact analysis of errors in radial and tangential orbital components on the estimation of satellite clocks and PPP through theoretical study and experimental evaluation. The relationship between the compensation of the orbital errors by the satellite clocks and the satellite-station geometry is discussed in details. Based on the satellite clocks estimated with regional station networks of different sizes (∼100, ∼300, ∼500 and ∼700 km in radius), results indicated that the orbital errors compensated by the satellite clock estimates reduce as the size of the network increases. An interesting regional PPP mode based on the broadcast ephemeris and the corresponding estimated satellite clocks is proposed and evaluated through the numerical study. The impact of orbital errors in the broadcast ephemeris has shown to be negligible for PPP users in a regional network of a radius of ∼300 km, with positioning RMS of about 1.4, 1.4 and 3.7 cm for east, north and up component in the post-mission kinematic mode, comparable with 1.3, 1.3 and 3.6 cm using the precise orbits and the corresponding estimated clocks. Compared with the DGPS and RTK positioning, only the estimated satellite clocks are needed to be disseminated to PPP users for this approach. It can significantly alleviate the communication burdens and therefore can be beneficial to the real time applications.     

高精度星载SAR多普勒参数估算方法

提出了一种基于空间坐标转换,利用卫星位置、速度参数精确估算星载SAR(Synthetic Aperture Radar)全观测带多普勒参数的方法.利用卫星速度、位置,通过星载SAR空间几何模型和坐标转换关系,建立SAR图像中斜距同卫星下视角之间的四次方程,解出下视角并进一步计算出该斜距处的多普勒参数值.仿真结果表明,该方法在无卫星位置、速度误差情况下估算精度达到0.02Hz(多普勒中心频率)和2×10^-4Hz/s(多普勒调频率);存在卫星位置测量误差(300m)以及速度测量误差(0.3m/s)的情况下,估算精度达到0.8Hz(多普勒中心频率)和0.07Hz/s(多普勒调频率).该方法适用于单星SAR以及分布式SAR高精度多普勒参数的估算.     

中轨道Walker导航星座在轨备份方案优化设计

针对中轨道Walker导航星座在轨备份方案的优化问题，首先提出了在轨备份星轨位优化设计方法，考虑星座运行期间在轨备份星与工作卫星存在共同提供服务的情况，选取PDOP值和可见卫星数作为轨位优化指标，建立了轨位优化模型，并基于NSGA-Ⅱ算法对不同轨位下在轨备份星对导航星座服务性能的提升效果进行仿真分析；然后，基于在轨备份星轨位的优化结果，建立在轨备份星替换的轨道机动模型，并综合考虑速度增量和替换时间，确定了以替换时间最少为优化目标的在轨备份星替换方案。结果表明，提出的在轨备份方案能有效满足备份星的设计需求，增强导航星座的服务性能，实现故障卫星的快速替换，可为导航星座备份星的建设提供借鉴。     

Accuracy of predicted earth's artificial satellite orbits

The accuracy of predicted satellite positions depends on many factors: the accuracy of initial conditions (observations), the accuracy of the adopted force model, the accuracy of the theory used in calculations and the time interval. The paper proposes an analytical algorithm for the inverse transformation from osculating to mean orbital elements and presents some results of an analysis of the influence of geopotential errors on predicted satellite position errors.     

基于线性协方差方法的交会对接误差分析

将线性协方差分析方法和蒙特卡罗仿真相结合,按交会任务和飞行特征把交会过程分为变轨飞行、自由飞行和中途速度修正三种特征段,研究了状态误差的传播规律和交会过程中各种误差对交会对接精度的影响。在变轨飞行段,分析了追踪航天器的姿态误差、控制系统性能状态估计误差,以及目标航天器轨道摄动对状态误差传播的影响。在自由飞行段,分析了追踪航天器估计状态误差的先验值和测轨误差对状态误差传播的影响。在中途速度修正段,分析了追踪航天器姿态误差和控制系统性能误差对状态误差传播的影响。仿真结果表明,误差分析方法设计合理,可以指导交会对接的轨道设计工作,能对已经设计好的交会策略进行误差分析和设计验证。     

Simulation of the land surface temperature from moon-based Earth observations

The land surface temperature (LST) is a key parameter for the Earth’s energy balance. As a natural satellite of the Earth, the orbital of the moon differs from that of current Earth observation satellites. It is a new way to measure the land surface temperature from the moon and has many advantages compared with artificial satellites. In this paper, we present a new method for simulating the LST measured by moon-based Earth observations. Firstly, a modified land-surface diurnal temperature cycle (DTC) method is applied to obtain the global LST at the same coordinated universal time (UTC) using the Moderate Resolution Imaging Spectroradiometer (MODIS) LST products. The lunar elevation angles calculated using the ephemeris data (DE405) from the Jet Propulsion Laboratory (JPL) were then applied to simulate the Earth coverage observed from the moon. At the same time, the modified DTC model was validated using in situ data, MODIS LST products, and the FengYun-2F (FY-2F) LST, respectively. The results show that the fitting accuracy (root-mean-square error, RMSE) of the modified DTC model is not greater than 0.72?°C for eight in situ stations with different land cover types, and the maximum fitting RMSE of the modified model is smaller than that of current DTC models. By the comparison of the simulated LST with MODIS and FY-2F LST products, the errors of the results were feasible and accredited, and the simulated global LST has a reasonable spatiotemporal distribution and change trend. The simulated LST data can therefore be used as base datasets to simulate the thermal infrared imagery from moon-based Earth observations in future research.     

Time synchronization of new-generation BDS satellites using inter-satellite link measurements

Autonomous satellite navigation is based on the ability of a Global Navigation Satellite System (GNSS), such as Beidou, to estimate orbits and clock parameters onboard satellites using Inter-Satellite Link (ISL) measurements instead of tracking data from a ground monitoring network. This paper focuses on the time synchronization of new-generation Beidou Navigation Satellite System (BDS) satellites equipped with an ISL payload. Two modes of Ka-band ISL measurements, Time Division Multiple Access (TDMA) mode and the continuous link mode, were used onboard these BDS satellites. Using a mathematical formulation for each measurement mode along with a derivation of the satellite clock offsets, geometric ranges from the dual one-way measurements were introduced. Then, pseudoranges and clock offsets were evaluated for the new-generation BDS satellites. The evaluation shows that the ranging accuracies of TDMA ISL and the continuous link are approximately 4?cm and 1?cm (root mean square, RMS), respectively. Both lead to ISL clock offset residuals of less than 0.3?ns (RMS). For further validation, time synchronization between these satellites to a ground control station keeping the systematic time in BDT was conducted using L-band Two-way Satellite Time Frequency Transfer (TWSTFT). System errors in the ISL measurements were calibrated by comparing the derived clock offsets with the TWSTFT. The standard deviations of the estimated ISL system errors are less than 0.3?ns, and the calibrated ISL clock parameters are consistent with that of the L-band TWSTFT. For the regional BDS network, the addition of ISL measurements for medium orbit (MEO) BDS satellites increased the clock tracking coverage by more than 40% for each orbital revolution. As a result, the clock predicting error for the satellite M1S was improved from 3.59 to 0.86?ns (RMS), and the predicting error of the satellite M2S was improved from 1.94 to 0.57?ns (RMS), which is a significant improvement by a factor of 3–4.     

基于星间距离测量的高精度自主导航

研究利用地球卫星和月球卫星之间的测距信息进行自主导航的方法.基于三体摄动轨道动力学方程和星间测距信息,可以同时确定参与导航的地球卫星和月球卫星的绝对位置;但是在初始位置误差较大的情况下,导航系统的定位性能会受到影响.为了解决这一问题,提出基于"星间测距+紫外导航敏感器"的组合导航方法.采用该导航方法,能够在初始位置误差和紫外导航敏感器测量误差较大的情况下实现高精度自主导航.基于Cramer-Rao下界(CRLB)分析了组合导航系统的性能,并通过数学仿真验证了该导航方法的有效性.     

“龙江2号”月球轨道微卫星定轨分析

微纳卫星深空探测任务中,通常所分配的测控资源有限,因此有必要对有限测控资源条件下微纳卫星的定轨精度进行分析。以微纳卫星深空探测为背景,采用"龙江2号"微卫星的轨道测量数据对其定轨精度进行了分析。"龙江2号"微卫星只有USB轨道测量数据,且环月段测控资源相对紧张,每天有两站跟踪,共约3～4 h的轨道测量数据。首先介绍了"龙江2号"微卫星飞行任务及其飞行过程中影响测定轨的因素;其次给出了定轨的动力学模型,对微卫星地月转移段的定轨精度进行了分析;最后通过分析摄动力、动量轮卸载以及数据弧段长度的影响,给出了微卫星环月阶段所使用的定轨策略,并通过重叠弧段比较的模式,给出了微卫星环月段的定轨精度。研究结论可以为后续微纳卫星深空探测任务提供有益参考。     

Orbital error propagation considering atmospheric density uncertainty

Due to the influence of various errors, the orbital uncertainty propagation of artificial celestial objects while orbit prediction is required, especially in some applications such as conjunction analysis. In the orbital error propagation of artificial celestial objects in low Earth orbits (LEOs), atmospheric density uncertainty is one of the important factors that require special attention. In this paper, on the basis of considering the uncertainties of position and velocity, the atmospheric density uncertainty is also taken into account to further investigate the orbital error propagation of artificial celestial objects in LEOs. Artificial intelligence algorithms are introduced, the MC Dropout neural network and the heteroscedastic loss function are used to realize the correction of the empirical atmospheric density model, as well as to provide the quantification of model uncertainty and input uncertainty for the corrected atmospheric densities. It is shown that the neural network we built achieves good results in atmospheric density correction, and the uncertainty quantization obtained from the neural network is also reasonable. Moreover, using the Gaussian mixture model - unscented transform (GMM-UT) method, the atmospheric density uncertainty is taken into account in the orbital uncertainty propagation, by adding a sampled random term to the corrected atmospheric density when calculating atmospheric density. The feasibility of the GMM-UT method considering atmospheric density uncertainty is proved by the further comparison of abundant sampling points and GMM-UT results (with and without considering atmospheric density uncertainty).