pidsmir1.m

来自「基于TKGK实验台的实验数据编码」· M 代码 · 共 15 行

function [u,s,a0]=pidsmir1(y,s,a0,u,ac,bc,cc,dc)

if isempty(y);
   y=zeros(1,5);
end

[u1,x]=lsim(ac,bc,cc,dc,[s(2),s(3);s(1),y(3)],[0:1],a0);
a0=x(2,:);
u(2)=u(1);
u(1)=u1(2,1);
s(2)=s(1);
s(3)=y(3);