[发明专利]一种航天器高精度自主目标预报方法

说明书

本发明公开了一种航天器自主高精度目标预报方法，首先对目标航天器各长期项轨道摄动进行建模和理论分析，并结合星载计算机的处理能力，对摄动力进行简化建模。由初始的平均根数，结合摄动模型，求取目标当前时刻的平根；建立短周期项摄动方法，由当前的平根，结合短周期摄动，求其当前的瞬时根数。不同于常规的轨道根数递推，不能忽略偏心率的影响。该方法有效实现了目标在轨高精度预报，为实现航天器自主控制创造条件。

技术领域

本发明涉及一种航天器自主目标预报技术，具体是一种基于地面上注的初始轨道根数，采用星上自主轨道递推算法，对目标航天器轨道进行预报的方法。

背景技术

随着空间交会对接试验和对非合作目标空间机动和接近观测试验的增多，对目标航天器高精度导航的需求增加，如何在无跟瞄单机的情况下实现对目标航天器高精度预报成为一项关键技术。

对于常规的航天器轨道理论，轨道摄动建模比较复杂，星载计算机处理能力满足复杂运算。本发明对对轨道摄动进行分析和简化，即满足导航精度要求又适应计算机处理要求。

发明内容

本发明针对现有技术中存在的上述不足，提供了一种航天器高精度自主目标预报方法。

本发明是通过以下技术方案实现的。一种航天器高精度自主目标预报方法，包括如下步骤：

步骤1，对一阶长期项轨道摄动进行建模，由初始地平均根数，求取目标航天器当前时刻平均轨道根数：

步骤2，由当前时刻平均根数，计算地球非球形摄动一阶短周期摄动项影响下轨道要素的变化，进而计算当前时刻的瞬时轨道根数；

步骤3，由当前时刻的瞬时轨道根数结合根数到位置速度转化关系，计算目标航天器位置和速度。

优选的，所述步骤1具体包括如下步骤：

(1.1)对于低轨航天器，考虑主要误差源及一阶长期项，确定轨道要素在轨道圈内的平均变化率所述的主要误差源为非球形引力摄动；

考虑大气阻力影响，主要改变航天器大小，确定轨道要素在轨道圈内的平均变化率；

(1.2)结合(1.1)中的结果，确定在地球非球形引力摄动及大气阻力共同作用下，航天器轨道根数平均变化率；

(1.3)由航天器初始的平均根数结合(1.2)中的轨道根数平均变化率，得到当前t时刻平均轨道根数。

优选的，考虑主要误差源及一阶长期项，轨道要素在轨道圈内的平均变化率为：

其中，a为航天器半长轴、e为航天器偏心率、i为航天器轨道倾角、Ω为航天器升交点赤经、ω为航天器近地点幅角、M为航天器平近点角、n为航天器轨道角速度、R_e为平均赤道半径；J₂为地球引力场二阶带谐系数。

优选的，考虑大气阻力影响下轨道要素在轨道圈内的平均变化率为：

其中，da为半长轴变化率。

优选的，所述的半长轴变化率da通过统计多轨航天器半长轴变化量计算得到。