轨道方程计算中Adams_Cowell方法外推改进.pdfVIP

中 国 空 间 科 学 技 术 2006 年 2 月 22 C H IN ESE SPA CE SCIEN CE AND T EC HN OL O GY 　　　 　　　第 　1 　期 轨道方程计算中A dam sCowell 方法的外推改进 付兆萍 　李红 (华中科技大学 , 武汉 430074) ( 摘要 　在人造地球卫星精密定轨中 , 有摄星历的计算常采用 AC 方法 即 A dam s ) Cowell 方法 , 针对这种经典的线性多步法 , 应用外推技术将计算过程进行了加速 , 从而 ( ) 提出了一种更高精度的 WAC 方法 即 A dam sCowell 方法的外推改进 。最后应用典型算 例将 AC 方法 、KS G 积分器 、改进的 AC 方法及文中的 WAC 方法进行数值验证 , 并比 较 4 种方法中卫星轨道沿迹误差的结果 , 结果表明采用 WAC 方法的精度有明显提高 。 主题词 　轨道计算 　运动方程 　计算方法 　航天器 1 　引言 在卫星轨道理论中 , 描述航天器轨道运动的微分方程相当复杂 , 除二体问题等少数几种情况 外 , 都不可能给出严格解 , 故计算常常采用数值解法 , 其中包括单步法和线性多步法 。单步法以 RungeKut t a 方法较为常用 , 它使用起来方便 , 且积分起步容易 , 同时又可以变步长计算 ; 但是随 ( ) 着方法阶数的增高 , 计算量也相应增大 右函数的计算 。由于人造地球卫星定轨问题需要处理的资 料多、计算量大 , 因此往往采用多步方法 , 多步法中常用的有 AC 方法 。Kro gh Shamp ine Gor don 在 AC 方法基础上稍作改进给出了一种名为 KS G 积分器的计算方法 , KS G 积分器校正过程的一大 特点是计算简单 , 不但可以改善一些舍入误差的状况 , 而且明显提高了计算效率 , 仅为 AC 方法 的3/ 4 , 但是 KSF 积分器计算方法随着方法阶数的增高反而使校正效果变差 。将 KS G 积分器采用 的校正过程 , 用于 AC 方法 , 就有了改进的 AC 方法 。改进后的 AC 方法 , 计算过程与 KS G 积 分器同样简单 , 耗时少 , 舍入误差小 , 但是对于提高精度确无益 , 特别对高阶方法更为不利 。于是 对于多步法 , 实际应用表明 , 还是采用一般的 AC 方法更好 。除了以上常用的线性多步法以外 , [ 1 ] [2 ] 还有近些年来研究的对 Hamilto n 系统数值计算新建立的辛算法 和线性对称多步法 ; 但是对称 算法在某些时候的确表现了数值不稳定性 。 综上所述知卫星轨道方程计算中已发展了多种数值方法 , 每种方法都有自己的优点和缺点 , 所 以弥补算法的不足提高算法的精度就成了重要的研究方向[3 ]