管理決策模型與方法---實驗報告一(共5頁)

1、 魯東大學(dxu)信息與電氣工程學院管理決策模型(mxng)與方法實驗報告（一） 實驗(shyn)名稱：MATLAB基本統(tǒng)計命令及其應用 專 業(yè)： 班 級： 姓 名： 學 號： 一、實驗目的：掌握matlab統(tǒng)計工具箱中的基本統(tǒng)計命令及其應用。二、實驗要求：將相關程序（M文件）和運行結果粘貼到word文檔中，說明操作過程中出現(xiàn)的問題及改正方法，以及后續(xù)學習應注意的問題。三、實驗(shyn)內容(nirng)： 1、某校60名學生(xu sheng)的一次考試成績如下: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)計算均值、標準差、極差、偏度、峰度，畫出直方圖； 2)檢驗分布的正態(tài)性； 3)若檢驗符合正態(tài)分布，估計正態(tài)分布的參數(shù)并檢驗參數(shù). 1）建立存儲數(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);%標準差range=max(a)-min(a);%極差，也可以用range(a)skewness=skewness(a);%偏度kurtosis=kurtosis(a);%峰度load aa hist(a, 10) %直方圖2) 分布的正態(tài)性檢驗：:load aanormplot(a)3)參數(shù)估

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

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）分別用兩個月的數(shù)據(jù)驗證這種說法的可靠性； 2）分別給出1月和2月汽油價格的置信區(qū)間； 3）給出1月和2月汽油價格差的置信區(qū)間.1)建立存儲數(shù)據(jù)文件,驗證這兩個月的數(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)驗證兩個月數(shù)據(jù)的均值是否都為115：load bbh,sig,ci = ttest(y1,115) 檢驗結果: 1. 布爾變量h=0, 表示不拒絕零假設. 說明提出的假設汽油價格均值115是合理的。2. 95%的

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

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

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

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

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