幅频特性曲线的MATLAB模拟
机自1206 吴文昊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 on
title('幅频特性')
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')