[論文レビュー] Penalty Method for Obliquely Reflected Diffusions

Charles Amponsah, Andrey Sarantsev|arXiv (Cornell University)|Sep 6, 2015

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

ひとこと要約

本稿では、滑らかな領域における確率的微分方程式の斜交反射を近似するためのペナルティ法を提案する。反射項を境界に近づくと内側に押し戻すようなドリフトに基づくペナルティ項に置き換えることで、解を近似する。広範なペナルティ関数のクラスに対して一般化された弱収束結果を確立し、数値シミュレーションにより性能を検証することで、標準的な法線反射に限らない柔軟な代替手法を提供する。

ABSTRACT

We approximate of solutions of stochastic differential equations driven by Brownian motion in a smooth domain and obliquely (not necessarily normally) reflected at the boundary. This approximation is done by replacing reflection term with a penalty function: an additional drift term at the boundary pointing inside the domain. Most literature on this topic is devoted to the half-line (one-dimensional case) or to the normal reflection. We propose a large class of penalty functions and show very general weak convergence results. We perform numerical simulations to assess approximation quality for various penalty functions.

研究の動機と目的

• 滑らかな領域における斜交反射拡散過程に対する一般化されたペナルティベースの近似法を開発すること。これは、標準的な法線反射のケースを超えるものである。
• 特に境界に対して法線でない反射が行われる状況において、多次元設定における斜交反射の収束理論の欠如を是正すること。
• 真の反射拡散過程へのペナルティ過程の弱収束を保証する広範なペナルティ関数のクラスを提案すること。
• 実用的状況における異なるペナルティ関数の近似品質を数値的に評価・比較すること。
• ペナルティ法が正確な反射の代替手段として計算的に実行可能であることを理論的および実証的に裏付けること。

提案手法

• SDEにおける反射項を、境界に近づいた際にパスを領域内に押し戻す力として作用するペナルティドリフト項に置き換える。
• ペナルティ関数は滑らかで内向きを指向するものとし、境界に近づくほどその大きさが増加するように構築する。
• ペナルティ関数および領域の滑らかさに関する一般条件の下で、ペナルティ付き過程の弱収束が確立される。
• マーティングル問題の定式化と弱収束技術を用いて、分布収束の証明を行う。
• 近似の精度と収束速度を評価するために、さまざまなペナルティ関数を用いた数値的シミュレーションを実施する。
• 非法線反射角を許容する多次元拡散過程への応用を実施する。

実験結果

リサーチクエスチョン

• RQ1 斜交反射を持つ反射SDEの解への弱収束を保証するペナルティ法を構築できるか？
• RQ2 ペナルティ関数の選択が、反射拡散過程の近似品質にどのように影響するか？
• RQ3 多次元斜交反射設定において、収束を保証するペナルティ関数の一般的条件は何か？
• RQ4 数値実験において、異なるタイプのペナルティ関数の性能はどのように比較されるか？
• RQ5 ペナルティ法は半直線を越える領域や非法線反射方向へ拡張可能か？

主な発見

• 提案されたペナルティ法は、領域およびペナルティ関数の弱い正則性条件のもとで、ペナルティ付き過程が真の斜交反射拡散過程に弱収束することを示した。
• 多項式的および指数的成長を示す境界近傍におけるペナルティ関数を含む広範な関数クラスが収束をもたらすことが示された。
• 数値的シミュレーションにより、ペナルティ強度が増加するにつれて近似誤差が減少し、収束速度がペナルティ関数の滑らかさに依存することが明らかになった。
• 非法線方向を含むさまざまな反射角においても、本手法は安定して高い性能を示し、その一般性が裏付けられた。
• 収束結果は多次元領域に対しても成り立つことが確認され、一次元または法線反射に限られた先行研究を拡張した。
• 特に解析的解が得にくい複雑な幾何構造において、正確な反射の代替手段として、本手法は実用的かつ柔軟な選択肢を提供する。