统计预测与决策论文!

1、各种预测模型在全国能源消耗总量中的应用091124104丁灯摘要：能源影响着我国社会经济的稳定持续发展，对未来能源消耗的准确预测具有重要的意义。本文以我国1978-2008年的全国能源消耗总量数据为基础，建立了arima预测模型、灰色预测模型、三次多项式预测模型和基于这三种模型的组合模型，并进行了精度比较，最后选择最优的组合预测模型对2009-2011年的全国能源消耗总量进行预测。关键词：arima模型；灰色预测模型；三次多项式；组合模型；能源消耗1 引言：能源是国民经济发展和人民生活水平提高的重要物质基础，能源短缺曾经长期制约我国经济的发展。近几年由于能源工业的发展，短缺局面虽然得到了缓解，

2、但从长远来看能源供需形势仍然非常严峻，因此做好未来能源消费预测分析，为能源规划及政策的制定提供科学的依据，对于保持我国社会经济健康、持续、稳定发展具有重要的理论与现实意义。本文利用中国统计年鉴得到31期全国能源消耗总量y的时间序列如下表一所示：表一：全国能源消耗总量（单位：万吨标准煤）年份197819791980198119821983198419851986y571445858860275594476206766040709047668280850年份198719881989199019911992199319941995y8663292997969349870310378310917011

3、5993122737131176年份199619971998199920002001200220032004y138948137798132214133831138552.6143199.2151797.3174990.3203226.7年份2005200620072008y2246822462702655832850002 预测方法介绍2.1 arima模型的基本原理arima模型是box和jenkins1970年提出的以随机理论为基础的时间序列分析方法，又称为“box-jenkins模型”，这以模型在经济领域的预测分析中得到了广泛的应用。时间序列是依赖时间t的一组随机变量，构成该时序的单个

4、序列值虽然具有不确定性，但对整个时间序列来说，它的变化却有一定的规律性，可以用相应的数学模型来近似描述。arima模型有三种基本类型：自回归模型、移动平均模型、单整自回归移动平均模型。单整是指将一个时间序列有非平稳性变为平稳性所要经过的差分的次数，这是对非平稳时间序列进行时间序列分析的必经步骤。假设一个随机过程含有d个单位根，其经过d次差分之后可以变换为一个平稳的自回归移动平均过程。则该随机过程称为单整自回归移动平均模型。模型中ar称为自回归分量，p为自回归分量的阶数；ma为移动平均分量，q为移动平均分量的阶数；i为差分，d为使时间序列具有平稳性所需要的差分次数。p阶自回归过程ar(p)的一般

5、表达式为： 其中白噪声过程。q阶的移动平均过程ma（q）可以表示为：，为白噪声过程。arima( p,d, q)模型一般表达式为：2.2 灰色预测法灰色预测法是一种对含有不确定因素的系统进行预测的方法。一般是利用时间序列数据，通过建立gm（1,1）模型进行预测。灰色预测模型的预测步骤如下：（1） 首先对原始时间序列数据，做一次累加生成，得到新的序列（2） 利用一次累加生成序列拟合微分方程：，得到参数和（3） 解微分方程得到预测模型函数：（4） 将得到的序列进行一次累减得到预测序列（5） 利用历史数据对数据模型进行精度检验，若通不过检验，则利用残差对原模型进行修正。（6） 通过预测方程进行预测。

6、2.3组合预测模型不同的预测方法根据相同的信息，往往会提供不同的结果，如果简单的将误差较大的一些方法舍弃掉，将会丢弃一些有用的信息，使得模型的精度不高。组合预测法是指通过建立一个组合预测模型，把多种预测方法所得到的预测结果进行综合。由于组合模型能够较大限度地利用各种预测样本信息，所以它比单项预测模型考虑问题更系统、更全面，因而能够有效地减少单个预测模型受随机因素的影响，可以提高预测的精度和稳定性。3 全国能源消耗总量的实证分析3.1 建立arima模型3.11平稳化处理用arima模型拟合的时间序列必须是平稳的，如果序列不平稳，则要通过差分或序列变换等先将序列平稳化。绘制原始序列的时序图得到图

7、形如图一所示：图一：y时序图由图可从直观上看出原始序列存在明显的长期递增趋势，原始序列不平稳。利用软件eviews6.0，运用单位根检验方法对序列进行平稳性检验发现原始序列确定不平稳，因此本文先对该序列取对数，令，然后对yl进行差分，差分两次之后得到平稳的序列ylii。单位根检验结果如下表二所示：表二：单位根检验结果adf值p值临界值a=1%临界值a=5%临界值a=10%原始序列y1.7585820.9995-3.68919-2.97185-2.62512取对数后yl0.6275170.988-3.68919-2.97185-2.62512一阶差分yli-2.958250.0529-3.724

8、07-2.98623-2.6326二阶差分ylii-4.765430.0007-3.69987-2.97626-2.62742由表可知，取对数后再做两次一步差分之后的序列yliiadf检验的p值为0.0007，小于0.05，因此拒绝序列非平稳的原假设，可以接受序列为平稳的备择假设。在此可知arima模型定阶为d=2。也可由图标直观说明序列的平稳性，作出ylii的时序图如图一所示：图二：ylii时序图时序图也显示ylii序列平稳，结果与单位根检验相符。3.12模型定阶模型定阶的方法有多中，本文选择基于自相关函数和偏自相关函数的定阶方法确定模型的阶数。首先，考察平稳序列ylii的自相关图和偏自相关

9、的性质，为拟合模型定阶，自相关函数(acf)和偏自相关函数（pacf）图形如图三所示;图三：由图形可知，对处理后的序列ylii可以选择建立arima(1,2,1)、arima(1,2,2)、arma(2,2,1)、arima(2,2,2)4种模型。分别拟合这四种模型得到结果如表三所示：表三：4种模型拟合结果模型变量估计系数t统计量伴生概率p值调整aicscarima(1,2,1)ar(1)-0.37510 -1.22811 0.23040 0.17179 -4.13804 -4.04288 ma(1)0.74524 3.23898 0.00330 arima(1.2,2)ar(1)0.3263

10、2 1.53493 0.13740 0.27357 -4.23696 -4.09422 ma(1)-0.25253 -1.74264 0.09370 ma(2)-0.73255 -5.34026 0.00000 arima(2,2,1)ar(1)-0.25552 -0.95215 0.35050 0.18863 -4.15745 -4.01347 ar(2)-0.23472 -1.22733 0.23160 ma(1)0.59424 2.30532 0.03010 arima(2,2,2)ar(1)-0.76684 -5.24699 0.00000 0.28900 -4.25798 -4.06

11、601 ar(2)-0.59005 -3.94857 0.00060 ma(1)1.24653 20.40994 0.00000 ma(2)0.92659 20.03634 0.00000 其中只有arima(2,2,2)模型的各系数通过了显著性检验，而且其模型的调整是四个模型中最大的，虽然它的aic，sc的绝对值不是最小的，arima(1,2,1)模型的aic最小，arima(2,2,1)的sc最小，但它们的系数都没有通过显著性检验，而且调整也较小，所以本文选择arima(2,2,2)模型，模型表达式为：3.13 模型检验用arima(2,2,2)模型做拟合得到残差序列，对残差序列进行自相关

12、和偏自相关分析，得到结果如图三所示：图三由图可知。残差序列p值几乎都是大于0.05的，说明残差序列近于白噪声，基本没有可提取的信息了，模型已经提取了有规律的信息，说明模型拟合效果较好。3.14模型预测利用1978-2008年的时间序列建立的arma(2,2,2)模型：来预测2004-2010年的能源消耗总量，结果如表四所示：表四：2004-2010年全国能源消耗总量arima(2,2,2)模型预测值年份2004.00 2005.00 2006.00 2007.00 2008.00 2009.00 2010.00 2011.00 实际值203226.68 224682.00 246270.00

13、265583.00 285000.00 预测值202971.88 233243.50 268974.47 311049.59 358189.95 413129.35 477104.00 549931.36 绝对相对误差百分比(%)0.13 3.81 9.22 17.12 25.68 3.2 灰色模型预测根据历史数据序列，做一次累加得到生成序列，对于微分方程，构造数据矩阵b和数据向量y，解该微分方程，得到其中： 进行矩阵运算得到发展灰数a=-0.055673，内生控制灰数=45508.908，得到预测模型为：。3.21残差检验将得到的序列进行一次累减生成预测序列，将预测值与真实值比较得到绝对误差

14、序列为：=0、8516.84、7337.171、3478.379、2894.069、3479.306、4761.587、6752.808、6917.228、8466.436、10356.31、9561.971、6328.756、6120.155、5915.771、6827.269、7321.323、9152.556、9938.481、1402.44、11990.47、18629.48、22636.52、27218.27、28376.96、15499.22、1831.269、11756.32、21153.92、27578.59、33369.39计算相对误差后发现其中有些较大，甚至大于10%，最后

15、五项的平均相对误差为7.36319507%。残差检验没有通过，下面进行关联度检验。当=0.5时，关联度大于0.6时就可以通过关联度检验了。由残差的绝对序列可知=0, =33369.39。然后根据公式计算每个序列值的关联系数，再求平均得到关联度r=0.650237。本文取取=0.5，r0.6,因此模型通过了关联度检验。虽然模型通过了关联度检验，但模型没有通过残差检验，精度不够理想，所以本文对该模型进行残差修正。3.22模型修正根据模型得到的残差序列，去掉第一项得到新的序列，然后进行累加得到序列，在此基础上建立相应的模型：，模型两边求导得到残差修正项为：。原预测模型加上此修正项得到修正后的预测模型

16、为：，其中为修正系数。最后进行累减得到原始序列预测模型：。根据所得数据，利用excel软件，得到=-0.06924，=2893.212。根据修正项的计算公式，可得，这样经过残差修正后的模型为：其中。在此序列的基础上进行累减得到2004-2011年的全国能源消耗总量的预测值如下表五所示：表五：2004-2011年全国能源消耗总量灰色模型预测值年份20042005200620072008200920102011实际值203226.7224682246270265583285000预测值201439.6212966.946225154.6238040.3251664.1266068.2281297.

17、3297398.4相对误差百分比(%)0.875.218.5710.3711.693.3三次多项式当预测对象依时间变化呈现某种上升或下降的趋势，并且无明显的季节波动，又能找到一条合适的函数曲线反应这种变化趋势时，就可用时间t为自变量，时序数值y为因变量，建立趋势模型：。 三次多项式预测模型是趋势外推预测模型中的一种比较常用的模型。为了探究时序数据变化的趋势，作出全国能源消费总量y的散点图如图四所示：图四：散点图由图可知全国能源消耗总量没有明显的季节波动，而且大致呈抛物线或指数的趋势上升，但经过精度比较之后，三次多项式的拟合的最好，所以本文采用三次多项式模型进行拟合。利用软件eviews 进行回

18、归得到结果如下表六所示：表六：三次多项式模型结果变量系数标准差t统计量p值常数项34209.538549.9254.0011500.0004t10136.882277.3004.4512700.0001t2-646.4941163.9748-3.9426440.0005t318.582983.3716555.5115310.0000其中模型的f统计量为345.9252，相伴概率为0.00，小于0.05，所以拒绝原假设，即模型是显著的。而且各系数的显著性检验的p值都小于0.05，通过了显著性检验。模型的调整的为0.971825，接近于1，说明模型拟合效果好。因此，建立三次多项式模型如下：用该模型

19、对2004-2011年的全国能源消费总量进行预测得到结果如表七所示：表七：2004-2011年全国能源消费总量三次多项式预测值年份20042005200620072008200920102011真实值203226.7224682246270265583285000预测值202380.4219124.89237698.4258212.3280778.2305507.6332511.9361902.6相对误差百分比（%）0.412.473.482.771.483.4组合模型以上三种方法有些精度较高，有些精度较低，但由于各自理论的缺陷，会使预测产生系统偏差。为了更充分的提取样本信息，提高预测精度，本

20、文将以上三种预测模型进行组合得到一个组合的预测模型。组合模型的权重有多种，可以等概率取权重，可以采用最小二乘估计法取权重，也可以使权重，其中为残差标准差，m为单一模型的个数，或者使权重，其中为第i个单一模型误差平方和，m为单一模型个数。本文采用第三种取取权重的方法，即。利用excel软件求得权重分别为：=0.466858，=0.171932824，=0.361209557根据预测模型：对2004-2011年的全国能源消耗总量进行预测结果如下表八所示：表八：2004-2011年全国能源消耗总量组合预测值年份20042005200620072008200920102011真实值203226.68224682246270265583285000预测值202142.34219328.12238544.92260010.57283571.66309765.5338872.9370931.5相对误差百分比