基于ARIMA模型对我国能源需求的预测(2)
时间:2026-01-25
时间:2026-01-25
基于ARIMA模型对我国能源需求的预测
60
实证研究2008年
三阶差分的自相关图和偏自相关图如下:1.3纯随机性检验
利用LB统计量对已经平稳的三阶差分序列进行
纯随机性检验,检验结果显示:能以95%的把握拒绝序列纯随机的原假设。因而可以认为差分以后的能源消费总量序列不属于纯随机波动。该序列不仅可以视为是平稳的,而且还蕴含着值得提取的信息,可以用来建立模型。
2
数据建模
2.1
长自回归模型初探
鉴于三阶差分的阶数偏高,首先尝试对原始序列
进行长自回归的建模。由于差分过程本身就是一个自回归的过程,从上面的分析来看,序列有可能建立长自回归模型。
在SPSS软件中试做AR(6)模型,结果显示残差的白噪声检验通过,但是滞后4、5、6期的参数不显著。经过分析得出可能是近三期的系统响应对当前期的响应影响较大的缘故,于是改作AR(3)模型。同样模型的显著性很好,但是只有延迟一期的参数通过了显著性检验。
综合前面的结果,发现对原始序列进行长自回归拟合的方法行不通。用三阶差分的序列来建模。
2.2模型定阶
三阶差分的自相关图和偏自相关图均显示了延
迟2阶的自相关系数和偏自相关系数在2倍标准差范围之外,其他阶数的系数都在2倍标准差范围内波动,而且,两图都显示了一定的拖尾性。根据这个特点可以判断该序列需拟合ARIMA(?,0,?)模型,即ARMA(?,?),所未知的是p、q阶数。可以尝试建立多种模型,直到模型通过显著性检验,以及模型参数通过显著性检验;在检验通过的模型中,选择易于解释的模型,同时兼顾拟合优度统计量。
“对多数数据取2或以下都能满足拟合混合模型
的需要。”[5]故本文尝试对序列分别建立ARMA(1,0)、
ARMA(1,1)、ARMA(1,2)、ARMA(2,0)、ARMA(2,1)、ARMA(2,2)、ARMA(0,1)、ARMA(0,2)模型,然后分别对其进行参数的显著性检验以及残差的白噪声检验。结果整理在下表中:
由上面建立的8个时序模型发现,p、q不应太小,否则模型不能通过显著性检验;q对模型的显著性贡献很小,从上表可以看出θi的显著性均很差;在模型的自回归部分,延迟一阶的参数不显著,延迟两阶的参数
非常显著。所以可以考虑建立疏系数模型AR(0,2)。
模型
拟合结果
参数的显著性检验
残差的白噪声检验
ARMA(1,0)!1不显著不通过ARMA(1,1)!1θ1不显著不通过ARMA(1,2)!1θ1θ2不显著通过ARMA(2,0)!1不显著!2显著通过ARMA(2,1)!1θ1不显著!2显著通过ARMA(2,2)!1!2θ1θ2不显著通过ARMA(0,1)θ1不显著不通过ARMA(0,2)
θ1θ2不显著
通过
2.2.1ARIMA建模
将三阶差分的序列记为X3,作为因变量,将其滞后两阶的序列记为LagsX3,作为自变量,进行回归(即做ARIMA(0,0,0)拟合模型)。由于滞后两期,所以
LagsX3序列中产生了缺省值,本文进行了如下处理:将X3序列删除了前两个样本值。回归结果如下:
!3xt=16.406-0.773B2!3xt+εtSE(854.919)(0.171)t
(0.019)(-4.509)
可见,回归系数已通过显著性检验。其中最大的绝对误差为8170.396。进一步检验模型的显著性,结果如下表所示:
基于ARIMA模型对我国能源需求的预测
第9期杜雨潇:基于ARIMA
模型对我国能源需求的预测
61
延迟阶数
LB统计量检验
LB统计量值
P值63.9180.688127.7980.80116
10.623
0.832
上表显示残差序列已经是一个白噪声序列,没有任何有价值的统计信息了。至此,模型已通过全部检验。分析结果表明,此模型是比较优良的。
2.2.2ExpertModeler建模
出于进一步优化模型的目的,将X3作为因变量,将其滞后两阶的序列记为LagsX3作为自变量,利用
SPSS的ExpertModeler模块进行回归,让软件自动选择最优模型进行拟合。结果显示最优的模型为ARIMA
(2,0,0)。结果如下:
!3
xt=-0.752B2
!3
xt+εt
SE(0.167)t
(-4.505)
其中,模型的最大绝对误差为8253.345。同样,回归系数已通过显著性检验,需进一步检验模型的显著性。
残差的白噪声检验如下:
延迟阶数
LB统计量检验
LB统计量值
P值63.7730.7071211.8150.46116
13.958
0.602
可见,该模型也通过全部检验,显著性也比较强。
2.2.3两模型的对比
将两个模型的部分信息统 …… 此处隐藏:948字,全部文档内容请下载后查看。喜欢就下载吧 ……