[论文解读] 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。
• 提供一种模拟框架，通过利用基于网格的跳跃动力学，避免使用复杂的数值格式。

提出的方法

• 扩散过程由二阶椭圆微分算子 A(x) = ∑ᵢⱼ ∂ᵢaᵢⱼ(x)∂ⱼ + bᵢ(x)∂ᵢ 定义，表示非对称广义扩散。
• 在 ℝᵈ 中的均匀且各向同性网格 Gₙ ⊂ ℝᵈ 上构造马尔可夫跳跃过程 Xₙ(·)，以近似扩散过程。
• 跳跃过程的生成元被显式推导，以确保在网格分辨率提高的极限下，与扩散过程的无穷小生成元一致。
• 该构造利用了网格的均匀性和各向同性，简化了转移率计算，并确保分布收敛。
• 对于 d ≥ 3，当扩散张量 {aᵢⱼ(x)} 满足额外结构条件时，该方法适用，从而实现近似。
• 模拟直接基于网格上的跳跃事件进行，避免了随机微分方程的离散化或路径的数值积分。

实验结果

研究问题

• RQ1 在规则网格上定义的马尔可夫跳跃过程能否弱收敛于 ℝᵈ 中的非对称广义扩散过程？
• RQ2 跳跃过程的生成元必须具有何种显式形式，才能确保收敛到目标扩散过程？
• RQ3 网格的各向同性和均匀性如何简化近似构造与模拟过程？
• RQ4 在何种扩散张量条件下，该方法可从 d = 2 扩展至 d ≥ 3？
• RQ5 该方法在代码中可实现的效率如何，尤其针对 d = 2？

主要发现

• 当网格间距趋于零时，马尔可夫跳跃过程 Xₙ(·) 在分布上收敛到目标非对称广义扩散过程 X(·)。
• 跳跃过程的生成元被显式构造，使得可直接模拟而无需求解随机微分方程。
• 对于 d = 2，由于网格结构和转移规则的简单性，该方法在计算机代码中实现非常直接。
• 对于 d ≥ 3，该方法在扩散张量 {aᵢⱼ(x)} 满足额外条件时依然有效，该条件确保了收敛所需的对称性和正则性。
• 该方法提供了一种数学上严谨且计算高效的模拟框架，特别适用于非对称扩散过程的路径模拟。
• 使用各向同性且均匀的网格使得转移率设计具有一致性，简化了实现与分析。