%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