exp4_10.m

matlab与控制系统仿真的PPt多媒体教程和源码

%exp4_10.m
%求典型二阶系统的波特图（阻尼比变化）
clear
close all
clc
wn=6;
kosi=0.1:0.2:1.0;
w=logspace(-2,2,100);% 在对数空间上，行矢量w的值从10^(-2)到10^(2),等距离取100个数据
num=wn^2;
for ko=kosi
den=[1 2*ko*wn wn^2];
[mag,pha,w1]=bode(num,den,w);
   % 注意mag的单位不是分贝，若需要分贝表示
   % 可以通过20*log10(mag)进行转换
subplot(221);
semilogx(w1,mag)
hold on
grid on
title('bode plot')
xlabel('frequency(rad/s)')
ylabel('amplitude')
subplot(222)
semilogx(w,mag) %绘制以x轴为对数分度的半对数坐标图
hold on
subplot(223);
semilogx(w1,pha)
hold on
grid on
xlabel('frequency(rad/s)')
ylabel('phase deg')
subplot(224)
semilogx(w,pha)
hold on
end