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）;???????????????????????%構(gòu)造累加矩陣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?