《灰色系统理论讲》.ppt

  1. 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
  2. 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  3. 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
灰色系统分析方法;主要内容;1 有关概念;; 该理论以“部分信息已知,部分信息未知”的小样本、信息不确定性系统为研究对象,主要通过对“部分”已知信息的生成、开发,提取出有价值的信息,从而实现对系统运行行为、演化规律的正确描述和有效监控。 目前,灰色系统理论应用范围已拓展到工业、农业、社会、经济、能源、地质、石油等众多科学领域,成功地解决了生产、生活和科学研究中的大量实际问题,取得了显著成果。 ; ;灰数的白化;灰色关联分析;;;;多因子的情况;;;灰色生成??列;;累减生成列;均值生成数列;;GM(1,1)模型;;;;;GM(1,N)的定义;;;;;2 灰色系统理论的应用范畴? ;3、建立灰色系统 GM(1,1)模型的步骤;2、建立模型GM(1,1);1、灰色预测方法;;例题 设有原始数据序列 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进行模拟;第3步:对参数列 进行最小二乘估计,得 =[-0.0372 3.06536] 第4步:确定模型为 时间响应式为 第5步:求X1的模拟值;第6步:还原出 的模拟值,由 得 对比原数据;3、检验预测值;4 如何利用灰色系统实现预测?;灰色灾害预测应用;定义:设原始序列X=(x(1),x(2),…,x(n)),给定上限异常值(灾变值)a,称X的子序列 Xa(x(q(1)), x(q(2)),…… x(q(m)))={x(q(i))=a, i=1,2,…,m}为上灾变序列。同理,可定义下灾变序列。二者统一称为灾变序列。;定义:设X为原始序列,称 Q0=(q(1),q(2),……q(m))为灾变日期序列。 灾变预测就是要通过对灾变日期序列的研究,寻找其规律性,预测以后若干次灾变发生的日期,灰色系统的灾变预测是通过对灾变日期序列建立GM(1,1)模型实现的。;例 某地区最近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时认为旱灾发生,试根据上述数据预测下一次旱灾发生在几年后? ;解:取灾变值为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,28,43.5);设q(k)+az1(k)=b,易知B,Y,由最小二乘法得[a,b]=[-0.253661 6.258339] 故灾变日期序列的GM(1,1)序号响应式为 即 ;由此可得Q0的模拟序列为 由 得残差序列为 ;再由相对误差序列 由此可计算出平均相对误差为;平均相对精度为1- =97.81%,故可用 进行预测, 即从最近一次旱灾发生的日期算起,5年以后,可能发生旱灾。 为了提高预测的可靠程度,可以取若干个不同的异常值,建立多个模型进行预测。 ;;;计算的 MATLAB 程序如下: clc,clear a=[390.6,412,320,559.2, 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=B\Y y=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) %为提

文档评论(0)

清风老月 + 关注
官方认证
内容提供者

该用户很懒,什么也没介绍

认证主体阳春市惠兴图文设计有限公司
IP属地广东
统一社会信用代码/组织机构代码
91441781MA53BEWA2D

1亿VIP精品文档

相关文档