[論文レビュー] Bivariate Lagrange interpolation at the node points of Lissajous curves - the degenerate case

Wolfgang Erb|arXiv (Cornell University)|Mar 3, 2015

Advanced Numerical Analysis Techniques参考文献 10被引用数 4

ひとこと要約

本稿は、非退化Lissajous曲線のノード点を一般化する退化Lissajous曲線における2変数ラグランジュ補間の理論的枠組みを確立する。定義された多項式空間内での一意な補間、明示的なラグランジュ基底多項式の導出、およびリーマン定数の対数的増加を示し、非等方的領域における連続関数に対する安定性と収束性を確認する。

ABSTRACT

In this article, we study bivariate polynomial interpolation on the node points of degenerate Lissajous figures. These node points form Chebyshev lattices of rank $1$ and are generalizations of the well-known Padua points. We show that these node points allow unique interpolation in appropriately defined spaces of polynomials and give explicit formulas for the Lagrange basis polynomials. Further, we prove mean and uniform convergence of the interpolating schemes. For the uniform convergence the growth of the Lebesgue constant has to be taken into consideration. It turns out that this growth is of logarithmic nature.

研究の動機と目的

• 既存の非退化Lissajousフレームワークではカバーされていない非退化Lissajous曲線へのPadua点の補間理論の拡張。
• 非退化Lissajous曲線のノード点をランク1のチェビシェフ格子として特徴づけ、既知の点集合と統一的に取り扱えるようにする。
• 補間スキームの解析的収束結果（平均収束および一様収束）を確立し、リーマン定数の挙動に注目する。
• 非等方的長方形領域におけるパラメータ $ p $ の影響を調査し、補間が方向的滑らかさに適合する役割を示す。
• ラグランジュ基底多項式の明示的公式を提供し、滑らかでない関数を含む数値実験により収束性を検証する。

提案手法

• 自然数 $ n, p \in \mathbb{N} $ で $ \gcd(n, n+p) = 1 $ を満たす条件下で、非退化Lissajous曲線 $ \gamma_{n,p}(t) = (\cos(nt), \cos((n+p)t)) $ を定義し、$ t_k = \frac{\pi k}{n(n+p)} $（$ k = 0, \dots, n(n+p) $）でサンプリングすることでノード集合 $ \mathrm{LD}_{n,p} $ を得る。
• ノード集合 $ \mathrm{LD}_{n,p} $ をランク1のチェビシェフ格子として特徴づけ、Padua点のような既知の補間点集合と関連付ける。
• 補間が $ \mathrm{LD}_{n,p} $ で一意に解ける多項式空間 $ \Pi_{n,p}^{2,L} $ を定義し、基本ラグランジュ基底多項式の明示的公式を導出する。
• $ L^r $-平均収束を、リーマン定数の有界性を用いて連続関数に対して証明する。
• リーマン定数 $ \Lambda_{n,p} $ の増加を分析することで一様収束を確立し、対数的増加、具体的には $ \mathcal{O}(\ln^2 n) $ を示し、Padua点の場合と一致することを確認する。
• $ [-1,1]^2 $、$ [0,1]^2 $、$ [0,2] \times [0,1] $ で数値実験を実施し、滑らかさの異なるテスト関数を用い、異なる $ p_n $ 値における性能を比較する。

実験結果

リサーチクエスチョン

• RQ1 非退化Lissajous曲線のノード点において、2変数ラグランジュ補間が一意に成立するか。
• RQ2 これらのノード点で補間が一意に解ける多項式空間の構造は何か。
• RQ3 補間スキームのリーマン定数は漸近的にどのように振る舞い、収束を保証するか。
• RQ4 パラメータ $ p $ の選択が非等方的長方形領域における近似品質にどのように影響するか。
• RQ5 方向的滑らかさを持つ関数に対して、補間スキームが最適収束率を達成できるか。また、$ p $ はその影響をどのように及ぼすか。

主な発見

• 非退化Lissajous曲線のノード点 $ \mathrm{LD}_{n,p} $ は、多項式空間 $ \Pi_{n,p}^{2,L} $ 内で一意な2変数ラグランジュ補間を可能とし、これは $ \Pi_n^2 $ の真部分空間である。
• Padua点の場合と類似した明示的公式がラグランジュ基底多項式に導出され、補間子の効率的計算が可能になる。
• リーマン定数 $ \Lambda_{n,p} $ は $ \mathcal{O}(\ln^2 n) $ のオーダーで増加し、対数的増加を確認することで補間スキームの安定性と収束性が保証される。
• $ L^r $-ノルムにおける平均収束が連続関数に対して証明され、収束速度は関数の方向的滑らかさに依存する。
• 数値実験では、$ [0,1]^2 $ のような等方的領域では $ p_n = 1 $（Padua点）が最良の結果をもたらすが、$ [0,2] \times [0,1] $ のような非等方的領域では $ p_n $ を大きく（例：$ p_n = n+1 $）することで性能が向上する。
• 滑らかさが混合する関数に対しては、$ p_n $ の選択を調整することで方向的モジュラスの連続性をバランスさせられ、特に $ p_n = \lfloor \sqrt{n} \rfloor n+1 $ が1変数で滑らかでない関数に対して最良の結果をもたらす。