- 1、本文档共6页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
灰色系统预测GM(1,1)模型及其Matlab实现
预备知识
(1)灰色系统
白色系统是指系统内部特征是完全的;黑色系统是指系统内部信息完全未知的;而灰色系统是介于白色系统和黑色系统之间的一种系统,灰色系统其内部一局部信息,另一局部信息未知或不确定。
(2)灰色预测
灰色预测,是指对系统行为特征值的开展变化进行的预测,对既含有信息又含有不确定信息的系统进行的预测,也就是对在一定范围内变化的、与时间序列有关的灰过程进行
预测。尽管灰过程中所显示的现象是随机的、杂乱无章的,但毕竟是有序的、有界的,因此得到的数据集合具备潜在的规律。灰色预测是利用这种规律建立灰色模型对灰色系统进行预测。
目前使用最广泛的灰色预测模型就是关于数列预测的一个变量、一阶微分的GM(1,1)模型。它是基于随机的原始时间序列,经按时间累加后所形成的新的时间序列呈现的规律可用一阶线性微分方程的解来逼近。经证明,经一阶线性微分方程的解逼近所揭示的原始时间序列呈指数变化规律。因此,当原始时间序列隐含着指数变化规律时,灰色模型GM(1,1)的预测是非常成功的。
1灰色系统的模型GM(1,1)
1.1GM(1
设有变量X(0)={X(0)(i),i=1,2,...,n}为某一预测对象的非负单调原始数据列,为建立灰色预测模型:首先对X(0)进行一次累加(1—AGO,AcumulatedGeneratingOperator)生成一次累加序列:
X(1)={X(1)(k),k=1,2,…,n}
其中
X(1)(k)=X(0)(i)
=X(1)(k-1)+X(0)(k)(1)
对X(1)可建立下述白化形式的微分方程:
十=u(2)
即GM(1,1)模型。
上述白化微分方程的解为(离散响应):
(1)(k+1)=(X(0)(1)-)+(3)
或
(1)(k)=(X(0)(1)-)+(4)
式中:k为时间序列,可取年、季或月。
1.2辩识算法
记参数序列为,=[a,u]T,可用下式求解:
=(BTB)-1BTYn(5)
式中:B—数据阵;Yn—数据列
B=(6)
Yn=(X(0)(2),X(0)(3),…,X(0)(n))T〔7〕
1.3预测值的复原
由于GM模型得到的是一次累加量,k{n+1,n+2,…}时刻的预测值,必须将GM模型所得数据(1)(k+1)(或(1)(k))经过逆生成,即累减生成(I—AGO)复原为(0)(k+1)(或(0)(k))才能用。
(1)(k)=(0)(i)
=(0)(i)+(0)(k)
(0)(k)=(1)(k)-(0)(i)
因为(1)(k-1)=(0)(i),所以(0)(k)=(1)(k)-(1)(k-1)。
2应用举例
取某高校1998年~2003年的某专业招生数据建模,见表1。
表1某高校专业招生数据表
年
招生人数
2000
132
2001
92
2002
118
2003
130
2004
187
2005
207
以表1中的数据构造原始数据列X(0),即
X(0)={X(0)(1),X(0)(2),X(0)(3),X(0)(4),X(0)(5),X(0)(6)}
={132,92,118,130,187,207}
对X(0)进行一次累加(1—AGO),生成数列:
X(1)(k)=X(0)(i)即
X(1)={X(1)(1),X(1)(2),X(1)(3),X(1)(4),X(1)(5),X(1)(6)}
={132,224,342,472,659,866}
和数据阵B、数据列Yn
B=,Yn=(92,118,130,187,207)T
由式(5)得
=[a,u]T
=
由式(4)得灰色预测模型GM(1,1)为
(1)(k)=(X(0)(1)-)+
=(132+277.0137483)-277.0137483
=409.0137483-277.0137483
预测值及预测精度见表2。
表2某高校专业招生预测值及预测精度表
年
GM(1,1)模型
计算值
1—AGO
还
原
值
实
际
值
误
差
拟合相对
误差(%)
2000
132
132
132
132
0
0
2001
225.0608796
224
93
92
1
-1
2002
339.2954418
342
114
118
4
3.38
2003
479.52123
472
140
130
-10
-7.69
2004
651.6519
659
172
187
15
8
2005
862.946
您可能关注的文档
- 三年级分数加减法.ppt
- 下8.10(1)第一框正义是人类良知的“声音”.ppt
- 专有名词-连词研学稿.ppt
- 猕猴桃试验报告.doc
- 环境艺术设计.ICAD模拟试题和答案.doc
- 万有引力与航天复习课件.ppt
- 牛津版七年级英语上册:-unit-5同步测试题(二).doc
- 点击查看下载-妇产科教案-第十八章.doc
- 物理第十二章和第十三章测试卷.doc
- 煤与瓦斯共采技术现状及展望.doc
- 10《那一年,面包飘香》教案.docx
- 13 花钟 教学设计-2023-2024学年三年级下册语文统编版.docx
- 2024-2025学年中职学校心理健康教育与霸凌预防的设计.docx
- 2024-2025学年中职生反思与行动的反霸凌教学设计.docx
- 2023-2024学年人教版小学数学一年级上册5.docx
- 4.1.1 线段、射线、直线 教学设计 2024-2025学年北师大版七年级数学上册.docx
- 川教版(2024)三年级上册 2.2在线导航选路线 教案.docx
- Unit 8 Dolls (教学设计)-2024-2025学年译林版(三起)英语四年级上册.docx
- 高一上学期体育与健康人教版 “贪吃蛇”耐久跑 教案.docx
- 第1课时 亿以内数的认识(教学设计)-2024-2025学年四年级上册数学人教版.docx
文档评论(0)