基于卫星观测数据的弹道导弹初始状态估计

利用卫星探测数据对大气层外飞行的弹道导弹弹道进行估计。在极大似然估计中,由于卫星探测导弹运动的可观测性较差,使初始点估计变得十分病态,传统的Gauss Newton迭代算法不再适用,而采用较为复杂Levenberg Marquardt算法,分别对单颗卫星和两颗卫星的情况进行了Monte Carlo仿真,得到了各自的归一化残差测量误差(代价函数),结果表明该算法有效。     

基于深度增强学习的卫星姿态控制方法

针对卫星在执行丢弃载荷或捕获目标等复杂任务时遭遇的姿态突然发生变化的问题,采用深度增强学习方法对卫星姿态进行控制,使卫星恢复稳定状态。具体来说,首先搭建飞行器的姿态动力学环境,并将连续的控制力矩输出离散化,然后采用Deep Q Network算法进行卫星自主姿态控制训练,以姿态角速度趋于稳定作为奖励获得离散行为的最优智能输出。仿真试验表明,面向空间卫星姿态控制的深度增强学习算法能够在卫星受到突发随机扰动后稳定卫星姿态,并能有效解决传统PD控制器依赖被控对象质量参数的难题。所提出的方法采用自主学习的方式对卫星姿态进行控制,具有很强的智能性和一定的普适性,在未来卫星执行复杂空间任务中的智能控制方面有着很好的应用潜力。     

一种应用于实时深空干涉测量的电离层时延修正方法

针对"嫦娥4号"中继星任务S频段信标信号的高精度实时干涉测量需求,结合深空测控干涉测量系统采用的稀疏标校工作模式,研究验证了一种面向测控模式实时干涉测量的电离层时延修正方法。首先分析了电磁波经电离层传播的延迟机理及特性;基于深空站历史观测数据,通过自相关函数分析验证了天顶向TEC的周日特性;在此基础上,结合深空干涉测量中心数据处理设备软件系统,讨论了电离层时延修正方法;通过任务期间的实测数据处理分析,验证了所提方法可以将实时测量精度提升1～3 ns,对低仰角跟踪弧段,该技术方法优势更为明显。该方法为后续深入推进深空测控干涉测量系统在任务中的实时应用提供了技术储备。     

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

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

Ionospheric TEC disturbances over China during the strong geomagnetic storm in September 2017

From September 7 to 8, 2017, a G4-level strong geomagnetic storm occurred, which seriously impacted on the Earth’s ionosphere. In this work, the global ionospheric maps released by Chinese Academy of Sciences are used to investigate the ionospheric responses over China and its adjacent regions during the strong storm. The prominent TEC enhancements, which mainly associated with the neutral wind and eastward prompt penetration electric field, are observed at equatorial ionization anomaly crests during the main phase of the storm on 8 September 2017. Compared with those on 8 September, the TEC enhancements move to lower-latitude regions during the recovery phase on 9 September. A moderate storm occurred well before the start of the strong storm causes similar middle-latitude TEC enhancements on 7 September. However, the weak TEC depletion is observed at middle and low latitude on 9–10 September, which could be associated with the prevailing westward disturbance electric field or storm-time neural composition changes. In addition, the storm-time RMS and STD values of the ionospheric TEC grids over China increase significantly due to the major geomagnetic storm. The maximum of the RMS reaches 12.0 TECU, while the maximum of the STD reaches 8.3 TECU at ~04UT on 8 September.     

The propagation features of LF radio waves at topside ionosphere and their variations possibly related to Wenchuan earthquake in 2008

The artificial low frequency (LF) signals recorded at topside ionosphere from DEMETER satellite were analyzed in this paper, and the typical diurnal and seasonal variations were illustrated around 162?kHz in electric field spectra. The larger power spectrum density (PSD) values in electric field at local nighttime and in winter season all demonstrate the correlation feature of lower plasma content to higher penetration of LF waves into ionosphere. Around Wenchuan earthquake, the comparison of signal-noise-ratio (SNR) values in electric field with each half month during January to May in 2008 and the same half month in May from 2005 to 2007 revealed their lowest values and small covering area around the preparation region of Wenchuan earthquake in 2008. Combined with other researches in VLF radio waves and geochemical observations from satellite, the interaction of ion accumulation and upward movement from gas-water release at surface might be a key factor to disturb the ionospheric plasma density, and then possibly leading to the decrease of low energy penetration of LF radio waves from the artificial transmitted source at ground.     

A holistic approach to SIM platform and its application to early-warning satellite system

This study proposes a new simulation platform named Simulation Integrated Management (SIM) for the analysis of parallel and distributed systems. The platform eases the process of designing and testing both applications and architectures. The main characteristics of SIM are flexibility, scalability, and expandability. To improve the efficiency of project development, new models of early-warning satellite system were designed based on the SIM platform. Finally, through a series of experiments, the correctness of SIM platform and the aforementioned early-warning satellite models was validated, and the systematical analyses for the orbital determination precision of the ballistic missile during its entire flight process were presented, as well as the deviation of the launch/landing point. Furthermore, the causes of deviation and prevention methods will be fully explained. The simulation platform and the models will lay the foundations for further validations of autonomy technology in space attack-defense architecture research.     

Equatorial anomaly according to the Interkosmos-19 data and IRI model: A comparison

The comparison of the IRI model with the foF2 distribution in the equatorial anomaly region obtained by topside sounding onboard the Interkosmos-19 satellite has been carried out. The global distribution of foF2 in terms of LT-maps was constructed by averaging Intercosmos-19 data for summer, winter, and equinox. These maps, in fact, represent an empirical model of the equatorial anomaly for high solar activity F10.7 ~ 200. The comparison is carried out for the latitudinal foF2 profiles in the characteristic longitudinal sectors of 30, 90, 210, 270, and 330°, as well as for the longitudinal variations in foF2 over the equator. The largest difference between the models (up to 60%) for any season was found in the Pacific longitudinal sector of 210°, where there are a few ground-based sounding stations. Considerable discrepancies, however, are sometimes observed in the longitudinal sectors, where there are many ground-based stations, for example, in the European or Indian sector. The discrepancies reach their maximum at 00 LT, since a decay of the equatorial anomaly begins before midnight in the IRI model and after midnight according to the Interkosmos-19 data. The discrepancies are also large in the morning at 06 LT, since in the IRI model, the foF2 growth begins long before sunrise. In the longitudinal variations in foF2 over the equator at noon, according to the satellite data, four harmonics are distinguished in the June solstice and at the equinox, and three harmonics in the December solstice, while in the IRI model only two and one harmonics respectively are revealed. In diurnal variations in foF2 and, accordingly, in the equatorial anomaly intensity, the IRI model does not adequately reproduce even the main, evening extremum.     

Orbit determination and prediction of GEO satellite of BeiDou during repositioning maneuver

In order to establish a continuous GEO satellite orbit during repositioning maneuvers, a suitable maneuver force model has been established associated with an optimal orbit determination method and strategy. A continuous increasing acceleration is established by constructing a constant force that is equivalent to the pulse force, with the mass of the satellite decreasing throughout maneuver. This acceleration can be added to other accelerations, such as solar radiation, to obtain the continuous acceleration of the satellite. The orbit determination method and strategy are illuminated, with subsequent assessment of the orbit being determined and predicted accordingly. The orbit of the GEO satellite during repositioning maneuver can be determined and predicted by using C-Band pseudo-range observations of the BeiDou GEO satellite with COSPAR ID 2010-001A in 2011 and 2012. The results indicate that observations before maneuver do affect orbit determination and prediction, and should therefore be selected appropriately. A more precise orbit and prediction can be obtained compared to common short arc methods when observations starting 1 day prior the maneuver and 2 h after the maneuver are adopted in POD (Precise Orbit Determination). The achieved URE (User Range Error) under non-consideration of satellite clock errors is better than 2 m within the first 2 h after maneuver, and less than 3 m for further 2 h of orbit prediction.