方程比未知数多的时候
四段弧、五个纬度,要定的只有两个数。多出来的那几条怎么办——他让误差的平方和最小。
1805 年那本讲彗星轨道的小书末尾附了八页东西,题目叫《论最小二乘法》。开头一句把问题说尽了:从观测里定几个量,往往得到一组形如 E = a + bx + cy + fz + … 的方程,方程比未知数多,误差不可能同时清零。怎么分摊?他说,在所有可以提出的原则里,没有比「让误差的平方和最小」更普遍、更准确、也更容易用的了——这样误差之间会建立起一种平衡,极端的那一个压不住其余的。接着他做了一件很现代的事:把方法用在一笔真账上。德朗布尔与梅尚刚量完敦刻尔克到蒙特惠奇的子午线,四段弧、五个纬度。地球若是椭球,一段弧长与两端纬度之间有一条关系,里面只有两个未知数:扁率与四十五度那一度的长。四条方程配两个未知数,配不上。他把五个纬度各自的改正数当成误差,求平方和最小,解得扁率 1/148、那一度长 28497.78 个模块(一模块合两托瓦兹)。1/148 离谱得很——今天的值约 1/298——他自己也说不能采信,因为摆长和天文给的是 1/320;于是把扁率钉死在 1/320 再解一次,得 28504.09,与米所依据的那个数 28504.10 只差 0.01。他顺手记了一句:若只认法国这一段的测量,米还得缩短大约四千五百分之一。这套办法就这样第一次印在了纸上。四年后高斯在《天体运动论》里用了它,并称自己 1795 年就在用。
拖动扁率与那一度的长,看五个纬度的改正数怎么跟着动;他自己解出的那两点各标一面旗
地球若是椭球,一段子午线弧的长度与两端的纬度之间有一条关系,里面只有两个未知数:扁率 α 与四十五度那一度的长。四段弧给出四条方程, 两个未知数配不上,于是他把五个纬度各自的改正数当成误差,让误差的平方和最小。 解出来的扁率是 1/148——离谱(今天的值约 1/298),他自己也说不能采信; 把扁率钉在摆长给出的 1/320 再解一次,那一度是 28504.09 modules, 与米所依据的 28504.10 只差 0.01。代价是误差变大:最大的一个落在埃沃,5.83 角秒。 他把这几秒归给局部的引力异常,并由此说了一句很有分量的话—— 定普遍度量单位,子午线弧不如摆长可靠。
