数字滤波器设计课件.ppt

第第 4 章章 数字滤波器设计主 要 内 容è滤波器概述滤波器概述 •典型模拟滤波器的设计典型模拟滤波器的设计•无限冲激响应无限冲激响应(IIR)数字滤波器的设计数字滤波器的设计•有限冲激响应有限冲激响应(FIR)数字滤波器的设计数字滤波器的设计 F 滤波器，是对输入信号起到滤波的作用 提取有用信号，抑制不需要的干扰；提取有用信号，抑制不需要的干扰； 把信号分成不同的频率分量；把信号分成不同的频率分量； 信号检测及信号估计信号检测及信号估计经典滤波现代滤波：4.1 滤波器概述一、滤波器基本原理 1 滤波器：对输入信号起滤波作用h(n)x(n)y(n)则则则则LTI(Linear Time Invariant)LTI(Linear Time Invariant)系统的输出为：系统的输出为：系统的输出为：系统的输出为：Ø假定输入信号x(n)中的有用成分和希望去除的成分，各自占有不同的频带当x(n）经过一个线性系统（即滤波器）后即可将欲去除的成分有效地去除X(ejw)|wwc无用wc|H(ejw)||Y(ejw)|wwcv 经典滤波器从功能上分又可分为：低通滤波器（低通滤波器（LPAF/LPDF):Low pass filterLPAF/LPDF):Low pass filter带通滤波器带通滤波器(BPAF/BPDF):Bandpass filter(BPAF/BPDF):Bandpass filter高通滤波器高通滤波器(HPAF/HPDF):High pass filter(HPAF/HPDF):High pass filter带阻滤波器带阻滤波器(BSAF/BSDF):Bandstop filter(BSAF/BSDF):Bandstop filterv 每一种又可分为：数字(Digital)和模拟(Analog)滤波器。

2 分类n n模拟滤波器：模拟滤波器：RCRC滤波、滤波、LCLC滤波滤波n n数字滤波器：是一种计算程序，处理离数字滤波器：是一种计算程序，处理离 散信号模拟滤波器的理想幅频特性模拟滤波器的理想幅频特性LPAFHPAFBPAFBSAF数字滤波器的理想幅频特性数字滤波器的理想幅频特性LPDFHPDFBPDFBSDF…….…….…….…….3 技术要求 通带最大衰减：阻带最小衰减：二 数字滤波器vv数字滤波器是用来处理数字信号的专用处理系统例1：已知一个数字滤波器的传输函数实现上式的运算：移位、乘、加则其相应的差分方程为：n n所谓数字滤波实质是一种运算过程，用来描述离散系统输入与输出关系的差分方程或卷积和，给它提供一个简明的运算规则，使它完成对数据的处理n n数字滤波不同于模拟滤波，因为模拟滤波是对函数进行处理例2：已知一个数字滤波器其传输函数为则根据传输函数幅度特性：实现：实现：&& 数字滤波器是一种数字式的处理设备，它的设计是建立在对离散系统分析的基础上但是某些滤波器的设计过程中，为了保留模拟滤波器的良好特性，往往按给定的指标先设计一个模拟滤波器，然后进行数字化，所以了解模拟滤波器的特性和设计方法是很有必要的。

主 要 内 容•滤波器概述滤波器概述 è典型模拟滤波器的设计典型模拟滤波器的设计•无限冲激响应无限冲激响应(IIR)数字滤波器的设计数字滤波器的设计•有限冲激响应有限冲激响应(FIR)数字滤波器的设计数字滤波器的设计 FAF的衰减（衰耗）特性工作衰耗取决于幅度平方函数或称模方函数理想的低通滤波器：通带过渡带阻带实际的低通滤波器：4.2 典型模拟滤波器设计由模方函数求AF的传递函数结论：模方函数与传递函数关系密切，可以从给定 的模方函数求出传递函数一、巴特沃思低通(Butterworth)滤波器v 基本性质 从给定的指标设计模拟滤波器，其中心是如何寻找一个恰当的近似函数来逼近理想特性 巴特沃思滤波器是以巴特沃思近似函数作为滤波器的传输函数，该函数以最高阶泰勒级数的形式来逼近理想矩形特性vv模方函数：随着阶次N的增加，通带越平坦，越接近理想的特性， ，逼近矩形巴特沃斯滤波器幅度特性性质：性质：Ø对对不同的不同的N N，，总存在总存在 ，幅度随，幅度随 着着 的升高单调下降。

的升高单调下降Ø当当 这说明在这说明在 处最处最“ “平坦平坦” ”——最大平坦近似最大平坦近似Ø频率远离 ，在高频范围内，频率每增加一倍，衰耗增加6NdBv 设计过程uu按按给定指标确定给定指标确定N N ，通带最大衰耗 ，阻带最小衰耗uu由由模方函数求传输函数模方函数求传输函数H(s)H(s) 巴特沃思函数极点分布规律极点分布图极点分布图N=3ResjImsS平面N=4ResjImsS平面N阶巴特沃思滤波器的传输函数阶巴特沃思滤波器的传输函数二、切比雪夫(Chebyshev)滤波器Ø切比雪夫滤波器是一种全极型滤波器，它是由切比雪夫多项式的正交函数 推导出来，采用了在通带内等波动，在通带外衰减单调递增的准则去逼近理想特性Ø它比巴特沃斯滤波器在通带内有更均匀的特性，是在所有全极型滤波器中过渡带最窄的最优滤波器 11切比雪夫滤波器幅度特性主 要 内 容•滤波器概述滤波器概述 •典型模拟滤波器的设计典型模拟滤波器的设计è无限冲激响应无限冲激响应(IIR)数字滤波器的设计数字滤波器的设计•有限冲激响应有限冲激响应(FIR)数字滤波器的设计数字滤波器的设计 4.3 无限冲激响应数字滤波器设计一、无限冲激响应数字滤波器一、无限冲激响应数字滤波器 (Infinite impulse response, IIR) 若bi全为零，则有vv具有(1)式形式的滤波器简称为FIR滤波器 (Finite Impulse Response)‚ b bi i中至少有一项不为零，并且分母中至少有中至少有一项不为零，并且分母中至少有一个根不为分子所抵消一个根不为分子所抵消vv具有(2)式形式的滤波器简称为IIR滤波器 (Infinite Impulse Response)二、设计二、设计IIR滤波器的基本条件滤波器的基本条件Ø可实现性（因果的，稳定的）Ø实现从模拟到数字的转换 在一定条件下用离散的变量Z代替连续变量S，则模拟滤波器的传输函数H(S)就有可能变换成IIR的传输函数H(Z)? 结论：结论：l l在实际设计在实际设计在实际设计在实际设计IIRIIRIIRIIR时，常根据指标先设计一个模拟时，常根据指标先设计一个模拟时，常根据指标先设计一个模拟时，常根据指标先设计一个模拟滤波器，然后将其数字化。

滤波器，然后将其数字化滤波器，然后将其数字化滤波器，然后将其数字化l l从模拟到数字的映射要满足两个条件：从模拟到数字的映射要满足两个条件：从模拟到数字的映射要满足两个条件：从模拟到数字的映射要满足两个条件：• •为了保持模拟滤波器频率特性，必须为了保持模拟滤波器频率特性，必须为了保持模拟滤波器频率特性，必须为了保持模拟滤波器频率特性，必须将将将将S S S S平面的复频平面的复频平面的复频平面的复频率轴映射到率轴映射到率轴映射到率轴映射到Z Z Z Z平面的单位圆上平面的单位圆上平面的单位圆上平面的单位圆上• • 为了保持模拟滤波器的稳定性，将为了保持模拟滤波器的稳定性，将为了保持模拟滤波器的稳定性，将为了保持模拟滤波器的稳定性，将S S S S左半平面映射到左半平面映射到左半平面映射到左半平面映射到Z Z Z Z平面的单位圆内平面的单位圆内平面的单位圆内平面的单位圆内三、双线性变换三、双线性变换1 基本思路基本思路 按给定指标设计AF得到H(s) H(z)压缩频带Z变换S平面RezZ平面2. 基本问题&&映射 ：FZ变换：vv为使模拟滤波器某一频率与数字滤波器的任一为使模拟滤波器某一频率与数字滤波器的任一频率有对应关系，引入系数频率有对应关系，引入系数 c c3 变换常数变换常数c的选择的选择•某一特定频率严格相对应：某一特定频率严格相对应：•低频处有较确切的对应关系：低频处有较确切的对应关系：特定频率处频率响应严格相等，可以较准确地特定频率处频率响应严格相等，可以较准确地特定频率处频率响应严格相等，可以较准确地特定频率处频率响应严格相等，可以较准确地控制截止频率位置。

控制截止频率位置控制截止频率位置控制截止频率位置4 逼近情况逼近情况1 1））s s平面虚轴平面虚轴z z平面单位圆平面单位圆2 2））左半平面左半平面单位圆内单位圆内s s平面平面z z平面平面右半平面右半平面单位圆外单位圆外虚轴虚轴单位圆上单位圆上5 优缺点优缺点• • 优点：优点：避免了频率响应的混迭现象避免了频率响应的混迭现象s s 平面与平面与 z z 平面为单值变换平面为单值变换•缺点：缺点： 除了零频率附近，除了零频率附近， 与与 之间严重非线性之间严重非线性2 2）要求模拟滤波器的幅频响应为分段常数型，不然）要求模拟滤波器的幅频响应为分段常数型，不然会产生畸变会产生畸变1 1）线性相位模拟滤波器）线性相位模拟滤波器 非线性相位数字滤波器非线性相位数字滤波器分段常数型模拟滤波器经变换后仍为分段常数型数字滤波器，但临界频率点产生畸变F 预畸变预畸变给定数字滤波器的截止频率给定数字滤波器的截止频率 ，则，则按 设计模拟滤波器，经双线性变换后，即可得到 为截止频率的数字滤波器6、模拟滤波器的数字化方法、模拟滤波器的数字化方法四、利用频率变换设计四、利用频率变换设计IIR滤波器滤波器vv方法1：vv方法2：设计模拟低通滤波器从模拟到模拟频率变换实现滤波器数字化脉冲响应不变法双线性变换脉冲响应不变法从数字到数字频率变化实现低通滤波器数字化设计模拟低通滤波器双线性变换数字域频带变换法数字域频带变换法归一化模拟低通数字低通数字低通、高通、带通、带阻数字域频带变换或双线性变换冲激响应不变法要求：阶数：阶数：N N极点：极点：零点：零点：则：则：即：即：为全通函数为全通函数u 数字低通数字低通 — 数字低通数字低通得：得：由主 要 内 容•滤波器概述滤波器概述 •典型模拟滤波器的设计典型模拟滤波器的设计•无限冲激响应无限冲激响应(IIR)数字滤波器的设计数字滤波器的设计è有限冲激响应有限冲激响应(FIR)数字滤波器的设计数字滤波器的设计 FIR数字滤波器的设计方法IIR数字滤波器：可以利用模拟滤波器设计可以利用模拟滤波器设计但相位非线性但相位非线性FIR数字滤波器： 可以严格线性相位，又可任意幅度特性可以严格线性相位，又可任意幅度特性因果稳定系统因果稳定系统可用可用FFTFFT计算计算但阶次比但阶次比IIRIIR滤波器要高得多滤波器要高得多 FIR FIR滤波器的单位冲激响应：滤波器的单位冲激响应：系统函数：系统函数：在在 z z 平面有平面有N N –1–1 个零个零点点在在 z z = 0 = 0 处是处是N N –1–1 阶极点阶极点 4.4 有限长脉冲响应(FIR)数字滤波器设计h h( (n n) )为实序列时，其频率响应：为实序列时，其频率响应：1、线性相位条件即群延时即群延时 是常数是常数第二类线性相位：第二类线性相位：第一类线性相位：第一类线性相位：线性相位是指线性相位是指 是是 的线性函数的线性函数一、 FIR滤波器的特性第一类线性相位：第一类线性相位：n n第一类线性相位第一类线性相位 的充要条件：的充要条件：n n = (= (N N – 1) /2 1) /2 为为h h( (n n) )的偶对称中的偶对称中心心n n第二类线性相位第二类线性相位 的充要条件：的充要条件：n n = (= (N N – 1) /2 1) /2 为为h h( (n n) )的奇对称中心的奇对称中心2、线性相位FIR滤波器频率响应的特点系统函数：系统函数：由n n频率响应：频率响应：1 1））））h h( (n n) )偶对称偶对称偶对称偶对称为第一类线性相位为第一类线性相位相位函数：相位函数：n n频率响应：频率响应：2 2））））h h( (n n) )奇对称奇对称奇对称奇对称相位函数：相位函数：为第二类线性相位为第二类线性相位3、幅度函数的特点1 1））h h( (n n) )偶对称，偶对称，N N为奇数为奇数幅度函数：幅度函数：其中：其中：其中：其中：2 2））））h h( (n n) )偶对称，偶对称，偶对称，偶对称，N N为偶数为偶数为偶数为偶数幅度函数：幅度函数：其中：其中：其中：其中：n n n n 故不能设计成高通、带阻滤波器故不能设计成高通、带阻滤波器 n n 3 3））））h h( (n n) )奇对称，奇对称，奇对称，奇对称，N N为奇数为奇数为奇数为奇数幅度函数：其中：其中：n n n n 4 4））））h h( (n n) )奇对称，奇对称，奇对称，奇对称，N N为偶数为偶数为偶数为偶数幅度函数：其中：其中：其中：其中：n n n n n nh h( (n n) )为奇对称时，有为奇对称时，有90900 0相移，适用于微分器和相移，适用于微分器和90900 0移相器，而选频滤波器采用移相器，而选频滤波器采用h h( (n n) )为偶对称时为偶对称时二、窗函数设计法 1、设计方法w(n)：窗函数序列要选择合适的形状和长度以低通滤波器为例讨论：以低通滤波器为例讨论：以低通滤波器为例讨论：以低通滤波器为例讨论：线性相位理想低通滤波器的频率响应：线性相位理想低通滤波器的频率响应：其理想单位抽样响应：其理想单位抽样响应：中心点为中心点为 α α 的偶对称无限长非因果序列的偶对称无限长非因果序列取矩形窗：取矩形窗：则则FIRFIR滤波器的单位抽样响应：滤波器的单位抽样响应： 加窗处理后对频率响应的影响：加窗处理后对频率响应的影响：加窗处理后对频率响应的影响：加窗处理后对频率响应的影响：时域乘积相当于频域卷积时域乘积相当于频域卷积而矩形窗的频率响应：而矩形窗的频率响应：理想滤波器的频率响应：理想滤波器的频率响应：其幅度函数：其幅度函数：则则FIRFIR滤波器的频率响应：滤波器的频率响应： n n n n n n n n n n n n 幅度函数：加窗函数的影响：加窗函数的影响：加窗函数的影响：加窗函数的影响：n n不连续点处边沿加宽形成过渡带，其宽度（两肩不连续点处边沿加宽形成过渡带，其宽度（两肩峰之间的宽度）等于窗函数频率响应的主瓣宽度。

峰之间的宽度）等于窗函数频率响应的主瓣宽度n n在在 处出现肩峰值，两侧形成起伏振处出现肩峰值，两侧形成起伏振 荡，振荡的幅度和多少取决于旁瓣的幅度和多少荡，振荡的幅度和多少取决于旁瓣的幅度和多少n n改变改变N N只能改变窗谱的主瓣宽度，但不能改变主只能改变窗谱的主瓣宽度，但不能改变主瓣与旁瓣的相对比例其相对比例由窗函数形状瓣与旁瓣的相对比例其相对比例由窗函数形状决定，称为吉布斯（决定，称为吉布斯（GibbsGibbs）效应）效应2、各种窗函数、各种窗函数n n窗函数的要求：窗函数的要求：n n窗谱主瓣尽可能窄以获得较陡的过渡带窗谱主瓣尽可能窄以获得较陡的过渡带n n尽量减少窗谱最大旁瓣的相对幅度尽量减少窗谱最大旁瓣的相对幅度以减小肩峰和波纹以减小肩峰和波纹n n矩形窗矩形窗主瓣宽度最窄：主瓣宽度最窄：旁瓣幅度大旁瓣幅度大窗谱：窗谱：幅度函数： 主瓣宽度宽：主瓣宽度宽：旁瓣幅度较小旁瓣幅度较小窗谱：窗谱：幅度函数： n n三角形（三角形（BartlettBartlett）窗）窗n n汉宁（汉宁（HanningHanning）窗）窗（升余弦窗）（升余弦窗）主瓣宽度宽：主瓣宽度宽：旁瓣幅度小旁瓣幅度小幅度函数： 主瓣宽度宽：主瓣宽度宽：旁瓣幅度更小旁瓣幅度更小幅度函数： n n海明（海明（HammingHamming）窗）窗 （改进的升余弦窗）（改进的升余弦窗）主瓣宽度最宽：主瓣宽度最宽：旁瓣幅度最小旁瓣幅度最小幅度函数： n n布莱克曼布莱克曼BlackmanBlackman））窗（二阶升余弦窗）窗（二阶升余弦窗） ：第一类变形零阶：第一类变形零阶 贝塞尔函数贝塞尔函数n n凯泽（凯泽（KaiserKaiser）窗）窗窗函数窗函数窗谱性能指标窗谱性能指标加窗后滤波器性能指标加窗后滤波器性能指标旁瓣峰值旁瓣峰值/dB/dB主瓣宽度主瓣宽度过渡带宽过渡带宽阻带最小衰减阻带最小衰减/dB/dB矩形窗矩形窗三角形窗三角形窗汉宁窗汉宁窗海明窗海明窗布拉克曼窗布拉克曼窗凯泽窗凯泽窗-13-13-25-25-31-31-41-41-57-57-57-572 24 44 44 46 60.90.92.12.13.13.13.33.35.55.55 5-21-21-25-25-44-44-53-53-74-74-80-80阻带最小衰减只由窗形状决定过渡带宽则与窗形状和窗宽N都有关3、窗函数法的设计步骤、窗函数法的设计步骤n n给定理想的频率响应函数给定理想的频率响应函数及技术指标及技术指标n n求出理想的单位抽样响应求出理想的单位抽样响应n n根据阻带衰减选择窗函数根据阻带衰减选择窗函数n n计算频率响应计算频率响应 ，验算指标是否满足要求，验算指标是否满足要求n n根据过渡带宽度确定根据过渡带宽度确定N N值值n n求所设计的求所设计的FIRFIR滤波器的单位抽样响应滤波器的单位抽样响应公式法：公式法：IFFTIFFT法：法：计算其IFFT，得：对对 M M点等间隔抽样：点等间隔抽样：4、线性相位、线性相位FIR低通滤波器的设计低通滤波器的设计解解：：1 1）求数字频率）求数字频率例：设计一个线性相位例：设计一个线性相位FIRFIR低通滤波器，低通滤波器，给定抽样频率给定抽样频率 ，通带截止频率通带截止频率 ，阻带起始频率为阻带起始频率为 ，，阻带衰减不小于阻带衰减不小于-50-50dBdB，幅度特性如图所示，幅度特性如图所示2 2）求）求h hd d( (n n) )4 4）确定）确定N N 值值3）选择窗函数：由 确定海明窗（-53dB）5 5）确定）确定FIRFIR滤波器的滤波器的h h( (n n) )6 6）求）求 ，验证，验证若不满足，则改变N或窗形状重新设计IIR和和FIR数字滤波器的比较数字滤波器的比较IIRIIR滤波器滤波器FIRFIR滤波器滤波器n nh h( (n n) )无限长无限长n nh h( (n n) )有限长有限长n n极点位于极点位于z z平面任意位置平面任意位置n n滤波器阶次低滤波器阶次低n n非线性相位非线性相位n n递归结构递归结构 n n不能用不能用FFTFFT计算计算n n可用模拟滤波器设计可用模拟滤波器设计n n用于设计规格化的选频用于设计规格化的选频滤波器滤波器n n极点固定在原点极点固定在原点n n滤波器阶次高得多滤波器阶次高得多n n可严格的线性相位可严格的线性相位n n一般采用非递归结构一般采用非递归结构n n可用可用FFTFFT计算计算n n设计借助于计算机设计借助于计算机n n可设计各种幅频特性和可设计各种幅频特性和相频特性的滤波器相频特性的滤波器4.5 与本章内容有关的MATLAB函数w buttord [N，，Wn] = buttord (Wp，，Ws，，Rp，，Rs) [N，，Wn] = buttord (Wp，，Ws，，Rp，，Rs，，‘s’)w buttap [z，，p，，k] = buttap(N) wlp2lplp2lp，，，，lp2hplp2hp，，，，lp2bplp2bp，，，，lp2bslp2bswbilinearbilinear [B，，A] = lp2lp (b，，a，，W0)[B，，A] = lp2hp (b，，a，，W0) [B，，A] = lp2bp (b，，a，，W0，，Bw) [B，，A] = lp2bp (b，，a，，W0，，Bw)[ Bz，，Az] = bilinear (B，，A，，Fs) w wbutterbutterw wfir1fir1 w wfir2fir2 [B，，A] = butter(N，，Wn)[B，，A] = butter(N，，Wn，，’high’)[B，，A] = butter(N，，Wn，，’stop’)[B，，A] = butter(N，，Wn，，’s’)b = fir1 (N，，Wn)b = fir1 (N，，Wn，，’high’)b = fir1 (N，，Wn，，’stop’)b = fir2 (N，，F，，M)Ø滤波器的定义、基本原理，以及滤波器的分滤波器的定义、基本原理，以及滤波器的分类与技术指标；类与技术指标；Ø巴特沃斯和切比雪夫巴特沃斯和切比雪夫Ⅰ型模拟低通滤波器的型模拟低通滤波器的设计方法；设计方法；Ø无限冲激响应数字滤波器和有限冲激响应数无限冲激响应数字滤波器和有限冲激响应数字滤波器的设计方法字滤波器的设计方法小 结。