gm（1，n）.m

来自「灰色系统理论预测 灰色系统理论预测」· M 代码 · 共 27 行

A=[];
[m,n]=size(A);
k=A(1,1);
Y=A(1,:);
Y(1)=[];
for i=1:(n-1)
    E(i)=A(1,i)/2+A(1,i+1)/2;
end
B=cumsum(A,2);
B(1,:)=[];
B(:,1)=[];
C=[-E;B];
D=inv(C*C')*C*Y;
e=D(1);
D(1)=[];
F=D';
for j=1:(n-1)
    G(j)=F*B(:,j);
end
H=[];   H(1)=k;
for j=2:n
    H(j)=(k+1/e*G(j-1))/exp(e*j)+1/e*G(j-1);
end
P(1)=k;
for i=2:n
    P(i)=H(i)-H(i-1);
end