[论文解读] Extreme Events for Fractional Brownian Motion with Drift: Theory and Numerical Validation

Maxence Arutkin, Benjamin Walter+1 more|arXiv (Cornell University)|Aug 28, 2019

Stochastic processes and financial applications被引用 1

一句话总结

该论文为带漂移的分数布朗运动（fBm）的首达时间、最大值分布和吸收概率构建了微扰理论，基于标准布朗运动展开的ε = H − 1/2展开。该研究推导出了一阶解析修正项，并通过使用自适应二分法的高精度数值模拟进行了验证，结果表现出极佳的一致性，且精度超过基于标度的估计方法。

ABSTRACT

We study the first-passage time, the distribution of the maximum, and the absorption probability of fractional Brownian motion of Hurst parameter $H$ with both a linear and a non-linear drift. The latter appears naturally when applying non-linear variable transformations. Via a perturbative expansion in $\epsilon = H-1/2$, we give the first-order corrections to the classical result for Brownian motion analytically. Using a recently introduced adaptive bisection algorithm, which is much more efficient than the standard Davies-Harte algorithm, we test our predictions for the first-passage time on grids of effective sizes up to $N_{ m eff}=2^{28}\approx 2.7 imes 10^{8}$ points. The agreement between theory and simulations is excellent, and by far exceeds in precision what can be obtained by scaling alone.

研究动机与目标

• 解决带漂移的分数布朗运动（fBm）在首达时间与极值统计方面缺乏解析结果的问题。
• 通过在ε = H − 1/2处对标准布朗运动进行微扰展开，将已知的H = 1/2时的标准布朗运动结果推广至一般H ≠ 1/2的Hurst参数。
• 同时包含线性和非线性漂移项，其动机源于非线性变换（如Cole-Hopf变换）。
• 通过高效的自适应二分法实现高精度数值模拟，验证理论预测。

提出的方法

• 以ε = H − 1/2为小参数，在标准布朗运动附近进行微扰展开，计算极值统计的一阶修正。
• 采用路径积分形式和图解展开系统地计算首达时间分布与吸收概率的修正项。
• 推导出编码漂移参数与Hurst参数H依赖关系的标度函数与辅助函数。
• 采用自适应二分法，实现对高达Neff ≈ 2.7 × 10^8个点的网格上首达时间的高效且精确的蒙特卡罗采样。
• 通过高分辨率模拟对理论预测进行数值验证，其精度超越了仅依赖标度的方法。
• 通过如yt = e^{zt}等变换引入非线性漂移，从而在变换后的过程中生成有效非线性漂移项。

实验结果

研究问题

• RQ1 对于H ≠ 1/2的fBm带漂移情况，其首达时间分布的一阶修正项是什么？相较于标准布朗运动情形，这些修正如何体现？
• RQ2 线性和非线性漂移如何影响fBm的极值统计特性，包括最大值分布与吸收概率？
• RQ3 以ε = H − 1/2为展开参数的微扰展开能否准确预测带漂移fBm的行为？其与数值模拟结果相比表现如何？
• RQ4 Hurst参数H在改变带漂移fBm的极值统计特性中起到何种作用，特别是在H < 1/2与H > 1/2时？
• RQ5 与Davies-Harte等标准方法相比，自适应二分法在首达时间模拟中如何提升精度与效率？

主要发现

• 以ε = H − 1/2为展开参数的微扰展开成功计算了带线性与非线性漂移的fBm的首达时间分布与吸收概率的一阶修正。
• 理论预测的首达时间分布与高分辨率数值模拟结果高度一致，即使在Neff ≈ 2.7 × 10^8个点的高分辨率下也表现优异。
• 通过如yt = e^{zt}等变换引入非线性漂移，可生成有效非线性漂移项，且这些项可被精确地纳入理论模型。
• 自适应二分法使模拟精度显著高于基于标度的估计方法，从而验证了理论框架的可靠性。
• 所推导的标度函数与辅助函数（如式91、94、101、105）完整刻画了首达时间分布与H、漂移参数及标度变量的关系。
• 研究结果确认，带漂移fBm的极值统计显著偏离标准布朗运动的预测，尤其在H < 1/2时更为明显，而微扰理论能准确捕捉这些偏离。