指数生成器有时会给出"weird"结果

Exponential number generator sometimes gives "weird" results

本文关键字：weird 结果 指数      更新时间：2023-10-16

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

如果突发的大小为 nPages，则

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

会超出它的界限，所以你很可能会胡说八道。您需要考虑边缘需要什么。

至于关于rand兰德被认为有害的评论值得一看。在您的情况下，将日志设置为 0 是不明智的。

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

相关文章：