%《MATLAB在電子信息課程中的應用》第七章例7.31程序q731
% 直接設計帶通數字橢圓濾波器
% 電子工業出版社出版  陳懷琛 吳大正 高西全合著 2001年10月
                       
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(3,2,1);
plot(w/pi,20*log10(abs(hw)));grid
axis([0,1,-80,5]);xlabel('ω/π');
ylabel('幅度（dB）')
subplot(3,2,3);plot(w/pi,angle(hw));
grid;axis([0,1,-pi,pi])
xlabel('ω/π');ylabel('相位（rad）')