灰色模型介绍及应用.docVIP

灰色模型介绍及应用

PAGE

PAGE13

第十章灰色模型介绍及应用（徐利艳天津农学院2.4万字）

10.1灰色理论基本知识

10.1.1概言

10.1.2有关名词概念

10.1.3GM建模机理

10.2灰色理论模型应用

10.2.1GM（1，1）模型的应用——污染物浓度问题

10.2.2GM（1，1）残差模型的应用——油菜发病率问题

10.2.3GM模型在复杂问题中的应用——SARS疫情问题

10.2.4GM（1，n）模型的应用——因素相关问题

本章小结

思考题

推荐阅读书目

第十章灰色模型介绍及应用

10.1灰色理论基本知识

10.1.1概言

客观世界的很多实际问题，其内部的结构、参数以及特征并未全部被人们了解，人们不可能象研究白箱问题那样将其内部机理研究清楚，只能依据某种思维逻辑与推断来构造模型。对这类部分信息已知而部分信息未知的系统，我们称之

利用常数变易法解得，通解为

若初始条件为，则可得到微分方程的特解为

或时间响应函数

其中白化微分方程中的项中的为的背景值，也称为初始值;为常数（有时也将写成）。

按白化导数定义有差分形式的微分方程，即

显然，当时间密化值定义为1，即当时，上式可记为

记为离散形式

这显然表明是一次累计生成，因此上述方程可改写为

这实际也表明，模型是以生成数（是以的一次累加）为基础的。

当足够小时，到不会发生突变，因此可取与的平均值作为时的背景值，因此，背景值便可记为

或

于是白化的微分方程可改写为

或

即

因此，上述方程可以改写为矩阵方程形式，即

引入下列符号，设

于是便有

令

则

解得

将求解得到的代入微分方程的解式（也称时间响应函数），则

由于，因此求导还原得

上述两式便为GM（1，1）的时间响应式，及灰色系统预测模型的基本算式，当然上述两式计算结果只是近似计算值。

为简记，一般可以将GM（1，1）的建模过程记为

10.2灰色理论模型应用

10.2.1GM（1，1）模型的应用——污染物浓度问题

GM（1，1）模型是灰色系统最基本的模型，下面以污染物浓度问题说明GM（1，1）模型的建立及求解过程。

例10.1某污染源中某种污染物质量浓度测量值如表10.1，试建立GM（1，1）模型

表10.1某污染物质量浓度测量值（mg/L）

年份

2001

2002

2003

2004

2005

2006

3.936

4.575

4.968

5.063

5.968

5.507

解：第一步，设原始数据为

第二步，对原始数据进行累加生成，即

因此累加生成数据为

第三步，构造矩阵

第四步，计算。

先求，即

根据逆矩阵的求解方法，得

再求的值，即

进而求得的值为

计算GM1_1的程序如下

function10toliti01(X0)

[m,n]=size(X0);

X1=cumsum(X0);

X2=[];

fori=1:n-1

X2(i,:)=X1(i)+X1(i+1);

end

B=-0.5.*X2;

t=ones(n-1,1);

B=[B,t];

YN=X0(2:end);

P_t=YN./X1(1:(length(X0)-1))

A=inv(B.*B)*B.*YN.;

a=A(1)

u=A(2)

B

b1=B.*B

b2=inv(B.*B)

b3=B.*YN.

b4=u/a

b5=X1(1)-b4

b6=-a*b5

第五步，将的值代入微分方程的时间响应函数，

令，得

第六步，求导还原得

第七步，对上述模型进行精度检验。

常用的方法是回代检验，即分别用模型求出各时刻值，然后求相对误差。

先利用时间响应函数模型求各时刻值（），并计算相对误差，结果如表10.2所示.

表10.2精度检验实测值、残差值表

GM计算值

实测值

残差

相对残差

8.6059

13.5344

18.7359

24.2254

30.0190

8.5110

13.479018.542024.510030.0170

-0.0949-0.0554-0.19390.2846

-0.0020

-0.0112

-0.0041

-0.0105

0.0116

-0.0001

再利用时间响应函数模型求各时刻值（），并计算相对误差，结果如表10.3所示