← 返回 插值与逼近 › 拉格朗日插值误差公式 定理 已证明
拉格朗日插值误差公式 命题陈述
设 f f f 在包含相异节点 x 0 , x 1 , … , x n x_0, x_1, \dots, x_n x 0 , x 1 , … , x n 和点 x x x 的区间上 ( n + 1 ) (n+1) ( n + 1 ) 次连续可微,P P P 是在这些节点处插值 f f f 的 n n n 次多项式。则存在该区间内的 ξ \xi ξ 使得 f ( x ) − P ( x ) = f ( n + 1 ) ( ξ ) ( n + 1 ) ! ∏ i = 0 n ( x − x i ) f(x) - P(x) = \dfrac{f^{(n+1)}(\xi)}{(n+1)!} \prod_{i=0}^{n} (x - x_i) f ( x ) − P ( x ) = ( n + 1 )! f ( n + 1 ) ( ξ ) ∏ i = 0 n ( x − x i ) 成立。
为什么成立?
插值多项式在节点处与f精确吻合,但对节点之间的f一无所知,因此剩余误差必须在每个节点处为零——这正是乘积项所强制的——并由一个衡量f弯曲程度(超出n次多项式所能捕捉范围)的剩余导数来调节大小。
证明思路 第一步(巧妙的辅助函数)。固定一个不是节点的点 x x x (若是节点,误差平凡地为 0 0 0 )。设 w ( t ) = ∏ i = 0 n ( t − x i ) w(t) = \prod_{i=0}^{n} (t - x_i) w ( t ) = ∏ i = 0 n ( t − x i ) ,并定义常数 c = f ( x ) − P ( x ) w ( x ) c = \dfrac{f(x)-P(x)}{w(x)} c = w ( x ) f ( x ) − P ( x ) (由于 w ( x ) ≠ 0 w(x)\ne0 w ( x ) = 0 故有良好定义)。定义辅助函数 g ( t ) = f ( t ) − P ( t ) − c w ( t ) g(t) = f(t) - P(t) - c\, w(t) g ( t ) = f ( t ) − P ( t ) − c w ( t ) 。
第二步(计数g的根)。在每个节点 x i x_i x i 处,由插值性质 f ( x i ) − P ( x i ) = 0 f(x_i)-P(x_i)=0 f ( x i ) − P ( x i ) = 0 ,且由 w w w 的定义 w ( x i ) = 0 w(x_i)=0 w ( x i ) = 0 ,故在全部 n + 1 n+1 n + 1 个节点处都有 g ( x i ) = 0 g(x_i)=0 g ( x i ) = 0 。此外,由 c c c 的选取,g ( x ) = f ( x ) − P ( x ) − c w ( x ) = f ( x ) − P ( x ) − [ f ( x ) − P ( x ) ] = 0 g(x) = f(x)-P(x) - c\,w(x) = f(x)-P(x) - [f(x)-P(x)] = 0 g ( x ) = f ( x ) − P ( x ) − c w ( x ) = f ( x ) − P ( x ) − [ f ( x ) − P ( x )] = 0 。故 g g g 有 n + 2 n+2 n + 2 个互异的根:n + 1 n+1 n + 1 个节点加上 x x x 本身。
第三步(反复应用罗尔定理)。在 g g g 的每对相邻根之间(n + 2 n+2 n + 2 个根之间共有 n + 1 n+1 n + 1 个这样的间隔),罗尔定理给出一点使 g ′ g' g ′ 为零,故 g ′ g' g ′ 至少有 n + 1 n+1 n + 1 个根。对 g ′ , g ′ ′ , … g', g'', \dots g ′ , g ′′ , … 重复此论证,每求一次导就减少一个根,故经过 n + 1 n+1 n + 1 次应用后,g ( n + 1 ) g^{(n+1)} g ( n + 1 ) 在该区间内至少有一个根 ξ \xi ξ 。
第四步(求导并解出误差)。由于 P P P 的次数不超过 n n n ,其 ( n + 1 ) (n+1) ( n + 1 ) 阶导数为 0 0 0 ;而 w w w 是次数为 n + 1 n+1 n + 1 的首一多项式,故对一切 t t t 都有 w ( n + 1 ) ( t ) = ( n + 1 ) ! w^{(n+1)}(t) = (n+1)! w ( n + 1 ) ( t ) = ( n + 1 )! 。对 g g g 求导得 g ( n + 1 ) ( t ) = f ( n + 1 ) ( t ) − 0 − c ( n + 1 ) ! g^{(n+1)}(t) = f^{(n+1)}(t) - 0 - c\,(n+1)! g ( n + 1 ) ( t ) = f ( n + 1 ) ( t ) − 0 − c ( n + 1 )! ,令 g ( n + 1 ) ( ξ ) = 0 g^{(n+1)}(\xi)=0 g ( n + 1 ) ( ξ ) = 0 得 c = f ( n + 1 ) ( ξ ) ( n + 1 ) ! c = \dfrac{f^{(n+1)}(\xi)}{(n+1)!} c = ( n + 1 )! f ( n + 1 ) ( ξ ) 。回想 c = f ( x ) − P ( x ) w ( x ) c=\dfrac{f(x)-P(x)}{w(x)} c = w ( x ) f ( x ) − P ( x ) 并解出误差,恰好得到 f ( x ) − P ( x ) = f ( n + 1 ) ( ξ ) ( n + 1 ) ! ∏ i = 0 n ( x − x i ) f(x) - P(x) = \dfrac{f^{(n+1)}(\xi)}{(n+1)!} \prod_{i=0}^{n} (x - x_i) f ( x ) − P ( x ) = ( n + 1 )! f ( n + 1 ) ( ξ ) ∏ i = 0 n ( x − x i ) 。