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

我正在用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;
    }
}

最佳答案

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

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

double x = dist(gen);