f_lsim.m

控制系统计算机辅助设计——MATLAB语言与应用(源代码）

function y=f_lsim(a,eta,b,gamma,u,t)
ueq=0; for i=1:length(b), ueq=ueq+b(i)*glfdiff(u,t,gamma(i)); end
h=t(2)-t(1); D=sum(a./[h.^eta]); W=[]; y=zeros(length(t),1); w=y;
for i=1:length(a)
    w=1; for j=2:length(t), w(j)=w(j-1)*(1-(eta(i)+1)/(j-1)); end, W=[W; w];
end
for i=2:length(t)
   for j=1:length(a), A(j)=W(j,2:i)*[y(i-1:-1:1)]; end
   y(i)=(ueq(i)-sum(A.*a./[h.^eta]))/D;
end