数学灰色系统理论.pptx

《数学灰色系统理论.pptx》由会员分享，可在线阅读，更多相关《数学灰色系统理论.pptx（38页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、什么是灰色系统？什么是灰色系统？从信息的完备性与模型的构建上看，工程技术等系统具有较充足的信息量，其发展变化规律明显，定量描述较方便，结构与参数较具体，人们称之为白色系统；一个系统的内部特性全部未知，则称之为黑色系统。对另一类系统诸如社会系统、农业系统、生态系统等，人们无法建立客观的物理原型，其作用原理亦不明确，内部因素难以辨识或之间关系隐蔽，人们很难准确了解这类系统的行为特征，因此对其定量描述难度较大，带来建立模型的困难。这类系统称之为灰色系统。灰色系统理论就是研究在信息大量缺乏或紊乱的情况下，如何对实际问题进行分析和解决。第1页/共38页 区别白色系统与灰色系统的重要标志是系统内各因素之间

2、是否具有确定关系。运动学中物体运动的速度、加速度与其所受到的外力有关，其关系可用牛顿定律以明确的定量来阐明，因此，物体的运动便是一个白色系统。当然，白、灰、黑是相对于一定的认识层次而言的，因而具有相对性。某人有一 天去他朋友家做客，发现当外面的汽车开过来时，他朋友家的狗就躲到屋角里瑟瑟发抖。他对此莫名其妙。但对他朋友来讲，狗的这种行为是可以理解的，因为他知道，狗在前不久曾被汽车撞伤过。显然，同样对于“狗的惧怕行为”，客人因不知内情而面临一个 黑箱，而主人则面临一个灰箱。第2页/共38页灰色系统的应用范畴大致分为以下几方面：（1 1）灰色关联分析；（2 2）灰色预测：人口预测；初霜预测；灾变预测

3、.等等；（3 3）灰色决策；（4 4）灰色预测控制。灰色系统理论的应用范畴？灰色系统理论的应用范畴？第3页/共38页 灰色系统的分析方法？灰色系统的分析方法？进行关联分进行关联分析析 客观世界中的事物往往现象复杂，因素繁多。我们往往需要对系统进行因素分析，这些因素中哪些对系统来讲是主要的，哪些是次要的，哪些需要发展，哪些需要抑制，哪些是潜在的，哪些是明显的。一般来讲，这些都是我们极为关心的问题。事实上，因素间关联性如何、关联程度如何量化等问题是系统分析的关键和起点。因素分析的基本方法过去主要采取回归分析等办法，但是这种方法需要大量数据作为基础，计算量大。而灰色系统理论采用的关联分析方法可以克服

4、这个弊端。灰色系统理论进行关联分析的两种方法：一 根据数据的几何关系分析法；二 利用关联公式分析法第4页/共38页生成数的生成方法生成数的生成方法生成方法生成方法应用相关应用相关公式公式一次累加时间一次累减时间均值生成时间+其他第5页/共38页问题引入问题引入某地区最近17年来年度平均降雨量数据（单位：mm）序列为：X=(390.6，412.0，320.0，559.2，380.8，542.4，553.0，310.0，561.0，300.0，632.0，540.0，406.2，313.8，576.0，586.6，318.5)如果将年平均降雨量低于320mm时认为旱灾发生，试根据上述数据预测下一次

5、旱灾发生在几年后？第6页/共38页如何建立灰色系统如何建立灰色系统 GMGM（1 1，1 1）模型？模型？第7页/共38页GM(1,1)模型G表示grey（灰色）M表示model（模型）GM（1，1）表示1阶的、1个变量的模型。第8页/共38页设 其中 则称 为GM(1,1)模型的基本形式。定义1.1第9页/共38页定理1.1 设有非负序列：为 的1-AGO（即一次累加）序列：其中 ；为 的紧邻均 值生成序列：其中 第10页/共38页若 为参数列，且 则GM(1,1)模型 的最小二乘估计参数列满足 第11页/共38页定义1.2设 为非负序列，为 的1-AGO（即一次累加）序列，为 的紧邻均值生

6、成序列，则称微分方程 为（，）模型（灰色方程）的 白化方程，也叫影子方程。第12页/共38页定理.2设 如定理1.1中所述，其中 ,则.白化方程 的解（也称时间响应函数）为.GM(1,1)模型 的时间响应序列为第13页/共38页.还原值第14页/共38页2.例题 设有原始数据序列 X0=(x0(1),x0(2),x0(3),x0(4),x0(5)=(2.874,3.278,3.337,3.390,3.679)试用GM(1,1)模型对X0进行模拟第15页/共38页第一步：对X0作 1-AGO，得 X1=(X1(1),X1(2),X1(3),X1(4),X1(5)=(2.874,6.152,9.4

7、89,12.897,16.558 )第二步：对X1作紧邻均值生成。令得Z1=(Z1(2),Z1(3),Z1(4),Z1(5)=(4.513,7.820,11.184,14.718 )可得B,Y第16页/共38页第三步：对参数列 进行最小二乘估计，得 =-0.0372 3.06536第四步：确定模型为时间响应式为第17页/共38页第六步：还原出 的模拟值，由得对比原数据 X0=(x0(1),x0(2),x0(3),x0(4),x0(5)=(2.874,3.278,3.337,3.390,3.679 )第五步：求X1的模拟值第18页/共38页如何利用灰色系统实现预测？第19页/共38页 灰色预测是

8、指利用 GM 模型对系统行为特征的发展变化规律进行估计预测，同时也可以对行为特征的异常情况发生的时刻进行估计计算，以及对在特定时区内发生事件的未来时间分布情况做出研究等等。这些工作实质上是将“随机过程”当作“灰色过程”，“随机变量”当作“灰变量”，并主要以灰色系统理论中的GM(1,1)模型来进行处理。灰色预测在工业、农业、商业等经济领域，以及环境、社会和军事等领域中都有广泛的应用。特别是依据目前已有的数据对未来的发展趋势做出预测分析。第20页/共38页灰色预测的步骤1.数据的检验与处理 首先，为了保证建模方法的可行性，需要对已知数据序列做必要的检验处理。2.建立模型 按照上节的方法建立灰色模型

9、GM(1,1),则可以得到预测值第21页/共38页3检验预测值4预测预报 由模型 GM(1,1)所得到的指定时区内的预测值，根据实际问题的需要，给出相应的预测预报。第22页/共38页定义 设原始数据序列 相应的预测模型模拟序列:残差序列:相对误差序列:第23页/共38页第24页/共38页第25页/共38页第26页/共38页第27页/共38页灰色灾害预测应用灰色灾害预测实质上是异常值预测，什么样的值算作异常值，往往是人们凭经验主观确定。灰色灾害预测的任务是给出下一个或几个异常值出现的时刻，以便人们提前准备，采取对策。第28页/共38页定义：设原始序列 X=(x(1),x(2),x(n)，给定上限

10、异常值（灾变值）a,称X的子序列 Xa(x(q(1),x(q(2),x(q(m)=x(q(i)=a,i=1,2,m 为上灾变序列。同理，可定义下灾变序列。二者统一称为灾变序列。第29页/共38页定义：设X为原始序列，称 Q0=(q(1),q(2),q(m)为灾变日期序列。灾变预测就是要通过对灾变日期序列的研究，寻找其规律性，预测以后若干次灾变发生的日期，灰色系统的灾变预测是通过对灾变日期序列建立GM(1,1)模型实现的。第30页/共38页例 某地区最近17年来的年度平均降雨量数据（单位：mm）序列为X=(390.6,412.0,320.0,559.2,380.8,542.4,553.0,310

11、.0,561.0,300.0,632.0,540.0,406.2,313.8,576.0,586.6,318.5)如果将年平均降雨量低于320mm时认为旱灾发生，试根据上述数据预测下一次旱灾发生在几年后？第31页/共38页解：取灾变值为a=320,得下限灾变序列为Xa=(x(3),x(8),x(10),x(14),x(17)=(320.0,310.0,300.0,313.8,318.5)与之对应的灾变日期序列为Q0=(q(1),q(2),q(3),q(4),q(5)=(3,8,10,14,17 )其1-AGO序列为 Q1=(3,11,21,35,52)的紧邻均值生成序列为 Z1=(7,16,2

12、8,43.5)第32页/共38页设q(k)+az1(k)=b,易知B,Y,由最小二乘法得a,b=-0.253661 6.258339故灾变日期序列的GM(1,1)序号响应式为即第33页/共38页由此可得Q0的模拟序列为由得残差序列为第34页/共38页再由相对误差序列由此可计算出平均相对误差为第35页/共38页平均相对精度为1-=97.81%，故可用进行预测，即从最近一次旱灾发生的日期算起，5年以后，可能发生旱灾。为了提高预测的可靠程度，可以取若干个不同的异常值，建立多个模型进行预测。第36页/共38页计算的 MATLAB 程序如下：clc,cleara=390.6,412,320,559.2,

13、380.8,542.4,553,310,561,300,632,540,406.2,313.8,576,587.6,318.5;t0=find(a=320);t1=cumsum(t0);n=length(t1);B=-0.5*(t1(1:end-1)+t1(2:end),ones(n-1,1);Y=t0(2:end);r=BYy=dsolve(Dy+a*y=b,y(0)=y0);y=subs(y,a,b,y0,r(1),r(2),t1(1);yuce1=subs(y,t,0:n+1)digits(6),y=vpa(y)%为提高预测精度，先计算预测值，再显示微分方程的解yuce=diff(yuce1);yuce=t0(1),yuce第37页/共38页感谢您的观看。第38页/共38页