Statistics\(U\)を一様分布\(U(0,1)\)に従う乱数とする。このとき、ある確率密度関数\(f(x)\)に対する累積分布関数\(F(x)\)の逆関数\(F^{-1}\)を用いて、乱数\(X\)を\(X=F^{-1}(U)\)と定義すると\(X\)は分布\(f(x)\)に従う。0,\qquad\int_{-\infty}^{+\infty}f(t)\mathrm{d}t=1\tag{1}\]を満たしていて、その累積分布関数\(F(x)\)は次のようになる。x]=F(x)=\int_{-\infty}^xf(t)\mathrm{d}t\tag{2}\]この定義による関数\(F(x)\)は必ず単調増加関数になるので累積分布関数には必ず逆関数\(F^{-1}\)がある。これを踏まえて、一様分布\(U(0,1)\)に従う確率変数\(U\)と何らかの累積分布関数の逆関数\(F^{-1}\)によって定義される確率変数\(X=F^{-1}(U)\)がどのような分布に従うかを計算してみると次のようになる。b)\)は適当な実数、\(f_u(u)\)は一様分布\(U(0,1)\)の確率密度関数で\(u\in[0,1]\)のとき\(f_u(u)=1\)、それ以外の時\(f_u(u)=0\)になる。累積分布関数は\([0,1]\)の範囲の値しかとらないからこのように積分する区間を変えずに変形できる。次に\(u=F(t)\)と変数変換する。\int_{a}^{b}f(t)\mathrm{d}t\tag{4}\]となって、これは\(X\)が確率密度分布\(f(x)\)で表される分布に従うことの定義そのもの、つまり\(X\)は分布\(f(x)\)に従うことが示された。0)\tag{5}\]に従う乱数を生成する。x}\tag{6}\]となる。\log(1-u)\tag{7}\]より、一様分布に従う乱数\(U\)として\(X=-(1/\lambda)\log(1-U)\)とすれば\(X\)は指数分布\(f_\lambda(x)\)に従う。U\)とした方が僅かに効率がよい。PostAbout個人的メモサイト的電子ノート。誰かに見られる可能性のある場所に置いておく事に意味があることもある。実際に誰か見ているかどうかはそんなに気にしないし後から見たら間違ったことも書いている。ノートだし消すこともないかということで置いたまま悪しからず。JsonSitemapXMLContact万が一何か言いたいことがあればコチラまで。