[发明专利]延长地球静止轨道通信卫星两行星历预报时段方法

权利要求书

1.一种延长地球静止轨道通信卫星两行星历预报时段方法，其特征在于，包括：

获取地球静止轨道通信卫星的两行轨道数据，其中，所述两行轨道数据包括平均运动圈数及开普勒轨道参数；

提取所述平均运动圈数和所述开普勒轨道参数；

修正所述平均运动圈数，根据修正的平均运动圈数及所述开普勒轨道参数以修正所述地球静止轨道通信卫星的轨道参数；

利用修正的地球静止轨道通信卫星的轨道参数计算所述地球静止轨道通信卫星的位置。

2.根据权利要求1所述的方法，其特征在于，所述提取所述平均运动圈数和开普勒轨道参数，包括：

轨道周期T满足：

T＝86400/N

其中，N为所述平均运动圈数。

3.根据权利要求1所述的方法，其特征在于，所述修正所述平均运动圈数，包括：修正所述平均运动圈数为1.00273896。

4.根据权利要求1所述的方法，其特征在于，所述根据修正的平均运动圈数及所述开普勒轨道参数以修正所述地球静止轨道通信卫星的轨道参数，包括：根据修正的平均运动圈数以修正轨道半长轴a、平近点角M、偏近点角E和真近点角f。

5.根据权利要求4所述的方法，其特征在于，所述轨道半长轴a、平近点角M、偏近点角E和真近点角f满足如下关系：

在时刻t的轨道半长轴a满足：

在时刻t的平近点角M满足：

以及，所述偏近点角E和所述平近点角M满足：

E＝M+e sin E，

所述真近点角f和所述偏近点角E满足：

其中，T＝86400/N，N为所述平均运动圈数，μ为万有引力常数μ＝3.986005×10¹⁴m³·s^-2；M₀为轨道历元时刻的平近点角，t_0e为轨道历元时刻，e为轨道偏心率，其中，所述轨道历元时刻的平近点角M₀、轨道历元时刻t_0e、轨道偏心率e从所述两行轨道数据中直接提取。

6.根据权利要求1所述的方法，其特征在于，所述利用修正的地球静止轨道通信卫星的轨道参数计算所述地球静止轨道通信卫星的位置，包括：

在时刻t满足时，所述地球静止轨道通信卫星在地球坐标系下的三维坐标为：

其中，为所述地球静止轨道通信卫星在地球坐标系下的三维坐标，为地心惯性坐标系下的坐标，GAST为t时刻的格林尼治恒星时，R₃(GAST)为三维坐标旋转矩阵，为转换矩阵，r为所述地球静止轨道通信卫星到地心的距离，f为真近点角。

7.根据权利要求6所述的方法，其特征在于，所述距离r满足：

r＝a(1-e cos E)，

其中，a为轨道半长轴，在时刻t的轨道半长轴a满足：

E为偏近点角，所述偏近点角E和平近点角M满足：

E＝M+e sin E，

在时刻t的平近点角M满足：

所述真近点角f和所述偏近点角E满足：

其中，T＝86400/N，N为所述平均运动圈数，μ为万有引力常数μ＝3.986005×10¹⁴m³·s^-2；M₀为轨道历元时刻的平近点角，t_0e为轨道历元时刻，e为轨道偏心率，其中，所述轨道历元时刻的平近点角M₀、轨道历元时刻t_0e、轨道偏心率e从所述两行轨道数据中直接提取。