[论文解读] Inter-arrival time distribution for the non-homogeneous Poisson process

G. Yakovlev, John B. Rundle+2 more|arXiv (Cornell University)|Jul 27, 2005

Advanced Frequency and Time Standards被引用 11

一句话总结

本文推导出在有限或无限时间区间内非时齐泊松过程（NHPP）的到达时间间隔分布的精确解析表达式。通过利用均值函数 Λ(t) 并对到达时间密度进行积分，作者得到了一个封闭形式的概率密度函数（pdf），该函数考虑了观测窗口效应，显示出由于过程中存在有限时间记忆而导致的与指数分布的偏离。

ABSTRACT

We derive an analytical expression of the inter-arrival time distribution for a non-homogeneous Poisson process (NHPP). This expression is exact and is applicable to any time interval, finite or infinite. As an illustration, we present simulation results for three different intensity functions.

研究动机与目标

• 推导非时齐泊松过程（NHPP）在有限或无限时间区间内的精确解析表达式，用于到达时间间隔分布。
• 解决在有限观测窗口下，NHPP 到达时间间隔分布缺乏精确解的问题。
• 量化有限观测区间如何引入记忆效应，导致与齐次泊松过程中所见的指数分布产生偏离。
• 通过三种不同强度函数（恒定、衰减和增长）的数值模拟，验证所推导公式的正确性。

提出的方法

• 通过在所有可能的事件数 n 和到达时间 t_k 上求和，利用不相交区间增量的独立性，推导出生存函数 G(x) = P(S > x)。
• 在数值实验中使用反变换法模拟到达时间间隔，确保与 NHPP 强度函数 λ(t) 一致。
• 应用均值函数 Λ(t) = ∫₀ᵗ λ(u)du，将时间映射为累积事件率，从而实现对时间区间内事件概率的精确计算。
• 通过对生存函数 G(x) 求导，推导出到达时间 pdf，得到涉及 λ(t) 和 Λ(t) 的双重积分表达式。
• 通过为首个到达时间间隔增加单独项，考虑有限区间中的边界效应，确保分布的归一化和完整性。
• 通过将解析解与 10⁸ 次实现的模拟结果进行比较，验证解析解的正确性，强度函数分别为：λ(t) = 1/c、c/(1+t) 和 ct。

实验结果

研究问题

• RQ1 在有限时间区间 T 内观测的非时齐泊松过程，其精确到达时间间隔分布是什么？
• RQ2 与无限时间情况相比，有限观测窗口 T 如何影响到达时间间隔分布？
• RQ3 有限区间引入的记忆效应在多大程度上导致与齐次泊松过程中所见的指数分布产生偏离？
• RQ4 能否推导出一个同时考虑时变强度函数和观测窗口边界约束的解析表达式？

主要发现

• 有限时间区间 T 的精确到达时间概率密度函数（pdf）为 pdf(x) = (1/Λ(T)) ∫₀ᵀ⁻ˣ λ(y)λ(y+x)exp{−[Λ(y+x)−Λ(y)]}dy + [λ(x)exp{−Λ(x)}]/Λ(T)，该表达式考虑了边界效应。
• 在齐次情形 λ(t) = 1/c 下，所推导的 pdf 简化为 pdf(x) = (T + c − x)/(cT) exp(−x/c)，显示出由于有限 T 而导致的与指数分布的偏离。
• 模拟结果证实了在 λ(t) = 1/c、c/(1+t) 和 ct 情况下，理论 pdf 与数值结果高度一致，所有测试强度函数下均表现出良好吻合。
• 有限观测窗口引入了过程中的记忆效应，因为在时间 y 之后无法观测到超过 T−y 的到达时间间隔，即使在恒定 λ(t) 下，分布也非指数分布。
• 所推导的表达式是正确归一化的，满足 G(0) = 1 且 G(T) = 0，与生存函数的定义一致。
• 该方法成功地将齐次泊松过程的指数到达时间间隔分布推广至具有任意确定性强度函数的非时齐情形。