samp5_2.m

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

%Samp5_2
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
   Rp=1;      %设置通带波纹为1dB
   [z,p,k]=cheb1ap(N,Rp);   %设计Chebyshev I型滤波器
   [b,a]=zp2tf(z,p,k);    %转换为传递函数形式
   [H,w]=freqs(b,a,n);    %求得传递函数的频率特性
   magH2=(abs(H)).^2;     %求得传递函数的幅频响应
posplot=['2,2,' num2str(ii)];%将数字ii转换为字符串，与’2,2,’合并并赋给posplot。
   subplot(posplot);
   plot(w,magH2);
   title(['N=' num2str(N)]);%将数字N转换为字符串与’N=’合并作为标题。
   xlabel('w/wc');
   ylabel('Chebyshev I |H(jw)|^2');
grid on
end