一种星敏感器-陀螺组合定姿的实时在轨标定方法

研究了一种星敏感器一陀螺组合定姿方式中的姿态敏感器误差的实时在轨标定方法。首先,选择直观的欧拉角作为姿态描述参数,根据星敏感器和陀螺的测量原理建立星敏感器一陀螺在轨标定的测量方程和状态方程,并以此建立数学模型。其次,采用简单高效的EKF（ExtendedKalmanFilter,扩展卡尔曼滤波）作为估值算法,进行了在轨标定数值仿真。对于航天器姿态定向中出现的姿态角和星敏感器安装角之间的耦合问题,通过在特定姿态通道上施加简单姿态机动实现了解耦。数值结果表明,该实时在轨标定方法,尤其是所提出的姿态角和星敏感器安装角解耦策略,可以实现对航天器姿态的实时精确估计以及对星敏感器安装误差、陀螺常值漂移和相关漂移等误差的实时在轨标定。该方法可用于航天器姿态测量设备的实时在轨标定和航天器姿态的高精度实时确定。     

基于北斗的车载接收机航向测量技术研究

研究一种BD2双天线的双频基线测量算法,实现了在车载动态条件下利用BD2卫星进行载体航向实时测量,满足军用车辆、无人机、船舶舰艇等载体高精度、快速实时、连续稳定长时间工作、低成本的要求,对算法实现进行了试验验证,取得了理想的成效,具有较强的工程实用价值。     

DORIS geodesy: A dynamic determination of geocentre location

Geoscience Australia contributed a multi-satellite, multi-year weekly time series to the International DORIS Service combined submission for the construction of International Terrestrial Reference Frame 2008 (ITRF2008). This contributing solution was extended to a study of the capability of DORIS to dynamically estimate the variation in the geocentre location. Two solutions, comprising different constraint configurations of the tracking network, were undertaken. The respective DORIS satellite orbit solutions (SPOT-2, SPOT-4, SPOT-5 and Envisat) were verified and validated by comparison with those produced at the Goddard Space Flight Center (GSFC), DORIS Analysis Centre, for computational consistency and standards. In addition, in the case of Envisat, the trajectories from the GA determined SLR and DORIS orbits were compared. The results for weekly dynamic geocentre estimates from the two constraint configurations were benchmarked against the geometric geocentre estimates from the IDS-2 combined solution. This established that DORIS is capable of determining the dynamic geocentre variation by estimating the degree one spherical harmonic coefficients of the Earth’s gravity potential. It was established that constrained configurations produced similar results for the geocentre location and consequently similar annual amplitudes. For the minimally constrained configuration Greenbelt–Kitab, the mean of the uncertainties of the geocentre location were 2.3, 2.3 and 7.6 mm and RMS of the mean uncertainties were 1.9, 1.2 and 3.5 mm for the X, Y and Z components, respectively. For GA_IDS-2_Datum constrained configuration, the mean of the uncertainties of the geocentre location were 1.7, 1.7 and 6.2 mm and RMS of the mean uncertainties were 0.9, 0.7 and 2.9 mm for the X, Y and Z components, respectively. The mean of the differences of the two DORIS dynamic geocentre solutions with respect to the IDS-2 combination were 1.6, 4.0 and 5.1 mm with an RMS of the mean 21.2, 14.0 and 31.5 mm for the Greenbelt–Kitab configuration and 4.1, 3.9 and 4.3 mm with an RMS 8.1, 9.0 and 28.6 mm for the GA_IDS-2_Datum constraint configuration. The annual amplitudes for each component were estimated to be 5.3, 10.8 and 11.0 mm for the Greenbelt–Kitab configuration and 5.3, 9.3 and 9.4 mm for the GA_IDS-2_Datum constraint configuration. The two DORIS determined dynamic geocentre solutions were compared to the SLR determined dynamic solution (which was determined from the same process of the GA contribution to the ITRF2008 ILRS combination) gave mean differences of 3.3, −4.7 and 2.5 mm with an RMS of 20.7, 17.5 and 28.0 mm for the X, Y and Z components, respectively for the Greenbelt–Kitab configuration and 1.1, −5.4 and 4.4 mm with an RMS of 9.7, 13.3 and 24.9 mm for the GA_IDS-2_Datum configuration. The larger variability is reflected in the respective amplitudes. As a comparison, the annual amplitudes of the SLR determined dynamic geocentre are 0.9, 1.0 and 6.8 mm in the X, Y and Z components. The results from this study indicate that there is potential to achieve precise dynamically determined geocentre from DORIS.     

Dynamical evolution of high area-to-mass ratio debris released into GPS orbits

A large set of simulations, including all the relevant perturbations, was carried out to investigate the long-term dynamical evolution of fictitious high area-to-mass ratio (A/M) objects released, with a negligible velocity variation, in each of the six orbital planes used by Global Positioning System (GPS) satellites. As with similar objects discovered in near synchronous trajectories, long lifetime orbits, with mean motions of about 2 revolutions per day, were found possible for debris characterized by extremely high area-to-mass ratios. Often the lifetime exceeds 100 years up to A/M ∼ 45 m²/kg, decreasing rapidly to a few months above such a threshold. However, the details of the evolution, which are conditioned by the complex interplay of solar radiation pressure and geopotential plus luni-solar resonances, depend on the initial conditions. Different behaviors are thus possible. In any case, objects like those discovered in synchronous orbits, with A/M as high as 20–40 m²/kg, could also survive in this orbital regime, with semi-major axes close to the semi-synchronous values, with maximum eccentricities between 0.3 and 0.7, and with significant orbit pole precessions (faster and wider for increasing values of A/M), leading to inclinations between 30° and more than 90°.     

DPOD2005: An extension of ITRF2005 for Precise Orbit Determination

For Precise Orbit Determination of altimetry missions, we have computed a data set of DORIS station coordinates defined for specific time intervals called DPOD2005. This terrestrial reference set is an extension of ITRF2005. However, it includes all new DORIS stations and is more reliable, as we disregard stations with large velocity formal errors as they could contaminate POD computations in the near future. About 1/4 of the station coordinates need to be defined as they do not appear in the original ITRF2005 realization. These results were verified with available DORIS and GPS results, as the integrity of DPOD2005 is almost as critical as its accuracy. Besides station coordinates and velocities, we also provide additional information such as periods for which DORIS data should be disregarded for specific DORIS stations, and epochs of coordinate and velocity discontinuities (related to either geophysical events, equipment problem or human intervention). The DPOD model was tested for orbit determination for TOPEX/Poseidon (T/P), Jason-1 and Jason-2. Test results show DPOD2005 offers improvement over the original ITRF2005, improvement that rapidly and significantly increases after 2005. Improvement is also significant for the early T/P cycles indicating improved station velocities in the DPOD2005 model and a more complete station set. Following 2005 the radial accuracy and centering of the ITRF2005-original orbits rapidly degrades due to station loss.     

空间目标轨道确定专家系统

空间目标轨道确定专家系统的面向空间监测需求开发的大型交互式应用软件。它的功能面向空间监测和信息分析中广泛的业务需求，基础是轨道计算软件，数据库支持和一些辅助支持软件。本文详细地介绍了空间目标轨道确定专家系统软件的设计方案。     

基于连续同伦算法的天基仅测角初定轨方法

要进一步提高天基短弧初定轨的精度,在观测资料精度较高的情况下,仅考虑二体问题是不够的,还应考虑轨道摄动的影响。因此,基于无摄初轨的单位矢量法原理和矢量斜分解方法,给出了考虑摄动的天基仅测角初定轨单位矢量法。针对天基仅测角观测条件方程组求解过程中易出现迭代不收敛或收敛到平凡解的问题,引入连续同伦算法求解观测条件方程组,提出了单星观测方式下的空间目标天基仅测角初定轨方法,并通过数值仿真算例验证了该算法在较大范围的收敛性和数值稳定性。     

锚固站数量及布局对导航星座自主定轨影响分析

针对导航星座自主定轨中的星座整体旋转问题,采用增设少量地面锚固站的方法可有效解决该问题。通过推导星地距离对卫星轨道升交点赤经的偏导数,证明了星地距离对卫星轨道升交点赤经可观。仅考虑在我国大陆范围内布设锚固站的条件下,仿真分析了锚固站数量以及布局对导航星座自主定轨精度的影响。仿真实验结果表明:采用3个以上的锚固站,即可有效控制星座整体旋转,在14d的仿真时段内卫星自主定轨精度保持4m以内;锚固站数量越多,自主定轨精度越高,但随着锚固站数量的增加,自主定轨精度改善程度越来越小;在保持4个锚固站的情形下,采用不同的锚固站布局方案,自主定轨精度并无明显差别。     

基于粒子群算法的极短弧定轨

针对经典的初轨计算方法在极短弧定轨中不适用的情况,建立了一种基于粒子群算法的极短弧(TooShort-Arc,TSA)定轨的计算方法。该方法将问题转化为两个三变量的分层优化问题,采用(a,e,M)作为优选变量,在保持问题维数较低的同时,实现了计算结果和观测资料的解耦。由于实测资料处理中的野值剔除方法不适用于粒子群算法,所以,采用稳健估计法,通过在适值函数中使用最小中值二乘准则,实现了稳健的极短弧计算方法。同时,应用MATLAB计算软件,选用缺省参数实现该算法,以进行数据验证。基于实测数据的数值验证表明,方法对于近圆轨道目标30s以下的弧段仍可以获得有效的结果,10s弧段误差仅为16km。此精度满足后续处理的需要,且方法稳健,具有很高的崩溃点。     

简单二体引力模型用于空间目标轨道面交线确定的效果分析

在空间目标碰撞预警分析中，准确地计算出空间目标的轨道面交线是进行地心距筛选、时间差筛选的前提，目前较多使用的快速确定轨道面交线的方法为简单二体引力模型。深入分析该模型，比较了其计算的空间目标轨道面交线与STK计算结果的差异，指出了简单二体引力模型在计算空间目标轨道面交线时的局限性，认为轨道摄动是影响轨道面交线计算准确性的主要原因，应该采用更能反映空间目标实际运动规律的改进的二体引力模型的方法。