c++ - 指数生成器有时会给出“奇怪”的结果

Question

我正在用 C++ 构建一个模拟，并且我有一个指数生成器来生成进程的突发时间。

通常它会返回如下值：3.14707,1.04998。但经常出现 1/10 次这样的数字：2.64823e-307

这是生成器的代码（我srand ( time(NULL) );在程序开头使用）：

double exponential(float u)
{
   double x,mean;
   mean = 10;
   // generate a U(0,1) random variate
   x = rand();
   u = x / RAND_MAX;
   return (-mean * log(u));
}

这就是我分配值的方式。里面的while部分是我努力摆脱这些价值观，但没有奏效：

for (int i = 0; i < nPages; i++)
{
    index[i] = i;
arrival[i]= poisson(r);
burst[i]=exponential(u);
while (burst[i]<1 || burst[i]>150)
    {
        cout<<"P"<<i<<endl;
        burst[i]=(burst[i-1]+burst[i+1])/2;
    }
}

Answer 1

为什么你使用 C 库而不是 C++ 库？

std::random_device rd;
std::default_random_engine gen(rd());
std::exponential_distribution<double> dist(lambda);

double x = dist(gen);

Answer 2

如果burst的大小是nPages，那么

for (int i = 0; i < nPages; i++)
{
    //...
    burst[i]=(burst[i-1]+burst[i+1])/2;
}

将超出其界限，因此您很可能会胡说八道。您需要考虑边缘需要什么。

至于关于randgo rand 被认为有害的评论值得一看。在您的情况下，记录 0 是不明智的。

Answer 3

使用exponential逐字复制的函数，我无法重现您描述的错误。PRNG 出现的问题要么出现 0，要么RAND_MAX应该只出现一次RAND_MAX，而不是 10% 的时间。我怀疑编译器有问题，或者您共享的不是产生所描述问题的实际代码。