samp5_4.m

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

%Samp5_4
n=0:0.01:2;   %设置频率点
for ii=1:4
   switch ii
   case 1,N=2;
   case 2,N=3;
   case 3,N=4;
   case 4,N=5;
   end
   Rp=1; Rs=15;   %设置通带波纹为1dB，阻带衰减为15dB.
   [z,p,k]=ellipap(N,Rp,Rs);   %设计椭圆滤波器
   [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); ylim([0 1])       %绘出幅度平方函数
   title(['N=' num2str(N)]);
   xlabel('w/wc');
   ylabel('椭圆|H(jw)|^2');
   grid on;
end