[论文解读] Two hierarchies of spline interpolations. Practical algorithms for multivariate higher order splines

Cristian C. Lalescu|ArXiv.org|May 21, 2009

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

一句话总结

本文提出在规则矩形网格上使用泰勒展开和中心差分法导出的多项式基函数，构建了两种多变量高阶样条插值层次结构——Hermite样条与网格样条。提供了网格样条的显式算法，确保单元间导数连续，并针对足够光滑的函数实现对平滑逼近的收敛。

ABSTRACT

A systematic construction of higher order splines using two hierarchies of polynomials is presented. Explicit instructions on how to implement one of these hierarchies are given. The results are limited to interpolations on regular, rectangular grids, but an approach to other types of grids is also discussed.

研究动机与目标

• 为规则矩形网格上的高阶样条插值开发系统化且实用的算法。
• 针对仅在网格节点处已知的函数构建平滑插值函数，避免存储导数带来的高内存开销。
• 通过多项式基函数确保在网格单元之间达到指定阶数的导数连续性。
• 通过局部计算和最小化进程间通信，实现高效的并行化计算。
• 将该方法推广至非矩形网格，同时在共享节点处保持光滑性。

提出的方法

• 采用两种层次结构：基于网格节点处函数值与导数值的Hermite样条，以及仅基于函数值并通过中心差分法构造的网格样条。
• 采用满足特定在网格点处导数条件的多项式基函数 α(n,l)i(x)，以确保 Cm 连续性。
• 利用二项式系数与下降阶乘，推导出基多项式 α(n,l)0(x) 与 α(n,l)1(x) 的显式闭式表达式。
• 将样条计算重新表述为函数值及其导数的线性组合，从而实现高效计算，无需求解大型线性系统。
• 通过一维基函数的张量积将方法推广至 D 维网格，并利用域分解实现并行计算。
• 使用中心差分法近似导数，降低内存需求，并确保单元边界处导数的连续性。

实验结果

研究问题

• RQ1 如何系统地构建多变量函数在规则网格上的高阶样条插值？
• RQ2 在内存受限条件下，如何高效计算高阶样条而不必存储所有导数？
• RQ3 能否推导出一种多项式基，使得样条计算仅依赖于函数值及其有限差分？
• RQ4 所得到的样条层次结构如何在保证高阶精度的同时，确保在网格单元之间具有 Cm 连续性？
• RQ5 当插值阶数增加时，样条逼近的收敛行为如何，特别是对于具有有限傅里叶表示的函数？

主要发现

• 通过使用中心差分法近似导数，网格样条层次结构实现了单元间的 Cm 连续性，从而无需显式存储所有导数。
• 推导出基多项式 α(n,l)0(x) 与 α(n,l)1(x) 的显式闭式表达式，涉及二项式系数以及 x 和 (1−x) 的幂次。
• 基多项式在 x=0 和 x=1 处满足所需的导数条件，确保样条在函数值及至阶数 m=(n−1)/2 的导数处与原函数匹配。
• 误差 ||s(n+1)−s(n)|| 随 n 增大而减小，对于足够光滑的函数趋于零，表明逼近收敛于真实函数。
• 通过将进程间通信限制为每个处理器仅与两个相邻域交换数据，并使用取模运算进行索引，该方法支持高效的并行化。
• 该方法可推广至非矩形网格，前提是共享节点处的泰勒展开保持一致，从而维持光滑性。