格型滤波器和简单整系数数字滤波器课件.ppt
- 【下载声明】
1. 本站全部试题类文档,若标题没写含答案,则无答案;标题注明含答案的文档,主观题也可能无答案。请谨慎下单,一旦售出,不予退换。
2. 本站全部PPT文档均不含视频和音频,PPT中出现的音频或视频标识(或文字)仅表示流程,实际无音频或视频文件。请谨慎下单,一旦售出,不予退换。
3. 本页资料《格型滤波器和简单整系数数字滤波器课件.ppt》由用户(晟晟文业)主动上传,其收益全归该用户。163文库仅提供信息存储空间,仅对该用户上传内容的表现方式做保护处理,对上传内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知163文库(点击联系客服),我们立即给予删除!
4. 请根据预览情况,自愿下载本文。本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
5. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007及以上版本和PDF阅读器,压缩文件请下载最新的WinRAR软件解压。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 滤波器 简单 系数 数字滤波器 课件
- 资源描述:
-
1、8.2 格型滤波器 8.2.1 全零点(FIR)格型滤波器 一个M阶的FIR滤波器的系统函数H(z)可写成如下形式:()01()()1MMiiiiMiiH zB zbzb z 其中,b(i)M表示M阶FIR滤波器的第i个系数,并假设首项系数b0=1。H(z)对应的格型结构如图 8.2.1 所示。图图 8.2.1 8.2.1 全零点格型滤波器网络结构全零点格型滤波器网络结构x(n)e0r0z1k1k1r1z1r2k2k2e2z1rM 1z1kM 1kM 1eM 1kMkMy(n)e1图图 8.2.2 8.2.2 全零点格型结构全零点格型结构=基本单元基本单元rm1kmkmrmz1em1em 图图
2、 8.2.2 所示基本格型单元的输入、所示基本格型单元的输入、输出关系如下式:输出关系如下式:em(n)=e m-1(n)+r m-1(n-1)km (8.2.2a)rm(n)=e m-1(n)km+rm-1(n-1)(8.2.2b)且 e0(n)=r0(n)=x(n)(8.2.2c)y(n)=em(n)(8.2.2d)格型结构网络系数格型结构网络系数kiki的公式推导:的公式推导:rm1kmkmrmz1em1em k km m与滤波器系数与滤波器系数b b(m)(m)m m之递推关系:之递推关系:()()()11()()()()1121mmiim immmmmmmim iimmmmmbkbb
3、k bkbbk bbk 例例 8.2.1 FIR滤波器由如下差分方程滤波器由如下差分方程给定:给定:1351()()(1)(2)(3)2483y nx nx nx nx n求其格型结构系数,求其格型结构系数,并画出格型结构图。并画出格型结构图。3()123331(1)(2)(3)333(3)331351()()1124831351,248313iiiH zB zb zzzzbbbkb 解:解:对差分方程两边进行对差分方程两边进行Z Z变换的变换的H(z)=BH(z)=B3 3(z)(z):(1)(2)(1)33 3223(2)(1)(2)33 3223(2)22(1)(1)(1)222122(
4、1)111353242481891121211414bk bbkbk bbkkbbk bbkkb图图 8.2.3 H(z)8.2.3 H(z)的格型结构流图的格型结构流图1/4z11/4z11/21/21/31/3z1y(n)x(n)IIR滤波器的格型结构受限于全极点系统滤波器的格型结构受限于全极点系统函数,函数,可以根据可以根据FIR格型结构开发。格型结构开发。设一个设一个全极点系统函数由下式给定全极点系统函数由下式给定:()111()()1MiiMiH zA za z图图 8.2.4 8.2.4 全极点全极点(IIR)(IIR)滤波器格型结构滤波器格型结构k2y(n)x(n)eMrMeM
5、1z 1rM 1kM 1 kM 1z 1e1 k2r2z 1r1e0k1 k1z 1r0kM kMe28.2.2 8.2.2 全极点全极点(IIR)(IIR)格型滤波器格型滤波器H(z)为为A(z)的逆的逆系统系统 例例 8.2.2 设全极点设全极点IIR滤波器系统函数滤波器系统函数为为 求其格型结构网络系数,求其格型结构网络系数,并画出格并画出格型结构。型结构。1231()135112483H zzzz解:解:312()31(1)(2)(3)333135()()1124813513,2483iiMMiBzAzzzb zMbbb 由例由例 8.2.1 所求所求FIR格型结构网络系数:格型结构网
6、络系数:123111,423kkk图图 8.2.5 8.2.5 例例 8.2.2 8.2.2 中的中的IIRIIR格型结构格型结构z11/4z11/21/21/3z1x(n)y(n)1/31/4按逆系统规则,得出:按逆系统规则,得出:在数字信号处理中,格型(在数字信号处理中,格型(Lattice)Lattice)网络起着网络起着重要的作用,它对有限寄存器长度效应敏感度低,重要的作用,它对有限寄存器长度效应敏感度低,在功率谱估计、语音处理、自适应滤波、线性预在功率谱估计、语音处理、自适应滤波、线性预测和逆滤波等方面已经得到广泛应用。测和逆滤波等方面已经得到广泛应用。MATLABMATLAB函数:
7、函数:latc2tf,tf2latclatc2tf,tf2latc8.3 简单整系数数字滤波器 前面介绍的前面介绍的IIRIIR、FIRFIR数字滤波器设计方法可以数字滤波器设计方法可以给出滤波性能相当好的滤波器,但其系数一般给出滤波性能相当好的滤波器,但其系数一般为非整数。在实际应用中,特别是实时信号处为非整数。在实际应用中,特别是实时信号处理场合,有时对滤波器性能要求并不是很高,理场合,有时对滤波器性能要求并不是很高,但对处理速度要求较高,且要求设计方法简单但对处理速度要求较高,且要求设计方法简单易行,这时,简单整系数数字滤波器是最好的易行,这时,简单整系数数字滤波器是最好的选择。选择。概
8、念:概念:指滤波器网络中的乘法支路增益均为指滤波器网络中的乘法支路增益均为 整数的滤波器。整数的滤波器。特点:乘法运算速度快,实际实现的运算单元只特点:乘法运算速度快,实际实现的运算单元只 有少量的移位和相加操作单元。左移一位可实现有少量的移位和相加操作单元。左移一位可实现乘乘2 2运算,左移一位,再加上移位前的原数数据运算,左移一位,再加上移位前的原数数据就可实现乘就可实现乘3 3的运算。其它整数相乘的实现可依的运算。其它整数相乘的实现可依次类推。次类推。整系数数字滤波器整系数数字滤波器 1.1.多项式拟合的基本概念多项式拟合的基本概念 设序列设序列x(n)x(n)中的一组数据为中的一组数据
9、为x(i),i=-x(i),i=-M,:,0,:,MM,:,0,:,M,我们可以构造一个我们可以构造一个p p阶多阶多项式项式f fi i来拟和这一组数据来拟和这一组数据x(i):x(i):2020220,2()()ppkiipkkpMMkikiMiMkfaaia ia ia ipMEfx ia ix i 总的拟合误差为总的拟合误差为(8.3.1)(8.3.2)8.3.1 8.3.1 建立在多项式拟合基础上的简单整系数滤波器建立在多项式拟合基础上的简单整系数滤波器 为了使拟合满足最小均方误差准则,令E对各系数的导数为零,即令002()0,1,()()pMkrkiMkrpMMk rrkkiMiM
10、MrriMMk rk riMEa ix i irpaaix i iFx i isi 则(8.3.3)式可写成如下形式:0prkk rkFa s(8.3.3)(8.3.4)2.最佳拟合模板与简单整系数FIR滤波器的单位脉冲响应h(n)在实际应用中,并不将fi的p+1 个系数全求出来,而是只求出a0,就可实现对x(n)的最佳拟合。由(8.3.1)式可知,例如,当M=2,p=2 时,为五点二次(抛物线)多项式拟合。据(8.3.4)式,并考虑当k+r=奇数时sk+r=0,有 00i ifa(8.3.5)00220204224022020 42s as aFs as aFs Fs Fas ss 其中,代
11、入上式可得 22202402225,10,34,(),()iisssFx i Fx i i03(2)12(1)17(0)12(1)3(2)35xxxxxa(8.3.6)3,12,17,12,3()351724cos6cos2()35jh nH e(8.3.7)(8.3.8)图 8.3.1 低通滤波器幅频特性 (a)M=2,p=2;(b)M=3,p=30H(e j)00.51.00.10.20.50.30.4/(a)0H(e j)00.51.00.10.20.50.30.4(b)2/2 如前所述,如前所述,在单位圆上等间隔在单位圆上等间隔分布分布N个零点,个零点,则构成则构成“梳状滤波器梳状滤波
12、器”。如果在如果在z=1 处再设置一个极点,处再设置一个极点,对消该对消该处的零点,处的零点,则构成低通滤波器,则构成低通滤波器,其系统其系统函数和频率响应函数分别为函数和频率响应函数分别为1(1)/21()11sin(/2)()1sin(/2)NLpj Njj NLpjzHzzeNHeee(8.3.9a)(8.3.9b)8.3.2 8.3.2 建立在零极点对消基础上的简单整系数滤波器建立在零极点对消基础上的简单整系数滤波器 图 8.3.2 低通滤波器零、极点分布及幅频特性(N=10)(a)(8.3.9a)式的零、极点分布图;(b)(8.3.9b)式的幅频特性ImzRez(a)0H(e j)0
13、5100.10.20.50.30.4(b)2 /10/2 基于同样的思想,基于同样的思想,在在z=-1 z=-1 处设处设置一个极点对消该处的零点,置一个极点对消该处的零点,则构成高则构成高通滤波器,通滤波器,其系统函数及频率响应函数其系统函数及频率响应函数分别为分别为 1(1)/2/21()11sin(/2)()1cos(/2)NHPj NjjNHPjzHzzeNHeee(8.3.10a)(8.3.10b)图 8.3.3 高通滤波器零、极点分布及幅频特性 (a)(8.3.10a)式零、极点分布;(b)幅频特性0H(e j)05100.10.20.50.30.4(b)ImzRez(a)2 /1
14、0/2 假设我们要求假设我们要求带通滤波器带通滤波器的中心的中心频率为频率为0,0 0,应当在应当在z=ej0和和z=e-j0处设置一对共轭极点,处设置一对共轭极点,则带通滤则带通滤波器的系统函数和频响函数为波器的系统函数和频响函数为0011120(2)/2/2011()(1)(1)12cos()sin(/2)()coscosNNBPjjjjNBPzzHzezezzzNHee(8.3.11a)(8.3.11b)图 8.3.4 带通滤波器零、极点分布及幅频特性(N=12,0=/6)(a)(8.3.11a)式的零、极点分布;(b)幅频特性曲线0H(e j)070.10.20.50.30.4(b)I
15、mzRez(a)/6654321/2 例如,取理想全通滤波器频响为 HAP(e j)=ce-jm,m为正整数,c为常数 要从HBP(ej)中减去带通滤波器HBP(ej)时,二者的相位特性必须一致。为此,HB P(z)取为如下形式(若取(8.3.11a)式,存在一常数相移/2):00112011()(1)(1)12cosNNBPjjzzHzezezz(8.3.12a)相应的频响函数为(2)/20cos(/2)()coscosjj NBPNHee(8.3.12b)带阻滤波器:带阻滤波器:全通滤波器减去带通滤波器全通滤波器减去带通滤波器 取HAP(ej)中的m=N/2-1即可满足相位特性一致条件,带
16、阻滤波器的系统函数和频响函数分别为(1)2120(2)/201()()()12cos()cos(/2)()()()coscosNNBSAPBPjjjj NBSAPBPzHzHzHzCzzzNHeHeHeCe(8.3.13a)(8.3.13b)11()()1NkLPzHzz(8.3.14)例 8.3.1 设计一个简单整系数低通滤波器,要求f60 Hz时,衰减不大于 3 dB,阻带最大衰减s=40 dB,采样频率fs=1200 Hz。解 由(8.3.9b)和(8.3.14)式知道sin(/2)()sin(/2)kjLPNHe(8.3.15)式中有两个未知数N和k。由已知条件可知:通带边界频率fp=
展开阅读全文