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

《数学建模灰色预测模型重制.pptx》由会员分享，可在线阅读，更多相关《数学建模灰色预测模型重制.pptx（74页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、预测问题：（可用灰色系统分析法）预测问题：（可用灰色系统分析法）CUMCM2003A SARS的传播CUMCM2005A 长江水质的评价和预测CUMCM2006A 出版社的资源配置 CUMCM2007A 中国人口增长预测MCM2016E 水资源预测 第1页/共75页CUMCM2003A SARS的传播SARS（Severe Acute Respiratory Syndrome，严重急性呼吸道综合症,俗称：非典型肺炎）是21世纪第一个在世界范围内传播的传染病。SARS的爆发和蔓延给我国的经济发展和人民生活带来了很大影响，我们从中得到许多重要的经验和教训，认识到定量地研究传染病的传播规律、为预测和

2、控制传染病蔓延创造条件的重要性。请你们对SARS 的传播建立数学模型，具体要求如下：（1）对附件1所提供的一个早期的模型，评价其合理性和实用性。第2页/共75页CUMCM2003A SARS的传播（2）建立你们自己的模型，说明为什么优于附件1中的模型；特别要说明怎样才能建立一个真正能够预测以及能为预防和控制提供可靠、足够的信息的模型，这样做的困难在哪里？对于卫生部门所采取的措施做出评论，如：提前或延后5天采取严格的隔离措施，对疫情传播所造成的影响做出估计。附件2提供的数据供参考。（3）收集SARS对经济某个方面影响的数据，建立相应的数学模型并进行预测。附件3提供的数据供参考。第3页/共75页

3、水是人类赖以生存的资源，保护水资源就是保护我们自己，对于我国大江大河水资源的保护和治理应是重中之重。专家们呼吁：“以人为本，建设文明和谐社会，改善人与自然的环境，减少污染。”长江是我国第一、世界第三大河流，长江水质的污染程度日趋严重，已引起了相关政府部门和专家们的高度重视。2004年10月，由全国政协与中国发展研究院联合组成“保护长江万里行”考察团，从长江上游宜宾到下游上海，对沿线21个CUMCM2005A 长江水质的评价和预测第4页/共75页 重点城市做了实地考察，揭示了一幅长江污染的真实画面，其污染程度让人触目惊心。为此，专家们提出“若不及时拯救，长江生态10年内将濒临崩溃”（附件），并发

4、出了“拿什么拯救癌变长江”的呼唤（附件2）。附件3给出了长江沿线17个观测站（地区）近两年多主要水质指标的检测数据，以及干流上个观测站近一年多的基本数据（站点距离、水流量和水流速）。通常认为一个观测站（地区）的水质污染主要来自于本地区的排污和上游的污水。CUMCM2005A 长江水质的评价和预测第5页/共75页 一般说来，江河自身对污染物都有一定的自然净化能力，即污染物在水环境中通过物理降解、化学降解和生物降解等使水中污染物的浓度降低。反映江河自然净化能力的指标称为降解系数。事实上，长江干流的自然净化能力可以认为是近似均匀的，根据检测可知，主要污染物高锰酸盐指数和氨氮的降解系数通常介于0.10

5、.5之间，比如可以考虑取0.2(单位：1/天)。附件4是“19952004年长江流域水质报告”给出的主要统计数据。下面的附表是国标(GB3838-2002)给出的地表水环境质量标准中4个主要项目标准限值，其中、类为可饮用水。CUMCM2005A 长江水质的评价和预测第6页/共75页CUMCM2005A 长江水质的评价和预测请你们研究下列问题：（1）对长江近两年多的水质情况做出定量的综合评价，并分析各地区水质的污染状况。（2）研究、分析长江干流近一年多主要污染物高锰酸盐指数和氨氮的污染源主要在哪些地区?（3）假如不采取更有效的治理措施，依照过去10年的主要统计数据，对长江未来水质污染的发展趋势做

6、出预测分析，比如研究未来10年的情况。（4）根据你的预测分析，如果未来10年内每年都要求长江干流的类和类水的比例控制在20%以内，且没有劣类水,那么每年需要处理多少污水？（5）你对解决长江水质污染问题有什么切实可行的建议和意见。第7页/共75页灰色系统分析方法灰色系统分析方法灰色系统分析方法灰色系统分析方法8 灰色预测的概念；灰色生成数列；案例：SARS疫情对某些经济指标影响。灰色模型GM；第8页/共75页 灰色系统是由华中科技大学的邓聚龙教授80年代初所创立，在短短的三十年里已得到了长足的发展。特点是模型使用的不是原始数据序列，而是生成数据序列。其核心是灰色模型，即对原始数据作累加生成（或其

7、他方法生成）得到的近似指数规律再进行建模的方法。灰色系统研究的是“部分信息明确，部分信息未知”的“小样本，贫信息”不确定性问题，并依据信息覆盖，通过序列算子的作用探索事物运动的现实规律。其特点是“少数据建模”。着重研究“外延明确，内涵不明确”的对象。灰色系统介绍第9页/共75页1、不需要很多的数据。一般只需要4个数据，就能解决历史数据少、序列的完整性及可靠性低的问题。2、能将无规律的原始数据进行生成得到规律性较强的生成序列，运算简便，易于检验，不考虑分布规律，不考虑变化趋势。3、定量分析结果与定性分析结果不会不一致。4、能利用微分方程来充分挖掘系统的本质，精准度高。灰色模型的优点第10页/共7

8、5页1、只适用于中短期预测2、只适合指数增长的预测 灰色模型的缺点第11页/共75页 灰色预测的概念：灰色系统、白色系统和黑色系统 白色系统是指一个系统的内部特征是完全 已知的，即系统的信息是完全充分的。黑色系统是指一个系统的内部信息对外界 来说是一无所知的，只能通过它与外界的 联系来加以观测研究。灰色系统内的一部分信息是已知的，另一部分信息是未知 的，系统内各因素间有不确定的关系。定义：系统作为一个包含若干相互关联、相互制约的任意种类元素组成的具有某种特定功能的整体。系统内部存在有物质流、信息流、能量流。系统（根据信息明确程度）黑色系统（信息毫无所知或知之甚少）灰色系统（既含有已知信息又有未

9、知信息）白色系统（信息完全明确）第12页/共75页 灰色预测法是一种对含有不确定因素的系统进行预测的方法。灰色预测是对既含有已知信息又含有不确定信息的系统进行预测，就是对在一定范围内变化的、与时间有关的灰色过程进行预测。灰色预测通过鉴别系统因素之间发展趋势的相异程度，即进行关联分析，并对原始数据进行生成处理来寻找系统变动的规律，生成有较强规律性的数据序列,然后建立相应的微分方程模型，从而预测事物未来发展趋势的状况。灰色预测法用等时距观测到的反映预测对象特征的一系列数量值构造灰色预测模型，预测未来某一时刻的特征量，或达到某一 特征量的时间。灰色预测法：与时间密切相关！很好理解第13页/共75页

10、灰色预测的四种常见类型 灰色时间序列预测：即用观察到的反映预测对象特征的时间序列来构造灰色预测模型，预测未来某一时刻的特征量，或达到某一特征量的时间。（常见）畸变预测：即通过灰色模型预测异常值出现的时刻，预测异常值什么时候出现在特定时区内。系统预测：通过对系统行为特征指标建立一组相互关联的灰色预测模型，预测系统中众多变量间的相互协调关系的变化。拓扑预测：将原始数据做曲线，在曲线上按定值寻找该定值发生的所有时点，并以该定值为框架构成时点数列，然后建立模型预测该定值所发生的时点。第14页/共75页公理1、差异信息原理。“差异”是信息，凡信息必有差异。公理2、解的非唯一性原理。信息不完全、不确定的解

11、是非唯一的。该原理是灰色系统理论解决实际问题所遵循的基本法则。公理3、最少信息原理。灰色系统理论的特点是充分开发利用已占有的“最少信息”。公理4、认知根据原理。信息是认知的根据。公理5、新信息优先原理。新信息对认知的作用大于老信息。公理6、灰性不灭原理。“信息不完全”是绝对的。灰色系统的基本原理第15页/共75页灰色系统的描述：灰色系统用灰色参数、灰色方程、灰色矩阵、灰色度等综合描述，其中灰数是灰色系统的基本单元。第16页/共75页 只知道大概范围而不知道其确切值的数称为灰数（只知道部分数学特征，而不知道具体数值的参数），通常记为：“”。例如：1.头发的多少才算是秃子。应该是个区间范围。模糊

12、2.多少层的楼房算高楼，中高楼，低楼。3.多么重才算胖子？。灰数及其运算第17页/共75页例如：“某人的身高约为170cm、体重大致为60kg”，这里的“（约为）170（cm）”、“60”都是灰数，分别记为 、。又如，“那女孩身高在157160cm之间”，则关于身高的灰数 。记为灰数的白化默认数，简称白化数。在灰色系统理论中，把随机变量看成灰数，即是在指定范围内变化的所有白色数的全体。如代购一件价格为100元左右的衣服，100可作为预购衣服价格的白化值。灰数有离散灰数（属于离散集）和连续灰数（属于某一区间）。第18页/共75页灰色代数方程含有灰色系数的代数方程如：灰色微分方程为含有灰色导数或灰

13、色微分的方程，如 灰色矩阵行列数确知而含有灰元的矩阵 若在A的m*n个元素中，有N个灰色元素，则可以用d表示这一矩阵的灰色度第19页/共75页灰色生成数列灰色生成数列 灰色系统理论认为，尽管客观表象复杂，但总是有整体功能的，因此必然蕴含某种内在规律。关键在于如何选择适当的方式去挖掘和利用它。灰色系统是通过对原始数据的整理来寻求其变化规律的，这是一种就数据寻求数据的现实规律的途径，即为灰色序列的生成。一切灰色序列都能通过某种生成弱化其随机性，显现其规律性。数据生成的常用方式有累加生成、累减生成和加权累加生成。第20页/共75页灰色序列生成是一种通过对原始数据的挖掘、整理来寻求数据变化的现实规律的

14、途径，简称灰生成。灰生成特点在保持原序列形式的前提下，改变序列中数据的值与性质。一切灰色序列都能通过某种生成弱化其随机性，显现其规律性。灰生成的作用1.统一序列的目标性质，为灰决策提供基础。2.将摆动序列转换为单调增长序列，以利于灰建模。3.揭示潜藏在序列中的递增势态，变不可比为可比序列。灰生成技术第21页/共75页常见的几种灰生成类型：1.累加生成算子（AGO）2.逆累加生成算子（IAGO）3.均值生成算子（MEAN）4.级比生成算子第22页/共75页1.累加生成算子（AGO）定义 它是对原序列中的数据依次累加以得到生成序列。令 为原序列我们说 是 的AGO序列，并记为当且仅当并满足第23页

15、/共75页例1 摆动序列为：通过AGO可以加工成单调增序列：第24页/共75页累加生成 把数列各项（时刻）数据依次累加的过程称为累加生成过程（AGO）。由累加生成过程所得的数列称为累加生成数列。设原始数列为 ，令称所得到的新数列为数列 的1次累加生成数列。类似地有称为 的r次累加生成数列。第25页/共75页累加生成计算示例累加生成计算示例例：x(0)=(x(0)(k)k=1，2，3，4，5)=x(0)(1)，x(0)(2)，x(0)(3)，x(0)(4)，x(0)(5)=(3.2，3.3，3.4，3.6，3.8)求 x(1)(k)解：第26页/共75页第27页/共75页累加生成的特点累加生成的

16、特点 一般经济数列都是非负数列。累加生成能使任意非负数列、摆动的与非摆动的，转化为非减的、递增的。原始数列作图 1AGO作图 某市的汽车销售量 递增的规律 第28页/共75页原始数列作图 1AGO作图 有明显的指数关系的规律 某钢厂产量某地区作物产量 s型变化规律 第29页/共75页2.逆累加生成算子（IAGO）定义 它是对AGO生成序列中相邻数据依次累减，又称累减生成。令 为原序列称 是 的IAGO序列，并记为当且仅当并 满足第30页/共75页例2令原始序列 为这表明第31页/共75页累减生成对于原始数据列依次做前后相邻的两个数据相减的运算过程称为累减生成过程IAGO。如果原始数据列为令 称

17、所得到的数列 为 的1次累减生成数列。注：从这里的记号也可以看到，从原始数列 ，得到新数列 ，再通过累减生成可以还原出原始数列。实际运用中在数列 的基础上预测出 ，通过累减生成得到预测数列 。第32页/共75页累减生成计算示例第33页/共75页3.均值生成算子（MEAN）定义 它是将AGO序列中前后相邻两数取平均数，以获得生成序列。令 为 的AGO序列称 为 的MEAN序列，并记为当且仅当并且每个 满足下述关系第34页/共75页例3 对于 ，有 第35页/共75页4.级比生成算子定义 设序列 ，则称为序列 的级比。检验准则 设序列 的级比满足时，序列 可做GM(1,1)建模。第36页/共75页

18、加权邻值生成设原始数列为称 为数列 的邻值。为后邻值，为前邻值，对于常数 ，令 由此得到的数列 称为数列 在权 下的邻值生成数，权 也称为生成系数。特别地，当生成系数 时，则称为均值生成数，也称等权邻值生成数。第37页/共75页灰色系统理论的主要方法关联度分析法最基本的方法（一个由众多因素构成的系统中哪些因素对系统的影响大/中/小？）基于白化权函数的灰色统计和灰色聚类法。灰色预测法（如GM（1，1）。灰色决策。灰色优化技术（如灰色规划等）。第38页/共75页GM(1,1)的符号含义：G M (1,1)Grey灰色Model模型1阶方程1个变量第39页/共75页灰色理论认为系统的行为现象尽管是朦

19、胧的，数据是复杂的，但它毕竟是有序的，是有整体功能的。灰数的生成，就是从杂乱中寻找出规律。同时，灰色理论建立的是生成数据模型，不是原始数据模型。因此，灰色预测的数据是通过生成数据的GM(1,1)模型所得到的预测值的逆处理结果。第40页/共75页灰色预测模型GM(m,n)灰色预测是应用灰色模型GM对灰色系统进行分析、建模、求解、预测的过程。由于灰色建模理论应用数据生成手段，弱化了系统的随机性，使紊乱的原始序列呈现某种规律，规律不明显的变得较为明显，建模后还能进行残差辨识，即使较少的历史数据，任意随机分布，也能得到较高的预测精度。因此，灰色预测在社会经济、管理决策、农业规划、气象生态等各个部门和行

20、业都得到了广泛的应用 第41页/共75页 灰色系统理论是基于关联空间、光滑离散函数等概念定义灰导数与灰微分方程，进而用离散数据列建立微分方程形式的动态模型，即灰色模型是利用离散随机数经过生成变为随机性被显著削弱而且较有规律的生成数，建立起的微分方程形式的模型，这样便于对其变化过程进行研究和描述。vG表示grey（灰色），M表示model（模型）灰色模型(1,1)第42页/共75页设 为原始数列，其1次累加生成数列为 ，其中定义 的灰导数为令 为数列 的邻值生成数列，即于是定义GM（1，1）的灰微分方程模型为第43页/共75页即或 （1）在式（1）中，称为灰导数，a称为发展系数，称为白化背景值，

21、b称为灰作用量。将时刻表 代入（1）式有引入矩阵向量记号：第44页/共75页于是GM（1，1）模型可表示为现在问题归结为求a,b在值。用一元线性回归，即最小二乘法求它们的估计值为注：实际上回归分析中求估计值是用软件计算的，有标准程序求解，如matlab等。GM（1，1）的白化型对于GM（1，1）的灰微分方程（1），如果将灰导数 的时刻 视为连续变量t，则 视为时间t函数 ，于是 对应于导数量级 ，白化背景值 对应于导数 。于是GM（1，1）的灰微分方程对应于的白微分方程为 （2）第45页/共75页GM（1，1）灰色预测的步骤）灰色预测的步骤1.数据的检验与处理为了保证GM（1，1）建模方法的可

22、行性，需要对已知数据做必要的检验处理。设原始数据列为了 ，计算数列的级比如果所有的级比都落在可容覆盖区间 内，则数据列 可以建立GM（1，1）模型且可以进行灰色预测。否则，对数据做适当的变换处理，如平移变换：取C使得数据列的级比都落在可容覆盖内。第46页/共75页2.建立GM（1，1）模型 不妨设 满足上面的要求，以它为数据列建立GM（1，1）模型用回归分析求得a,b的估计值，于是相应的白化模型为 解为 （3）于是得到预测值从而相应地得到预测值：第47页/共75页3.检验预测值（1）残差检验：计算相对残差如果对所有的 ，则认为达到较高的要求：否则，若对所有的 ，则认为达到一般要求。（2）级比偏

23、差值检验：计算如果对所有的 ，则认为达到较高的要求；否则若对所有的 ，则认为达到一般要求。第48页/共75页灰色预测计算实例例 北方某城市19861992 年道路交通噪声平均声级数据见表6表6 市近年来交通噪声数据dB(A)序号年份 eq L 1 1986 71.12 1987 72.43 1988 72.44 1989 72.1 5 1990 71.4 6 1991 72.07 1992 71.6第一步:级比检验建立交通噪声平均声级数据时间序列如下：=(71.1,72.4,72.4,72.1,71.4,72.0,71.6)第49页/共75页（1）求级比(k)=(0.982,1,1.0042,

24、1.0098,0.9917,1.0056)（2）级比判断由于所有的(k)0.982,1.0098，k=2,3,7，故可以用x(0)作满意的GM（1，1）建模。第二步:GM（1，1）建模（1）对原始数据 作一次累加，即 =(71.1,143.5,215.9,288,359.4,431.4,503）（2）构造数据矩阵B 及数据向量Y第50页/共75页（3）计算u 于是得到a=0.0023，b=72.6573。（4）建立模型第51页/共75页求解得（5）求生成数列值 及模型还原值 ：令k=1,2,3,4,5,6,由上面的时间响应函数可算得 ，其中取由 取k=2,3,4,7，得 (71.1,72.4,

25、72.2,72.1,71.9,71.7,71.6)第52页/共75页第三步:模型检验模型的各种检验指标值的计算结果见表 7.表7 GM(1,1)模型检验表序号 年份 原始值 模型值 残差 相对误差 级比偏差1 1986 71.1 71.1 0 02 1987 72.4 72.4 -0.0057 0.01%0.00233 1988 72.4 72.2 0.1638 0.23%0.02034 1989 72.1 72.1 0.0329 0.05%-0.00185 1990 71.4 71.9 -0.4984 0.7%-0.00746 1991 72.0 71.7 0.2699 0.37%0.010

26、77 1992 71.6 71.6 0.0378 0.05%-0.0032经验证，该模型的精度较高，可进行预测和预报。第53页/共75页案例：案例：SARS疫情对某些经济指标影响疫情对某些经济指标影响1 问题的提出2003年的SARS疫情对中国部分行业的经济发展产生了一定的影响，特别是对帮分疫情较严重的省市的相关行业所造成的影响是明显的，经济影响主要分为直接经济影响和间接影响直接经济影响涉及到商品零售业、旅游业、综合服务等行业很多方面难以进行定量地评估，现仅就SARS疫情较重的某市商品零售业、旅游业和综合服务业的影响进行定量的评估分析究竟SARS疫情对商品零售业、旅游业和综合服务业的影响有多大

27、，已知该市从1997年1月到2003年10月的商品零售额、接待旅游人数和综合服务收入的统计数据如下表1、表2、表3第54页/共75页表表1 商品的零售额商品的零售额(单位：亿元单位：亿元)年代 1月 2月 3月 4月 5月 6月 7月 8月 9月 10月 11月 12月1997 830 798 781 851 866 882 903 867 933 925 909 9691998 1017 851 878 916 934 945 974 995 1042 1023 1010 12351999 922 1140 933 1010 1035 1052 1095 1092 1096 1112 121

28、7 13132000 1050 1257 1066 1160 1176 1180 1217 1187 12O2 1278 1218 12192001 1393 1295 1225 1245 1357 13O8 1387 1337 1368 1389 1296 13372002 1375 1353 1330 1334 1428 1416 1429 1473 1596 1621 1535 15592003 1632 1597 1584 1452 124 1441 157 1626 1718 1807 1735 1765第55页/共75页表表2 接待海外旅游人数接待海外旅游人数(单位：万人单位：万人

29、)年代 1月 2月 3月 4月 5月 6月 7月 8月 9月 1O月 l1月 12月1997 94 113 168 198 2O3 188 2O9 249 247 243 194 1861998 96 117 158 199 195 178 178 233 214 245 2O1 1591999 1O1 129 177 210 210 2O4 219 258 293 298 236 1652000 114 260 196 259 276 243 230 278 273 285 328 1852001 115 264 2O4 261 289 280 252 3O8 287 281 222 2O7

30、2002 137 297 231 289 290 274 260 322 314 326 292 2292003 154 171 235 l16 178 261 88 162 2O1 249 265 218第56页/共75页表表3 综合服务业累计数据综合服务业累计数据(单位：亿元单位：亿元)年代 2月 3月 4月 5月 6月 7月8月 9月 1O月 11月 12月l997 96 l44 l94 276 383 466 554 652 747 832 9721998 111 169 235 400 459 565 695 805 881 1011 1139l999 151 238 335 425

31、 541 641 739 866 975 1087 12382000 164 263 376 531 600 711 913 1O38 1173 1296 14972001 182 318 445 576 708 856 1000 1 145 1292 1435 16672002 216 361 504 642 818 979 1142 1305 1479 1644 19202003 241 404 584 741 923 1114 1298 1492 1684 1885 2218试根据这些历史数据建立预测评估模型，评估2003年SARS疫情给该市的商品零售业、旅游业和综合服务业所造成的影响第

32、57页/共75页2 模型的分析与假设模型的分析与假设根据所掌握的历史统计数据可以看出，在正常情况下，全年的平均值较好地反映了相关指标的变化规律，这样可以把预测评估分成两部分：(i)利用灰色理论建立灰微分方程模型，由19972002年的平均值预测2003年平均值；(ii)通过历史数据计算每个月的指标值与全年总值的关系，从而可预测出正常情况下2003年每个月的指标值，再与实际值比较可以估算出SARS疫情实际造成的影响 第58页/共75页2 模型的分析与假设模型的分析与假设给出下面两条假设：(1)假设该市的统计数据都是可靠准确的；(2)假设该市在SARS疫情流行期间和结束之后，数据的变化只与SARS

33、疫情的影响有关，不考虑其他随机因素的影响第59页/共75页3 建立灰色预测模型建立灰色预测模型GM(1,1)由已知数据，对于19972002年某项指标记为矩阵 计算每年的年平均值，记为 (3)并要求级比 对 作一次累加，则 记 (4)取 的加权均值，则 为确定参数,于是GM(1,1)的白化微分方程模型为 (5)其中a是发展灰度,b是内生控制灰度.第60页/共75页3 建立灰色预测模型建立灰色预测模型GM(1,1)由于 ，取 为灰导数，为 背景值，则建立灰微分方程为:或 其矩阵形式为 ，其中 用最小二乘法求得参数的估计值为 (6)第61页/共75页3 建立灰色预测模型建立灰色预测模型GM(1,1

34、)则灰微分方程模型(4)的解为则 (7)由(7)式可以得到2003年的平均值为 ，则预测2003年的总值为 根据历史数据，可以统计计算出2003年第i个月的指标值占全年总值的比例为 ，即 (8)则 于是可得2003年每一个月的指标值为Y=X u 第62页/共75页4 模型的求解模型的求解 (i)商品零售额由数据表1，计算可得每年月平均值、一次累加值分别为 =(876167，985000，1084750，1184167，1328083，1454083)，=(876167，1861167，2945917，4130083，5458167，6912250)显然 的所有级比都在可行域内经检验，在这里取参

35、数 比较合适，则有 (1270167，2295067，3419583，4661317，6039800)第63页/共75页(i)商品零售额商品零售额由最小二乘法求得a=-00993，b=855985可得2003年的月平均值为 =1628826亿元；年总值为=19546亿元 由(8)式得每月的比例为u=(00794，00807，00749，00786，00819，00818，00845，00838，00872，00886，00866，00920)故2003年112月的预测值为Y=uX=(1552，1578，1464，1536，1601，1599，1652，1638，17o5，1732，1693，1

36、799)(亿元)第64页/共75页将预测值与实际统计值进行比较如下表4所示 表4 2003年商品的零售额(单位：亿元)月份 1月 2月 3月 4月 5月 6月 7月 8月 9月 10月l1月 12月预测值 1552 1578 1464 1536 1601 1599 1652 1638 1705 1732 1693 1799实际值 1632 1597 1584 1452 1240 1441 1570 1626 1718 18O7 1735 1765第65页/共75页clc,clearhan1=83.0 79.8 78.1 85.1 86.6 88.2 90.3 86.7 93.3 92.5 90

37、.9 96.9101.7 85.1 87.8 91.6 93.4 94.5 97.4 99.5 104.2 102.3 101.0 123.592.2 114.0 93.3 101.0 103.5 105.2 109.5 109.2 109.6 111.2 121.7 131.3105.0 125.7 106.6 116.0 117.6 118.0 121.7 118.7 120.2 127.8 121.8 121.9139.3 129.5 122.5 124.5 135.7 130.8 138.7 133.7 136.8 138.9 129.6 133.7137.5 135.3 133.0

38、133.4 142.8 141.6 142.9 147.3 159.6 162.1 153.5 155.9163.2 159.7 158.4 145.2 124.0 144.1 157.0 162.6 171.8 180.7 173.5 176.5;han1(end,:)=;m=size(han1,2);x0=mean(han1,2);x1=cumsum(x0);alpha=0.4;n=length(x0);z1=alpha*x1(2:n)+(1-alpha)*x1(1:n-1)Y=x0(2:n);B=-z1,ones(n-1,1);ab=BYk=6;x7hat=(x0(1)-ab(2)/ab

39、(1)*(exp(-ab(1)*k)-exp(-ab(1)*(k-1)z=m*x7hatu=sum(han1)/sum(sum(han1)v=z*u第66页/共75页x1=87.6167 186.1167 294.5917 413.0083 545.8167 691.2250z1=127.0167 229.5067 341.9583 466.1317 603.9800ab=-0.0993 85.5985x7hat=162.8793 z=1.9546e+003u=0.0794 0.0807 0.0749 0.0786 0.0819 0.0818 0.0845 0.0838 0.0872 0.08

40、86 0.0866 0.0920v=155.2152 157.7365 146.4023 153.5421 160.1400 159.8337 165.0649 163.7924 170.5317 173.1473 169.3064 179.8394第67页/共75页(ii)接待海外旅游人数接待海外旅游人数由数据表2，计算每年月平均值 和一次累加值 取参数 ，可得加权平均值 进而求得a=-00938，b=162670，=302649，X=3631788，以及u=(00407，00732，00703，00878，00907，00848，00836，01022，01010，01041，00914，

41、00701)于是可得到2003年的接待海外旅游人数的预测值，并与实际值比较如下表5所示第68页/共75页 表5 2003年接待海外旅游人数(单位：万人)月份 1月 2月 3月 4月 5月 6月 7月8月 9月 10月 11月 12月预测值 14.8 26.6 25.5 31.9 33.0 30.8 30.4 37.1 36.7 37.8 33.2 25.5实际值 15.4 17.1 23.5 11.6 1.78 2.61 8.8 16.2 20.1 24.9 26.5 21.8(iii)综合服务收入类似处理(略)第69页/共75页5 模型的结果分析模型的结果分析根据该市的统计报告显示，2003

42、年4、5、6三个月的实际商品零售额分别为1452，124、1441亿元在这之前，根据统计部门的估计4、5、6三个月份SARS疫情对该市的商品零售业的影响最为严重，这三个月估计大约损失62亿元左右从我们的模型预测结果来计算，4、5、6三个月的损失为603亿元，这个数据基本与专家的估算值相符，8月份基本恢复正常，这也说明了模型的正确性和可靠性对于旅游业来说是受影响最严重的行业之一，最严重的4、5、6、7四个月就损失100多万人，按最新统计数据，平均每人消费1002美元计算，大约损失1O亿美元全年大约损失160万人，约合16亿美元，到年底基本恢复正常 第70页/共75页对于综合服务业中的部分行业影响较大，如航空交通运输、宾馆餐饮等，但有些行业影响不大，如电信、通讯等，总平均来看，影响还不算太大，5、6、7、8四个月大约损失7O亿元从预测结果可以看出，虽然下半年没有发生疫情，但人们一直担心SARS会卷土重来，所以，对这些行业还是有一定的影响，即SARS影响的延续性的作用该模型虽是就某经济指标的发展规律进行评估预测而建立的，但类似的也适用于其他方面的一些 数据规律的评估预测问题，即该模型具有很广泛的应用性第71页/共75页四种灰色预测模型对比第72页/共75页四种灰色预测模型对比第73页/共75页谢谢大家！预祝大家取得好成绩！第74页/共75页