积分器和程序设计.doc

紫金山天文台台刊稿件 PMOE 2003行星历表框架 （II）积分器和程序设计 李广宇 田兰兰 （ 中国科学院紫金山天文台，210008） 摘要: 积分器是实现历表数学模型的必要工具. 程序设计则是实现这两者，进行历表研究和完成实用历表的最终手段. 本文从历表研究需要和理论与实践结合的角度，统一整理了有关的知识和方法，详细介绍了张家祥积分器的理论基础，推导了实用公式，给出了有关数据结构、算法和程序设计的原则与方法，并对积分器的精度进行了分析. 关键词：积分器，数据结构，程序设计 数值算法 我们在文[1]中叙述了PMOE 2003精密行星历表框架的数学模型，对框架所用的积分器和程序结构及设计问题，只作了概要介绍. 作为历表框架的重要组成部分，积分器是实现历表数学模型的必要工具. 程序设计则是实现这两者，进行历表研究和完成实用历表的最终手段. 积分器和程序设计的质量直接影响历表的实现和质量. 有关这两部分内容的叙述散见于天体力学[2-5] 、计算方法[6,7]和程序设计[8-10]的文献中. 张家祥先生长期从事太阳系动力学模型数值积分的研究，积累了处理各种问题的丰富经验. 他改进Cowell方法建立的积分器具有快速高效和大幅减少舍入误差的显著优点，我们已在文[1]的引言中概要介绍了他的有关研究工作. 张先生关于积分器的工作在文[11,12]中已有介绍，但对学习者尚觉简略，并且没有涉及