§4 曲线拟合的最小二乘法
1 最小二乘法及其计算
在函数的最佳平方逼近中 f (x) C如[a果, b],
f (x)
只在一组离散点集 {xi , i 0,1,上给, m定},这就是科 学实验中经常见到的实验数据 {(xi , yi ),i 0的,1, , m} 曲线拟合.
1
问题为利用 yi f (xi ),i 求0出,1一,个, m函,数 y S * (x) 与所给数据{(xi , yi ),i 0,拟1,合. , m}
7
Ga d ,
其中 a (a0 , a1, , an )T , d (d0 , d1, , dn )T ,
(0 ,0 ) (0 ,1)
G
(1
,
0
)
(1 , 1 )
(0 ,n )
(1
,
n
)
.
n (n ,0 ) (n ,1) (n ,n )
(k , j )a j dk (k 0,1, , n).
6
若记
m
( j ,k ) ( xi ) j ( xi )k ( xi ), i0 m
( f ,k ) (xi ) f (xi )k (xi ) dk i0 (k 0,1, , n).
上式可改写为
n
(k , j )a j dk
j 0
(k 0,1, , n).
这个方程称为法方程,可写成矩阵形式
要使j法0 方程有唯一解, 就要求矩阵 非奇G异,
而 0 (x),1(x),在 ,n上(x线) 性[无a,关b]不能推出
矩阵 G非奇异,必须加上另外的条件.
8
例如
(0 x) sin x,(1 x) sin 2x,显然(0 x),(1 x)线性无关。 取X {x0,...,x4} {0, ,2 ,3 ,4}, (0 x j ) (1 x j ) 0; j 0,...,4