形式陈述
一个窗口内的计数服从 Poisson 分布,并不能决定相邻窗口如何关联。Poisson 过程需要哪些过程级条件?
速率 λ > 0 的齐次 Poisson 过程是右连续的计数型随机过程 公理库 随机过程 Stochastic process · Random process 由同一随机实验产生、按时间或空间指标组织的一族随机变量。 N ( t ) ,满足 N ( 0 ) = 0 ,不交时间区间的增量相互独立 公理库 独立性 Statistical independence 从概率表理解独立性,区分两两、相互和条件独立,并用可计算反例澄清零协方差与条件均值的限度。 ,且
N ( t ) − N ( s ) ∼ Pois ( λ ( t − s ) ) , 0 ≤ s < t . 其中 λ 的单位是“每单位时间的次数”,λ ( t − s ) 才是Poisson 分布 公理库 Poisson 分布 Poisson distribution · 泊松分布 以强度参数描述固定窗口内稀有事件计数的离散分布。 的无量纲参数。该过程几乎必然每次只跳一格,有限区间内只有有限个跳点。
等价构造是取独立同分布的指数间隔 公理库 指数分布 Exponential distribution 具有恒定失效率和连续无记忆性的非负等待时间分布。 X n ∼ Exp ( λ ) ,令 S 0 = 0 、S n = ∑ k = 1 n X k ,再取 N ( t ) = max { n ≥ 0 : S n ≤ t } 。因此它是更新过程 公理库 更新过程 Renewal process 用独立同分布的正间隔构造到达时刻,再由到达时刻反演更新次数。 的特殊情形,但更新过程的一般间隔不产生独立增量。它也是连续时间 Markov 链 公理库 连续时间 Markov 链 Continuous-time Markov chain · CTMC 以指数停留时间和跳转概率构造连续时间链,并将局部速率与非爆炸条件分开。 的纯生特例:从 n 以速率 λ 转到 n + 1 ,没有向下跳变;独立未来增量保证给定当前计数后无需更早历史。
直觉
极短时间 h 内,一次到达的概率是 λ h + o ( h ) ,两次及以上是 o ( h ) 。加上不交区间独立,意味着把时间切得更细不会改变模型:每段都有同样的单位时间机会,没有隐含的班次相位或共同高峰环境。
指数间隔与计数的对应可以从零次事件看出:P ( X 1 > t ) = P ( N ( t ) = 0 ) = e − λ t 。反方向从独立指数间隔出发,任意确定时刻 t 所在间隔的剩余部分仍为新指数时间,并独立于已观察到的历史;这由指数无记忆性与后续间隔独立性得到。因而每次从确定时刻重新观察,未来计数与历史独立且规律相同。对 n ≥ 1 ,累积间隔 S n 的密度为 λ n s n − 1 e − λ s / ( n − 1 ) ! ;与下一段在 t − s 内不结束的概率相乘并积分,得到
Pr ( N ( t ) = n ) = ∫ 0 t λ n s n − 1 ( n − 1 ) ! e − λ t d s = e − λ t ( λ t ) n n ! . n = 0 则直接由首段生存函数给出。这验证了构造具有所定义的全部增量规律。反之,独立 Poisson 增量确定所有有限维分布;右连续计数路径由有理时刻的值确定,所以这两个描述给出同一过程律。
例子与边界
设两条独立输入流速率分别为 λ 1 、λ 2 。在任意长度 h 的窗口,两者计数相加的生成函数为
e λ 1 h ( z − 1 ) e λ 2 h ( z − 1 ) = e ( λ 1 + λ 2 ) h ( z − 1 ) , 不同窗口的总计数也独立,因此合并流仍为速率 λ 1 + λ 2 的 Poisson 过程。这不是只对一个窗口使用分布可加性,而是同时核验了整个过程的增量结构。
反过来,给速率 λ 的每次到达贴上标签,标签相互独立且与原过程独立,以概率 p ∈ ( 0 , 1 ) 分到第一类。长度 h 窗口两类计数的联合生成函数为
E [ s N 1 z N 2 ] = e λ h ( p s + ( 1 − p ) z − 1 ) = e p λ h ( s − 1 ) e ( 1 − p ) λ h ( z − 1 ) . 这个因式分解也可逐项核对:对非负整数 k , l ,
Pr ( N 1 = k , N 2 = l ) = e − λ h ( λ h ) k + l ( k + l ) ! ( k + l k ) p k ( 1 − p ) l = [ e − p λ h ( p λ h ) k k ! ] [ e − ( 1 − p ) λ h ( ( 1 − p ) λ h ) l l ! ] . 所以两类窗口计数独立;不交窗口的原计数及所用标签也独立,这进一步给出两个过程的独立性。其速率分别为 p λ 和 ( 1 − p ) λ 。若 p = 0 或 1 ,其中一路恒为空,可称为退化的零速率过程;它不需要定义参数为零的指数间隔。标签若根据最近的拥堵状态选择,就不再是这里的独立稀疏化。
对 T > 0 、整数 n ≥ 1 ,给定 N ( T ) = n ,到达时刻 ( S 1 , … , S n ) 的联合分布等于 n 个独立 Unif ( 0 , T ) 样本排序后的分布。因此即使已知窗口内恰有两次到达,它们也不是固定落在三等分点;应在 0 < s 1 < s 2 < T 的三角域上均匀,密度为 2 / T 2 。
若一开始随机抽一个速率 Λ ,之后按该速率运行,条件于 Λ 是 Poisson 过程,但不同窗口会共享 Λ 而相关。非齐次 Poisson 过程则保留独立增量,将参数换成 ∫ s t λ ( u ) d u ,但失去平稳增量。批量到达又需要复合计数模型,不能保留“一次只跳一格”的假设。
极端超越的 Poisson 点过程极限 公理库 极端超越的 Poisson 点过程极限 Poisson point-process limit for extremes · Extreme exceedance point process 正则变化独立样本的极端观测保留发生时间与归一化高度后,形成具有幂型高度强度的 Poisson 点过程极限。 从独立重尾样本推导时间—高度平面中的稀有点模型:强度随高度变化,只有远离高度零的区域具有有限点数,因而比固定速率时间计数保留更多极端幅度信息。
推论与应用
Poisson 到达与指数服务使系统当前人数足以预测下一次变化,由此得到M/M/1 队列 公理库 M/M/1 排队模型 M/M/1 queue 从单服务台的生灭速率推出几何稳态、负载条件与平均等待,区分队内和系统内指标。 。它也可用作统一的候选跳变时钟,在均匀化算法 公理库 连续时间链的均匀化 Uniformization · Randomization method for CTMC · 均匀化 用共同 Poisson 时钟和自环表示有界速率连续时间链,并以 Poisson 尾控制数值截断误差。 中表示一般有界速率的连续时间链。真实输入是否满足独立、恒定速率和单个到达,应在使用这些模型以前检查。
参考资料