www.pudn.com > controlsystem.ZIP > EX3005.M

% Example 3.5 
% 
for alph=[-1:1] 
 alph 
 num=[1,alph]; 
 den=[1 10 27 18]; 
 [a,b,c,d]=tf2ss(num,den) 
 cam=ctrb(a,b); 
 rcam=rank(cam) 
 oam=obsv(a,c); 
 coam=rank(oam) 
end 
 
 
% The result 
 
alph = 
    -1 
a = 
   -10   -27   -18 
     1     0     0 
     0     1     0 
b = 
     1 
     0 
     0 
c = 
     0     1    -1 
d = 
     0 
rcam = 
     3 
coam = 
     3 
 
 
alph = 
     0 
a = 
   -10   -27   -18 
     1     0     0 
     0     1     0 
b = 
     1 
     0 
     0 
c = 
     0     1     0 
d = 
     0 
rcam = 
     3 
coam = 
     3 
 
 
alph = 
     1 
a = 
   -10   -27   -18 
     1     0     0 
     0     1     0 
b = 
     1 
     0 
     0 
c = 
     0     1     1 
d = 
     0 
rcam = 
     3 
coam = 
     2 
»