〈灰色系統理論及其應用〉第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 是灰色作用量。

a^=11+1x

Σnk=2z(1)(k)x(0)(k)