资源描述：
《灰色系统预测GM(1,1)模型及其Matlab实现.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、灰色系统预测GM(1,1)模型及其Matlab实现预备知识(1)灰色系统白色系统是指系统内部特征是完全已知的；黑色系统是指系统内部信息完全未知的；而灰色系统是介于白色系统和黑色系统之间的一种系统，灰色系统其内部一部分信息已知，另一部分信息未知或不确定。(2)灰色预测灰色预测，是指对系统行为特征值的发展变化进行的预测，对既含有已知信息又含有不确定信息的系统进行的预测，也就是对在一定范围内变化的、与时间序列有关的灰过程进行预测。尽管灰过程中所显示的现象是随机的、杂乱无章的，但毕竟是有序的、有界的，因此得到的数据集合具备潜在的规律。灰色预测是利用这种规律建立灰色模型对灰色系统进行预测
2、。目前使用最广泛的灰色预测模型就是关于数列预测的一个变量、一阶微分的GM(1,1)模型。它是基于随机的原始时间序列，经按时间累加后所形成的新的时间序列呈现的规律可用一阶线性微分方程的解来逼近。经证明，经一阶线性微分方程的解逼近所揭示的原始时间序列呈指数变化规律。因此，当原始时间序列隐含着指数变化规律时，灰色模型GM(1,1)的预测是非常成功的。1灰色系统的模型GM(1,1)1.1GM(1,1)的一般形式设有变量X(0)＝{X(0)(i)，i=1,2，...，n}为某一预测对象的非负单调原始数据列，为建立灰色预测模型：首先对X(0)进行一次累加(1—AGO,AcumulatedG
3、eneratingOperator)生成一次累加序列： 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,可用下式求解
4、:＝(BTB)-1BTYn(5)式中:B—数据阵；Yn—数据列B＝(6)Yn＝(X(0)(2),X(0)(3),…,X(0)(ｎ))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)(ｋ)(0)(ｋ)＝(1)(k)－(0)(i)因为(1)(k-1)＝(0)(i)，所以(0)(ｋ)＝(1)(k)－(1)(k-1)。2应用举例取某高校1998年～2
5、003年的某专业招生数据建模,见表1。表1某高校专业招生数据表年招生人数20001322001922002118200313020041872005207以表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,6
6、59,866}和数据阵B、数据列YnB＝，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还原值实际值误差拟合相对误差(％)2000132132132132002001225.060879622493921－12002339.295441834211411843.382003
7、479.52123472140130－10－7.692004651.65196591721871582005862.9466129866211207－4－1.9320061122.316167259252－7－2.78由表2知预测精度较高。2006年某专业招生人数预测值为259人。由于人数为整数，所以结果取整数部分。GM(1,1)是一种长期预测模型，在没有大的市场波动及政策性变化的前提下，该预测值应是可信的。如前所述，影响招生人数的因素很多且难以预测。因此，在采用灰色系统理论进行定量预测时，
