%《MATLAB在電子信息課程中的應(yīng)用》第八章例8.9程序q809
% 連續(xù)和離散系統(tǒng)的多種輸出響應(yīng)曲線
% 電子工業(yè)出版社出版  陳懷琛 吳大正 高西全合著 2001年10月

clear,clf
[a,b,c,d]=rmodel(4);
s1=ss(a,b,c,d);Ts=0.2;
sd1=c2d(s1,Ts,'t')
t=0:Ts:25;
u=sin(0.5*t);
for i=1:2
   if i==1 s=s1;
   else s=sd1; end
   figure(i)
   subplot(2,2,1)
   impulse(s,5);grid
   subplot(2,2,2)
   lsim(s,u,t);
   subplot(2,2,3)
   x0=[1,-1,0,2];
   initial(s,x0,5),grid
   subplot(2,2,4)
   pzmap(s),grid
end