%〈灰色系统理论及其应用〉第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 是灰色作用量。