计算方法插值方法课堂教学版.pptVIP

* 四．拉格朗日插值 2．抛物插值 l0(x1)=0 l0(x2)=0 l0(x)=c(x-x1)(x-x2) l0(x0)=1 插值基函数 有 位有效数字 * 例3 解： 四．拉格朗日插值 2．抛物插值 4 x0=100，y0=10 x1=121，y1=11 x2=144，y2=12 x=115 代入(7)式 y=10.7228 10.723805… * 定义1 若n次多项式li(x)(i=0,1,…,n)在n+1个节点x0x1…xn上满足条件 就称这n+1个n次多项式l0(x)，l1(x)，…，ln(x)为节点x0,x1,…,xn上的n次插值基函数。 四．拉格朗日插值 3．一般情形 定义1 若n次多项式li(x)(i=0,1,…,n)在n+1个节点x0x1…xn上满足条件 就称这n+1个n次多项式l0(x)，l1(x)，…，ln(x)为节点x0,x1,…,xn上的n次插值基函数。 * 四．拉格朗日插值 3．一般情形 li(x) ? n = - = = j j ? i j i i x x C x l 0 ) ( ) ( Ci(x-x0)…(x-xi-1) (x-xi+1)…(x-xn) ? - = = j ? i j i i i i x x C x l ) ( 1 1 ) ( 每个 li 有 n 个根 x0, …,xi-1, xi+1, …, xn 与 有关，而与 无关 节点 f 拉格朗日插值多项式 * 四．拉格朗日插值 4．算法框图 * 定理3 设x0,x1,…,xn ∈[a,b]，f(x)在[a,b]内有连续的直到n+1阶导数，且f(xi)=yi(i=0,1,…,n)，则当x∈[a,b]时，对于式(9)给出的pn(x)，存在ξ ∈[a,b] ，使得： 五．插值余项 1．拉格朗日余项定理 定理3 设x0,x1,…,xn ∈[a,b]，f(x)在[a,b]内有连续的直到n+1阶导数，且f(xi)=yi(i=0,1,…,n)，则当x∈[a,b]时，对于式(9)给出的pn(x)，存在ξ ∈[a,b] ，使得： * 证明： 五．插值余项 Rolle’s Theorem: 若 充分光滑， ，则 存在 使得 。 推广：若 使得 使得 存在 使得 Rn(x) 至少有 个根 n+1 任意固定 x ? xi (i = 0, …, n), 考察 ? = - ? = n i i x t x K t Rn t 0 ) ( ) ( ) ( ) ( φ ?(t)有 n+2 个不同的根 x0 … xn x ! ) 1 ( ) ( ) ( ) 1 ( + - + n x K R n n ξ 注意这里是对 t 求导 ! ) 1 ( ) ( ) ( ) 1 ( + = + n f x K n ξ = + - ! ) 1 )( ( n x K 1．拉格朗日余项定理 即对于任意的x，都有一个ξ满足余项公式 * 只能在理论上说明ξ存在； 五．插值余项 2．补充说明 若存在 ，可得误差估计： 为插值区间。如果插值点x位于插值区间内，这种插值过程称为内插(interpolation)，否则称为外推（extrapolation )，插值的外推过程是不可靠的； f(x)要足够光滑，否则插值效果不好。当f(x)为任一个次数? n 的多项式时， , 可知 ，即n次插值多项式对于次数? n 的多项式是精确的。 Hn为次数不大于n的多项式的集合。 * 例4 已给sin0.32=0.314567,sin0.34=0.333487,sin0.36=0.352274，用线性插值及抛物插值计算sin0.3367的值并估计截断误差。 五．插值余项 3．应用实例 解：x0=0.32，y0=0.314567，x1=0.34，y1=0.333487，x2=0.36，y2=0.352274，x=0.3367。 取x0和x1： 取x1和x2： 截断误差限： 截断误差限： 线性插值： 外推 内插 * 例4 已给sin0.32=0.314567,sin0.34=0.333487,sin0.36=0.352274，用线性插值及抛物插值计算sin0.3367的值并估计截断误差。 五．插值余项 3．应用实例 解：x0=0.32，y0=0