灰色预测代码

clear 
syms a b
c=[a,b]';
A=[89677,99215,109655,120333,135823,159878,182321,209407,246619,300670];
B=cumsum(A);           %原始数据累加
n=length(B);
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);c
E
%预测后续数据
F=[];
F(1)=A(1);
for i=2:(n+10)
    F(i)=(A(1)-b/a)/exp(a*(i-1))+b/a;
end
G=[];G(1)=A(1);
for i=2:(n+10)
    G(i)=F(i)-F(i-1);                           %得到预测出来的数据
end
t1=1999:2008;
t2=1999:2018;
G
plot(t1,A,'o',t2,G)                              %原始数据和预测的数据做比较





%结果：：：：；

G =

  1.0e+006 *

  Columns 1 through 18

    0.0897    0.0893    0.1034    0.1196    0.1385    0.1602    0.1854    0.2146    0.2483    0.2873    0.3325    0.3847    0.4452    0.5152    0.5962    0.6899    0.7984    0.9239

  Columns 19 through 20

    1.0691    1.2371