插值.pptVIP

编程框图如图2.3.3，可用二重循环来完成 值的计算，先通过内循环，即先固定 ，令 从0到 累乘求得 然后再通过外循环，即令 从0到 ，累加得出插值结果 2、由误差公式有 差分与函数值的关系 等距节点插值多项式 其截断误差得 其中 ，因 f(x)=sinx，f”(x)= -sinx， 可取 ，于是 ?R1(0.3367)?=?sin 0.3367 –L1(0.3367)? ?1/2(0.3335)(0.0167)(0.0033)?0.92?10–5， 若取x1=0.34，x2=0.36为节点，则线性插值为 其截断误差为 其中 于是 用抛物插值计算 sin0.3367时，可得 这个结果与六位有效数字的正弦函数表完全一样， 这说明查表时用二次插值精度已相当高了。 其截断误差得 其中 于是 第三节 Newton插值 本节介绍Newton插值多项式，该公式是另一种具有承袭性的插值多项式，且在后续章节有重要应用。为了使Newton插值多项式具有承袭性，令 式中c0,c1,…,cn 为插值多项式系数。 为便于表示Nn(x), 引出差商概念. 差商及其性质 定义1　给定一个函数表 一般的, f(x)关于xi,xi+1,…,xi+k的k 阶差商记做 f[xi,xi+1,…,xi+k ] 例： 定理1: 差商具有如下性质 (1)差商与函数值的关系为 (2)差商与结点排列顺序无关 由插值条件 Nn(xi)=f(xi) i=0,1,…,n 导出 Nn(x0) =c0=f(x0) Newton 插值公式 依次类推，得： cn=f[x0,x1,…,xn] 因此，每增加一个结点，Newton插值多项式只增加一项，克服了 Lagrange插值的缺点。 差商表 定理2：Newton插值多项式的余项为 Rn(x)= f(x)-Nn(x)=f[x,x0,x1,…,xn] ?n+1（x） 其中?n+1 (x）=(x-x0)(x-x1)(x-x2)…(x-xn) 证明： 过n+1个点 x0 ,x1 ,…, xn 作f(x)的n次Newton插值多项式Nn(t)，设x与x0 ,x1 ,…, xn 互异, 再取n+2个节点 x ,x0 ,x1 ,… xn, 构造f(x)的n+1次Newton插值多项式 Nn+1(t)=Nn(t)+f[x0 ,x1 ,… ,xn,x]wn+1(t) (a) 满足 Nn+1(xi)=f(xi) (b) Nn+1(x)=f(x) (c) 利用Newton插值余项，可证明差商的性质3） 事实上，由Newton插值余项与Lagrange插值 余项的等价性，有 消掉 ，即得性质(3）。 由(a), (c)有 f(x)=Nn(x)+ f[x0 x1 … xn ,x]wn+1(x) Rn(x)=f(x)-Nn(x)= f[x0 x1 …xn,x]wn+1(x) 例：设f(x)=10x3-100x+1, 计算f[x0 , x1 , x2 , x3]， f[x0 , x1 , x2 , x3 , x4] ? Newton插值实际计算过程为 f ( x0 ) f ( x1 ) f ( x2 ) … f (xn?1) f ( xn ) f [x0, …, xn] f [x0 , x1] f [x1, x2] … f [xn-1, xn] f [x0, x1, x2] … f [xn?2, xn?1, xn] x0 x1 x2 … xn?1 xn 差商表 2 6 9 2 2 3 1 1 0 二阶均差 一阶均差 y