佳木斯湛栽影视文化发展公司

主頁 > 知識(shí)庫 > 從np.random.normal()到正態(tài)分布的擬合操作

從np.random.normal()到正態(tài)分布的擬合操作

熱門標(biāo)簽:服務(wù)器配置 智能手機(jī) 銀行業(yè)務(wù) 美圖手機(jī) 鐵路電話系統(tǒng) 網(wǎng)站文章發(fā)布 檢查注冊(cè)表項(xiàng) 呼叫中心市場(chǎng)需求

先看偉大的高斯分布(Gaussian Distribution)的概率密度函數(shù)(probability density function):

對(duì)應(yīng)于numpy中:

numpy.random.normal(loc=0.0, scale=1.0, size=None)

參數(shù)的意義為:

loc:float

此概率分布的均值(對(duì)應(yīng)著整個(gè)分布的中心centre)

scale:float

此概率分布的標(biāo)準(zhǔn)差(對(duì)應(yīng)于分布的寬度,scale越大越矮胖,scale越小,越瘦高)

size:int or tuple of ints

輸出的shape,默認(rèn)為None,只輸出一個(gè)值

我們更經(jīng)常會(huì)用到的np.random.randn(size)所謂標(biāo)準(zhǔn)正態(tài)分布

對(duì)應(yīng)于np.random.normal(loc=0, scale=1, size)。

采樣(sampling)

# 從某一分布(由均值和標(biāo)準(zhǔn)差標(biāo)識(shí))中獲得樣本
mu, sigma = 0, .1
s = np.random.normal(loc=mu, scale=sigma, size=1000)

也可使用scipy庫中的相關(guān)api(這里的類與函數(shù)更符合數(shù)理統(tǒng)計(jì)中的直覺):

import scipy.stats as st
mu, sigma = 0, .1
s = st.norm(mu, sigma).rvs(1000)

校驗(yàn)均值和方差:

>>> abs(mu  np.mean(s))  .01
True
>>> abs(sigma-np.std(s, ddof=1))  .01
True
            # ddof,delta degrees of freedom,表示自由度
            # 一般取1,表示無偏估計(jì),

擬合

我們看使用matplotlib.pyplot便捷而強(qiáng)大的語法如何進(jìn)行高斯分布的擬合:

import matplotlib.pyplot as plt
count, bins, _ = plt.hist(s, 30, normed=True)
        # normed是進(jìn)行擬合的關(guān)鍵
        # count統(tǒng)計(jì)某一bin出現(xiàn)的次數(shù),在Normed為True時(shí),可能其值會(huì)略有不同
plt.plot(bins, 1./(np.sqrt(2*np.pi)*sigma)*np.exp(-(bins-mu)**2/(2*sigma**2), lw=2, c='r')
plt.show()

或者:

s_fit = np.linspace(s.min(), s.max())
plt.plot(s_fit, st.norm(mu, sigma).pdf(s_fit), lw=2, c='r')

np.random.normal()的含義及實(shí)例

這是個(gè)隨機(jī)產(chǎn)生正態(tài)分布的函數(shù)。(normal 表正態(tài))

先看一下官方解釋:

有三個(gè)參數(shù)

loc:正態(tài)分布的均值,對(duì)應(yīng)著這個(gè)分布的中心.代表下圖的μ

scale:正態(tài)分布的標(biāo)準(zhǔn)差,對(duì)應(yīng)分布的寬度,scale越大,正態(tài)分布的曲線 越矮胖,scale越小,曲線越高瘦。 代表下圖的σ

size:你輸入數(shù)據(jù)的shape,例子:

下面展示一些 內(nèi)聯(lián)代碼片。

// An highlighted block
a=np.random.normal(0, 1, (2, 4))
print(a)
輸出:
[[-0.29217334  0.41371571  1.26816017  0.46474676]
 [ 1.33271487  0.80162296  0.47974157 -1.49748788]]

看這個(gè)圖直觀些:

以下為官方文檔:

以上為個(gè)人經(jīng)驗(yàn),希望能給大家一個(gè)參考,也希望大家多多支持腳本之家。

您可能感興趣的文章:
  • python numpy之np.random的隨機(jī)數(shù)函數(shù)使用介紹
  • Numpy中np.random.rand()和np.random.randn() 用法和區(qū)別詳解
  • Numpy之random函數(shù)使用學(xué)習(xí)

標(biāo)簽:河南 上海 新疆 長(zhǎng)治 紅河 沈陽 樂山 滄州

巨人網(wǎng)絡(luò)通訊聲明:本文標(biāo)題《從np.random.normal()到正態(tài)分布的擬合操作》,本文關(guān)鍵詞  ;如發(fā)現(xiàn)本文內(nèi)容存在版權(quán)問題,煩請(qǐng)?zhí)峁┫嚓P(guān)信息告之我們,我們將及時(shí)溝通與處理。本站內(nèi)容系統(tǒng)采集于網(wǎng)絡(luò),涉及言論、版權(quán)與本站無關(guān)。
  • 相關(guān)文章
  • 收縮
    • 微信客服
    • 微信二維碼
    • 電話咨詢

    • 400-1100-266
    房产| 博兴县| 兰考县| 遂川县| 南安市| 东乡县| 阿克苏市| 金堂县| 阳东县| 英超| 济阳县| 收藏| 黑水县| 临江市| 卫辉市| 紫阳县| 六安市| 旺苍县| 惠水县| 泾源县| 淳化县| 阿荣旗| 山东| 东丰县| 忻州市| 通榆县| 萝北县| 法库县| 海林市| 海原县| 汝城县| 奉化市| 离岛区| 汕尾市| 舟山市| 西和县| 定远县| 隆昌县| 屏东市| 突泉县| 宿迁市|