专利名称：：一种星上实时计算环月轨道的方法
技术领域：
：本发明涉及航天器实时计算环月轨道
技术领域：
。技术背景目前星上实时轨道计算一般只在近地轨道卫星上使用，所使用的方法是平根数轨道计算方法。国内航天器上采用注入19根数进行计算，国外采用SGP4模型进行计算，这两种方法类似，都是采用平均轨道根数及其长期项系数、长周期项系数、短周期项系数计算卫星瞬时轨道根数。由于地球引力场的J2项摄动和大气阻力摄动比其他项大的多，在使用平根数计算时只需要考虑地球引力的J2项和大气阻力项，因而公式较为简单，计算量不大。而月球引力场的J2项摄动的量级与其他高阶项接近，如果采用平根数法进行计算，则需要考虑20阶月球引力和地球引力摄动，即使经过简化后计算量仍然非常大，不适合星上计算机进行计算。在月球卫星环月期间，定姿需要卫星在月球坐标系内实时的轨道根数，对卫星轨道计算精度提出的要求是位置误差小于化m。如果采用传统的平根数法，要做到较高精度的轨道计算，由于月球引力场的复杂性，计算时采用的月球引力模型就需要取到20阶左右的较高阶次，另外在较高精度要求下还需要考虑地球引力摄动的影响，从而导致计算量过大，不适合在星上计算。而数值解的方法尽管精度较高，但也存在计算量较大的问题。
发明内容本发明的技术解决问题是克服现有技术的不足之处，提供一种地面运算与星上运算相结合的方式，使得环月卫星能够在星上实施计算出环月轨道。本发明的方法的技术解决方案是一种星上实时计算环月轨道的方法，包括以下步骤(1)以定轨结果作为初值，用RKF78数值方法计算出实际轨道；(2)选取一条二体轨道作为参考轨道；(3)计算实际轨道相对于参考轨道的运动轨迹，采用三角函数来拟合实际轨道相对于参考轨道的运动；(4)根据公知的两条二体轨道之间的相对运动规律，计算出一条相对道；步骤(1)至步骤(4)均在地面完成，将一天分成12个时间段，分别重复利用步骤(1)至步骤(4)计算相应的二体轨道，然后将计算结果注入到星上；(5)卫星采用二体轨道计算方法，根据步骤(4)中算得的二体轨道实时计算出卫星的轨道根数。本发明与现有技术相比的有益效果是本发明通过地上运算与星上运算相结合的方式，特别是选取了一条二体轨道作为参考轨道，使得星上在计算量很小的情况下能够实施计算出卫星轨道，同时不需要地面和卫星保持实时联系，这种轨道计算方法的位置误差小于1km。图1为本发明星上计算轨道的水平面内误差情况；图2为本发明星上计算轨道的轨道面内误差情况。具体实施方式(1)在地面通过高精度的数值运算计算出月球卫星在未来一段时间内的月心惯性系轨道根数，以一定的步长h保存(hK60s),轨道计算的初值由定轨所得。因为摄动原因得到的结果随时间而变化，以^表示此根数，并认为此根数表示了卫星的实际位置。计算过程需要考虑的摄动包括月球引力非球形引力摄动，其阶次应取到30阶以上，同时还需考虑地球引力、太阳引力等对月球卫星的摄动。计算方法采用高精度的数值解法。例如，取卫星定轨的结果为起始时间(北京时)t0=2007年11月7日12:0:0.000卫星轨道根数为半长轴a=1934km;偏心率e=0.00237;倾角i=87.41度；升交点赤经0=265.51度；近地点幅角w=281.585;平近点角M=202。精确轨道计算采用RKF78数值方法，考虑月球引力场为70x70阶、考虑地球引力摄动、太阳引力摄动。(2)选取一条二体轨道作为参考轨道，使得上述计算出的轨道相对于该二体轨道在小范围内运动。可以选择开始时刻对应的二体轨道作为参考轨道，其轨道根数记为cr。。参考轨道选择为半长轴a(m)偏心率e倾角i(度)升交点赤经Q(度)近地点幅角U)(度)平近点角M(度)19340.0023787.41265.51281.585202(3)计算实际轨道相对于参考轨道的运动轨迹。利用坐标变换可以得到tr,相对于o"。的运动轨迹<formula>formulaseeoriginaldocumentpage5</formula>其中<formula>formulaseeoriginaldocumentpage5</formula>0凡，n表示旋转矩阵下标t表示(1)中计算出的数值解轨道A，下标o表示参考轨道，、r'分别表示^、A的月心距。(4)对得到的相对运动结果^。，^。做最小二乘拟合，对x方向和y方向的相对运动采用的拟合函数分别为y="sinM+&2cosM其中，在一个拟合时间段内(^确定的周期7;始，以//"。=、%的步长增加.2;r)M取值从0开解法方程A"Ac-A、得到拟合公式的系数c-(c。，Cl，C2，c3f，其中sinOcosO、1sinM,cosM,X=IXlVX2S。，…，X"S(ja=1M2sinM2cosM2、1sinAC,cosM—u解法方程A"Ab=A、得到拟合公式的系数b=(6l'62)r，其中，、in0cosO、/、rcosilyy=(ms0,>^0，…，Ja=sinA/2cosM2cosily上面A矩阵中的Mt(k=1，2…n-1),"o，由拟合系数计算相对轨道根数a0，w00=w0+M00其中Moo为参考轨道CT。在一个拟合时间段初始时刻的平近点角M(c2cosw00+c3sint/00)+1=2。0(c2sinm00—c3cosw00)2o0由相对轨道根数求拟合轨道根数:<formula>formulaseeoriginaldocumentpage7</formula>每组轨道的使用时间H取为1或2个平均轨道周期，对12个拟合时间段分别使用以上(1)至(4)中所述的方法进行拟合，得到12组轨il^艮数，在加上轨道参数初始时刻对应的星时to、使用时间段对应的轨道半长轴a0两个参数，最终形成上行注入的74个参数，分别用以下符号表示<table>tableseeoriginaldocumentpage7</column></row><table><table>tableseeoriginaldocumentpage8</column></row><table>具体计算结果为<table>tableseeoriginaldocumentpage8</column></row><table>(5)卫星采用二体轨道刮-算方法，根据(4)中算得的二体4九道实时计算出卫星的轨道根数。计算公式为<formula>formulaseeoriginaldocumentpage8</formula><formula>formulaseeoriginaldocumentpage9</formula>惯性系至轨道坐标系转换矩阵记为<formula>formulaseeoriginaldocumentpage9</formula>AOI用作卫星相对轨道坐标系进行姿态确定的输入。如图1、2所示，为星上使用地面注入的二体轨道进行计算与地面高精度计算结果相比的水平面内误差情况和轨道面内误差情况。图中x轴表示卫星飞行方向，y轴表示轨道负法线方向，z轴表示月心方向，三个轴的单位都为km，从图中可见星上计算卫星x方向误差小于200m,y方向误差小于400m，z方向误差小于300m。权利要求1.一种星上实时计算环月轨道的方法，其特征在于包括以下步骤(1)以定轨结果作为初值，用数值方法计算出实际轨道；(2)选取一条二体轨道作为参考轨道；(3)计算实际轨道相对于参考轨道的运动轨迹；(4)根据公知的两条二体轨道之间的相对运动规律，计算出一条相对于参考轨道运动轨迹与实际轨道相对于参考轨道运动轨迹相接近的二体轨道；步骤(1)至步骤(4)均在地面完成，然后将计算结果注入到星上；(5)卫星采用二体轨道计算方法，根据步骤(4)中算得的二体轨道实时计算出卫星的轨道根数。2.根据权利要求1所述的星上实时计算环月轨道的方法，其特征在于所述步骤(1)中的数值方法为RKF78数值方法。3.根据权利要求1所述的星上实时计算环月轨道的方法，其特征在于所述步骤(3)中计算实际轨道相对于参考轨道的运动轨迹，是采用三角函数来拟合实际轨道相对于参考轨道的运动。4.根据权利要求1所述的星上实时计算环月轨道的方法，其特征在于:将一天分成12个时间段，分别重复利用步骤(1)至步骤(4)计算相应的二体轨道，再注入星上计算相应的轨道根数。全文摘要一种星上实时计算环月轨道的方法，涉及航天器实时计算环月轨道
技术领域：
，包括以下步骤(1)以定轨结果作为初值，用数值方法计算出实际轨道；(2)选取一条二体轨道作为参考轨道；(3)计算实际轨道相对于参考轨道的运动轨迹；(4)根据公知的两条二体轨道之间的相对运动规律，计算出一条相对于参考轨道运动轨迹与实际轨道相对于参考轨道运动轨迹相接近的二体轨道；步骤(1)至步骤(4)均在地面完成，然后将计算结果注入到星上；(5)卫星采用二体轨道计算方法，根据步骤(4)中算得的二体轨道实时计算出卫星的轨道根数。本发明使得星上在计算量很小的情况下能够实施计算出卫星轨道，同时不需要地面和卫星保持实时联系。文档编号G01C21/24GK101226062SQ20071030174公开日2008年7月23日申请日期2007年12月26日优先权日2007年12月26日发明者红宗,张洪华,李铁寿,哲王,寨王,王大轶,王淑一,綦艳霞,冬韩申请人:北京控制工程研究所