graymodel.m

GM(1,1)灰色预测的matlab的.m文件

y=[-0.0682
-0.1014
-0.2746
0.1889
-0.3378
0.0077
-0.3535
0.0991
-0.0444
0.0229
0.152
0.4003
-0.0273
0.3693
-0.0799
-0.0556
0.1679
0.2836
0.429
0.244
0.5578
0.418
0.1376
0.2311
0.3712
0.6631
0.2187
0.5713
0.8368
0.6998
0.5325
0.7552
0.8597
0.7773
0.7296
0.9798
0.84
1.0162
0.8029
];
n=length(y);
yy=ones(n,1);
yy(1)=y(1);
for i=2:n
    yy(i)=yy(i-1)+y(i);
end
B=ones(n-1,2);
for i=1:(n-1)
    B(i,1)=-(yy(i)+yy(i+1))/2;
    B(i,2)=1;
end
BT=B';
for j=1:n-1
    YN(j)=y(j+1);
end
YN=YN';
A=inv(BT*B)*BT*YN;
a=A(1);
u=A(2);
t=u/a;
t_test=input('请输入需要预测个数：');
i=1:t_test+n;
yys(i+1)=(y(1)-t).*exp(-a.*i)+t;
yys(1)=y(1);
for j=n+t_test:-1:2
    ys(j)=yys(j)-yys(j-1);
end
x=1:n;
xs=2:n+t_test;
yn=ys(2:n+t_test);
plot(x,y,'^r',xs,yn,'*-b');
det=0;
for i=2:n
    det=det+abs(yn(i)-y(i));
end
det=det/(n-1);
disp(['百分绝对误差为：',num2str(det),'%']);
disp(['预测值为： ',num2str(ys(n+1:n+t_test))]);