x0=[213.9166189 235.879359 364.8587163 227.6641233 303.7526842 362.882818];
n=6;
for i=1:n
x1(i)=sum(x0(1:i));
end
for k=2:n % 紧邻均生成 z
z(k)=0.5*x1(k)+0.5*x1(k-1);
end
for i=1:n-1
b(i,1)=-z(i+1);
y(i)=x0(i+1);
end
b(:,2)=1;y=y'; % 转置为列向量
au=b\y; % 作矩阵除法,计算a uyc1(1)=x0(1);
for k=1:n
c=x0(1)-au(2)/au(1);
yc1(k+1)=c*exp( -au(1)*k)+au(2)/au(1);
end
yc0(1)=x0(1);
for k=1:n
yc0(k+1)=yc1(k+1)-yc1(k);
end
disp(uint16(yc0(2:1:n+1)))
显示结果
??? Error using ==> mldivide
Matrix dimensions must agree.
Error in ==> Untitled at 14
au=b\y; % 作矩阵除法,计算a uyc1(1)=x0(1);