[論文レビュー] On Wavelet-Galerkin methods for Semilinear Parabolic Equations with Additive Noise

Mihály Kovács, Stig Larsson+1 more|arXiv (Cornell University)|Aug 2, 2012

Image and Signal Denoising Methods参考文献 13被引用数 6

ひとこと要約

本稿では、加法的ノイズを伴う半線形確率的放物型方程式を解くための新しいウェーブレット・ガレルキン法を提案する。問題は、非適応的ウェーブレット・ガレルキン法で解く線形確率的方程式と、適応的ウェーブレット・ガレルキン法で解く非線形確率的発展方程式に分解される。主な貢献は、ノイズおよび初期データの正則性に関する仮定の下で、平均二乗誤差推定値が得られ、収束次数 $ \tau^{\beta/2} $ が最適であることを示したことである。$ 0 \leq \beta < 1 $。

ABSTRACT

We consider the semilinear stochastic heat equation perturbed by additive noise. After time-discretization by Euler's method the equation is split into a linear stochastic equation and a non-linear random evolution equation. The linear stochastic equation is discretized in space by a non-adaptive wavelet-Galerkin method. This equation is solved first and its solution is substituted into the nonlinear random evolution equation, which is solved by an adaptive wavelet method. We provide mean square estimates for the overall error.

研究の動機と目的

• 加法的ノイズによって駆動される半線形確率的放物型方程式の数値解法を開発すること。特に、ノイズが粗いまたはマルコフ的でない場合を想定する。
• 線形確率的部と非線形決定的部を分離することで、このような方程式を効率的に解く課題に対処すること。
• 線形確率的成分には非適応的ウェーブレット・ガレルキン法を、非線形発展には適応的ウェーブレット・ガレルキン法を組み合わせ、異なる正則性特性を活用すること。
• 時間的および空間的離散化誤差を含む、全体の数値スキームに対する厳密な平均二乗誤差推定値を確立すること。
• ノイズおよび初期データの正則性をパrameter $ \beta $ で定量化したヒルバート＝シュミットノルム条件 $ \|A^{(\beta-1)/2}Q^{1/2}\|_{\mathrm{HS}} < \infty $ を用いて、収束率を導出すること。

提案手法

• ローテの方法を適用：まず時間に後退オイラー法を用いて確率的PDEを離散化し、解を $ u^n = v^n + w^n $ に分割する。ここで $ w^n $ は線形確率的問題を、$ v^n $ は非線形問題をそれぞれ解く。
• 線形確率的方程式 $ w^n $ を、次数 $ m $ の多重分解空間 $ S_J $ に基づく非適応的ウェーブレット・ガレルキン法で空間に離散化し、安定性と収束性を保証する。
• 非線形方程式 $ v^n $ を、解の特徴を効率的に捉えるために空間基底を動的に細分化する適応的ウェーブレット・ガレルキン法で離散化する。
• SPDEの弱解表現を用い、解を確率的畳み込み $ w(t) = \int_0^t e^{-(t-s)A} dW(s) $ と $ v $ のための確率的発展方程式に分割することで、逐次的解法を可能にする。
• 半群 $ e^{-tA} $ の近似に、時間に依存するリゾルベント作用素 $ r^n(\tau A) $ を用い、ウェーブレット・ガレルキン射影の安定性および近似性質を活用する。
• イトの等長性とスペクトル関数計算を用いて、確率的畳み込みの誤差を評価し、ヒルバート＝シュミットノルム条件 $ \|A^{(\beta-1)/2}Q^{1/2}\|_{\mathrm{HS}} < \infty $ を活用する。

実験結果

リサーチクエスチョン

• RQ1 ウェーブレット・ガレルキン法は、加法的ノイズを伴う半線形確率的放物型方程式に対して最適収束率を達成できるか？
• RQ2 線形確率的部と非線形決定的部の分離が、全体の誤差および数値的安定性にどのように影響するか？
• RQ3 ノイズおよび初期データの正則性を考慮した場合、時間ステップ $ \tau $ におけるウェーブレット・ガレルキンスキームの最適収束率は何か？
• RQ4 このクラスのSPDEに対して、非適応的および適応的ウェーブレット・ガレルキン手法は、誤差制御および計算効率の観点でどのように比較できるか？
• RQ5 共分散作用素 $ Q $ および初期データにどのような条件下で、数値スキームが平均二乗収束を維持するか？

主な発見

• ノイズおよび初期データの正則性に関する仮定 $ \|A^{(\beta-1)/2}Q^{1/2}\|_{\mathrm{HS}} < \infty $ の下で、全体の平均二乗誤差は $ \mathbb{E}\|e^n\|^2 \leq C\tau^\beta \mathrm{e}^{Ct_n} $ を満たす。$ 0 \leq \beta < 1 $。
• 平均二乗ノルムにおける収束率は $ \tau^{\beta/2} $ であり、与えられたノイズおよび初期データの正則性仮定の下で最適である。
• 条件 $ \tau L_f \leq 1/2 $ の下で、非線形部の誤差 $ \bar{v}^n - v^n $ は、確率的畳み込みのウェーブレット近似誤差 $ w_J^n - w^n $ の定数倍で抑えられる。
• 時間離散化された確率的畳み込みの誤差は $ \mathbb{E}\|e_{21}\|^2 \leq C\tau^\beta \|A^{(\beta-1)/2}Q^{1/2}\|_{\mathrm{HS}}^2 $ で抑えられ、ノイズの正則性が収束に与える影響が明確になる。
• 半群の時間離散化誤差は、有理近似 $ r^n(\tau A) $ の安定性により制御され、離散グロワールド補題を用いて導かれた境界が得られる。
• 解析により、非線形部における適応的ウェーブレット法が有効であることが確認された。誤差伝播は、確率的成分の誤差に線形に依存しており、安定性および収束性が保証される。