samp5_3.m

MATLAB7.x数字信号处理 光盘内容

%Samp5_3
n=0:0.01:2;  %设置频率点
for ii=1:4
   switch ii
   case 1,N=2;
   case 2,N=4;
   case 3,N=6;
   case 4,N=8;
   end
   Rs=16;     %阻带衰减设置为16dB
   [z,p,k]=cheb2ap(N,Rs);   %设计Chebyshev II型模拟原型滤波器
   [b,a]=zp2tf(z,p,k);    %转化为传递函数
   [H,w]=freqs(b,a,n);    %求出滤波器的频率响应
   magH2=(abs(H)).^2;     %频率响应的幅度平方
   posplot=['2,2,' num2str(ii)];
   subplot(posplot);
   plot(w,magH2);    %绘出幅频平方响应函数
   title(['N=' num2str(N)]);
   xlabel('w/wc');
   ylabel('Chebyshev II |H(jw)|^2');
   grid on;
end