GM(1,1).m %建立符号变量a(发展系数)和b(灰作用量) syms a b; c = [a b]'; %原始数列 A A = [174, 179, 183, 189, 207, 234, 220.5, 256, 270, 285];%填入已有的数据列! n = length(A); %对原始数列 A 做累加得到数列 B B = cumsum(A); %对数列 B 做紧邻均值生成 for i = 2:n C(i) = (B(i) + B(i - 1))/2; end C(1) = [];
function SGrey X0 = input('请输入原始负荷数据:'); %输入原始数据 n = length(X0); %原始n年数据 %累加生成 X1 = zeros(1,n); for i = 1:n if i == 1 X1(1,i) = X0(1,i); else X1(1,i) = X0(1,i) + X1(1,i-1); end end X1 %计算数据矩阵B和数据向量Y B = zeros(n-1,2); Y = zeros(n-1,1); for i = 1:n-1 B