p178_1.m

来自「matlab仿真的文件包然后写出其具体功能(至少要20个字)。尽量不要让站长把时」· 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)