数学建模灰色预测模型.pptx

灰色预测模型（Gray Forecast Model）是通过少量的、不完全的信息，建立数学模型并做出预测的一种预测方法.当我们应用运筹学的思想方法解决实际问题，制定发展战略和政策、进行重大问题的决策时，都必须对未来进行科学的预测.预测是根据客观事物的过去和现在的发展规律，借助于科学的方法对其未来的发展趋势和状况进行描述和分析，并形成科学的假设和判断.第1页/共82页灰色系统理论是研究解决灰色系统分析、建模、预测、决策和控制的理论.灰色预测是对灰色系统所做的预测.目前常用的一些预测方法（如回归分析等），需要较大的样本.若样本较小，常造成较大误差，使预测目标失效.灰色预测模型所需建模信息少，运算方便，建模精度高，在各种预测领域都有着广泛的应用，是处理小样本预测问题的有效工具.第2页/共82页1 灰色系统的定义和特点灰色系统的模型Sars 疫情 4 销售额预测5 城市道路交通事故次数的灰色预测6 城市火灾发生次数的灰色预测7灾变与异常值预测第3页/共82页1 灰色系统的定义和特点灰色系统的定义和特点第4页/共82页灰色系统的定义和特点 灰色系统理论是由华中理工大学邓聚龙教授于19821982年提出并加以发展的。二十几年来，引起了不少国内外学者的关注，得到了长足的发展。目前，在我国已经成为社会、经济、科学技术在等诸多领域进行预测、决策、评估、规划控制、系统分析与建模的重要方法之一。特别是它对时间序列短、统计数据少、信息不完全系统的分析与建模，具有独特的功效，因此得到了广泛的应用.在这里我们将简要地介绍灰色建模与预测的方法.第5页/共82页一、灰色系统的定义和特点1.1.灰色系统的定义 灰色系统是黑箱概念的一种推广。我们把既含有已知信息又含有未知信息的系统称为灰色系统.作为两个极端，我们将称信息完全未确定的系统为黑色系统；称信息完全确定的系统为白色系统.区别白色系统与黑色系统的重要标志是系统各因素之间是否具有确定的关系。第6页/共82页1灰色系统的定义和特点2.灰色系统的特点（1）用灰色数学处理不确定量，使之量化.（2）充分利用已知信息寻求系统的运动规律.（3）灰色系统理论能处理贫信息系统.第7页/共82页1灰色系统的定义和特点常用的灰色预测有五种：（1）数列预测，即用观察到的反映预测对象特征的时间序列来构造灰色预测模型，预测未来某一时刻的特征量，或达到某一特征量的时间。（2）灾变与异常值预测，即通过灰色模型预测异常值出现的时刻，预测异常值什么时候出现在特定时区内。（3）季节灾变与异常值预测，即通过灰色模型预测灾变值发生在一年内某个特定的时区或季节的灾变预测。（4）拓扑预测，将原始数据作曲线，在曲线上按定值寻找该定值发生的所有时点，并以该定值为框架构成时点数列，然后建立模型预测该定值所发生的时点。（5）系统预测.通过对系统行为特征指标建立一组相互关联的灰色预测模型，预测系统中众多变量间的相互协调关系的变化。第8页/共82页2 灰色系统的模型灰色系统的模型第9页/共82页 在灰色系统理论中，把一切随机变量都看作灰色数，即使在指定范围内变化的所有白色数的全体，对灰数处理主要是利用数据处理的方法去寻求数据间的内在规律，通过对已知数据列中的数据进行处理而产生新的数据列，以此来研究寻求数据的规律性，这种方法称为数据的生成。常用的方法有：累加生成累减生成均值生成第10页/共82页1）累加生成 把数列各时刻数据依次累加的过程称为累加生成过程，记为AGO，有累加生成过程所得到的新数列称为累加生成数列。设原始数列为 ，令 则称 为数列 的1次累加生成，数列 称为数列 的1次累加生成数列。类似有第11页/共82页称之为 的r次累加生成，记 称之为 的r次累加生成数列第12页/共82页2）累减生成 对于原始数据列依次做前后相邻的两个数据相减的运算过程称为累减生成过程，记为IAGO，设原始数列为 ，令则称 为数列 的1次累减生成 一般地，对于r次累加生成数列则称为数列 的r次累减生成2）累减生成第13页/共82页3）均值生成设原始数列则称 与 为数列 的邻值，为后邻值，为前邻值。对于常数 ，则称为有数列 的邻值在生成系数（权）下的邻值生成数特别地，当生成系数 时则称为邻均值生成数，即等权邻值生成数第14页/共82页2 灰色系统的模型 通过下面的数据分析、处理过程，我们将了解到，有了一个时间数据序列后，如何建立一个基于模型的灰色预测。1.1.数据的预处理 首先我们从一个简单例子来考察问题.【例7.17.1】设原始数据序列第15页/共82页7.2 灰色系统的模型对数据累加 于是得到一个新数据序列第16页/共82页7.2 灰色系统的模型 归纳上面的式子可写为归纳上面的式子可写为 称此式所表示的数据列为原始数据列的一次累加称此式所表示的数据列为原始数据列的一次累加生成，简称为生成，简称为一次累加生成一次累加生成.显然有显然有 将上述例子中的 分别做成图、图7.2.可见图上的曲线有明显的摆动，图呈现逐渐递增的形式，说明原始数据的起伏已显著弱化.可以设想用一条指数曲线乃至一条直线来逼近累加生成数列 第17页/共82页7.2 灰色系统的模型图7.2 图为了把累加数据列还原为原始数列，需进行后减运算或称相减生成，它是指后前两个数据之差，如上例中第18页/共82页7.2 灰色系统的模型归纳上面的式子得到如下结果：一次后减其中第19页/共82页第20页/共82页第21页/共82页第22页/共82页第23页/共82页第24页/共82页第25页/共82页第26页/共82页第27页/共82页白化定义第28页/共82页第29页/共82页第30页/共82页第31页/共82页7.2 灰色系统的模型3.3.精度检验 (1)(1)残差检验：分别计算第32页/共82页7.2 灰色系统的模型（3 3）预测精度等级对照表，见表7.1.7.1.第33页/共82页7.2 灰色系统的模型注：由于模型是基于一阶常微分方程建立的，故称为一阶一元灰色模型，记为GM(1,1).GM(1,1).须指出的是，建模时先要作一次累加，因此要求原始数据均为非负数.否则，累加时会正负抵消，达不到使数据序列随时间递增的目的.如果实际问题的原始数据列出现负数，可对原始数据列进行“数据整体提升”处理.注意到一阶常微分方程是导出GM(1,1)GM(1,1)模型的桥梁，在我们应用GM(1,1)GM(1,1)模型于实际问题预测时，不必求解一阶常微分方程。第34页/共82页7.2 灰色系统的模型4.GM(1,1)4.GM(1,1)的建模步骤 综上所述，GM(1,1)GM(1,1)的建模步骤如下：第35页/共82页第36页/共82页第37页/共82页第38页/共82页第39页/共82页第40页/共82页第41页/共82页第42页/共82页第43页/共82页第44页/共82页第45页/共82页第46页/共82页第47页/共82页销售额预测销售额预测第48页/共82页7.3 销售额预测 随着生产的发展、消费的扩大，市场需求通常总是增加的，一个商店、一个地区的销售额常常呈增长趋势.因此，这些数据符合建立灰色预测模型的要求。【例7.2】表列出了某公司19992003年逐年的销 售额.试用建立预测模型，预测2004年的销售额，要求作精度检验。第49页/共82页7.3 销售额预测 表7.2 7.2 逐年销售额（百万元）年份年份1999199920002000200120012002200220032003 序号序号1 12 23 34 45 5 2.8742.8743.2783.2783.3373.3373.3903.3903.6793.679【例7.2】表列出了某公司19992003年逐年的销 售额.试用建立预测模型，预测2004年的销售额，要求作精度检验。第50页/共82页7.3 销售额预测 解（1 1）由原始数据列计算一次累加序列 ，结果见表7.3.7.3.表7.3 7.3 一次累加数据年份年份1999199920002000200120012002200220032003序号序号1 12 23 34 45 52.8742.8743.2783.2783.3373.3373.3903.3903.6793.6792.8742.8746.1526.1529.4899.48912.87912.87916.55816.558第51页/共82页7.3 销售额预测（2 2）建立矩阵：第52页/共82页7.3 销售额预测第53页/共82页7.3 销售额预测第54页/共82页7.3 销售额预测第55页/共82页7.3 销售额预测第56页/共82页7.3 销售额预测下面我们用用GMGM预测软件求解例7.2.7.2.参考附录B B（1 1）调用GMGM预测软件.见图7.3.7.3.图第57页/共82页 7.3 销售额预测（2 2）在“文件”菜单中打开“新建问题”，见到数据输入界面.见图7.4.7.4.第58页/共82页7.3 销售额预测（3 3）输入题目名称及元素个数后，点击“下一步”键，得到原始数据序列 的输入表格.见图7.5.第59页/共82页7.3 销售额预测（4 4）点击“运行”键，输出分析数据如下：题目:123:123原始数列(5(5个，预测结果如下：1dx/dt+ax=u1dx/dt+ax=u：，22时间响应方程：33残差E(k)E(k)：(1)0.00000000 (2)0.04596109(1)0.00000000 (2)0.04596109 (3)-0.01754976(4)-0.09170440 (5)0.06532115 (3)-0.01754976(4)-0.09170440 (5)0.06532115 44第一次累加值:(1)2.874000 (2)6.152000:(1)2.874000 (2)6.152000 (3)9.489000 (4)12.879000 (5)16.558000 (3)9.489000 (4)12.879000 (5)16.558000 55相对残差e(k)e(k)：(1)0.00000000 (2)0.01402108(1)0.00000000 (2)0.01402108 第60页/共82页7.3 销售额预测66原数据均值avg(x)avg(x)：3.31160000 3.31160000 77原数据方差 S(1)S(1)：88残差的均值avg(E)avg(E)：99残差的方差 S(2)S(2)：1010后验差比值:C:C：1111小误差概率 P P：1212模型计算值X(k)X(k)：(3)3.35454976 (4)3.48170440 (5)3.61367885(3)3.35454976 (4)3.48170440 (5)3.61367885 1313预测的结果X*(k)X*(k)：(6)4.51747233(6)4.51747233 预测精度等级：好！第61页/共82页7.4 城市道路交通事故次数城市道路交通事故次数 的灰色预测的灰色预测第62页/共82页 7.4 城市道路交通事故次数的灰色预测灰色理论以“部分信息已知、部分信息未知”的“小样本”、“贫信息”的不确定问题为研究对象,通过对“部分”已知的信息的生成开发,提取有价值的信息,构造生成序列的手段来寻求现实现象中存在的规律。交通事故作为一个随机事件,其本身具有相当大的偶然性和模糊性,如果把某地区的道路交通作为一个系统来看，则此系统中存在着一些确定因素(灰色系统称为白色信息),),如道路状况、信号标志,同时也存在一些不确定因素(灰色系统称为灰色信息)如车辆状况、气候因素、驾驶员心理状态等等,具有明显的不确定性特征。因此可以认为一个地区的道路交通安全系统是一个灰色系统,可以利用灰色系统理论进行研究。第63页/共82页 7.4 城市道路交通事故次数的灰色预测【例7.37.3】某市20042004年1-61-6月的交通事故次数统计见表7.5.7.5.试建立灰色预测模型.表7.5 7.5 交通事故次数统计解 利用GMGM预测软件计算，输出分析数据如下：原始数列(元素共6 6个):83):83，9595，130130，141141，156156，185185预测结果如下：第64页/共82页 7.4 城市道路交通事故次数的灰色预测1dx/dt+ax=u1dx/dt+ax=u：，22时间响应方程：33残差 E(k)E(k)：(1)0.00000000(2)-8.71441263(1)0.00000000(2)-8.71441263 (3)10.22065739(4)2.66733676 (3)10.22065739(4)2.66733676 (5)-3.75981586(6)0.49405494 (5)-3.75981586(6)0.49405494 44第一次累加值:(1)83.000000(2)178.000000 (3)308.000000:(1)83.000000(2)178.000000 (3)308.000000(4)449.00000(5)605.000000 (6)790.000000(4)449.00000(5)605.000000 (6)790.000000 55相对残差e(k)e(k)：第65页/共82页 7.4 城市道路交通事故次数的灰色预测66原数据均值avg(x)avg(x)：77原数据方差 S(1)S(1)：88残差的均值avg(E)avg(E)：99残差的方差 S(2)S(2)：1010后验差比值 C C：1111小误差概率 P P：1212模型计算值X(k)X(k)：(1)83.00000000(2)(1)83.00000000(2)103.71441263(3)119.77934261 (4)138.33266324 103.71441263(3)119.77934261 (4)138.33266324(5)159.75981586(6)184.50594506(5)159.75981586(6)184.50594506 1313预测的结果X*(k)X*(k)：(1)213.08514646(2)(1)213.08514646(2)246.09114698(3)284.20963932 (4)328.23252716 246.09114698(3)284.20963932 (4)328.23252716(5)379.07437672(6)437.79141674(7)(5)379.07437672(6)437.79141674(7)505.60348139 505.60348139 预测精度等级：好！这表明：如果该市不采取更有效的管制措施，7 7月的交通事故次数将上升至213213次.第66页/共82页7.5 城市火灾发生次数城市火灾发生次数 的灰色预测的灰色预测第67页/共82页7.5 城市火灾发生次数的灰色预测【例【例7.47.4】某市某市2001200520012005年火灾的统计数据见年火灾的统计数据见表表7.7.7.7.试建立模型，并对该市试建立模型，并对该市20062006年的火灾发年的火灾发生状况做出预测。生状况做出预测。表表7.7 7.7 某市某市2001200120052005年火灾数据年火灾数据年份年份2001200120022002200320032004200420052005 火灾火灾(起起)87879797120120166166161161第68页/共82页7.5 城市火灾发生次数的灰色预测解 利用GMGM预测软件计算，输出分析数据如下：原始数列(元素共5 5个):87):87，9797，120120，166166，161161预测结果如下：1dx/dt+ax=u：，2时间响应方程：3残差 E(k)：(1)0.00000000 (2)-7.05165921 第69页/共82页7.5 城市火灾发生次数的灰色预测4 4 第一次累加值:(1)87.000000 (2)184.000000 (3)304.000000 :(1)87.000000 (2)184.000000 (3)304.000000 (4)470.000000 (5)631.000000 (4)470.000000 (5)631.000000 5 5 相对残差e(k)e(k)：(1)0.00000000 (2)-0.07269752 (3)-0.02437316(1)0.00000000 (2)-0.07269752 (3)-0.02437316 (4)0.12517381 (5)-0.06560050 (4)0.12517381 (5)-0.06560050 6 6 原数据均值avg(x)avg(x)：7 7 原数据方差 S(1)S(1)：8 8 残差的均值avg(E)avg(E)：9 9 残差的方差 S(2)S(2)：10 10 后验差比值 C C：11 11 小误差概率 P P：12 12 模型计算值X(k)X(k)：(1)87.00000000 (2)104.05165921 (1)87.00000000 (2)104.05165921 (3)122.92477940 (4)145.22114789 (5)171.56168104 (3)122.92477940 (4)145.22114789 (5)171.56168104 13 13 预测的结果X*(k)X*(k)：(1)202.67991837(2)239.44245045 (1)202.67991837(2)239.44245045 (3)282.87305194 (4)334.18119203 (5)394.79571611(6)(3)282.87305194 (4)334.18119203 (5)394.79571611(6)466.40463669 466.40463669 n预测精度等级：合格！n结果表明：如果该市不采取更有效的防火措施，20062006年的火灾事故次数约为 203 次.第70页/共82页7.6 灾变与异常值预测灾变与异常值预测第71页/共82页7.6 灾变与异常值预测灰色灾变与异常值预测指运用灰色动态模型，对系统变化过程中某个异常数值在未来什么时间还会出现进行的预测.由于这个异常值的出现经常对人类产生不利的影响，即造成灾害，如：某年降雨量低于300mm,300mm,便形成旱灾，使粮食生产歉收；某年发生蝗灾，农作物就要减产；破坏性地震、特大洪水、台风与海啸等自然灾害的发生，更是给人们的生活和生产带来巨大的损失.因此，对这一类事件发生的时间和程度进行预报，是很有实际意义的.第72页/共82页7.6 灾变与异常值预测1.1.灾变预的数学原理与特征 灾变预测与数据预测的不同点，在于它不是预测序列数据的量的变化，而是预测异常值或“灾变”点出现的时间，它是应用灰色区间（间隔）的预测而进行的。所以，灾变预测的基本要求是“定量求时”。灾变预测的数学原理描述如下：第73页/共82页7.6 灾变与异常值预测第74页/共82页7.6 灾变与异常值预测第75页/共82页7.6 灾变与异常值预测3.3.实际问题旱灾预测【例7.57.5】某地年降水量原始数据序列如表所示，根据多年的时间观测，每当年降水量小于430430440mm440mm时，该地区将发生旱灾.所以，选择阈值=435mm,=435mm,利用GM(1,1)GM(1,1)模型进行旱灾预报.第76页/共82页7.6 灾变与异常值预测表7.9 某地年降水量（mm）原始数据第77页/共82页7.6 灾变与异常值预测第78页/共82页7.6 灾变与异常值预测第79页/共82页7.6 灾变与异常值预测第80页/共82页7.1 灰色系统的定义和特点7.2 灰色系统的模型7.3 销售额预测7.4 城市道路交通事故次数的灰色预测7.5 城市火灾发生次数的灰色预测7.6 灾变与异常值预测本章内容回顾第81页/共82页感谢您的观看！第82页/共82页