[论文解读] On the first hitting times of one dimensional elliptic diffusions

Noufel Frikha, Arturo Kohatsu‐Higa+1 more|arXiv (Cornell University)|Sep 29, 2016

Markov Chains and Monte Carlo Methods参考文献 8被引用 5

一句话总结

本文使用参数化方法，研究了一维均匀椭圆扩散过程中首次 hitting 时间与吸收扩散过程的转移密度的存在性、正则性及高斯上界。在漂移和扩散系数具有较弱光滑性假设下，推导了密度的显式级数展开，并提供了可用于无偏蒙特卡洛模拟的概率表示。

ABSTRACT

In this article, we obtain properties of the law associated to the first hitting time of a threshold by a one-dimensional uniformly elliptic diffusion process and to the associated process stopped at the threshold. Our methodology relies on the parametrix method that we apply to the associated Markov semigroup. It allows to obtain explicit expressions for the corresponding transition densities and to study its regularity properties up to the boundary under mild assumptions on the coefficients. As a by product, we also provide Gaussian upper estimates for these laws and derive a probabilistic representation that may be useful for the construction of an unbiased Monte Carlo path simulation method, among other applications.

研究动机与目标

• 分析一维均匀椭圆扩散过程中首次 hitting 某一阈值的分布。
• 研究关联吸收扩散过程在边界处停止的转移密度。
• 建立首次 hitting 时间密度及吸收过程密度的正则性性质与精确高斯上界。
• 通过参数化展开发展一种概率表示，以期用于无偏蒙特卡洛路径模拟。
• 利用基于半群的参数化方法，将现有首次通过时间密度结果推广至光滑性较低的系数情形。

提出的方法

• 将参数化方法应用于与吸收扩散过程相关的马尔可夫半群。
• 构造转移核 $ P_t h(u,x) $ 的无限级数展开，表示 $ (u + \tau^{u,x}_{u+t}, X^{u,x}_{u+\tau^{u,x}_{u+t}}) $ 的分布。
• 使用扰动技术处理由于在边界 $ L $ 处吸收导致的联合分布的奇异行为。
• 推导密度分量的显式表达式：$ p^{K,N}(x,t) $ 表示 hitting 时间，$ p^{D,N}_T(x,z) $ 表示吸收状态。
• 在紧子集上建立 $ p^{K,N}(x,t) $ 和 $ p^{D,N}_T(x,z) $ 级数的统一且绝对收敛性。
• 应用控制收敛定理并取极限 $ N \to \infty $，恢复极限密度 $ p_T(u,x,dt,dz) $，其高斯上界与 $ N $ 无关。

实验结果

研究问题

• RQ1 对于具有非光滑系数的一维椭圆扩散过程，首次 hitting 时间密度是否存在且保持正则？
• RQ2 参数化方法能否被调整以处理首次 hitting 时间与吸收过程的奇异联合分布？
• RQ3 首次 hitting 时间密度与吸收扩散过程密度的精确高斯上界是什么？
• RQ4 能否推导出支持无偏蒙特卡洛模拟的转移密度的概率表示？
• RQ5 密度的正则性与上界如何依赖于漂移和扩散系数的光滑性？

主要发现

• 在系数具有较弱光滑性假设下，过程 $ (u + \tau^{u,x}_{u+t}, X^{u,x}_{u+\tau^{u,x}_{u+t}}) $ 的转移密度存在且光滑。
• 密度具有统一的高斯上界：$ p^{K,N}(x,t) \leq C t^{-1/2} g(ct, L - x) $ 与 $ p^{D,N}_T(x,z) \leq C g(cT, z - x) $，其中常数 $ C, c > 1 $ 与 $ N $ 无关。
• 在 $ (0,T] \times (-\infty,L] \times \mathbb{R} $ 的紧子集上，$ p^{K,N}(x,t) $ 与 $ p^{D,N}_T(x,z) $ 的无限级数展开绝对且一致收敛。
• 极限密度 $ p_T(u,x,dt,dz) $ 作为 $ N \to \infty $ 的极限获得，且保持高斯上界。
• 半群 $ P_T h(u,x) $ 的概率表示有效，且其级数结构支持无偏模拟。
• 该方法可推广至非光滑系数情形，为退化或光滑性较低情形下的密度估计提供了一种稳健的替代方案，优于马尔可夫场微积分。