[論文レビュー] Space-time approximation of stochastic $p$-Laplace systems

Dominic Breit, Martina Hofmanová+1 more|arXiv (Cornell University)|Apr 5, 2019

Stochastic processes and financial applications参考文献 36被引用数 6

ひとこと要約

本稿では、時間正則性が低い問題に対処するための確率的 $p$-ラプラシアン系に対する時空間有限要素法を提示する。ランダムな時間グリッドを導入することで、空間で線形収束、時間で $α \in (0,\frac{1}{2})$ の収束率を達成する。自然な正則性仮定のもとで、数値実験によりバイアス補正済みの収束率が確認された。

ABSTRACT

We consider systems of stochastic evolutionary equations of the $p$-Laplace type. We establish convergence rates for a finite-element based space-time approximation, where the error is measured in a suitable quasi-norm. Under natural regularity assumptions on the solution, our main result provides linear convergence in space and convergence of order $α$ in time for all $α\in(0,\frac{1}{2})$. The key ingredient of our analysis is a random time-grid, which allows us to compensate for the lack of time regularity. Our theoretical results are confirmed by numerical experiments.

研究の動機と目的

• 時空間有限要素近似を、準ノルムにおける誤差を制御できるように開発すること。
• 確率的 $p$-ラプラシアン系における時間正則性の低さという課題に対処するため、ランダムな時間グリッド戦略を導入すること。
• 解の最小限の正則性仮定のもとで、時空間近似の収束率を確立すること。
• バイアス補正済みの収束率推定を用いた数値実験により、理論的結果を検証すること。

提案手法

• 確率的 $p$-ラプラシアン系の変分定式化を用いて、時空間有限要素法を定式化する。
• 誤差は非線形変換 $\mathbf{F}(\boldsymbol{\xi}) = (\kappa + |\boldsymbol{\xi}|)^{\frac{p-2}{2}}\boldsymbol{\xi}$ を含む準ノルムで測定され、これにより誤差構造が線形化される。
• 解の時間正則性の欠如に対応するため、ランダムな時間グリッドを導入し、より良い時間収束解析を可能にする。
• 収束解析は、$\mathbf{F}(\nabla \mathbf{u})$ および $\partial_t \mathbf{u}$ の正則性仮定に依存しており、特に $\mathbf{F}(\nabla \mathbf{u}) \in L^2(0,T;W^{1,2})$ および $\partial_t \mathbf{u} \in L^\infty(0,T;L^2)$ が含まれる。
• 理論的収束率の妥当性を検証するため、数値実験を実施し、真の収束次数を推定する際にバイアス補正回帰を用いる。
• 漸近解析から導かれるバイアス要因 $\beta_{\tau,\tilde{\tau}}(a)$ の関数形を用いて、回帰ベースの収束率推定に対するバイアス補正技術を適用する。

実験結果

リサーチクエスチョン

• RQ1 最小限の正則性仮定のもとで、確率的 $p$-ラプラシアン系の時空間有限要素近似が達成できる収束率は何か？
• RQ2 解の時間正則性の欠如は、収束解析においてどのように補償可能か？
• RQ3 ランダムな時間グリッドの導入が、時空間有限要素スキームの収束挙動に与える影響は何か？
• RQ4 確率的ノイズと離散サンプリングの存在下で、数値的収束率を信頼性高く推定し、バイアスを補正する方法は何か？
• RQ5 数値実験は、解析によって予測された理論的収束率をどの程度確認できるか？

主な発見

• 標準的な正則性仮定のもとで、空間で線形収束、時間で $\alpha \in (0,\frac{1}{2})$ の収束率を達成する。
• 理論的収束率は数値実験により確認され、バイアス補正済み推定値 $\alpha \approx 0.44 \pm 0.07$ が理論的上限と整合的である。
• ランダムな時間グリッド戦略は、時間正則性の低さの影響を効果的に軽減し、滑らかさの欠如にもかかわらず収束率の導出を可能にする。
• 数値実験では、対数-対数回帰における推定収束率が $\tilde{a} \approx 1.3 \pm 0.1$ に達するが、これは $g(\tau) - g(\tilde{\tau})$ を $g(\tau) - g(0)$ の代わりに使用しているため、上向きバイアスが生じている。
• バイアス補正手順により、物理的に意味のある収束次数 $a \approx 0.88 \pm 0.14$（対応する $\alpha \approx 0.44 \pm 0.07$）が回復され、理論的結果と一致する。
• 解析により、非線形確率的 $p$-ラプラシアン系における誤差の定量的評価に、$\mathbf{F}(\nabla \mathbf{u})$ を用いた準ノルムが適切であることが確認された。