曲线拟合的最小二乘法.ppt

3.1问题的提出函数解析式未知,通过实验观测得到的一组数据,即在某个区间[a,b]上给出一系列点的函数值yi=f(xi)3.2.曲线拟合的最小二乘法数据含有误差。节点上的函数值是由实验或观测得到的数据，不可避免地带有测量误差，如果要求所得的近似函数曲线精确无误地通过所有的点(xi,yi),就会使曲线保留着一切测试误差。当个别数据的误差较大时,插值效果显然是不理想的。数据量很大。由实验或观测提供的数据个数往往很多,如果用插值法,势必得到次数较高的插值多项式，这样是不可行的。为此,我们希望从给定的数据(xi,yi)出发,构造一个近似函数,不要求函数完全通过所有的数据点，只要求所得的近似曲线能反映数据的基本趋势，如图3.1所示。与函数插值问题不同,曲线拟合不要求曲线通过所有已知点,而是要求得到的近似函数能反映数据的基本关系。在某种意义上,曲线拟合更有实用价值。函数插值是插值函数P(x)与被插函数f(x)在节点处函数值相同,即而曲线拟合函数不要求严格地通过所有数据点,也就是说拟合函数在xi处的偏差(亦称残差）不都严格地等于零。但是,为了使近似曲线能尽量反映所给数据点的变化趋势,要求按某种度量标准最小。若记向量,即要求向量的某种范数最小,如的1-范数或∞-范数即或一般曲线拟合的最小二乘法的求法例1设有某实验数据如下：12341.361.371.952.2814.09416.84418.47520.963例2设某实验数据如下：123456012345521123解：将已给数据点描在坐标系中，可以看出这些点接近一条抛物线，因此设所求的多项式为（4）可化为线性拟合的非线性拟合对于一个实际的曲线拟合问题，一般先按观测值在直角坐标平面上描出散点图，看一看散点的分布同哪类曲线图形接近，然后选用合适的拟合函数。非线性拟合函数可以通过变量替换转化为线性拟合问题，按线性拟合解出后再还原为原变量所表示的曲线拟合方程。表3-4列举了几类经适当变换后化为线性拟合求解的曲线拟合方程及变换关系表3-4几种常见的数据拟合情况。图(a)数据接近于直线，故宜采用线性函数拟合；图(b)数据分布接近于抛物线可采拟合二次多项式拟合；图(c)：开始曲线上升较快随后逐渐变慢,宜采用双曲线型函数或指数型函数图(d)：开始曲线下降快,随后逐渐变慢,宜采用或或等数据拟合。例3设某实验数据如下:12345600.511.522.52.01.00.90.60.40.3则正规方程组为由得小结插值法中的拉格朗日插值多项式是研究数值微积分与微分方程数值解的重要工具。牛顿插值多项式是拉格朗日插值多项式的变形，具有承袭性，比拉格朗日插值多项式节省计算量。分段低次多项式插值由于具有良好的稳定性与收敛性，且算法简单，便于应用。特别是应用广泛的三次样条插值，不但有较好的稳定性和收敛性，而且具有较好的光滑性，从而满足了许多实际问题的要求。需对样条函数作进一步了解的读者可参阅有关文献曲线拟合的最小二乘法是处理实验数据的常用方法。本章主要介绍了最小二乘法的基本原理和线性最小二乘问题的求解方法。多项式拟合是线性最小二乘拟合问题的一种特殊情况,其特点是拟合多项式形式简单,但当n较大时,法方程组往往是病态的。用正交多项式进行曲线拟合，避免了法方程组病态所造成的麻烦。关于非线性最小二乘曲线拟合问题,一般求解比较困难,但对一些特殊情形,可以转换为线性最小二乘拟合问题。作业P89习题三：1Thankyouverymuch!