资源描述:
《最新5.曲线拟合-最小二乘法课件ppt.ppt》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、5.曲线拟合-最小二乘法y=f(x)y=p(x)xx0x1x2……xnyy0y1y2……yn插值求一条曲线,使数据点均在离此曲线的上方或下方不远处,所求的曲线称为拟合曲线,它既能反映数据的总体分布,又不至于出现局部较大的波动,更能反映被逼近函数的特性,使求得的逼近函数与已知函数从总体上来说其偏差按某种方法度量达到最小。拟合例某物质的溶解度y和温度x的关系经测定满足下面数据表,试建立该问题的数学模型.将(x,y)的数据点描在一坐标纸上,则如下图所示.7y与x近似成抛物线关系,数据点分布在一抛物线的两侧.从图中可见
2、,因此,可以猜测即有其中是待定常数,这就是本问题的数学模型.8确定了问题的数学模型后,中的待定常数.还需求出不妨假设与待定常数呈线性关系,即这里是线性无关函数系,为待定常数.9在例1中,设函数误差为我们希望猜想的数学模型应尽量接近观测数据,即使得误差带权平方和越小越好.10根据最小二乘原理,应取和使有极小值,故和应满足下列条件:(1)直线拟合设已知数据点,分布大致为一条直线。作拟合直线,该直线不是通过所有的数据点,而是使偏差平方和为最小。即得如下方程组例设有某实验数据如下:12341.361.371.952.2
3、814.09416.84418.47520.963用最小二乘法求以上数据的拟合函数解:把表中所给数据画在坐标纸上,将会看到数据点的分布可以用一条直线来近似地描述,设所求的拟合直线为记x1=1.36,x2=1.37,x3=1.95x4=2.28,y1=14.094,y2=16.844,y3=18.475,y4=20.963则正规方程组为其中将以上数据代入上式正规方程组,得解得即得拟合直线(2)多项式拟合有时所给数据点的分布并不一定近似地呈一条直线,这时仍用直线拟合显然是不合适的,可用多项式拟合。对于给定的一组数据
4、寻求次数不超过n(n<5、得所求的多项式为几种常见的数据拟合情况。图(a)表示数据接近于直线,故宜采用线性函数拟合;图(b)数据分布接近于抛物线。可采拟合;二次多项式拟合;(a)(b)图(c)的数据分布特点是开始曲线上升较快随后逐渐变慢,宜采用双曲线型函数或指数型函数图(d)的数据分布特点是开始曲线下降快,随后逐渐变慢,宜采用或或等数据拟合。(c)(d)例3.13设某实验数据如下:12345600.511.522.52.01.00.90.60.40.3用最小二乘法求拟合曲线解:将已给数据点描在坐标系中下图所示,可以看出这些点接近指数曲线
6、,因而可取指数函数作为拟合函数.对函数两边取对数得.令得则就得到线性模型则正规方程组为其中将以上数据代入上式正规方程组,得解得由得由得于是得到拟合指数函数为(4)超定方程组的最小二乘解设线性方程组Ax=b中,,b是m维已知向量,x是n维解向量,当m>n,即方程组中方程的个数多于未知量的个数时,称此方程组为超定方程组。一般来说,超定方程组无解(此时为矛盾方程组),这时需要寻求方程组的一个“最近似”的解.记,称使,即最小的解为方程组Ax=b的最小二乘解。定理是Ax=b的最小二乘解的充分必要条件为是的解.证明:充分性
7、若存在n维向量,使任取一n维向量,令,则且所以是Ax=b的最小二乘解。必要性:r的第i个分量为,,记由多元函数求极值的必要条件,可得即由线性代数知识知,上式写成矩阵形式为它是关于的线性方程组,也就是我们所说的正规方程组或法方程组。可以证明如果A是列满秩的,则方程组(5.48)存在惟一解例求超定方程组的最小二乘解,并求误差平方和。解:方程组写成矩阵形式为正规方程组为即解得此时误差平方和为我们已经讨论了最小二乘意义下的曲线拟合问题,由于方程比较简单,实际中应用广泛,特别是因为任何连续函数至少在一个较小的邻域内可以用
8、多项式任意逼近,因此用多项式作数据拟合,有它的特殊重要性。从而在许多实际问题中,不论具体函数关系如何,都可用多项式作近似拟合,但用多项式拟合时,当n较大时(n≥7),其法方程的系数矩阵的条件数一般较大,所以往往是病态的,因而给求解工作带来了困难。蜀黍田:造就了我们平凡而充实的人生中国是个面食大国,从古至今,面食一直在老百姓的生活中起到了不可替代的作用。种类繁多,十分美味,驰名中外。面食