% Example 6.35

num = [1 2 16];
den = [1 4 8 0];
H = tf(num,den);
t = 0:10/300:10;

% response to exponential
x = exp(-2*t);
y=lsim(H,x,t);
ya = exp(-2*t).*cos(2*t)-.5*exp(-2*t).*sin(2*t) + 1 - 2*exp(-2*t);
plot(t,y,t,ya,'--')  % ya is actual response from example 6.30
xlabel('Time (sec)')
ylabel('y(t)')
title('Response to exp(-t)')
text(2,.6,'----- Result computed from lsim')
text(2,.5,'- - - Analytical result from Ex. 6.30')
pause

% step response
y =step(H,t);    
plot(t,y)
xlabel('Time (sec)')
ylabel('y(t)')
title('Step Response')
pause
% impulse response
y=impulse(H,t);
plot(t,y)
xlabel('Time (sec)')
ylabel('y(t)')
title('Impulse Response')