[論文レビュー] Shannon wavelet approximations of linear differential operators

Erwan Deriaz|ArXiv.org|Jan 16, 2007

Image and Signal Denoising Methods参考文献 16被引用数 5

ひとこと要約

本稿では、線形微分作用素の近似にシャノンウェーブレットを用いたウェーブレットベースのアルゴリズムの収束条件を確立し、特にレーライ射影作用素と陰的ラプラシアンに注目している。スペクトル分解とカンタロヴィッチ不等式を用いて収束を証明し、シャノンウェーブレットでは理論的収束率 ≤ 0.56、ウェーブレットパッケージでは ≤ 0.17 を達成しており、ナビエ＝ストークス方程式を含むPDEの効率的数値解法が可能になる。

ABSTRACT

Recent works emphasized the interest of numerical solution of PDE's with wavelets. In their works, A.Cohen, W.Dahmen and R.DeVore focussed on the non linear approximation aspect of the wavelet approximation of PDE's to prove the relevance of such methods. In order to extend these results, we focuss on the convergence of the iterative algorithm, and we consider different possibilities offered by the wavelet theory: the tensorial wavelets and the derivation/integration of wavelet bases. We also investigate the use of wavelet packets. We apply these extended results to prove in the case of the Shannon wavelets, the convergence of the Leray projector algorithm with divergence-free wavelets.

研究の動機と目的

• シャノンウェーブレットを用いた線形微分作用素のウェーブレット近似の一般的な収束条件を確立すること。
• テンソリアル分解とウェーブレットパッケージを含めたウェーブレットベースの手法を拡張し、スペクトル局在化を向上させること。
• PDEを解く反復アルゴリズムの理論的基盤を提供すること、特にレーライ射影作用素と陰的ラプラシアンに焦点を当てる。
• 周波数領域におけるスペクトル解析とカンタロヴィッチ不等式を用いて収束率を最適化すること。
• 非圧縮性流体動力学、特にナビエ＝ストークス方程式に対して、効率的で発散のないウェーブレット・ソルバーを実現すること。

提案手法

• 周波数空間における乗法的表現を可能にするフーリエ変換とシンボル理論を用いて、微分作用素を表現する。
• シャノンウェーブレットとウェーブレットパッケージを用い、周波数領域を二進区間で分割する。
• レマリエ＝リオー＝セーの双直交ウェーブレット導出技術を応用して、作用素近似を構築する。
• 周波数空間における射影行列 $M_\omega$ と $L_\omega$ を用いて、レーライ射影作用素の反復アルゴリズムを構築する。
• カンタロヴィッチ不等式を用いて、反復行列 $(Id - M_\omega - L_\omega)$ のスペクトル半径を用いて収束を分析する。
• 周波数区間比 $a/b$ に基づく収束境界を導出する。シャノンでは $b=2a$、ウェーブレットパッケージでは $b=1.5a$ である。

実験結果

リサーチクエスチョン

• RQ1 定数係数を有する線形微分作用素に対して、反復的ウェーブレットアルゴリズムが収束する条件は何か？
• RQ2 テンソリアルウェーブレットとウェーブレットパッケージは、標準的なマルチスケール解析に比べて、微分作用素の近似をどのように改善するか？
• RQ3 シャノンウェーブレットを用いたレーライ射影作用素アルゴリズムの理論的収束率は何か？
• RQ4 ウェーブレットパッケージによる周波数局在化を用いて、収束率を最適化できるか？
• RQ5 発散のないウェーブレット変換の文脈において、反復行列のスペクトル的挙動は収束にどのように影響するか？

主な発見

• シャノンウェーブレットを用いたレーライ射影作用素の反復アルゴリズムは、$L^2$ノルムにおいて収束し、理論的収束率が $\rho \leq 9/16 \approx 0.56$ で抑えられている。
• シャノンウェーブレットパッケージを用いるとかなり収束率が向上し、$\rho \leq 25/144 \approx 0.17$ にまで低下する。これは周波数局在化の向上によるものである。
• 収束解析は、周波数区間の端点比 $a/b$ に依存する反復行列のスペクトル半径にカンタロヴィッチ不等式を適用することで成り立つ。
• 反復的射影により、ベクトル場の発散のない成分と勾配成分を効果的に分離でき、ナビエ＝ストークス方程式の正確な解法が可能になる。
• 理論的枠組みは、ウェーブレットパッケージがPDEの解法に有望なツールであることを支持しているが、実装面ではまだ主に予測的である。
• 結果として、数値PDEソルバーにおける効率的で発散のない前処理行列を構築する目的でシャノンウェーブレットが有効であることが検証された。