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

% Example 3.26 
% 
a=[0 1 0 0; 0 0 1 0; 0 0 0 1; -62.5 -213.8 -204.2 -54]; 
b=[0 0 0 1]'; c=[1562 1875 0 0]; d=0; 
w=logspace(-1,1); 
figure(1) 
bode(a,b,c,d,1,w) 
subplot(2,1,1) 
title('Bode Plot') 
figure(2) 
nyquist(a,b,c,d,1,w) 
title('Nyquist Plot') 
[z,p,k]=ss2zp(a,b,c,d); 
p 
figure(3) 
[a1,b1,c1,d1]=cloop(a,b,c,d); 
impulse(a1,b1,c1,d1) 
title('Impulse Response')