[論文レビュー] Markov Jump Processes Approximating a Nonsymmetric Generalized Diffusion

N. Limić|arXiv (Cornell University)|Apr 5, 2008

Advanced Mathematical Modeling in Engineering参考文献 2被引用数 3

ひとこと要約

本稿では、ℝᵈ における非対称一般化拡散過程を、均一で等方的なグリッド上のマコフ跳躍過程を用いて近似する手法を提案する。これらの跳躍過程の生成作用素は明示的に構成されており、目的の拡散過程への弱収束を保証する。この手法により、特に d = 2 の場合に高効率なサンプルパスのシミュレーションが可能となり、d ≥ 3 の場合にも追加の条件下で適用可能である。

ABSTRACT

Consider a nonsymetric generalized diffusion X(·) in R d generated by the differential operator A(x) = ∑ ∑ ∂iaij(x)∂j + bi(x)∂i. In this paper the diffusion process is ij i approximated by Markov jump processes Xn(·) in homogeneous and isotropic grids Gn ⊂ R d which converge in distribution to diffusion. The generators of Xn(·) are constructed explicitly. Due to the homogeneity and isotropy of grids the proposed method for d ≥ 3 can be applied to processes for which the diffusion tensor {aij(x)} dd 11 fulfills an additional condition. The proposed construction offers a simple method for simulation of sample paths of nonsymetric generalized diffusion. Simulations are carried out in terms of jump processes Xn(·). For d = 2 the construction can be easily implemented into a computer code.

研究の動機と目的

• ℝᵈ における非対称一般化拡散過程のサンプルパスをシミュレートする計算的に実行可能な手法の開発。
• 目的の拡散過程に弱収束する構造的グリッド上のマコフ跳躍過程の構築。
• 特に d = 2 の場合に明示的かつ実装可能であるように構成を保証し、d ≥ 3 の場合の拡張は、拡散テンソルに関する追加条件のもとで可能とする。
• SDE の離散化やパスの数値積分を避けることで、複雑な数値スキームを回避するシミュレーションフレームワークの提供。

提案手法

• 拡散過程は、2階の楕円型微分作用素 A(x) = ∑ᵢⱼ ∂ᵢaᵢⱼ(x)∂ⱼ + bᵢ(x)∂ᵢ で定義され、非対称一般化拡散を表す。
• 均一かつ等方的なグリッド Gₙ ⊂ ℝᵈ 上に、マコフ跳躍過程 Xₙ(·) を構築して、拡散過程を近似する。
• 跳躍過程の生成作用素は、グリッド解像度が向上する極限において、元の拡散過程の無限小生成作用素と一致するように明示的に導出される。
• グリッドの均一性と等方性を活用することで、遷移率の計算が簡素化され、分布収束が保証される。
• d ≥ 3 の場合、拡散テンソル {aᵢⱼ(x)} が追加の構造的条件を満たす場合にのみ、この手法が適用可能である。
• シミュレーションは、SDE の離散化やパスの数値積分を避けて、グリッド上の跳躍イベントそのものに基づいて直接実行される。

実験結果

リサーチクエスチョン

• RQ1 正則なグリッド上でのマコフ跳躍過程は、ℝᵈ における非対称一般化拡散過程を弱収束的に近似可能か？
• RQ2 跳躍過程の生成作用素が、目的の拡散過程に収束するためには、どのような明示的形をとる必要があるか？
• RQ3 グリッドの等方性と均一性は、近似の構築とシミュレーションをどのように簡素化するか？
• RQ4 d = 2 から d ≥ 3 への拡張は、拡散テンソルにどのような条件下で可能か？
• RQ5 この手法は、特に d = 2 の場合に、コード実装においてどの程度効率的に実現可能か？

主な発見

• グリッド間隔がゼロに近づくにつれて、マコフ跳躍過程 Xₙ(·) は、目的の非対称一般化拡散過程 X(·) に分布収束する。
• 跳躍過程の生成作用素は明示的に構成されており、確率微分方程式を解くことなく直接シミュレーションが可能である。
• d = 2 の場合、グリッド構造と遷移ルールの単純さのおかげで、コンピュータコードへの実装が容易である。
• d ≥ 3 の場合、拡散テンソル {aᵢⱼ(x)} が追加の条件を満たす限り、この手法は依然として有効であり、収束に必要な対称性と正則性を保証する。
• 本手法は、数学的に厳密でありながら、計算的にも効率的なシミュレーションフレームワークを提供し、非対称拡散過程のパスシミュレーションに特に適している。
• 等方的かつ均一なグリッドの使用により、均一な遷移率設計が可能となり、実装と解析が簡素化される。