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（發展係數）。