第4章_插值法与最小二乘拟合教材教学课件.pptVIP

第4章 插值法与最小二乘拟合;4.1 多项式插值;从几何上看; 用Pn表示所有次数不超过n的多项式函数类,若Pn(x)? Pn ,则; 结论 给定n+1个互异节点x0,x1,…,xn上的函数值y0,y1,…,yn,则满足插值条件(4.2)的n次插值多项式Pn(x)是存在且唯一的.; Rn(x)=K(x)?n+1(x);注意;x 0 y ; y 1 0 xk xk+1 x;y 1 0 x;用二次插值计算ln11.25的近似值,并估计误差.; 称?(xj)-?(xi)与xj-xi(i?j)的比值为?(x)关于节点xi, xj的一阶差商,并记为?[xi ,xj],即; 差商具有以下性质:;; 例 给出函数y=?(x)的函数表; ?(x)=?(x0)+(x-x0)?[x0,x];则有; Nk+1(x)=Nk(x)+?k+1(x)f [x0,x1,…,xk+1] k=1,2,…,n-1;4.1.4 Hermite插值;带导数的两点插值（重要特例);证明思路与不带导数插值类似。;方法：基函数法，待定参数法，重节点差商表;解：基函数法，待定参数法，重节点差商表; 对?(x)=(1+25x2)-1,在区间[-1,1]上取等距节点xi=-1+ih, i=0,1,…,10,h=0.2,作?(x)关于节点xi(i=0,1,…,10)的10次插值多项式L10(x), 如图所示; 1. 分段线性插值;此时,在小区间[xi-1,xi]上有四个插值条件,故能构造一个三次多项式;则,H3(i)(x)就是满足条件;满足插值条件H3(xi)=yi ,H?3(xi)=y?i(i=0,1,…,n)的分段三次多项式H3(x)为;4.3.3 三次样条插值; S(x)在区间[xi-1,xi]上是三次多项式,S(x)=aix3+bix2+cix+di,有4个待定系数,要确定S(x)共有4n个待定系数. ; 若假设S?(xi)=mi ,i=0,1,…,n,利用分段Hermite插值多项式,当x?[xi-1,xi]时,有 ;于是有;由连续性条件S??(xi-0)=S??(xi+0)可得 ;(2)若边界条件为第二类:S ??(x0)=y??0 , S ??(xn)=y??n ,则有 ; (3) 若边界条件为周期性边界条件，由S?(x0+0)=S?(xn-0) ,和S??(x0+0)=S??(xn-0),有m0=mn，和 ;于是有 ;例 设?(0)=1,?(1)=0 ,?(2)= -1 ,?(3)=0 , ??(0)=1, ??(3)=0, 试求?(x)在区间[0,3]的三次样条插值函数S(x).; 三次样条函数S(x)也可以利用在节点处的二阶导数为参数来表示,设S ??(xi)=Mi,i=0,1,…, n,则对x?[xi-1,xi]有 ;于是有;若边界条件为: M0=y??0 ,Mn=y??n ,可得 ;若边界???件为周期性边界条件，由S?(x0+0)=S?(xn-0) ,和S??(x0+0)=S??(xn-0),有 ;于是有 ;于是有 ;三次样条插值函数的误差估计;4.4 最小二乘拟合;拟合的标准 ;线性拟合（简单例子）;当系数行列式不等于0，用克莱姆方法解出方程的唯一解; ?(x)=a0?0(x)+a1?1(x)+…+am?m(x); 就是在函数类?=Span??0(x),?1(x),…,?m(x)?中找一个函数y=?*(x),使误差向量e*的2-范数的平方达到最小值,即 ;引进向量;最小二乘求解;(是否是最优的)？; 由于正规方程组的系数矩阵是对称正定矩阵,可由平方根法或SOR法求得唯一解a0*,a1*,…,am*,于是得到拟合函数:;举例;最小二乘问题中，如何选择数学模型很重要，即如何选取函数空间 ? = span{?0, ?1, ?, ?m} ，通常需要根据物理意义，或所给数据的分布情况来选取合适的数学模型。;例：求下面数据表的二次最小二乘拟合多项式 ;正交多项式拟合;Schmit正交化过程;