形式陈述
在区间 [ a , b ] 上固定互异节点 x 0 , … , x n ,令 ℓ j 为相应 Lagrange 基,即在第 j 个节点取 1 、其余节点取 0 的基多项式。插值 理路 多项式插值问题 Polynomial interpolation 由互异节点上的有限数据唯一确定次数受限的插值多项式,并区分对象存在性与具体表示算法。 算子
I n : C [ a , b ] ⟶ P n , ( I n f ) ( x ) = ∑ j = 0 n f ( x j ) ℓ j ( x ) 是从一致范数空间到次数不超过 n 的多项式空间的线性算子。Lebesgue 函数和 Lebesgue 常数定义为
Λ n ( x ) = ∑ j = 0 n | ℓ j ( x ) | , Λ n = max x ∈ [ a , b ] Λ n ( x ) . 在输入、输出都采用无穷范数时,Λ n = ‖ I n ‖ ∞ 。若节点数据受到扰动 | δ y j | ≤ ε ,插值结果的变化满足
| δ p ( x ) | ≤ Λ n ( x ) ε , ‖ δ p ‖ ∞ ≤ Λ n ε . 这个上界能被最坏扰动达到:固定查询点 x ,令每个 δ y j 的符号与 ℓ j ( x ) 相同,所有贡献就同向相加。若从连续函数空间看算子范数,取节点值为 sgn ℓ j ( x ) ,再用分段线性函数连接,就能保持输入一致范数不超过 1 ,而输出在该点达到 Λ n ( x ) 。因此这里得到的是精确的算子范数,不只是一个方便的粗上界。
所以 Λ n 是固定节点插值问题的绝对放大因子 理路 问题条件性与条件数 Conditioning of a problem · Condition number 度量问题真解对输入微扰的局部敏感性,并把这种固有敏感性与算法误差分开。 ,与采用单项式、Newton 或重心公式求值无关。算法还会产生额外舍入误差,但不能改变同一节点集合定义的算子范数。
对任意 q ∈ P n ,利用 I n q = q 可写
f − I n f = ( f − q ) − I n ( f − q ) . 取范数并对 q 取下确界,得到准最优界
‖ f − I n f ‖ ∞ ≤ ( 1 + Λ n ) inf q ∈ P n ‖ f − q ‖ ∞ . 右侧同时包含最佳多项式逼近误差和插值算子放大,两项缺一不可。这里的下确界由最佳一致逼近存在性 理路 最佳一致多项式逼近 Best uniform polynomial approximation · Minimax polynomial approximation · 极小极大多项式逼近 在固定次数内最小化整区间最大误差,证明最优解存在,并区分插值、积分投影与一致最佳解。 保证能够达到;若要证明某个具体候选已达到它,可用交错点证书 理路 Chebyshev 交错定理 Chebyshev alternation theorem · Equioscillation theorem · 等振荡定理 用交替达到最大误差的点认证最佳多项式,完整证明必要、充分、唯一性及近最佳下界。 。
直觉
Lagrange 基在节点上像一组精确的开关:第 j 个基在 x j 为 1 、在其余节点为 0 。离开节点后,这些基可能出现很大的正负摆动;数据中很小的变化乘上这些摆动后,会在区间某处叠成较大偏差。Lebesgue 函数画出这种局部放大,常数则记录最坏位置。
它是问题条件性 理路 问题条件性与条件数 Conditioning of a problem · Condition number 度量问题真解对输入微扰的局部敏感性,并把这种固有敏感性与算法误差分开。 ,不是目标函数平滑度,也不是某段求值程序的评分。重心公式 理路 重心 Lagrange 插值 Barycentric Lagrange interpolation · Barycentric interpolation 预计算重心权后以线性成本求值 Lagrange 插值多项式,并说明节点命中、权重尺度与数值稳定性的条件。 可以更稳定地算出既定插值多项式 理路 多项式插值问题 Polynomial interpolation 由互异节点上的有限数据唯一确定次数受限的插值多项式,并区分对象存在性与具体表示算法。 ,却不会让 Λ n 变小。
例子与边界
取 [ − 1 , 1 ] 上的三个节点 − 1 , 0 , 1 ,其基函数为
ℓ 0 ( x ) = x ( x − 1 ) 2 , ℓ 1 ( x ) = 1 − x 2 , ℓ 2 ( x ) = x ( x + 1 ) 2 . 在 x = 1 / 2 处,三个值为 − 1 / 8 , 3 / 4 , 3 / 8 ,和为 1 ,绝对值之和却为 5 / 4 。因此若三个数据分别扰动 ( − ε , ε , ε ) ,查询值就改变 5 ε / 4 ;若三个数据都增加 ε ,查询值却只增加 ε 。最坏放大取决于扰动方向,不能把上界理解成每次都会发生。进一步在 0 ≤ x ≤ 1 上求得 Λ 2 ( x ) = 1 + x − x 2 ,利用左右对称可知全区间最大值确为 5 / 4 。
随着次数增加,固定区间上的等距节点 Lebesgue 常数会快速、最终呈指数级增长;Chebyshev 型节点的增长仅为对数级。端点附近的基函数摆动是两者差异最明显的区域,因此节点向端点聚集并不是装饰,而是在重新控制插值算子的几何放大。
Runge 函数的完整节点 sweep 与误差数据集中放在Chebyshev 多项式与节点 理路 Chebyshev 多项式与节点 Chebyshev polynomials and nodes · Chebyshev nodes 以余弦表示和端点聚集节点控制区间上一致逼近、插值条件性与高次振荡。 页面。本页只解释机制:等距节点同时面临较差的节点乘积和迅速增长的算子范数,换成稳定求值公式不能改变这两个问题。
小 Λ n 也不自动保证对每个连续函数快速收敛。准最优界中仍有最佳逼近误差;若函数缺乏光滑性或具有靠近区间的奇异性,低次多项式本身可能逼近得慢。Lebesgue 常数只控制“插值比当前最佳多项式差多少倍”。
推论与应用
节点设计应同时考虑目标函数类、误差范数与 Lebesgue 常数。Chebyshev 点在一致范数多项式插值中提供近最优条件性,但换成加权范数、周期问题或特殊采样约束后,合适的节点也可能改变。
实验比较若要说明节点条件性,应对两组节点使用同一种求值实现,并报告相同的一致误差度量;否则公式误差、节点误差和采样网格误差会混在一起。这个原则也适用于噪声数据:必须给出扰动尺度,不能只展示一条看似平滑的曲线。
参考资料
Lloyd N. Trefethen, Approximation Theory and Approximation Practice , extended ed., SIAM, 2019, Chs. 15–16.
Jean-Paul Berrut and Lloyd N. Trefethen, “Barycentric Lagrange Interpolation ,” SIAM Review 46(3), 2004,§9.1,Lebesgue 常数与算子范数。
NIST Digital Library of Mathematical Functions, §3.3 Interpolation , accessed 2026.