資源描述:
《SAS時間序列分析上機報告》由會員上傳分享,免費在線閱讀,更多相關內容在工程資料-天天文庫。
1、實驗序號:1240703020-1班級計算1001班姓名王傳雷學號20102346成績.基礎實驗模擬AR模型【實驗目的】直觀了解AR模型的平穩(wěn)性要求的含義,熟悉各種AR模型的樣本自相關函數(shù)和偏自相關函數(shù)的特點,為理論學習提供直觀印象?!緦嶒瀮热荨侩S機模擬各種AR模型實驗編號:B124070012班級姓名成績.【實驗要求】觀察各種AR序列的圖形及其樣本自相關函數(shù)和偏自相關函數(shù),體會AR序列的平穩(wěn)性含義及樣本自相關函數(shù)和樣本偏自相關函數(shù)的特點,總結如何判定時間序列的平穩(wěn)性。3、4.在editor窗中輸入
2、如下程序:模擬時間序列模型:dataa;x1=0.5;n=-50;doi=-50to1000;a=rannor(32565);*隨機正態(tài),32565是種子;x=a+0.4*x1;x1=x;n=n+1;ifi>0thenoutput;end;run;procprintdata=a;varx;procgplotdata=a;symboli=splinec=red;plotx*n;run;得到的結果如圖一所示:圖一時間序列5.在editor窗中輸入如下程序:(1)實驗序號:1240703020-1班級計算
3、1001班姓名王傳雷學號20102346成績.dataa;x1=0.5;n=-50;doi=-50to1000;a=rannor(32565);*隨機正態(tài),32565是種子;x=a+1.1*x1;x1=x;n=n+1;ifi>0thenoutput;end;run;procprintdata=a;varx;procgplotdata=a;symboli=splinec=red;plotx*n;run;得到的結果如下圖二所示:圖二時間序列dataa;x1=0.5;n=-50;doi=-50to1000
4、;a=rannor(32565);*隨機正態(tài),32565是種子;x=a-1.1*x1;x1=x;n=n+1;ifi>0thenoutput;end;run;procprintdata=a;varx;procgplotdata=a;symboli=splinec=red;plotx*n;run;(2)模擬時間序列模型:dataa;實驗序號:1240703020-1班級計算1001班姓名王傳雷學號20102346成績.x1=0.5;n=-50;doi=-50to1000;a=rannor(32565);
5、*隨機正態(tài),32565是種子;x=a-1.1*x1;x1=x;n=n+1;ifi>0thenoutput;end;run;procprintdata=a;varx;procgplotdata=a;symboli=splinec=red;plotx*n;run;得到的結果如下圖三所示:圖三時間序列6.從上面的圖形我們可以看到圖一的所有值均在0左右波動,整個時間序列是平穩(wěn)的,而圖二和圖三都是隨著時間的推遲,越發(fā)偏離初始值,是不平穩(wěn)的。根據(jù),只有,模型才是平穩(wěn)的,顯然后面兩個模型不符合條件。7.模擬自回歸
6、模型的程序如下:dataa;x1=0.5;x2=0.5;n=-50;doi=-50to1000;a=rannor(32565);*隨機正態(tài),32565是種子;x=a+0.5*x1+0.3*x2;x2=x1;x1=x;n=n+1;ifi>0thenoutput;end;實驗序號:1240703020-1班級計算1001班姓名王傳雷學號20102346成績.run;procprintdata=a;varx;procgplotdata=a;symboli=splinec=red;plotx*n;run;得
7、到的結果如下圖四所示:圖四模擬圖形8.針對第八題,我們利用用第七題的程序產生數(shù)據(jù),在根據(jù)第八題的相關程序觀察自相關函數(shù)和偏自相關函數(shù)的圖像。程序如下:dataa;x1=0.5;x2=0.5;n=-50;doi=-50to1000;a=rannor(32565);*???ú?yì?£?32565ê???×ó;x=a+0.5*x1+0.3*x2;x2=x1;x1=x;n=n+1;ifi>0thenoutput;end;run;procprintdata=a;varx;procgplotdata=a;sy
8、mboli=splinec=red;plotx*n;run;procarimadata=a;identifyvar=xnlag=10outcov=exp1;run;procgplotdata=exp1;symboli=needlewidth=6;plotcorr*lag;run;procgplotdata=exp1;實驗序號:1240703020-1班級計算1001班姓名王傳雷學號20102346成績.symboli=needlewidth=6;plotpartcor