matlab-第10講數(shù)據(jù)統(tǒng)計分析(參考版)
2025-05-17 14:25本頁面
【正文】 4 5 6],2,3) 結(jié)果為： M= To MATLAB（ liti6） 此命令產(chǎn)生了 2?3的正態(tài)分布隨機數(shù)矩陣，各數(shù)分別服從 N(1,), N(2,22), N(3, 32), N(4,), N(5, 22),N(6, 32) 返回 2021/6/14 33 給出數(shù)組 data的 頻數(shù)表 的命令為： [N,X]=hist(data,k) 此命令將區(qū)間 [min(data),max(data)]分為 k個小區(qū)間（缺省為 10），返回數(shù)組 data落在每一個小區(qū)間的頻數(shù) N和每一個小區(qū)間的中點 X. 描繪數(shù)組 data的 頻數(shù)直方圖 的命令為： hist(data,k) 四、頻 數(shù) 直 方 圖 的 描 繪 返回 2021/6/14 34 五、參數(shù)估計 正態(tài)總體的參數(shù)估計 設(shè)總體服從正態(tài)分布，則其點估計和區(qū)間估計可同時由以下命令獲得： [muhat,sigmahat,muci,sigmaci] = normfit(X,alpha) 此命令在顯著性水平 alpha下估計數(shù)據(jù) X的參數(shù)（ alpha缺省時設(shè)定為 ），返回值 muhat是 X的均值的點估計值，sigmahat是標準差的點估計值 , muci是均值的區(qū)間估計 ,sigmaci是標準差的區(qū)間估計 . 2021/6/14 35 其它分布的參數(shù)估計 有兩種處理辦法 : 一 .取容量充分大的樣本 （ n50） ， 按中心極限定理 ， 它近似地 服從正態(tài)分布； 二 .使用 Matlab工具箱中具有特定分布總體的估計命令 . （ 1） [muhat, muci] = expfit(X,alpha) 在顯著性水平 alpha下，求 指數(shù)分布的數(shù)據(jù) X的 均值的點估計及其區(qū)間估計 . （ 2） [lambdahat, lambdaci] = poissfit(X,alpha) 在顯著性水平 alpha下，求 泊松分布的數(shù)據(jù) X 的參數(shù)的點估計及其區(qū)間估計 . （ 3） [phat, pci] = weibfit(X,alpha) 在顯著性水平 alpha下，求 Weibull分布的數(shù)據(jù) X 的參數(shù)的點估計及其區(qū)間估計 . 返回 2021/6/14 36 六、假設(shè)檢驗 在總體服從正態(tài)分布的情況下，可用以下命令進行假設(shè)檢驗 . 總體方差 sigma2已知時，總體均值的檢驗使用 z檢驗 [h,sig,ci] = ztest(x,m,sigma,alpha,tail) 檢驗數(shù)據(jù) x 的關(guān)于均值的某一假設(shè)是否成立，其中 sigma 為已知方差， alpha 為顯著性水平，究竟檢驗什么假設(shè)取決于 tail 的取值： tail = 0，檢驗假設(shè)“ x 的均值等于 m ” tail = 1，檢驗假設(shè)“ x 的均值大于 m ” tail =1，檢驗假設(shè)“ x 的均值小于 m ” tail的缺省值為 0， alpha的缺省值為 . 返回值 h 為一個布爾值， h=1 表示可以拒絕假設(shè)， h=0 表示不可以拒絕假設(shè)， sig 為假設(shè)成立的概率， ci 為均值的 1alpha 置信區(qū)間 . 2021/6/14 37 例 7 Matlab統(tǒng)計工具箱中的數(shù)據(jù)文件 1993年一月份和二月份的汽油平均價格（ price1,price2分別是一，二月份的油價，單位為美分），它是容量為 20的雙樣本 .假設(shè)一月份油價的標準偏差是一加侖四分幣（ s=4），試檢驗一月份油價的均值是否等于 115. 解 作假設(shè)： m = 115. 首先取出數(shù)據(jù)，用以下命令： load gas 然后用以下命令檢驗 [h,sig,ci] = ztest(price1,115,4) 返回： h = 0， sig = ， ci = [ ]. 檢驗結(jié)果 : 1. 布爾變量 h=0, 表示不拒絕零假設(shè) . 說明提出的假設(shè)均值 115 是合理的 . 2. sig值為 , 遠超過 , 不能拒絕零假設(shè) 3. 95%的置信區(qū)間為 [, ], 它完全包括 115, 且精度很 高 . . To MATLAB（ liti7） 2021/6/14 38 總體方差 sigma2未知時，總體均值的檢驗使用 t檢驗 [h,sig,ci] = ttest(x,m,alpha,tail) 檢驗數(shù)據(jù) x 的關(guān)于均值的某一假設(shè)是否成立，其中alpha 為顯著性水平，究竟檢驗什么假設(shè)取決于 tail 的取值： tail = 0，檢驗假設(shè)“ x 的均值等于 m ” tail = 1，檢驗假設(shè)“ x 的均值大于 m ” tail =1，檢驗假設(shè)“ x 的均值小于 m ” tail的缺省值為 0， alpha的缺省值為 . 返回值 h 為一個布爾值， h=1 表示可以拒絕假設(shè)， h=0 表示不可以拒絕假設(shè)， sig 為假設(shè)成立的概率， ci 為均值的 1alpha 置信區(qū)間 . 2021/6/14 39 返回： h = 1， sig = ， ci =[ ]. 檢驗結(jié)果 : 1. 布爾變量 h=1, 表示拒絕零假設(shè) . 說明提出的假 設(shè)油價均值 115是不合理的 . 2. 95%的置信區(qū)間為 [ ], 它不包括 115, 故不能接受假設(shè) . 3. sig值為 , 遠小于 , 不能接受零 假設(shè) . To MATLAB（ liti8） 例 8 試檢驗例 8中二月份油價 Price2的均值是否等于 115. 解 作假設(shè)： m = 115， price2為二月份的油價，不知其方差，故用以下命令檢驗 [h,sig,ci] = ttest( price2 ,115) 2021/6/14 40 兩總體均值的假設(shè)檢驗 使用 t檢驗 [h,sig,ci] = ttest2(x,y,alpha,tail) 檢驗數(shù)據(jù) x ， y 的關(guān)于均值的某一假設(shè)是否成立 ， 其中alpha 為顯著性水平 ， 究竟檢驗什么假設(shè)取決于 tail 的取值： tail = 0， 檢驗假設(shè) “ x 的均值等于 y 的均值 ” tail = 1， 檢驗假設(shè) “ x 的均值大于 y 的均值 ” tail =1， 檢驗假設(shè) “ x 的均值小于 y 的均值 ” tail的缺省值為 0， alpha的缺省值為 . 返回值 h 為一個布爾值， h=1 表示可以拒絕假設(shè)， h=0 表示不可以拒絕假設(shè)， sig 為假設(shè)成立的概率， ci 為與 x與 y均值差的的 1alpha 置信區(qū)間 . 2021/6/14 41 返回： h = 1， sig = ， ci =[,]. 檢驗結(jié)果 :1. 布爾變量 h=1, 表示拒絕零假設(shè) . 說明提出的
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1