资源描述：
《数值分析之最小二乘法与最佳一致逼近.ppt》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、§4曲线拟合的最小二乘法1最小二乘法及其计算在函数的最佳平方逼近中如果只在一组离散点集上给定，这就是科学实验中经常见到的实验数据的曲线拟合.1记误差则的各分量分别为个数据点上的误差.问题为利用求出一个函数与所给数据拟合.2设是上线性无关函数族，在中找一函数，使误差平方和这里3这个问题称为最小二乘逼近，几何上称为曲线拟合的最小二乘法.用最小二乘求拟合曲线时,首先要确定的形式.确定的形式问题不仅是数学问题,还与问题的实际背景有关.通常要用问题的运动规律及给定的数据进行数据描图,确定的形式,然后通过实际计算选出较好的结果.4为了使问题的提法更有一般
2、性，通常在最小二乘法中考虑加权平方和这里是上的权函数，它表示不同点处的数据比重不同.就是次多项式.若是次多项式，的一般表达式为线性形式.5这样，最小二乘问题就转化为求多元函数的极小点问题.用最小二乘法求拟合曲线的问题，就是在中求一函数，由求多元函数极值的必要条件，有使误差取得最小.6若记上式可改写为这个方程称为法方程，可写成矩阵形式7其中要使法方程有唯一解，就要求矩阵非奇异，而在上线性无关不能推出矩阵非奇异，必须加上另外的条件.8显然在任意个点上满足哈尔条件.哈尔条件，则法方程的系数矩阵非奇异，如果在上满足函数的最小二乘解为定义10设的任意线
3、性组合在点集上至多只有个不同的零点，则称在点集上满足哈尔(Haar)条件.方程存在唯一的解从而得到于是9这样得到的，对任何的都有故确是所求最小二乘解.10一般可取，但这样做当时，通常对的简单情形都可通过求法方程得到给定的离散数据，求解法方程时将出现系数矩阵为病态的问题，我们在下面考虑用正交多项式的方法解决。11例1已知一组实验数据如下，求它的拟合曲线.12解从图中看到各点在一条直线附近，故可选择线性函数作拟合曲线，将所给数据在坐标纸上标出，见图3-4.图3-413令这里故14解得可得方程组于是所求拟合曲线为15关于多项式拟合，Matlab中有
4、现成的程序其中输入参数为要拟合的数据，为拟合多项式的次数，输出参数为拟合多项式的系数.利用下面的程序，可在Matlab中完成上例的多项式拟合.16x=[11233345];f=[444.566688.5];aa=poly(x,f,1);y=polyval(aa,x);plot(x,f,’r+’,x,y,’k’)xlabel(‘x’);ylabel(‘y’);gtext(‘y=s1(x)’）17结果如下：18有时根据给定数据图形，其拟合函数表面上不是线性模型的形式，但通过变换仍可化为线性模型.例如，，若两边取对数得此时，若令这样就变成了线性模型
5、.19例2设数据由表3-1给出，用最小二乘法确定及.解表中第4行为通过描点可以看出数学模型为它不是线性形式.用给定数据描图可确定拟合曲线方程为两边取对数得20若令先将转化为为确定，根据最小二乘法，取则得数据表见表3-1.得21故有法方程解得于是得最小二乘拟合曲线为22利用下面的程序，可在Matlab中完成曲线拟合.x=[1.001.251.501.752.00];y=[5.105.796.537.458.46];y1=log(y);aa=poly(x,y1,1);a=aa(1);b=exp(aa(2));y2=b*exp(a*x);plo
6、t(x,y,’r+’,x,y2,’k’)xlabel(‘x’);ylabel(‘y’);gtext(‘y=a*exp(bx))’;23结果如下：242用正交多项式做最小二乘拟合如果是关于点集用最小二乘法得到的法方程组，其系数矩阵是病态的.带权正交的函数族，即（5.6）25则方程的解为且平方误差为26接下来根据给定节点及权函数构造带权正交的多项式.注意，用递推公式表示，即这里是首项系数为1的次多项式，根据的正交性，得27下面用归纳法证明这样给出的是正交的.28假定对及要证对均成立.有由的表达式，有均成立，29而，当时，另外，是首项系数为1的次
7、多项式，它可由由归纳法假定，当时的线性组合表示.由归纳法假定又有30由假定有再考虑利用表达式及以上结果，得31至此已证明了此多项式组成一个关于点集的正交系.用正交多项式的线性组合作最小二乘曲线拟合，只要根据公逐步求的同时，相应计算出系数最后，由和的表达式（5.11）有32并逐步把累加到中去，最后就可得到所求的用这种方法编程序不用解方程组，只用递推公式，并且当逼近次数增加一次时，只要把程序中循环数加1，其余不用改变.这里可事先给定或在计算过程中根据误差确定.拟合曲线33以上述例1为例先求正交系34例3设X={1.00,1.25,1.50,1.
8、75,2.00},在X上定义内积5(f,g)=∑xif(xi)g(xi)i=11)在函数系{1,x2}中求一个X上的正交函数系.2)用最小二乘法求一个形如y=a+b
