perl : 正規分布の乱数を得る | luxの日記 | スラド
正規乱数(正規分布の乱数)を作る方法です。 ボックス・ミューラーの正規乱数発生法より $sigma は標準偏差、$m は、中央値です。 3.14159265359 はπ(円周率)です。 sub randn { my ($m, $sigma) = @_; my ($r1, $r2) = (rand(), rand()); while ($r1 == 0) { $r1 = rand(); } return ($sigma * sqrt(-2 * log($r1)) * sin(2 * 3.14159265359 * $r2)) + $m; } print &randn(5, 3); # 平均値5、標準偏差3
■ - 泥船
N($mean, $deviation^2)なる正規分布に従った乱数生成。一様乱数はメルセンヌツイスタ。 任意の確率分布の乱数をしたい場合は、、、確か逆関数法とか棄却法とかで頑張るんだった気が。 use strict; use warnings; use Math::Random::MT qw(rand srand); use Math::Trig; srand(1); my $mean = 10.0; my $deviation = 1.0; my $numData = 10000; my @data; for(my $i = 0; $i < $numData; $i++){ my $rnd = normalRnd($mean, $deviation); push(@data, $rnd); print "$rnd\n"; } sub normalRnd{ my ($mean, $dev
1
リリース、障害情報などのサービスのお知らせ
最新の人気エントリーの配信
処理を実行中です
j次のブックマーク
k前のブックマーク
lあとで読む
eコメント一覧を開く
oページを開く
Math-Random-MT-1.17 - The Mersenne Twister PRNG - metacpan.org
