灰色预测模型理论与应用

. . .

灰色预测模型理论及其应用

灰色系统理论认为对既含有已知信息又含有未知或非确定信息的系统进行预测，就是对在一定方位内变化的、与时间有关的灰色过程的预测. 尽管过程中所显示的现象是随机的、杂乱无章的，但毕竟是有序的、有界的，因此这一数据集合具备潜在的规律，灰色预测就是利用这种规律建立灰色模型对灰色系统进行预测.

灰色预测模型只需要较少的观测数据即可，这和时间序列分析，多元回归分析等需要较多数据的统计模型不一样. 因此，对于只有少量观测数据的项目来说，灰色预测是一种有用的工具.本文主要围绕灰色预测GM(1,1)模型及其应用进行展开。

一、灰色系统及灰色预测的概念

1.1灰色系统

灰色系统产生于控制理论的研究中。 若一个系统的内部特征是完全已知的，即系统的信息是充足完全的，我们称之为白色系统。

若一个系统的内部信息是一无所知，一团漆黑，只能从它同外部的联系来观测研究，这种系统便是黑色系统。

灰色系统介于二者之间，灰色系统的一部分信息是已知的，一部分是未知的。 区别白色和灰色系统的重要标志是系统各因素间是否有确定的关系。 特点：灰色系统理论以“部分信息已知、部分信息未知”的 “小样本”、“贫信息”不确定型系统的研究对象。 1.2灰色预测

灰色系统分析方法是通过鉴别系统因素之间发展趋势的相似或相异程度，即进行关联度分析，并通过对原始数据的生成处理来寻求系统变动的规律。生成数据序列有较强的规律性，可以用它来建立相应的微分方程模型，从而预测事物未来的发展趋势和未来状态。灰色预测是用灰色模型GM(1,1)来进行定量分析的，通常分为以下几类：

(1) 灰色时间序列预测。用等时距观测到的反映预测对象特征的一系列数量（如产量、销量、人口数量、存款数量、利率等）构造灰色预测模型，预测未来某一时刻的特征量，或者达到某特征量的时间。

(2) 畸变预测（灾变预测）。通过模型预测异常值出现的时刻，预测异常值什么时候出现在特定时区内。

(3) 波形预测，或称为拓扑预测，它是通过灰色模型预测事物未来变动的轨迹。 (4) 系统预测，是对系统行为特征指标建立一族相互关联的灰色预测理论模型，在预测系统整体变化的同时，预测系统各个环节的变化。

上述灰预测方法的共同特点是： （1）允许少数据预测；

（2）允许对灰因果律事件进行预测，比如 灰因白果律事件：在粮食生产预测中，影响粮食生产的因子很多，多到无法枚举，故为灰因，然而粮食产量却是具体的，故为白果。粮食预测即为灰因白果律事件预测。

白因灰果律事件 ：在开发项目前景预测时，开发项目的投入是具体的，为白因，而项目的效益暂时不很清楚，为灰果。项目前景预测即为灰因白果律事件预测。

..........

. . .

（3）具有可检验性，包括：建模可行性的级比检验（事前检验），建模精度检验（模型检验），预测的滚动检验（预测检验）。

二、GM（1,1）模型

2.1GM(1,1)模型

GM(1,1)模型是基于灰色系统的理论思想，将离散变量连续化，用微分方程代替差分方程，按时间累加后所形成的新的时间序列呈现的规律可用一阶线性微分方程的解来逼近，用生成数序列代替原始时间序列，弱化原始时间序列的随机性，这样可以对变化过程作较长时间的描述，进而建立微分方程形式的模型. 其建模的实质是建立微分方程的系数，将时间序列转化为微分方程，通过灰色微分方程可以建立抽象系统的发展模型. 经证明，经一阶线性微分方程的解逼近所揭示的原始时间数列呈指数变化规律时，灰色预测GM(1,1)模型的预测将是非常成功的.

2.2GM(1,1)模型的建立

GM(1,1)模型是指一阶，一个变量的微分方案预测模型，是一阶单序列的线性动态模型，用于时间序列预测的离散形式的微分方程模型.

模型符号含义为

G M （1， 1） Grey Model 1阶方程 1个变量

设时间序列X?0?有n个观察值，X?0??x?0??1?,x?0??2?,?,x?0??n?，为了使其成为

?有规律的时间序列数据，对其作一次累加生成运算，即令

x从而得到新的生成数列X?1??1??t???x?0??n?

n?1t，X?1??x?1??1?,x?1??2?,?,x?1??n?,称

?x(0)(k)?ax(1)(k)?b

为GM(1,1)模型的原始形式。 新的生成数列X?1?一般近似地服从指数规律. 则生成的离散形式的微分方程具体的形式为

dx?ax?u dt即表示变量对于时间的一阶微分方程是连续的. 求解上述微分方程，解为

x(t)?ce?a(t?1)?u a当t=1时，x(t)?x(1)，即c?x(1)? ..........

u，则可根据上述公式得到离散形式微分方程的具体a. . .

形式为

u?u?x?t???x?1???e?a?t?1??

a?a?其中，ax项中的x为

dx的背景值，也称初始值；a，u是待识别的灰色参数，a为发展系dt数，反映x的发展趋势；u为灰色作用量，反映数据间的变化关系.

按白化导数定义有

dxx(t?t)?x(t)?lim t?0dtt显然，当时间密化值定义为1时，当t?1时，则上式可记为

dx?lim(x(t?t)?x(t)) dtt?1这表明

dx是一次累减生成的，因此该式可以改写为 dtdx?x(1)(t?1)?x(1)(t) dt当t足够小时，变量x从x(t)到x(t?t)是不会出现突变的，所以取x(t)与x(t?t)的

(1)平均值作为当t足够小时的背景值，即x成序列）将其值带入式子，整理得

1(1)(1)??x(t)?x(t?1)??（紧邻均值（MEAN）生2?1(1)(1)x(0)(t?1)??a?x(t)?x(t?1)??u（GM(1,1)模型的均值形式） ??2由其离散形式可得到如下矩阵：

1(1)??(1)???x(1)?x(2)???2??x(0)(2)???1?(0)?(1)(1)???x(2)?x(3)?x(3)??a??2?????u

?????????x(0)(n)????1(1)?(1)?????x(n?1)?x(n)???2?(0)(0)令 Y???x(2),x(3),,x(n)??

(0)T ..........

. . .

1(1)?(1)1????x(1)?x(2)???2???1???x(1)(2)?x(1)(3)??1?? B??2??????1(1)?(1)?????x(n?1)?x(n)?1??2????au?

称Y为数据向量，B为数据矩阵，?为参数向量. 则上式可简化为线性模型：

Y?B?

由最小二乘估计方法得

?1T?a?T??????BB?BY

?u?T上式即为GM(1,1)参数a,u的矩阵辨识算式，式中BB义逆矩阵.

将求得的a，u值代入微分方程的解式，则

?T??1BTY事实上是数据矩阵B的广

u??a(t?1)u(1)?(1)(t)??xx(1)?? ??ea?a?其中，上式是GM(1,1)模型的时间响应函数形式，将它离散化得

u??a(t?1)u(0)?(1)(t)??xx(1)?? ??eaa???对序列x?1??t?再作累减生成可进行预测. 即

?(0)(t)?x?(1)(t)?x?(1)(t?1)x u??(0)a?a(t?1)??x(1)???1?e?ea??上式便是GM(1,1)模型的预测的具体计算式.

2.3 GM(1,1)模型的检验

GM(1,1)模型的检验包括残差检验、关联度检验、后验差检验三种形式. 每种检验对应不同功能：残差检验属于算术检验，对模型值和实际值的误差进行逐点检验；关联度检验属于几何检验范围，通过考察模型曲线与建模序列曲线的几何相似程度进行检验，关联度越大模型越好；后验差检验属于统计检验，对残差分布的统计特性进行检验，衡量灰色模型的精度. ? 残差检验

残差大小检验，即对模型值和实际值的残差进行逐点检验.

设模拟值的残差序列为e(t)，则

(0) ..........