example6_21.m

来自「<MATLAB 7辅助控制系统设计与仿真>书籍内的源代码.」· M 代码 · 共 29 行

k1=16.7/0.125;
z1=[0];
p1=[-1.25 -4 -16];
[num1,den1]=zp2tf(z1,p1,k1);
[num,den]=cloop(num1,den1);
[z,p,k]=tf2zp(num,den);
p
figure(1)
nyquist(num,den);
title('Nyquist 曲线图');
xlabel('实数轴');
ylabel('虚数轴');
figure(2)
[num2,den2]=cloop(num,den);
impulse(num2,den2);
title('冲激响应');
xlabel('时间');
ylabel('振幅');