灰色预测MATLAB对比代码,灰色马尔可夫模型代码

%%发展系数a灰色作用量b

syms a b；

%%原始数列e列

a1=xls read (c : (用户(0223 )文档) MATLAB )劳动人口. xlsx )、(‘C1:C11 )；

a=a1’；

A=A.*0.01；

n=长度(a；

%%累计一次

b(1)=a ) 1；

for i=2 : n

b(I )=b ) I-1 ) a ) I；

结束

%%平均生成

for i=2 : n

z(I )=(b ) I ) b(I-1 ) ) *0.5；

结束

z(1)=[]；

%%构建Y B矩阵

Y=A； y(1)=[]； y=y’；

B=[-Z； ONES(1，n-1 )；

%%计算a b参数的值

q=sum(z；

w=sum(a )-a ) 1；

e=0；

for i=2 : n

e=ez(I-1 ) a ) I；

结束

r=0；

for i=2 : n

r=rz(I-1 ) z (I-1 )；

结束

a=() qw ) (n-1 ) e )/() n-1 ) r-) QQ )；

b=(wr )-(QE ) )/(n-1 ) r-) ) QQ )；

%%预测了过去四年的数据

yuce=[]；

YUCE(1)=a )；

forI=2:(n4 ) ) )

YUCE(I )=(a )1)- b/a ) )/exp ) a* (I-1 ) ) b/a；

结束

恢复%%数据

环宇安=[ ]；

环屿安(1)=a )；

forI=2:(n4 ) ) )

环屿安(I )=Yuce(I )-Yuce ) I-1；

结束

%%模型检查

cancha zhi=ABS (0.784-0.789691929 )/0.784