灰色预测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