第五章曲线拟合与最小二乘法课件.pptVIP

第五章 曲线拟合与最小二乘法;为此,我们希望从给定的数据(xi,yi)出发,构造一个近似函数 ,不要求函数 完全通过所有的数据点，只要求所得的近似曲线能反映数据的基本趋势，如图5-1所示。 ; 与函数插值问题不同,曲线拟合不要求曲线通过所有已知点,而是要求得到的近似函数能反映数据的基本关系。在某种意义上,曲线拟合更有实用价值。 在对给出的实验(或观测)数据 作曲线拟合时,怎样才算拟合得最好呢？一般希望各实验(或观测)数据与拟合曲线的偏差的平方和最小,这就是最小二??原理。 两种逼近概念: 插值: 在节点处函数值相同. 拟合: 在数据点处误差平方和最小; 函数插值是插值函数P(x)与被插函数f(x)在节处函数值相同,即 而曲线拟合函数 不要求严格地通过所有数据点 ,也就是说拟合函数 在xi 处的偏差(亦称残差） 不都严格地等于零。但是,为了使近似曲线能尽量反 映所给数据点的变化趋势,要求 按某种度量标准 最小。若记向量 ,即要求向量 的某种范数 最小,如 的1-范数 或∞-范数 即 ;或 ;??? （1）直线拟合 设已知数据点 ,分布大致为一条直线。作拟合直线 ,该直线不是通过所有的数据点 ,而是使偏差平方和;即得如下正规方程组 ;拟合直线为 记x1=1.36, x2=1.37, x3 =1.95 x4 =2.28, y1 =14.094, y2= 16.844, y3=18.475, y4=20.963 则正规方程组为 ;（2）多项式拟合 有时所给数据点的分布并不一定近似地呈一条直线,这时仍用直线拟合显然是不合适的,可用多项式拟合。对于给定的一组数据 寻求次数不超过m (mN ) 的多项式， ;由于Q可以看作是关于 ( j=0,1,2,…, m)的多元函数, 故上述拟合多项式的构造问题可归结为多元函数的极值问题。令;这是关于系数 的线性方程组，通常称为正规方程组。可以证明，正规方程组有惟一解。 ;解：将已给数据点描在坐标系中，可以看出这些点 接近一条抛物线，因此设所求的多项式为 ;（3）可化为线性拟合的非线性拟合 有些非线性拟合曲线可以通过适当的变量替换转化为线性曲线，从而用线性拟合进行处理，对于一个实际的曲线拟合问题，一般先按观测值在直角坐标平面上描出散点图，看一看散点的分布同哪类曲线图形接近，然后选用相接近的曲线拟合方程。再通过适当的变量替换转化为线性拟合问题，按线性拟合解出后再还原为原变量所表示的曲线拟合方程。 表5-1列举了几类经适当变换后化为线性拟合求解的曲线拟合方程及变换关系 ;表5-1;几种常见的数据拟合情况。图 ( a ) 表示数据接近于 直线，故宜采用线性函数 拟合；图(b)数 据分布接近于抛物线。可采拟合；二次多项式 拟合； ;图 ( c ) 的数据分布特点是开始曲线上升较快随后逐 渐变慢,宜采用双曲线型函数 或指数型函 数 图 ( d ) 的数据分布特点是开始曲线下降快,随后逐渐变慢,宜采用 或 或 等数据拟合。;例5. 3 设某实验数据如下: 1 2 3 4 5 6 0 0.5 1 1.5 2 2.5 2.0 1.0 0.9 0.6 0.4 0.3;则正规方程组为 ; （4）超定方程组的最小二乘解 设线性方程组Ax=b中， ,b 是m维已知向量，x是n维解向量，当m＞n，即方程组中方程的个数多于未知量的个数时，称此方程组为超定方程组。一般来说，超定方程组无解（此时为矛盾方程组),这时需要寻求方程组的一个“最近似” 的解.记