p178_1.m

来自「反馈控制系统分析与设计,虽然有点抽象」· M 代码 · 共 11 行

wn=1; zet=[0:0.1:1,2,3,5]; 
w=logspace(-1,1,100); mm=[]; pp=[];
for i=1:length(zet)
   m=wn^2./sqrt((wn^2-w.^2).^2+4*zet(i)^2*wn^2*w.^2);
   p=atan2(-2*zet(i)*wn*w,wn^2-w.^2);
   pp=[pp,180*p'/pi]; mm=[mm,20*log10(m')];
end
subplot(211); semilogx(w,mm), 
subplot(212); semilogx(w,pp)