在 Ruby 中生成高斯（正态分布）随机数的代码

在 ruby​​ 中生成正态分布随机数的代码是什么？

（注：我回答了我自己的问题，但我会等几天再接受，看看是否有人有更好的答案。）

EDIT:

为了搜索这个，我查看了两次搜索产生的 SO 上的所有页面：

+“正态分布”红宝石

and

+高斯 +随机红宝石

蟒蛇的随机.高斯() http://www.opensource.apple.com/source/python/python-3/python/Lib/random.py和升压的正态分布 http://www.boost.org/doc/libs/1_43_0/boost/random/normal_distribution.hpp两者都使用Box-Muller 变换 http://en.wikipedia.org/wiki/Box%E2%80%93Muller_transform，所以这对于 Ruby 来说也应该足够好了。

def gaussian(mean, stddev, rand)
  theta = 2 * Math::PI * rand.call
  rho = Math.sqrt(-2 * Math.log(1 - rand.call))
  scale = stddev * rho
  x = mean + scale * Math.cos(theta)
  y = mean + scale * Math.sin(theta)
  return x, y
end

该方法可以封装在一个类中，该类一一返回样本。

class RandomGaussian
  def initialize(mean, stddev, rand_helper = lambda { Kernel.rand })
    @rand_helper = rand_helper
    @mean = mean
    @stddev = stddev
    @valid = false
    @next = 0
  end

  def rand
    if @valid then
      @valid = false
      return @next
    else
      @valid = true
      x, y = self.class.gaussian(@mean, @stddev, @rand_helper)
      @next = y
      return x
    end
  end

  private
  def self.gaussian(mean, stddev, rand)
    theta = 2 * Math::PI * rand.call
    rho = Math.sqrt(-2 * Math.log(1 - rand.call))
    scale = stddev * rho
    x = mean + scale * Math.cos(theta)
    y = mean + scale * Math.sin(theta)
    return x, y
  end
end

(CC0) http://creativecommons.org/publicdomain/zero/1.0/

在法律允许的范围内，安托纳科斯 https://stackoverflow.com/users/513598/antonakos已放弃所有版权以及相关或邻接权RandomGaussian红宝石类。本作品出版自：丹麦。