Cowell数值积分器的变步长与自起步方法  

The Method of Step Size Change and Self-startup of Cowell Integrator

作　　者：夏炎 田波[1] XIA Yan;TIAN Bo(College of Data Science,Tongren University,Tongren 554300,China;Purple Mountain Observatory,Chinese Academy of Sciences,Nanjing 210008,China)

机构地区：[1]铜仁学院大数据学院,铜仁554300 [2]中国科学院紫金山天文台,南京210008

出　　处：《科学技术与工程》2020年第17期6742-6747,共6页Science Technology and Engineering

基　　金：国家自然科学基金(61741214);贵州省教育厅项目(黔教合KY字[2016]051);铜仁学院博士科研启动基金(trxyDH1913)。

摘　　要：Cowell多步法积分器简单易用、计算高效,但也存在着步长变换困难和积分起步需借助其他工具等缺点,限制了其应用范围。为了使Cowell积分器更加方便实用,可以处理更多的问题,对积分器实现了以下3个方面的改进:首先实现了无精度损失的变步长操作,经实际验算,即便经过几十万次的反复步长变换,计算精度也没有损失;其次利用逐阶迭代的方式实现了积分器的自起步,摆脱了对其他工具的依赖;最后推荐了一种增加和分值存储字长的方法来减小舍入误差,误差减小效果显著,积分109圈以后相对误差仍小于1%。The Cowell multistep integrator with equal-space is easy to use and efficient,but difficult to change the step size and realize the self-start process.The disadvantage limits its application.In order to make Cowell integrator more convenient and handle more problems,the integrator was improved in the following three aspects.Firstly,a method was provided to change the step size without any additional error created.Through practical calculation,even after hundreds of thousands of times of repeated step transformation,there was no loss of calculation accuracy.Secondly,the self-starting of integrator was realized by the way of order-by-order iteration,and the dependence on other tools was eliminated.At last,a method to reduce rounding error was recommended by increasing and storing word length with significant effect on error reduction,and the relative error was still less than 1%after integrating 109 cycles.

关 键 词：数值积分 Cowell积分器 多步法 变步长 自起步 轨道计算 

分 类 号：P138.2[天文地球—天体力学]