管理決策模型與方法---實(shí)驗(yàn)報(bào)告一(共5頁(yè))

1、 魯東大學(xué)(dxu)信息與電氣工程學(xué)院管理決策模型(mxng)與方法實(shí)驗(yàn)報(bào)告（一） 實(shí)驗(yàn)(shyn)名稱：MATLAB基本統(tǒng)計(jì)命令及其應(yīng)用 專 業(yè)： 班 級(jí)： 姓 名： 學(xué) 號(hào)： 一、實(shí)驗(yàn)?zāi)康模赫莆誱atlab統(tǒng)計(jì)工具箱中的基本統(tǒng)計(jì)命令及其應(yīng)用。二、實(shí)驗(yàn)要求：將相關(guān)程序（M文件）和運(yùn)行結(jié)果粘貼到word文檔中，說(shuō)明操作過(guò)程中出現(xiàn)的問(wèn)題及改正方法，以及后續(xù)學(xué)習(xí)應(yīng)注意的問(wèn)題。三、實(shí)驗(yàn)(shyn)內(nèi)容(nirng)： 1、某校60名學(xué)生(xu sheng)的一次考試成績(jī)?nèi)缦?93 75 83 93 91 85 84 82 77 76 77 95 94 89 91 88 86 83 96 81 79

2、97 78 75 67 69 68 84 83 81 75 66 85 70 94 84 83 82 80 78 74 73 76 70 86 76 90 89 71 66 86 73 80 94 79 78 77 63 53 551)計(jì)算均值、標(biāo)準(zhǔn)差、極差、偏度、峰度，畫出直方圖； 2)檢驗(yàn)分布的正態(tài)性； 3)若檢驗(yàn)符合正態(tài)分布，估計(jì)正態(tài)分布的參數(shù)并檢驗(yàn)參數(shù). 1）建立存儲(chǔ)數(shù)據(jù)文件aa：a1=93 75 83 93 91 85 84 82 77 76;a2=77 95 94 89 91 88 86 83 96 81;a3=79 97 78 75 67 69 68 84 83 81;a4=75

3、 66 85 70 94 84 83 82 80 78;a5=74 73 76 70 86 76 90 89 71 66;a6=86 73 80 94 79 78 77 63 53 55;a=a1 a2 a3 a4 a5 a6;save aa a; mean=mean(a);%均值std=std(a);%標(biāo)準(zhǔn)差range=max(a)-min(a);%極差，也可以用range(a)skewness=skewness(a);%偏度kurtosis=kurtosis(a);%峰度load aa hist(a, 10) %直方圖2) 分布的正態(tài)性檢驗(yàn)：:load aanormplot(a)3)參數(shù)估

4、計(jì)：load aamuhat,sigmahat,muci,sigmaci = normfit(a)假設(shè)檢驗(yàn):load aah,sig,ci = ttest(a,80.1) %依據(jù)參數(shù)估計(jì)得到的均值的點(diǎn)估計(jì)值為80.1檢驗(yàn)結(jié)果: 1. 布爾變量(binling)h=0, 表示(biosh)不拒絕零假設(shè). 說(shuō)明提出(t ch)的假設(shè)成績(jī)均值80.1是合理的.2. 95%的置信區(qū)間為77.5915 82.6085, 它完全包括80.1, 且精度很高.3. sig值為1, 遠(yuǎn)超過(guò)0.5, 不能拒絕零假。2、據(jù)說(shuō)某地汽油的價(jià)格是每加侖115美分，為了驗(yàn)證這種說(shuō)法，一位學(xué)者開(kāi)車隨機(jī)選擇了一些加油站，得到某

5、年1月和2月的數(shù)據(jù)如下： 1月：119 117 115 116 112 121 115 122 116 118 109 112 119 112 117 113 114 109 109 1182月：118 119 115 122 118 121 120 122 128 116 120 123 121 119 117 119 128 126 118 1251）分別用兩個(gè)月的數(shù)據(jù)驗(yàn)證這種說(shuō)法的可靠性； 2）分別給出1月和2月汽油價(jià)格的置信區(qū)間； 3）給出1月和2月汽油價(jià)格差的置信區(qū)間.1)建立存儲(chǔ)數(shù)據(jù)文件,驗(yàn)證這兩個(gè)月的數(shù)據(jù)都服從正態(tài)分布：y1=119 117 115 116 112 121 115

6、 122 116 118 109 112 119 112 117 113 114 109 109 118;y2=118 119 115 122 118 121 120 122 128 116 120 123 121 119 117 119 128 126 118 125;save bb y1 y2 hist(y1,7）hist(y2,7)normplot(y1)normplot(y2)驗(yàn)證兩個(gè)月數(shù)據(jù)的均值是否都為115：load bbh,sig,ci = ttest(y1,115) 檢驗(yàn)結(jié)果: 1. 布爾變量h=0, 表示不拒絕零假設(shè). 說(shuō)明提出的假設(shè)汽油價(jià)格均值115是合理的。2. 95%的

7、置信區(qū)間為113.3388 116.9612, 它完全(wnqun)包括115, 且精度(jn d)很高。3. sig值為0.8642, 遠(yuǎn)超過(guò)0.5, 不能拒絕(jju)零假。load bbh,sig,ci = ttest(y2,115) 檢驗(yàn)結(jié)果: 1. 布爾變量h=1, 表示拒絕零假設(shè). 說(shuō)明提出的假設(shè)汽油價(jià)格均值115是不合理的。2. 95%的置信區(qū)間為119.0129 122.4871, 它不包括115, 故不能接受假設(shè)。3. sig值為1.3241e-006, 遠(yuǎn)小于0.5, 不能接受零假。2）1月份汽油價(jià)格參數(shù)估計(jì)：load bbmuhat,sigmahat,muci,sigma

8、ci = normfit(y1)均值的點(diǎn)估計(jì)值muhat =115.1500標(biāo)準(zhǔn)差的點(diǎn)估計(jì)值sigmahat =3.8699均值的區(qū)間估計(jì)muci = 113.3388 116.9612標(biāo)準(zhǔn)差的區(qū)間估計(jì)sigmaci = 2.9430 5.6523假設(shè)檢驗(yàn)h,sig,ci = ttest(y1,115.15) 1. 布爾變量h=0, 表示不拒絕零假設(shè). 說(shuō)明提出的假設(shè)汽油價(jià)格均值115.15是合理的。2. 95%的置信區(qū)間為113.3388 116.9612, 它完全包括115, 且精度很高。3. sig值為1, 遠(yuǎn)超過(guò)0.5, 不能拒絕零假。2月份汽油價(jià)格參數(shù)估計(jì)：load bbmuhat,

9、sigmahat,muci,sigmaci = normfit(y2)均值的點(diǎn)估計(jì)值muhat = 120.7500標(biāo)準(zhǔn)差的點(diǎn)估計(jì)值sigmahat = 3.7116均值的區(qū)間估計(jì)muci = 119.0129 122.4871標(biāo)準(zhǔn)差的區(qū)間估計(jì)sigmaci = 2.8227 5.4211假設(shè)檢驗(yàn)h,sig,ci = ttest(y2,120.75) 1. 布爾變量h=0, 表示不拒絕零假設(shè). 說(shuō)明提出的假設(shè)汽油價(jià)格均值120.75是合理的。2. 95%的置信區(qū)間為119.0129 122.4871, 它完全包括115, 且精度很高。3. sig值為1, 遠(yuǎn)超過(guò)0.5, 不能拒絕零假。3）lo

10、ad bbh,sig,ci = ttest2(y1,y2)檢驗(yàn)結(jié)果:1. 布爾變量h=1, 表示拒絕零假設(shè). 說(shuō)明提出的假設(shè)“油價(jià)均值相同”是不合理的。2. 95%的置信區(qū)間為-8.0273 -3.1727,說(shuō)明(shumng)一月份油價(jià)比二月份油價(jià)約低3.2至8.1美分。3. sig值為3.6952e-005, 遠(yuǎn)小于0.5, 不能接受(jishu)“油價(jià)均相同(xin tn)”假設(shè). 3、從某電工器材廠生產(chǎn)的一批保險(xiǎn)絲中抽取10根，測(cè)試其融化時(shí)間，得到數(shù)據(jù)如下：42 65 75 78 71 59 57 68 55 54設(shè)這批保險(xiǎn)絲的融化時(shí)間服從正態(tài)分布，檢驗(yàn)總體方差是否等于122？x=42

11、 65 75 78 71 59 57 68 55 54;var0=12.2;h,p,ci=vartest(x,var0)運(yùn)行結(jié)果h = 0p = 0.8530ci = 57.6362 406.0155說(shuō)明：1、布爾變量h=0說(shuō)明在顯著水平為0.05（缺省狀態(tài)下為0.05）下，不拒絕零假設(shè)（方差為144），說(shuō)明提出的假設(shè)是合理的；2. 95%的置信區(qū)間為57.6362 406.0155, 它完全包括144,說(shuō)明假設(shè)合理。3. p值為0.8530, 假設(shè)成立的概率大于0.5,接受總體方差等于122的假設(shè).對(duì)于版本較低的matlab，統(tǒng)計(jì)包中可能不包含vartest函數(shù)包，可使用如下語(yǔ)句編寫（依據(jù)概率論課本中的正態(tài)總體方差的假設(shè)檢驗(yàn)）alpha=0.05;%顯著水平var0=12.2;%已知的方差x=42 65 75 78 71 59 57 68 55 54;n=length(x);%樣本容量var=var(x);%樣本方差chi2=(n-1)*var/var0 %卡方分布統(tǒng)計(jì)量的觀測(cè)值lambda1=chi2inv(alpha/2,n-1);%給定水平下卡方統(tǒng)計(jì)的臨界值lambda2=chi2inv(1-alpha/2,n-1);la