matlab求解－a（发展系数）的方法。

<灰色系统理论及应用>刘思峰第7版，P165~166，求解－a（发展系数）的方法。b:灰色作用量。

clear all;clc;close all

a=-0.1;
X0=[];
for i=0:5
X0= [X0,exp(-a*i)];
end
X1=cumsum(X0);
Z1=[0];
for i=2:6
    Z1=[Z1,(X1(i)+X1(i-1))/2];
end

a=(sum(Z1(2:6))*sum(X0(2:6))-(6-1)*sum(Z1(2:6).*X0(2:6)))/((6-1)*sum(Z1(2:6).^2)-sum(Z1(2:6))^2);

b=( sum(X0(2:6))*sum(Z1(2:6).^2) - sum(Z1(2:6))* sum(Z1(2:6).*X0(2:6))) /((6-1)*sum(Z1(2:6).^2)-sum(Z1(2:6))^2)​

求解出来：a是个正值 ，前面加个负号，变成－a（发展系数）。