q731.m

MATLAB及在电子信息课程中的应用,电子工业出版社出版 陈怀琛 吴大正 高西全合著

%《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）')