〈灰色系统理论及其应用〉第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\hat a= \cfrac{1}{1+\cfrac{1}{x}}

Σnk=2z(1)(k)x(0)(k) {\underset{k=2} {\overset{n}{\Large \Sigma}}z^{(1)}(k)x^{(0)}(k)}