MATLAB參數(shù)估計(jì)與假設(shè)檢驗(yàn)

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

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

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運(yùn)行結(jié)果標(biāo)準(zhǔn)差區(qū)間估計(jì) （4.4863，6.6926）標(biāo)準(zhǔn)差點(diǎn)估計(jì) 5.3707均值區(qū)間估計(jì)（171.1777, 174.2303）均值點(diǎn)估計(jì) 172.7040第二部分假設(shè)檢驗(yàn)總體均值的假設(shè)檢驗(yàn) 總體方差的假設(shè)檢驗(yàn) 兩總體的假設(shè)檢驗(yàn) 0-1分布總體均值的假設(shè)檢驗(yàn) 總體分布正態(tài)性檢驗(yàn) 假設(shè)檢驗(yàn)的MATLAB實(shí)現(xiàn)假設(shè)檢驗(yàn)MATLAB的實(shí)現(xiàn)MATLAB命令使用說(shuō)明輸入?yún)?shù)x是樣本（n維數(shù)組），mu是H0中的µ0，sigma是總體標(biāo)準(zhǔn)差，alpha是顯著性水平（缺省時(shí)設(shè)定為0.05），tail是對(duì)雙側(cè)檢驗(yàn)和兩個(gè)單側(cè)檢驗(yàn)的標(biāo)識(shí)，用備選假設(shè)H1確定：H1為µµ0時(shí)令

5、tail=0（可缺省）；H1為µ>µ0時(shí)令tail=1；H1為µ<µ0時(shí)令tail=-1。輸出參數(shù)h=0表示接受H0,h=1表示拒絕H0,sig標(biāo)示對(duì)假設(shè)的接受和拒絕程度。ci給出的置信區(qū)間，zval是樣本統(tǒng)計(jì)量z的值。例題用N(5,1)隨機(jī)數(shù)產(chǎn)生n=100的樣本，在總體方差未知的情況下分別取=0.05和=0.01檢驗(yàn)總體均值µ 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計(jì)算機(jī)結(jié)果m = 5.0111,h1,sig1,ci1 = 1 0.0343 -Inf 5.1815h2,sig2,ci2 = 0 0.0343 -Inf 5.2537可知在=0.05下拒絕H0（此時(shí)sig1<），