head

    由于我们要处理的是很大的距离和很长的时间单位,因此在这种情况下用天文学上常用的天文单位(AU)来表示距离,高斯年(Year)来表示时间会比较方便()。 这时公式(2)可写为 :
             
    2)四阶龙格-库塔法解微分方程
    四阶龙格-库塔法[2]是解带初值的常微分方程最普遍的数值方法,与欧拉方法相比,它的精度更高,相应的推导过程也比较复杂,在这里我们只把结果介绍给大家。
在物理学中,多数的运动学问题都能归纳成以下形式:
            
,则可得到两个一阶方程:
              
用四阶龙格-库塔方法求得方程组(4)的解为:
            
其中:
            


注:1976年,国际天文学联会把一天文单位定义为一颗质量可忽略、公转轨道不受干扰而且公转周期为365.2568983日(即一高斯年)的粒子与一个质量相等约一个太阳的物体的距离。所以一天文单位即1.4959787066 × 1011米

 

bottom