[论文解读] Cubature formula and interpolation on the cubic domain

Huiyuan Li, Jiachang Sun+1 more|ArXiv.org|Aug 1, 2008

Mathematical functions and polynomials参考文献 19被引用 4

一句话总结

本文提出了一种基于离散傅里叶分析与格点密铺的新型三重积分公式，用于三维立方体 $[-1,1]^3$，实现了度数为 $2n-1$ 的公式仅需 $n^3/4 + \mathcal{O}(n^2)$ 个节点，显著少于乘积高斯积分公式所需的 $n^3$ 个节点。该方法推导出基本插值多项式的紧凑公式，并证明了Lebesgue常数保持在 $\mathcal{O}((\log n)^3)$，从而实现了高维多项式插值的稳定性。

ABSTRACT

Several cubature formulas on the cubic domains are derived using the discrete Fourier analysis associated with lattice tiling, as developed in \cite{LSX}. The main results consist of a new derivation of the Gaussian type cubature for the product Chebyshev weight functions and associated interpolation polynomials on $[-1,1]^2$, as well as new results on $[-1,1]^3$. In particular, compact formulas for the fundamental interpolation polynomials are derived, based on $n^3/4 +\CO(n^2)$ nodes of a cubature formula on $[-1,1]^3$.

研究动机与目标

• 开发适用于立方体域 $[-1,1]^3$ 的高效积分公式，其节点数显著少于标准乘积高斯规则。
• 提出一种基于离散傅里叶分析与格点密铺的新方法，用于构造类高斯积分公式，避免依赖正交多项式根的求解。
• 基于积分节点，推导出基本插值多项式的紧凑显式公式。
• 分析所得到插值过程的Lebesgue常数，并确认其在维度上呈对数增长。

提出的方法

• 利用与格点密铺相关的离散傅里叶分析，使用两种不同的格点：一种定义立方体域（$\mathbb{Z}^3$），另一种确定三角函数空间。
• 通过菱形十二面体作为谱集，在 $[-1,1]^3$ 上构造一个度数为 $2n-1$ 的积分公式，得到 $n^3/4 + \mathcal{O}(n^2)$ 个节点。
• 通过变量替换 $s = \cos 2\pi x$，将环面上的三角函数插值转化为 $[-1,1]^3$ 上的代数多项式插值。
• 定义插值算子 $\mathcal{L}_n f(x) = \sum_{k \in \Xi_n} f(k/(2n)) \ell_k(x)$，其中 $\ell_k(x)$ 由对偶函数 $\Phi_n^*$ 导出。
• 通过 $\mathcal{P}f(x) = \frac{1}{8} \sum_{\varepsilon \in \{-1,1\}^3} f(\varepsilon_1 x_1, \varepsilon_2 x_2, \varepsilon_3 x_3)$ 的对称性平均，确保函数为偶函数并实现节点插值。
• 将插值空间 $\Pi_n^*$ 定义为单项式 $s_1^{k_1}s_2^{k_2}s_3^{k_3}$ 的张成空间，其中 $k_i + k_j \leq n$，对应菱形十二面体的正卦限。

实验结果

研究问题

• RQ1 能否使用超越正交多项式根求解的方法，在 $[-1,1]^3$ 上构造一个度数为 $2n-1$ 的积分公式，且其节点数显著少于 $n^3$？
• RQ2 此类积分公式的节点所对应的插值空间结构为何？
• RQ3 能否为三元变量的基底插值多项式推导出紧凑且显式的公式？
• RQ4 该插值方案的Lebesgue常数的增长率如何？
• RQ5 是否可能在三维立方体上实现类高斯积分公式的理论最低节点数 $n^3/6 + \mathcal{O}(n^2)$？

主要发现

• 本文构造了一个度数为 $2n-1$ 的积分公式，应用于 $[-1,1]^3$，仅使用 $n^3/4 + \mathcal{O}(n^2)$ 个节点，优于乘积高斯积分公式的 $n^3$ 个节点。
• 所构造节点的插值空间被识别为 $\Pi_n^* = \mathrm{span}\{s_1^{k_1}s_2^{k_2}s_3^{k_3} : k_i + k_j \leq n\}$，对应菱形十二面体的正卦限。
• 从对偶函数 $\Phi_n^*$ 推导出基本插值多项式 $\ell_k^*(s)$ 的紧凑公式，实现了高效计算。
• 插值过程的Lebesgue常数被控制在 $\mathcal{O}((\log n)^3)$ 以内，确保了对连续函数的稳定性和收敛性。
• 该方法将二维中基于格点密铺与离散傅里叶分析的方法推广至三维，为高维积分与插值提供了新框架。
• 数值测试验证了 $n^3/4 + \mathcal{O}(n^2)$ 的节点数量与 $\mathcal{O}((\log n)^3)$ 的Lebesgue常数，证实了理论边界的正确性。