离散时间滤波器的设计

《离散时间滤波器的设计》由会员分享，可在线阅读，更多相关《离散时间滤波器的设计（62页珍藏版）》请在金锄头文库上搜索。

1、8.1 技术指标 8.2 由连续时间滤波器设计IIR滤波器8.2.1 脉冲响应不变法8.2.2 双线性变换法 8.3 窗函数法设计FIR滤波器8.3.1 设计思想8.3.1 布莱克曼窗族8.3.2 凯泽窗族 8.4 FIR滤波器的等波纹最佳逼近设计离散时间滤波器的设计滤波器的设计与实现包括以下几个步骤：(1)根据需要制定滤波器的技术指标即频率响应，包括幅度和 相位两方面的要求；(2)用因果稳定的离散时间系统逼近这些技术指标。对于FIR滤 波器就是找到一个频率响应满足指标要求的因果有限长的 单位脉冲响应；而对IIR滤波器则是找到一个频率响应满足 指标要求的因果稳定的有理系统函数；(3)为实现该滤

2、波器选择一有限精度的软件算法或硬件结构；(4)验证在选定的实现算法和精度情况下频响特性是否符合指 标要求，如果产生了偏差则对步骤（3）的结果进行修正 ；(5)实现该滤波器。8.1 技术指标等效的模拟系统通带容差阻带容差通带截止频率阻带截止频率3dB/半功率 截止频率幅度响应容限图若给出的是等效模拟系统的幅度响应,即横坐标是 , 给 出的截止频率是 ,则需转换得到数字的截止频率 :对数增益: 3dB截止频率对应的对数增益 :在工程实际中，相对指标最受欢迎。8.2 由连续时间滤波器设计IIR滤波器首先设计一个连续时间滤波器（称为原型连续时间滤波 器），得到系统函数Hc(s) ； 再将系统函数Hc(

3、s)转换成离散时间滤波器的系统函数 H(z) 。Hc(s)到H(z)的转换过程实际上是s平面到z平面的映射 过程，该映射需要满足两个要求：一是要求所得到的离散时间滤波器的频率响应保持连续 时间滤波器频率响应的基本特征，即s平面的虚轴映射 到z平面的单位圆，且离散时间滤波器与连续时间滤波 器的频率响应的类型保持不变；二是要求将一个因果稳定的连续时间滤波器转换成一个 因果稳定的离散时间滤波器，即将在s左半平面的极点 全部转换成在z平面单位圆内的极点。虚轴映射 到单位圆左半平面的极点映 射到单位圆内S平面映射 到Z平面8.2.0 连续时间滤波器设计简介区别： （1）波动不同 （2）相同阶数 ， 性能

4、增加 （3）设计复杂度增3种常用的连续时间滤波器的幅度响应：巴特沃斯滤模拟低通波器的设计公式：幅度平方函数：计算阶数和3 分贝截止频率 ：OR计算极点:写出系统函数：向上取整全极点型N,Wc=buttord(2000*pi,4000*pi,1,15, s) Bs,As=butter(N,Wc, s) H,W=freqs(Bs,As); plot(W/2/pi,20*(log10(abs(H) axis(1000,2000,-16,0) grid on设计模拟低通滤波器:举例输出: N = 4 Wc = 8.1932e+003 Bs = 1.0e+015 * 0 0 0 0 4.5063 As

5、= 1.0e+015 * 0.0000 0.0000 0.0000 0.0014 4.5063幅度响应设计模拟低通chebyI型滤波器:N,Wc=cheb1ord(2000*pi,4000*pi,1,15, s) Bs,As=cheby1(N,1,Wc, s) H,W=freqs(Bs,As); plot(W/2/pi,20*(log10(abs(H) axis(0,4000,-30,0) grid on举例设计模拟低通chebyII 型滤波器:N,Wc=cheb2ord(2000*pi,4000*pi,1,15, s) Bs,As=cheby2(N,15,Wc, s) H,W=freqs(B

6、s,As); plot(W/2/pi,20*(log10(abs(H) axis(0,4000,-30,0) grid on举例设计模拟高通滤波器:N,Wc=buttord(4000*pi,2000*pi,1,15, s) Bs,As=butter(N,Wc, high, s) H,W=freqs(Bs,As); plot(W/2/pi,20*(log10(abs(H) axis(0,4000,-16,0) grid on举例带通带阻滤波器的通带和阻带 截止频率分别是2维,带阻用参 数stop8.2.1 脉冲响应不变法1.设计思想 ：使离散时间滤波器的单位脉冲响应为原型连续时间滤波 器的单位脉

7、冲响应的取样 :就可以根据原型连续时间滤波器的系统函数确定离 散时间滤波器的系统函数:先将连续时间滤波器的系统函数分解成部分分式：2.系统函数转换方法 再写出离散时间滤波器的系统函数:证明见课堂笔记3. 因果稳定性极点的映射关系:4频率响应的映射特点如果连续时间滤波器是带限的，即 则混迭小的情况下，近似认为频响是线性关系。离散时间滤波器与原型连续时间滤波器的频率响应的自变 量之间的映射关系是： s 平面z 平面虚轴映射到单位圆是一对多，频率轴每绕一圈频响累 加一次，所以频响有混迭。(1) 离散时间滤波器的指标 原型连续时间滤波器的指标 (2) 设计(3)注意：如果给出的等效模拟指标，则必须 首

8、先按照比例T线性转换成数字指标。5.设计步骤 6. 关于 Td若从离散时间指标开始，则Td的取值大小不影响混迭程度 和离散时间系统。所以可以任意取值。 为了简单起见，常令 优点:频率的线性映射;缺点：频率响应有混迭，应用只限于带限滤波器的设计， 应用于高通或带阻滤波器设计时会产生严重的频响混叠。7.特点wp=0.2*pi; ws=0.4*pi ap=1;as=12 Td=1;Wp=wp/Td;Ws=ws/Td N,Wc=buttord(Wp,Ws, ap , as, s) Bs,As=butter(N,Wc, s) Bz,Az=impinvar(Bs,As,1/Td) H,W=freqs(Bs

9、,As); plot(W/pi,20*(log10(abs(H), r*) hold on H,w=freqz(Bz,Az); plot(w/pi,20*(log10(abs(H) axis(0.2,0.4,-20,0) grid举例输出： Bs = 0 0 0 0.5150 As =1.0000 1.6031 1.2849 0.5150 Bz = 0 0.1453 0.0855 0 Az =1.0000 -1.4782 0.9106 -0.2013手算：N,Wc=buttord(0.6*pi,0.5*pi, 1 ,12, s) Bs,As=butter(N,Wc, high, s) Bz,A

10、z=impinvar(Bs,As,1) H,W=freqs(Bs,As); plot(W/pi,20*(log10(abs(H), r) figure H,w=freqz(Bz,Az); plot(w/pi,20*(log10(abs(H)举例思考本例的问题出在哪里。8.2.2 双线性变换法1.设计思想s平面的整个虚轴非线性地一对一地映 射到z平面的单位圆上。2.系统函数转换方法如果 0, 则|z|1 ，即s左半平面的所有点都映射到z平面的单 位圆内，包括极点，所以如果原型连续时间滤波器因果稳定则 离散时间滤波器也一定因果稳定 。3.因果稳定性|z|=1. 即 虚轴映射到单位圆上 (一对一映射

11、). 将 带入上式得到 如果 则4.频率响应的映射特点频率响应间的关系需要采用预畸变确定原型连续时间滤波器的技术指标:1) 离散时间滤波器指标 原型连续时间滤波器指标2)3)5.设计步骤 注意：如果给出的等效模拟指标，则必须 首先用比例T线性转换成数字指标。6. 关于 Td 可任取，一般取1优点：避免了脉冲响应不变法中的频率响应混叠问题； 缺点:引入了频率失真，频率失真问题可以通过预畸变加以解 决。只适用于滤波器具有近似理想的分段恒定幅度响应的情 况，并且它无法实现线性的幅度或相位映射 7. 特点例如模拟微分器（幅度响应线性）不能通过双线性变换法 得到数字微分器。wp=0.2613*pi;ws

12、=0.4018*pi;ap=0.75; as=20; Td=1;Ws=2/Td*tan(ws/2); Wp=2/Td*tan(wp/2) N,Wc=buttord(Wp,Ws,ap,as,s) Bs,As=butter(N,Wc, s) Bz,Az=bilinear(Bs,As,1/Td) H,W=freqs(Bs,As); plot(W/pi,20*(log10(abs(H),Rx) hold on H,w=freqz(Bz,Az); plot(w/pi,20*(log10(abs(H) ylabel(虚线：模拟滤波器幅度dB 实线：数字滤波器幅度dB) xlabel(虚线：模拟角频率*弧度

13、/秒 实线:数字角频率*弧度) axis(0.25,0.5,-20,-0.45) grid举例或着 wp=0.2613*pi;ws=0.4018*pi; ap=0.75;as=20 N,wc=buttord(wp/pi,ws/pi,ap,as) Bz,Az=butter(N,wc) H,w=freqz(Bz,Az); plot(w/pi,20*(log10(abs(H) N,Wc=buttord(0.6,0.5,1,12); b,a=butter(N,Wc, high); H,w=freqz(b,a); plot(w/pi,20*(log10(abs(H) grid举例设计离散时间带通滤波器举

14、例N,wc=buttord(0.45 0.55,0.4 0.6,3,10) %双线性变换法 B,A=butter(N,wc) H,w=freqz(B,A); plot(w/pi,20*(log10(abs(H) ylabel(20log|H(ej)| dB) xlabel(数字角频率*弧度) axis(0.4,0.6,-10,0);grid onOutput: N = 2 wc = 0.4410 0.5590 B = 0.0271 0 -0.0541 0 0.0271 A = 1.0000 0 1.4838 0 0.59208.2节 总结 1. 设计步骤2.冲击响应不变法：频率轴线性多对一映射,频响有混迭,不适用于高通等.双线性变换法：频率轴有畸变一对一映射,频响无混迭,不适用于微分器.8.3 窗函数法设计FIR滤波器 8.3.1 设计思想（1）先对具有广义线性相位的理想频率响应作傅里叶反变换 ，得到无限长相对M/2偶对称的单位脉冲响应 （2）然后对hdn加长度为M的相对M/2偶对称的窗进行截短， 得到有限长的相对M/2偶对称的序列hn就是一个逼近理想频率响应的因果FIR滤波器的单位脉冲 响应，其频率响应是：加窗的方式（包括窗形状和窗长）决定了滤波器的频率响应对 理想频率响应的逼近程度。 窗函数的主瓣宽（与窗 长成反比）越大，则滤 波器的过渡带宽越宽； 窗函数的旁瓣面积与主