c++中的poisson到达码给出了相同的结果

poisson arrival code in c++ gives same result

本文关键字:结果 poisson 中的 c++      更新时间:2023-10-16

你好,我正在尝试模拟一些作业调度算法,并尝试为请求创建一个泊松到达函数。我在维基百科上找到了泊松到达的算法,我实现并运行了它,但我一直得到相同的结果(例如l=15->返回14,l=1/15返回0)

#include <cmath>
#include <iostream>
#include <cstdlib>
using namespace std;
main (){
for (int i=0;i<1000;i++)
{
float l=25;
float L=exp(-l);
float k=0;
float p=1;
//begin while
do{ 
k=k+1;
// Generate a random number between 0 and 1
// return a uniform number in [0,1].
double u = rand() / (double)RAND_MAX;   
p=p*u;
}   while (p>L);    
return k-1;
}}

这是我用来创建这个的算法

algorithm poisson random number (Knuth):
init:
     Let L ← e−λ, k ← 0 and p ← 1.
do:
     k ← k + 1.
     Generate uniform random number u in [0,1] and let p ← p × u.
while p > L.
return k − 1.

提前感谢

您的问题有点不清楚,但我认为您的问题是每次运行程序时都会期望不同的随机数。请注意,每次运行程序时,rand()都会生成相同的随机数序列,除非您为其提供某种种子值。

这样做的标准方法是使用当前时间为随机数生成器播种,如下所示。

#include <iostream>
#include <ctime>
#include <cstdlib>
int main() {
    srand(time(NULL));
    std::cout << rand();
    return 0;
}

这将每次产生不同的随机序列。如果删除srand(时间(NULL))行,那么每次运行程序时都会得到相同的结果。

您只需要在程序开始时调用srand(时间(NULL))一次。

编辑:我也要编辑我的帖子,提到C++11中有一个专用的"随机"标头,我还没有使用过,所以我不能对此发表太多评论,但你可以在这里阅读

http://www.cplusplus.com/reference/random/