体问题的几种数值算法比较.pdfVIP

[文章编号】 1001．246x(2006105_0599_05 Ⅳ体问题的几种数值算法比较 杨远玲1， 聂清香1， 吴晓梅2， 徐顺福1 (1．山东师范大学物理与电子科学学院，山东济南250014；2．泰山学院，山东泰安271021) [摘要】 中应用时引起的能量误差、半长径和偏心率的变化进行比较．结果发现：RKF法精度最高，但长时间内有误差 积累；辛算法无人工耗散，能较好保持能量误差的稳定性；厄米算法虽然误差较大，但构造简单，耗机时较少． [关键词】 H蛐iIton系统；RKF法；辛算法；厄米算法；二体问题；Ⅳ体问题 [中图分类号】 P132 【文献标识码] A 0 引言 在天体力学中，行星系统的起源和演化问题一直是大家关注的重点．对这类问题，需要解Ⅳ体运动的微 分方程．但这些方程一般都相当复杂，除了二体问题等几种少数情况外，都不能得到严格的分析解．因此，数 值积分在天体力学中格外重要． 目前，数值积分已被广泛地应用在天体力学的Ⅳ体问题中．从人造卫星的定轨问题，到太阳系的演化问 题，再到对新发现的太阳系的Kuiper带天体和外太阳系行星系统的研究，都离不开数值模拟的方法．在天体 力学的Ⅳ体问题进行数值积分的过程中，通常遇到2个主要困难，一是积分误差的积累较大，二是积分时间 积分，如果积分方法选择不当，只考虑适当的机时，保证不了精度，那么误差的积累会使结果失真；若只考虑 积分精度，积分时间又相当长，从而使计算几乎不能实现．因此，选择什么样的积分方法，才能既保证积分结 果的精度又节省机时，是Ⅳ体问题数值积分中首先要解决的问题． 达到较高的精度，故在制定天文历表和人造卫星的精密定轨上经常使用．但由于其长时间的积分不能保持自 治系统的能量守恒bo，使得能量误差随时间线性增长，从而导致天体的精度误差随时间呈平方增长，所以在 多粒子的Ⅳ体演化问题中不很适用． Ⅳ体问题大部分都是Hamilton系统，而传统的算法却不能保持这个系统的一些特性．鉴于此，一种新的 算法——辛算法(symplectic 系统的辛几何结构，自动维持系统的能量，不存在人工耗散，在研究行星系统演化的收敛性和稳定性问题上 显示出优越性№j．但其计算精度较传统算法要低，尤其在短时间段内不能显示其优越性．所以，辛算法经常用 于多行星系统的演化模拟计算．针对截断误差较大的问题，许多研究者对其进行了改进b’9叫1|，使得辛算法 在计算上得到了简化b]，在较长积分步长和较长积分时间数值计算中的精度也得到了明显提高归1． 厄米算法(Hemlite 加速度的差值多项式中的一种四阶预估修正算法．由于普通的预估修正厄米算法在应用于长时间的模拟中 精度不高，EiichimKokubo等提出了一种时间对称的厄米算法，它的主要特征是没有长期的误差积累¨3|． 本文对上述可用的几种方法进行了试验分析比较，分别将它们应用于二体问题和Ⅳ体问题，并进行了 长时间计算，从能量误差、轨道半长径和偏心率变化等方面进行比较，不但确定了研究Kuiper带天体所需要 的最佳方法，也为其他有关的研究工作提供了可行的参考方法． [收稿日期】2005一05—24；[修回日期】2005—12一16 【作者简介]杨远玲(1979，)，女，山东省青州，研究生，从事天体物理研究，山东省济南市文化东路88号250014 万方数据万方数据 计 算 物 理 第23卷 1算法原理 1．1 RKF法 设天体运动的常微分方程为 警=以川)，戈∈彤，t∈尺， (1) 泰勒展开，用在积分区间[t。，t川]上若干个点的右函数值厂的线性组合来代替，的导数，然后用泰勒展开确 定相应的系数．常用的7阶RKF法(即8阶改正，简记RKF7(8))的积分公式为㈨ 菇州=菇。+^∑ciI|}；+D(^3)， ；州=¨^擎后t+o(^9)，f=0 ＼二，(2) 后o=，(￡。，戈。)， |i}一八f