月球轨道计算计算机可精确地求解星球运动微分方程。本算法适用于用计算机求解任意多体问题的星球运动微分方程，也能精确地确认万有引力公式本身和实际情况的相差程度。一．计算机计算星球轨道原理。以地球绕日运行为例，如图1-1。地球受力如图1-2。 月轨图改1-1,2.gif (2.3 KB)2007-5-11 23:00为叙述方便，以下用矢量形式表达。 文字图1.gif (3.12 KB)2007-5-11 23:00根据a=dv/dt,v=dr/dt,我们可按等时间间隔dt（即等步长），以微分形式从地球的初值点逐点向下推算。设t=0时，地球的初值点为r0,v0和a0,于是，地球经dt时间从初值点到达第一点,递推式为, 月轨图补7.gif (743 Bytes)2007-5-11 23:00由于dt是人为设定的，是已知的，因此地球到达1点的近似值v,r和a可由上式算出，算出1点值后，可把1点值作为初值，按步长dt继续推算出下一点的值，如此，可推算到第n点。由于dt值取得越小，递推的精度越高，我们可据此来控制计算误差。设要计算地球在t=T时的r值，要求计算误差为e，t=0时的初值r0,a0,v0为已知。我们可将0到T的时间间隔划分为n个dt，即令计算步长dt=T/n,然后根据上述，按步长dt从t=0时的初值点推算到t=n•dt=T时的r值。然后将dt二分，即令计算步长dt1=dt/2,再按此新步长值dt1从t=0时的初值点算到t=2•n•dt1=T时的r值为r2，比较一下二分前后的r值，即看一看是否满足条件r2 - r

r1,r2的初值为r10x,r10y,r10z,r20x,r20y,r20z,r的初值为r0x,r0y,r0z,地球的初速度为v10x,v10y,v10z,月球的初始速度为v0x,v0y,v0z。这样，地球和月球从各自的初值点同时出发，经dt时间后，地球就到达了它的下一点r1x,r1y,r1z,月球到达了它的下一点rx,ry,rz,于是可得如下递推式： 三维月轨推导式5.gif (9.92 KB)2007-5-11 23:00控制计算误差的6个条件为：|rx2-rx|

一帖一图吧，这样看太费力了

围观4楼

力学很好嘛，为什么要被淘汰。我所做的工作就是和力学打交道。现在造的房屋、桥梁，都是各种力学分析的结果。如果被淘汰了，那我们住的房子会不会塌？

都好厉害哈~

顶```

兴趣中……

用相对论公式算一下，看看结果如何？

刚高中 不懂这个 ...