最小二乘法的来源:
最小二乘法(Least Squares method),又称最小平方法,起源于十八世纪的大航海探索时期,发展于天文领域和航海领域[1]。法国科学家勒让德于1805年首先提出了最小二乘法,后来在1822年由高斯证明了最小二乘法的优越性,也就是高斯-马尔可夫定理。现在最小二乘法已被普遍应用于误差估计、曲线拟合、参数辨识、预测、预报等领域。
理解最小二乘法:
为了帮助大家直观理解最小二乘法,我们以一个简单的例子,介绍最小二乘法的思想。
假设下面这20个数据是你用同一把尺子测量20次水杯高度的数据,每次的测量都是有误差的,那么水杯的高度到底是多少呢?
表120次水杯高度的测量值
一般情况下我们可以用数据的平均值作为水杯的高度:
利用上面的数据可求得高度的平均值为10.0015cm。
让我们用最小二乘法来思考刚才的问题。先将上面的数据画在坐标系中,如下图所示,其中红色虚线为水杯高度的估计值。
三轴加速度计的标定:
三轴加速度计有两种常用的标定方法,一种借助转台设备进行标定,该方法操作简单,标定准确,但需要额外购买设备,成本较高;另一种方法则不需要外部的标定设备,但精度略微差。本部分将介绍如何利用最小二乘法对加速度计进行数据的标定以提高其测量准确性。
三轴加速度计的测量精度受非正交误差、零漂、尺度因子等多源误差的影响,在标定时需要首先建立包含这些误差的输出方程。
1)首先考虑非正交误差,示意图如图2所示,相应的误差矩阵为式(5)。
这里最小二乘法辨识利用的原理是:在静止状态下,无论角度值如何变化,加速度计三个轴测得的重力加速度矢量与当地重力加速度矢量的模值相等[4]。
根据这个原理,我们可以求得的最小二乘估计:
为当地重力加速度值,为观测值的数量,一般情况下应远大于待辨识参数的个数。
至此,利用最小二乘法标定三轴加速度计的思路已经介绍完毕,式(10)的求解可以使用最速下降法、高斯牛顿法、遗传算法等最优求解的方法进行计算。