GM(1,1)

[meatball1982]

function [out]=gm(x0)
n=length(x0);
%由数列x0生成数列x1;
s=0;
for i=1:n
s=s+x0(i);
x1(i)=s;
end
%由数列x1生成矩阵B和Y；
for i=1n-1);
Q(i)=-(x1(i)+x1(i+1))/2;
end
M=ones(n-1,1);
B=[Q' M];
for i=1n-1);
C(i)=x0(i+1);
end
Y=C';
beta=inv(B'*B)*B'*Y;
%计算出a和u;
a=beta(1);
u=beta(2);
out=(x0(1)-u/a)*(1-exp(a))*exp(-a*3)

这是我写的一个GM11的函数程序，定义了一个输出。你的函数没有输出。

function [pre]=fun_GM11(X0,pre_num)
X1=cumsum(X0)
n=length(X0);
G(:,1)=-(X1(2:n)+X1(1:n-1))/2;
G(:,2)=ones(n-1,1);
Y(:,1)=X0(2:n);
pre=Y;
a1=inv(G'*G)*G'*Y;
a=a1(1);
u=a1(2);
X_pre(1)=X0(1);
for k=1:n-1+pre_num
X_pre(k+1)=(1-exp(a))*(X0(1)-u/a)*exp(-a*k);
end
pre=X_pre;