clear
x10=
x20=
x30=
x40=
%生成列
x11=cumsum(x10)
x21=cumsum(x20)
x31=cumsum(x30)
x41=cumsum(x40)
%求参数，建立模型
for k=2:length(x11)
    z1=(1/2)*(x11(k)+x11(k-1))
end
z1
B=[-(z1(2:end))',(x21(2:end))',(x31(2:end))',(x41(2:end))'];
Y=(x10(2:end))';
C=inv(B'*B)*B'*Y
自己根据书编出的代码，用matlab做出来的结果是
C =

   -0.0072
    0.1476
   -0.1282
这个是GM（1，4）不是应该有四个参数结果出来么？为什么我的是三个？？？？还有书上给的最后的公式是i
可是，我做的分析的参考论文里给出的最后的动态模型是

我用上面的结果尅直接代入下面这个模型中么？

问题整理，1.为什么我的结果出现了三个系数，而不是四个？是哪里出问题了？2.可以如果有四个系数的结果直接代入最下面的模型里么？下面的模型便于分析和简洁！