[发明专利]一种地球静止卫星非球形引力切向漂移加速度计算方法

权利要求书

1.一种地球静止卫星非球形引力切向漂移加速度计算方法，其特征在于，包括以下步骤：

步骤1、建立地球非球形引力半长轴摄动方程，得到半长轴变化率；建立的地球非球形引力半长轴摄动方程为：

式(1)中，J₂₂、J₃₁、J₃₃、J₄₂、J₄₄、l₂₂、l₃₁、l₃₃、l₄₂、l₄₄为地球引力场4阶摄动系数，r_e为地球半径，a为轨道半长轴，l＝ω+M+Ω-Θ，l为平经度，ω为近地点幅角，M为平近点角，Ω为升交点赤经，Θ为格林尼治恒星时角，μ为地球引力常数；

步骤2、建立漂移加速度计算方程；具体包括：

步骤2.1、设定漂移率公式：

式(2)中，ω_e为地球自转角速度，a_c为地球静止卫星标称半长轴；

步骤2.2、式(2)对时间t求导，得到漂移加速度计算公式：

步骤3、将步骤1中得到的半长轴变化率代入步骤2中建立的漂移加速度计算方程，进而通过代入地球静止卫星当前平经度解析漂移加速度。

2.如权利要求1所述的一种地球静止卫星非球形引力切向漂移加速度计算方法，其特征在于，所述式(3)中轨道半长轴a设定为标称值，即a＝a_c。

3.如权利要求2所述的一种地球静止卫星非球形引力切向漂移加速度计算方法，其特征在于，所述步骤3具体包括将式(1)代入式(3)得到漂移加速度计算方程为：

式(4)中，代入地球静止卫星当前平经度l，解析计算出漂移加速度D_acc。