数学建模+灰色预测模型+MATLAB

§12.5 灰色预测

我们通常所说的系统是指：由客观世界中相同或相似的事物和因素按一定的秩序相互关联、相互制约而构成的一个整体．例如：工程技术系统、社会系统、经济系统等．如果一个系统中具有充足的信息量，其发展变化的规律明显、定量描述方便、结构与参数具体，则这种系统通常称为白色系统．如果一个系统的部特征全部是未知的，则称此系统为黑色系统．如果系统部信息和特征是部分已知的，另一部分是未知的，这种系统称为灰色系统．例如：社会系统、农业系统、经济系统、气象系统、生物系统等．对于这类系统，部因素难以辨识，相互之间的关系较为隐蔽，人们难以准确了解这类系统的行为特征．因此，对于这类问题进行定量描述，即建立模型难度较大．区别白色系统与灰色系统的重要标志是系统各因素之间是否具有确定的关系．

灰色系统分析方法主要是根据具体灰色系统的行为特征数据，充分利用数量不多的数据和信息寻求相关因素自身与各因素之间的数学关系，建立相应的数学模型．目前，灰色系统理论在实际中已得到了广泛的应用，例如：在工程技术、经济管理、气象预报以及政治、社会、工业、农业等领域都取得了一定的应用成果．

我们往往要对农业问题、商业问题等做未来的预测工作，另外，进行军事战争以及治理生态环境也需对未来的发展情形做一可靠的分析，这就产生了灰色预测．灰色预测是对灰色系统问题进行未来的预测，实际问题中，应用最多的灰色预测模型是以GM(1,1)（即GM(1,N)当N=1时的特例）模型为基础的．

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

设X(0)=(X(0)(1),X(0)(2),…,X(0)(n))，做1－AGO，得

则GM(1,1)模型相应的微分方程为：

（1）

式中：a称为发展灰数；μ称为生控制灰数．

设=(a,μ)T，按最小二乘法得到

（2）

其中

易求得，方程（1）的解为

（3）

例4 100m成绩预测

1983～1990年世界男子和中国女子100m最好成绩如表6．

表6 各年度最好成绩

记世界男子100m成绩的原始数列为

建立GM(1,1)模型，即按式（1）、（2）、（3）得到预测模型为

由预测模型得预测值为

年份 模型预测值/s

1991 9.92

1992 9.91

2000 9.85

记中国女子的原始数列为

同样建立GM(1,1)模型，得到预测模型为

从而得到中国女子100m成绩的预测值

年份 模型预测值/s

1991 11.30

1992 11.24

2000 10.85

12.5.2 模型检验

灰色预测检验一般有残差检验、关联度检验和后验差检验．

（1）残差检验

按预测模型计算，并将累减生成，然后计算原始序列X(0)(i)

与的绝对误差序列及相对误差序列．

（2）关联度检验

定义1 选取参考数列

其中k表示时刻．假设有m个比较数列

则称

（1）

为比较数列Xi对参考数列X0在k时刻的关联系数，其中ρ∈[0,1]为分辨系数，一般取ρ=0.5．称式（1）中| X0(k)-Xi(k)|、 | X0(k)-Xi(k)|分别为两级最小差和两级最大差．由（1）式易看出，ρ越大，分辨率越大；ρ越小，分辨率越小．

式（1）定义的关联系数是描述比较数列与参考数列在某时刻关联程度的一种指标，由于各个时刻都有一个关联数，因此信息显得过于分散，不便于比较，为此我们给出以下定义

定义2 称

（2）

为数列Xi对参考数列X0的关联度．

由式（2）易看出，关联度是把各个时刻的关联系数集中为一个平均值，也就是把过于分散的信息集中处理．根据前面所述关联度计算方法计算出与原始序列X(0)(i)的关联系数，然后计算出关联度，根据经验，当ρ=0.5时，关联度大于0.6便满意了．

（3）后验差检验

1.计算原始序列标准差：

2.计算绝对误差序列的标准差：

3.计算方差比：

4.计算小误差概率：

令则．

表7 检验标准

若残差检验、关联度检验和后验差检验都能通过，则可以用所建模型进行预测；若用原始时间序列X(0)建立的GM(1,1)模型检验不合格或精度不理想时，这时要对建立的GM(1,1)模型进行修正或提高模型的预测精度．其修法如下：

设原始时间序列X(0)建立的GM(1,1)模型为

可获得生成序列X(1)的预测值，即对于，有预测序列，定义残差为

若取j=i,i+1,…,n,则与X(1)及对应的残差序列为

为便于计算上式改写为

e(0)的累加生成序列为

e(1)可建立相应的GM(1,1)模型：

的导数加上修正，得修正模型：

其中为修正系数．

最后给出经过残差修正的原始序列预测模型：

§12.6 灰色预测模型案例

一、问题描述

表8给出了市1991年－1996年国生产总值总消费资料．

生产决定消费，国生产总值总消费决定了居民的消费水平，为此很有必要对国生产总值总消费进行科学预测，分析国生产总值总消费发展趋势，为宏观经济政策的制定提供重要的参考．

表8 国生产总值总消费 单位：亿元

试根据表8的资料，建立市国生产总值总消费的灰色预测模型GM(1,1)，并预测市1998年国生产总值总消费．

二、模型的建立及求解

1.令X(0)(1),X(0)(2),…,X(0)(6)对应于原始序列数据．

第一步，构造累加生成序列：

第二步，构造数据矩阵B和数据向量Y1：

第三步，计算BTB, (BTB)-1, BTY1：

，

，

即

第四步，得出预测模型：

三、模型检验

第五步，进行关联度检验：

（1）计算：

（2）计算关联系数：

，是基本满足时，r>0.57的．

所以关联度检验通过．

第六步，后验差检验：

（1）计算：

（2）计算残差的均值：

残差的标准差：

，

所有ek都小于S0，故P=1,C<0.35．

所以后验差检验通过．

第七步，残差检验：

（1）计算

（2）累减生成序列：

（3）计算绝对误差序列及相对误差序列：

绝对误差序列△(0)={0,53.87,26.28,69.82,95.37,33.48}

相对误差序列Ф={0,11.3%,3.87%,7.99%,8.79%,2.67%}

相对误差序列中有的相对误差很大，所以要对原模型进行残差修正以提高精度．

（4）利用残差对原模型进行修正：

取 e(0)={53.87,26.28,69.82,95.37,33.48}

e(1)={53.87,80.15,149.94,245.34,278.82}

得

最后得修正模型为：

其中．

表9 修正后的残差计算表

因此，可用上述经过残差修正后的模型来预测市1998年国生产总值总消费：

故市1998年国生产总值总消费预测值为：

（亿元）

注：灰色预测是通过对原始数据的处理和灰色模型的建立，发现、掌握系统发展规律，对系统的未来状态作出科学的定量预测．一个模型要经过多种检验才能判定其是否合理有效，只有通过检验的模型才能用作预测．

本文来源：https://www.2haoxitong.net/k/doc/ff8abfb40a1c59eef8c75fbfc77da26924c5962e.html