資源描述:
《matlab統(tǒng)計(jì)筆記》由會(huì)員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在應(yīng)用文檔-天天文庫(kù)。
1、X=[-3:0.5:3];f=normpdf(x,0,1);(其中normpdf為正態(tài)分布的Matlab分布實(shí)現(xiàn)函數(shù),可由以下介紹的函數(shù)代替。)功能:可選分布的概率密度函數(shù)。格式:Y=pdf(‘name’,X,A1,A2,A3)說(shuō)明:‘name’為特定分布的名稱,如‘Normal’,’Gamma’等。X為分布函數(shù)的自變量X的取值矩陣,而A1,A2,A3分別為相應(yīng)的分布參數(shù)值。Y給出結(jié)果,為概率密度值矩陣。舉例:p=pdf(‘Normal’,-2:2,0,1)給出標(biāo)準(zhǔn)正態(tài)分布在-2到2的分布函數(shù)值。
2、而p=pdf(‘Poisson’,0:4,1:5)給出Poisson分布函數(shù)。---random功能:產(chǎn)生可選分布的隨機(jī)數(shù)。格式:y=random(‘name’,A1,A2,A3,m,n)說(shuō)明:random函數(shù)產(chǎn)生統(tǒng)計(jì)工具箱中任一分布的隨機(jī)數(shù)?!畁ame’為相應(yīng)分布的名稱。A1,A2,A3為分布參數(shù),意義同pdf參數(shù)。m,n確定了結(jié)果y的數(shù)量,如果分布參數(shù)A1,A2,A3為矢量,則m,n是可選的,但應(yīng)注意,它們給出的長(zhǎng)度或矩陣行列數(shù)必須與分布參數(shù)的長(zhǎng)度相匹配。舉例:rn=random(‘Norma
3、l’,0,1,2,4)基本思想:已知一組觀測(cè)值,給定這組值出自的某類分布中,求得最有可能出現(xiàn)這組值的一個(gè)分布.調(diào)用方法:[phat,pci]=mls[‘dist’,data,alpha]phat為參數(shù)估計(jì)結(jié)果,pci為置信區(qū)間計(jì)算結(jié)果dist為用戶給定的分布名稱,data為數(shù)據(jù)列表,(1-alpha)置信區(qū)域.正態(tài)分布的負(fù)對(duì)數(shù)似然函數(shù)調(diào)用方法L=normlike(params,data)Params為正態(tài)分布參數(shù):params(1)為μ,params(2)為σβ分布的負(fù)對(duì)數(shù)似然函數(shù)logL=bet
4、alike(params,data)γ分布的負(fù)對(duì)數(shù)似然函數(shù)logL=gamlike(params,data)威布爾分布的負(fù)對(duì)數(shù)似然函數(shù)logL=weiblike(params,data)參數(shù)設(shè)置與正態(tài)分布的負(fù)對(duì)數(shù)似然函數(shù)類似,不加冗述.?功能:樣本的幾何均值。?格式:m=geomean(X)?說(shuō)明:幾何均值的定義為m=(1.4.1)geomean函數(shù)計(jì)算樣本的幾何均值。X若為矢量,它返回X中元素的幾何均值;X若為矩陣,它的結(jié)果為一個(gè)行矢量,每個(gè)元素為X對(duì)應(yīng)列元素的幾何均值。舉例:x=exprnd(
5、1,10,6);geometric=geomean(X);average=mean(X);?樣本數(shù)據(jù)的調(diào)和均值。?格式:m=harmmean(X)Median樣本數(shù)據(jù)的中值?剔除極端數(shù)據(jù)的樣本均值。?格式:m=trimmean(X,percent)?說(shuō)明:函數(shù)計(jì)算剔除觀測(cè)量中最高百分比和最低百分比數(shù)據(jù)后的均值。函數(shù)中percent代表百分比。在Matlab里,有關(guān)散布度量計(jì)算的函數(shù)為:1:計(jì)算樣本的內(nèi)四分位數(shù)間距的iqr(X).2:求樣本數(shù)據(jù)的平均絕對(duì)偏差的mad(X).3:計(jì)算樣本極差的rang
6、e(X).4:計(jì)算樣本方差的var(X,w).5:求樣本的標(biāo)準(zhǔn)差的std(X).6:求協(xié)方差矩陣的cov(X).這些函數(shù)的詳細(xì)說(shuō)明可以參見Matlab的幫助文檔。