[論文レビュー] A new and efficient method for the computation of Legendre coefficients

Enrico De Micheli, G. A. Viano|arXiv (Cornell University)|Jun 2, 2011

Digital Filter Design and Implementation参考文献 7被引用数 4

ひとこと要約

本稿では、関数 f(x) のアベル型積分変換のフーリエ係数を計算する問題に変換することで、Legendre係数を計算するための新規な O(N log N) アルゴリズムを提示する。Dirichlet–Murphy 積分表現を用い、1回の高速フーリエ変換（FFT）を適用することで、標準的な O(N²) の数値積分法よりも著しく高速かつ高精度な性能を達成する。

ABSTRACT

An efficient procedure for the computation of the coefficients of Legendre expansions is here presented. We prove that the Legendre coefficients associated with a function f(x) can be represented as the Fourier coefficients of an Abel-type transform of f(x). The computation of N Legendre coefficients can then be performed in O(N log N) operations with a single Fast Fourier Transform of the Abel-type transform of f(x).

研究の動機と目的

• 近似理論、偏微分方程式（PDEs）、スペクトル法において不可欠な Legendre係数を、より高速かつ高精度に計算するアルゴリズムの開発。
• 急速に振動する被積分関数を伴う高インデックス Legendre係数の計算において、標準的 quadrature 法の O(N²) の計算コストを克服すること。
• アベル型積分変換を通じて Legendre 展開とフーリエ級数の関係を活用し、効率的な計算を可能にすること。
• 標準的 quadrature 法でアベル変換を評価した後、1回の高速フーリエ変換（FFT）を用いることで、O(N log N) の計算量を達成すること。

提案手法

• Legendre多項式の Dirichlet–Murphy 積分表現を用い、Legendre係数を f(x) のアベル型変換のフーリエ係数に関連付ける。
• f(x) と平方根特異性を含む積分を用いて変換関数 f̂(y) を定義し、弱い特異性に対処するために非線形変数変換を適用する。
• 係数 cₙ は cₙ = (2n+1) × aₙ として表され、ここで aₙ は [−π, π] 上での f̂(y) の非負のフーリエ係数である。
• フーリエ係数 aₙ は1回の高速フーリエ変換（FFT）により効率的に計算され、全体の計算量が O(N log N) に削減される。
• アベル型積分は、Gauss-Legendre quadrature と非線形変数変換を用いて、代数的特異性を適切に扱うことで高精度に評価される。
• 実装はGNU科学ライブラリ（GSL）を用いた倍精度で行われ、高精度な基準値との照合による検証が実施された。

実験結果

リサーチクエスチョン

• RQ1 スペクトル法を用いて、O(N²) よりも高速に Legendre係数を計算できるか？
• RQ2 Legendre展開とフーリエ級数の間に、より高速な計算を可能にする関係が存在するか？
• RQ3 f(x) のアベル型積分変換を高精度かつ効率的に計算できるか？ これにより O(N log N) のアルゴリズムが実現可能か？
• RQ4 提案手法は、標準的 quadrature 法および既存の O(N(log N)²) アルゴリズムと比較して、精度と速度の両面で優れているか？

主な発見

• 本手法は O(N log N) の計算量を達成し、標準的 O(N²) の quadrature 法と比べて顕著な改善を示す。
• f(x) = |x|³/² の場合、計算された Legendre係数は基準値と 10⁻¹⁰ から 10⁻¹³ の誤差範囲で一致し、高精度であることが確認された。
• アルゴリズムは N / log N の割合に比例した性能向上比を示し、N が大きくなるほどその重要性が増す。
• 正確なアベル変換の評価後に1回の FFT を適用することで、速度と精度の両方を実現でき、多項式、指数関数、有理関数を含む複数のテスト関数で検証された。
• 本手法は一般化可能でスケーラブルであり、多次元 Legendre 変換や Gegenbauer 多項式展開への応用が有望である。