[论文解读] On the time of first level crossing and inverse Gaussian distribution

V. K. Malinovskiǐ|arXiv (Cornell University)|Aug 29, 2017

Stochastic processes and statistical mechanics参考文献 4被引用 5

一句话总结

本文提出了一种针对复合更新过程 $V_s - cs$ 穿过水平 $u$ 的首爆时分布的新近似方法，利用逆高斯分布。该方法在临界漂移 $c^* = \mathbb{E}[Y]/\mathbb{E}[T]$ 附近提供了具有竞争力的、渐近精确的近似，误差界为 $\underline{O}\left(\frac{\ln(u+cv)}{u+cv}\right)$，当 $u+cv \to \infty$ 时成立，且在矩和有界密度条件下成立。

ABSTRACT

We propose a new approximation for the distribution of the time of the first level $u$ crossing by the random process $\homV{s}-cs$, where $\homV{s}$, $s>0$, is compound renewal process and $c>0$. It is competitive with respect to existing approximations, particularly in the region around the critical point $c=\cS$ which separates processes with positive and negative drifts. This approximation is tightly related to inverse Gaussian distributions.

研究动机与目标

• 为复合更新过程中的首爆时分布 $\Upsilon = \inf\{s > 0 : V_s - cs > u\}$ 开发一种新的、解析可处理的近似方法。
• 改进临界区域 $c = c^*$ 附近的准确性，即漂移从正到负转变的区域，此区域现有近似方法常失效。
• 建立首爆时分布与逆高斯分布之间的紧密联系，从而通过标准正态分布累积分布函数实现闭式表达。
• 在跳跃和更新时间分布的弱正则性条件下，推导出近似误差的统一界 $\underline{O}\left(\frac{\ln(u+cv)}{u+cv}\right)$。

提出的方法

• 该方法基于逆高斯分布的累积分布函数（CDF），构建一个依赖于时间的近似 $\mathcal{M}_t(u,c,v)$，其参数由底层更新过程的均值和方差导出。
• 通过涉及 $\mu = 1/(1 - cM)$ 和 $\lambda = (u + cv)/(c^2 D^2)$ 的变换正态 CDF 表达式，对条件概率 $\mathbb{P}\{v < \Upsilon \leq t \mid T_1 = v\}$ 建模，其中 $M = \mathbb{E}[T]/\mathbb{E}[Y]$。
• 当 $c \leq c^*$ 时，近似使用标准逆高斯 CDF；当 $c \geq c^*$ 时，应用指数倾斜以反映负漂移行为。
• 推导基于局部中心极限定理框架，结合对涉及 $\exp(-y^2/2)$ 和有理函数的积分进行严格的渐近分析，通过误差项的精确界进行验证。
• 利用关于 $\exp(-y^2/2)$ 和 $1/(y+R)$ 的积分渐近行为的引理，控制近似中的误差项，尤其在 $u+cv$ 较大时。
• 最终通过将主积分分解为两部分，并应用引理 4.9–4.11，证明两部分均为 $\underline{O}(U^{-1})$，从而得出所述收敛速率。

实验结果

研究问题

• RQ1 如何在临界漂移 $c^*$ 附近对复合更新过程的首爆时分布实现高精度近似？
• RQ2 能否在一般矩和密度条件下，基于逆高斯分布导出首爆时 $\Upsilon$ 的闭式近似？
• RQ3 当 $u+cv \to \infty$ 时，此类近似的渐近误差率是多少？能否在 $t$ 和 $v$ 上实现统一有界？
• RQ4 与现有近似方法相比，该方法在漂移接近零的临界区域中的准确性如何？
• RQ5 该近似能否推广至一般初始更新时间分布 $T_1$（超越指数分布），同时保持误差界？

主要发现

• 在 $\mathbb{E}[T^3] < \infty$、$\mathbb{E}[Y^3] < \infty$ 和有界密度条件下，所提出的近似 $\mathcal{M}_t(u,c,v)$ 对 $\mathbb{P}\{v < \Upsilon \leq t \mid T_1 = v\}$ 的误差界为 $\underline{O}\left(\frac{\ln(u+cv)}{u+cv}\right)$，当 $u+cv \to \infty$ 时成立。
• 该方法通过在 $c \geq c^*$ 情况下对逆高斯 CDF 进行倾斜，统一处理了 $c \leq c^*$ 和 $c \geq c^*$ 的情形，确保了在不同漂移区域内的连续性和准确性。
• 在 $T_1$ 为指数分布的情况下，通过将近似对指数密度进行积分，可获得完整的生存函数 $\mathbb{P}(\Upsilon \leq t)$，并得到可处理的表达式。
• 该近似在性能上与现有方法具有竞争力，尤其在 $c \approx c^*$ 的临界区域中，传统近似方法常失效。
• 误差分析依赖于对涉及 $\exp(-y^2/2)$ 和有理函数的积分的精确渐近估计，通过引理 4.9–4.11 建立了误差项的 $\underline{O}(R^{-1})$ 衰减。
• 收敛速率 $\underline{O}(\ln(u+cv)/(u+cv))$ 是紧的且在 $t$ 和 $v$ 上一致，使该方法适用于风险与排队论中的实际应用。