[論文レビュー] Random walk approximation of BSDEs with H{\\"o}lder continuous terminal condition

Christel Geiß, Céline Labart+1 more|arXiv (Cornell University)|Jun 19, 2018

Stochastic processes and financial applications参考文献 34被引用数 7

ひとこと要約

本稿は、局所 Hölder 連続な端末条件を有するマルコフ型後向きストキャスティック微分方程式（BSDE）を解くためのランダムウォークに基づく数値スキームの $L_2$-収束速度を確立する。関連する PDE 解の2階微分特性を分析することにより、$Y$-成分の誤差が $h^{\alpha/4}$ の速度で収束し、$Z$-成分が $h^{\alpha/4}/\sqrt{T - v}$ の速度で収束することを証明する。ここで $\alpha$ は Hölder 指数、$h$ は時間刻みサイズである。

ABSTRACT

In this paper we consider the random walk approximation of the solution of a Markovian BSDE whose terminal condition is a locally H{\\"o}lder continuous function of the Brownian motion. We state the rate of the L 2-convergence of the approximated solution to the true one. The proof relies in part on growth and smoothness properties of the solution u of the associated PDE. Here we improve existing results by showing some properties of the second derivative of u in space.

研究の動機と目的

• BSDE の古典的ランダムウォーク近似法の収束速度を確立することで、文献における空白を埋めること。これは、基礎的であるが、長年の未解決問題であった。
• 既存のランダムウォークスキームの結果を拡張するため、関連する PDE の解の正則性、特に2階空間微分に焦点を当てた分析を行うこと。
• Hölder 連続な端末条件の下で、BSDE 解の $Y$ および $Z$ 成分の収束速度を定量的に特定すること。
• 解の導関数に対する鋭い推定値を導出し、$L_2$-ノルムにおける明確な誤差境界を可能にする。
• 数値シミュレーションを通じて理論的結果を検証し、導出された収束速度の実用的関連性を示すこと。

提案手法

• 時間刻みサイズ $h$ に基づく離散時間ランダムウォーク近似 $B^n_t = \sqrt{h}\sum_{i=1}^{[t/h]} \varepsilon_i$ を用い、BSDE 内のブラウン運動を置き換える。ここで $\varepsilon_i$ は i.i.d. のラデマッハ変数である。
• 近似解 $(Y^n, Z^n)$ は、時間点 $t_k = kh$ における後向き再帰的定義により定義され、条件付き期待値と生成子 $f$ を用いる。
• Malliavin 微積分の枠組みを用いて、$Z$ および $Z^n$ をMalliavin重みを介して表現し、$Z$-成分の正則性と誤差構造の分析を可能にする。
• 特に、端末条件 $g$ が局所 Hölder 連続で指数 $\alpha$ を持つという仮定の下で、関連する PDE の解 $u$ の1階および2階空間微分に対する重要な推定値を導出する。
• 誤差伝播を扱うために、特異核 $(w - r)^{-1/2}$ を含む、修正されたGronwall型不等式を適用する。これは、誤差ダイナミクスの非マルコフ性に対処するためである。
• 尾部の制御のため、切り詰めと局所化技術を用い、$Z$-成分の可積分性を保証し、最終的な収束速度の導出を可能にする。

実験結果

リサーチクエスチョン

• RQ1 端末条件が局所 Hölder 連続である場合、BSDE のランダムウォーク近似スキームの $L_2$-収束速度は何か？
• RQ2 関連する PDE の解の正則性、特に2階空間微分の性質が、数値スキームの収束速度にどのように影響するか？
• RQ3 離散近似におけるMalliavin重みと条件付き期待値の構造を活用することで、収束速度を向上させることは可能か？
• RQ4 収束速度は、端末条件の Hölder 指数 $\alpha$ にどのように依存するか？
• RQ5 時間 $v$ が終端時刻 $T$ に近づくにつれて、$Z$-成分の誤差はどのように振る舞い、その振る舞いは定量的に記述可能か？

主な発見

• $Y$-成分の $L_2$-誤差は、時間刻みサイズ $h$ と端末条件の Hölder 指数 $\alpha \in (0,1]$ の下で、$h^{\alpha/4}$ の速度で収束する。
• $Z$-成分の $L_2$-誤差は、$h^{\alpha/4}/\sqrt{T - v}$ の速度で収束し、これは終端時刻 $T$ に近づくと誤差が発散する傾向を示しており、Malliavin 重み表現における特異性と整合的である。
• 証明は、関連する PDE の解 $u$ の2階空間微分に対する新しい推定値に依拠しており、Zhang [36, Theorem 3.2] を局所 Hölder 連続な端末条件の状況に拡張したものである。
• 収束速度は、Malliavin 微積分、PDE の正則性理論、および特異核を含む精密なGronwall型不等式の組み合わせにより導出された。
• 数値シミュレーションにより理論的収束速度が確認され、$Y$ および $Z$ 成分の両方において $h^{\alpha/4}$ スケーリングが妥当であることが示された。
• 解析は、端末条件 $g$ の多項式成長および Hölder 連続性を考慮しており、広範な実用的 BSDE のクラスに適用可能であることを保証している。