河南工程学院课程设计《时间序列分析课程设计》学生姓名学号：学院：理学院专业班级：专业课程：时间序列分析课程设计指导教师：2017年 6 月 2 日目录1. 实验一澳大利亚常住人口变动分析..... 错误!未定义书签。
实验目的............................................... 错误!未定义书签。
实验原理............................................... 错误!未定义书签。
实验内容............................................... 错误!未定义书签。
实验过程............................................... 错误!未定义书签。
2. 实验二我国铁路货运量分析........... 错误!未定义书签。
实验目的............................................... 错误!未定义书签。
实验原理............................................... 错误!未定义书签。
实验内容............................................... 错误!未定义书签。
实验过程............................................... 错误!未定义书签。
3. 实验三美国月度事故死亡数据分析...... 错误!未定义书签。
实验目的............................................... 错误!未定义书签。
实验原理............................................... 错误!未定义书签。
实验内容............................................... 错误!未定义书签。
实验过程............................................... 错误!未定义书签。
课程设计体会 ............................ 错误!未定义书签。
1.实验一澳大利亚常住人口变动分析1971年9月—1993年6月澳大利亚常住人口变动（单位：千人）情况如表1-1所示（行数据）。
表1-1（1）判断该序列的平稳性与纯随机性。
（2）选择适当模型拟合该序列的发展。
（3）绘制该序列拟合及未来5年预测序列图。
实验目的掌握用SAS软件对数据进行相关性分析，判断序列的平稳性与纯随机性,选择模型拟合序列发展。
实验原理（1）平稳性检验与纯随机性检验对序列的平稳性检验有两种方法，一种是根据时序图和自相关图显示的特征做出判断的图检验法；另一种是单位根检验法。
（2）模型识别先对模型进行定阶，选出相对最优的模型，下一步就是要估计模型中未知参数的值，以确定模型的口径，并对拟合好的模型进行显著性诊断。
（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）判断该序列的平稳性与纯随机性图1-1 1971年9月-1993年6月澳大利亚季度常住人口变动序列时序图时序图显示澳大利亚季度常住人口围绕在52千人附近随机波动，没有明显趋势或周期，基本可视为平稳模式。
图1-2序列自相关图自相关图显示该序列的自相关系数一直都比较小，始终控制在2倍的标准差范围以内，故认为该序列是平稳序列。
图1-3 序列的单位根检验结果根据第五列、第六列输出的结果我们可以判断，当显著性水平 取时，序列非平稳，但当消除线性趋势之后序列平稳。
图1-4 白噪声检验输出结果可以看到延迟6阶、12阶的检验P值均小于，故拒绝原假设，认为该序列为非白噪声序列（非纯随机序列）。
（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年我国铁路货运量。
实验目的掌握用SAS软件对数据进行相关性分析，掌握对非平稳时间序列的随机分析，选择合适模型，拟合序列发展。
实验原理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 建模流程实验内容由于ARMA 模型是ARIMA 模型的一种特例，所以在SAS 系统中这两种模型的拟合都放在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；实验过程（1）判断序列的平稳性图2-2 我国1949—2008年每年铁路货运量时序图通过分析可知，该时序图有明显的上升趋势，所以为非平稳序列。
在此，对该序列进行1阶差分运算。
difhuoyunliang-30000-20000-10000100002000030000timeJAN1945JAN1950JAN1955JAN1960JAN1965JAN1970JAN1975JAN1980JAN1985JAN1990JAN1995JAN2000JAN2005JAN2010图2-3 1阶差分后序列时序图图2-4 1阶差分后序列自相关图通过分析可知，时序图显示差分后序列没有明显的非平稳特征；自相关图显示序列有很很强的短期相关性，所以可认为1阶差分后序列平稳。
对平稳的1阶查分序列进行白噪声检验，检验结果如图图2-5 1阶差分后序列白噪声检验默认显著性水平为的条件下，由于延迟6阶、12阶的P值为和，小于，所以该差分后序列不能视为白噪声序列，即差分后的序列还蕴含着不容忽视的相关信息可供提取。
（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模型对该序列建模成功。