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

% Example 3.3 
% 
a1=[0,1;-1,-2]; 
b1=[0;1]; 
c1=[1,3];d1=[1]; 
a2=[0,1;-1,-3]; 
b2=[0;1]; 
c2=[1,4];d2=0; 
[a,b,c,d]=series(a1,b1,c1,d1,a2,b2,c2,d2) 
[a,b,c,d]=parallel(a1,b1,c1,d1,a2,b2,c2,d2) 
[a,b,c,d]=feedback(a1,b1,c1,d1,a2,b2,c2,d2) 
[a,b,c,d]=feedback(a1,b1,c1,d1,a2,b2,c2,d2,+1) 
 
% the result 
 
a = 
     0     1     0     0 
    -1    -2     0     0 
     0     0     0     1 
     1     3    -1    -3 
b = 
     0 
     1 
     0 
     1 
c = 
     0     0     1     4 
d = 
     0 
a = 
     0     1     0     0 
    -1    -2     0     0 
     0     0     0     1 
     0     0    -1    -3 
b = 
     0 
     1 
     0 
     1 
c = 
     1     3     1     4 
d = 
     1 
a = 
     0     1     0     0 
    -1    -2    -1    -4 
     0     0     0     1 
     1     3    -2    -7 
b = 
     0 
     1 
     0 
     1 
c = 
     1     3    -1    -4 
d = 
     1 
a = 
     0     1     0     0 
    -1    -2     1     4 
     0     0     0     1 
     1     3     0     1 
b = 
     0 
     1 
     0 
     1 
c = 
     1     3     1     4 
d = 
     1