应用时间序列分析 第5章

《应用时间序列分析 第5章》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、佛山科学技术学院应用时间序列分析实验报告实验名称第五章非平稳序列的随机分析一、上机练习通过第4章我们学习了非平稳序列的确定性因素分解方法，但随着研究方法的深入和研究领域的拓宽，我们发现确定性因素分解方法不能很充分的提取确定性信息以及无法提供明确有效的方法判断各因素之间确切的作用关系。第5章所介绍的随机性分析方法弥补了确定性因素分解方法的不足，为我们提供了更加丰富、更加精确的时序分析工具。5.8.1拟合ARIMA模型【程序】dataexample5_1;inputx@@;difx=dif(x);t=_n_;cards;1.05-0.84-1.4

2、20.202.816.725.404.385.524.462.89-0.43-4.86-8.54-11.54-16.22-19.41-21.61-22.51-23.51-24.49-25.54-24.06-23.44-23.41-24.17-21.58-19.00-14.14-12.69-9.48-10.29-9.88-8.33-4.67-2.97-2.91-1.86-1.91-0.80;procgplot;plotx*tdifx*t;symbolv=starc=blacki=join;procarima;identifyvar=x(1);e

3、stimatep=1;estimatep=1noint;forecastlead=5id=tout=out;procgplotdata=out;plotx*t=1forecast*t=2l95*t=3u95*t=3/overlay;symbol1c=blacki=nonev=star;symbol2c=redi=joinv=none;symbol3c=greenI=joinv=none;12run;【结果及分析】1、输出序列x时序图：由图1-1显示出这是一个典型的非平稳序列，考虑对该序列进行1阶差分。图1-1序列x时序图2、序列difx时序图

4、：如图1-2所示，时序图显示差分后序列difx没有明显的非平稳特征。图1-2序列difx时序图3、序列difx白噪声检验：图1-3所示，由结果可知拒绝原假设，1阶差分后序列difx为平稳非白噪声序列。12图1-3序列difx白噪声检验结果4、模型定阶：由图1-4自相关图和图1-5偏自相关图我们可以看到，序列difx具有显著的自相关系数不截尾，偏自相关系数1阶截尾的性质，对1阶差分后序列difx拟合AR(1)模型。图1-4序列difx自相关图图1-5序列difx偏自相关图5、参数估计及模型检验：图1-6参数估计结果由图1-6我们看到，参数估计结

5、果显示常数项不显著，消除常数项再拟合AR(1)模型，结果如图1-7。图1-7消除常数项后的参数估计结果12由图1-7结果知，参数显著。再看序列difx的残差自相关检验结果，如图1-8所示：图1-8残差自相关检验结果图1-8结果显示，延迟各阶的LB统计量的P值均显著大于，所以该拟合模型显著有效，拟合结果如图1-9所示：图1-9模型拟合结果6、模型口径：根据前面1-5的分析以及输出结果，我们可以得到序列x的拟合模型为，模型记为：。7、短期预测：利用拟合模型对序列x进行5期预测，预测结果和拟合效果图如下。图1-10预测结果图1-11拟合效果图125

6、.8.2拟合Auto-Regressive模型【程序】dataexample5_2;inputx@@;lagx=lag(x);t=_n_;cards;3.038.4610.229.8011.962.838.4313.7716.1816.8419.5713.2614.7824.4828.1628.2732.6218.4425.2538.3643.7044.4650.6633.0139.9760.1768.1268.8478.1549.8462.2391.49103.20104.53118.1877.8894.75138.36155.68157.

7、46177.69117.15;procgplotdata=example5_2;plotx*t=1;symbol1c=blacki=joinv=star;run;procautoregdata=example5_2;modelx=t/dwprob;procautoregdata=example5_2;modelx=t/nlag=5backstepmethod=ml;outputout=outp=xppm=trend;procautoregdata=example5_2;modelx=t/nlag=5backstepmethod=mlnoint

8、;outputout=outp=xppm=trend;procgplotdata=out;plotx*t=2xp*t=3trend*t=4/overlay;symb