gm（11）模型的matlab程序

function gm（x0）;
clc
format long;
if length（x0（:1））==1    %对输入矩阵进行判断，如不是一维列矩阵，进行转置变换
    x0=x0‘;
end
n=length（x0）;
%由数列x0生成数列x1;
s=0;
for i=1:n
    s=s+x0（i）;
    x1（i）=s;
end
%由数列x1生成矩阵C和A； 
for i=1:（n-1）; 
A（i）=-（x1（i）+x1（i+1））/2; 
end 
A=[A‘ ones（n-11）]; 
for i=1:（n-1）; 
C（i）=x0（i+1）; 
end 
C=C‘;
beta=inv（A‘*A）*A‘*C;
%计算出a和u;
a=beta（1）;