时间序列分析第三章平稳时间序列分析1816.pdfVIP

创作时间：二零二一年六月三十日 应 用 时 间 序 列 分 析 实 验 报 告之巴公井开创作 创作时间：二零二一年六月三十日 实验名称 第三章 平稳时间序列分析 一、上机练习 data example3_1; input x@@; time=_n_; cards; ; procgplotdata=example3_1; plot x*time=1; symbolc=red i=join v=star; run; 建立该数据集, 绘制该序列时序图得： 根据所得图像, 对序列进行平稳性检验.时序图就是一个平面二维坐标图, 通常横轴暗示时间, 纵轴暗 示序列取值.时序图可以直观地帮手我们掌握时间序列的一些基天职布特征. 根据平稳时间序列均值、方差为常数的性质 , 平稳序列的时序图应该显示出该序列始终在一个常数值附 近随机摆荡, 而且摆荡的范围有界的特点.如果观察序列的时序图, 显示出该序列有明显的趋势性或周期 性, 那它通常不是平稳序列.从图上可以看出, 数值围绕在 0 附近随机摆荡, 没有明显或周期, 其本可以视 为平稳序列, 时序图显示该序列摆荡平稳. procarimadata=example3_1; identifyvar=x nlag=8; run; 创作时间：二零二一年六月三十日 创作时间：二零二一年六月三十日 图一 图二样本自相关图 图三 样本逆自相关图 图四 样本偏自相关图 图五 纯随机检验图 实验结果分析： （1）由图一我们可以知道序列样本的序列均值为-0.06595, 标准差为1.561613, 观察值个数为84 个. （2 ）根据图二序列样本的自相关图我们可以知道该图横轴暗示自相关系数, 综轴暗示延迟时期数, 用水平方向的垂线暗示自相关系数的年夜小.我们发现样本自相关图延迟 3 阶之后, 自相关系数都落入2 倍标准差范围以内, 而且自相关系数向0.03 衰减的速度非常快, 延迟5 阶之后自相关系数即在0.03 值附 近摆荡.这是一个短时间相关的样本自相关图.所以根据样本自相关图的相关性质, 可以认为该序列平稳. （3 ）根据图五的检验结果我们知道, 在各阶延迟下LB 检验统计量的P 值都非常小（0.0001 ）, 所 创作时间：二零二一年六月三十日 创作时间：二零二一年六月三十日 以我们可以以很年夜的掌控（置信水平99.999% ）判定该序列样本属于非白噪声序列. procarimadata=example3_1; identifyvar=x nlag=8minicp= (0:5) q=(0:5); run; IDENTIFY 命令输出的最小信息量结果 某个观察值序列通过序列预处置, 可以判定为平稳非白噪声序列, 就可以利用ARMA 模型对该序列 建模.建模的基本步伐如下： A ：求出该观察值序列的样本自相关系数（ACF ）和样本偏自相关系数（PACF ）的值. B ：根据样本自相关系数和偏自相关系数的性质, 选择适本地ARMA （p, q ）模型进行拟合. C ：估计模型中未知参数的值. D ：检验模型有效性.如果拟合模型欠亨过检验, 转向步伐B, 重新选择模型再拟合. E ：模型优化.如果拟合模型通过检验, 仍然转向步伐B, 充沛考虑各种可能, 建立多个拟合模型, 从所 有通过检验中选择最优模型. F ：利用拟合模型, 预测序列的将来走势. 为了尽量防止因个人经验缺乏招致的模型识别问题, SAS 系统还提供了相对最优模型识别.最后一条 信息显示, 在自相关延迟阶数小于即是5, 移动平均延迟阶数也小于即是5 的所有ARMR （p,q ）模型中, BIC 信息量相对最小的是ARMR （0,4 ）模型, 即MA （4 ）模型. 需要注意的是, MINIC 只给出一定