灰色预测模型matlab预测20个数据 怎么matlab灰色模型预测这组数据的下一...

clc

clearall

%本程序主要用来计算根据灰色理论建立的模型的预测值。

%应用的数学模型是GM(1,1)。

%原始数据的处理方法是一次累加法。

y=[1662.872163.41965.352472.482900.663034.932755.532073462];%已知数据

n=length(y);

yy=ones(n,1);

yy(1)=y(1);

fori=2:n

yy(i)=yy(i-1)+y(i);

end

B=ones(n-1,2);

fori=1:(n-1)

B(i,1)=-(yy(i)+yy(i+1))/2;

B(i,2)=1;

end

BT=B';

forj=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=4;%需要预测个数

i=1:t_test+n;

yys(i+1)=(y(1)-t).*exp(-a.*i)+t;

yys(1)=y(1);

forj=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;

fori=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))]);

输出结果：

百分绝对误差为：228.3113%

预测值为：3710.1523978.21424265.64424573.8413