Impact of covariance information of kinematic positions on orbit reconstruction and gravity field recovery

Gravity missions are equipped with onboard Global Positioning System (GPS) receivers for precise orbit determination (POD) and for the extraction of the long wavelength part of the Earth’s gravity field. As positions of low Earth orbiters (LEOs) may be determined from GPS measurements at each observation epoch by geometric means only, it is attractive to derive such kinematic positions in a first step and to use them in a second step as pseudo-observations for gravity field determination. The drawback of not directly using the original GPS measurements is, however, that kinematic positions are correlated due to the ambiguities in the GPS carrier phase observations, which in principle requires covariance information be taken into account. We use GRACE data to show that dynamic or reduced-dynamic orbit parameters are not optimally reconstructed from kinematic positions when only taking epoch-wise covariance information into account, but that essentially the same orbit quality can be achieved as when directly using the GPS measurements, if correlations in time are taken into account over sufficiently long intervals. For orbit reconstruction covariances have to be considered up to one revolution period to avoid ambiguity-induced variations of kinematic positions being erroneously interpreted as orbital variations. For gravity field recovery the advantage is, however, not very pronounced.     

基于CAD二次开发的飞行器隐身性能计算方法

在飞行器同雷达对抗时的不同威胁情况下,分析了在雷达仰视照射下飞行器的隐身性能。通过绘制雷达仰角图,得到了目标飞行器在某典型雷达照射下的暴露与隐身区域,从而量化了目标飞行器的隐身性能。通过CATIA二次开发技术实现了CAD模型的便捷修改,为优化外形隐身特性提供支持。最后通过一个具体算例验证种方法的效果,并说明方法具有实用意义。     

Bayesian linking of geosynchronous orbital debris tracks as seen by the Large Synoptic Survey Telescope

We describe a Bayesian sampling model for linking and constraining orbit models from angular observations of “streaks” in optical telescope images. Our algorithm is particularly suited to situations where the observation times are small fractions of the orbital periods of the observed objects or when there is significant confusion of objects in the observation field. We use Markov Chain Monte Carlo to sample from the joint posterior distribution of the parameters of multiple orbit models (up to the number of observed tracks) and parameters describing which tracks are linked with which orbit models. Using this algorithm, we forecast the constraints on geosynchronous (GEO) debris orbits achievable with the planned Large Synoptic Survey Telescope (LSST). Because of the short 15 s exposure times, preliminary orbit determinations of GEO objects from LSST will have large and degenerate errors on the orbital elements. Combined with the expected crowded fields of GEO debris it will be challenging to reliably link orbital tracks in LSST observations given the currently planned observing cadence.     

月球探测软着陆与采样返回段弹道确定

针对月球探测中软着陆与采样返回段弹道计算问题,提出用数值逼近弹道确定方法。通过B样条对探测器状态进行建模,进而综合全弧段数据进行统计定轨的方法。由于样条法良好的数值逼近性能,使得该方法对探测器弹道异常复杂情况下的状态确定较为有效。对嫦娥三号探测器动力软着陆弧段进行了仿真与实测数据处理。分析了采样返回段的基本动力学与控制特征,为后续的嫦娥五号探测器的软着陆及其采样返回提供初步的可行弹道计算方法。在嫦娥三号探测器动力落月段实测数据处理中,通过评估,该段弹道确定精度优于100 m,其弹道末点与NASA的月球勘测轨道器（LRO）给出的结果差异优于50 m,证实了文章提出的软着陆弹道确定方法的有效性。     

航天器自主交会对接的视觉相对导航方法

针对航天器自主交会对接在最终逼近阶段的相对导航问题,研究了基于特征光标的航天器视觉相对导航方法。首先建立了小角度假设条件下的相对导航模型。给出了航天器相对位置和采用欧拉角描述的相对姿态的近似解析算法。其次利用相对姿态的解析解作为观测信息,采用降阶的无迹卡尔曼滤波方法对目标航天器的姿态角和惯性角速度进行了估计,间接得到了相对姿态和相对姿态角速度的估计值。数值仿真表明提出的解析算法能够有效获得相对位置和相对姿态测量信息,通过降阶滤波可明显提高相对姿态确定精度。     

基于BP神经网络的自主定轨自适应Kalman滤波算法

针对Sage Husa自适应滤波方法存在的窗函数开窗大小选择问题,提出一种基于BP神经网络学习估计系统协方差矩阵的自适应Kalman滤波算法。该算法以Kalman滤波预测残差向量作为网络输入,通过网络分段离线学习确定预测残差向量与预测残差协方差矩阵间的非线性关系,自适应地估计Kalman滤波系统协方差矩阵。将其应用到自主定轨系统,仿真结果表明利用本文算法自主定轨60天星座平均URE误差小于1.9米,且能够快速跟踪到系统噪声的突变,较Kalman滤波方法和Sage Husa自适应滤波方法具有更好的性能。     

嫦娥二号探测器轨道确定与支持

针对嫦娥二号探测器直接进入地月转移轨道、距月面100km高度捕获月球、完成既定任务后飞往日地第二拉格朗日平动点等飞行轨道方面的新特点,分析了定轨预报策略,利用事后精密轨道,全面评估了关键变轨点定轨预报和变轨后快速定轨的精度,其中,近月制动前3h定轨预报至近月点的位置误差为1km,速度误差为 0.3m/s 。利用不同月球引力场模型进行环月轨道精密定轨,根据实测数据残差分析和精密星历比对的结果,采用SGM100h引力场模型的定轨残差均方根最大。此外,针对嫦娥二号扩展任务,分析了不同测轨条件下的定轨精度,测量数据残差分析结果表明,在扩展任务中途修正前的定轨弧段内,测距、时延和时延率数据的残差分别为5m,5ns和1ps/s。     

典型过失速机动运动规律建模研究

对典型过失速机动的运动特性进行了分析,根据其运动特点将过失速机动的运动过程分解为定直飞行、快速俯仰和绕速度矢量旋转三种基本运动的组合。选取基本运动的关键参数,并针对关键参数建立过失速机动的运动模型。模型计算结果和飞行仿真结果的比较表明,所建立的运动模型在过失速机动中的应用是可行的。     

一种改进的IMU加表标定模型及快速标定方法

针对MEMS产品工程应用的需求,提出了一种改进的IMU加表标定模型,更清晰地分离出与安装误差有关的投影关系阵以及与器件误差有关的刻度转化阵。并且,在此模型的基础上给出了不需转台的快速标定方法。投影关系阵不变,该方法减少了观测量,只需6位置来估计加速度计的零位和标度因数,同时进行了现场误差补偿。     

卫星编队飞行的相对轨道自主确定研究

为了研究卫星编队飞行相对轨道的自主确定,基于相对轨道根数建立编队卫星间的相对运动方程,利用测量所得到的星间距离和方位信息作为观测量。不同于目前广泛采用的扩展卡尔曼滤波算法,设计Unscented Kalman Filter(UKF)算法实现卫星编队飞行的相对轨道自主确定。仿真结果表明这种相对轨道自主确定方案能获得较高的定轨精度。