摘要
倾角函数是天体力学分析理论中一种常用的函数.当把摄动方程展开成时间和根数的形式时需要用到.历史上提出了很多经典的倾角函数递推算法,并在双精度平台下开发了Fortran程序.进行了1次四精度计算倾角函数的试验,结果表明:L平面递推方法的四精度计算精度可达10^(-22),计算速度比双精度Jacobi方法快6倍.
In the analytical theories of celestial mechanics, the inclination function is a kind of frequently-used function when we expand the perturbation function as a form of time and orbit elements. A lot of classic recursion methods were brought out in history, and their Fortran-code were developed on the double-precision computer platform. While in this paper, a test of quad-precision calculation are reported. The results show that the quad-precision calculation precision of the L-plane recursion method can reach 10-22, and its computation speed is six times faster than the double-precision calculation of Jacobi method.
出处
《天文学报》
CSCD
北大核心
2013年第3期235-244,共10页
Acta Astronomica Sinica
基金
国家自然科学基金项目(11033009
11003049)资助
关键词
天体力学
方法
数值
celestial mechanics, methods: numerical