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

如果爆发的大小是nPages，那么

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

将走出它的界限，所以你可能最终会胡说八道。你需要考虑边缘需要什么。

关于rand go rand considered harmful的评论值得关注。在你的情况下，记录0是不明智的。

Answer 3

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

指数生成器有时会产生“奇怪”的结果

3 个答案: