Ch2_3Hermite插值解说.ppt

§2.6 埃尔米特插值 /* Hermite Interpolation */ §5.2 Hermite Interpolation 一般地，已知 x0 , …, xn 处有 y0 , …, yn 和 y0’ , …, yn’ ，求 H2n+1(x) 满足 H2n+1(xi) = yi ， H’2n+1(xi) = yi’。 例：设 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 分析完全类似 ? + = n i ) ( ) ( ) ( = 0 i i x ? x ? yi x H2n+1 ? n = 0 i yi’ 其中 ? i(xj) = ?ij , ?i’(xj) = 0, (xj) = 0, ’(xj) = ?ij ? ? 解：设 ? + = n i ) ( ) ( ) ( = 0 i i x ? x ? yi x H2n+1 ? n = 0 i yi’ 其中 ? i(xj) = ?ij , ?i’(xj) = 0, (xj) = 0, ’(xj) = ?ij ?i ?i(x) 有根 x0 , …, xi , …, xn且都是2重根 ? ? 由余下条件 ?i (xi) = 1 和 ?i’(xi) = 0 可解Ai 和 Bi ? (x) ?i 有根 x0 , …, xn, 除了xi 外都是2重根 ? ?i ) ( ) ( i i li2(x) x x C x - = 又 ?’i (xi) = 1 ? Ci = 1 ?i ) ( x ) ( i li2(x) x x - = 设 则 这样的Hermite 插值唯一 ?i ) ( ) ( ) ( 2 x l B x A x i i i + = ?i 类似的， Th. 设f(x)?C 2n+2[a,b],则?? ? [a,b], s.t. 满足下面插值条件 ? 求Hermite多项式的基本步骤： ? 写出相应于条件的?i(x)、 ?i(x) 的组合形式； ? 对每一个?i(x)、 ?i(x)找出尽可能多的条件给出的根； ? 根据多项式的总阶数和根的个数写出表达式； ? 根据尚未利用的条件解出表达式中的待定系数； ? 最后完整写出H2n+1(x)。 注：待定系数法仍适用，但插值节点多时比较麻烦。 * * * *