数值分析第6章-2.pptVIP

* 设函数f(x)在[a,b]上有定义,xi,xj?[a,b],且 xi≠xj ,则称 第三节 牛顿插值 拉格朗日插值法最大的弱点是,当插值节点增加后,所有的插值基函数都要重新构造,增加节点前的所有计算结果将毫无作用,造成计算的浪费.牛顿插值却克服了此弱点. * 结束 一 差商 1 差商概念 为f(x)在xi,xj的一阶差商。设f(xi,xj )与f(xj,xk )为f(x)在点xi,xj与点xj,xk的一阶差商,则称 结束 为f(x)在xi, xj , xk的二阶差商。 设f (x0,x1,?,xk-1) 与f (x1,x2,?,xk) 为f(x)在点x0,x1,?,xk-1 与点x1,x2,?,xk的k-1阶差商,则称 为f(x)在点x0,x1,?,xk的k阶差商。 2 差商的基本性质 (1) n次多项式的k阶差商,当k?n时,为一个n-k次多项式,而当kn时恒为零。 * 结束 例如一阶差商 (2) 差商可以表示为节点处函数值的线性组合： 其中?′(xi)=(xi- x0)? (xi- xi-1)(xi- xi+1)?(xi- xk) * 结束 二阶差商 * 结束 类似地 (4) 若f (x)在含x0,x1,?,xk的区间D上k阶可导。则至少存在一点?? D,使得 3 差商表 点x0,x1,x2,x3,?,的所有差商可以列成如下的一张差商表: (3) 差商关于节点对称,与节点的排列次序无关。 例如一阶差商 * 结束 二 牛顿插值公式 设x0,x1 ,?,xn为插值节点,x为插值区间内任意一点,以x,x0,x1 ,?,xn为节点,依次作一阶、二阶,? , n+1阶差商,可得 * 结束 * 依次将后式代入前式得 结束 其中 ? (x)=(x - x0)(x - x1)?(x - xn) * 结束 记 则 其中 n次多项式Nn(x)称为n次牛顿插值多项式,Rn(x)称为插值余项。 ∵ ∴有 由于插值多项式的唯一性,牛顿插值多项式与拉格朗日插值多项式应该是同一多项式,只是写法不同而已。所以其余项也应该相等,即 * 结束 视n为k-1,x为xk,即可得差商的性质（4） 即当插值节点增加后,仅在原有计算结果的基础上增加一些项,增加节点前的所有计算结果均起作用,不会造成计算的浪费。克服了拉格朗日插值的弱点。 因为牛顿插值公式是按x的方次作升幂排列,所以有递推公式 显然,k次牛顿插值公式仅在k-1阶牛顿插值公式后面增加一项(x- x0)?(x- xk-1) f(x0,x1 ,?,xk),作为对k-1阶牛顿插值公式的补偿或修正,提高了插值精度。 * 结束 例6 已知 ,求 ,?=0.001。 解 先构造差商表 三 终止计算原则 * 结束 （ 的准确值是0.6428）,三阶牛顿插值已相当准确。 * 结束 * 6.4 埃尔米特(Hermite)插值 在对函数f(x)进行插值时,有时不仅要求插值多项式在节点处的值等于被插函数在这些点的值,还要求插值多项式的导数在这些点的值也等于函数f(x)的导函数在这些点的值,即带指定导数值的插值,这便是埃尔米特插值. 6.4.1 埃尔米特插值多项式 一般情形为：已知函数f(x)在n+1个互异点xi , (i=0,1,…,n)上对应的函数值yi , (i=0,1,…,n)及导数值y′i , (i=0,1,…,n)求一个2n+1次多项式H2n+1(x),使其满足： 可借鉴拉格朗日插值法构造插值基函数的思想, 令 其中Hi(x)与hi(x)均为2n+1次多项式.且 结束 * 先求Hi (x),由于Hi (xj) =?ij,这一点与拉格朗日插值基函数li(x)相同.所以Hi (x)含有li (x)作为因子,又因为H′i(x) =0,所以li(x)的根都是Hi(x)的二重根,所以Hi(x)含有因子li2(x),显然li2(x)是2n次多项式,因此可令Hi (x) =(ax+b) li2 (x) 由Hi (xi) = 1推得axi+b=1, 结束 * 又Hi′(x) =a li2(x) +2(ax+b) li (x) li′(x) 显然Hi′(xj) =0, (j? i) 而Hi′(xi) =a+2(a xi +b) li′(xi) =0′ 所以a=- li′(xi) , b=1+2 xi li′(xi) 所以 同理可得 这就是带完全的一阶导数的埃尔米特插值公式. 结束 * 当仅两个节点x0,x1时,埃尔米特插值公式为： (6.14)称为两点三次埃尔米特插值,它在后面的样