数值插值方法.ppt

H3(x)应用四个插值基函数表示。设H3(x)的插值基函数为α0(x)，α1(x)，β0(x)，β1(x)，则其中第61页，共93页，星期日，2025年，2月5日可知x1是α0(x)的二重零点，即可假设由可得第62页，共93页，星期日，2025年，2月5日Lagrange插值基函数同理可得第63页，共93页，星期日，2025年，2月5日将以上结果代入得两个节点的三次Hermite插值公式：第64页，共93页，星期日，2025年，2月5日二、三次Hermite插值的余项定理：设f(x)在区间[a,b]上有定义，f(x)在(a,b)内有4阶导数，H3(x)是满足插值条件(j=0,1)的三次Hermite插值函数，则对任意的x∈[a,b]，H(x)的插值余项为证明：由(i=0,1)第65页，共93页，星期日，2025年，2月5日可知，x0，x1均为R3(x)的二重零点，因此可设其中K(x)待定构造辅助函数i=0,1因此φ(t)至少有5个零点。连续4次使用Rolle定理可得，至少存在一点ξ∈[x0,x1]，使得第66页，共93页，星期日，2025年，2月5日第29页，共93页，星期日，2025年，2月5日5.3分段低次插值高次插值的病态性质：对于一个确定的区间，如果插值节点之间的距离较小，自然插值节点就增多，如果用一个多项式插值，自然次数就会升高，也就是说要用高次多项式插值。但是否次数越高，插值多项式的逼近效果越好呢？20世纪初，Runge就给出了一个等距节点插值多项式不收敛的例子。第30页，共93页，星期日，2025年，2月5日Runge反例:(-5≤x≤5)它在[-5,5]上各阶导数均存在，在该区间上取n+1个等距节点：构造拉格朗日插值多项式为：令则第31页，共93页，星期日，2025年，2月5日20.1379310.759615-0.62168440.066390-0.3568260.42321660.0544630.607879-0.55341680.049651-0.8310170.880668100.0470591.578721-1.531662120.045440-2.7550002.800440140.0443345.332743-5.288409160.043530-10.17386710.217397180.04292020.123671-20.080751200.042440-39.95244939.994889下表列出了n=2,4,…,20的Ln(xn-1/2)和R(xn-1/2)的值：第32页，共93页，星期日，2025年，2月5日从表中可以看出，随着n的增加，R(xn-1/2)的绝对值几乎成倍地增加，这说明当n－∞时Ln在[-5,5]上不收敛。Runge证明了，存在一个常数c≈3.63，使得当|x|≤c时，lim(Ln(x))=f(x)(n－∞)；而当|x|c时，Ln(x)发散。下图给出当n=10时，y=L10(x)及f(x)=1/(1+x2)在[-5,5]上的图形。第33页，共93页，星期日，2025年，2月5日取xk=-5+k计算:f(xk)(k=0,1,…,10)构造L10(x).取:tk=-5+0.05k(k=0,1,…,200),计算:L10(tk)L10(t)f(t)f(x)第34页，共93页，星期日，2025年，2月5日一、分段线性Lagrange插值构造Lagrange线性插值1.分段线性插值的构造设插值节点为xi，函数值为yi，i=0,1,2,……,nhi=xi+1-xi，i=0,1,2,……,n-1，任取两个相邻的节点xk，xk+1，形成一个插值区间[xk，xk+1],k=0,1,2,…,n-1第35页，共93页，星期日，2025年，2月5日显然我们称由上式构成的插值多项式L1(x)为分段线性Lagrange插值多项式。i=0,1,2,…,n第36页，共93页，星期日，2025年，2月5日内插外插外插第37页，共93页，星期日，2025年，2月5日故也称折线插值，如右图：但曲线的光滑性较差，且在节点处有尖点。如果增加节点的数量，减小步长，会改善插值效果。因此则第38页，共93页，星期日，2025年，2月5日由前述余项定理可知，n次Lagrange插值多项式的余项为：2.分段线性插值的误差估计则分段线性插值L1(x)的余项为第39页，共93页，星期日，2025年，2月5日二、分段二次La