本文举了4种情况的Nyquist图的MATLAB程序,供参考
1.惯性环节
clear
clc
k=24;
nunG1=k;
% denG1=conv(conv([0,1],[0,1]),conv([1 2],[1 3]));
denG1=[1 1];
[re,im]=nyquist(nunG1,denG1);
plot(re,im);
grid
图1
2.比例+积分+惯性
clear
clc
k=24;
nunG1=k;
% denG1=conv(conv([0,1],[0,1]),conv([1 2],[1 3]));
denG1=conv([1 0],[1 3]);
[re,im]=nyquist(nunG1,denG1);
plot(re,im);
grid
3.比例+2个积分+2个惯性
clear
clc
k=24;
nunG1=k;
denG1=conv(conv([1,0],[1,0]),conv([1 2],[1 3]));
图3
4.比例+积分+导前+惯性
clear
clc
k=24;
nunG1=k*[9 1];
denG1=conv([1,0],[2 1]);
%t1要比t2大
[re,im]=nyquist(nunG1,denG1);
plot(re,im);
grid
[re,im]=nyquist(nunG1,denG1);
plot(re,im);
grid
图4