[论文解读] Exact Simulation of One-dimensional Stochastic Differential Equations involving the local time at zero of the unknown process

Pierre Étoré, Miguel Martínez|arXiv (Cornell University)|Feb 13, 2011

Stochastic processes and financial applications参考文献 12被引用 4

一句话总结

本文将一维扩散过程的精确模拟方法扩展至包含在零点处局部时的扩散过程，通过Girsanov变换将目标过程与带 drift 的偏斜布朗运动联系起来。核心贡献在于提出了一种模拟算法，通过利用带 drift 的偏斜布朗运动的可处理桥接分布，实现轨迹的精确抽样，从而克服了 SDE 动力学中局部时项带来的挑战。

ABSTRACT

In this article we extend the exact simulation methods of Beskos et al. to the solutions of one-dimensional stochastic differential equations involving the local time of the unknown process at point zero. In order to perform the method we compute the law of the skew Brownian motion with drift. The method presented in this article covers the case where the solution of the SDE with local time corresponds to a divergence form operator with a discontinuous coefficient at zero. Numerical examples are shown to illustrate the method and the performances are compared with more traditional discretization schemes.

研究动机与目标

• 开发一种针对包含零点局部时项的一维 SDE 的精确模拟方法，此类 SDE 未被经典精确模拟技术覆盖。
• 解决此类 SDE 由于局部时分量的存在，与 Wiener 测度不绝对连续所带来的挑战。
• 将 Beskos 等人的精确模拟框架扩展至在零点发生偏斜的过程，特别是那些由具有不连续系数的散度型算子所生成的过程。
• 实现具有不连续漂移和零点局部时的扩散过程的模拟，此类过程在建模不连续介质时具有实际意义。
• 为未来扩展至多个偏斜点或无界漂移提供基础。

提出的方法

• 使用 Girsanov 变换将目标过程的分布表示为相对于带 drift 的偏斜布朗运动的绝对连续测度，从而在 Radon-Nikodym 导数中消除局部时项。
• 利用带 drift 的偏斜布朗运动的已知转移密度，计算 Girsanov 变换中的似然比。
• 应用 Pitman-Yor (1991) 提出的时间反转与点蚀分解技术，刻画带 drift 的偏斜布朗运动的桥接分布。
• 利用最大值与首次 hitting 时间的联合分布，将桥接路径分解为独立片段，从而实现精确路径重构。
• 通过条件化于最大值与 hitting 时间，利用带 drift 的偏斜布朗运动下首次通过时间的已知分布，模拟桥接路径。
• 通过结合 Girsanov 测度变换与参考过程——带 drift 的偏斜布朗运动桥的精确模拟，构建完整的精确模拟算法。

实验结果

研究问题

• RQ1 能否将精确模拟方法扩展至包含零点局部时项的一维 SDE，而这些 SDE 超出了经典精确模拟方法的适用范围？
• RQ2 在 Girsanov 变换中，如何处理零点处的局部时项，以确保参考测度具有可处理性？
• RQ3 带 drift 的偏斜布朗运动的桥接分布是什么？能否实现其精确模拟以支持路径重构？
• RQ4 是否可以放松精确模拟算法中对漂移函数有界性的假设，如同在经典 SDE 中所做的那样？
• RQ5 偏斜参数 β ≠ 0 对具有不连续系数的扩散过程模拟有何影响？

主要发现

• 带零点局部时的 SDE 解的分布相对于带 drift 的偏斜布朗运动的分布是绝对连续的，从而支持基于 Girsanov 的精确模拟。
• 带 drift 的偏斜布朗运动的转移密度被显式推导，并用于计算 Girsanov 变换中的似然比。
• 带 drift 的偏斜布朗运动的桥接分布基于 Pitman-Yor (1991) 的结果，通过最大值与首次 hitting 时间的联合分布进行刻画。
• 通过时间反转与点蚀分解，桥接路径可被分解为独立片段，从而实现完整轨迹的精确模拟。
• 数值示例验证了该方法的准确性，并显示其性能优于标准的欧拉离散化方案等传统方法。
• 该方法适用于对应于在零点具有不连续系数的散度型算子的 SDE，例如在生态学、地球物理学和生物医学建模中出现的情形。