用Mathematica8计算原数据中因海王星的轨扁心率输入有误,行星轨道牛顿引力摄动重新计算后所得最后结果为:
近日点每百年进动值(括符内为观测值):水星528.769"(573.57"),金星10.934"(-108.80"),地球1143.22(1198.28"),火星1593.54"(1560.78"),土星641.59"(839.93"),木星1660.66"( -1948.89"),天王星291.435"(1312.56"),海王星82.626(-844.43")。
升交点每百年进动值(括符内为观测值):水星-451.484"(-446.30"),金星-1000.55"( -996.89"),地球-5792934.37"(-18228.25"),火星-1062.37"(-1020.19"),土星645.147"(1212.17"),木星-939.924"(-1591.05"),天王星270.448"(-168140"),海王星 -21.530(-151.25")。
因金星与地球轨道半径比值较大,用拉格朗日摄动公式展开到30级的计算精度还不够,同时有7重级数展开计算量太大电脑速度太慢,故采用高斯摄动公式20次迭代运算,所以重新计算轻松得到了足够高精度的计算结果。
相关人员不应再沉迷于水星问题的43"了,看看其它行星问题又如何?不过这里还有太阳扁率摄动和银河系自转数据未加入修正,是否还有其它更大因素对所有行星的影响不得而知。
@天涯树叶