[论文解读] Compact Finite Differences and Cubic Splines

Robert M. Corless|arXiv (Cornell University)|May 19, 2018

Advanced Numerical Analysis Techniques参考文献 25被引用 3

一句话总结

本文在均匀网格上建立了四阶紧凑有限差分与三次样条之间的新联系，表明基于Padé的紧凑有限差分格式可产生三次样条所要求的精确导数，以实现二阶连续可微性。本文提出一种新型紧凑三次插值函数，在边界处具有更高阶精度，并证明在均匀或轻微非均匀网格下，其对应的矩阵为完全非负且正定。

ABSTRACT

In this paper I uncover and explain---using contour integrals and residues---a connection between cubic splines and a popular compact finite difference formula. The connection is that on a uniform mesh the simplest Padé scheme for generating fourth-order accurate compact finite differences gives extsl{exactly} the derivatives at the interior nodes needed to guarantee twice-continuous differentiability for cubic splines. %I found this connection surprising, because the two problems being solved are different. I also introduce an apparently new spline-like interpolant that I call a compact cubic interpolant; this is similar to one introduced in 1972 by Swartz and Varga, but has higher order accuracy at the edges. I argue that for mildly nonuniform meshes the compact cubic approach offers some potential advantages, and even for uniform meshes offers a simple way to treat the edge conditions, relieving the user of the burden of deciding to use one of the three standard options: free (natural), complete (clamped), or "not-a-knot" conditions. Finally, I establish that the matrices defining the compact cubic splines (equivalently, the fourth-order compact finite difference formulæ) are positive definite, and in fact totally nonnegative, if all mesh widths are the same sign.

研究动机与目标

• 揭示四阶紧凑有限差分格式与均匀网格上三次样条导数之间此前未被注意到的数学联系。
• 开发一种类似样条的新插值函数——称为紧凑三次插值函数——其在网格边界处的精度优于标准三次样条。
• 提供一个统一的框架以处理边界条件，无需用户在自然、夹紧或非打结条件之间进行选择。
• 证明当网格宽度符号一致时，紧凑三次样条的系统矩阵为正定且完全非负。

提出的方法

• 使用围线积分与留数计算，从分析上推导并解释紧凑有限差分公式与三次样条所需导数值之间的等价性。
• 通过Padé逼近推导紧凑有限差分格式，确保内部节点具有四阶精度。
• 在边界处引入一种改进的四节点紧凑有限差分公式，即使在轻微非均匀网格上也能实现$O(h^4)$精度的导数计算。
• 采用重心形式以实现对所得分段三次插值函数的稳定且高效的求值。
• 证明当所有网格宽度同为正或同为负时，系数矩阵为完全非负且正定。
• 将新提出的紧凑三次插值函数与经典样条及pchip进行比较，强调其在节点处导数精度的提升以及边界处理的简化。

实验结果

研究问题

• RQ1 为何在均匀网格上，基于Padé逼近的四阶紧凑有限差分格式能产生与三次样条完全相同的导数值？
• RQ2 能否构造一种新插值函数，既保持三次样条的光滑性，又在边界处实现更高阶精度？
• RQ3 在非均匀网格下，紧凑三次插值函数的稳定性和条件数与经典样条方法相比如何？
• RQ4 紧凑三次样条的系统矩阵具有哪些代数性质（例如正定性）？
• RQ5 紧凑有限差分与三次样条之间观察到的等价性是否为巧合，还是存在更深层的解析解释？

主要发现

• 在均匀网格上，通过Padé逼近的四阶紧凑有限差分格式所产生的导数值与三次样条所要求的完全一致，证实了精确的数学等价性。
• 所提出的紧凑三次插值函数在端点处的导数精度达到$O(h^4)$，优于经典样条中端点条件的$O(h^3)$精度。
• 当所有网格宽度符号相同时，紧凑三次样条的系统矩阵被证明为正定且完全非负，确保了数值稳定性。
• 通过围线积分与留数方法，解释了紧凑有限差分与三次样条之间的联系，为这一看似意外的数值巧合提供了复分析层面的解释。
• 基于样条分解为纯三次项与$|x - x_0|^3$项的另一种解释也能说明该等价性，但留数法因其更广泛的理论意义而更受青睐。
• 该方法已在MATLAB中实现为`compactcubic.m`，提供了一种实用、稳定且用户友好的传统样条构造替代方案，具备自动边界处理功能。