2014統(tǒng)計計算課程設計.doc

統(tǒng)計計算課程設計報告學院 專業(yè) 姓名 學號 評語：分數(shù) 二一四年五月2014統(tǒng)計計算課程設計題型題型一：若產(chǎn)生總體，其中未知，請設計一個隨機模擬實驗，要求從該總體產(chǎn)生一個容量為100的樣本，考慮置信水平分別取0.95和0.5時，對上述過程重復1000次，統(tǒng)計有多少個區(qū)間包含均值5，要求畫出置信水平分別取0.95和0.5時均值的置信區(qū)間圖，并給出實驗總結(jié)。問題分析： 實驗要求在置信水平分別取0.95和0.5的情況下，從總體中產(chǎn)生一個容量為100的樣本，由于未知，在分析中先暫且把它設置為。由于方差未知，估計正態(tài)總體均值的置信區(qū)間時使用公式。如果均值5在置信區(qū)間內(nèi)，那么符合條件的區(qū)間數(shù)加1。該過程重復1000次，統(tǒng)計最終符合條件的區(qū)間的頻數(shù)為多少，對應的頻率為多少。SAS結(jié)果：圖1 輸出結(jié)果：符合條件的區(qū)間數(shù)累加結(jié)果以及頻率圖2 置信水平取0.95時均值的置信區(qū)間圖圖3 置信水平分別取0.5時均值的置信區(qū)間圖 由于對SAS作圖操作的了解程度有限，尚未能掌握畫出標準的置信區(qū)間圖的方法。圖2、圖3中，橫軸表示置信區(qū)間的編號，縱軸表示總體均數(shù)；星號表示的是該編號的置信區(qū)間上限，點表示的是該編號的置信區(qū)間下限；中間是總體均值等于5的參考線，方便觀察對比得出結(jié)論。結(jié)論：在置信水平取時，如果從同一總體中重復抽取1000份樣本含量相同（本實驗樣本容量為100）的獨立樣本，每份樣本分別計算1個置信區(qū)間，在這1000個置信區(qū)間中將大約有個置信區(qū)間覆蓋總體均數(shù)，大約有個置信區(qū)間并不覆蓋總體均數(shù)。所以，對于某一次估計的置信區(qū)間，我們平時總是宣稱這個區(qū)間覆蓋了總體均數(shù)，但不一定是真的覆蓋了總體均數(shù)，于是，我們補充一句：置信度為。題型二：在實際觀察中，已知腐蝕深度與腐蝕時間有線性關(guān)系，設給定腐蝕時間X時腐蝕深度Y的總體均數(shù)E(Y|X)與X的關(guān)系滿足方程E(Y|X)=70+0.6X，且腐蝕時間，腐蝕深度?，F(xiàn)隨機抽取該總體20對腐蝕深度與腐蝕時間的關(guān)系，構(gòu)成一份樣本，做一次回歸分析；重復抽取相同樣本量的10份樣本，分別進行回歸，得到10條直線，觀察它們的圖形，得出結(jié)論。要求：（1）給出隨機樣本表；（2）10條回歸重疊圖形；（3）實驗結(jié)論。問題分析： 實驗要求從總體，中隨機抽取20對和的關(guān)系。然后根據(jù)這20對樣本做一次回歸分析。該過程重復10次，并畫出這10條回歸直線，觀察并得出相應結(jié)論。SAS結(jié)果：圖4 隨機樣本表圖5 10條回歸直線重疊圖形結(jié)論：觀察圖5可以發(fā)現(xiàn)，10條回歸直線的趨勢大致相同，但是具體每條直線的截距和斜率都存在著差異。同時可以比較10個模型的回歸結(jié)果和樣本的來源(截距為70，斜率為0.6)，相差也很大而且不穩(wěn)定。綜上所述，這10個回歸模型的擬合效果并不理想，造成這一現(xiàn)象的主要原因是樣本量不夠大。在一元線性回歸中，有。顯然越大，越小。所以，要想使的估計值更穩(wěn)定，在收集數(shù)據(jù)時，樣本量應盡可能大一些，樣本量大小時，估計量的穩(wěn)定性肯定不會太好。題型三：設有一個由兩個服務臺串聯(lián)組成的服務機構(gòu)（雙服務太串聯(lián)排隊系統(tǒng)）。顧客在第一個服務臺接受服務后進入第二個服務臺，服務完畢后離開。假定顧客達到第一個服務臺的時間間隔是均值為1分鐘的指數(shù)分布，顧客在第一個和第二個服務臺的服務時間分別是均值為0.7分和0.9分的指數(shù)分布。請模擬這種雙服務臺串聯(lián)排隊系統(tǒng)（分別模擬600分和1000分的系統(tǒng)）；并估計出顧客在兩個服務臺的平均逗留時間和排隊中的顧客平均數(shù)。問題分析：首先引入幾個記號：顧客到達第一個服務臺的時刻顧客到達第二個服務臺的時刻顧客在第一個服務臺的服務時間顧客在第二個服務臺的服務時間顧客在第一個服務臺的等待時間顧客在第二個服務臺的等待時間在第一個服務臺排隊的顧客數(shù)在第二個服務臺排隊的顧客數(shù)顧客離開第一個服務臺的時刻顧客離開第二個服務臺的時刻模擬時鐘從分開始，產(chǎn)生指數(shù)分布隨機數(shù)，比如得；在第一個服務臺的服務時間，產(chǎn)生隨機數(shù)比如得；在第二個服務臺的服務時間，產(chǎn)生隨機數(shù)比如得。分時，第一個顧客到達第一個服務臺，記為，因沒有人排隊，馬上接受服務，即，此時；第一個顧客在第一個服務臺接受服務時間為1分，計算分；接著進入第二個服務臺，記；因沒有人排隊，馬上接受服務，即，此時；第一個顧客在第二個服務臺接受服務時間為0.2分，計算分，即第一個顧客于開門后1.5分離開（即分時離開）。分時，第六個顧客到達第一個服務臺，記為，而根據(jù)前面的計算， ，即；此時在第一個服務臺的排隊中，第四個和第五個顧客仍在（因為，都大于，即第六個顧客到達時他們都還沒走），所以；第六個顧客在第一個服務臺接受服務時間為0.3分，計算分；接著進入第二個服務臺，記；而根據(jù)前面的計算，即；此時在第二個服務臺的排隊中，第四個和第五個顧客仍在（因為，都大于，即第六個顧客到達時他們都還沒走），所以；第六個顧客在第二個服務臺接受服務時間為2.0分，計算分，即第六個顧客于開門后6.2分離開（即分時離開）。一直按這個過程循環(huán)直至模擬時鐘的時間到達600或者1000分。表1、表2列出模擬600分系統(tǒng)試驗的部分結(jié)果。表1 模擬過程（輸入過程）顧客序號輸入過程到達間隔在第一個服務臺的服務時間在第二個服務臺的服務時間10.31.00.220.90.60.430.30.30.340.40.31.550.10.10.160.40.32.0表2 模擬過程（輸出結(jié)果）顧客序號模擬試驗過程的輸出結(jié)果第一個服務臺第二個服務臺到達時刻服務時間等待時間離開時刻排隊的顧客數(shù)到達時刻服務時間等待時間離開時刻排隊的顧客數(shù)10.31.00.01.301.30.20.01.5021.20.60.11.911.90.40.02.3031.50.30.42.212.20.30.12.6141.90.30.32.512.51.50.14.1152.00.10.52.622.60.11.54.2162.40.30.22.922.92.01.36.22思路框圖： 下面以模擬600分系統(tǒng)為例，畫出流程圖。否是是否是否是置初始狀態(tài)產(chǎn)生eie(1),s1ie(1/0.7), s2ie(1/0.9) x1i=x1i+eix1ic1i? d1i=0d1i=c1i-x1iENDa1=a1+1st1i=d1i+s1ic1i=x1i+st1iy(a1)=c1i在第二個服務臺的思路流程是跟在第一個服務臺一樣的，這里由于篇幅原因就沒有具體給出。T=c2i T600?ENDk=1X1iy(k)?n1i=n1i+1k+1k=ai?sti=st1i+st2iSAS結(jié)果：圖6 顧客在兩個服務臺的平均逗留時間和排隊中的顧客平均數(shù)的估計結(jié)果結(jié)論：在模擬600分的系統(tǒng)中，顧客在第一個服務臺的平均逗留時間為2.59分，排隊中的顧客平均數(shù)為2人；在第二個服務臺的平均逗留時間為6.43分，排隊中的顧客平均數(shù)為7人；平均每個顧客在服務機構(gòu)逗留的總時間為9.02分。在模擬1000分的系統(tǒng)中，顧客在第一個服務臺的平均逗留時間為2.15分，排隊中的顧客平均數(shù)為2人；在第二個服務臺的平均逗留時間為6.05分，排隊中的顧客平均數(shù)為6人；平均每個顧客在服務機構(gòu)逗留的總時間為8.20分。程序：題型一：data ex1;array x(100) x1-x100; /*定義數(shù)列的元素*/do alpha=0.95,0.5; /*置信水平取0.95、0.5時各執(zhí)行程序一次*/t=tinv(1-(1-alpha)/2,99); /*計算*/sigma=1; /*設置方差*/m=0; /*初始化符合條件的區(qū)間數(shù)*/do j=1 to 1000;/*重復1000次實驗*/do i=1 to 100;/*每次實驗產(chǎn)生100個樣本*/r=rannor(32789);x(i)=5+r*1; /*產(chǎn)生服從，正態(tài)分布的隨機數(shù)*/if i=100 then do; /*當輸出第100個樣本時*/mean=mean(of x1-x100); /*計算均值*/s=std(of x1-x100); /*計算標準差*/delta=t*s/sqrt(100); /*計算*/lcl=mean-delta; /*計算置信區(qū)間下限*/ucl=mean+delta; /*計算置信區(qū)間上限/if lcl5&5=time); /*工作時間小于設定的time*/n1=0;n2=0; /*各服務臺排隊顧客初始化*/e=round(-log(ranuni(seed),0.1); /*到達時間間隔*/s1=round(-0.7*log(ranuni(seed),0.1); /*顧客在第一服務臺的服務時間*/s2=round(-0.9*log(ranuni(seed),0.1); /*顧客在第二服務臺的服務時間*/x1=x1+e; /*顧客到達第一個服務臺時刻*/if x1c1 then d1=0; /*如果顧客到達第一個服務臺時刻大于上一名顧客離開第一個服務臺時刻，就不需要等待*/else do; /*否則，在第一個服務臺等待時間=上一名顧客離開時刻-這名顧客的到達時刻*/d1=c1-x1;do k=1 to a1;if x1c2) then d2=0; /*如果顧客到達第二個服務臺時刻大于上一名顧客離開第二個服務臺時刻，就不需要等待*/else do; /*否則，在第二個服務臺等待時間=上一名顧客離開時刻-這名顧客的到達時刻*/d2=c2-x2;do k=1 to a2;if x2z(k) then n2=n2+1; /*在第二個服務臺前排隊的顧客數(shù)*/end;end;a2=a2+1; /*累計到達第二個服務臺的顧客數(shù)*/st2=d2+s2; /*顧客在第二服務臺的逗留時間*/c2=x2+st2; /*顧客離開第二服務臺的時刻*/z(a2)=c2; /*將離開時間賦值到數(shù)列z*/st=st1+st2; /*顧客在