符合泊松分布的事件模拟到达时间生成

符合泊松分布的事件模拟到达时间生成

如果某个事件的到达符合泊松分布，那么如何模拟这些事件的到达时间呢
首先要明白如果事件到达符合泊松分布的话，那么事件到达时间的间隔符合指数分布。
事件的达到时间可以用Knuth的算法生成
至于原理，可以参考博文，主要借助了CDF（累积分布函数）的思想。生成的代码如下：

//生成泊松随机数，其中lambda是到达率
public int nextPoisson() {      
    double L = Math.exp(-lambda);
    double p = 1.0;
    int k = 0;
    Random r = new Random();

    do {
        k++;
        p *= r.nextDouble();
    } while (p > L);
    return k - 1;
}
//生成事件的到达时间
public void generateArrival(){
    int[] arrival_time;
    arrival_time[0]=r.nextPoisson();
    for(int i=1;i<arrival_time.length;i++){
        arrival_time[i]=arrival_time[i-1]+r.nextPoisson();
    }
}

这里需要注意lambda的取值，假设事件平均每秒达到2个，那么lambda取值为0.5；如果事件每2秒到达一个，那么lambda取值为2。