ボックス＝ミューラー法

よしコクのある乱数の実際のサンプルをアップした。こうやって並べてみると、味の違いは一目瞭然。 https://t.co/esbYaYWYLh

— 深津 貴之 (@fladdict) November 4, 2016

ボックスミューラー法の養殖された乱数とはｗｗ？

10000.times {
  x = rand()
  y = rand()

  z1 = Math.sqrt(-2 * Math.log(x)) * Math.cos(2 * Math::PI * y)
  z2 = Math.sqrt(-2 * Math.log(x)) * Math.sin(2 * Math::PI * y)

  puts "#{z1}\t#{z2}"
}

ruby bm.rb | pbcopy 

googleスプレッドシートに貼り付けてヒストグラムのグラフを作る（重い…）

rand - Google スプレッドシート

1万回程度ではあんまり綺麗なグラフにならないすね

---

おまけ

• Ruby gem → GitHub - s-osa/random_bell
• Javaクラスライブラリ → Random (Java Platform SE 8 )
• npm → gaussian

乱数にコクを出す？

togetter.com

コクとはｗｗ？

count = (ARGV[0] || 1).to_i

10000.times do
  puts Array.new(count) { rand() }.inject(&:+) / count 
end

rand.rbとして保存

$ ruby rand.rb 1 | pbcopy

これをgoogleスプレッドシートでグラフにする

rand - Google スプレッドシート

なるほど正規分布っぽくなる

たとえば分子のエネルギーなどは正規分布なので自然界の挙動をエミュレートするには単純な乱数よりいいかもしれませんね

---

よしコクのある乱数の実際のサンプルをアップした。こうやって並べてみると、味の違いは一目瞭然。 https://t.co/esbYaYWYLh

— 深津 貴之 (@fladdict) November 4, 2016

なるほどなるほど

中心極限定理 - Wikipedia