%exp4_10.m
%求典型二階系統(tǒng)自然振蕩頻率固定，阻尼比變化時的波特圖
clear
close all
clc
wn=6;
kosi=[0.1:0.1:1.0];
%在對數(shù)空間上生成從10^(-1)到10^1共100個數(shù)據(jù)的橫坐標
w=logspace(-1,1,100);
num=wn^2;
for kos=kosi
   den=[1 2*kos*wn wn^2];
   [mag,pha,w1]=bode(num,den,w);
   % 注意mag的單位不是分貝，若需要分貝表示
   % 可以通過20*log10(mag)進行轉(zhuǎn)換
   subplot(221);
   hold on;
   semilogx(w1,mag)
   % 注意在所繪制的圖形窗口會發(fā)現(xiàn)x軸并沒有取對數(shù)分度
   subplot(222)
   hold on;
   semilogx(w,mag)
   subplot(223);
   hold on;
   semilogx(w1,pha)
   subplot(224)
   hold on;
   semilogx(w,pha)
end
subplot(221)
grid on
title('bode plot')
xlabel('frequency(w)')
ylabel('amplitude')
text(6.2,5,'kosi=0.1')
text(2,0.5,'kosi=1.0')

subplot(223)
grid on
xlabel('frequency(w)')
ylabel('phase deg')
text(5,-20,'kosi=0.1')
text(2,-85,'kosi=1.0')
hold off