%《MATLAB及在電子信息課程中的應用》第七章例7.31程序q731
% 直接設計帶通數字橢圓濾波器
% 電子工業出版社出版  陳懷琛 吳大正 高西全合著 2001年10月初版，2003年7月第二版
%
clear;close all
Wp=[0.25,0.45];
Ws=[0.15,0.55];
Rp=0.1;
Rs=60;
[N,wc]=ellipord(Wp,Ws,Rp,Rs)            % 計算橢圓濾波器階數
[b,a]=ellip(N,Rp,Rs,wc)                 % 直接計算橢圓濾波器的系數
[hw,w]=freqz(b,a);                      % 檢驗該濾波器的頻率特性
subplot(2,1,1);
set(gcf,'color','w')                    % 設置圖形背景色為白色
plot(w/pi,20*log10(abs(hw)));grid on    % 畫出該濾波器的幅頻特性
axis([0,1,-80,5]);xlabel('ω/π');
ylabel('幅度（dB）')
subplot(2,1,2);plot(w/pi,angle(hw));    % 畫出該濾波器的相頻特性
grid on;axis([0,1,-pi,pi])
xlabel('ω/π');ylabel('相位（rad）')