Luk*_*ser 3 random numpy scipy
我试图使用numpy/scipy从对数正态分布创建随机数.
平均值为2000,西格玛为800.
如果我使用numpy.random.lognormal(mean = 2000,sigma = 800,size = 10000)创建我的随机值,我得到的是非常高或inf数字.
有办法解决这个问题吗？
注意:mean和sigma参数对应于分发日志的lognormal分布; 分布的实际算术平均值是exp(mean + sigma**2/2),inf当mean=2000和时,它在标准双精度浮点数时计算sigma=800.
见 http://docs.scipy.org/doc/numpy/reference/generated/numpy.random.lognormal.html#numpy.random.lognormal 和 http://en.wikipedia.org/wiki/Log-normal_distribution 更多细节.