[Math]服从高斯分布的随机生成器

最新推荐文章于 2022-02-14 20:23:39 发布

weixin_30500105

最新推荐文章于 2022-02-14 20:23:39 发布

阅读量170

点赞数

原文链接：http://www.cnblogs.com/jimmyqwy/archive/2010/10/25/1860857.html

版权

组里的同事突然问我如何生成服从高斯的随机生成器。

想当然的就感觉给定高斯函数，然后利用rand()的均匀分布产生的随机数，来模拟。

总感觉有点问题。

然后调查了一番，感觉还是有点东西需要记录，算是再学习吧。

1）按照中心极限定理，对于独立同分布的随机变量之和逼近于正态分布。

这样就可以用rand()生成的随机数，有限次相加就可以逼近正态分布了。

加上平均值和偏差也就可以模拟服从高斯分布的随机数了。

2）由于上述需要多次相加，计算量比较大。由此有另外一种方法就是

使用Box-Muller Transformation.

将两个独立生成的服从均匀分布的随机数U1，U2（0,1之间）。进行如下转换

$z_0 = \sqrt{-2 \ln U_1} \cos(2 \pi U_2) = \sqrt{-2 \ln s} \left(\frac{u}{\sqrt{s}}\right) = u \cdot \sqrt{\frac{-2 \ln s}{s}}$

$z_1 = \sqrt{-2 \ln U_1} \sin(2 \pi U_2) = \sqrt{-2 \ln s}\left( \frac{v}{\sqrt{s}}\right) = v \cdot \sqrt{\frac{-2 \ln s}{s}}.$

即服从两个独立的正态分布。

u,v为极坐标，s为平方和开根号（半径）

s与U1对应， $\scriptstyle \theta/(2 \pi)$ 与U2对应。

为啥经过这个变化就可以直接得到服从正态分布的数呢。

根据Marsaglia polar method。

If u is uniformly distributed in the interval 0 ≤ u < 1,
then the point (cos(2πu), sin(2πu)) is uniformly distributed on the unit circumference x² + y² = 1,
and multiplying that point by an independent random variable ρ whose distribution is

$\Pr(\rho<a)=\int_0^a re^{-r^2/2}\,dr$

will produce a point

$\left(\rho\cos(2\pi u),\rho\sin(2\pi u)\right)$

whose coordinates are jointly distributed as two independent standard normal random variables.

此处的u即为上述的服从均匀分布的U2

而rho，正好是上述U1的chi-2距离。正好服从上述分布

$\Pr(\rho<a)=\int_0^a re^{-r^2/2}\,dr$

那么就说明，产生的两个点z0,z1服从独立的正态分布。

而至于Marsaglia polar method方法的证明，我没有直接找到。

经过一些相关的资料，可以发现，使用Inverse transform sampling的概念，

生成点函数的反函数的Jacobian行列式为高斯密度函数，可以推导出此函数即服从于高斯分布。

（可能理解还有偏差，如有不对请高人指点）

具体证明请参照下一篇续-修改@2010.11.2

转载于:https://www.cnblogs.com/jimmyqwy/archive/2010/10/25/1860857.html

weixin_30500105

关注

0
点赞
踩
0

收藏

觉得还不错? 一键收藏
0
评论
复制链接

分享到 QQ

分享到新浪微博

扫一扫

评论

被折叠的条评论为什么被折叠?

到【灌水乐园】发言

查看更多评论

添加红包

成就一亿技术人!

hope_wisdom

发出的红包

实付元

使用余额支付

点击重新获取

扫码支付

钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。