[论文解读] On the backward stability of the second barycentric formula for interpolation

Walter F. Mascarenhas|arXiv (Cornell University)|Oct 9, 2013

Mathematical functions and polynomials被引用 11

一句话总结

该论文在 Lebesgue 常数较小的条件下，建立了第二类重心插值公式的后向稳定性。论文为节点和权重扰动引起的后向误差提供了严格的上下界，表明数值计算结果等价于在轻微扰动的节点和函数值下进行的精确计算，误差界按 ϵn²log n（Salzer 权重）或 ϵn log n（数值权重）缩放。

ABSTRACT

We present a new stability analysis for the second barycentric formula for interpolation, showing that this formula is backward stable when the relevant Lebesgue constant is small.

研究动机与目标

• 分析在节点和权重存在实际扰动的情况下，第二类重心插值公式的后向稳定性。
• 通过证明当 Lebesgue 常数较小时稳定性依然成立，解决 Higham（2004）在最坏情况下的不稳定性结论之间的矛盾。
• 为后向误差提供紧致的上下界，同时考虑节点和权重的扰动。
• 证明在实际情形（如使用标准权重计算的切比雪夫节点）下，该公式仍保持后向稳定性。
• 通过扰动理论框架，统一多项式与有理函数重心插值的分析。

提出的方法

• 将第二类重心公式形式化为线性算子 Ix,w，将函数值 y 映射为有理插值结果。
• 将 Lebesgue 常数 Λx−,x+,x,w 定义为算子 Ix,w 的算子范数，用于衡量误差的最大放大程度。
• 引入相对误差度量 δjk、δ−j、δ+j 表示节点扰动，ζk 表示权重扰动，以量化数值不精确性。
• 利用 Higham（2002）提出的误差反向计数技术，通过扰动分析推导出后向误差的上界，同时考虑节点和权重误差。
• 建立扰动节点与原始节点之间的双射 χ，将计算值映射到附近点的精确插值结果。
• 运用后向误差分析理论，证明计算值等于在扰动节点和函数值下精确插值的结果，且相对误差较小。

实验结果

研究问题

• RQ1 在何种条件下，第二类重心公式即使在节点和权重存在扰动时仍保持后向稳定性？
• RQ2 节点和权重的相对误差如何共同影响插值公式的后向误差？
• RQ3 为何在实际应用中（如切比雪夫节点）公式保持稳定，而 Higham（2004）声称在最坏情况下存在不稳定性？
• RQ4 使用第二类切比雪夫点进行拉格朗日插值时，后向误差的精确界是什么？
• RQ5 不同的权重计算策略（Salzer 的闭式表达 vs. 数值计算）如何影响后向误差？

主要发现

• 当 Lebesgue 常数 Λx−,x+,x,w 较小时，第二类重心公式即使在节点和权重存在扰动时仍保持后向稳定性。
• 对于第二类切比雪夫节点，使用 Salzer 权重时后向误差界为 O(ϵn²log n)，使用数值计算权重时为 O(ϵn log n)。
• 对于拉格朗日插值，后向误差界在 log n 因子范围内是紧致的，如第 3 节所给出的下界所示。
• 分析同时考虑了节点和权重的扰动；若忽略其中一项，将导致对真实后向误差的低估。
• 计算值 fl(q(ˆx; ˆx, y, ˆw)) 等于在附近点 x 和扰动函数值 ˜yk = yk(1 + αk) 下的精确插值结果 q(x; x, ˜y, w)，其中 ∥α∥∞ 受 δ 和 Lebesgue 常数的函数有界。
• 节点扰动的界为 max{∥x − ˆx∥∞, |x− − ˆx−|, |x+ − ˆx+|}，确保计算点接近原始点。