典型的二阶环节传函：G(s）=wn*wn/(s^2+2*sigma*wn*s+wn*wn);自然震荡角频率wn=0.7,绘制不同sigma取值下的Bode图和Nyquist图如下：
Bode图：
程序代码如下：
>> w=[0,logspace(-2,2,200)]; %w为~之间对数等间距分布的200个数
wn=0.7; %自然震荡角频率
sigma=[0.1,0.4,0.707,1,3,5,7]; %取不同的阻尼比
for j=1:7
num=[wn*wn];den=[1,2*sigma(j)*wn,wn*wn];
sys=tf(num,den); %建立系统传函
bode(sys,w);hold on; %绘制波特图
end
>> gtext('sigma=0.1'); gtext('sigma=0.4'); gtext('sigma=0.707'); gtext('sigma=1'); gtext('sigma=3'); gtext('sigma=5'); gtext('sigma=7'); %放置sigma取不同值时的文字注释
显示图形如下：
Nyquist图：
程序代码如下：
>> w=[0,logspace(-2,2,200)];
>> wn=0.7;
>> sigma=[0.1,0.4,0.707,1,3,5,7]; %取不同的阻尼比
>> for j=1:7
num=[wn*wn];den=[1,2*sigma(j)*wn,wn*wn];
nyquist(num,den,w);hold on;
end
>> grid
显示图形：
