6 插值.pptVIP

第6章 插值 §6.1 插值多项式 定理 6-1 定理 6-1 §6.2 拉格朗日多项式 例1: n=2：抛物插值 例2: 插值余项 例3: 例4． §6. 3 牛顿插值 ? 差商 差商的性质： 例： 例5： 牛顿插值多项式： 例6: 例7： §6.4 分段低次插值 ? 分段线性插值 §6.5 埃尔米特插值 例: §6.6 三次样条 插值法小结 3 3 0 1 -1 -1/3 -2 -3/2 -1/6 1/24 1 4 7 8 6 1 2 3 4 5 0 1 2 3 4 一阶差商 二阶差商 三阶差商 四阶差商 f(xk) xk k . . 6 , 8 , 7 , 4 , 1 ) ( , 5 , 4 , 3 , 2 , 1 插值多项式 求四次牛顿 时 设当 = = i i x f x 可造插值多项式： /* piecewise polynomial approximation */ 例：在[?5, 5]上考察 的Ln(x)。 - 5 - 4 - 3 - 2 - 1 0 1 2 3 4 5 - 0.5 0 0.5 1 1.5 2 2.5 取 记: 分段低次插值 n 越大，端点 附近抖动越大， 称为Runge 现象 Ln(x) ? f (x) ? /* piecewise linear interpolation */ 在每个区间 上，用1阶多项式 (直线) 逼近 f (x): 记 ， 失去了原函数的光滑性。 易证：当 时， 一致 ? 分段Hermite插值 /* Hermite piecewise polynomials */ 给定 在 上利用两点的 y 及 y’ 构造3次Hermite函数 导数一般不易得到。 不仅要求函数值重合，而且要求若干阶导数也重合。 即：要求插值函数 ? (x) 满足? (xi) = f (xi), ?’ (xi) = f ’ (xi), …, ?(mi) (xi) = f (mi) (xi). 注：? N 个条件可以确定 阶多项式。 N ? 1 ?要求在1个节点 x0 处直到m0 阶导数都重合的插值多项式即为Taylor多项式 其余项为 ?一般只考虑 f 与f ’的值。 设 x0 ? x1 ? x2, 已知 f(x0)、 f(x1)、 f(x2) 和 f ’(x1), 求多项式 P(x) 满足 P(xi) = f (xi)，i = 0, 1, 2，且 P’(x1) = f ’(x1), 并估计误差。 模仿 Lagrange 多项式的思想，设 解：首先，P 的阶数 = 3 ?  + = 2 1 3 ) ( ) ( ) ( ) ( ) ( = 0 i i i x h x1 f ’ x h x f x P ? h0(x) 有根 x1, x2，且 h0’(x1) = 0 ? x1 是重根。 ) ( ) ( ) ( 2 2 1 0 0 x x x x C x h - - = 又: h0(x0) = 1 ? C0 h2(x) h1(x) 有根 x0, x2 ? ) )( )( ( ) ( 2 0 1 x x x x B Ax x h - - + = 由余下条件 h1(x1) = 1 和 h1’(x1) = 0 可解 与h0(x) 完全类似 (x) ? h1 有根 x0, x1, x2 ? ? h1 ) )( )( ( ) ( 2 1 0 1 x x x x x x C x - - - = ? h1 又: ’(x1) = 1 ? C1 可解。 其中 hi(xj) = ?ij , hi’(x1) = 0, (xi) = 0, ’(x1) = 1 ? h1 ? h1 与 Lagrange 分析完全类似 一般地，已知 x0 , …, xn 处有 y0 , …, yn 和 y0’ , …, yn’ ，求 H2n+1(x) 满足 H2n+1(xi) = yi ， H’2n+1(xi) = yi’。 解：设 ? + = n i ) ( ) ( ) ( = 0 i i x h x h yi x H2n+1 ? ? n = 0 i yi’ 其中 hi(xj) = ?ij , hi’(