%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