%
% Example 8.10 and Figures 8.3 and 8.4
%
A=[0 1; -2 -3]
B=[0;1]
C=[1 0]
T=0.1
[Ad,Bd]=c2d(A,B,T)
Dd=0
Cd=C
[ys,xs]=dstep(Ad,Bd,Cd,Dd)
[yi,xi]=dimpulse(Ad,Bd,Cd,Dd)
figure (1)
xiL=length(xi)
k=0:1:xiL-1
plot(k,xi(1:xiL,1:1),'*',k,xi(1:xiL,2:2),'o')
xlabel('discrete time k')
ylabel('system impulse responses')
text(10,0.03,'xi1=yi')
text(30,-0.01,'xi2')
print -deps fig8_3.eps
figure (2)
xsL=length(xs)
k=0:1:xsL-1
plot(k,xs(1:xsL,1:1),'*',k,xs(1:xsL,2:2),'o')
xlabel('discrete time k')
ylabel('system step responses')
text(30,0.43,'xs1=ys')
text(30,0.07,'xs2')
print -deps fig8_4.eps