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

% Example 3.14 
% 
a1=[-1.6,-0.9,0,0; 0.9,0,0,0 
 0.4,0.5,-5.0,-2.45; 0,0,2.45,0]; 
b1=[1;0;1;0]; 
c1=[1,1,1,1]; 
d1=[0]; 
t=0.5; 
[a,b,c,d]=c2dm(a1,b1,c1,d1,t,'tustin'); 
figure(1) 
subplot(2,2,1) 
dstep(a,b,c,d) 
title('Discrete Step Response') 
subplot(2,2,2) 
dimpulse(a,b,c,d) 
title('Discrete Impulse Response') 
subplot(2,2,3) 
x0=[1;1;1;-1]; 
dinitial(a,b,c,d,x0) 
axis([0 6 -0.5 2.5]) 
title('Discrete Initial Response') 
subplot(2,2,4) 
[z,p,k]=ss2zp(a,b,c,d,1); 
zplane(z,p) 
title('Discrete Pole-Zero Map')