poisson distから乱数を作成します。メルセンヌ・ツイスタ

Question

を使用してこんにちは、私は私が簡単に均一で、正規分布からの乱数を取得している私が走ってるのシミュレーションがあります。

#include <iostream> 
#include "MersenneTwister.h" 
using namespace std; 
int main() 
{ 
    MTRand mtrand1; 
    double r1,r2; 
    r1 = mtrand.rand(); // from a uninform dist. 
    r2 = mtrand1.randNorm(); //from a normal dist. 
} 

私はAから乱数を取得するために、この乱数ジェネレータを使用したいが平均 'A'を有するポアソン分布。

MersseneTwisterコードを使用してこの手順を実装する方法については何か考えてください。 コードは https://gcc.gnu.org/bugzilla/attachment.cgi?id=11960にあります。広く使用されています。

Answer 1

あなたはstd::random_deviceで播種するunlikely to be satisfactoryであることに注意してください標準ライブラリ

#include<random> 

double mean = 3.1415926; 
std::mt19937 mt{std::random_device{}()}; 
std::poisson_distribution<> pd{mean}; 

auto n = pd(mt); // get a number 

を使用することができます。