[論文レビュー] Probabilistic Representation and Approximation for Coupled Systems of Variational Inequalities

Romuald Élie, Idris Kharroubi|arXiv (Cornell University)|Sep 5, 2009

Stochastic processes and financial applications参考文献 5被引用数 4

ひとこと要約

本稿は、ジャンプを伴う1次元制約付き後向きSDEを用いて、結合された準線形変分不等式系の確率的表現を導入し、新たなモンテカルロベースの数値スキームを可能にする。この手法は、レジームスイッチングジャンプ過程を活用し、緩い条件下でも得られるアルゴリズムの収束保証を提供する。

ABSTRACT

Our study is dedicated to the probabilistic representation and numerical approximation of solutions to coupled systems of variational inequalities. The dynamics of each component of the solution is driven by a different linear parabolic operator and suffers a non-linear dependence in all the components of the solution. This dynamics is combined with a global structural constraint between all the components of the solution including the practical example of optimal switching problems. In this paper, we interpret the unique viscosity solution to this type of coupled systems of variational inequalities as the solution to one-dimensional constrained BSDEs with jumps introduced recently in [6]. In the spirit of [3], this new representation allows for the introduction of a natural entirely probabilistic numerical scheme for the resolution of these systems.

研究の動機と目的

• 最適スイッチングおよび確率的制御問題に生じる高次元の結合された準線形変分不等式系を解く数値的課題に対処する。
• レジーム数 $ m $ が大きい場合に古典的PDE手法が計算的に非現実的になるのを克服する。
• 粘性解とジャンプを伴う制約付き後向きSDEを結びつける、確率的フェインマン＝カース型表現を構築する。
• 実用的な数値的解決のため、この表現に基づいた収束性が保証された実装可能なモンテカルロアルゴリズムを設計する。
• 反射BSDEおよびジャンプ付きBSDEに関する先行研究を、レジーム依存のダイナミクスと相関する制約を持つ系へ一般化する。

提案手法

• 結合された変分不等式系の解を、レジーム空間 $ \mathcal{I} $ 上の純ジャンプ過程によって駆動される1次元制約付き後向きSDEの $ Y $-成分として表現する。
• 前向き過程 $ (I_s, X_s) $ を、レジーム $ I_s $ におけるジャンプと $ X_s $ における拡散を伴うSDEとしてモデル化し、ブラウン運動およびポアソン測度によって駆動する。
• ドライバがジャンプ成分 $ U_s(j) = v_j - v_{I_{s-}} $ に依存するBSDEのダイナミクスを導出する。制約は $ h(I_{s-}, j, X_s, Y_{s-}, Y_{s-}+U_s(j), Z_s) \geq 0 $ を満たす。
• 時間分割 $ h $ を用いた離散時間近似スキームを構築し、時間間隔 $ [t_k, t_{k+1}) $ で後向き再帰を実行する。条件付き期待値は最小二乗回帰によって近似する。
• 状態変数 $ (I_{t_k}, X^{h}_{t_k}) $ の基底関数を用いたOLS回帰により、条件付き期待値を近似する。$ N $ 本のシミュレートパスと $ L $ 個の基底関数を用いる。
• 離散時間収束および統計的誤差制御に基づき、$ h \to 0 $、$ N, L \to \infty $ の下で、全アルゴリズムの収束を確立する。

実験結果

リサーチクエスチョン

• RQ1 レジーム依存の生成子と相関する制約を持つ、結合された準線形変分不等式系に対して、確率的表現を確立できるか？
• RQ2 ジャンプを伴う制約付きBSDEの定式化は、PDE系に対して一意な粘性解を与えるか？
• RQ3 この表現に基づいて、数値的に安定かつ収束性が保証されたモンテカルロスキームを構築できるか？
• RQ4 レジーム過程にジャンプを含めることで、拡散のみまたはジャンプなしのモデルと比較して、表現および近似がどのように改善されるか？
• RQ5 時間離散化、パスシミュレーション、および回帰に基づく条件付き期待値近似を組み合わせた全アルゴリズムの収束挙動はいかなるものか？

主な発見

• 結合された変分不等式系の一意な粘性解は、ジャンプを伴う1次元制約付きBSDEの $ Y $-成分として表現される。
• BSDEの定式化は、拡散成分とジャンプ成分を両方含み、ジャンプ部分がレジーム間での価値関数成分の相関を符号化する。
• 時間ステップ $ h \to 0 $ の下で、$ \mathbf{S}^2 $、$ \mathbf{L}^2_W $、$ \mathbf{L}^2_{\tilde{\mu}} $ のノルムにおいて、緩い正則性条件のもとでアルゴリズムが収束する。
• OLS回帰による条件付き期待値の近似は、パス数 $ N $ および基底関数数 $ L $ が無限大に近づく際に収束する。
• 反射BSDEおよびジャンプ付きBSDEに関する先行研究を、レジームスイッチングダイナミクスと非マルコフ的制約を持つ系へ一般化する。
• このアプローチにより、古典的PDE手法が計算的に非現実的になる高次元スイッチング問題の効率的数値的解決が可能になる。