火星探测直接转移轨道精确求解研究

针对火星探测直接转移轨道精确求解问题, 提出了一种快速微分修正算法. 基于二体模型建立控制参数和目标参数偏差关系的数学模型, 并由此求解系统的偏导数矩阵; 采用B平面参数作为中间变量, 对算法进行两层迭代设计, 有效地减少了求解过程中的积分运算次数. 以2018年火星探测机会为例对算法进行了验证, 仿真结果表明, 利用圆锥曲线拼接法得到的轨道初值, 求解一条标准轨道只需6~9次积分迭代. 通过STK对计算结果进行了对比和验证.      

一种图像辅助火星着陆段自主导航方法

针对未来深空探测软着陆高精度实时导航的需求,提出了一种图像辅助的自主导航方案。首先通过下降段图像与落点区域地形匹配,获得着陆器相对于落点的位置和姿态;然后基于误差Kalman模型估计着陆器的状态误差,修正惯性导航的结果;在图像信息不可用的情况下,只进行惯性递推导航;这种方案既提高了导航的精度,也能保证实时性的要求。 数值仿真验证了该方案的可行性,该方案对未来实际工程任务具有很大的参考价值。     

卫星星座整网轨道确定分析与仿真

编队飞行的卫星或卫星星座对轨道确定自主性和精度提出了较高要求。针对这个问题,通过建立星座的轨道动力学模型和星间观测的测量模型,将星座中的星间观测数据和地面观测数据融合起来,将待估的卫星轨道参数和部分动力学参数进行适当的分类,研究卫星星座整网轨道确定的新方法,并在理论上分析了整网定轨方法能提高定轨精度的原因;最后采用自主开发的卫星星座整网轨道确定软件进行了仿真计算。计算表明,该方法能有效地减少对地面站的依赖,并较大幅度提高定轨时卫星绝对位置和相对位置精度。
     

月球卫星精密定轨

随着空间应用需求的日益增大,深空探测已成为现实,而月球显然是人类走向深空的首选目标。发射月球探测器通常分3个阶段,其运动状态分别对应3种不同类型的轨道:近地停泊轨道、地月转移轨道和绕月轨道。月球是1个慢自转天体且无大气,就轨道解而言这些因素导致环月卫星的运动与地球卫星有所差别。本文针对月球探测任务的特点,从月球与地球的差别入手,在仔细分析月球卫星的受力状况前提下,着重阐述月球探测器在环月段精密定轨的方法原理和具体实现过程。     

基于多站数据融合的参数精估计方法

针对侦察设备处于星载SAR副瓣照射范围,从而导致截获信号湮没于强噪声背景这个问题,本文提出一种基于多站接收机之间的数据融合方法。在信号形式未知的情况下,通过此方法可以检测出淹没在噪声中的微弱信号,进行信号的分类和时频域参数的精估计。首先,将参考接收机与其他接收机之间进行互相关处理,得到峰值信息,根据峰值信息的位置得到信号与参考信号之间的延迟位置,进行延迟校准;其次,各个接收机分别进行粗步长的分数阶傅里叶变换（Fractional Fourier Transform,FrFT）,记录峰值信息为精估计做准备,根据峰值角度和分数阶傅里叶反变换恢复出原始信号;最后,判定是否存在信号,若信号存在实现多站原始信号功率比的加性融合,根据多站峰值信息限定旋转角度范围,采用精步长的分数阶傅里叶变换估计出调频率和中心频率;利用联合互相关谱实现信号能量的累积,采用自适应门线和边界波谷连续取小方法,找到信号存续状态中的左右边界,估计出带宽和中心频率,计算脉宽,实现时频域信号的精估计。仿真实验表明：该方法可以在低信噪比的高斯白噪声和有色噪声背景下,对线性调频信号（Chirp）的时频参数进行有效的精估计。     

Comparison of convergence time and positioning accuracy among BDS,GPS and BDS/GPS precise point positioning with ambiguity resolution

Precise point positioning (PPP) usually takes about 30?min to obtain centimetre-level accuracy, which greatly limits its application. To address the drawbacks of convergence speed and positioning accuracy, we develop a PPP model with integrated GPS and BDS observations. Based on the method, stations with global coverage are selected to estimate the fractional cycle bias (FCB) of GPS and BDS. The short-term and long-term time series of wide-lane (WL) FCB, and the single day change of narrow-lane (NL) FCB are analysed. It is found that the range of GPS and BDS non-GEO (IGSO and MEO) WL FCB is stable at up to a 30-day-time frame. At times frame of up to 60?days, the stability is reduced a lot. Whether for short-term or long-term, the changes in the BDS GEO WL FCB are large. Moreover, BDS FCB sometimes undergoes a sudden jump. Besides, 17 and 10 stations were used respectively to investigate the convergence speed and positioning errors with six strategies: BDS ambiguity-float PPP (Bfloat), GPS ambiguity-float PPP (Gfloat), BDS/GPS ambiguity-float PPP (BGfloat), BDS ambiguity-fixed PPP (Bfix), GPS ambiguity-fixed (Gfix), and BDS/GPS ambiguity-fixed (BGfix). The average convergence speed of the ambiguity-fixed solution is greatly improved compared with the ambiguity-float solution. In terms of the average convergence time, the Bfloat is the longest and the BGfix is the shortest among these six strategies. Whether for ambiguity-float PPP or ambiguity-fixed PPP, the convergence reduction time in three directions for the combined system is the largest compared with the single BDS. The average RMS value of the Bfix in three directions (easting (E), northing (N), and up (U)) are 2.0?cm, 1.5?cm, and 5.9?cm respectively, while those of the Gfix are 0.8?cm, 0.5?cm, and 1.7?cm. Compared with single system, the BDS/GPS combined ambiguity-fixed system (BGfix) has the fastest convergence speed and the highest accuracy, with average RMS as 0.7?cm, 0.5?cm, and 1.9?cm for the E, N, U components, respectively.     

Real-time zenith tropospheric delays in support of numerical weather prediction applications

The Geodetic Observatory Pecný (GOP) routinely estimates near real-time zenith total delays (ZTD) from GPS permanent stations for assimilation in numerical weather prediction (NWP) models more than 12 years. Besides European regional, global and GPS and GLONASS solutions, we have recently developed real-time estimates aimed at supporting NWP nowcasting or severe weather event monitoring. While all previous solutions are based on data batch processing in a network mode, the real-time solution exploits real-time global orbits and clocks from the International GNSS Service (IGS) and Precise Point Positioning (PPP) processing strategy. New application G-Nut/Tefnut has been developed and real-time ZTDs have been continuously processed in the nine-month demonstration campaign (February–October, 2013) for selected 36 European and global stations. Resulting ZTDs can be characterized by mean standard deviations of 6–10 mm, but still remaining large biases up to 20 mm due to missing precise models in the software. These results fulfilled threshold requirements for the operational NWP nowcasting (i.e. 30 mm in ZTD). Since remaining ZTD biases can be effectively eliminated using the bias-reduction procedure prior to the assimilation, results are approaching the target requirements in terms of relative accuracy (i.e. 6 mm in ZTD). Real-time strategy and software are under the development and we foresee further improvements in reducing biases and in optimizing the accuracy within required timeliness. The real-time products from the International GNSS Service were found accurate and stable for supporting PPP-based tropospheric estimates for the NWP nowcasting.     

Influence of time variable geopotential models on precise orbits of altimetry satellites,global and regional mean sea level trends

During the last decade a significant progress has been reached in the investigation of the gravity field of the Earth. Besides static, also time variable geopotential models have been recently created. In this paper we investigate the impact of the recent time variable geopotential models on altimetry satellite orbits and such altimetry products based on these orbits, as global and regional mean sea level trends. We show that the modeling of time variable gravity improves the orbit solutions, at least for the GRACE period where time variable gravity is sufficiently accurately observed by this mission. Our analysis includes six geopotential models jointly developed by GFZ German Research Centre for Geosciences and Space Geodesy Research Group (CNES/GRGS) Toulouse: the stationary model EIGEN-GL04S, a stationary version of EIGEN-6S (EIGEN-6S_stat), a corrected version of EIGEN-6S and three enhanced versions of EIGEN-6S called EIGEN-6S2, EIGEN-6S2A and EIGEN-6S2B. By “stationary” we mean “containing periodic parameters such as annual and semi-annual variations, but no secular (drift) terms”. We computed precise orbits for the radar altimetry satellites ERS-1, ERS-2, TOPEX/Poseidon, and Envisat over 20 years between 1991 and 2011. The orbit, single-mission and multi-mission altimetry crossover analyses show that the time variable models EIGEN-6S_corrected, EIGEN-6S2 and its two precursors EIGEN-6S2A/B perform notably better than the stationary models for the GRACE period from 2003 onwards. Thus, using EIGEN-6S2 and EIGEN-6S2A/B we have got 3.6% smaller root mean square fits of satellite laser ranging observations for Envisat, as when using EIGEN-GL04S. However, for the pre-GRACE period 1991–2003, the stationary geopotential models EIGEN-GL04S and EIGEN-6S_stat as well as EIGEN-6S2 having no drift terms for degree 3–50 at this time interval perform superior compared to EIGEN-6S_correct and EIGEN-6S2A/B which contain drifts for this period. We found, that the time variable geopotential models have a low (0.1–0.2 mm/yr) impact on our results for the global mean sea level trend. However, we found strong East/West differences up to 3 mm/yr in the regional mean sea level trends when using orbits of all four satellites based on time variable and stationary geopotential models. We show that these differences are related to the relative drifts of the centers-of-origin between the orbit solutions based on the time variable and stationary geopotential models. From the results of our detailed study, we conclude that the final version of the time variable gravity field model EIGEN-6S2 performs best for the four satellites tested. This model provides the most reliable and mission-consistent sea level estimates for the whole time period from 1992 to 2010. This model is of maximum spherical harmonic degree and order 260 and contains time series for drifts as well as annual and semiannual variations of the spherical harmonic coefficients for degree 2–50.     

频率捷变时间精密测试系统

较详细地介绍了捷变时间的定义以及６种常用测量方法的原理、数据处理过程、及测量结果等内容。６种方法分别为检相－示波器法、鉴频－示波器法、检相－时间间隔计数器法、鉴频一时间间隔计数器法，ＨＰ５３７２Ａ测频和ＨＰ５３７２Ａ测相。给出了这几种测量方法的适用范围和优缺点。     

高精度连续时频信号产生和控制系统的分析

具体分析了NTSC和LORAN-C的两套时频基准系统。结果表明,NTSC时频基准系统可以产生高精度的时间频率信号,但是没有主标备份,不利于连续性的保持;而LORAN-C的时频基准系统产生控制高精度时频信号,同时主标发生故障时自动切换到副标,不停止信号的输出,而且切换主标时不会发生信号瞬间中断或相位突变,这些又保证了信号的瞬间连续性。它可广泛应用于守时、授时及导航系统中。