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

% Example 3.28 
% 
k1=16.7/.0125; 
z1=[0]; 
p1=[-1.25 -4 -16]; 
[num1,den1]=zp2tf(z1,p1,k1); 
[num,den]=cloop(num1,den1); 
[z,p,k]=tf2zp(num,den); 
p 
figure(1) 
nyquist(num,den) 
title('Nyquist Plot') 
figure(2) 
[num2,den2]=cloop(num,den); 
impulse(num2,den2) 
title('Impulse Response')