第六章无限脉冲响应数字滤波器的设计.ppt

程序清单[bz,az]=impinvar(b,a,fs)T=1;Fs=1/T;wp=0.2*pi/T;ws=0.35*pi/T;rp=1;rs=10;[N1,wc1]=buttord(wp,ws,rp,rs,s);[B1,A1]=butter(N1,wc1,s);fk=0:Fs/2/512:Fs/2;wk=2*pi*fk;Hk=freqs(B1,A1,wk);subplot(2,2,1);plot(fk,20*log10(abs(Hk)));gridon；xlabel(‘f(Hz)’);ylabel(‘幅度(dB))axis([0,0.5,-50,5]);title((a));第93页，共133页，星期日，2025年，2月5日程序清单[bz,az]=impinvar(b,a,fs)[Bz1,Az1]=impinvar(B1,A1,Fs);Hz=freqz(Bz1,Az1,wk);subplot(2,2,3);plot(wk/pi,20*log10(abs(Hz)));gridon;xlabel(‘\omega/\pi’);ylabel(‘幅度(dB))axis([0,1,-50,5]);title((b));第94页，共133页，星期日，2025年，2月5日③设计归一化低通G(p)。采用巴特沃斯滤波器④求模拟高通第61页，共133页，星期日，2025年，2月5日巴特沃斯高通滤波器的MATLAB实现（1）实现函数：LP2HP调用格式：[bh,ah]=lp2hp(b,a,wph)wp=1;ws=4;Rp=0.1;As=40;[N,wc]=buttord(wp,ws,Rp,As,s);[B,A]=butter(N,wc,s);wph=2*pi*4000;[BH,AH]=lp2hp(B,A,wph);例6.3.2程序：第62页，共133页，星期日，2025年，2月5日例6.3.2程序结构图第63页，共133页，星期日，2025年，2月5日（2）直接用buttord和butter实现巴特沃斯高通滤波器的MATLAB实现例6.3.2程序：wp=2*pi*4000;ws=2*pi*1000;Rp=0.1;As=40;[N,wc]=buttord(wp,ws,Rp,As,s);[BH,AH]=butter(N,wc,high,s);实际应用中常用该种方法第64页，共133页，星期日，2025年，2月5日二、模拟低通到模拟带通的变换Ωpu：通带上截止频率Ωpl：通带下截止频率Ωsl：阻带下截止频率Ωsu：阻带上截止频率Bw=Ωpu-Ωpl:带通滤波器的带宽通带中心频率第65页，共133页，星期日，2025年，2月5日模拟低通到模拟带通的频率变换关系第66页，共133页，星期日，2025年，2月5日模拟低通到模拟带通的频率变换关系模拟带通滤波器的系统函数为：频率变换关系函数变换关系第67页，共133页，星期日，2025年，2月5日(1)确定模拟带通滤波器的技术指标Ωpu，Ωpl，Ωsu，Ωsl模拟带通的设计步骤(2)确定归一化低通滤波器的技术指标(3)设计归一化低通G(p)。(4)将G(p)转换成带通H(s)。实现函数：lp2bp第68页，共133页，星期日，2025年，2月5日设计巴特沃斯模拟带通滤波器，要求通带上、下边界频率分别为4kHz和7kHz，阻带上、下边界截止频率分别为2kHz和9kHz，通带内最大衰减αp=1dB，阻带最小衰减αs=20dB。例6.3.3解(见书P175）程序如下：wp=2*pi*[4000,7000];ws=2*pi*[2000,9000];Rp=1;As=20;[N,wc]=buttord(wp,ws,Rp,As,s);[BB,AB]=butter(N,wc,s);第69页，共133页，星期日，2025年，2月5日例6.3.3结果图：第70页，共133页，星期日，2025年，2月5日三、模拟低通到模拟带阻的变换第71页，共133页，星期日，2025年，2月5日模拟低通到模拟带阻的频率变换关系模拟带阻滤波器的系统函数为:频率变换关系函数变换关系第72页，共133页，星期日，2025年，2月5日设计巴特沃斯、椭圆模拟带阻滤波器，要求阻带上、下边界频率分别为4kHz和7kHz，通带上、下边界截止频率分别为2kHz和9kHz，通带内最大衰减αp=1dB，阻带最小衰减αs=20dB。例6.3.4见书P177实现函数：lp2bs第73页