E=∑(yi-F(xi))^2
它是一个多元函数,有an共n个未知量,现在要求的是最小值。所以必然满足对各变量的偏导等于0,于是得到n个方程:
n个方程确定n个未知量为常量是理论上可以解出来的。用这种误差分析的方法进行回归方程的方法就是最小二乘法。
2、最小二乘与投影
我这个人不喜欢看这些理论,公式推导,而更喜欢用例子来展示算法的思想。例如,在二维坐标系中,有三点,(1,1),(2,2),(3,2),那如何用一条直线来拟合这些点呢?
首先,我们可以假设直线表达式如下所示:
然后计算误差函数:
在求得误差函数E对系数a,b的偏导,并使之为0:
由上式得到系数a,b的值,并得到拟合直线表达式:
通过最小二乘法得到的曲线如下:
拟合的目的就是使得数据点都满足上述函数表达式,即:
用矩阵形式表示如下:
Aa=p是肯定有解的,因为p在矩阵A的列空间中。要使得e向量的长度最短,当且仅当p为b在矩阵列空间上的投影!有上一篇《正交投影》中投影矩阵的通式可得:
那么将p代入公式Aa=p可得:
将具体数值代入得:
则可以得到:
b,p,e向量分别可以表示如下:
p,b在图中的表示如下:
原文地址:http://blog.csdn.net/tengweitw/article/details/41745923