灰色预测-MATLAB程序

预测2013年到2017年措施实施的最坏结果的无偏灰色预测程序
%GM(1,1) MATLAB 程序 最坏预测结果
clc
clear
syms a b;
A=[202 219 221 325 347 333 360 407 441 476 479 458]; %原始数据
B=cumsum(A); %生成序列
n=length(A);
for i=1:(n-1)
C(i)=(B(i)+B(i+1))/2;
end
D=A;
D(1)=[];
D=D'; %数据向量
E=[-C;ones(1,n-1)]; %数据矩阵
c=inv(E*E')*E*D;
c=c';
a=c(1); %发展系数
b=c(2); %灰作用量
a1=log((2-a)/(2+a));
A1=2*b/(2+a);
G=[];
G(1)=A(1);
for i=2:(n+6)
G(i)=A1*exp(a1*(i-1)); %原始数据的估计值
end
t1=2000:2011;
t2=2000:2017;
plot(t1,A,'o',t2,G)
for i=1:n
d(i)=abs(A(i)-G(i)); %绝对误差
o(i)=d(i)/A(i); %相对误差
end
disp(['绝对残差： ',num2str(d(1:n))])
disp(['相对残差： ',num2str(o(1:n))])<