我現(xiàn)在是在英國留學(xué),有一門課要用到MATLAB,我現(xiàn)在有個(gè)問題想讓老師輔導(dǎo)一下,就是如何生成隨機(jī)數(shù),計(jì)算和繪制pdf,然后使用概率分布對(duì)象計(jì)算多項(xiàng)式分布的描述性統(tǒng)計(jì)?
同學(xué)提到的這個(gè)問題是概率分布中的離散分布問題,也就是將概率分布擬合到樣本數(shù)據(jù),評(píng)估概率函數(shù)(如pdf和cdf),計(jì)算匯總統(tǒng)計(jì)數(shù)據(jù)(如平均值和中值),可視化樣本數(shù)據(jù),生成隨機(jī)數(shù)。具體步驟如下:
1、定義分布參數(shù)
創(chuàng)建一個(gè)向量p,包含每個(gè)結(jié)果的概率。結(jié)果1的概率是1/2,結(jié)果2的概率是1/3,結(jié)果3的概率是1/6。每個(gè)實(shí)驗(yàn)的試驗(yàn)次數(shù)n為5次,實(shí)驗(yàn)的重復(fù)次數(shù)為8次。
p = [1/2 1/3 1/6];
n = 5;
reps = 8;
2、創(chuàng)建一個(gè)多項(xiàng)式概率分布對(duì)象
使用為probability參數(shù)指定的值p創(chuàng)建一個(gè)多項(xiàng)概率分布對(duì)象。
pd = makedist('Multinomial','Probabilities',p)
3、生成一個(gè)隨機(jī)數(shù)
從多項(xiàng)式分布中生成一個(gè)隨機(jī)數(shù),這是單次試驗(yàn)的結(jié)果。
rng('default') % For reproducibility
r = random(pd)
4、生成隨機(jī)數(shù)矩陣
可以從多項(xiàng)分布中生成一個(gè)隨機(jī)數(shù)矩陣,該矩陣報(bào)告多個(gè)實(shí)驗(yàn)的結(jié)果,每個(gè)實(shí)驗(yàn)都包含多個(gè)試驗(yàn)。生成一個(gè)矩陣,包含n = 5次試驗(yàn)和reps = 8次重復(fù)的實(shí)驗(yàn)結(jié)果。
r = random(pd,reps,n)
5、計(jì)算并繪制pdf
計(jì)算分布的pdf。
x = 1:3;
y = pdf(pd,x);
bar(x,y)
xlabel('Outcome')
ylabel('Probability Mass')
title('Trinomial Distribution')
6、計(jì)算描述性統(tǒng)計(jì)
計(jì)算分布的平均值、中位數(shù)和標(biāo)準(zhǔn)差。
m = mean(pd)
med = median(pd)
s = std(pd)
同學(xué)如果還有其他MATLAB相關(guān)問題需要老師輔導(dǎo)的話,可以直接聯(lián)系英國課程輔導(dǎo)老師,我們的老師非常熟悉這方面的知識(shí),同學(xué)只管提問就行。