example6_9.m

<MATLAB 7辅助控制系统设计与仿真>书籍内的源代码.

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('虚部')