GM(1,1)模型应用及残差修正.doc

一．GM(1,1)预测模型应用举例 灰色预测是基于GM(1,1)预测模型的预测，按其应用的对象可有四种类型： 数列预测。这类预测是针对系统行为特征值的发展变化所进行的预测。 灾变预测。这类预测是针对系统行为的特征值超过某个阙值的异常值将在何时出现的预测。 季节灾变预测。若系统行为的特征有异常值出现或某种事件的发生是在一年中的某个特定的时区，则该预测为季节性灾变预测。 拓扑预测。这类预测是对一段时间内系统行为特征数据波形的预测。 例1（数列预测）：设原始序列 试用GM(1,1)模型对进行模拟和预测，并计算模拟精度。 解：第一步：对进行一次累加，得 第二步：对作准光滑性检验。由 得。 当k3时准光滑条件满足。 第三步：检验是否具有准指数规律。由 得 当k3时，，准指数规律满足，故可对建立GM(1,1)模型。 第四步：对作紧邻均值生成，得 于是 第五步：对参数列进行最小二乘估计。得 第六步：确定模型 及时间响应序列 第七步：求的模拟值 第八步：还原求出的模拟值。由 得 第九步：检验误差。由下表可算出残差平方和： 误差检验表 序号 实际数据 模拟数据 残差 - 相对误差 2 3 4 5 3.278 3.337 3.390 3.679 3.2300 3.3545 3.4817 3.6136 0.0460 -0.0175 -0.0917 0.0654 1.40% 0.52% 2.71% 1.78% 平均相对误差 1.6025% 第十步：预测 例2 （灾变预测）：某企业生产用原料属受自然灾害影响较大的农产品。一般来说，自然灾害的发生有其偶然性，但对历史数据的整理，仍可发现一定的规律性。为尽量减少生产不受自然灾害的影响，该企业希望了解影响原料供应的规律性并提前做好原料储备，所收集数据见下表，并规定每亩平均收获量小于320千克时为欠收年份，将影响原料的正常供应，现应用灰色灾变预测来预测下次发生欠收的年份。 原料收获统计表 年份 1991 1992 1993 1994 1995 1996 1997 1998 1999 收获量(千克) 390.6 412 320 559 380 542 553 310 561 年份 2000 2001 2002 2003 2004 2005 2006 2007 收获量(千克) 300 632 540 406.2 314 576 587 318 第一步：将上表中年份用序号替换，并找出收获量小于320千克的年份序号形成初始序列。 本例初始序列： 一次累加生成序列： 的紧邻均值生成序列： 第二步：按建GM(1,1)模型。 第三步：预测 当t=6时， 因此，下次发生收获量小于320千克的年份为：2011年至2012年，即四至五年后将出现欠收年份。 其他预测类型见参考书。残差GM(1,1)模型 当GM(1,1)模型精度不符合要求时，可使用残差序列建立GM(1,1)模型，对原来模型进行修正，以提高精度。 定义4 设 其中，-为的残差序列。若存在k0，满足 1. 2.，则称 为可建模残差尾段，仍记为 命题1 设为可建模残差尾段，其一次累加序列 的GM(1,1)模型的时间响应式为 则残差尾段的模拟序列为 其中 定义5 若用修正则称修正后的时间响应式 为残差修正GM(1,1)模型，简称残差GM(1,1)。其中残差修正值 的符号应与残差尾段的符号保持一致。 定义6 若则相应的残差修正时间响应式 称为累减还原式的残差修正模型。 例题 湖北省云梦县油菜发病率数据为 建立GM(1,1)模型，得时间响应式为 作累减还原，得 检验其精度：列出误差检验表 误差检验表 序号 实际数据 模拟数据 残差 - 相对误差 2 3 4 5 6 7 8 9 10 11 12 13 20 40 25 40 45 35 21 14 18 15.5 17 15 35.6704 33.4303 31.3308 29.3682 27.5192 25.7901 24.1719 22.6534 21.2307 19.8974 18.6478 17.4768 -15.6704 6.5697 -6.3308 10.6318 17.4808 9.2099 -3.1719 -8.65