c - 一様分布からポアソン分布への変換

Question

乱数の一様分布 (0 から 1 など) をポアソン分布に変換する C プログラムを作成する必要があります。誰でも助けることができますか？

Answer 1

Gnu Scientific Library である GSL を使用します。gsl_ran_poissonという関数があります:

この関数は、平均 mu のポアソン分布からランダムな整数を返します。ポアソン変量の確率分布は、p(k) = {\mu^k \over k!} \exp(-\mu) for k >= 0 です。

それ以外の場合は、コードを見てアイデアをコピーしてください。

Answer 2

U(0,1) の乱数が与えられた場合に、ポアソン分布から乱数をサンプリングできる C プログラムを書きたいと想定しています。

通常、これは U(0,1) から数値の逆 CDF を取得することによって行われます。ポアソンのような離散分布の場合、最初に CDF 関数が整数点間で滑らかであると仮定して連続分布に変換し、次に適切な近似 (フロア関数) を適用します。

Numerical Recipes in C++ (第 3 版)という本には、完全な説明と C++ コードもあります。秒 7.3.12、372 ページ。