%输入原始序列X0
X0=[4747 6162 7375 10661 11648 13224 22310];
%对原始序列X0一次累加生成序列X1
[n,m]=size(X0);
for i=1:n   c=0;
   for j=1:i
       c=X0(j,m)+c;
   end
   X1(i,m)=c;
end
disp('X1=');
disp(X1);
%计算数据矩阵B
for i=2:n
     Z1(i,m)=(X1(i-1,m)+X1(i,m))/2;   
end
disp('Z1=');
disp(Z1);
B=[-Z1 ones(n,1)];
B(1,:)=[];
disp('B=');
disp(B);
%计算数据矩阵Y及a,b的值
X0(1,:)=[];
Y=X0;
disp('Y=');
disp(Y);
Q=inv(B'*B)*B'*Y;
disp('Q=');
disp(Q);
X1=
   1.0e+004 *
     0.0011          0          0          0          0     2.2310     2.2310
     0.0024          0          0          0          0          0          0
     0.0039          0          0          0          0          0          0
     0.0058          0          0          0          0          0          0
     0.0082          0          0          0          0          0          0
     0.0122          0          0          0          0          0          0
     0.0168          0          0          0          0          0          0
Z1=
          0
    17.8100
    31.8150
    48.7600
    69.8750
   101.6600
   144.6600
??? Error using ==> horzcat
CAT arguments dimensions are not consistent.

小弟菜鸟....有急用啊！！！各位高手帮忙！！！