灰色模型1.docx

上传人:王** 文档编号:62432 上传时间:2022-12-12 格式:DOCX 页数:19 大小:95.88KB
下载 相关 举报
灰色模型1.docx_第1页
第1页 / 共19页
灰色模型1.docx_第2页
第2页 / 共19页
灰色模型1.docx_第3页
第3页 / 共19页
灰色模型1.docx_第4页
第4页 / 共19页
灰色模型1.docx_第5页
第5页 / 共19页
灰色模型1.docx_第6页
第6页 / 共19页
灰色模型1.docx_第7页
第7页 / 共19页
灰色模型1.docx_第8页
第8页 / 共19页
灰色模型1.docx_第9页
第9页 / 共19页
灰色模型1.docx_第10页
第10页 / 共19页
亲,该文档总共19页,到这儿已超出免费预览范围,如果喜欢就下载吧!
资源描述

《灰色模型1.docx》由会员分享,可在线阅读,更多相关《灰色模型1.docx(19页珍藏版)》请在优知文库上搜索。

1、灰色预测模型灰色预测是就灰色系统所做的预测.所谓灰色系统是介于白色系统和黑箱 系统之间的过渡系统,其具体的含义是:如果某一系统的全部信息已知为白色系 统,全部信息未知为黑箱系统,部分信息已知,部分信息未知,那么这一系统就 是灰箱系统.一般地说,社会系统、经济系统、生态系统都是灰色系统.灰色系统理论认为对既含有已知信息又含有未知或非确定信息的系统进行 预测,就是对在一定方位内变化的、与时间有关的灰色过程的预测.尽管过程中 所显示的现象是随机的、杂乱无章的,但毕竟是有序的、有界的,因此这一数据 集合具备潜在的规律,灰色预测就是利用这种规律建立灰色模型对灰色系统进行 预测.灰色预测模型只需要较少的观

2、测数据即可,这和时间序列分析,多元回归分 析等需要较多数据的统计模型不一样.因此,对于只有少量观测数据的项目来 说,灰色预测是一种有用的工具.一、GM(Ij)模型多色系统理论是邓聚龙教授在1981年提出来的,是一种对含有不确定因索 彳进行预测的方法.通.过鉴别系统因素之间发展趋势的相异程度,进行关联分 析,并通过对原始数据进行生成处理来寻找系统的变化规律,生成较强规律性数 据字歹U,然后建立相应微分方程模型,从而预测事物未来的发展&势和未来状态.GM(IjI:M,,.代忤差分方拧,按时间累加后所形成的新的时间序列呈现的规律可用一阶线性微 分方程的解来逼近,用生成数序列代替原始时间序列,弱化原始

3、时间序列的随机 性.这样可以对变化过程作较长时间的描述,进而建立微分方程形式的模型.其 建模的实质是建立微分方程的系数,将时间序列转化为微分方程,通过灰色微分 方程口以建立抽象系统的发展模型.经证明,经阶线性微分方程的解逼近所揭.始时间数.J .功的.1.1 GM(Ll)模型的建立成为有规律的时间序列数据,对其作一次累加生成运算,即令XSSn=l从而得到新的生成数列X,X=卜(1),*(2),数列 则生成的离散形式的微分方程具体的形式为即表示求解上述微分方程,解为x(t) = cea+-a当尸1时,x(r) = x(l), KPc = X(I)-,则可根据上述公式得到离散形式微分方程 a的具体

4、形式为x(r) = pl)-,)+其中,G项中的X为公的背景值,也称初始值;a, 是待识别的灰色参数, 出按白化导数定义有dx 1. x(t+t)-x(t)=Iimdt 3,显然,当时间密化值定义为1时,当41时,则上式可记为dx=lm(x( +,)- x(f)dt Af这表明虫是一次累减生成的,因此该式可以改写为dt= (r + l)-(r)dt当,足够小时,变量JV从Xa)到、(,+/)是不会出现突变的,所以取力,)与+,)的平均值作为当,足够小时的背景值,即/)=;/() +“(+ 1),匚1: 值带入式子,整理得?0)(r 1) = -g 力 + X (Z +1) + w由其离散形式可

5、得到如下矩阵:二x+ X f(2)2 L令y = ()(2),x (3),0n)r-扪+X一如(I)+ x S) 1 = ( P称Y为数据向量,8为数据矩阵,为参数向量.则上式可简化为线性模型:Y = Ba由最小二乘估计方法得a = pzl = (rl B上式即为GM(Ij)参数出的矩阵辨识算式,式中(3%厂/丫事实上是数据矩阵 8的广义逆矩阵.将求得的。,值代入微分方程的解式,则,r) = wi)-Vfl-0+-a a其中,上式是GM(IJ)模型的时间响应函数形式,将它离散化得xm(Z) = f0i(l)-L-0+-(0)对序列外)(。再作累减生成可进行预测.即(O = Jf(O-iO-l)

6、=卜用(IY上式便是GM(1,1)模型的预测的具体计算式.或对x。) = ce + 2求导还原得 aa1.2 GM(1,1)模型的检验GM(1,1)模型的检验包括残差检验、关联度检验、后验差检验三种形式.每种检验对应不同功能:残差检验属于算术检验,对模型值和实际值的误差 进行逐点检验;关联度检验属于几何检验范围,通过考察模型曲线与建模序列曲 线的几何相似程度进行检验,关联度越大模型越好;后验差检验属于统计检验, 对残差分布的统计特性进行检验,衡量灰色模型的精度.残差检验残差大小检验,即对模型值和实际值的残差进行逐点检验.设模拟值的残差序列为O,则e(o)(r) = x(o)(r)-x(o)(r

7、)令(/)为残差相对值,即残差百分比为) =x(o)(O-x(o)(O令A为平均残差,X = -f(t). /=I1 2设残差的方差为s;,则=一2卜(力-可.故后验差比例C为c = S2s, n /=I误差频率 P为 P= Pet)-e0.950.800.700.70C0.350.500.65表1灰色预测精确度检验等级标准一般要求V 20%,最好是V10%,符合要求.关联度检验关联度是用来定量描述各变化过程之间的差别.关联系数越大,说明预测值 和实际值越接近.设X(o)(r) = x(o)(l),x(o)(2),并5)X(o)(r) = x(o)(l),x(o)(2),.,xo)()序列关联

8、系数定义为min |x0) (r) -x(0) (r)| + max |x(0) (r) -x(0) (r)4(r) = |i(0) (r) - x(0) (r) + max |i(0) (r) - x(0) (r)1, r = 0式中,W()7()(U为第/个点X和髀的绝对误差,a)为第,个数据的关联系 数,P称为分辨率,即取定的最大差百分比,O0l,一般取0 = 0.5.X和土的关联度为精度等级关联度均方差比值小误差概率好(1级)0.900350.95合格(2级)0.800.500.80勉强(3级)0.700.650.70不合格(4级)0.65后验差检验后验差检验,即对残差分布的统计特性进

9、行检验.检验步骤如下:1、计算原始时间数列X()=U,9)(2),小)(研的均值和方差仁(U2 tl1 ro,=-yo) , s; /-12、计算残差数列*)=忖,e的均值0和方差s;e=-e(t) ,2n 1n /=I其中dQ) = ()一小)(f) , f = l,2,为残差数歹U3、计算后验差比值C = Sj S4、计算小误差频率P=P 1|e(0) (Z) - e I 0,当Co,当P时.,称模型为小残差概率合格模型.PC模型精度0.950.800.700.65勉强合格0.65不合格表3后验差检验判别参照表1.3 残差GM(1,1)模型当原始数据序列X建立的GM(1,1)模型检验不合格

10、时,可以用GM(1,1)残 差模型来修正.如果原始序列建立的GM(1,1)模型不够精确,也可以用GM(1,1) 残差模型来提高精度.若用原始序列X建立的GM(1,1)模型x(r+l) = 0l)-Kr+-a a可获得生成序列X的预测值,定义残差序列*)/)=/伏)-/(左).若取k=t, 什1,则对应的残差序列为) = e ,e22),.,e()()计算其生成序列e伏),并据此建立相应的GM(1,1)模型e(t + 1) = ?0(1) - I + 殳ae 4得修正模型x( + l)= xo)(l)- e-ak + - + (k-t)(-ae) e(l) 4 eL。L ae其中次ZT)=F为修

11、正参数.0 kl时,模拟精度 低于70%;发展系数-a1.5时,模拟精度低于50%.进一步对预测误差进行考虑,当发展系数-av.3时,1步预测精度在98% 以上,2步和5步预测精度都在90%以上,10步预测精度亦高于80%;当发展系 数- 0.8时,1步预测精度已低于70%.通过以上分析,可得下述结论:1、当-av.3时,GM(1,1)可用于中长期预测;2、当0.3-q0.5时,GM(1,1)可用于短期预测,中长期预测慎用;3、当O.5v-aO.8时,GM(1,1)作短期预测应十分谨慎;4、当0.8-4l时,应采用残差修正GM(1,1)模型;5、当-心1时,不宜采用GM(1,1)模型.1.5

12、GM(U)模型实例分析例:预测学生后悭个学期的成绩.学期1成绩 学期2成绩 学期3成绩 学期4 丽某学生 7974.82574.2976.98 则该学生成绩时间序列如下:x(o)=()(i), )(2), 0)(3), 04) = (79,74.825,74.29,76.98)对X作一次累加后的数列为X=(X,D(2), X(3), ”(4) =(79,153.825,228.115,305.095)对X做紧邻均值生成.令Z6) = 0.5均”(A)+ 0.5Xa)(A-1),得Z =(Z(,)(2),Z(,)(3),Z(1,(4) = (116.4125,151.47,150.1925)则数据矩阵B及数据向量y为-Z(2) 11 -l 16.4125 10j(2)- 74.825B= -Z 1 = -151.471 , Y= x(o)(3) = 74.29-Z(4) 1-150.1925 1J0,(4) |_76.98对参数列3 = ,勿进行最小二乘估计,得a = (B B)i BY = BY= =ayur-0.0144即 = ).0

展开阅读全文
相关资源
猜你喜欢
相关搜索

当前位置:首页 > 论文 > 大学论文

copyright@ 2008-2023 yzwku网站版权所有

经营许可证编号:宁ICP备2022001189号-2

本站为文档C2C交易模式,即用户上传的文档直接被用户下载,本站只是中间服务平台,本站所有文档下载所得的收益归上传人(含作者)所有。装配图网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。若文档所含内容侵犯了您的版权或隐私,请立即通知装配图网,我们立即给予删除!