Skip to content

模型Model

多项式插值问题

施工中 · 内容逐步修订

Polynomial interpolation

由互异节点上的有限数据唯一确定次数受限的插值多项式,并区分对象存在性与具体表示算法。

形式陈述 ​

给定非负整数 n 及域 F 上两两互异的节点

x0,x1,…,xn

以及数据 y0,…,yn∈F,多项式插值问题要求寻找 p∈F[x],使

deg⁡p≤n,p(xi)=yi(0≤i≤n).

这样的 p 存在且唯一。存在性可由 Lagrange 基

ℓj(x)=∏0≤k≤nk≠jx−xkxj−xk

直接验证:ℓj(xi)=δij,所以 p=∑j=0nyjℓj 满足全部条件。唯一性则只需观察:若 p 与 q 都符合条件,次数不超过 n 的差多项式 p−q 有 n+1 个互异根,只能恒为零。这是多项式对象的结论,不依赖选用哪组基;单项式、Lagrange 与 Newton 表示在这里的任意域上都可用,数值计算中常见的 Chebyshev 基则放在实数或复数系数空间中讨论。

若写成单项式系数 p(x)=∑k=0nakxk,条件等价于 Vandermonde 系统 Va=y,其中 Vik=xik,且

det⁡V=∏0≤i<j≤n(xj−xi)≠0.

在域上,行列式非零等价于方阵可逆,因此再次证明唯一性,却不说明直接求解该系统具有良好数值性质。存在唯一的数学对象、选择何种表示,以及怎样在浮点数中计算,是三个不同问题。

直觉

n+1 个互异节点给次数至多 n 的多项式施加 n+1 个独立线性条件,恰好确定这个 (n+1) 维空间中的一个元素。插值关注“逐点穿过数据”;逼近关注在某个范数下整体接近目标函数;最小二乘允许不穿过所有点,以换取对噪声或超定数据的整体平衡。若指定最大误差范数,最佳一致逼近还给出独立的存在性与最优性证书。

插值先确定函数,再选择保存它的坐标。单项式系数便于代数识别,Newton 形式适合逐点增加数据,重心 Lagrange 形式适合求值;这些表示在精确算术中相同。在实数或复数域上,固定节点并把 yj 改为 yj+Δyj 时,有

Δp(x)=∑j=0nΔyjℓj(x),|Δp(x)|≤maxj|Δyj|∑j=0n|ℓj(x)|.

在节点上只有一个基函数非零,数据扰动不会被额外放大;在节点之间,基函数绝对值之和可能很大。稳定的求值表示能减少计算误差,却不能消除由节点决定的数据敏感性。

例子与边界

在实数域上,三组数据

(−1,1),(0,0),(1,1)

唯一确定次数不超过 2 的多项式 p(x)=x2。这个简单例子中单项式系数很好辨认,但定理只保证对象唯一;节点数增多或尺度悬殊时,Vandermonde 系统即使可逆也可能非常病态。

若两个普通插值节点重合而给出不同函数值,约束彼此矛盾;若给出相同函数值,它们又是重复条件,不能唯一确定原来的次数上限。Hermite 插值通过在重合节点补充导数数据建立另一套独立条件,但那不是把普通差分公式中的零分母直接保留下来。

即使数据没有噪声,精确穿过节点也不保证区间上逼近良好。经典 Runge 函数 f(x)=1/(1+25x2) 在 [−1,1] 的高次等距节点插值会在端点附近产生越来越明显的振荡;问题的插值多项式仍然唯一,失败发生在节点选择和逼近稳定性,而不是存在唯一性。若数据本身还带噪,精确穿过每个观测又可能把噪声当成信号,此时回归、最小二乘或平滑样条对应的是另一种任务目标。

推论与应用

重心 Lagrange 插值、差商与 Newton 形式都计算本页定义的同一多项式。在域上的精确算术中,乘积树插值还可先计算根积导数权重,再自底向上合并,直接输出同一插值多项式的系数;互异节点保证所需差值都可逆。插值余项研究当数据来自函数采样时,插值多项式离原函数多远;Lebesgue 常数研究节点集合怎样放大采样扰动。

有限域上的秘密共享与 Reed–Solomon 译码也使用多项式插值的唯一性,但其误差模型和计算语境不同。本页以实、复数值计算为主要背景,不把有限域的精确代数步骤误写成浮点算法。

解析函数的快速逼近还能量化到具体几何速率。解析函数的谱收敛用复平面最近奇点及 Bernstein 椭圆控制 Chebyshev 系数,并区分级数截断、节点插值和微分误差。

如果目标是对任意连续函数提供可验证的一致逼近,Bernstein 多项式使用等距样本的非负加权平均,并不要求穿过中间节点。例如对 f(x)=x2,它给出 Bnf=x2+x(1−x)/n,内点误差明确非零,却受 1/(4n) 控制。插值条件与全域逼近保证是两个不同的接口。

参考资料
关系图谱28 个相邻概念 · 4 类关系

拖动节点调整位置。

显示关系

显示:依赖

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