A1=[-1.5,-0.8,0,0;0.8,0,0,0;0.3,0.4,-4.0,-1.25;0,0,-1.25,0];
B1=[1,0,1,0]';
C1=[1,2,1,2];
D1=0;
t=0.5;
[A,B,C,D]=c2dm(A1,B1,C1,D1,t,'tustin');
subplot(2,2,1)
dstep(A,B,C,D)
grid on;
title('離散階躍響應')
xlabel('時間')
ylabel('振幅')
subplot(2,2,2)
dimpulse(A,B,C,D)
grid on;
title('離散沖激響應')
xlabel('時間')
ylabel('振幅')
subplot(2,2,3)
x0=[1,1,1,1];
dinitial(A,B,C,D,x0)
grid on;
title('零輸入響應')
xlabel('時間')
ylabel('振幅')
axis([0 6 -0.5 2.5])
subplot(2,2,4)
[z,p,k]=ss2zp(A,B,C,D,1)
zplane(z,p)
grid on;
title('離散零極點圖')
xlabel('實部')
ylabel('虛部')