灰色预测模型matlab代码

clear
syms a b;
c=[a b]';              %待定系数
m=xlsread('1.xls');     %取Excel中的数据
d=m(:,3)                %取数据中的列向量
A=d.';
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);
F=[];F(1)=A(1);
for i=2:(n+1)
  F(i)=(A(1)-b/a)/exp(a*(i-1))+b/a;       %求出X（t+1）1
end
G=[];G(1)=A(1);
for i=2:(n+1)
G(i)=F(i)-F(i-1);                            %X（t+1）0
end
for i=1:(n-1)
   o=G(i)-A(i);
end
 err=sqrt(o^2/(n-1))
t1=2008:2015;
t2=2008:2016;
G
a
b