pidsmir2.m

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

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(1,1)=u(2,1);
u(2,1)=u1(2,1);
u(1,2)=u(2,2);
u(2,2)=u1(2,2);

s(2)=s(1);
s(3)=y(3);