河南工程学院课程设计《时间序列分析课程设计》学生姓名学号：学院：理学院专业班级：专业课程：时间序列分析课程设计指导教师：2017年6 月2 日目录1. 实验一澳大利亚常住人口变动分析 (1)1.1 实验目的 (1)1.2 实验原理 (1)1.3 实验内容 (2)1.4 实验过程 (3)2. 实验二我国铁路货运量分析 (8)2.1 实验目的 (8)2.2 实验原理 (8)2.3 实验内容 (9)2.4 实验过程 (10)3. 实验三美国月度事故死亡数据分析 (14)3.1 实验目的 (14)3.2 实验原理 (15)3.3 实验内容 (15)3.4 实验过程 (16)课程设计体会 (19)1.实验一澳大利亚常住人口变动分析1971年9月—1993年6月澳大利亚常住人口变动（单位：千人）情况如表1-1所示（行数据）。
表1-1（1）判断该序列的平稳性与纯随机性。
（2）选择适当模型拟合该序列的发展。
（3）绘制该序列拟合及未来5年预测序列图。
1.1 实验目的掌握用SAS软件对数据进行相关性分析，判断序列的平稳性与纯随机性,选择模型拟合序列发展。
1.2 实验原理（1）平稳性检验与纯随机性检验对序列的平稳性检验有两种方法，一种是根据时序图和自相关图显示的特征做出判断的图检验法；另一种是单位根检验法。
（2）模型识别先对模型进行定阶，选出相对最优的模型，下一步就是要估计模型中未知参数的值，以确定模型的口径，并对拟合好的模型进行显著性诊断。
（3）模型预测模型拟合好之后，利用该模型对序列进行短期预测。
1.3 实验内容（1）判断该序列的平稳性与纯随机性时序图检验，根据平稳时间序列均值、方差为常数的性质，平稳序列的时序图应该显示出该序列始终在一个常识值附近波动，而且波动的范围有界。
如果序列的时序图显示该序列有明显的趋势性或周期性，那么它通常不是平稳序列。
对自相关图进行检验时，可以用SAS 系统ARIMA 过程中的IDENTIFY 语句来做自相关图。
而单位根检验我们用到的是DF 检验。
以1阶自回归序列为例：11t t t x x φε-=+该序列的特征方程为：0λφ-=特征根为：λφ=当特征根在单位圆内时：11φ<该序列平稳。
当特征根在单位圆上或单位圆外时：11φ≥该序列非平稳。
对于纯随机性检验，既白噪声检验，可以用SAS 系统中的IDENTIFY 语句来输出白噪声检验的结果。
（2）选择适当模型拟合该序列的发展先对模型进行定阶，选出相对最优的模型，下一步就是要估计模型中未知参数的值，以确定模型的口径，并对拟合好的模型进行显著性诊断。
ARIMA过程的第一步是要IDENTIFY命令对该序列的平稳性和纯随机性进行识别，并对平稳非白噪序列估计拟合模型的阶数。
使用命令如下：proc print data=example3_20;IDENTIFY VAR =people nlag=8 minic p= (0:5) q =(0:5);run;（3）绘制该序列拟合及未来5年预测序列图模型拟合好之后，利用该模型对序列进行短期预测。
预测命令如下：forecast lead=5 id=time out=results;run;其中，lead指定预期数；id指定时间变量标识；out指定预测后期的结果存入某个数据集。
利用存储在临时数据集RESULTS里的数据，我们可以绘制拟合预测图，相关命令如下：proc gplot data=results;plot people*time=1 forecast*time=2 l95*time=3 u95*time=3/overlay;symbol1 c=red i=none v=star;symbol2 c=black i=join v=none;symbol3 c=green i=join v=none l=32;run;1.4 实验过程按照实验的过程运行程序，对程序结果的分析如下：（1）判断该序列的平稳性与纯随机性图1-1 1971年9月-1993年6月澳大利亚季度常住人口变动序列时序图时序图显示澳大利亚季度常住人口围绕在52千人附近随机波动，没有明显趋势或周期，基本可视为平稳模式。
图1-2序列自相关图自相关图显示该序列的自相关系数一直都比较小，始终控制在2倍的标准差范围以内，故认为该序列是平稳序列。
图1-3 序列的单位根检验结果根据第五列、第六列输出的结果我们可以判断，当显著性水平 取0.05时，序列非平稳，但当消除线性趋势之后序列平稳。
图1-4 白噪声检验输出结果可以看到延迟6阶、12阶的检验P值均小于0.05，故拒绝原假设，认为该序列为非白噪声序列（非纯随机序列）。
（2）选择适当模型拟合该序列的发展图1-5 IDENTIFY命令输出的最小信息量结果最后一条信息显示，在自相关延迟阶数也小于等于5的所有ARMA(p，q)模型中，BIC信息量相对于最小的是ARMA（1，3）模型。
图1-6 ESTIMATE命令输出的未知参数结果图1-7 ESTIMATE命令输出的拟合统计量结果图1-8 ESTIMATE 命令输出的系数矩阵图1-9 ESTIMATE 命令输出的残差自相关检验结果从输出结果可以看出由于延迟各阶的LB 统计量的P 值均显著大于α（0.05α≥），所以该拟合模型显著成立。
图1-10 ESTIMATE 命令输出的拟合模型形式该输出形式等价于：23(10.62415B 0.253693B 0.2953B )t t x ε=-++或记为：1230.624150.2536930.2953t t t t t x εεεε---=-++（3）绘制该序列拟合及未来5年预测序列图图1-11 FORECAST 命令输出的5年预测结果拟合效果图如图1-11：图1-12 拟合效果图2.实验二我国铁路货运量分析我国1949—2008年每年铁路货运量（单位：万吨）数据如表2-1所示。
表2-1请选择适当的模型拟合该序列，并预测2009—2013年我国铁路货运量。
2.1 实验目的掌握用SAS软件对数据进行相关性分析，掌握对非平稳时间序列的随机分析，选择合适模型，拟合序列发展。
2.2 实验原理ARIMA模型的预测和ARMA模型的预测方法非常类似。
(p,d,q)ARIMA模型的一般表示方法为：(B)(B)d t t x φε∇=Θ同时可以简记为：(B)(B)d t t x εΘ∇=Φ 式中，{}t ε 为零均值白噪声序列。
我们可以从上式看出，ARIMA 模型的实质就是差分与ARMA 模型的组合，这说明任何非平稳序列如果能通过适当阶数的差分实现差分后平稳，就可以对差分后序列进行ARMA 模型拟合。
（1）对差分平稳后的序列可以使用ARIMA 模型进行拟合，ARIMA 建模操作流程如图2-1所示。
图2-1 建模流程2.3 实验内容由于ARMA 模型是ARIMA 模型的一种特例，所以在SAS 系统中这两种模型的拟合都放在ARMA 过程中。
先利用时序图分析模型是否平稳，可以运用实验一的程序来实现。
再对该序平稳性检白噪声检分析结束 通过 差分运算拟合ARMA 模型未通过平稳不平稳获得观察值序列列进行1阶差分运算，同时考虑差分后序列的平稳性，添加如下命令：difhuoyunliang=dif(huoyunliang);命令“difhuoyunliang=dif(huoyunliang);”是指令系统对变量进行的1阶差分后的序列值赋值给变量difhuoyunliang，其中dif()是差分函数。
利用差分函数得出平稳模型。
再对模型进行定阶和进行预测。
模型定阶：identify var=difhuoyunliang(1) nlag=8 minic p=(0:5) q=(0:5)；模型预测：forecast lead=5 id=time；2.4 实验过程（1）判断序列的平稳性图2-2 我国1949—2008年每年铁路货运量时序图通过分析可知，该时序图有明显的上升趋势，所以为非平稳序列。
在此，对该序列进行1阶差分运算。
图2-3 1阶差分后序列时序图图2-4 1阶差分后序列自相关图通过分析可知，时序图显示差分后序列没有明显的非平稳特征；自相关图显示序列有很很强的短期相关性，所以可认为1阶差分后序列平稳。
对平稳的1阶查分序列进行白噪声检验，检验结果如图图2-5 1阶差分后序列白噪声检验默认显著性水平为0.05的条件下，由于延迟6阶、12阶的P值为0.0012和0.0098，小于0.05，所以该差分后序列不能视为白噪声序列，即差分后的序列还蕴含着不容忽视的相关信息可供提取。
（2）对平稳非白噪声查分序列进行拟合图2-6 IDENTIFY命令输出的最小信息量结果最后一条信息显示，在自相关延迟阶数也小于等于5的所有(p,q)ARMA模型中，BIC信息量相对于最小的是(1,0)ARMA模型。
考虑到前面已经进行的1阶差分运算，实际上是用(1,1,0)ARIMA模型拟合原序列。
图2-7 ESTIMATE命令输出的未知参数结果图2-8 ESTIMATE命令输出的拟合统计结果图2-8 ESTIMATE命令输出的残差自相关检验结果α≥），所以显然，拟合检验统计量的P值均显著大于显著性水平α（0.05可以认为改残差序列即为白噪声序列，显著性检验显示两参数均显著，这说明ARIMA模型对该序列建模成功。
(1,1,0)图2-10 ESTIMATE 命令输出的拟合模型形式输出结果显示，序列t x 的拟合模型为(1,1,0)ARIMA ,模型口径为：10.51983tt x Bε∇=-等价记为：121.519830.51983t t t t x x x ε--=-+利用拟合模型对序列做5期预测，结果如图2-10：图2-11 2009-2013我国铁路货运量预测3.实验三美国月度事故死亡数据分析据美国国家安全委员会统计，1973—1978年美国月度事故死亡数据如表3-1所示。
请选择适当模型拟合该序列的发展。
3.1 实验目的掌握用SAS软件对数据进行相关性分析，掌握对非平稳时间序列的随机分析，选择合适模型，拟合序列发展。
3.2 实验原理在SAS 系统中有一个AUTOREG 程序，可以进行残差自相关回归模型拟合。
残差自回归模型的构思是首先通过确定性因素分解方法提取序列中主要的确定性信息：t t t t x T S ε=++ （1）式中，t T 为趋势效应拟合；t S 为季节效应拟合。
考虑到因素分解方法对确定性信息的提取可能不够充分，因而需要进一步检验残差序列{}t ε的自相关性。
如果检验结果显示残差序列的自相关性不显著说明确定性回归模型（1）对信息的提取比较充分，可以停止分析。
如果检验结果显示残差序列的自相关显著，说明确定性回归模型（1）对信息的提取不充分，这时可以考虑对残差序列拟合自回归模型，进一步提取相关信息：11t t p t p t a εφεφε--=+⋅⋅⋅++这样构造的模型:t t t t x T S ε=++11t t p t p t a εφεφε--=+⋅⋅⋅++(a )0t E =,2(a )t Var σ=,(a ,a )0t t i Cov -=,1i ∀≥这就是自回归模型。