Matlab累積分布函數(shù)cdf與概率密度函數(shù)pdf

原文:http://blog.sina.com.cn/s/blog_7054a1960102vy7x.html

累積分布函數(shù)cdf (Cumulative Distribution Function)

背景知識:http://www.lifelaf.com/blog/?p=746

語法

y = cdf('name',x,A)

y = cdf('name',x,A,B)

y = cdf('name',x,A,B,C)

y = cdf(pd,x)

y = cdf(___,'upper')

描述

y = cdf('name',x,A)?計算某種分布(由'name'定義衷掷,如'Normal'正態(tài), 'Poisson'泊松, 'T' t分布…)下,x值處的累計分布基显,A,B,C等為'name'函數(shù)的參數(shù)

y = cdf(pd,x)?直接計算概率分布函數(shù)pd(probability distribution)?脱衙,在x處的累計分布侥猬,實際上,這里的pd?已被'name', A定義好捐韩,舉栗如下:

%?定義一個正態(tài)分布函數(shù)pd,?均值mu = 0,?標準差sigma = 1.

mu = 0;

sigma = 1;

pd = makedist('Normal',mu,sigma);

%?定義x

x = [-2,-1,0,1,2];

%?計算x值處的累計分布

y = cdf(pd,x)

y =

?0.0228?0.1587?0.5000?0.8413?0.9772

用第一種語句表達相同內(nèi)容為:

y2 = cdf('Normal',x,mu,sigma)?%正態(tài)分布退唠,x值處,均值為0荤胁,標準差為1

y2 =

?0.0228?0.1587?0.5000?0.8413?0.9772

http://nl.mathworks.com/help/stats/cdf.html

t分布累積分布函數(shù)tcdf (Student'stcumulative distribution function)

%?事實上就是y = cdf('T',x,A)函數(shù)

語法

p = tcdf(x,nu)

p = tcdf(x,nu,'upper')

描述

計算t分布在x值處的累積分布瞧预,nu是t分布的自由度

再舉個栗子

mu = 1;?% Population mean

sigma = 2;% Population standard deviation

n = 100;?% Sample size

x = normrnd(mu,sigma,n,1);% Random sample from population

xbar = mean(x);% Sample mean

s = std(x);?% Sample standard deviation

t = (xbar - mu)/(s/sqrt(n)) %?這里t分布出現(xiàn)了,正態(tài)分布總體與樣本均值的差符合t分布

t =

?1.0589

p = 1-tcdf(t,n-1)?% Probability of larger t-statistic

p =

?0.1461

該p值(即t函數(shù)的累積分布就是t檢驗在相同x值處的概率ptest)

[h,ptest] = ttest(x,mu,0.05,'right')

h =

?0

ptest =

?0.1461

http://nl.mathworks.com/help/stats/tcdf.html

概率密度函數(shù)pdf (Probability density functions)

搞懂了累積分布函數(shù)cdf仅政,這個就沒什么需要多說了

語法

y = pdf('name',x,A)

y = pdf('name',x,A,B)

y = pdf('name',x,A,B,C)

y = pdf(pd,x)

舉例

%?定義一個正態(tài)分布函數(shù)pd,?均值mu = 0,?標準差sigma = 1.

mu = 0;

sigma = 1;

pd = makedist('Normal',mu,sigma);

%?定義x值

x = [-2 -1 0 1 2];

%?計算x值處的概率密度(cdf是累計分布)

y = pdf (pd,x)

y =

0.0540?0.2420?0.3989?0.2420?0.0540

同樣垢油,另一種表達

y = pdf(pd,x)

y =

?0.0540?0.2420?0.3989?0.2420?0.0540

http://nl.mathworks.com/help/stats/pdf.html

t分布概率密度函數(shù)tpdfStudent's t probability density function

語法

y = tpdf(x,nu)

舉例

tpdf(0,1:6)

ans =

?0.3183?0.3536?0.3676?0.3750?0.3796?0.3827

http://nl.mathworks.com/help/stats/tpdf.html

相反,還可以通過p求t分布的t值

tinv (Student's t inverse cumulative distribution function)

語法

x = tinv(p,nu)

舉例

% the 99th percentile of the Student's t distribution for one to six degrees of freedom

percentile = tinv(0.99,1:6)

percentile =

?31.8205?6.9646?4.5407?3.7469?3.3649?3.1427

http://nl.mathworks.com/help/stats/tinv.html

有一個問題圆丹,Matlab有一個inv矩陣求逆函數(shù)滩愁,不知與tinv什么關系,莫非tinv是在t分布下調(diào)用了inv計算程序辫封?但p并不等是t的逆矩陣啊(即t*p = E)跋跬鳌?求解答

inv是矩陣求逆的意思倦微。具體用法A=inv(B)妻味,其中B是輸入的可逆矩陣,輸出A就是B的逆矩陣璃诀,逆矩陣滿足性質(zhì) AB=BA=E (E是單位陣)弧可。如果輸入的是不可逆矩陣會彈出警告,并返回inf。

調(diào)用舉例:

>> inv([1 0;0 0])

警告: 矩陣為奇異工作精度棕诵。

ans =

Inf Inf

Inf Inf

>> inv(rand(2))

ans =

-13.0929 5.2640

12.0501 -3.3159

另附官方英文解釋(輸入doc inv也可以自己查看):

Y = inv(X) returns theinverse of the square matrix X. A warning messageis printed if X is badly scaled or nearly singular.

In practice, it is seldom necessary to form the explicit inverseof a matrix. A frequent misuse of inv arises whensolving the system of linear equations Ax = b.One way to solve this is with x = inv(A)*b.A better way, from both an execution time and numerical accuracy standpoint,is to use the matrix division operator x = A\b.This produces the solution using Gaussian elimination, without formingthe inverse. See mldivide (\)for further information.

http://zhidao.baidu.com/link?url=dTjjN7fsj2EBX7zwCm1_TZA2LDv4Abtmgnq0mwfb3pyCLpZ26g1RWFGFlKOiADzVsnpeUh-bu4o9QvP_e5T5q

?著作權歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末裁良,一起剝皮案震驚了整個濱河市,隨后出現(xiàn)的幾起案子校套,更是在濱河造成了極大的恐慌价脾,老刑警劉巖,帶你破解...
    沈念sama閱讀 221,548評論 6 515
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件笛匙,死亡現(xiàn)場離奇詭異侨把,居然都是意外死亡,警方通過查閱死者的電腦和手機妹孙,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 94,497評論 3 399
  • 文/潘曉璐 我一進店門秋柄,熙熙樓的掌柜王于貴愁眉苦臉地迎上來,“玉大人蠢正,你說我怎么就攤上這事骇笔。” “怎么了嚣崭?”我有些...
    開封第一講書人閱讀 167,990評論 0 360
  • 文/不壞的土叔 我叫張陵笨触,是天一觀的道長。 經(jīng)常有香客問我雹舀,道長芦劣,這世上最難降的妖魔是什么? 我笑而不...
    開封第一講書人閱讀 59,618評論 1 296
  • 正文 為了忘掉前任说榆,我火速辦了婚禮虚吟,結(jié)果婚禮上,老公的妹妹穿的比我還像新娘娱俺。我一直安慰自己稍味,他們只是感情好,可當我...
    茶點故事閱讀 68,618評論 6 397
  • 文/花漫 我一把揭開白布荠卷。 她就那樣靜靜地躺著模庐,像睡著了一般。 火紅的嫁衣襯著肌膚如雪油宜。 梳的紋絲不亂的頭發(fā)上掂碱,一...
    開封第一講書人閱讀 52,246評論 1 308
  • 那天,我揣著相機與錄音慎冤,去河邊找鬼疼燥。 笑死,一個胖子當著我的面吹牛蚁堤,可吹牛的內(nèi)容都是我干的醉者。 我是一名探鬼主播绵咱,決...
    沈念sama閱讀 40,819評論 3 421
  • 文/蒼蘭香墨 我猛地睜開眼奢浑,長吁一口氣:“原來是場噩夢啊……” “哼!你這毒婦竟也來了?” 一聲冷哼從身側(cè)響起柳爽,我...
    開封第一講書人閱讀 39,725評論 0 276
  • 序言:老撾萬榮一對情侶失蹤惭每,失蹤者是張志新(化名)和其女友劉穎损谦,沒想到半個月后朦乏,有當?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體,經(jīng)...
    沈念sama閱讀 46,268評論 1 320
  • 正文 獨居荒郊野嶺守林人離奇死亡粒竖,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點故事閱讀 38,356評論 3 340
  • 正文 我和宋清朗相戀三年颅崩,在試婚紗的時候發(fā)現(xiàn)自己被綠了。 大學時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片蕊苗。...
    茶點故事閱讀 40,488評論 1 352
  • 序言:一個原本活蹦亂跳的男人離奇死亡沿后,死狀恐怖,靈堂內(nèi)的尸體忽然破棺而出朽砰,到底是詐尸還是另有隱情得运,我是刑警寧澤,帶...
    沈念sama閱讀 36,181評論 5 350
  • 正文 年R本政府宣布锅移,位于F島的核電站,受9級特大地震影響饱搏,放射性物質(zhì)發(fā)生泄漏非剃。R本人自食惡果不足惜,卻給世界環(huán)境...
    茶點故事閱讀 41,862評論 3 333
  • 文/蒙蒙 一推沸、第九天 我趴在偏房一處隱蔽的房頂上張望备绽。 院中可真熱鬧,春花似錦鬓催、人聲如沸肺素。這莊子的主人今日做“春日...
    開封第一講書人閱讀 32,331評論 0 24
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽倍靡。三九已至,卻和暖如春课舍,著一層夾襖步出監(jiān)牢的瞬間塌西,已是汗流浹背。 一陣腳步聲響...
    開封第一講書人閱讀 33,445評論 1 272
  • 我被黑心中介騙來泰國打工筝尾, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留捡需,地道東北人。 一個月前我還...
    沈念sama閱讀 48,897評論 3 376
  • 正文 我出身青樓筹淫,卻偏偏與公主長得像站辉,于是被迫代替她去往敵國和親。 傳聞我的和親對象是個殘疾皇子,可洞房花燭夜當晚...
    茶點故事閱讀 45,500評論 2 359

推薦閱讀更多精彩內(nèi)容