是的,看到这里,天文专业的同学应该发现了一个问题:
某度小编把冥王星的数据计算成了77。2——这特么是太阳系内边界的距离……
实际上呢。
在计算过程中,由于k次多项式存在的缘故,冥王星和海王星是共用n=8来计算的。
所以根据提丢斯-波得定则计算,冥王星的误差率是2%,而非200%。
这是天体物理以及天体测量第二学期就会明确标注在课本上的内容,作为一个百科栏目居然会犯这种错误,也是挺无奈的……
上辈子徐云恰好有某段情节正好用到了提丢斯-波得定则,在骚扰……咳咳,咨询某位在凤凰山观测站工作的朋友时,对方一度对百科表达了某些极其亲切的问候与祝福。
当然了。
造成这种情况的很大部分因素要归结于知识的冷门,提丢斯-波得定则本身就是个小众知识,更别说冥王星这个小众中的小众了。
总而言之。
后世对于提丢斯-波得定则在数学计算的数值方面基本是没意见的。
它的主要争议在于物理意义模糊,是一个纯粹的经验公式,很难从原理上进行解释。
像an+1∶an=β之类的其他测定方式,基本上也都是数学方面精准,但物理意义不明的情况。
随后徐云又写下了两个个公式,也就是k次多项式的函数和最小误差值:
f(x)≈g(x)=a0+a1x+a2x2+a3x3+……+akxk。
loss=i=0∑10(g(i)-f(i))2。
这样一来。
只要找到合适的系数,就能令误差值最小了。
而就在徐云优化函数的同时。
其他人也没闲着,各自按着预定好的计划在行事。
例如老汤正和来自格林威治天文台的技术人员拍摄着今天的星图,高斯则整理起了布莱德雷家族留下来的独门观测记录:
“0。00066045……0。01072261……0。12684538……0。43146853……”
众所周知。
如果是需要仅仅通过数学来计算行星轨道数据,那么必然会用到开普勒行星三定律:
第一定律:
每一个行星都沿各自的椭圆轨道环绕太阳,而太阳则处在椭圆的一个焦点中。
第二定律:
在相等时间内,太阳和运动着的行星的连线所扫过的面积都是相等的。
也就是Sab=Scd。
第三定律则是:
各个行星绕太阳公转周期的平方,和它们的椭圆轨道的半长轴的立方成正比。
即T^2a^3=K,T为行星周期,K为常数。
另外还需要用到笛卡尔坐标系下的椭圆曲线,即:
Ax^2+Bxy+Cy^2+Dx+Ey+F=0。
有了这些,只要在加上某个工具就能进行计算了。
后世科技发达,计算轨道的工具一般是numpy,几秒钟就能计算出结果。
眼下虽然没有numpy协助,但这玩意儿的计算逻辑实际上就是最小二乘法。
而最小二乘法的发明者不是别人,正是高斯……
“g(x)=-0。43146853+0。12684538x-0。01072261x^2+0。00066045x^3……”
“下一组是0。31468531……0。21538462……0。12960373……”