典型的二階環節傳函：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
顯示圖形：
