建立同步卫星精密星历表的最佳数值方法.pdf

总第32期 I993年 I2月 人造卫星观测与研究 建立同步卫星精密星历表的 最佳数值方法 徐 劲 王-L{1 (中国科学院陕西天文台，西i安临i套．710600) 主题谒： 同步卫星 精街星历表 数值方法 _ — — 一 — — — — — — — — 一 - - — _ 一 — 提 要 本文简要舟绍了当前走体力学 中常用的各种敷值计算方法，蛄合同步卫星运动方程 的特点 和轨道解的性质。指出了各种敷值计算方法在 同步卫星情况下使用的优劣。最后通过有效的数 值实验 。确定了一次和分形式的Cowcll方法是建立同步卫星精密星历表 的最佳方法。这种方法 能适用于 同步卫星厘米叛定位精度的需要。妾验姑果表明： 在 Ikm 的精度要求下．瘟采用 I3阶方法。 在 lOOm～lm 的精度要求下．应采用 14阶方法． 在厘米级精度要求下，应考虐采用 l5阶方法 1 天体力学中常用的数值方法 随着天体力学研究的 日益广泛和深人。它所涉及到的计算手段也在不断发展．近年来 篙 又ffJ现了一些新 的数侬 计算方法．日前在国际天体力学研究中使用的主要数值方法可归结 2 为表 I(见 p51)。 由表 l可以看 出：天体力学计算中已发展了多种数值方法，大多数方法都有若干种演 变形式，以适应 同的要求或使计算更有效．每种方法都有各 自的优点和缺点，究竟哪一 种方法最适合同步卫星精密星历表的计算呢?这就需要对同步卫星运动方程的特点，轨道 解的性质 以及各种方法的优劣属性作一具体分析． 2 同步卫星运动方程与轨道解的特点 ¨步 I丛运动方程叮表示为 ： 一 十 第 2坤j 棣 劲：建 训司步。l￡精 帐 裘的最代数fi『 让 F 是摄 动 力，它包括 以 F几项： a．地球引力场搬动 ； b．LI爿搬动； c．光 搬动； d．附加搬动 。 当只求同步一ll星轨道 一阶锯时，可 以4考虑附 u搬动 “。出方程 (1)IIl以肴 III1 步：lJ星运动方程是二阶常微分方程组。当然也呵化为一阶常微分方程组。其右函数非常复 杂， 所包 括的蛸种_I}{动力通常情形 F都nT视 为保守力，H有当I司步1J 进地影时，光 搬动 能视为保守力，但这种情况H发生存奋分和秋分前后很短一段时阀内，对 ’J『星运动 产生的影响很小，以F不妨将砌 为保守力。因此同步卫星运动方程是一个保守系统，右 函数不显含；’。南同步_J『星的作用可知其轨道解是近阐轨道，偏心率 租小． 袭 I 常用敷疽方法分类 3 最佳方法的选择 根据前面对同步 1J型运动方程及轨道解的讨论 ，可知对同步卫星的数位计算4：需要变 阶、变步长的处理，所选取的数值方法应适用于问步卫星右函数较为复杂的汁算，在此前 提 F，数值方法的效率越高越好 。因此 MDMS方法 ‘’与级数展开方法都 适于埘同步 的数值训算，其原因是复杂的右函数的高阶导数递归公式推导幽难．在其余的各种方 法中，究竟选取哪一种方法最好呢?应主要考察各种方法的效率．当然考察应建立在方法 52． A造]J 观剐与蚵f究