[論文レビュー] The Fourier Transforms of the Chebyshev and Legendre Polynomials

A. S. Fokas, S. A. Smitheman|arXiv (Cornell University)|Nov 21, 2012

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

ひとこと要約

本稿では、チェビシェフ多項式およびルジャンドル多項式の有限フーリエ変換に対する明示的な解析的表現を導出し、半整数階ベッセル関数の新たな表現を可能にする。この手法は、凸多角形における偏微分方程式を解くための統一変換法に応用され、フーリエ空間における最適化された配置点を用いたルジャンドル多項式を用いることで、数値的安定性が向上し、スペクトル精度を達成することが示された。

ABSTRACT

Analytic expressions for the Fourier transforms of the Chebyshev and Legendre polynomials are derived, and the latter is used to find a new representation for the half-order Bessel functions. The numerical implementation of the so-called unified method in the interior of a convex polygon provides an example of the applicability of these analytic expressions.

研究の動機と目的

• チェビシェフ多項式およびルジャンドル多項式の有限フーリエ変換に対する閉形式の解析的表現を導出すること。
• これらの変換を、凸多角形における線形楕円型PDEを解くための統一法に応用すること。
• 特定の配置点を用いたルジャンドル多項式を用いることで、他の基底と比較してより良好な条件数を持つ線形系が得られることを示すこと。
• ルジャンドル多項式変換を用いて、半整数階ベッセル関数の新たな明示的表現を提供すること。
• 正方形領域における数値的検証を通じて、スペクトル精度と行列の条件数の向上を示すこと。

提案手法

• 変数変換 $ x = \cos w $ を用いて、チェビシェフ多項式の有限フーリエ変換を導出し、積分を三角関数形式に変換する。
• 補助関数 $ K_m(z) $ の再帰的関係を確立し、フーリエ変換係数の再帰的計算を可能にする。
• 部分積分と三角恒等式を用いて、一般形の変換を $ i\lambda $ の逆数の累乗の和として導出し、係数 $ \alpha_n^m $ を二項係数および積項で定義する。
• 統一法におけるグローバル関係にルジャンドル多項式変換を適用し、未知の境界値をルジャンドル関数の展開として表現する。
• 先行研究で提案されたように、複素 $ \lambda $-平面に沿った線に沿った配置点を選択し、半ブロック巡回構造を維持するとともに、条件数の改善を図る。
• 過剰決定線形系を $ N $、$ 3N/2 $、$ 2N $ 個の配置点を用いて数値的に解き、行と列をそれぞれ $ l^1 $-ノルムで正規化することで安定性を維持する。

実験結果

リサーチクエスチョン

• RQ1 チェビシェフ多項式およびルジャンドル多項式の有限フーリエ変換に対して、明示的な解析的表現を導出できるか？
• RQ2 これらの変換を用いて、半整数階ベッセル関数の新たな表現を構築できるか？
• RQ3 統一変換法を用いて凸多角形におけるPDEを解く際、ルジャンドル多項式を他の基底と比較して使用した場合の影響は何か？
• RQ4 異なる配置点の選択が、得られる線形系の条件数および精度に与える影響は何か？
• RQ5 ルジャンドル多項式と最適化された配置点を用いることで、PDEの数値解法においてスペクトル精度が達成できるか？

主な発見

• ルジャンドル多項式の有限フーリエ変換により、半整数階ベッセル関数の新たな明示的表現が得られた。
• チェビシェフ多項式の変換は、再帰的部分積分により導出された係数 $ \alpha_n^m $ を持つ $ \lambda^{-n} $ 項の有限和として表現される。
• 統一法においてルジャンドル多項式を基底関数として用いることで、フーリエ基底や他の基底と比較して、より良好な条件数を持つ線形系が得られた。
• 数値実験の結果、未知数 $ N $ 個に対して配置点を $ 2N $ 個に増やす（過剰決定を2倍にすること）ことで、優れた行列の条件数が保証された。
• 本手法はスペクトル精度を達成しており、相対誤差 $ E_\infty $ が $ N $ の増加に従い急速に減少することが、図2のプロットで確認された。
• 修正されたヘルムホルツ方程式に対しても、ラプラス方程式と同様に、配置行列の半ブロック巡回構造が保持され、効率的な解法戦略が可能となった。