灰色预测理论详解实用教案

1、目录(ml) 灰色系统理论简介以及发展 灰色预测理论 灰色预测简介 灰色预测类型 灰色预测模型 灰色预测检验(jinyn) 案例以及软件实现第1页/共59页第一页，共60页。灰色系统理论(lln)简介 灰色系统理论是由我国著名学者邓聚龙教授于1982年创立的“以部分信息已知，部分信息未知的小样本、贫信息”不确定系统为研究对象的一门(y mn)系统科学新学科，具有原创性的科学意义，是我国对系统科学的新贡献，目前已受到国内外学术界的广泛重视，并在农业科学、经济管理、环境科学、医药卫生、矿业工程、教育科学、水利水电、图像信息、生命科学、控制科学、航空航天等众多领域中得到了广泛的应用，解决了许多过去难

2、以解决的实际问题。 第2页/共59页第二页，共60页。灰色系统理论(lln)的内容 灰色系统理论经过20多年(du nin)的发展，已基本建立起一门新兴学科的结构体系。其主要内容包括以灰色朦胧集为基础的理论体系，以灰色关联空间为依托的分析体系，以灰色序列生成为基础的方法体系，以灰色模型(GM )为核心的模型体系，以系统分析、评估、建模、预测、决策、控制、优化为主体的技术体系。第3页/共59页第三页，共60页。 灰色朦胧集、灰色代数系统、灰色方程、灰色矩阵等是灰色系灰色朦胧集、灰色代数系统、灰色方程、灰色矩阵等是灰色系统理论的基础，从学科体系自身的优美、完善出发，这里有许多统理论的基础，从学科体

3、系自身的优美、完善出发，这里有许多问题值得进一步研究。灰色系统分析除灰色关联问题值得进一步研究。灰色系统分析除灰色关联(gunlin)(gunlin)分分析外，还包括灰色聚类和灰色统计评估等内容。灰色序列生成通析外，还包括灰色聚类和灰色统计评估等内容。灰色序列生成通过序列算子的作用来实现，序列算子主要包括缓冲算子过序列算子的作用来实现，序列算子主要包括缓冲算子( (弱化算弱化算子、强化算子子、强化算子) )、均值生成算子、级比生成算子、累加生成算子、均值生成算子、级比生成算子、累加生成算子和累减生成算子等。灰色模型按照五步建模思想构建，通过灰色和累减生成算子等。灰色模型按照五步建模思想构建，通

4、过灰色生成或序列算子的作用弱化随机性，挖掘潜在规律，经过灰色差生成或序列算子的作用弱化随机性，挖掘潜在规律，经过灰色差分方程与灰色微分方程之间的互换实现了利用离散的数据序列建分方程与灰色微分方程之间的互换实现了利用离散的数据序列建立连续的动态微分方程的新飞跃。灰色预测是基于立连续的动态微分方程的新飞跃。灰色预测是基于GM GM 模型作出模型作出的定量预测，按照其功能和特征可分为数列预测、区间预测、灾的定量预测，按照其功能和特征可分为数列预测、区间预测、灾变预测、季节灾变预测、波形预测和系统预测等几种类型。灰色变预测、季节灾变预测、波形预测和系统预测等几种类型。灰色决策包括灰靶决策、灰色关联决策

5、包括灰靶决策、灰色关联(gunlin)(gunlin)决策、灰色统计、聚决策、灰色统计、聚类决策、灰色局势决策和灰色层次决策等。灰色控制的主要内容类决策、灰色局势决策和灰色层次决策等。灰色控制的主要内容包括本征性灰色系统的控制问题和以灰色系统方法为基础构成的包括本征性灰色系统的控制问题和以灰色系统方法为基础构成的控制，如灰色关联控制，如灰色关联(gunlin)(gunlin)控制和控制和GM (1GM (1，1)1)预测控制等。预测控制等。灰色优化技术包括灰色线性规划、灰色非线性规划、灰色整数规灰色优化技术包括灰色线性规划、灰色非线性规划、灰色整数规划和灰色动态规划等。划和灰色动态规划等。灰色

6、系统(xtng)理论的内容第4页/共59页第四页，共60页。基础知识1.灰色系统是指“部分信息已知,部分信息未知”的“小样本”,“贫信息”的不确定性系统。 信息不完全包含： 1、系统因素不完全明确； 2、因素关系不完全清楚； 3、系统结构不完全知道； 4、系统作用原理不完全明了。2.白色系统、灰色系统、黑色系统 白色系统：一个系统的内部特征是完全已知的，即系统的信息是完全 充分的。如：存取款系统，存款金额明确，利息固定则最终取款金额就已知。 灰色系统：一个系统的内部特征是不完全已知的系统。人体是一个系统，人的身高、体温、血压等都是已知的，可是，人体内部在结构及部位功能(gngnng)上还有许多

7、问题尚未可知。 黑色系统：一个系统的内部信息对外界来说是一无所知的，只能通过它与外界的联系来加以观测研究。如：观测到的星体。第5页/共59页第五页，共60页。灰色系统分析法、数理统计(sh l tn(sh l tn j)j)法及模糊法对比灰色系统数理统计方法模糊法内涵小样本不确定大样本不确定界限不确定依据信息覆盖概率统计隶属度函数手段生成统计边界取值特点少数据少数据多数据经验（数据）要求允许任意分布允许任意分布要求典型分布函数目标现实规律历史统计规律认知表达信息准则最少信息无限信息经验信息第6页/共59页第六页，共60页。灰色灰色(hus)(hus)预测预测 灰色预测是对既含有已知信息又含有不

8、确定信息的系统进行预测，就是对在一定范围内变化的、与时间有关的灰色过程进行预测。 通过对原始数据的生成处理来和灰色模型的建立，挖掘、发现、掌握寻求系统变动的规律。生成数据序列有较强的规律性，可以用它来建立相应的微分方程模型，从而预测事物未来(wili)的发展趋势和未来(wili)状态，对系统的未来(wili)状态做出科学的定量分析。第7页/共59页第七页，共60页。常用(chn yn)的灰色预测有五种（1）数列预测，即用观察到的反映预测对象特征的时间序列来构造灰色预测模型，预测未来某一时刻的特征量，或达到某一特征量的时间。（2）灾变与异常值预测，即通过灰色模型预测异常值出现的时刻，预测异常值什

9、么时候出现在特定(tdng)时区内。（3）季节灾变与异常值预测，即通过灰色模型预测灾变值发生在一年内某个特定(tdng)的时区或季节的灾变预测。（4）拓扑预测，将原始数据作曲线，在曲线上按定值寻找该定值发生的所有时点，并以该定值为框架构成时点数列，然后建立模型预测该定值所发生的时点。（5）系统预测，通过对系统行为特征指标建立一组相互关联的灰色预测模型，预测系统中众多变量间的相互协调关系的变化。第8页/共59页第八页，共60页。预测(yc)模型单序列(xli)灰色预测模型GM（1，1）模型DGM（1，1）模型GM（1，N）模型Verhulst模型区间灰数预测模型基于几何坐标法的区间灰数预测模型I

10、GPM-G(1,1)基于信息分解法的区间灰数预测模型IGPM-P(1,1)基于灰色属性法的区间灰数预测模型IGPM-D(1,1)第9页/共59页第九页，共60页。单序列(xli)灰色预测模型灰色系统理论认为：系统的行为现象尽管朦胧，数据尽管复杂，但它必然是有序的，都存在着某种内在规律。不过这些规律被纷繁复杂的现象所掩盖,人们很难直接从原始数据中找到某种内在的规律.灰色生成：建立灰色模型之前，需要对原始时间序列按照某种要求进行预处理，得到有规律的时间序列数据生成列。即对原始数据的生成就是企图从杂乱无章的现象中去发现(fxin)内在规律.常用的灰色系统生成方式有: 累加生成,累减生成,均值生成,级

11、比生成等,下面对这几种生成做简单介绍： 第10页/共59页第十页，共60页。1.累加生成 通过数列间时刻各数据的依个累加以得到新的数据与数列，累加所得的新数列叫做累加生成数列。 具体地，记原始数列为 X(0)=(x(0)(1), x(0)(2),x(1)(n) 累加生成序列 X(i)=(x(i)(1), x(i)(2),x(i)(n) 一次累加生成关系 累加生成的作用：通过累加生成可以看出灰量积累过程的发展态势,使离乱的原始数据中蕴含的积分( jfn)特性或规律加以显化。 2.累减生成 对数列求相邻两数值的差，是累加生成的逆运算。 记原始序列为 X(1)=(x(1)(1), x(1)(2),)

12、,x(1)(n) 一次累减生成序列为 X(0)=(x(0)(1), x(0)(2),x(0)(n) 其中，x(0)(k)=x(1)(k)-x(1)(k-1) 累减生成的作用 累减生成可将累加生成还原为非生成数列,在建模方 程用来获得增量信息。(1)(0)(1)(0)0( )( )(1)( )kixkxixkxk第11页/共59页第十一页，共60页。3.3.均值生成 设原始序列为 X(0 X(0）=x(0)(1), x(0)(2),x(0)(n) =x(0)(1), x(0)(2),x(0)(n) ，x(0)(k-1)x(0)(k-1)与x(0)(k)x(0)(k)为数列X(0)X(0)的一对（

13、紧）邻值，则称x(0)(k-1)x(0)(k-1)为值，x(0)(k)x(0)(k)称为后值。 对于常数(0,1)(0,1) 则称 z(0)(k)=x(0)(k) + (1-)x(0)(k z(0)(k)=x(0)(k) + (1-)x(0)(k 1)1) 为由数列x(0)x(0)的邻值在生成系数（权）下的邻值生成数。 特别(tbi)(tbi)地，当生成系数为0.50.5时，则称z(0)(k) = 0.5x(0)(k) + 0.5x(0)(k z(0)(k) = 0.5x(0)(k) + 0.5x(0)(k 1)1) 为（紧）邻均值生成数，即等权邻值生成数。4.4.级比生成 级比生成是一种常用

14、的填补序列端点空穴的方法。对数列端点值的生成，我们无法采用均值生成填补空缺，只能采用级比生成。 级别生成在建模中可以获得较好的灰指数律。级比生成是级比(k)(k)与光滑比(k)(k)生成的总称。 设原始序列为 X(0 X(0）=x(0)(1), x(0)(2),x(0)(n)=x(0)(1), x(0)(2),x(0)(n)，称 (k) (k)为级比，(k)(k)为光滑比，其表达式为 (k)=x(0)(k)/x(0)(k (k)=x(0)(k)/x(0)(k 1)1) (k)=x(0)(k)/x(1)(k-1) (k)=x(0)(k)/x(1)(k-1) 第12页/共59页第十二页，共60页。

15、GM（1.1）模型(mxng)l 模型符号含义模型符号含义l GM(1，1) Grey Model(1阶方程，阶方程，1个变量个变量)l GM(1，1)建模过程建模过程l 令令X（0）为）为GM(1，1)为原始建模序列：为原始建模序列：l X(0)=(x(0)(1), x(0)(2),x(0)(n)l 其中其中x(0)(k)0,k=1,2,.,nl X（1） 为为X（0）累加生成序列）累加生成序列l X(1)=(x(1)(1), x(1)(2),x(1)(n)l 令令Z（1）为）为X（1）的紧邻均值生成序列）的紧邻均值生成序列l Z（1）=（ z(1)(1 ), z(1)(2), z(1)(k

16、 )）l z(1)(k)=o.5 x(1)(k)+0.5 x(1)(k-1)l GM(1，1)的灰微分方程模型为的灰微分方程模型为l x(0)(k+1)+a z(1)(k)=bl 式中称为式中称为-a发展系数，发展系数，b为灰色为灰色(hus)作用量作用量(1)(0)(1)(0)0( )( )(1)( )kixkxixkxk第13页/共59页第十三页，共60页。v准备知识准备知识v一阶微分方程模型一阶微分方程模型 dx/dt+ax=bv导数的定义导数的定义 v当当t很小并取很小的单位很小并取很小的单位1时时vx(t+1)-x(t)=x/tv则离散形式可写为则离散形式可写为vx/t=x(1)(k

17、+1)-x(1)(k)=x(0)(k+1)v由由dx/dtx/tx(1)(k+1)-x(1)(k)，在，在x(1)(k)，x(1)(k+1)范围内，由于范围内，由于v很短时间内背景值（即很短时间内背景值（即x值）不会发生突变，则取均值值）不会发生突变，则取均值v z(1)(k+1)=o.5 x(1)(k)+0.5 x(1)(k+1)作为作为(zuwi)x的值。的值。v则得到灰微分方程为则得到灰微分方程为vx(0)(k+1)+a z(1)(k)=bv则可得矩阵方程则可得矩阵方程vx(0)(k+1)=-a z(1)(k)+bvYn=B&0()( )limtdxx ttx tdtt 第14页

18、/共59页第十四页，共60页。 设&为待估参数向量，即&=（a,b）T，则灰微分方程(wi fn fn chn)的最小二乘估计参数列满足&= 其中， B= Y n= 则称 为灰色微分方程(wi fn fn chn)x(0)(k)+az(1)(k)=b的白化方程，也叫影子方程。 将上面所求参数代入白化方程，求得其解为 还原到原始数据 (1)(1)(1)(2)1(3)1.( )1zzzn(0)(0)(0)(2)(3).( )xxxnYBBBTT1(1 )(1 )d xa xbd tabeabxkak) 1 () 1(x)1(1)() 1() 1(x) 1 (0)kxkxk第

19、15页/共59页第十五页，共60页。DGM（1，1）模型(mxng)第16页/共59页第十六页，共60页。1 11 11 11）1）- -(n(nX X(1)(1)X X(1)(1)X XB B(1)(1)(1)(1)(1)(1)第17页/共59页第十七页，共60页。第18页/共59页第十八页，共60页。DGM(1DGM(1，1)1)模型和GM(1GM(1，1)1)是完全等价的。DGM(1DGM(1，1)1)模型全面符合灰色预测模型的建模机理. .是一种新的灰色预测模型. .或者说是灰色预测模型的一种新形式。其中DGM(1DGM(1，1)1)模型更能够精确模拟齐次指数序列。对于非指数增长序列和

20、震荡(zhndng)(zhndng)序列，应选择微分，差分混合形态的GM(1GM(1，1)1)；对于接近齐次指数序列的非指数增长序列和震荡(zhndng)(zhndng)序列，应优先选择DGM(1DGM(1，1)1)模型。DGMDGM模型可以看做是GMGM模型的精确形式，当GMGM模型中的a a取值很小时，二者可替换。GM(1,1)和DGM(1,1)的关系(gun x)第19页/共59页第十九页，共60页。GM（1，N）模型(mxng)如果考虑的系统由若干个相互影响的因素组成，设 为系统特征数据(shj)序列，而)(),.,2(),1 ()0(1)0(1)0(1)0(1nxxxX)(),.,2

21、(),1 ()0(2)0(2)0(2)0(2nxxxX. . )(),.,2(),1 ()0()0()0()0(nxxxXNNNN为相关因素序列。 为 的1-AGO序列 ， 为 的紧邻(jnln)生成序列，则称 )1(iX)0(iXNi,.,2 , 1)1(1Z)1(1XNiiikxbkazkx2)1()1(1)0(1)()()(为GM(1,N)灰色微分方程。第20页/共59页第二十页，共60页。定义(dngy)TNbbaa.2为GM(1,N)灰色微分方程(wi fn fn chn)的参数列，根据最小二乘法可以得出：YBBBaTT1式中)(.)()(.)3(.)3()3()2(.)2()2()

22、1()1(2)1(1)1()1(2)1(1)1()1(2)1(1nxnxnzxxzxxzBNNNTnxxxY)(.)3()2()0(1)0(1)0(1)1 ()1 (33)1 (22)1 (1)1 (1.NNxbxbxbaxdtdx为GM(1,N)灰色(hus)微分方程（2）的白化方程，也称影子方程。第21页/共59页第二十一页，共60页。于是(ysh)，我们有1）白化方程（3）的解为NiNiatiiiiatNiNiiiatiiatdtetxbxbtxedtxbxdtetxbetx22)1()1()1(122)1()1(1)1()1(1)()0()0()0()0()()(2）当 变化幅度很小时

23、(xiosh)，可视 为灰常量，这样，GM(1,N)灰色微分方程（3）的近似时间响应式为),.,2 , 1()1(NiXiNiiikxb2)1()(NiiiakNiiiikxbaekxbaxkx2) 1 (2)() 1 (1) 1 (1) 1(1) 1(1) 0() 1(其中(qzhng) 取为 。)0()1(1x) 1 ()0(1x3）累减还原式为)() 1() 1()1(1)1(1)0(1kxkxkx (5) (4) 第22页/共59页第二十二页，共60页。灰色(hus)Verhulst(hus)Verhulst模型灰色系统理论自1982年诞生以来，它被广泛应用于农业生产、经济、管理和工程

24、技术等领域。灰色预测理论是灰色系统理论的主要内容之一，而GM（1，1）模型又是灰色预测理论中最核心的模型。而对于具有s型序列，则不适宜用GM（1，1）模型预测，适合用灰色Verhulst模型等进行预测。灰色Verhulst模型避免了传统的Verhulst模型建模的大样本要求，因此，该模型得到了较广泛的使用。Verhulst模型主要用来描述具有饱和状态的过程,即S形过程,常用于人口预测,生物生长,繁殖预测和产品(chnpn)经济寿命预测等.由Verhulst方程的解可以看出,当t时,若a0,则x(1)(t) 0;若at, x(1)(k+1) 与x(1)(k) 充分接近,此时x(0)(k+1) 0

25、 ,系统趋于死亡.第23页/共59页第二十三页，共60页。设X(0)为原始数据序列(xli), X(1)为X(0)的1-AGO序列(xli), Z(1)为X(1)的紧邻均值生成序列(xli),则称为GM(1,1)幂模型. 称为GM(1,1)幂模型的白化方程.GM(1,1)幂模型之白化方程的解为)()1()1()1(xbaxdtdx)()1()1()0(ZbaZX11)1 ()1 () 1 ()1()(CdtbeetxataataVerhulst模型(mxng)第24页/共59页第二十四页，共60页。 设 ，有则GM(1,1)幂模型(mxng)参数列 的最小二乘估计为)1()1()0(,ZXX

26、)() 3 () 2()0()0()0(nxxxY )()()3 () 3 ()2 () 2 () 1 () 1 () 1 () 1 () 1 () 1 (nznzzzzzBTbaa, YBBBbaaTTT1)(,当a=2时,称为灰色Verhulst模型(mxng).称 为灰色Verhulst模型(mxng)的白化过程.2)1()1()0()(ZbaZX2) 1 () 1 () 1 ()(xbaxdtdx第25页/共59页第二十五页，共60页。 Verhulst白化方程的解为灰色Verhulst模型(mxng)的时间响应式为atatatatatebxabxaxebxaeaxeabxetx)0(

27、)0()0()1)(0()0()1 ()0(11)()1()1()1()1()1()1()1(akebxabxaxkx)0() 0() 0() 1() 1 () 1 () 1 () 1 (第26页/共59页第二十六页，共60页。区间(q jin)灰数预测模型 区间灰数：定义 既有下界 ,又有上界 的灰数称为区间灰数，记为 ，其中 。 白化权函数：是用来描述一个区间灰数 在其取值范围 内对不同数值(shz)的“偏爱”程度的函数，称为 的白化权函数。 区间灰数的“核”和“灰度”是区间灰数的两个重要属性，是研究区间灰数代数运算法则以及建立区间灰数预测模型误差检验方法的基础。区间灰数的“核”，是在充分

28、考虑已知信息的条件下，最有可能代表区间灰数“白化值”的实数；区间灰数的“灰度”，反映了人们对灰色系统认识的不确定程度。我们通常的实数则是灰度为零且核为本身的特殊区间灰数。kakb kkkbat,kakb kkkbat,kkb,a kt第27页/共59页第二十七页，共60页。区间灰数序列(xli)的白化方法 区间灰数序列： 其中 。 根据区间灰数序列的几何特点，信息特征以及灰色属性(shxng)，通常将区间灰数序列的白化方法分为几何坐标法，信息分解法以及灰色属性(shxng)法三类，将区间灰数序列转化为相关的实数序列。 ,)(),.,(),(21ntttXnkbatkkk,.,2 , 1,)(,

29、第28页/共59页第二十八页，共60页。转换方法结果转换公式信息等价性几何坐标法面积序列S=（s(1),s(2),.,s(n-1)） s(p)=等价坐标序列W=（w(1),w(2),.,w(n-1)）w(p)=信息分解法实部序列R=等价灰部序列H=等价灰色属性法核序列K=灰度序列G= 2bb11ppppaa 4bb11ppppaanaa,.,2, 1anhhh,.,21, nt,.,t,t21)(g),.,(g),(g21ntttkkkhat)(kkkab h2)(kkkbatababtkkk)(g第29页/共59页第二十九页，共60页。附图-基于几何(j h)坐标法的区间灰数序列白化方法注释

30、基于灰色(hus)属性法的区间灰数序列白化方法其中所有区间灰数产生的背景或论域为 ，也就是区间灰数的最大边界(binji)，为已知。且 在其范围内的取值可能性均等，即区间灰数 的白化权函数为矩形。ba,第30页/共59页第三十页，共60页。 基于几何坐标法的区间灰数预测模型IGPM-G(1,1) 即分别构建基于面积序列和坐标序列的DGM(1,1)模型，并在此基础上通过几何坐标法的反向推导，实现对区间灰数上界(shngji)及下届的模拟，进而实现区间灰数预测模型的构建。 具体模型建立参照前边的DGM(1,1)模型，可以得到面积序列和坐标序列的的预测模型， 从而得到11211) 1 (1sksk

31、112111w) 1(wkk 221121131kk221121131kk1114ab1112ababCabCkkkskkks第31页/共59页第三十一页，共60页。 解方程组，可得到区间灰数 上界和下界的模拟(mn)及预测公式， 其中)(kt22113121131221131211311211b1211abFFaFFkkkskkskkkkskkwk 212111w11s11,1112,12wCwsCCFCFsws第32页/共59页第三十二页，共60页。 基于信息分解法的区间灰数预测模型(mxng)IGPM-P(1,1) 即分别构建基于实部（也称“白部”）序列和灰部序列的DGM(1,1)模型(

32、mxng)，并在此基础上通过推导区间灰数上界及下界的模拟表达式，实现区间灰数预测模型(mxng)的构建。 具体模型(mxng)建立参照前边的DGM(1,1)模型(mxng)，可以得到白部序列和序列的的预测模型(mxng)， 11211112111k112111k1k1k1k112111k112111k11b1aabh1h1akkkkkahaha可得，由第33页/共59页第三十三页，共60页。 基于灰色属性法的区间灰数预测模型IGPM-D(1,1) 即通过“核”序列为基础建立DGM(1,1)预测模型，实现对未来区间灰数“核”的预测；然后以“灰度不减公理”为理论依据，以“核”为中心拖展得区间灰数的

33、上界和下界，在不破坏区间灰数独立性和完整性的前提下，实现区间灰 数的模拟和预测。 核序列具体模型建立参照(cnzho)前边的DGM(1,1)模型，可以得到核序列的预测模型， 灰度不减原理：两个灰度不同的区间灰数进行和，差，积，商运算时，运算结果的灰度不小于灰数较大的区间灰数的灰度，即 1121111)(kktt第34页/共59页第三十四页，共60页。 可得 1111112112baab,.,2 , 1)(g)(g)(g),.(g),(gmax)(gkkkxxkkxknktabnxtttttt由则其中21)(21)(112111112111xxkkxxkkabtbabta第35页/共59页第三十

34、五页，共60页。单序列预测(yc)(yc)模型检验 灰色模型是灰色系统理论的重要组成部分，在各个领域的广泛应用说明：灰色理论具有很多独特的优越性。灰色模型应用的广泛，其精度成为了一个很受关注的问题。 灰色模型的误差a.a.抽样误差或试验误差它是原始数据记录或统计过程中出现的误差b.b.舍入误差或计算(j sun)(j sun)误差它是建模计算(j sun)(j sun)过程中由于运算数据的舍入引起的误差C.C.模型表达误差由于灰色系统以“近似地”表达为描述问题和就决问题的指导思想，所以原始数据和模拟值之间就必然存在微小的差距 。第36页/共59页第三十六页，共60页。误差分析 前两种误差带有一

35、定的客观性和必然性，是模型不可避免的误差；后一种误差可以通过不断的改进使模型的误差逐渐降低。 由于人们用精准性理论建立符合复杂事物发展规律的模型几乎是不可能的，因此建立背景值为灰的近似模型也就不可避免地出现微小误差，但是这种很小的误差并不影响人们利用灰模型研究事物发展的内在规律，并以此指导(zhdo)(zhdo)生产实践，只有当模型误差超出了一定的限制范围时，模型的实用性才会降低。第37页/共59页第三十七页，共60页。模模型型选选定定之之后后, ,一一定定要要经经过过检检验验才才能能判判定定其其是是否否合合理理, ,只只有有通通过过检检验验的的模模型型才才能能用用来来作作预预测测, ,灰灰色

36、色模模型型的的精精度度检检验验一一般般有有三三种种方方法法: :相相对对误误差差大大小小检检验验法法, ,关关联联度度检检验验法法和和后后验验差差检检验验法法. .下下面面对对这这三三种种方方法法做做个个简简单单介介绍绍. .(1)(1)(0)(1.1),GMXXX设设按按建建模模法法已已求求出出并并将将做做一一次次累累减减转转化化为为即即1.1.1 相对误差(xin du w ch)检验法(0)(0)(0)(0)(1),(2), ,( )(2 31)Xxxxn (1-1)第38页/共59页第三十八页，共60页。计计算算残残差差得得(0)(0) (1), (2), ( )(2 32)Eeee

37、nXX (0)(0), ( )( )( ),1,2,e kxkxkkn其其中中计计算算相相对对误误差差得得(0)( )( )100% ,1,2, ,(2 33)( )ekrel kknxk 计计算算平平均均相相对对误误差差得得11( ),(2 34)nkrelrel kn (1-2)第39页/共59页第三十九页，共60页。计计算算后后验验差差比比为为21/(2 36)C S S 计计算算小小误误差差概概率率为为 1( )0.6745(2 37)p P e keS 1212,.CpCCSSSSC指指标标 和和 是是后后验验差差检检验验的的两两个个重重要要指指标标. .指指标标 越越小小越越好好越

38、越小小表表示示 大大而而越越小小大大表表示示原原始始数数据据方方差差大大, ,即即原原始始数数据据离离散散程程度度大大. . 小小表表示示残残方方差差小小, ,即即残残差差离离散散程程度度小小. . 小小就就表表明明尽尽管管原原始始数数据据很很离离散散, ,而而模模型型所所得得计计算算值值与与实实际际值值之之差差并并不不太太离离散散. .第40页/共59页第四十页，共60页。1.1.2 后验差检验法( 0 )( 0 )2212(1 .1)(23 1),(23 2 ),G MXXESS 设设 按按建建 模模 法法 所所 求求 出出 的的如如所所 示示残残差差 如如所所 示示原原 始始 序序 列列

39、及及 残残 差差 序序 列列的的 方方 差差分分 别别 为为和和则则2(0)21122211( )1 ( )(2 35)nknkSxkxnSe ken (0)1111,( ),( )nnkkxxk ee knn其其中中第41页/共59页第四十一页，共60页。1,0.6745,.ppC p指指标标 越越大大越越好好越越大大 表表明明残残差差与与残残差差平平均均值值之之差差小小于于给给定定值值的的点点较较多多 即即拟拟合合值值( (或或预预测测值值) )分分布布比比较较均均匀匀. .按按两两个个指指标标 可可综综合合评评定定预预测测模模型型的的精精度度 模模型型的的精精度度由由后后验验差差和和小小

40、误误差差概概率率共共同同刻刻划划. .一一般般地地, ,将将模模型型的的精精度度分分为为四四级级, ,见见表表2 2- -1 1模型精度级模型精度级均方差比值均方差比值C小误差概率小误差概率p1级（好）级（好）C=0.350.95=p2级（合格）级（合格）0.35C=0.50.80=p0.953级（勉强）级（勉强）0.5C=0.650.70=p0.804级（不合格）级（不合格） 0.65CP0.70 ,MaxpC 模模型型的的精精度度级级别别的的级级别别于于是是的的级级别别表21 精度(jn d)检验等级参照表第42页/共59页第四十二页，共60页。1.1.3 关联度检验法关联度检验，即通过考

41、察模型值曲线和建模序列(xli)曲线的相似程度进行检验，是一种几何检验。关联度计算方法，先计算出 与 原始序列(xli)的关联系数，然后算出关联度，根据经验，关联度大于0.6便是满意的。 )()0(ix)()0(ix如果 在所有(suyu)关联度中最大，则第i种灰色建模方法为所建模型中最好的模型。),.,2 , 1(miri第43页/共59页第四十三页，共60页。 设， 则关联系数定义(dngy)为： nXXXkX0000,.,2,1 nXXXkX0000,.,2,1 kXtXkXkXkXkXkXkXk00000000maxmaxmaxmaxminmin)( kXkX00 0X 0X kXkX00minmin kXkX00maxmax式中:为第k个点与的绝对误差(ju du w ch)；为两级最小差； 为两级最大差； 第44页/共59页第四十四页，共60页。称为分辨率，01,一般取=0.5；对单位不一，初值不同的序列，在计算相关系数前应首先进行初始化，即将(jjing)该序列所有数据分别除