samp5_1.m

来自「数字信号处理的Matlab实现的随书光盘」· M 代码 · 共 25 行

%Samp5_1
n=0:0.01:2;    %频率点
for ii=1:4     %取4种滤波器
   switch ii
   case 1,N=2;
   case 2,N=5;
   case 3,N=10;
   case 4,N=20;
   end
   [z,p,k]=buttap(N);    %设计Butterworth滤波器
   [b,a]=zp2tf(z,p,k);   %将零点极点增益形式转换为传递函数形式
   [H,w]=freqs(b,a,n);   %按n指定的频率点给出频率响应
   magH2=(abs(H)).^2;    %给出传递函数幅度平方
   hold on;
   plot(w,magH2);     %绘制传递函数幅度平方
end
xlabel('w/wc');
ylabel('|H(jw)|^2');
title('Butterworth模拟原型滤波器');
text(1.5,0.18,'n=2')   %作必要的标记
text(1.3,0.08,'n=5')
text(1.16,0.08,'n=10')
text(0.93,0.98,'n=20')
grid on;