程序: f=-5:0.001:5;y=1./sqrt(1+(2*pi*f).^2);plot(f,y)xlabel('f')ylabel('y')title('y=1/√[1+(2π f)^2] ')grid on 结果: