GM11模型matlab程序

function GM=gm11（dataN）
T=length（data）; 
X0=data;
for i=2:T
   X1（1）=X0（1）;
   X1（i）=X1（i-1）+X0（i）;                %用AGO生成一阶累加生成模块
end
for i=1:T-1
   M（i）=-（0.5*（X1（i）+X1（i+1）））;
end
B=zeros（T-12）;                       %构造累加矩阵B
for i=1:T-1
     for j=1:2
         if j<2
            B（ij）=M（i）;
         elseif j>1
            B（ij）=1;
         end
     end
end
for i=2:T