Question

我正在尝试使用 C ++ STD TechnicalReport1 扩展程序按照正态分布生成数字，但此代码（改编自this article）：

mt19937 eng;
eng.seed(SEED);

normal_distribution<double> dist;
// XXX if I use the one below it exits the for loop
// uniform_int<int> dist(1, 52);

for (unsigned int i = 0; i < 1000; ++i) {
  cout << "Generating " << i << "-th value" << endl;
  cout << dist(eng) << endl;
}

只打印1“Generating ...”日志消息，然后永远不会退出for循环！如果我使用我注释掉的发行版，它会终止，所以我想知道我做错了什么。有什么想法吗？

非常感谢！

Answer 1

我对最初发布的代码存在同样的问题，并调查了

的GNU实现

首先观察一些事情：使用g ++ - 4.4并使用代码挂起，使用g ++ - 4.5并使用-std = c ++ 0x（即不是TR1，但真实的东西）代码正常工作

恕我直言，关于随机数生成和随机数消耗之间的适配器，TR1和c ++ 0x之间发生了变化 - mt19937产生整数，normal_distribution消耗双精度

c ++ 0x自动使用自适应，g ++ TR1代码不

为了让您的代码使用g ++ - 4.4和TR1，请执行以下操作

std::tr1::mt19937 prng(seed);
std::tr1::normal_distribution<double> normal;
std::tr1::variate_generator<std::tr1::mt19937, std::tr1::normal_distribution<double> > randn(prng,normal);
double r = randn();

Answer 2

这肯定不会挂起程序。但是，不确定它是否真的符合您的需求。

 #include <random>
 #include <iostream>

 using namespace std;

 typedef std::tr1::ranlux64_base_01 Myeng; 

 typedef std::tr1::normal_distribution<double> Mydist; 

 int main() 
 { 
      Myeng eng; 
      eng.seed(1000);
      Mydist dist(1,10); 

      dist.reset(); // discard any cached values 
      for (int i = 0; i < 10; i++)
      {
           std::cout << "a random value == " << (int)dist(eng) << std::endl; 
      }

 return (0); 
 }

Answer 3

如果您的TR1随机数生成实现有问题，您可以通过编写自己的普通生成器来避免TR1，如下所示。

使用您信任的任何随机生成器生成两个均匀（0,1）随机样本u和v。然后让r = sqrt（-2 log（u））并返回x = r sin（2 pi v）。（这被称为Box-Mueller方法。）

如果您需要具有平均μ和标准差sigma的正常样本样本，请返回sigma * x + mu而不是x。

Answer 4

虽然这似乎是一个错误，但快速确认将传递默认的0.0,1.0参数。 normal_distribution<double>::normal_distribution()应该等于normal_distribution<double>::normal_distribution(0.0, 1.0)

C ++ TR1：如何使用normal_distribution？

4 个答案: