计算方法-曲线拟合课件.ppt

《计算方法-曲线拟合课件.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、曲线拟合的最小二乘法直线拟合直线拟合多项式拟合一般最小二乘法的拟合应用由上述我们已经知到上述线性模型实际上是最小二乘法的推广，实际上也就是多项式逼近函数的问题。它不仅可以解决一元问题还可用于多元问题。除此外还可求解某些非线性问题。求解方法是将其通过一定的代数变换转换为可用线性模型求解的问题。比如对方程y=aebx取对数，得lny=lna+bx，令Y=lny,A=lna,B=b则问题转化为解Y=A+Bx的线性问题。类似的再如，对y=a+b/x拟和可对此方程取倒数，则新变量1/y于x成线性关系。线性模型引深及推广主页拟合与插值的关系函数插值与曲线拟合都是要根据一组数据构造一个函数作为

2、近似，由于近似的要求不同，二者的数学方法上是完全不同的。实例：下面数据是某次实验所得，希望得到X和f之间的关系？MATLAB(cn)问题：给定一批数据点，需确定满足特定要求的曲线或曲面解决方案：若不要求曲线（面）通过所有数据点，而是要求它反映对象整体的变化趋势，这就是数据拟合，又称曲线拟合或曲面拟合。若要求所求曲线（面）通过所给所有数据点，就是插值问题；曲线拟合问题最常用的解法——线性最小二乘法的基本思路第一步:先选定一组函数r1(x),r2(x),…rm(x),m

3、确定a1,a2,…am的准则（最小二乘准则）：使n个点（xi,yi)与曲线y=f(x)的距离i的平方和最小。记问题归结为，求a1,a2,…am使J(a1,a2,…am)最小。线性最小二乘法的求解定理：当RTR可逆时，超定方程组（3）存在最小二乘解，且即为方程组RTRa=RTy的解：a=(RTR)-1RTy所以，曲线拟合的最小二乘法要解决的问题，实际上就是求以下超定方程组的最小二乘解的问题。其中Ra=y（3）线性最小二乘拟合f(x)=a1r1(x)+…+amrm(x)中函数{r1(x),…rm(x)}的选取1.通过机理分析建立数学模型来确定f(x)；+++++++++++++++

4、+++++++++++++++f=a1+a2xf=a1+a2x+a3x2f=a1+a2x+a3x2f=a1+a2/xf=aebxf=ae-bx2.将数据(xi,yi)i=1,…n作图，通过直观判断确定f(x)：实例讲解某种合成纤维的强度与其拉伸倍数有直接关系，下表是实际测定的24个纤维样品的强度与相应拉伸倍数的记录。提示：将拉伸倍数作为x,强度作为y,在座标纸上标出各点，可以发现什么?数据表格从上图中可以看出强度与拉伸倍数大致成线形关系，可用一条直线来表示两者之间的关系。解：设y*=a+bxi，令δ=yi-y*i=yi-a-bxi，根据最小二乘原理，即使误差的平方和达到最小，也就

5、是令nQ=∑δi2i=1为最小，即求使(a,b)=有最小值的a和b的值。计算出它的正规方程得解得：a=0.15,b=0.859直线方程为：y*=0.15+0.859x曲线拟合问题的提法已知一组（二维）数据，即平面上n个点（xi,yi)i=1,…n,寻求一个函数（曲线）y=f(x),使f(x)在某种准则下与所有数据点最为接近，即曲线拟合得最好。+++++++++xyy=f(x)(xi,yi)ii为点（xi,yi)与曲线y=f(x)的距离多项式的最小二乘拟合的MATLAB函数文件agui_fit.m如下:最小二乘法拟合函数polyfit格式：a=polyfit(x,y,n)