[論文レビュー] Stability of densities for perturbed Diffusions and Markov Chains

Valentin Konakov, Anna Kozhina+1 more|arXiv (Cornell University)|Jun 29, 2015

Advanced Mathematical Modeling in Engineering参考文献 31被引用数 3

ひとこと要約

本稿は、非退化拡散過程およびマコフ連鎖の係数に対する摂動の下で、遷移密度の定量的安定性バインディングを確立する。パラメトリックス法とホルダー連続性の仮定を用いて、係数の差の観点から明示的な誤差推定値を導出し、収束速度はドリフトおよび拡散係数の正則性に依存する。特にホルダー連続性または多項式減衰条件の下で成立する。

ABSTRACT

We are interested in studying the sensitivity of diffusion processes or their approximations by Markov Chains with respect to a perturbation of the coefficients.

研究の動機と目的

• 非退化拡散過程およびマコフ連鎖の遷移密度がドリフトおよび拡散係数の摂動に対してどれほど感度を持つかを定量的に評価すること。
• 係数の正則性仮定を最小限に抑えた状況下で、真の密度と摂動付き密度の差に対する明示的な誤差バインディングを提供すること。
• 特にオイラースキームの弱誤差解析の文脈において、SDEの離散時間マコフ連鎖近似への安定性結果を拡張すること。
• モデル不適合、パラメータ推定、数値スキーム収束といった応用分野における係数摂動の影響を分析すること。
• ホルダー連続性および多項式減衰仮定の下で、密度差の収束速度を確立し、摂動の大きさおよび正則性パラメータに明示的な依存関係を示すこと。

提案手法

• 拡散過程およびマコフ連鎖の遷移密度を、ガウスカーネルおよび反復積分を含む級数展開で表現するため、パラメトリックス法を用いる。
• 時間に関する滑らかさの欠如に強く、離散時間マコフ連鎖にも適応可能な、McKean–Singerのパラメトリックス展開の手法を採用する。
• 特に2階および3階の導関数に注目し、パラメトリックス級数の項の摂動を分析することで、密度差のバインディングを導出する。
• ホルダー不等式およびガウス型バインディングを用いて、特にホルダー連続性仮定下でドリフトおよび拡散係数の差を制御する。
• 多項式減衰（例：$ q_m(t,x) riangleq c^d (t)^{-d/2} Q_{M-(d+5+ u)}(x / t^{1/2}) $）を持つ密度の畳み込み解析を導入し、畳み込み積分における非対角および対角領域を扱う。
• 帰納法およびガンマ関数の恒等式を用いて、パラメトリックス展開における反復畳み込みの成長を制御し、$ au^{rac{r heta}{2}} $ の観点から明示的な収束速度を導出する。

実験結果

リサーチクエスチョン

• RQ1 非退化拡散過程の遷移密度は、ドリフトおよび拡散係数の小さな摂動に対してどのように変化するか？
• RQ2 係数が摂動された場合、マコフ連鎖近似の密度が真の拡散密度にどれほど速く収束するか？
• RQ3 係数のホルダー連続性および多項式減衰条件が、摂動下での密度安定性にどのように影響するか？
• RQ4 不規則な係数を持つオイラースキームの弱誤差に対して、パラメトリックス法を用いて明示的な誤差バインディングを導出可能か？
• RQ5 高階導関数は密度差にどのように寄与するか？また、パラメトリックス展開における集中損失はどのように制御されるか？

主な発見

• 係数のホルダー連続性の下で、真の密度と摂動付き密度の差は、$ \frac{\tilde{\rho}_{\text{sup}}}{(t_j - t_k)^{1 - \frac{\nu}{2}}} $ で有界である。ここで $ \tilde{\rho}_{\text{sup}} $ は係数差のノルムの上界を表す。
• マコフ連鎖近似において、密度差は $ |\tilde{p}^h \boxtimes H^{h,(r)} - \tilde{p}_\nu^h \boxtimes H_\nu^{h,(r)}| \triangleq \text{error} \triangleq (r+1)\rho_{\nu,\nu,\text{sup}} \frac{\big\bracevert (1 \bigvee T^{(1-\nu)/2})c_1 \big\bracevert^{r+1} [\text{B}(\nu/2)]^r}{\text{B}(1 + r\nu/2)} \frac{c^d}{(t_j - t_i)^{d/2}} Q_{M-(d+5+\nu)}\big( \frac{y-x}{(t_j - t_i)^{1/2}/c} \big) (t_j - t_i)^{r\nu/2} $ を満たし、$ (t_j - t_i)^{r\nu/2} $ における明示的な収束が示される。
• 多項式減衰密度の畳み込みの解析から、$ I_{t_k}^1(t_i,t_j,x,y) \triangleq \big| \big( q_m(t_k - t_i, \bullet - x) * q_m(t_j - t_k, \bullet - y) \big)(y) \big| \triangleq \text{convolution} \triangleq \bar{c} q_m(t_j - t_i, y - x) $ が、$ m > d $ の下で $ i < k < j $ に対して一様に成立する。
• パラメトリックス展開における集中損失は、定義域を分割し、$ |w|^3 $ および $ |w|^4 $ 項を統合することで制御される。これにより、可積分性を保証するためには $ m > d+5+\nu $ が必要となる。
• 密度差の収束速度は $ O((t_j - t_i)^{r\nu/2}) $ であり、仮定されたホルダー正則性 $ \nu $ の下で最適であり、空間的・時間的に一様に成立する。
• 結果は $ \rho_{\nu,\nu,\text{sup}} $ の意味で係数摂動に対して頑健であり、これは $ (b, \nu) $ と $ (b_\nu, \nu_\nu) $ の差の上界ノルムを測る。バインディングはこの摂動に対して線形にスケーリングされる。