。
这样一来。
只要找到合适的系数,就能令误差值最小了。
而就在徐云优化函数的同时。
其他人也没闲着,各自按着预定好的计划在行事。
例如老汤正和来自格林威治天文台的技术人员拍摄着今天的星图,高斯则整理起了布莱德雷家族留下来的独门观测记录:
“0.000660450.010722610.126845380.43146853”
众所周知。
如果是需要仅仅通过数学来计算行星轨道数据,那么必然会用到开普勒行星三定律:
第一定律:
每一个行星都沿各自的椭圆轨道环绕太阳,而太阳则处在椭圆的一个焦点中。
第二定律:
在相等时间内,太阳和运动着的行星的连线所扫过的面积都是相等的。
也就是sab=scd。
第三定律则是:
各个行星绕太阳公转周期的平方,和它们的椭圆轨道的半长轴的立方成正比。
即t/a=k,t为行星周期, k为常数。
另外还需要用到笛卡尔坐标系下的椭圆曲线,即:
ax+bxy+cy+dx+ey+f=0。
有了这些,只要在加上某个工具就能进行计算了。
后世科技发达,计算轨道的工具一般是numpy,几秒钟就能计算出结果。
眼下虽然没有numpy协助,但这玩意儿的计算逻辑实际上就是最小二乘法。
而最小二乘法的发明者不是别人,正是高斯
“g(x)=0.43146853+0.12684538x0.01072261x+0.00066045x”
“下一组是0.314685310.215384620.12960373”
“0.053379950.017249420.32307692”(注:所有数据都来自nasa开放的数据库,非杜撰)
过了大概十多分钟。
负责最终计算的黎曼抹了把额头上的汗水,在纸上写下了一个数字:
0.4857342657342658。
虽然目前还无法知晓冥王星的具体位置,更不知道它的重量大小。
但此前曾经提及过。
天王星在扣除海王星的引力之后