基于学习的中段反导拦截时间和拦截点预测方法

弹道导弹实时、准确地预测拦截弹的拦截点与拦截时间，是实现中段突防的有效手段。针对弹道导弹中段突防中的拦截点坐标及拦截时间的预测问题，提出了一种基于监督学习的在线预测方法。以拦截弹的主动段关机参数和关机时刻为输入量，建立拦截时间和拦截点预测模型。在多层感知机神经网络的基础上构建有监督学习算法，通过攻防仿真获取拦截弹的参数制作训练数据集，在线下完成网络训练。仿真结果表明:神经网络能够有效在线预测拦截时间和拦截点坐标，预测结果的相对误差分别为0.124 3%和0.128 5%，拦截时间预测结果误差的平均值为0.224 0 s，拦截点预测结果距离误差平均值为2 016.48 m，均满足精度要求。      

拦截弹道快速设计方法

针对传统优化算法在解决多约束条件下拦截弹道设计问题时迭代计算耗时过长的缺点,提出了一种快速拦截弹道设计方法.该方法基于人工神经网络的非线性拟合技术,先通过普通优化算法进行离线样本计算,对样本中目标位置和优化得到的飞行程序参数进行输入-输出映射拟合训练,从而实现在给定目标位置的情况下对最优飞行程序参数的快速计算,达到提高拦截弹道设计效率的目的.经过数值仿真验证,将该方法应用于拦截弹道设计,在保证良好命中精度的同时,设计所需的计算时间大大减少,从而提高了弹道设计的时效性.      

脉冲推力轨道拦截可达性描述及求解方法

针对航天器单脉冲轨道拦截可达性分析问题,基于共面变轨、逆轨拦截假设,考虑能量、时间和交会角约束,提出了拦截航天器可拦截区和可发射区的概念。在航天器单脉冲空间可达范围的基础上,进一步考虑了目标轨道的约束,建立了目标轨道命中区的计算策略,对异面轨道交叉点为燃料最省点做出了解释。把拦截可达性相关的问题归纳为8个基本拦截问题,通过这8个问题的组合,描述了考虑能量、时间和交会角约束下拦截问题的可拦截区和可发射区的计算方法。采用圆轨道共面变轨、逆轨拦截场景进行了仿真验证,结果表明该方法能够快速有效地计算出约束条件下航天器的拦截可达范围,能够用于分析特定任务情况下的拦截可达性。     

高精度快速趋近滑模变结构末端导引方法

针对反导拦截过程中拦截时间短时传统比例导引制导脱靶量大的问题,运用滑动模态变结构控制理论,提出了快速趋近滑模制导律及高精度制导信息快速提取算法,在典型交战条件下与比例导引、扩展比例导引及开关偏置比例导引方法进行了仿真对比,结果表明该制导律实现了滑模趋近的快速性,有效抑制了指令的发散,滤波算法能快速有效提取有色噪声条件下制导信息,大大提高了制导精度,为高精度反导拦截打下了基础.     

基于最优化问题的混合再入制导方法

在深入研究机动再入飞行器标准轨道和预测制导方法的基础上,针对标准轨道和预测制导相结合的制导方法,分析了由预测制导指令延迟引起的误差,提出了误差补偿策略.根据最优化理论,设计了相关性能指标函数,提出新的预测制导算法.并结合标准轨道法,通过对有限几个特征点位置上引入预测制导,形成了一种基于最优化问题的混合制导方法.有效减小了预测制导指令的解算时间,进一步降低了落点误差.仿真结果表明,该方法既具有较强的抗干扰能力和较高的制导精度,又具有一定的工程实用价值.     

高精度在轨实时轨道机动决策

为保证在轨机动实时性和高精度的要求,提出了一种基于机器学习的在轨实时机动决策方法。通过优化算法离线获得摄动下的精确解,减去二体解得到速度增量差,将其投影到轨道坐标系获得速度增量摄动修正项,以此作为神经网络输出,设计网络参数并训练得到摄动修正网络、组合应用摄动修正网络和二体解实现高精度的在轨实时轨道机动决策。仿真结果表明:卫星按照该决策机动完成后的终端位置偏差与按照优化算法给出的决策机动完成后终端位置偏差精度一致,且前者决策耗时仅为后者决策耗时的0.01%左右。所提轨道机动决策方法兼顾了精度与实时性,适用于星上决策。      

基于预测碰撞点的剩余飞行时间估计方法

分别针对顺轨与逆轨拦截飞行轨迹的特点,设计了相应的剩余飞行时间（TGO）估计方法。该方法通过对弹目碰撞点的预测,降低了发射条件差异对TGO估计精度的影响。首先对线性比例导引运动方程进行变形,得到拦截弹飞行轨迹关于弹目距离的一阶微分方程,基于预测碰撞点,对不同初始发射角造成的积分结果误差进行修正,得出了2种拦截模式下的TGO解析表达式。通过与3种现有估计方法对比分析,验证了提出方法的实时性和估计精度,且能够优化制导性能。      

基于T-H方程的多脉冲最优交会方法

针对椭圆参考轨道交会问题,采用T-H方程描述相对运动,提出了一种时间固定燃料最省的多脉冲最优交会方法,优化参数为交会脉冲及其施加时刻.当考虑到J2摄动或航天器初始相对距离较大时,用T-H方程进行状态预测其线性化误差一般不容忽略,而若用轨道积分预测则耗时较多,进而导致优化时间过长.针对此问题,提出一种采用前一优化脉冲节点的状态导出的轨道根数预测当前节点状态的预测方法.该方法简单实用,有效地加快了优化收敛速度.最后基于多脉冲优化解进行了数值轨道积分以验证交会精度.仿真结果表明,即使加上J2摄动,在初始相对距离为1 000 km时,该方法的终端位置精度仍能达到75 m.     

椭圆轨道航天器相对导航的多项式插值滤波算法

针对运行于椭圆轨道的航天器的相对自主导航问题,建立了惯性坐标系下的近距离相对运动方程,方程中没有圆轨道的假设,相对Clohessy-Wiltshire(C-W)方程具有更广的应用范围。将二阶多项式插值(Divided Difference 2,DD2)滤波应用于相对导航算法中,和扩展卡尔曼滤波(Extended Kalman Filter,EKF)相比,前者不需要对状态方程组进行微分处理,且对非线性函数统计特性近似到二阶程度。数值仿真比较了DD2滤波和EKF两种算法的性能,验证了利用惯性坐标系下的近距离相对运动方程实现相对导航的可行性,并分析了主动航天器的轨道和姿态确定误差对相对导航精度的影响。     

基于单站CCD漂移扫描光电设备的同步轨道目标测定轨

地基光电观测在同步轨道目标监测领域具有重要作用.为评估单站光电设备对同步轨道目标的实际测定轨能力,利用上海天文台佘山站1.56m望远镜,采用CCD漂移扫描光电技术,对3颗北斗同步卫星开展试验观测,基于卫星精密星历评估目标的测定轨外符精度.结果表明:同步轨道目标的天文定位在方位和俯仰方向上的外符精度均好于0.3";在单圈次观测情况下,尽管轨道预报精度较低,约为数千米量级,但是观测弧段内定轨精度可优于百米;在多圈次观测情况下,轨道改进效果显著,定轨精度优于50m,外推至4d的轨道预报精度为百米量级.此外,定量评估了每晚不同观测时间跨度下同步轨道目标的测定轨精度,为单站光电设备实际应用提供了参考.      

“嫦娥4号”中继星任务轨道确定问题初探

"嫦娥4号"任务将采用着陆器、巡视器和绕飞地月拉格朗日L2点中继星进行月球背面的探测,中继星已先期发射,进入环绕地月L2点的晕(Halo)轨道。在中继星使命轨道动力学模型的基础上,通过相关仿真工作,开展了中继星在Halo轨道上的摄动源量级及影响定轨预报的主要因素分析,结果表明:太阳光压摄动是其主要影响因素。为降低其相关影响,提高定轨精度,在太阳光压球模型的基础上,结合中继星在轨运行特点及其星体结构特点,提出了一种求解光压等效面积的方法。经仿真分析,使用修正后的太阳光压球模型进行定轨求解,速度精度可提升约一个量级。     

基于多项式近似的空间碎片群体轨道预报算法

快速准确地分析空间碎片群轨道演化行为对于其他在轨航天器碰撞规避至关重要。在各摄动力的作用下,空间碎片群演化运动呈现出复杂的非线性特征。空间碎片群体个体数量巨大,如果通过对空间碎片群中每个空间碎片进行轨道积分来分析群体预报的方法会导致计算量过大。针对该问题,提出一种基于多项式近似的轨道快速预报分析方法。该方法将空间碎片群分为少量的标称碎片和其他大量关联碎片。针对标称碎片的轨道预报采用数值积分求解保证预报精度;而针对其他大量的关联碎片轨道预报问题,采用多项式泰勒展开半解析方法求解,从而在保证预报精度的前提下有效减少空间碎片群轨道预报的计算量。为了验证方法的有效性,对不同空间碎片群进行了轨道预报仿真。仿真结果表明,当轨道预报精度设定在1m范围内时,多项式近似算法的计算量较蒙特卡洛方法计算效率提高了2.2~17.2倍,验证了所提出方法的有效性。     

大气阻力参数修正对低轨空间目标轨道预报精度的改进

对于低轨空间目标, 大气阻力是影响轨道预报精度的主要摄动力. 本文提出了一种 基于空间环境数据和神经网络模型的空间目标大气阻力参数修正方法, 基于目 标的历史两行元根数, 通过模拟得到外推一天轨道预报中预报结果与观测数据 符合最好的阻力调制系数, 分析表明其与太阳F_10.7指数和地磁Ap指数具有很好的相关性. 根据已有数据, 构建神经网络模型, 实现对阻力调制系数 的补偿计算, 从而改进低轨目标外推一天的轨道预报. 结果表明, 神经网络模 型相比两行元根数能够更及时地对空间环境变化进行响应. 将该方案应用于天 宫一号和国际空间站的外推一天轨道预报, 验证了方案的正确性和普适性, 对 地磁扰动引起的较大预报误差改进效果更好, 误差能够降低50%~60%; 平均而言, 预报精度可以提高约30%, 改进成功率达到80%左右.      

基于经验加速度的低轨卫星轨道预报新方法

研究将定轨过程中的经验加速度应用于地球低轨卫星轨道预报的新方法. 利用GPS伪距观测数据和简化动力学最小二乘批处理方法对地球低轨卫星定 轨, 其中卫星位置、速度及大气阻力系数和辐射光压系数可以直接用于轨道预报. 作为简化动力学最重要特征的经验加速度呈现准周期、余弦曲线特点, 可通过 傅里叶级数拟合建模. 确定性动力学模型与补偿大气阻力模型误差的切向经验 加速度级数拟合模型组成增强型动力学模型用于提高轨道预报精度. 应用 GRACE-A星载GPS伪距观测数据和IGS超快星历定轨并进行轨道预报, 结果表明 轨道预报初值位置精度达到0.2m, 速度精度达到1×10-4m·s-1, 预报3天位置精度优于60m, 比只利用确定性动力学模型进行预报精度平 均提高2.3倍. 先定轨后预报的模式可用在星上自主精确导航系统中.      

基于ARVM模型的液体火箭发动机试验台故障预测方法

针对液体火箭发动机试验台故障子样少,故障预测精度低,故障维修保障困难等问题,在分析标准RVM优缺点的基础之上,提出了一种自适应能力较强的故障预测模型——ARVM(Adaptive Relevance Vector Machine)。为测试该模型,以某型轨控发动机高空模拟试验台管路流量、燃烧室压力为输入参量对推力矢量进行了预测,预测结果表明,ARVM方法能够有效跟踪推力矢量参数的变化趋势,并且获得了较高的预测精度和模型稀疏性。该方法对于复杂系统的故障预测和维修保障具有一定的理论价值和工程应用意义。     

ERS-1 and ERS-2 operational and precise orbit determination

This paper presents the European Space Operations Centre's orbit determination and prediction systems for the ERS-1 mission. The routine operational orbit determination and prediction subsystem is discussed briefly, and statistics of the accuracy compared to the requirements are given. The precise orbit determination subsystem is then described, and the accuracy of its results are compared to those of the operational orbit system and to the D-PAF preliminary orbit solutions. Some geophysical results from the altimeter data, processed in these orbit determinations, are also presented. The ESOC/OAD ‘ERS-1 Orbit Report’ is introduced as a document providing this information on a monthly basis. Finally, this paper describes how the experience gained with the precise orbit determination will be exploited to further improve the accuracy of the routine system that will be used for ERS-2, and provides an estimate of this accuracy.     

星载观测平台的视线误差分析

为分析星载观测平台各项误差源对视线测量精度的影响,研究了一种基于多参量的视线测量误差建模与评价方法.不同于以往卫星、相机的空间测量方式,针对卫星、转台、相机的观测结构,构建了从惯性空间到光学传感器像平面的目标成像模型及星载观测平台视线测量模型.通过推导星载观测平台视线测量误差与观测中13项误差源的关系,提出一种基于灵敏度分析的误差评价方法,并在三个轴向上分析了各项误差源对星载观测平台视线测量精度的影响.利用蒙特卡罗仿真试验验证了理论模型的有效性.结果表明,卫星轨道误差、卫星姿态误差、载荷平台角振动误差、内外框架转动误差、像平面目标像点的位置量化误差是影响星载观测平台视线测量精度的主要因素.该方法能够科学评估各项误差源对星载观测平台视线测量精度的影响,对星载观测平台的总体设计具有重要的应用价值.      

长周期高精度回归轨道与脉冲轨道控制策略设计

针对长周期回归轨道设计和维持问题，研究了一种基于高阶Poincaré映射的高精度引力场中回归轨道优化设计和控制的通用性半解析方法，其中摄动包含大气阻力、太阳辐射压力和日月三体引力等因素。通过对Poincaré映射进行高阶展开并表示为多项式形式，可精确近似对一个或者多个回归周期内的轨道递推，从而在赤道升交点处施加脉冲推力，实现高精度的回归轨道设计和控制。提出了分别解决严格和宽松2种回归约束下问题的方法，并应用于实际在轨TerraSAR-X、Landsat-8、IRS-P6、SPOT-7和UoSAT-12任务的回归模式。所提方法具有计算效率和精度高的优点，可用于星上自主轨道递推和轨道控制。      

Orbit determination of BeiDou navigation satellite system maneuvered satellites based on instantaneous velocity pulses parameters

The BeiDou navigation satellite system (BDS) comprises geostationary earth orbit (GEO) satellites as well as inclined geosynchronous orbit (IGSO) and medium earth orbit (MEO) satellites. Owing to their special orbital characteristics, GEO satellites require frequent orbital maneuvers to ensure that they operate in a specific orbital window. The availability of the entire system is affected during the maneuver period because service cannot be provided before the ephemeris is restored. In this study, based on the conventional dynamic orbit determination method for navigation satellites, multiple sets of instantaneous velocity pulses parameters which belong to one of pseudo-stochastic parameters were used to simulate the orbital maneuver process in the orbital maneuver arc and establish the observed and predicted orbits of the maneuvered and non-maneuvered satellites of BeiDou regional navigation satellite system (BDS-2) and BeiDou global navigation satellite system (BDS-3). Finally, the single point positioning (SPP) technology was used to verify the accuracy of the observed and predicted orbits. The orbit determination accuracy of maneuvered satellites can be greatly improved by using the orbit determination method proposed in this paper. The overlapping orbit determination accuracy of maneuvered GEO satellites of BDS-2 and BDS-3 can improve 2–3 orders of magnitude. Among them, the radial orbit determination accuracy of each maneuvered satellite is basically better than 1 m. simultaneously, the combined orbit determination of the maneuvered and non-maneuvered satellites does not have a great impact on the orbit determination accuracy of the non-maneuvered satellites. Compared with the multi GNSS products (indicated by GBM) from the German Research Centre for Geosciences (GFZ), the impact of adding the maneuvered satellites on the orbit determination accuracy of BDS-2 satellites is less than 9 %. Furthermore, the orbital recovery time and the service availability period are significantly improved. When the node of the predicted orbit is traversed approximately 3 h after the maneuver, the accuracy of the predicted orbit of the maneuvered satellite can reach that of the observed orbit. The SPP results for the BDS reached a normal level when the node of the predicted orbit was 2 h after the maneuver.