MATLAB參數估計與假設檢驗

1、 MATLAB參數估計與假設檢驗課型：新授課 教具：多媒體教學設備，matlab教學軟件一、目標與要求掌握matlab統(tǒng)計工具箱中的基本統(tǒng)計命令及其應用。二、教學重點與難點本堂課教學的重點在于引導學生在編寫matlab程序時能夠熟練運用基本統(tǒng)計量的相關命令實現相應的功能。三、教學方法 本課程主要通過講授法、演示法、練習法等相結合的方法來引導學生掌控本堂課的學習內容。四、教學內容上機內容回顧一、基本的統(tǒng)計量命令二、常見概率分布函數新授課統(tǒng)計推斷： 通過對樣本的處理和分析，得出與總參數相關的結論。統(tǒng)計推斷包括參數估計和假設檢驗兩部分內容。示例：吸煙對血壓有影響嗎？對吸煙和不吸煙兩組人群進行24小時

2、動態(tài)監(jiān)測，吸煙組66人，不吸煙組62人，分別測量24小時收縮壓（ 24hSBP）和舒張壓（ 24hDBP），白天（ 6Am-10Pm）收縮壓（ dSBP）和舒張壓（ dDBP ）， 夜間（ 10Pm-6Am）收縮壓（ nSBP）和舒張壓（ nDBP）。然后分別計算每類的樣本均值和標準差問題： 1）任何一個考察的時段，吸煙和不吸煙群體的血壓的真值分別是多少？（參數估計）2）吸煙和不吸煙群體的血壓的真值是否有區(qū)別？（假設檢驗）概念：第一部分：一：點估計1 矩估計法2 似然函數法二、評價估計優(yōu)劣的標準1 無偏性2 有效性3一致性三、區(qū)間估計參數估計的MATLAB實現：例題：50名17歲城市男性學生身

3、高（單位： cm）：170.1 179.0 171.5 173.1 174.1 177.2 170.3 176.2 163.7 175.4163.3 179.0 176.5 178.4 165.1 179.4 176.3 179.0 173.9 173.7173.2 172.3 169.3 172.8 176.4 163.7 177.0 165.9 166.6 167.4174.0 174.3 184.5 171.9 181.4 164.6 176.4 172.4 180.3 160.5166.2 173.5 171.7 167.9 168.7 175.6 179.6 171.6 168.1

4、172.2運行結果標準差區(qū)間估計 （4.4863，6.6926）標準差點估計 5.3707均值區(qū)間估計（171.1777, 174.2303）均值點估計 172.7040第二部分假設檢驗總體均值的假設檢驗 總體方差的假設檢驗 兩總體的假設檢驗 0-1分布總體均值的假設檢驗 總體分布正態(tài)性檢驗 假設檢驗的MATLAB實現假設檢驗MATLAB的實現MATLAB命令使用說明輸入參數x是樣本（n維數組），mu是H0中的µ0，sigma是總體標準差，alpha是顯著性水平（缺省時設定為0.05），tail是對雙側檢驗和兩個單側檢驗的標識，用備選假設H1確定：H1為µµ0時令

5、tail=0（可缺?。籋1為µ>µ0時令tail=1；H1為µ<µ0時令tail=-1。輸出參數h=0表示接受H0,h=1表示拒絕H0,sig標示對假設的接受和拒絕程度。ci給出的置信區(qū)間，zval是樣本統(tǒng)計量z的值。例題用N(5,1)隨機數產生n=100的樣本，在總體方差未知的情況下分別取=0.05和=0.01檢驗總體均值µ 5.2。H 0 : µ 5.2, H1 : µ < 5.2MATLAB命令x = normrnd(5,1,100,1);m = mean(x),h1,sig1,ci1 = ttest(x,5.2,0.05,-1)h2,sig2,ci2 = ttest(x,5.2,0.01,-1)H 0 : µ 5.2, H1 : µ < 5.2計算機結果m = 5.0111,h1,sig1,ci1 = 1 0.0343 -Inf 5.1815h2,sig2,ci2 = 0 0.0343 -Inf 5.2537可知在=0.05下拒絕H0（此時sig1<），