新MATLAB辅助现代工程数字信号处理 第二版 教学课件 李益华 第1 5章 第4章.ppt

第4章 工程数字滤波器设计 4.1 数字滤波器的技术指标与设计方法4.1.1 数字滤波器的技术指标　　物理上可实现的线性时不变系统是因果的，具有频率响应 图 4.1 物理上可实现滤波器的幅度特性 　　从图4.1可知，数字低通滤波器的技术指标主要有: 　　(1) 通带截止频率ωp;　　(2) 阻带截止频率ωs;　　(3) 通带波纹δ1;　　(4) 阻带波纹δ2;　　(5) 通带内允许的最大衰减αp: 　　　　　αp=－20 lg(1－δ1)=－20 lg|H(ejωp)| 图 4.1 物理上可实现滤波器的幅度特性 4.1.2 数字滤波器的设计方法　　数字滤波器的一般设计步骤如下: 　　(1) 按照实际需要确定滤波器的性能要求;　　(2) 用一个因果稳定系统的H(z)或h(n)去逼近该性能要求，即求h(n)的表达式， 再确定系数ai、bi或零、极点ci、di，以使滤波器满足给定的性能要求; 　　(3) 用一个有限精度的运算去实现这个系统函数，包括选择运算结构(如级联型、并联型、卷积型、频率采样型、快速卷积(FFT)型等);　　(4) 选择合适的字长和有效的数字处理方法，用数字电路或在计算机上编写软件实现所设计的H(z)。 4.2 无限冲激响应(IIR)数字滤波器的设计4.2.1 IIR滤波器的经典设计　　从模拟滤波器设计IIR数字滤波器，就是按照一定的转换关系将s平面上的Ha(s)转换成z平面上的H(z)。这种映射变换应遵循如下两个基本原则: 　　(1) H(z)的频响要能模仿Ha(s)的频响，即s平面的虚轴应映射到z平面的单位圆ejω上。　　(2) Ha(s)的因果稳定性映射成H(z)后保持不变，即s平面的左半平面应映射到z平面的单位圆以内。 　　IIR数字滤波器经典设计法的步骤如下: 　　(1) 将给定的数字滤波器的技术指标转换为模拟滤波器的技术指标; 　　(2) 估计模拟滤波器的最小阶数和边界条件，可利用MATLAB的工具函数buttord、cheb1ord等; 　　(3) 根据转换后的技术指标设计模拟低通滤波器原型，可利用MATLAB的工具函数buttap、cheb1ap、ellipap等; 　　(4) 由频率变换将模拟低通滤波器原型转化为所需的模拟滤波器(低通、高通、带通等)，可利用MATLAB的工具函数lp2lp、lp2hp、lp2bp等; 　　(5) 按照一定规则将模拟滤波器转换为数字滤波器。 　　1. 脉冲响应不变法　　脉冲响应不变法是从滤波器的脉冲响应出发，使数字滤波器的单位脉冲响应序列h(n)正好等于模拟滤波器的冲激响应ha(t)的采样值，即　　　　　　　　　h(n)=ha(t)|t=nT　　　　　　　　　(4.2)式中，T为采样周期。　　因此，数字滤波器的系统函数H(z)可由下式求得:　　　　　　　H(z)=Z [h(n)]=Z [ha(nT)]　　　　　　(4.3)　　采用脉冲响应不变法将模拟滤波器变换为数字滤波器时，它所完成的s平面到z平面的变换，正是拉普拉斯变换到Z变换的标准变换关系z=esT。这种映射关系反映的是Ha(s)的周期延拓与H(z)的关系，而不是Ha(s)本身与H(z)的关系。 　　脉冲响应不变法的数字角频率ω和模拟角频率Ω满足线性变换关系： ω=ΩT。　　在MATLAB中，可用函数impinvar来实现脉冲响应不变法， 调用格式为 　　　[bz， az]= impinvar(b， a， fs) 　　　[bz， az]= impinvar(b， a)其中，b和a分别为模拟滤波器的分子和分母多项式系数向量; fs为采样频率，缺省值为1 Hz; bz和az分别为数字滤波器的分子和分母多项式系数向量。 　　【例 4.1】 用脉冲响应不变法设计一个巴特沃斯低通数字滤波器，设计指标为ωp=0.2 π,ωs=0.3 π,rp=3 dB, rs=15 dB，采样频率fs=10 kHz。　　MATLAB程序如下: 　　 　%MATLAB PROGRAM 4-1　　　wp=0.2*pi*10000; %数字指标转化为模拟指标　　　ws=0.3*pi*10000; 　　　rp=3; 　　　rs=15; 　　　fs=10000; 　　　Nn=256; 　　　[n, wn]=buttord(wp, ws, rp, rs,′s′); %计算阶数n和截止频率wn　　　[z，p，k]=buttap(n); 　　　%设计模拟低通滤波器原型　　　[bap，aap]=zp2tf(z，p，k); %传递函数形式　　　[b，a]=lp