gm(1,1).m

来自「灰色系统理论预测 灰色系统理论预测」· M 代码 · 共 26 行

clear
syms a b;
c=[a b]';
A=[174 179 183 189 207 234 220.5 256 270 285];
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);u=c(2);
F=[];F(1)=A(1);
for i=2:(n+10)
    F(i)=(A(1)-u/a)/exp(a*(i-1))+u/a ;
end
G=[];G(1)=A(1);
for i=2:(n+10)
    G(i)=F(i)-F(i-1);
end
x1=1995:1:2004;x2=2005:1:2014;
plot(x1,G(1:10));hold on;plot(x2,G(11:20));
xlabel('时间(年)'),ylabel('排污量(亿吨)'),title('近10年（1995-2004）的模拟情况和未来10年（2005-2014）污染趋势图');