【原创】八大行星的牛顿引力摄动轨道近日点进动的精确计算结果

轨道升交点进动的牛顿引力理论计算值是(括号内是观测值，单位是每百年角秒)：
水星 -451.4844(-446.30)， 金星 -994.2752(-996.89)， 地球 -5754713.4(-18228.25)，火星 -1059.5182(-1020.19)， 土星 645.1411(1217.17)， 木星 -928.1120(-1591.05)， 天王星 279.8878(-1681.40)， 海王星 -23.1862(-151.25)。

原来用Mathematica 软件进行数值积分很简单，根本就不需要作任何人为级数展开(展开后的项目太多太多计算量太大反而运算速度太慢了)。
八大行星的牛顿引力摄动轨道近日点进动和升交点进动的计算结果分别是(括符内是升交点进动值，单位是每百年角秒)：
水星 528.7693(-451.4845)， 金星 10.9329(-1000.7622)， 地球 1143.2093(-579295.583)， 火星 1593.5351(-1062.3763)， 土星 643.0834(645.1460)，木星 1659.1982(-939.9329)， 天王星 310.9990(272.9742)， 海王星 71.0442（-22.8706)

以上数据表明远不只水星轨道近日点进动值理论计算值与观测值偏差较大，而且所有行星的轨道近日点进动和升交点进动理论计算与观测值偏差都很大。其中，地球轨道升交点显然是观测位点不确定性造成的，而土星木星存在长周期轨道共振，其它所有行星难道不都与观测的不确定性有关？如果我们把所有观测数据都从太阳中心变换到太阳系质心为基点的结果又将如何？可惜我没有原始观测数据无法完成这种坐标系变换，但至少可以证明相对论是荒谬透顶的。

行星进动，源在太阳运动，太阳外空间，电磁力不均分布；地球与其它行星，也有自己的电磁空间不均分布，一年一度新轨道，新旧轨道起点变化是进动。进动，不是参考系观测所为的。

以上是拉格朗日摄动方程，用Mathematica 软件通过八大行星轨道坐标函数、与时间相关的数百个反函数和隐函数进行直接数值积分的分析数据图。图示的主波是三大内行星轨道长半轴在土星方向受摄结果，土星位于远日点方向的摄动值为正，近日点方向为负，这就是水星轨道近日点进动观测值与牛顿理论计算值不符的真相。

终于成功用平近点角设置隐函数方法完成计算。过去所有计算都是对轨道真近点角求平均值，而百年前的手工计算是用开普勒公式展开成平近点角计算。尽管今天得到的公式在Mathematica中进行数值积分结果与以前多次不同方法计算结果基本相同，但由于与八大行星的观测值都有较大差异，这里有必要再公布一次。
八大行星轨道近日点进动的牛顿引力摄动值分别是(括号内是观测值，单位是每百年角秒)：
水星 528.7694“(573.57)，金星 10.9334”(-108.80)，地球 1143.2093“(1198.28)，火星 1593.5351”(1560.78)，土星 643.0834“(839.93)，木星 1659.1982”(-1948.89)，天王星 310.9990“(1312.56)，海王星 71.0442”(-844.43)。
轨道升交点进动的牛顿引力理论计算值(括号内是观测值，单位是每百年角秒)：
水星 -451.4847(-446.30)， 金星 -1000.4750(-996.89)， 地球 -5792955.85(-18228.25)，火星 -1062.3763(-1020.19)， 土星 645.1460(1217.17)， 木星 -939.9329(-1591.05)， 天王星 272.9742(-1681.40)， 海王星 -22.8706(-151.25)。

@绝世圣人 彬彬还有时间上网吗？Mathematica 比MATLAB功能强大多了，用起来也得心应手，就是因为MATLAB用起来太死板我才改用Mathematica 的。

重新审核推导了我6年前论文初稿中的相关公式，根据我的引力场以太理论引力质量与惯性质量不等推断得到的对水星轨道近日点进动的贡献是正值，每百年约为14.6“，根据10楼计算结果仍有约30“的差额。
但问题是，如果以这些数据为依据，地球的却还有约50”的差额又去了哪？还有其它各大行有的差额可能更大，是否与观测过程中的平近点角和真近点角换算位点差异有关？反正地球春分点的平近点角计算值并不在其轨道真近点角位置上。