[論文レビュー] New method of computing prolate spheroidal wavefunctions and bandlimited extrapolation: A general approach to bandlimited Fredholm kernels

Vishal Vaibhav|arXiv (Cornell University)|Apr 8, 2018

Numerical methods in inverse problems被引用数 2

ひとこと要約

本稿では、球ベッセル関数を用いた新しい退化核法を導入し、帯域制限付きフリードホルム積分方程式を解くために、準楕円体付随関数（PSWFs）を計算する。核関数を球ベッセル関数の基底関数で離散化することにより、関連する固有値問題および帯域制限付き関数の逆問題の効率的解法が可能となり、実数直線上で高い精度を達成する。この手法は、一般の帯域制限付き、平方可積分核に対しても拡張可能である。

ABSTRACT

The paper deals with numerical solution of the Fredholm integral equation associated with the classical problem of extrapolating bandlimited functions known on $(-1,1)$ to the entire real line. The approach presented can be characterized as the degenerate kernel method using the spherical Bessel functions as basis functions. This discretization also facilitates the solution of the associated eigenvalue problem whose eigenfunctions are the prolate spheroidal wave functions of order zero, thus yielding a new method of computing these functions on the entire real line. These ideas are then extended to Fredholm integral equations whose kernel belongs to a class of bandlimited functions that are square integrable. The proposed discretization scheme is used to solve the associated eigenvalue problem as well as the inverse problem that arises in the estimation of object function from its image function.

研究の動機と目的

• 全実数直線上で、零次の準楕円体付随関数を数値的に安定して計算するための手法を開発すること。
• 区間 $(-1,1)$ を超えて帯域制限付き関数の外挿に現れるフリードホルム積分方程式を解くこと。
• 固有値問題および逆問題を解くために、一般の帯域制限付き、平方可積分核に対してもこの手法を拡張すること。
• 球ベッセル関数を用いた離散化を通じて、帯域制限付きフリードホルム核を体系的に取り扱うフレームワークを提供すること。

提案手法

• 本手法は、帯域制限付き核関数を球ベッセル関数を基底関数として展開することにより、退化核近似を採用する。
• 離散化により、連続的なフリードホルム積分方程式が行列固有値問題に変換され、準楕円体付随関数の数値的解法が可能になる。
• 得られた行列系の固有関数は、古典的な零次のPSWFsを近似することができ、全実数直線上での計算が可能になる。
• 同じ球ベッセル基底展開を適用することにより、任意の帯域制限付き、平方可積分核に対してもこの手法を一般化できる。
• 対象関数をその帯域制限付き核による像から推定する逆問題は、同じ離散化フレームワークを用いて解かれる。
• 関数空間内での球ベッセル関数の直交性および完全性を活用することで、高い精度が保証される。

実験結果

リサーチクエスチョン

• RQ1 零次の準楕円体付随関数を、$[-1,1]$ に限定されず、全実数直線上で効率的かつ正確に計算するにはどうすればよいか？
• RQ2 帯域制限付きフリードホルム核を離散化する最適な方法は何か？ これにより、スペクトル性質が保たれ、関連固有値問題の数値的解法が可能になる。
• RQ3 球ベッセル関数を用いた退化核法を、任意の帯域制限付き、平方可積分核に一般化できるか？
• RQ4 本手法は、PSWF計算および帯域制限付き外挿の既存数値手法と比較して、精度および安定性においてどのように優れているか？
• RQ5 このフレームワークは、信号処理における逆問題、例えば帯域制限付き核を用いたデコンボリューションに、どの程度応用可能か？

主な発見

• 提案手法により、零次の準楕円体付随関数が全実数直線上で正確に計算可能となり、従来の $[-1,1]$ に限定される制限を克服した。
• 球ベッセル関数を基底関数として用いることで、良好な条件数を持つ行列固有値問題が得られ、数値的安定性が向上した。
• 対象関数をその帯域制限付き像から再構成する逆問題を、一般の帯域制限付き核に対して効果的に解くことができ、ロバスト性を示した。
• 離散化スキームは核の帯域制限性を保ち、基底次元を増加させることで真の解に収束することが保証された。
• このフレームワークは、任意の平方可積分で帯域制限付き核に適用可能であり、信号処理および数理物理学におけるフリードホルム方程式の一般用途ツールを提供する。