〈灰色系统理论及其应用〉第7版，刘思峰 P153例7.3.1

%〈灰色系统理论及其应用〉第7版，刘思峰 P153例7.3.1

clear all;clc;close all
X0=[6,20,40,25,40,45,35,21,14,18,15.5,17,15];
X1=cumsum(X0);
Z1=[0];
for i=2:13
    Z1=[Z1,(X1(i)+X1(i-1))/2];
end
a=(sum(Z1(2:13))*sum(X0(2:13))-(13-1)*sum(Z1(2:13).*X0(2:13)))/((13-1)*sum(Z1(2:13).^2)-sum(Z1(2:13))^2) %公式见P165
b=( sum(X0(2:13))*sum(Z1(2:13).^2) - sum(Z1(2:13))* sum(Z1(2:13).*X0(2:13))) /((13-1)*sum(Z1(2:13).^2)-sum(Z1(2:13))^2)
disp(num2str(b/a))

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