进动是椭轨道长轴方向的转动现象，具体表现为行星绕行一周到达近日点位置时刻大于360°，与轨道速度是否加快或减慢没有直接关系。月球使是使地球的自转速度减慢，是因为引潮力的摩擦作用使地球转动能量减小的结果，同时摩擦力又把部份能量直接传递给月球，使月球获得能量而远离地球。月球引潮力对地球赤道隆起部份的另一种作用与行星轨道相似而产生进动现象，这就是岁差和章动。

    这个画图没用，全靠非常复杂而严密的数理逻辑，一环套一环来不得一丝差错，经常是一个符号输错就推倒重来(我划分为几百个各自独立的小部份分开计算)。我没去学Matlab软件，而是用Excel程序，根本就无法用数值方法计算。在Excel中只能显示9位有效数，而将勒让德级数展开到20级以上，恐怕输入16位有效数值的计算结果也会出差错，起初的近一个月时间我改用计算器将数值计算输入到12位以上根本查不出任何其它误差来源，全部改用原公式输入才消除了计算误差。现在符号输入也找到了系统方法，出差错的概率小了，主要时间花在逻辑组合的输入上。

http://tieba.baidu.com/f?z=908454022&ct=335544320&lm=0&sc=0&rn=30&tn=baiduPostBrowser&word=%CF%E0%B6%D4%C2%DB&pn=0>

老刘：

我发现一个问题。23楼CloudK>的推导有个问题。前几式都没问题，就是从（3）到（4）的过程有问题。

俞允强的《广义相对论引论》也是类似的能量守恒法。他的6.3.2式为：

1/2(dr/dt)^2=E+GM/r-L^2/(2r^2)。错了，应为E/m（E就是上面的C）。

利用dφ/dt=L/r^2可化为上面（3）式。接着可以化为：

1/2[d（1/r )/dφ]^2=E/L^2+GM/(rL^2)-1/(2r^2)……6.4.1式。少了个m。

接着就出现了同样的问题，把6.4.1式两边对φ求微商，就能得到（4）式。这根本不能成立吗。

d《[d（1/r )/dφ]^2》/dφ=d《2E/L^2+2GM/(rL^2)-1/r^2》/dφ会得到（4）式？？？

你说的乘以u^4，那能把（3）式变换成6.4.1式(E那项还少个m）。接下来就是将6.4.1式两边对φ求微商啊

我算过了不能化简。你算一下,最好给个过程。
有问题的东西才会向你确认一下。你也知道，我发生误判的情况也是较少的。

※※※※※※
天地之道，以阴阳二气造化万物。是故易有太极，是生两仪。两仪生四象，四象生八卦。

    对不起，还是你说的对。这既不是我的计算问题也不是参考系搞没搞清楚的问题，而是观测时所处黃经位置的归算问题！

    根据牛顿平方反比力的摄动计算，行星轨道的近日点进动值在不同黃经位点的观测差异非常之大，再把数据归算到某一纪元时刻也就必然会出计算误差。也许我得先直接计算实际近日点时刻的观测值，再用开普勒平近点角展开式归算到某一纪元时刻的理论值。现在我直接用2000J时刻数据计算得到的值分别是：水星 528.7439"，金星 83.1211"，地球 1190.3872"，火星 1571.149"。还有很多计算没有完成，相关研究将继续进行。