main1.m

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

%main1

[y0,s,filename]=indata;
%if filename=='pidsmir1';
if strcmp(filename,'pidsmir1')
[ac,bc,cc,dc]=linmod2('smith1c');
[b,a0]=smith1c([],[],[],0);
else
  a0=[0 0 0 0 0 0];
end

i=1;m=500;
u=[0 0];
%v=140;
s0=[0 0 1];
while i<m
   t0=now;
   [y,s,filename]=indata;
   if isempty(y)
      y=[0 0 0 0 0];
   end;
     s0(1)=s;
    if strcmp(filename,'pidsmir1')
      fm=[filename, '(y,s0,a0,u,ac,bc,cc,dc)'];
    else
      fm=[filename, '(y,s0,a0,u)'];
   end
   
   
   [u,s0,a0]=eval(fm);
   u
outdata(u(2));
%i=i+1;
t01=etime(datevec(now),datevec(t0));
 pause(abs(1-t01));  %Running time is 1 sec. for a cycle 
 end
fclose(f11);