[发明专利]一种星上实时计算环月轨道的方法

说明书

技术领域

本发明涉及航天器实时计算环月轨道技术领域。

背景技术

目前星上实时轨道计算一般只在近地轨道卫星上使用，所使用的方法是平根数轨道计算方法。国内航天器上采用注入19根数进行计算，国外采用SGP4模型进行计算，这两种方法类似，都是采用平均轨道根数及其长期项系数、长周期项系数、短周期项系数计算卫星瞬时轨道根数。由于地球引力场的J2项摄动和大气阻力摄动比其他项大的多，在使用平根数计算时只需要考虑地球引力的J2项和大气阻力项，因而公式较为简单，计算量不大。而月球引力场的J2项摄动的量级与其他高阶项接近，如果采用平根数法进行计算，则需要考虑20阶月球引力和地球引力摄动，即使经过简化后计算量仍然非常大，不适合星上计算机进行计算。

在月球卫星环月期间，定姿需要卫星在月球坐标系内实时的轨道根数，对卫星轨道计算精度提出的要求是位置误差小于1km。如果采用传统的平根数法，要做到较高精度的轨道计算，由于月球引力场的复杂性，计算时采用的月球引力模型就需要取到20阶左右的较高阶次，另外在较高精度要求下还需要考虑地球引力摄动的影响，从而导致计算量过大，不适合在星上计算。而数值解的方法尽管精度较高，但也存在计算量较大的问题。

发明内容

本发明的技术解决问题是：克服现有技术的不足之处，提供一种地面运算与星上运算相结合的方式，使得环月卫星能够在星上实施计算出环月轨道。

本发明的方法的技术解决方案是：一种星上实时计算环月轨道的方法，包括以下步骤：

(1)以定轨结果作为初值，用RKF78数值方法计算出实际轨道；

(2)选取一条二体轨道作为参考轨道；

(3)计算实际轨道相对于参考轨道的运动轨迹，采用三角函数来拟合实际轨道相对于参考轨道的运动；

(4)根据公知的两条二体轨道之间的相对运动规律，计算出一条相对于参考轨道运动轨迹与实际轨道相对于参考轨道运动轨迹相接近的二体轨道；

步骤(1)至步骤(4)均在地面完成，将一天分成12个时间段，分别重复利用步骤(1)至步骤(4)计算相应的二体轨道，然后将计算结果注入到星上；

(5)卫星采用二体轨道计算方法，根据步骤(4)中算得的二体轨道实时计算出卫星的轨道根数。

本发明与现有技术相比的有益效果是：本发明通过地上运算与星上运算相结合的方式，特别是选取了一条二体轨道作为参考轨道，使得星上在计算量很小的情况下能够实施计算出卫星轨道，同时不需要地面和卫星保持实时联系，这种轨道计算方法的位置误差小于1km。

附图说明

图1为本发明星上计算轨道的水平面内误差情况；

图2为本发明星上计算轨道的轨道面内误差情况。

具体实施方式

(1)在地面通过高精度的数值运算计算出月球卫星在未来一段时间内的月心惯性系轨道根数，以一定的步长h保存(h≤60s)，轨道计算的初值由定轨所得。因为摄动原因得到的结果随时间而变化，以σ_t表示此根数，并认为此根数表示了卫星的实际位置。计算过程需要考虑的摄动包括月球引力非球形引力摄动，其阶次应取到30阶以上，同时还需考虑地球引力、太阳引力等对月球卫星的摄动。计算方法采用高精度的数值解法。

例如，取卫星定轨的结果为：

起始时间(北京时)t0＝2007年11月7日12:0:0.000

卫星轨道根数为：

半长轴a＝1934km；

偏心率e＝0.00237；

倾角i＝87.41度；

升交点赤经Ω＝265.51度；

近地点幅角ω＝281.585；

平近点角M＝202。

精确轨道计算采用RKF78数值方法，考虑月球引力场为70×70阶、考虑地球引力摄动、太阳引力摄动。

(2)选取一条二体轨道作为参考轨道，使得上述计算出的轨道相对于该二体轨道在小范围内运动。可以选择开始时刻对应的二体轨道作为参考轨道，其轨道根数记为σ₀。参考轨道选择为：