Question

我使用以下函数生成高斯随机数：

double r_norm(double mean, double sigma){
    random_device rd;
    mt19937 gen(rd());
    normal_distribution<double> d(mean, sigma);
    return d(gen);
}

但是，当我用cout：

在main（）中调用它时

for (int k = 0; k < 10; k++){ 
    cout << r_norm(2,0.5) <<endl;
}

输出相同的数字10次。理想情况下，我需要能够在任何地方调用此函数，以便每次都接收新生成的数字。

更新：我设法通过将范围之外的随机设备和mersenne twister声明为全局变量来解决这个问题，但有没有更简洁的方法来做到这一点？

Answer 1

问题是您需要回收random_device。

如果你真的想保留完全相同的函数签名，最简单的方法是使用全局变量：

random_device rd;
mt19937 gen(rd());

double r_norm(double mean, double sigma){
    normal_distribution<double> d(mean, sigma);
    return d(gen);
}

话虽如此：stl随机分布是有状态的，所以如果你想要一个实际的有效分布，你还需要回收d。

此时，维护您的界面需要static std::map<pair<double, double>, normal_distribution>，以便您也可以正确回收它们。

c ++创建高斯随机数的函数

1 个答案: