数学建模 灰色预测方法.ppt

1灰色预测理论2GM(1,1)模型3GM(1,1)模型的改进,灰色预测法,4灰色预测实例,1灰色预测理论,一、灰色预测的概念,（1）灰色系统、白色系统和黑色系统,白色系统是指一个系统的内部特征是完全已知的，即系统的信息是完全充分的。,黑色系统是指一个系统的内部信息对外界来说是一无所知的，只能通过它与外界的联系来加以观测研究。,灰色系统内的一部分信息是已知的，另一部分信息是未知的，系统内各因素间有不确定的关系。,灰色预测法是一种对含有不确定因素的系统进行预测的方法。灰色预测是对既含有已知信息又含有不确定信息的系统进行预则，就是对在一定范围内变化的、与时间有关的灰色过程进行预测。,（2）灰色预测方法,灰色预测法用等时距观测到的反映预测对象特征的一系列数量值构造灰色预测模型，预测未来某一时刻的特征量，或达到某一特征量的时间。,（3）灰色预测数据的特点：,1）序列性：原始数据以时间序列的形式出现。,2）少数据性：原始数据序列可以少到只有4个数据。,（4）灰色预测的四种常见类型,灰色时间序列预测即用观察到的反映预测对象特征的时间序列来构造灰色预测模型，预测未来某一时刻的特征量，或达到某一特征量的时间。灾变预测即通过灰色模型预测异常值出现的时刻，预测异常值什么时候出现在特定时区内。,系统预测通过对系统行为特征指标建立一组相互关联的灰色预测模型，预测系统中众多变量间的相互协调关系的变化。拓扑预测（波形预测）将原始数据做曲线，在曲线上按定值寻找该定值发生的所有时点，并以该定值为框架构成时点数列，然后建立模型预测该定值所发生的时点。,二、灰色生成数列,对灰数的处理主要是利用数据处理方法去寻求数据间的内在规律，通过对已知数据列中的数据进行处理而产生新的数据列，以此来研究寻找数据的规律性，这种方法称为数据的生成。数据的生成方式有多种，常用的方法有累加生成、累减生成和加权累加生成等。,（1）累加生成设原始数列为,，令,称为数列,的-次累加生成数列。类似地有,累加的规则:,将原始序列的第一个数据作为生成列的第一个数据，将原始序列的第二个数据加到原始序列的第一个数据上，其和作为生成列的第二个数据，将原始序列的第三个数据加到生成列的第二个数据上，其和作为生成列的第三个数据，按此规则进行下去，便可得到生成列。,对非负数据，累加次数越多则随机性弱化越多，累加次数足够大后，可认为时间序列已由随机序列变为非随机序列。一般随机序列的多次累加序列，大多可用指数曲线逼近。,累加举例：设原始时间序列为,一次累加生成列为,的曲线是摆动的，起伏变化幅度较大，而已呈现明显的增长规律性。,（2）累减生成,将原始序列前后两个数据相减得到累减生成序列,累减是累加的逆运算，累减可将累加生成列还原为非生成列，在建模中获得增量信息。,一次累减的公式为：,一般地，对于r次累加生成数列,则称,（3）均值生成,设原始数列,的邻值生成数（或生成值）。,为紧邻均值生成数，即等权邻值生成数。类似地，可以定义非紧邻值生成数,由,而得的数列,称为紧邻均值生成数列。,2GM（1，1）模型,灰色模型是利用离散随机数经过生成变为随机性被显著削弱而且较有规律的生成数，建立起的微分方程形式的模型，这样便于对其变化过程进行研究和描述。,灰色预测模型称为GM模型，G为grey的第一个字母，M为model的第一个字母。,GM（1，1）表示一阶的，一个变量的微分方程型预测模型。GM（1，1）是一阶单序列的线性动态模型，主要用于时间序列预测。,一、GM（1，1）模型概述,有：,其中：,则GM(1,1)的灰微分方程模型为：,灰导数,发展系数,灰作用量,GM(1,1)的白化型：,控制灰数。,为灰导数，对应于,为白化背景值，对应于,则灰微分方程对应的白化方程为：,灰方程也可改写为：,按最小二乘法求解，有：,式中：,预测模型白化响应式（解）为：,注意：GM(1,1)白化型不是从定义推导出来的，是一种“借用”或“白化默认”，所以，一切从白化推导出来的结果，只在不与定义型有矛盾时才成立，否则无效。,也可由GM(1,1)模型推导出另一表达式内涵型表达式：,灰色预测的事前检验,般用序列的光滑比对作准光滑性检验；用累加序列的级比对作准指数规律性检验来判断满足建模条件,光滑比定义：,（1）若光滑比满足递减且,则称为准光滑序列。,级比定义：,（2）若级比满足：,则认为具有准指数规律。当（1）（2）都满足时可对建GM(1,1)模型。,若原始数据不适合建立GM(1,1)模型，则进行予处理。,注：GM(1，1)模型中发展系数a的取值范围,二、GM（1，1）的建模步骤,灰色预测检验一般有残差、关联度和后验差检验。,三、模型检验,（1）残差检验,按预测模型计算,并将,累减生成,然后计算原始序列,与,的绝对误差序列及相,对误差序列。,残差：,相对误差：,残差序列,称为残差合格模型。,（2）关联度检验,第三步计算最小差与最大差,最小差为：,最大差为：,第四步计算关联系数,分辨系数，一般取0.5,第五步计算关联度,样本个数。,(3)后验差检验,a.计算原始数列的均值,b.计算原始数列的方差,根据经验，当=0.5时，关联度大于0.6便满意了。,此外，也可计算与的绝对关联度【1】，若对于给定的，则称为关联度合格模型。,e.计算均方差比,f.计算小误差概率,注：对给定的称模型为均方差比合格模型。,注：对给定的称模型为小误差概率合格模型。,g.检验根据经验，对给定的一组取值，就确定了检验模型模拟精度的等级划分如下表。,通过以上检验，如果相对误差、关联度、均方差比值、小误差概率都在允许范围之内时，则可用所建模进行预测，否则应进行残差修正。,表预测精度等级划分,四、预测,GM(1,1)模型经以上检验合格后可用于预测，其预测公式为：,对于数列预测，要建立多个预测模型，得到多组预测值，然后进行分析，从中确定出一个合适的预测模型，以取定一组合适的预测值。,对于一组数列，要建立多个预测模型，是通过对原始数列进行不同的取舍，形成新的数列，即对数列中的数据用不同的组合方式和取舍方式派生出新的数列，对原始数列和派生出来的新数据都建立预测模型，这样就对一个数列建立了多个预测模型。,例如有下述原始数列：,再加上原始数列，就可建立五个GM（1，1）模型。,3GM(1,1)模型的改进,一、用残差模型进行修正,若用原始时间序列,检验不合格或精度不理想时，则可用GM（1，1）残差模型进行修正以提高原GM（1，1）模型的预测精度，从而达到改进目的。,建立的GM(1,1)模型,（1）GM（1，1）残差模型,的数据序号；,注意：,符号一致；,必须满足,将上述的残差尾段仍记为,对上式求导数得残差尾段的模拟序列,其中,二、对GM(1,1)模型的其它改进方法,（1）新信息GM（1，1）模型不断地补充新出现的信息，即在预测下一时刻的值时，将最新的信息加入。此模型随着时间推移，序列长度会越来越长；（2）新陈代谢GM（1，1）模型即新信息出现后，将老信息去掉，加入新信息，保持序列长度不变；（3）GM（1，1）模型群法用原始时间序列数据建立多个GM（1，1）模型，给出预测值的区间；,一、数列的预测实例,原始数据,(1)求原始序列的一届累加生成,4灰色预测实例,(2)对作准光滑性检验,(3)对作准指数规律性检验,(4)作的紧邻均值生成序列并且确定B,Y,(5)按最小二乘法确定a,b的估计值,(6)确定模型,其时间响应式,并得的模拟值,(7)检验误差,残差序列,相对误差序列,相对误差检验,(8)预测,二、灾变灰预测方法,由异常值构成的序列,指序列中有异常值，是异常值可能在未来某时区发生的预测,步骤一：原始序列、阈值,给出原始序列，指定阈值（正常值与异常值的界限）,步骤二：构造异常值序列,步骤三：时分布序列,由异常值序列得时分布序列,步骤四：对时分布序列作GM(1,1)建模,步骤五：预测,三、季节灾变灰预测,对发生在每年特定时区的事件作预测,如春雨；根据往年数据预测某时间出现春雨的年份。初霜、棉铃虫虫灾发生的日期、长江洪水、小麦干热风等，都是季节性事件,对季节灾变预测方法及实例参考文献【