o u 一 ∞ n 实证研究 2008生 三阶差分的自相关图和偏自相关图如下: DIFF( ̄,3) 2.2模型定阶 三阶差分的自相关图和偏自相关图均显示了延 ■Coefficient ——迟2阶的自相关系数和偏自相关系数在2倍标准差范 围之外,其他阶数的系数都在2倍标准差范围内波动, 而且,两图都显示了一定的拖尾性。根据这个特点可 Upper Confidence Limit Lower Confidence Limit —以判断该序列需拟合ARIMA(?,0,?)模型,即ARMA(?,?), 所未知的是P.q阶数。可以尝试建立多种模型,直到模 I 2 3 4 5 6 7 8 910l1l2l314l516 Lag Num ̄r DIFF(x,3) ■Coefficient ——‘ UpperConfidence Limit 一..- —Lower Confidence Limit 1 1 . I 2 3 4 5 6 7 8 910111213I4l5l6 Lag Number lI3纯随机性检验 利用LB统计量对已经平稳的三阶差分序列进行 纯随机性检验,检验结果显示:能以95%的把握拒绝 序列纯随机的原假设。因而可以认为差分以后的能源 消费总量序列不属于纯随机波动。该序列不仅可以视 为是平稳的,而且还蕴含着值得提取的信息,可以用 来建立模型。 2数据建模 2.1长自回归模型初探 鉴于三阶差分的阶数偏高,首先尝试对原始序列 进行长自回归的建模。由于差分过程本身就是一个自 回归的过程,从上面的分析来看,序列有可能建立长 自回归模型。 在SPSS软件中试做AR(6)模型,结果显示残差的 自噪声检验通过,但是滞后4、5、6期的参数不显著。经 过分析得出可能是近三期的系统响应对当前期的响 应影响较大的缘故,于是改作AR(3)模型。同样模型的 显著性很好,但是只有延迟一期的参数通过了显著性 检验。 综合前面的结果,发现对原始序列进行长自回归 拟合的方法行不通。用三阶差分的序列来建模。 型通过显著性检验,以及模型参数通过显著性检验; 在检验通过的模型中,选择易于解释的模型,同时兼 顾拟合优度统计量。 “对多数数据取2或以下都能满足拟合混合模型 的需要。”[51故本文尝试对序列分别建立ARMA(1,01、 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对模型的显著性贡献 很小,从上表可以看出Oi的显著性均很差;在模型的 自回归部分,延迟一阶的参数不显著,延迟两阶的参数 非常显著。所以可以考虑建立疏系数模型AR(0,2)。 模型 2.2.1 ARIMA建模 将三阶差分的序列记为X3,作为因变量,将其滞 后两阶的序列记为LagsX3,作为自变量,进行回归(即 做ARIMA(0,0,0)拟合模型)。由于滞后两期,所以 LagsX3序列中产生了缺省值,本文进行了如下处理: 将X3序列删除了前两个样本值。回归结果如下: =16.406—0.773B V 3x +oo SE(854.919)(0.171) t (0.019)(一4.509) 可见,回归系数已通过显著性检验。其中最大的 绝对误差为8170.396。进一步检验模型的显著性,结果 如下表所示: 维普资讯 http://www.cqvip.com 第9期 延迟阶数 杜雨潇:基于ARIMA模型对我国能源需求的预测 LB统计量检验 20,000 .——-——61 LB统计量值P值 l0,000 —.—-Observed Fit UCL LCL Forecast ・一鱼 o . 上表显示残差序列已经是一个白噪声序列,没有 、√ 、 、 V 任何有价值的统计信息了。至此,模型已通过全部检 验。分析结果表明,此模型是比较优良的。 2.2.2 Expert Modeler建模 lOL000 . 20,000 出于进一步优化模型的目的,将X3作为因变量, 将其滞后两阶的序列记为LagsX3作为自变量,利用 SPSS的Expert Modeler模块进行回归,让软件自动选 择最优模型进行拟合。结果显示最优的模型为ARIMA (2,0,01。结果如下: 3x产一0.752B 3x + (单位:万吨标准煤) SE t (0.167) (一4.505) 4结论及建议 依照本文构建的模型进行预测,到2010年我国的 能源消费总量将超过29亿吨标准煤。 我国未来经济的发展将会面临严峻的能源形势。 为了缓解巨大的能源需求压力,今后应联合各方努 力,在保证经济稳定发展的前提下,严格控制能源消 费,加强能源利用的监督体制。本文试从以下几方面 给出一些对策及建议: 其中,模型的最大绝对误差为8253.345。同样,回 归系数已通过显著性检验,需进一步检验模型的显著 性。 残差的白噪声检验如下: 延迟阶数 LB统计量检验 LB统计量值P值 第一,积极开发利用新能源。加强对太阳能、风能、 生物能、潮汐能以及地热能等的利用,统筹发展各种 可见,该模型也通过全部检验,显著性也比较 强。 2.2.3两模型的对比 新能源以降低全社会对传统能源的依赖;第二,扩大 节能技术的推广。节能是我们应当密切关注解决能源 问题的出路之一。通过提高能源的利用率,降低能耗, 达到节约资源的目的;第三,政府加强对能源利用的 监督机制。对于那些能源消耗过高的产业给以必要的 约束和管制。同时,在全社会推广节能,通过增强全体 社会成员的节约意识,达到能源的合理、可持续利用。 将两个模型的部分信息统计在下表进行比较: 从比较结果来看,无论从BIC准则还是R— 参考文献 squared方面来讲,2.2.2模型均优于2、2、1模型,所以 选择2、2、2模型为最终拟合的模型: V 3x =一0、752B + 【l】崔民选.2007中国能源发展报告[^删E京:社会科学文献, ̄2007(1) 2007 【2】陶爱元.中国经济增长对能源的依赖程度分析『J】.统计与决策. 下面是预测值与观测值在时序图中的分布: 【3】熊国强,刘海磊.我国能源消费的组合预测模型【J1.统计与决策. 2007(2) 3模型预测 利用上述拟合的时间序列模型对未来四年的能 源消费总量进行预测,预测结果如下: 【4 汪陶,41张意翔.湖北能源需求预测研究【J】.统计与决策.2007(3) 【5】宇传华.SPSS与统计分析【M】.北京:电子工业出版社,2007 【6 王燕.61应用时问序列分析【M E京:中国人民大学出版社,2005 (责任编辑李本钊)
因篇幅问题不能全部显示,请点此查看更多更全内容