大神帮帮忙。。。新手不懂。

clear;
syms a b;
c=[a b]';
A=[170 184 210 210 200];
B=cumsum(A);
n=length(A);
for i=1:(n-1)
C(i)=(B(i)+B(i+1)/2);
end
D=A;D(1)=[];
D=D';
E=[-C;ones(1,n-1)];
c=inv(E*E')*E*D;
c=c';
a=c(1);b=c(2);
F=[];F(1)=A(1);
for i=2:(n+5)
G(i)=F(i)-F(i-1);
end
t1=1:5;
t2=1:10;
G
plot(t1,A,'o',t2,G)