[論文レビュー] Pathwise approximation of SDEs by coupling piecewise abelian rough paths

Guy Flint, Terry Lyons|arXiv (Cornell University)|May 6, 2015

Stochastic processes and financial applications参考文献 74被引用数 10

ひとこと要約

本稿では、多変数ブラウン運動によって駆動される確率微分方程式（SDE）のための、新たなパスワイズ近似スキームを提案する。この手法は、区分的アーベル型の粗い経路とガウス近似を用いて、リーヴィ面積の増分を近似する。$ \gamma $-リプシッツなベクトル場（$ \gamma > 2 $）に対して、Wasserstein収束速度 $ O(h^{1-2/\gamma - \varepsilon}) $ を達成し、標準的なEulerスキームの $ O(\sqrt{h}) $ よりも優れている。リーヴィ面積の明示的シミュレーションを必要としない。

ABSTRACT

We present a new pathwise approximation scheme for stochastic differential equations driven by multidimensional Brownian motion which does not require the simulation of Lévy area and has a Wasserstein convergence rate better than the Euler scheme's strong error rate of $O(\sqrt{h})$, where $h$ is the step-size. By using rough path theory we avoid imposing any non-degenerate Hörmander or ellipticity assumptions on the vector fields of the SDE, in contrast to the similar papers of Alfonsi, Davie, and Malliavin et al. The scheme is based on the log-ODE method with the Lévy area increments replaced by Gaussian approximations with the same covariance structure. The Wasserstein coupling is achieved by making small changes to the argument of Davie, the latter being an extension of the Komlós-Major-Tusnády Theorem. We prove that the convergence of the scheme in the Wasserstein metric is of the order $O(h^{1-2/γ-\varepsilon})$ when the vector fields are $γ$-Lipschitz in the sense of Stein.

研究の動機と目的

• $ d \geq 2 $ の場合に計算的に困難なリーヴィ面積増分のシミュレーションを回避するパスワイズSDE近似スキームの開発。
• 特に高次元または退化した設定において、標準的な $ O(\sqrt{h}) $ のEulerスキームより優れたWasserstein収束速度を達成すること。
• 非退化ホルンダーや楕円型条件の必要性を排除し、より広範なSDEのクラスへの適用を可能にすること。
• 粗い経路理論とログ-ODE法を活用し、真のSDE解と近似解との間のカップリングを構築し、Wasserstein誤差を制御すること。
• Komlós-Major-Tusnády定理の枠組みを、SDEの文脈におけるWassersteinカップリングを通じて、多変数ブラウン運動によって駆動される確率過程へと拡張すること。

提案手法

• リーヴィ面積増分の近似に、独立なガウス確率変数 $ W^{(j)} \sim N(0, hI_d) $、$ z^{(j)} \sim N(0, 12^{-1}hI_d) $、および $ \lambda^{(j)} \sim N(0, 12^{-1}h^2I_{d(d-1)/2}) $ を用いる。
• 近似値 $ x_j^h $ は指数写像を用いて反復的に更新される：$ x_{j+1}^h = \exp\left(hV_0 + \sum_k W_k^{(j)} V_k + \sum_{k<l} B_{kl}^{(j)} [V_k, V_l] \right)(x_j^h) $。ここで $ B_{kl}^{(j)} $ は反復積分のガウス近似である。
• 非予測的でパスワイズな近似を扱うために、SDE解を粗い経路フレームワークに埋め込むログ-ODE法を採用する。
• 真のブラウン運動と近似ノイズを同一の確率空間上に一致させることでWassersteinカップリングを構築し、パスワイズ誤差の共同制御を実現する。
• 収束解析はKomlós-Major-Tusnády定理とその粗い経路設定への拡張に依存し、Baker–Campbell–Hausdorff公式のリー代数展開を用いる。
• 主な技術的革新は、リー代数レベルの摂動を用いてBCH展開の差をバインドすることであり、グローバルな群レベル誤差を制御するための補題11.3を活用する。

実験結果

リサーチクエスチョン

• RQ1 リーヴィ面積のシミュレーションを回避しつつ、$ O(\sqrt{h}) $ よりも優れた収束速度を達成できるパスワイズSDE近似を構築できるか？
• RQ2 退化または非楕円型のベクトル場を有するSDEに対して、Wasserstein距離における収束速度をEulerスキームの $ O(\sqrt{h}) $ よりも向上させられるか？
• RQ3 明示的なリーヴィ面積計算を一切行わず、ガウス確率変数のみを用いてWasserstein距離で高次の収束を得ることは可能か？
• RQ4 SDEの文脈において、多変数ブラウン運動によって駆動される確率過程にKomlós-Major-Tusnády定理をどのように適応できるか？
• RQ5 粗い経路に基づく近似スキームを用いた場合、$ \gamma $-リプシッツなベクトル場を有するSDEのWasserstein距離における最適収束速度は何か？

主な発見

• 提案スキームは、任意の $ \varepsilon > 0 $ に対して、$ \gamma > 2 $ であるベクトル場のホルダー・リプシッツ正則性を満たす場合、Wasserstein収束速度 $ O(h^{1-2/\gamma - \varepsilon}) $ を達成する。
• リーヴィ面積のシミュレーションを完全に回避し、ブラウン運動の反復積分と一致する共分散構造を持つガウス近似に依存する。
• 標準的なEulerスキームの $ O(\sqrt{h}) $ よりも収束速度が向上し、特に $ \gamma > 4 $ の場合、収束速度は $ O(h^{1-\varepsilon}) $ にまで向上する。
• 最小限の仮定の下で有効である：$ V_0 \in \text{Lip}^1(\mathbb{R}^q) $、$ V \in \text{Lip}^\gamma(\mathbb{R}^q) $ で $ \gamma > 2 $、非退化性やホルンダ型条件は必要ない。
• Wassersteinカップリングは明示的に同一の確率空間上に構築されており、真の解と近似解の同時分布的制御が保証される。
• 理論的基盤は、リー代数におけるBaker–Campbell–Hausdorff公式の精密解析にあり、リー代数の摂動による群元の差の新しいバインド（補題11.3）を用いる。