q732.m

MATLAB及在电子信息课程中的应用例题索引与程序

%《MATLAB在电子信息课程中的应用》第七章例7.32程序q732
% 用双线性变换法设计高通和带通巴特沃斯数字滤波器
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著 2001年10月

%用双线性变换法设计数字高通和带通滤波器
clear;close all
T=1;wch=pi/2;			% T为采样间隔，wch为数字高通3dB截止频率
wlc=0.35*pi;wuc=0.65*pi;		% wlc,wuc:数字带通3dB截止频率
B=1;A=[1,2.6131,3.4142,2.6131,1];
[h,w]=freqs(B,A,512);			% 求原归一化模拟滤波器的频率响应
subplot(3,2,1);plot(w,20*log10(abs(h)));	% 画模拟滤波器幅频特性
grid;axis([0,10,-90,0])
xlabel('ω/π');ylabel('模拟低通幅度(dB)')
%(1)设计高通
omegach=2*tan(wch/2)/T; 		% 预畸变求模拟高通3dB截止频率
[Bhs,Ahs]=lp2hp(B,A,omegach);		% 模拟域低通转换为高通系数
[Bhz,Ahz]=bilinear(Bhs,Ahs,1/T);	% 模拟转换为数字高通系数向量
[h,w]=freqz(Bhz,Ahz,512);			% 求并画数字滤波器幅频特性
subplot(3,2,3);plot(w/pi,20*log10(abs(h)));
grid;axis([0,1,-150,0])
xlabel('ω/π');ylabel('数字高通幅度(dB)')
%(2)设计带通
omegalc=2*tan(wlc/2)/T; 	% 预畸变求滤波器通带低端截止频率
omegauc=2*tan(wuc/2)/T;		% 预畸变求滤波器通带高端截止频率
wo=sqrt(omegalc*omegauc);Bw=omegauc-omegalc;
[Bbs,Abs]=lp2bp(B,A,wo,Bw);	%		% 模拟域低通转换为带通系数
[Bbz,Abz]=bilinear(Bbs,Abs,1/T);	% 模拟转换为数字带通系数向量
[h,w]=freqz(Bbz,Abz,512);			% 求并画数字滤波器幅频特性
subplot(3,2,4);plot(w/pi,20*log10(abs(h)));
grid;axis([0,1,-150,0])
xlabel('ω/π');ylabel('数字带通幅度(dB)')