samp5_5.m

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

%Samp5_5
clf
n=0:0.01:2;   %设置频率点
for ii=1:2
   switch ii
   case 1,pos=1;N=5;
   case 2,pos=3;N=10;
   end
   [z,p,k]=besselap(N);   %设计Bessel模拟滤波器
   [b,a]=zp2tf(z,p,k);   %将零点极点增益形式转化为传递函数形式
   [H,w]=freqs(b,a,n);   %求得传递函数的复数频率响应
   magH2=(abs(H)).^2;
   phaH=unwrap(angle(H));    %求得函数的相位角并进行解缠绕运算
   phaH=phaH*180/pi;       %将相位角由弧度转化为度
   posplot=['2,2,' num2str(pos)];   %设置绘图位置字符串
   subplot(posplot);
   plot(w,magH2); grid on       %绘出平方幅频响应
   title(['N=' num2str(N)]);
   xlabel('w/wc');   ylabel('Bessel |H(jw)|^2');
   grid on
   subplot(['2,2,' num2str(pos+1)]);
   plot(w,phaH);
   xlabel('w/wc');ylabel('Bessel相位/^o'); 
 title(['N=' num2str(N)]);
   grid on
end