曲线拟合最小二乘法一般最小二乘逼近曲线拟合最小二乘法.pptVIP

§5曲线拟合的最小二乘法 般的最小二乘逼近(曲线拟合的最小二乘 法)的一般提法是:对给定的一组数据(x,y) (=01…,m),要求在函数类9={93,…,}中找 个函数y=S(x),使误差平方和 |2=∑82=∑S(x)-y=min∑[S(x)-y2 其中S(x)=a090(x)+a192(x)+…+an9(x)(nm) 带权的最小二乘法: (x)S(x)-f(x 其中o(x)≥0是[a,b]上的权函数。  用最小二乘法求曲线拟合的问题,就是在S(x) 中求一函数y=S(x),使|2取的最小。它 转化为求多元函数 /(o,a1,…,an) o(x D∑a9(x1)-f(x) 的极小点(an,a1,…an)问题。由求多元函数极值 的必要条件,有 (x)∑a,(x)-f(x)2(x)=0  若记 (q,9)=∑(x)(x)92(x) f,q)=2o(x)f(x厘(x)≡dk(k=0,1 则上式可改写为∑(9,9)1=4(k=0 这个方程称为法方程,矩阵形式Ga=d 其中 (q0,90)(q,9)…(q,pn) G,)(,91)…(,9 (qn,90)(qn,91) 9n,9  由于9,1,…n线性无关,故(G≠0,方 程组存在唯一解 (k=0,1,…,n) 从而得到函数f(x)的最小二乘解为 S(x)=a09(x)+a191(x)+…+ann(x) 可证 O(xS(x)-f(x)2≤∑o(x)S(x)-f(x) 故S(x)是所求最小二乘解。  =74 (q,f)=O,f=4 (q,)=∑x=1455 由方程组 +22c1 2.77 22a0+74a1=145.5 1.13 所求拟合曲线为 S1(x)=2.77+1.13x  例:在某化学反应里,根据实验所得生成物的浓 度与时间关系如下表,求浓度y与时间t的拟合曲 线y=F(1 时间12345678910111213141516 t(分) 浓度4.6.8.8.9.9.9.9.10.10.10.10.10.10.10.10. ×10300400098022507086002032425055560 解:将数据标在坐标纸上,可发现数据符合双曲 线函数或指数函数。 1)双曲线函数拟合 双曲线型: 即 (at +b)  为了确定a,b,令y 由数据表1y生成数据表x,y.于是可用x的 线性函数S(x)=y=a+bx拟合数据(x,y;) (=1,…,16)。方法与上例一样解方程组 16a+3.38073b=1.8372×10 3.38073+1.58435b=0.52886×10, 得 a=80.6621.b=161.6822. (at+b) 从而有 (80.621+161.6822)≈F(, 其各点误差为 (1)=yi  2)指数函数拟合 拟合曲线形如y=ae两边取对数hy=ha+ 为了确定ab,令j=y,A=ma,x=1, 由(t,y)计算出x,),拟合数据的曲线仍为 S, (x)=y=A+bx 用上例的方法计算出 A=-4.48072,b=-1.0567, 从而 11.3253×10 最后求得y=113253×103c100=F2 各点误差为 y-F2(1)(i=1…,16)  3)两个模型的比较 本例经计算可得 nax6=0568×103,max2=027710 均方误差为 1.19×10 62)2=0.34×103 由此可知|2及=都比较小,所以用 y=F2(t)作拟合曲线较好。 确定拟合曲线的数学模型需要选择比较。