数学建模+灰色预测模型+MATLAB

1、文档大全12.5灰色预测我们通常所说的系统是指：由客观世界中相同或相似的事物和因素按一定的秩序相互关联、相互制约而构成的一个整体例如：工程技术系统、社会系统、经济系统等如果一个系统中具有充足的信息量，其发展变化的规律明显、定量描述方便、结构与参数具体，则这种系统通常称为白色系统如果一个系统的内部特征全部是未知的，则称此系统为黑色系统如果系统内部信息和特征是部分已知的，另一部分是未知的，这种系统称为灰色系统例如：社会系统、农业系统、经济系统、气象系统、生物系统等对于这类系统，内部因素难以辨识，相互之间的关系较为隐蔽，人们难以准确了解这类系统的行为特征因此，对于这类问题进行定量描述，即建立模型难度

2、较大区别白色系统与灰色系统的重要标志是系统内各因素之间是否具有确定的关系灰色系统分析方法主要是根据具体灰色系统的行为特征数据，充分利用数量不多的数据和信息寻求相关因素自身与各因素之间的数学关系，建立相应的数学模型目前，灰色系统理论在实际中已得到了广泛的应用，例如：在工程技术、经济管理、气象预报以及政治、社会、工业、农业等领域都取得了一定的应用成果我们往往要对农业问题、商业问题等做未来的预测工作，另外，进行军事战争以及治理生态环境也需对未来的发展情形做一可靠的分析，这就产生了灰色预测灰色预测是对灰色系统问题进行未来的预测，实际问题中，应用最多的灰色预测模型是以GM(1,1)(即GM(1,N)当N

3、=1时的特例)模型为基础的.12.5.1GM(1,1)模型的建立设X(o)=(X(o)(l),X(o)(2),X(o)(n),做1AGO,得X=(X(1)(1),X(1)(2),X(1)(n)=(X(1)(1),X(1)(1)+X(0)(2),X(1)(n-1)+X(0)(n)TOC o 1-5 h z则GM(1,1)模型相应的微分方程为：dX(1)+aX(1)=u(1)dt式中：a称为发展灰数；“称为内生控制灰数.设&=(a,)T,按最小二乘法得到a=(BtB)-1BtY(2)1其中-(X(i)(l)+X(1)(2)2-(X(1)(2)+X(1)(3)2-丄(X(1)(n-1)+X(1)(n

4、)1V2丿厂X(0)(2)、X(0)(3)1:VX(0)(n)丿易求得，方程(1)的解为X(1)(k+1)二(X(0)(1)-u)e-ak+-(3)aa例4100m成绩预测19831990年世界男子和中国女子100m最好成绩如表6.表6各年度最好成绩年份19831984198519861987198819891990男子9.939.969.989.959.939.929.949.93女子11.9511.6611.6311.6511.3511.3211.5811.32记世界男子100m成绩的原始数列为X(0)=(9.93,9.96,9.9&9.95,9.93,9.92,9.94,9.93)建立G

5、M(1,1)模型，即按式(1)、(2)、(3)得到预测模型为X(1)(k+1)=(9.93-13884.61)e-0.0007185266k+13884.61由预测模型得预测值为年份模型预测值/s19919.9219929.9120009.85记中国女子的原始数列为X(0)=(11.95,11.66,11.63,11.65,11.35,11.32,11.58,11.32)同样建立GM(1,1)模型，得到预测模型为X(1)(k+1)=(11.95一2602.187)e-0.00451067k+2602.187从而得到中国女子100m成绩的预测值年份模型预测值/s199111.30199211.2

6、4200010.8512.5.2模型检验灰色预测检验一般有残差检验、关联度检验和后验差检验(1)残差检验按预测模型计算乂(i)(i)，并将文(i)(i)累减生成文(o)(i)，然后计算原始序列X(o)(i)与X(O)(i)的绝对误差序列及相对误差序列.A(0)(i)=1X(0)(i)-X(0)(i)|i=1,2,n=號x100%i=卩，n(2)关联度检验定义1选取参考数列X=X(k)|k=1,2,n=(X(1),X(2),X(n)00000i=1,2,mTOC o 1-5 h z其中k表示时刻.假设有m个比较数列X=X(k)|k=1,2,n=(X(1),X(2),X(n)iiiii则称minm

7、inIX(k)一X(k)I+pmaxmaxIX(k)一X(k)Ig(k)=ik0iik0i(1)iIX(k)一X(k)I+pmaxmaxIX(k)一X(k)I0iik0i为比较数列X.对参考数列X0在k时刻的关联系数，其中pe0,1为分辨系数，一般取i0p=0.5.称式(1)中minminiX0(k)-Xi(k)Imaxmax丨X0(k)-Xi(k)I分别为两级最小差和两ikik级最大差由(1)式易看出，p越大，分辨率越大；p越小，分辨率越小.式(1)定义的关联系数是描述比较数列与参考数列在某时刻关联程度的一种指标，由于各个时刻都有一个关联数，因此信息显得过于分散，不便于比较，为此我们给出以下

8、定义定义2称2)r=g(k)inik=1为数列X,对参考数列X0的关联度.i0由式(2)易看出，关联度是把各个时刻的关联系数集中为一个平均值，也就是把过于分散的信息集中处理.根据前面所述关联度计算方法计算出文(0)(i)与原始序列X(0)(i)的关联系数，然后计算出关联度，根据经验，当j=0.5时，关联度大于0.6便满意了.3)后验差检验实用标准实用标准文档大全文档大全计算原始序列标准差：巨X(o)(i)X(0)2n1计算绝对误差序列的标准差：工A(0)(i)A(0)2计算方差比：计算小误差概率：P=plA(0)(i)A(0)l0.6745S1令e=lA(0)(i)A(0)l,S=0.6745

9、S,则P=pe0.950.800.700.65勉强合格0.65不合格若残差检验、关联度检验和后验差检验都能通过，则可以用所建模型进行预测；若用原始时间序列X建立的GM(1,1)模型检验不合格或精度不理想时，这时要对建立的GM(1,1)模型进行修正或提高模型的预测精度其修正方法如下：设原始时间序列X(0)建立的GM(1,1)模型为X(1)(i+1)二(X(0)(1)u)e-ai+-aa可获得生成序列X(1)的预测值幺，即对于X(1)=X(1)(1),X(1)(2),X(1)(n)，有预测序列X=X(1)(1),X，,X(1)(n)，定义残差为e(0)(j)=X(i)(j)-文(i)(j)若取j=

10、i,i+1,n,则与X(i)及文对应的残差序列为e(0)=e(0)(i),e(o)(i+1),e(o)(n)为便于计算上式改写为e(0)=e(o)(l),e(o)(2),e(0)(n)e(0)的累加生成序列为e(1)=e(1)(1),e(1)(2),e(1)(n)n=n-ie(1)可建立相应的GM(1,1)模型:e(1)(k+1)=(e(0)(1)-Ue)e-ak+Ueaaeeue(k+1)的导数e(k+1)=(-a)(e(o)(l)小匕d加上e(k+1)修正X(1)(k+1)，得eae修正模型：X(1)(k+1)=(X(0)(1)u)e-ak-+6(k1)(a)(e(o)(1)鋼e-ae(k

11、1)aaeae其中6(k1)=,:2为修正系数0,k0.57的所以关联度检验通过第六步，后验差检验：(1)计算：-1X(0=(386.(+6476+576+9.35+873.89+1085.73932.2：552.33)6S=1工%(0)(k)一X如2=341.0652)计算残差的均值：A=1(0+53.87+26.28+69.82+95.37+33.48)=46.47残差的标准差：工A(k)-A2s338438S=33.8438,C=-2=0.099232n1S341.0651S二0.6745S二0.6745x341.065二230.04801e=1A(k)-AI二46.47,7.4,20.

12、19,23.35,48.9,12.98k所有。点都小于S0,故戶=心0.35.所以后验差检验通过第七步，残差检验：(1)计算X(1)(1)=386.06,X(1)(2)=916.498734,X(1)(3)=1569.574052,X(1)(4)=2373.639335,X(1)(5)=3363.603222,X(1)(6)=4582.44517(2)累减生成序列：X(0)(1)=386.06,X(0)(2)=530.44,X(0)(3)=653.07X(0)(4)=804.07,X(0)(5)=989.96,X(0)(6)=1218.35计算绝对误差序列及相对误差序列：绝对误差序列(0=0,

13、53.87,26.28,69.82,95.37,33.48相对误差序列=0,11.3%,3.87%,7.99%,8.79%,2.67%相对误差序列中有的相对误差很大，所以要对原模型进行残差修正以提高精度利用残差对原模型进行修正：取e(0)=53.87,26.28,69.82,95.37,33.48e(1)=53.87,80.15,149.94,245.34,278.82得e(k+1)=848.1722e0.0589k一794.3022最后得修正模型为：X(1)(k+1)=2294.301108e0.207988k一1908.24118+5(k一1)49.9573e0.0589(k-1)其中5(k-1)=表9修正后的残差计算表序号kX(1)(k)x(1Xk)修正后误差相对误差1386.06386.06002862.63900.0137.384.33%31541.981622.5680.585.23%42451.872426.6310.760.45%53501.203419.20822.34%64753.334622