Skip to content

定理Theorem

多项式插值余项

施工中 · 内容逐步修订

Interpolation remainder · Polynomial interpolation error

在足够光滑条件下以高阶导数和节点乘积表示插值函数误差,并明确该表达式的适用边界。

形式陈述 ​

设 x0,…,xn 是区间 [a,b] 中互异节点,pn 是 f 在这些节点上的插值多项式。设 a<b,并假设 f 是实值函数,具有直到 n+1 阶的连续导数,记为 f∈Cn+1[a,b]。则对每个 x∈[a,b],存在依赖于 x 的点 ξx∈(a,b),使

f(x)−pn(x)=f(n+1)(ξx)(n+1)!ωn+1(x),ωn+1(x)=∏j=0n(x−xj).

当 x 等于节点时,两侧都为零;对非节点,可令

g(t)=f(t)−pn(t)−f(x)−pn(x)ωn+1(x)ωn+1(t).

这里减去的最后一项同时保留节点处的零点,并额外迫使 g(x)=0。将这 n+2 个互异零点从左到右排列,相邻零点之间各有一个导数零点,于是 g′ 至少有 n+1 个零点。继续对各阶导数重复这一过程,也就是反复应用Rolle 定理。最终存在 ξx 使 g(n+1)(ξx)=0;再用 pn(n+1)=0 和 ωn+1(n+1)=(n+1)! 即得公式。

若

Mn+1=maxt∈[a,b]|f(n+1)(t)|,

则有一致上界

‖f−pn‖∞≤Mn+1(n+1)!maxx∈[a,b]|ωn+1(x)|.

这是由光滑性和节点乘积得到的先验上界,不是对每个 x 的精确误差数值。

直觉

插值误差必须在所有节点处消失,所以它带有乘积 ωn+1(x);函数无法被次数 n 多项式捕捉的剩余弯曲,则由某处的 (n+1) 阶导数度量。公式把“目标函数有多弯”与“节点怎样分布”分成两个因子。

点 ξx 由中值定理保证,通常既未知又随查询点变化。余项适合推导界和阶数,不应被当作一个可直接计算的误差修正公式。

例子与边界

先看一个可以把实际误差与上界对照的例子。取 f(x)=x2,在 0,1 处采样得到 0,1,穿过两点的插值直线就是 p1(x)=x。余项中的 f″=2 为常数,因此 f(x)−p1(x)=x(x−1) 完全精确。在 x=1/2 处,直线给出 1/2、真值为 1/4,误差为 −1/4。

一般地,两个节点 a<b 的线性插值满足

f(x)−p1(x)=f″(ξx)2(x−a)(x−b).

若 f 在区间上凸,则 f″≥0,而内部点有 (x−a)(x−b)<0,所以弦线位于函数图像上方。这里余项不仅给大小,也在光滑与符号条件下解释误差方向。又因为 |(x−a)(x−b)| 在中点达到 (b−a)2/4,若 |f″|≤M2,便得到整段的误差界 M2(b−a)2/8。上述二次函数恰好达到此界;对一般函数它只是保证,而非每点都会达到的误差。

若 f 只有连续性,互异节点仍唯一确定插值多项式,但经典导数余项可能毫无意义;例如 f(x)=|x| 在节点 −1,1 上给出常数插值 p1=1,中点误差为 −1,而 f″(0) 不存在,不能套用二阶导数公式。反过来,高阶导数界很小也不能消除采样数据扰动被插值算子放大的问题,那属于节点条件性。

由于 ξx 随 x 变化,不能把

f(n+1)(ξx)(n+1)!ωn+1(x)

当成普通乘积直接对 x 求导。数值微分公式虽可通过对插值多项式求导获得,其导数误差必须另行证明,不能靠形式微分未知的 ξx。

推论与应用

在给定 Mn+1 的情况下,节点选择问题转化为控制 max|ωn+1|。Chebyshev 节点与首一多项式的一致范数极小性质由此进入插值误差分析;它解决的是节点因子,不替代对函数光滑性的检查。

余项还为有限差分和 Newton–Cotes 公式提供局部截断误差来源。应用时必须重新核对求导、积分或复合网格对余项的处理条件,而不是只把原公式中的幂次抄成算法阶数。

参考资料
关系图谱10 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

  1. 前置三跳
  2. 前置二跳
  3. 前置一跳
  4. 当前条目
  5. 后续一跳
  6. 后续二跳
  7. 后续三跳
文字版关系按与当前条目的最短距离分组
类型化关系

使用的工具