1. 线性拟合和二次拟合函数
最小二乘法本质上就是求一个事先定义一个函数,然后使用已知的采样点结果拟合函数的参数,使得所有采样点的均方误差最小。
即:
1. 线性拟合
我们假定拟合曲线为:
则有拟合误差为:
要使得拟合误差Q 最小,则我们有Q 对于参数a,b 的偏导均为0,因此,我们即有:
可以解得:
2. 二次拟合函数
类似的,我们可以得到二次拟合函数的最小二乘法的结果。
定义拟合函数为:
则有拟合误差:
同样,我们可以解得,在极值条件下,有:
亦即:
可解得:
类似的也可以推至n阶的情况,这里就不一一展开了。
3. a\cdot e^{bx}型函数形如a\cdot e^{bx} 的函数直接用最小二乘法倒是没法直接求解,不过可以通过一定的函数变换转换成n 阶函数形式,然后我们就可以仿照上述方式进行求解了。
比如这里的a\cdot e^{bx} 类型的函数,我们做函数变换,令z= ln(y) ,则有:
即可使用最小二乘法进行求解。
类似的,还可以有函数y = \frac{1}{ax+b} ,我们令z = \frac{1}{y} ,即可有:
2. 解矛盾方程组
书中这一小节事实上就是给前面最小二乘法的内容提供一些理论上的支持,没有啥更多的内容,因此,我们仅在这里摘录书中的定理如下:
定理3.1
(1) A 为m 行n 列的矩阵,b 为列向量,A^TAX = A^Tb 称为矛盾方程组AX = b 的法方程,法方程恒有解。
(2) X 是min||AX-b||_2^2 的解,当且仅当X 满足A^TAX = A^Tb ,即X 是法方程的解。