[論文レビュー] 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

ひとこと要約

本稿は、離散フーリエ解析と格子タイル張りを用いて、3次元立方体 $[-1,1]^3$ に対する新しい立方則を導入し、次数 $2n-1$ の公式で $n^3/4 + \mathcal{O}(n^2)$ 個のノードを達成した。これは積型ガウス立方則の $n^3$ 個のノードよりも顕著に少ない。基本的補間多項式の簡潔な公式を導出し、リーマン定数が $\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$ に対して、標準的な積型ガウス則よりもノード数を減らした効率的な立方則の構築を目的とする。
• 直交多項式の根を求める手法に依存せずに、離散フーリエ解析と格子タイル張りを用いてガウス型立方則を構築する新しい手法の確立。
• 立方則のノードに基づく基本的補間多項式の簡潔で明示的な公式の導出。
• 得られた補間プロセスのリーマン定数の解析を行い、次元に応じて対数的に増加することの確認。

提案手法

• 格子タイル張りに関連する離散フーリエ解析を活用し、2つの異なる格子を用いる：1つは立方体領域を定義する（$\mathbb{Z}^3$）、もう1つは三角関数空間を決定する。
• 菱形12面体をスペクトル集合として用い、$[-1,1]^3$ 上の次数 $2n-1$ の立方則を構築し、$n^3/4 + \mathcal{O}(n^2)$ 個のノードを達成した。
• トーラス上の三角関数補間を $[-1,1]^3$ 上の代数的多項式補間に変換するために、変数変換 $s = \cos 2\pi x$ を適用した。
• 補間作用素 $\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^*$ を、$k_i + k_j \leq n$ を満たす単項式 $s_1^{k_1}s_2^{k_2}s_3^{k_3}$ の線形包として定義し、これは菱形12面体の正の象限に対応する。

実験結果

リサーチクエスチョン

• RQ1 $[-1,1]^3$ 上の次数 $2n-1$ の立方則を、直交多項式の根を求める手法に依存せずに、$n^3$ より顕著に少ないノード数で構築することは可能か？
• RQ2 このような立方則のノードに対応する補間空間の構造は何か？
• RQ3 3変数における基本的補間多項式の簡潔で明示的な公式を導出することは可能か？
• RQ4 この補間スキームのリーマン定数の増加率は何か？
• RQ5 ガウス型立方則の理論的下界 $n^3/6 + \mathcal{O}(n^2)$ 個のノードを3次元立方体で達成することは可能か？

主な発見

• 本稿は、$n^3/4 + \mathcal{O}(n^2)$ 個のノードを用いて、$[-1,1]^3$ 上の次数 $2n-1$ の立方則を構築した。これは積型ガウス立方則の $n^3$ 個のノードよりも優れている。
• 構築されたノードの補間空間は、$\Pi_n^* = \mathrm{span}\{s_1^{k_1}s_2^{k_2}s_3^{k_3} : k_i + k_j \leq n\}$ として特定され、これは菱形12面体の正の象限に対応する。
• 双対関数 $\Phi_n^*$ から導かれた $\ell_k^*(s)$ のための簡潔な公式を導出し、効率的な評価が可能になった。
• 補間プロセスのリーマン定数は $\mathcal{O}((\log n)^3)$ で有界であり、連続関数に対して安定で収束することが保証された。
• 2次元でのアプローチを、格子タイル張りと離散フーリエ解析を用いて3次元に一般化し、高次元立方則と補間のための新しい枠組みを提供した。
• 数値的検証により、$n^3/4 + \mathcal{O}(n^2)$ 個のノード数と $\mathcal{O}((\log n)^3)$ のリーマン定数が理論的境界と整合していることが確認された。