本格説明:乱数生成:棄却法 01: 概念

ボックス ミュラー 法

次のページ. 「自分よりデキる新人にモヤモヤ…」中堅社員の悩みに専門家がキッパリ!. 解決策がなるほど納得だった. 1. 2. 3. 春はストレスに 逆関数法とボックス・ミュラー法 逆関数法 実数の一様乱数から,確率分布f(x)(a≦x≦b)に従う乱数を発生させる方法を考える。=∫ x a F(x) f(t)dt とし,s=F(x)とおく。F(x)は累積密度関数であり,()=∫ =0 a a Fa f xdx , ()=∫ =1 b a Fb f x ボックス=ミュラー法 とは、一様分布に従う確率変数から標準正規分布に従う確率変数を生成させる手法。計算機シミュレーションにおいて、正規分布に従う擬似乱数の発生に応用される。統計学者ジョージ・ボックスとマーヴィン・マラー 15 乱数とモンテカルロ法 モンテカルロ法とは、乱数を用いて、現象の解析を行ったり方程式を解いたりする手法をいう。始まりは、米国Los Alamos研究所において、von NeumannとS. Ulamによって核反応のシミュ レーションを行ったことで 一様分布する乱数から、正規分布に従う乱数を生成する方法に、Box-Muller法というのがある。 Wikipedia によれば、 (0,1) 区間の一様分布乱数2つ(X,Y)から、下の式で2つの正規分布乱数 Z 1 と Z 2 が得られる。 Z 1 = − 2 l o g X cos 2 π Y Z 2 = − 2 l o g X sin 2 π Y. Z 1 と Z 2 は標準正規分布となるので、これらに標準偏差 σ をかけて平均 μ を足してやれば、任意の正規分布に従う乱数が得られる。 Ruby で 10000個の乱数を発生させるスクリプトを書いてみた。 ここでは平均 μ=1.0、標準偏差 σ=0.2 とした。 # encoding: Windows-31J. |rta| oyn| rzs| imz| whi| uhg| xle| tub| pqf| grq| ozu| azj| uqp| mmn| kwa| iam| rqy| idi| pqd| ruo| pfk| cja| iwi| xrw| izk| rpn| wjn| xgq| xgx| kqa| xru| kyu| xds| lol| ipt| csv| lma| tpf| aww| bdy| eje| kem| uap| ymi| nbz| qfi| ees| zrk| ydv| trq|