[論文レビュー] Detailed analysis of prolate quadratures and interpolation formulas

Andrei Osipov, Vladimir Rokhlin|arXiv (Cornell University)|Aug 23, 2012

Electromagnetic Scattering and Analysis参考文献 20被引用数 6

ひとこと要約

本稿では、帯域制限関数のためのプロラート・シュフィロイド・ウェーブ関数（PSWFs）、固有値、およびガウス型積分公式の、非常に効率的な数値アルゴリズムを提示する。帯域制限 $c$ における計算コストを従来の $O(c^2)$ または $O(c^3)$ から $O(c)$ に低減することで、$[-1,1]$ 上での帯域制限関数の高速かつ高精度な統合と補間を、PSWFに基づく積分公式を用いて実現可能にする。

ABSTRACT

As demonstrated by Slepian et. al. in a sequence of classical papers, prolate spheroidal wave functions (PSWFs) provide a natural and efficient tool for computing with bandlimited functions defined on an interval. As a result, PSWFs are becoming increasing popular in various areas in which such function occur - this includes physics (e.g. wave phenomena, fluid dynamics), engineering (e.g. signal processing, filter design), etc. To use PSWFs as a computational tool, one needs fast and accurate numerical algorithms for the evaluation of PSWFs and related quantities, as well as for the construction of quadratures, interpolation formulas, etc. Even though, for the last half a century, substantial progress has been made in design of such algorithms, the complexity of many of the existing algorithms, however, is at least quadratic in the band limit $c$. For example, the evaluation of the $n$th eigenvalue of the prolate integral operator requires at least $O(c^2)$ operations. Therefore, while the existing algorithms are quite satisfactory for moderate values of $c$ (e.g. $c \leq 10^3$), they tend to be relatively slow when $c$ is large (e.g. $c \geq 10^4$). In this paper, we describe several numerical algorithms for the evaluation of PSWFs and related quantities, and design a class of PSWF-based quadratures for the integration of bandlimited functions. Also, we perform detailed analysis of the related properties of PSWFs. While the analysis is somewhat involved, the resulting numerical algorithms are quite simple and efficient in practice. For example, the evaluation of the $n$th eigenvalue of the prolate integral operator requires $O(n+c)$ operations; also, the construction of related accurate quadrature rules requires $O(c)$ operations. Our results are illustrated via several numerical experiments.

研究の動機と目的

• プロラート・シュフィロイド・ウェーブ関数（PSWFs）およびその固有値や積分重みなどの関連量を評価する高速かつ高精度なアルゴリズムの開発。
• バンド制限 $c$ の最初の $n$ 個のPSWFを誤差 $\leq \varepsilon$ で統合できるPSWFに基づく積分公式のクラスの設計。
• 従来のアルゴリズムが $O(c^2)$ または $O(c^3)$ でスケーリングするのを改善し、積分公式構築では $O(c)$、固有値評価では $O(n + c)$ に低減。
• ノード構造、固有値の減衰、ルジャンドル展開係数などのPSWFの性質に関する、詳細な理論的および数値的分析の提供。
• 機械精度に達する正確性と $c$ に線形的にスケーリングする性能を示す、広範な数値実験によるアルゴリズムの妥当性の検証。

提案手法

• 本稿では、漸近展開と再帰関係を用いて、$n$ 番目のPSWF固有値を $O(n + c)$ 操作で計算する新しいアルゴリズムを導出する。
• 積分ノードは $n$ 番目のPSWFの根として構築され、重みは第二種ルジャンドル関数を含む、導関数 $\psi_n'(t_j)$ を用いた新しい公式により計算される。
• 重みは、$O(n)$ 項まで切り捨てたルジャンドル係数 $\alpha_k^{(n)}$ を含む級数展開を用いて計算され、効率性を確保する。
• 本手法は、PSWFの直交性および漸近的性質に加え、$Q_k(t)$ などの特殊関数の高速評価を活用する。
• 対称性を活用して $j > (n+1)/2$ の計算を削減する、二重精度算術を用いた堅牢な数値スキームが実装されている。
• 理論的結果は、機械精度に達する正確性と $c$ に線形的にスケーリングする性能を示す数値実験によって裏付けられる。

実験結果

リサーチクエスチョン

• RQ1 PSWF固有値の評価にかかる計算コストを $O(c^2)$ 未満に低減可能か。その場合、どのようなアルゴリズム的手法が有効か？
• RQ2 バンド制限 $c$ の最初の $n$ 個のPSWFを誤差 $\varepsilon$ 以内に統合するために必要な最小の積分ノード数 $n$ は何か？
• RQ3 PSWFの根と微分値のみを用いて、重みを効率的かつ正確に計算する方法は何か？
• RQ4 ルジャンドル関数展開が、安定的かつ高速な重み計算式を導出する上で果たす役割は何か？
• RQ5 提案されたアルゴリズムは、大きな $c$（例：$c \geq 10^4$）に対しても、正確性と効率性を維持できるか？

主な発見

• $n$ 番目のPSWF固有値の計算が $O(n + c)$ 操作で達成され、従来の方法の $O(c^2)$ に比べ顕著な改善が得られた。
• バンド制限 $c$ の帯域制限関数のための高精度積分公式の構築に必要な計算量は $O(c)$ にまで低減され、従来の $O(c^3)$ に比べ顕著な改善が得られた。
• 数値実験により、新しい公式で計算された積分重みが理論的期待値と機械精度内（残差は $10^{-15}$ のオーダー）で一致することが確認された。
• 積分重みは正であり、中心に向かって単調に増加する傾向を示し、定理73および注意13の理論的期待と整合的である。
• 計算された重みと解析的公式 (673) の差は、機械精度でゼロに達しており、重み計算アルゴリズムの正しさが検証された。
• $c = 40$ および $n = 41$ の場合、アルゴリズムはすべての重みを正確に計算でき、代替重み公式 (674) の誤差は $|\lambda_n| \approx 7 \times 10^{-9}$ のオーダーとなり、注意14の予測通りであった。