幅频特性曲线的matlab模拟
幅频特性曲线的 MATLAB 模拟机自 1206 mldxs 41263010代码如下：s=0:0.01:5;y0=(200.^(-1))*(sqrt(((1-s.^2)).^2+(2*1*s).^2)).^(-1)y1=(200.^(-1))*sqrt(((1-s.^2)).^2+(2*0.5*s).^2).^(-1)y2=(200.^(-1))*sqrt(((1-s.^2)).^2+(2*0.375*s).^2).^(-1)y3=(200.^(-1))*sqrt(((1-s.^2)).^2+(2*0.25*s).^2).^(-1)y4=(200.^(-1))*sqrt(((1-s.^2)).^2+(2*0.125*s).^2).^(-1)y5=(200.^(-1))*(sqrt(((1-s.^2)).^2+(2*0.1*s).^2)).^(-1)plot(s,y0, -r ,s,y1, -r ,s,y2, -r ,s,y3, -r ,s,y4, -r ,s,y5, -r ),grid ontitle( 幅频特性 )ylabel( |H(jw)| ),xlabel( w/wn )gtext( 1 )gtext( 0.5 )gtext( 0.375 )gtext( 0.25 )gtext( 0.125 )gtext( 0.1 )gtext( ζ=c/cc )相频特性曲线的 MATLAB 模拟代码如下：s=0:0.001:5y0=atan((2*1*s).*(1-s.^2).^(-1))y1=atan((2*0.5*s).*(1-s.^2).^(-1))y2=atan((2*0.375*s).*(1-s.^2).^(-1))y3=atan((2*0.25*s).*(1-s.^2).^(-1))y4=atan((2*0.125*s).*(1-s.^2).^(-1))y5=atan((2*0.10*s).*(1-s.^2).^(-1))plot(s,y0, -r ,s,y1, -r ,s,y2, -r ,s,y3, -r ,s,y4, -r ,s,y5, -r ),grid ontitle( 相频特性 )ylabel( φ ),xlabel( w/wn )gtext( 1 )gtext( 0.5 )gtext( 0.375 )gtext( 0.25 )gtext( 0.125 )gtext( 0.1 )