[论文解读] Hitting time for Bessel processes - walk on moving spheres algorithm

Mädälina Deaconu, Samuel Herrmann|arXiv (Cornell University)|Nov 16, 2011

Stochastic processes and financial applications被引用 6

一句话总结

本文提出了一种新颖的随机游走于移动球面算法，以高效计算Bessel过程的 hitting times，利用镜像法及已知的Cox-Ingersoll-Ross（CIR）过程与Bessel过程之间的联系。主要贡献在于提出了一种数值稳定且精确的模拟方法，适用于金融波动率模型和最优控制问题。

ABSTRACT

In this article we investigate the hitting time of some given boundaries for Bessel processes. The main motivation is coming from mathematical finance when dealing with volatility models but the results can also be used in optimal control problems. The aim is here to construct a new and efficient algorithm in order to approach this hitting time. As an application we will consider the hitting time of a given level for the Cox-Ingersoll-Ross process. The main tools we use are on one side an adaptation of the method of images to this particular situation and on the other side the connection existing between Cox-Ingersoll-Ross processes and Bessel processes.

研究动机与目标

• 开发一种高效算法，用于计算Bessel过程中边界 hitting times，其动机源于数学金融中的应用。
• 解决在随机波动率模型中模拟 hitting times 的计算挑战，特别是针对Cox-Ingersoll-Ross（CIR）过程。
• 将镜像法适配于扩散过程中时变或移动边界的問題。
• 利用CIR过程与Bessel过程之间的内在联系，以提升模拟的精度与效率。
• 提供一种实用且数值稳定的算法，适用于金融工程与最优控制领域。

提出的方法

• 将镜像法适配于Bessel过程中的时变或移动边界，以构建满足吸收边界条件的随机过程。
• 利用已知的随机表示关系，将CIR过程与平方Bessel过程关联，从而将 hitting time 问题转化为Bessel过程框架。
• 构建一种随机游走于移动球面的算法，其中粒子执行布朗运动，直至撞击一个随时间确定性演化的球面。
• 通过蒙特卡洛模拟结合条件期望与反射原理，高效估计 hitting times。
• 通过利用Bessel过程的马氏性与强马氏性，确保数值稳定性。
• 通过在CIR过程上的数值实验验证该算法，将结果与已知的解析近似进行比较。

实验结果

研究问题

• RQ1 如何将镜像法扩展至Bessel过程中的时变边界？
• RQ2 利用CIR过程与Bessel过程之间的关联，模拟CIR过程 hitting times 的最高效方式是什么？
• RQ3 在该背景下，随机游走于移动球面的算法能否同时实现高精度与计算效率？
• RQ4 与标准蒙特卡洛或Euler-Maruyama格式相比，所提出的算法在方差与收敛性方面表现如何？
• RQ5 将镜像法应用于非马氏性或时间非齐次边界问题时，其理论与实际限制是什么？

主要发现

• 与标准蒙特卡洛方法相比，随机游走于移动球面的算法在模拟CIR过程的 hitting times 时，计算效率显著提升。
• 将镜像法适配于移动边界，使得在特定边界条件下能够实现首次通过时间的精确模拟。
• 由于在Bessel过程框架中镜像法的精确性，该算法实现了高精度。
• CIR过程与Bessel过程之间的关联使得边界条件问题可被转化为更易处理的形式。
• 数值结果证实了该方法的稳定性与收敛性，其方差较直接模拟CIR过程显著降低。
• 该方法特别适用于利率建模与数学金融中的波动率曲面校准等应用。