服從Gaussian分布的近似隨機數(shù)生成算法

中心極限定理是概率論中的一組定理,研究的是大量相互獨立的隨機變量之和在什么樣的條件下會收斂于正態(tài)分布。我們利用林德伯格-列維(Lindeberg-Levy)中心極限定理,可以通過服從均勻分布的隨機數(shù)生成服從高斯分布的隨機數(shù)生成算法

figu.png

林德伯格-列維中心極限定理

具有有限期望和方差(不為0)的獨立同分布的隨機變量之和,經(jīng)過標(biāo)準(zhǔn)化后,其分布收斂于標(biāo)準(zhǔn)正態(tài)分布

令隨機變量 X_1, X_2, … , X_n 獨立同分布,且具有有限的期望和方差 E(X_i) = \mu, D(X_i) = \sigma^2 \neq 0 \quad (i = 1,2, .., n)


Y = \frac{ \sum_{i=1}^{n} X_i - n\mu }{\sqrt{n} \sigma } \tag{1}

\lim_{n \to \infty} P(Y \leq a) = \Phi(a)
其中\Phi(a)是標(biāo)準(zhǔn)正態(tài)分布的分布函數(shù)

生成服從Gaussian分布的近似隨機數(shù)

標(biāo)準(zhǔn)正態(tài)分布

根據(jù)上文,我們可以得知,取多次指定分布下的隨機數(shù),然后經(jīng)過求和、標(biāo)準(zhǔn)化后,即可使得該隨機數(shù)滿足標(biāo)準(zhǔn)正態(tài)分布

這里我們?nèi)∪菀撰@得的均勻分布X_n \sim U(0,1)為例,則易知:

E(X_i) = 1/2, D(X_i) = 1/12

則根據(jù)(1)式可得,所需的服從標(biāo)準(zhǔn)正態(tài)分布的隨機變量Y:
Y = \sqrt(\frac{12}{n}) (\sum_{i=1}^{n}X_n - \frac{n}{2}) \sim N(0,1)
在有些算法的實現(xiàn)源碼可以看到n取值12,其原因就在這里根號正好可以消去。n越大,其精度越高

普通高斯分布

一般高斯分布Z \sim N( \mu_2, \sigma_2 )和標(biāo)準(zhǔn)正態(tài)分布Y \sim N(0, 1)的轉(zhuǎn)換公式如下:
Y = \frac{Z-\mu_2}{\sigma_2} \tag{2}
則可以用(2)式獲得指定的高斯分布

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
【社區(qū)內(nèi)容提示】社區(qū)部分內(nèi)容疑似由AI輔助生成,瀏覽時請結(jié)合常識與多方信息審慎甄別。
平臺聲明:文章內(nèi)容(如有圖片或視頻亦包括在內(nèi))由作者上傳并發(fā)布,文章內(nèi)容僅代表作者本人觀點,簡書系信息發(fā)布平臺,僅提供信息存儲服務(wù)。
禁止轉(zhuǎn)載,如需轉(zhuǎn)載請通過簡信或評論聯(lián)系作者。

相關(guān)閱讀更多精彩內(nèi)容

友情鏈接更多精彩內(nèi)容