Cのポアソン分布の計算

Question

kの値が720までのポアソン分布を計算するには、C関数が必要です。効率の良いソリューションが必要です。/K - （ラムダ）gsl_ran_poisson_pdf

Answer 1

はGSLしようとする与えます！
あなたは私はかなり遅れてあなたの緊急用思うよログ（N！）のログ（N）+ログ=（N-1を！）と動的計画

Answer 2

をライブラリ を使用するのではなく、それを自分で計算したい場合は、式を使用して、それを計算することができます。.. E^k個の*電子^：

Answer 3

ポアソン乱数生成器

int poissonRandom(double expectedValue) { 
    int n = 0; //counter of iteration 
    double limit; 
    double x; //pseudo random number 
    limit = exp(-expectedValue); 
    x = rand()/INT_MAX; 
    while (x > limit) { 
    n++; 
    x *= rand()/INT_MAX; 
    } 
    return n; 
} 

を使用することができますデマンド。

Answer 4

これは元々のリクエストでは遅すぎると思いますが、答えのいくつかが間違っていると思います。配布から乱数を生成したいとは思っていませんが、配布自体が必要です。ここでは、これを行う関数が大きくなる可能性のある階乗の計算を避ける です。

double poisson(int k, double mean) { 
    double p = std::exp(-mean); 
    double f = 1; 
    for (int i=0 ; i<k ; i++) f *= mean/(i+1);  
    return p*f; 
}