[論文レビュー] Weak backward error analysis for Langevin process

Marie Kopec|arXiv (Cornell University)|Oct 9, 2013

Markov Chains and Monte Carlo Methods参考文献 14被引用数 4

ひとこと要約

本稿では、確率的ランジュヴィン方程式に適用された陰的数値スキームに対する弱い後退誤差解析を提示する。数値解の生成子がステップサイズの高次項まで修正されたコルモゴロフ作用素を近似することを示し、漸近展開を用いることで修正された不変測度と数値不変測度が非常に近いことを示している。さらに、スキームの力学的性質が無視できる誤差の範囲で指数的混合性を示しており、長時間にわたる安定性と精度を保証している。

ABSTRACT

We consider numerical approximations of stochastic Langevin equations by implicit methods. We show a weak backward error analysis result in the sense that the generator associated with the numerical solution coincides with the solution of a modified Kolmogorov equation up to high order terms with respect to the stepsize. This implies that every measure of the numerical scheme is close to a modified invariant measure obtained by asymptotic expansion. Moreover, we prove that, up to negligible terms, the dynamic associated with the implicit scheme considered is exponentially mixing.

研究の動機と目的

• 確率的ランジュヴィン方程式に対する陰的数値スキームの長時間挙動を分析すること。
• このようなスキームのための弱い後退誤差解析フレームワークを確立すること。
• 数値解の不変測度が漸近展開によって得られる修正された不変測度に非常に近いかどうかを示すこと。
• 陰的スキームの力学的性質が無視できる項の範囲で指数的混合性を示すことを証明すること。
• エルゴディックな拡散過程をシミュレートする際の陰的スキームの精度と安定性に対する理論的根拠を提供すること。

提案手法

• 数値スキームの生成子を用い、ステップサイズの高次補正項を含む修正されたコルモゴロフ方程式と比較する。
• 漸近展開技術を用いて、摂動された生成子に対応する修正された不変測度を導出する。
• 弱収束解析を用いて、数値解と修正されたSDEの解との関係を結びつける。
• 生成子の性質と摂動論を活用して、コルモゴロフ方程式における誤差項を制御する。
• 適切な仮定の下で、修正された生成子の固有値的性質を解析することで指数的混合性を確立する。
• 特に、ドリフト項と拡散係数が与えられたSDEに対して、陰的時間ステッピングスキーム（例：陰的エイラー法）にこのフレームワークを適用する。

実験結果

リサーチクエスチョン

• RQ1 ランジュヴィンSDEのための陰的数値スキームの生成子は、どのように修正されたコルモゴロフ作用素に関連しているか？
• RQ2 漸近展開によって得られる修正された不変測度と、数値不変測度のずれはどの程度小さいか？
• RQ3 無視できる誤差の範囲で、陰的スキームの力学的性質が指数的混合性を示せるか？
• RQ4 ランジュヴィン過程の文脈において、陰的スキームの弱い後退誤差解析の次数は何か？
• RQ5 ステップサイズにおける高次項は、数値解の長時間統計的性質にどのように影響を与えるか？

主な発見

• 数値スキームの生成子は、ステップサイズの高次項まで、修正されたコルモゴロフ方程式の生成子と一致しており、高次の弱い精度を保証する。
• 数値不変測度は、生成子の漸近展開によって得られる修正された不変測度に任意に近い。
• 陰的スキームの力学的性質は指数的混合性を示しており、平衡への速やかな収束と長時間シミュレーションにおける安定性を意味する。
• この解析は一般の陰的スキームに適用可能であり、滑らかな係数を持つ広範なクラスの確率的ランジュヴィン方程式に適用可能である。
• 結果として、統計誤差を制御し、長時間精度を保証する観点から、エルゴディックな拡散過程のシミュレーションに陰的スキームを用いる根拠が得られる。
• このフレームワークは、SDEの不変測度からのサンプリングにおける陰的スキームの信頼性の理論的基盤を提供する。