[論文レビュー] Numerical Algorithms for 1-d Backward Stochastic Differential Equations: Convergence and Simulations

Shigē Péng, Mingyu Xu|arXiv (Cornell University)|Nov 28, 2006

Stochastic processes and financial applications参考文献 22被引用数 3

ひとこと要約

本稿では、1次元ブラウン運動のランダムウォーク近似を用いて、1次元後向きストキャスティック微分方程式（BSDE）および反射型BSDEを解くための、暗黙的および陽的数値スキームの提案と分析を行う。これらのスキームの収束性はリプシッツ条件のもとで確立され、特に1つまたは2つのバリアを持つペナルティ型および制約付きBSDEにおいて、安定性と精度がシミュレーションで示されている。

ABSTRACT

In this paper we study different algorithms for backward stochastic differential equations (BSDE in short) basing on random walk framework for 1-dimensional Brownian motion. Implicit and explicit schemes for both BSDE and reflected BSDE are introduced. Then we prove the convergence of different algorithms and present simulation results for different types of BSDEs.

研究の動機と目的

• 非線形生成子を有する1次元後向きストキャスティック微分方程式（BSDE）を解くための数値アルゴリズムの開発および分析を目的とする。
• これらのアルゴリズムを、1つまたは2つのバリアを有する反射型BSDEに拡張し、ペナルティ型および陽-暗黙スキームを含む。
• リプシッツ条件および可積分性条件のもとで、提案された数値スキームの理論的収束結果を確立すること。
• 数値シミュレーションによるアルゴリズムの実装と検証を通じて、実用的状況における安定性と精度を示すこと。
• 非専門家ユーザーが簡単に使用できる直感的な入出力インターフェースを備えた、ユーザーフレンドリーなソフトウェアパッケージの提供

提案手法

• 時間区間を、1次元ブラウン運動のランダムウォーク近似（ステップサイズ δ = T/n）を用いて離散化する。
• 条件付き期待値を、ベルヌーイ増分上に計算することで、Zを求める暗黙的および陽的スキームを適用して離散BSDE系を解く。
• 解Yまたは制御Zに対する制約を満たすためにペナルティ法を用い、制約違反の正の部分に比例するペナルティ項を導入する。
• 生成子およびペナルティ項から導かれる関数の逆関数を含む固定点方程式を解くことにより、暗黙スキームを実装する。
• 将来の値の条件付き期待値を用いて現在のYおよびZを逐次近似する、陽-暗黙スキームを導入する。
• 解の表面および軌道の動的シミュレーションを通じてスキームの妥当性を検証し、制約およびペナルティの効果を可視化する。

実験結果

リサーチクエスチョン

• RQ1 生成子に関するリプシッツ条件および終端条件の可積分性のもとで、1次元BSDEの暗黙的および陽的数値スキームは収束するか？
• RQ2 ペナルティ型スキームは、1つまたは2つのバリアを有する反射型BSDEをどの程度適切に近似できるか。また、ペナルティパラメータpを大きくし、時間刻みδを小さくする際に収束するか？
• RQ3 制約の強制が、特にp√δ > 1のとき、数値解の安定性および精度に与える影響は何か？
• RQ4 非線形生成子g(t,y,z) = -2|y+z| - 1を有する場合、提案されたスキームはどのように動作するか。解の表面およびプロセス軌道の挙動はいかなるものか？
• RQ5 非専門家ユーザーが簡単に使用できるインターフェースを備えた、ユーザーフレンドリーなソフトウェアパッケージとして、数値アルゴリズムを効果的に実装できるか？

主な発見

• リプシッツ条件および可積分性条件のもとで、n → ∞ および p → ∞ のとき、標準BSDEの暗黙的および陽的スキームはL²およびLβノルム（β < 2）で収束する。
• 1つのバリアを有する反射型BSDEに対しては、暗黙的および陽的反射スキームの両方が、時間刻みδ → 0およびペナルティパラメータp → ∞ のとき、真の解に収束する。
• YまたはZに制約があるBSDEに対するペナルティ型陽-暗黙スキームは、同じ条件下でL²およびLβノルム（β < 2）で収束する。
• シミュレーションの結果、p√δ > 1 であると、数値解が不安定になり、発散する可能性があることが示され、安定性の臨界閾値が存在することが示唆された。
• 解の表面Y^{p,n}(t, B_t^n)およびプロセスA_t^{p,n}（ペナルティ項）が明確に可視化され、YまたはZが許容領域から逸脱した場合に制約が効果的に強制されていることが確認された。
• ソフトウェアパッケージは、便利なインターフェースを備えており、非専門家でも複雑なBSDEおよび反射型BSDEを簡単にシミュレーションできることが確認された。