krotlica.m

H.S.Park, D-S.Kim, and W.H.Kwon "A Scheduling Method for Network-based Control Systems", IEEE Transa

A=[0 1;0 -0.7]
B=[0;1]
C=[1 0]
K=36

STOP=0;

while STOP==0,

T=input('Input Sampling Time = ');

[P,G]=c2d(A,B,T)

H0=[P -G*K;C*P -C*G*K]
H1=[P -G*K;C 0]

echo on all;

abs(eig(H0))
abs(eig(H1))

pause;

abs(eig(H0*H1))
abs(eig(H1*H0))

echo off;

pause;

Htemp=H0;
for i=1:10,
  Htemp=Htemp*H1;
  eee=abs(eig(Htemp));
  if max(eee) > 1,
    echo on all;
    % H0*H1*H1*H1 .....
    echo off
    i
    eee
    pause;
    break;
  end
end

Htemp=H0*H1*H0;
for i=1:10,
  Htemp=Htemp*H1;
  eee=abs(eig(Htemp));
  if max(eee) > 1,
    echo on all;
    % H0*H1*H0*H1*H1*H1 .....
    echo off
    i
    eee
    pause;
    break;
  end
end

Htemp=input('Input the combination of H0*H1*H0 .... ');

for i=1:50,
  Htemp=Htemp*H1;
  eee=abs(eig(Htemp));
  if max(eee) > 1,
    echo on all;
    % H0*H1*H0*H1*H1*H1 .....
    echo off
    i
    eee
    pause;
    break;
  end
end

keyboard;

end;