第-10-章-无限脉冲响应滤波器课件.ppt
- 【下载声明】
1. 本站全部试题类文档,若标题没写含答案,则无答案;标题注明含答案的文档,主观题也可能无答案。请谨慎下单,一旦售出,不予退换。
2. 本站全部PPT文档均不含视频和音频,PPT中出现的音频或视频标识(或文字)仅表示流程,实际无音频或视频文件。请谨慎下单,一旦售出,不予退换。
3. 本页资料《第-10-章-无限脉冲响应滤波器课件.ppt》由用户(晟晟文业)主动上传,其收益全归该用户。163文库仅提供信息存储空间,仅对该用户上传内容的表现方式做保护处理,对上传内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知163文库(点击联系客服),我们立即给予删除!
4. 请根据预览情况,自愿下载本文。本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
5. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007及以上版本和PDF阅读器,压缩文件请下载最新的WinRAR软件解压。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 10 无限 脉冲响应 滤波器 课件
- 资源描述:
-
1、 infinite impulse response filter(IIR)无限脉冲响应滤波器 bilinear transformation 双线性变换 prewarping equation 预扭曲方程 Butterworth filter 巴特沃斯滤波器 Chebyshev Type I filter 切比雪夫I 型滤波器 Chebyshev Type II filter 切比雪夫II 型滤波器 elliptic filter 椭圆滤波器 Impulse invariance method 脉冲响应不变法10.1 无限脉冲响应滤波器基础无限脉冲响应滤波器基础 滤波器新的输出,和过去的输出
2、及过去的输入和现在滤波器新的输出,和过去的输出及过去的输入和现在输入有关。输入有关。差分方程差分方程 akyn-k=bkxn-k Nk=0 Mk=0若若 a0=1,则,则 yn=-a1yn-1-aNyn-N +b0 xn+b1n-1+bMxn-M传输函数传输函数H(z)=b0+b1z-1+b2z-2+bMz-M 1+a1z-1+a2z-2+aNz-N 这种递归滤波器的极点由分母多项式确定,这种递归滤波器的极点由分母多项式确定,这就可能使滤波器不稳定,同时递归滤波器这就可能使滤波器不稳定,同时递归滤波器很难实现线性相位,也就是递归滤波器使相很难实现线性相位,也就是递归滤波器使相位失真,但实现某种
3、性能要求时比非递归滤位失真,但实现某种性能要求时比非递归滤波器所需要的系数少。波器所需要的系数少。设计递归滤波器的方法是选择具有待求特设计递归滤波器的方法是选择具有待求特性的原型模拟滤波器,然后将其转换为数字性的原型模拟滤波器,然后将其转换为数字滤波器。滤波器。返回返回10.2 低通模拟滤波器低通模拟滤波器 在连续域中,滤波器用在连续域中,滤波器用 S 而不用而不用 Z 来描述。来描述。简单低通模拟滤波器的传输函数为简单低通模拟滤波器的传输函数为 H(s)=1 s+1频率响应为频率响应为 H()=s用用 j代替代替 1 j+1幅度幅度|H()|=1 2+1|H()|0|H()|1具有低通特性具
4、有低通特性=1弧度弧度/秒秒|H()|=20log()=-3dB 12 12 即即 f=/2=1/2 Hz 为它的带宽为它的带宽为使滤波器更通用,传输函数变为为使滤波器更通用,传输函数变为H(s)=p1 s+p1H()=p1 j+p1 1 1+j/p1其幅度其幅度|H()|=1(/p1)2+1=p1 时,出现时,出现 3dB图图10.3 通常模拟滤波器类型包括:巴特沃斯、切比雪夫通常模拟滤波器类型包括:巴特沃斯、切比雪夫型、型、切比雪夫切比雪夫型和椭圆滤波器等。型和椭圆滤波器等。图图10.4返回返回10.3 双线性变换双线性变换 它为模拟滤波器和数字滤波器之间的转换提供了一种它为模拟滤波器和数
5、字滤波器之间的转换提供了一种方法。方法。s 2fs z 1 z+1 z用用 ej 代替可得到离散时间傅立叶变换代替可得到离散时间傅立叶变换2fs =2fs =2fs ej 1 ej+1 ej/2(ej/2-e-j/2)ej/2(ej/2+e-j/2)ej/2e-j/2 ej/2+e-j/2 sin=和和 cos=eje-j 2j ej+e-j 2将将 s 变为变为 j 则有则有 2fs tan(/2)预扭曲方程预扭曲方程 数字频率数字频率 范围范围 0 弧度弧度 模拟频率模拟频率 范围范围 0 弧度弧度/秒秒 逆双线性变换逆双线性变换 2tan-1(/2fs)图图 10.5 给出了数字频率给出
6、了数字频率 和模拟频率和模拟频率 之间的关系。之间的关系。双线性变换建立了数学域和模拟域之间的另一种联系:双线性变换建立了数学域和模拟域之间的另一种联系:数字滤波器的稳定区域是数字滤波器的稳定区域是 Z 平面单位圆内,模拟滤波器平面单位圆内,模拟滤波器的稳定区域在的稳定区域在 S 复平面的左半部复平面的左半部图图 10.6 例1:一阶模拟低通滤波器的传输函数 H(s)=wp1/(s+wp1)。滤波器的-3dB频率是2000弧度/秒。求出与此模拟滤波器相对应的数字滤波器的传输函数H(z),采样频率为1500Hz.预扭曲的数字滤波器幅度响应如图预扭曲的数字滤波器幅度响应如图 10.10 所示,可所
7、示,可看出预扭曲的滤波器的截止频率非常接近所要求的值看出预扭曲的滤波器的截止频率非常接近所要求的值 318.3 Hz。图图 10.10返回返回10.4 巴特沃斯滤波器设计巴特沃斯滤波器设计 巴特沃斯滤波器是巴特沃斯滤波器是 IIR 无限脉冲响应滤波器中最简单的。无限脉冲响应滤波器中最简单的。它的传输函数它的传输函数 H(s)=p1S+p1|H()|=1(/p1)2+1一阶模拟巴特沃斯滤波器一阶模拟巴特沃斯滤波器 n阶模拟巴特沃斯滤波器具有负责的传输函数阶模拟巴特沃斯滤波器具有负责的传输函数 H(s),其滤波器形状其滤波器形状|H()|=,与一阶相似。,与一阶相似。1(/p1)2n+1阶数越高,
8、滚降越陡,图阶数越高,滚降越陡,图10.12 显示显示 n=1,2,3,4时模拟时模拟巴特沃斯滤波器形状。巴特沃斯滤波器形状。图图10.12描述数字描述数字 IIR 滤波器的术语:滤波器的术语:(由图由图10.13(a)、(b)看看)图图10.13 p 通带波纹:通带内最大和最小增益之差通带波纹:通带内最大和最小增益之差 p 取取 3dB (巴特沃斯滤波器巴特沃斯滤波器)1 p:通带边缘增益(对应的频率通带边缘增益(对应的频率fp1)(或或20log(1 p)dB)s:阻带波纹,对应的阻带边缘增益,:阻带波纹,对应的阻带边缘增益,对应频率对应频率 fs1 阻带边缘频率。阻带边缘频率。阻带边缘的
9、增益为阻带边缘的增益为20logsdB 设计时阻带边缘频率和阻带衰减(设计时阻带边缘频率和阻带衰减(-20logs)或阻)或阻带增益(带增益(20logs)要给出。图)要给出。图 10.14 说明非线性相说明非线性相位响应,表明滤波后的信号在某种程度上有相位失真。位响应,表明滤波后的信号在某种程度上有相位失真。对于对于 n 阶巴特沃斯滤波器,满足通带和阻带设计要求阶巴特沃斯滤波器,满足通带和阻带设计要求所需的阶数所需的阶数 n n log(-1)2log()1s2s1 p1s1:模拟滤波器预扭曲阻带边缘频率。模拟滤波器预扭曲阻带边缘频率。p1:模拟滤波器预扭曲通带边缘频率。模拟滤波器预扭曲通带
10、边缘频率。设计低通巴特沃斯滤波器步骤:设计低通巴特沃斯滤波器步骤:1)确定待求通带边缘频率确定待求通带边缘频率 fp1 Hz、待求阻带边缘频率、待求阻带边缘频率 fs1 Hz 和待求阻带衰减和待求阻带衰减-20logsdB(或待求阻带增或待求阻带增 益益 20logsdB)。通带边缘频率必须对应。通带边缘频率必须对应 3dB增益。增益。4)由已给定的阻带衰减由已给定的阻带衰减-20logs(或增益或增益-20logs)确定阻带边缘增益确定阻带边缘增益 s。3)计算预扭曲模拟频率以避免双线性变化带来的失真。计算预扭曲模拟频率以避免双线性变化带来的失真。由由=2fs tan(/2)求得求得 p1和
11、和 s1,单位是弧度,单位是弧度/秒。秒。2)用式用式=2f/fs 把由把由 Hz 表示的待求边缘频率转成表示的待求边缘频率转成 由弧度表示的数字频率,得到由弧度表示的数字频率,得到 p1 和和s1。5)用下式计算所需滤波器的阶数:用下式计算所需滤波器的阶数:n log(-1)2log()1s2s1 p1 n 取整数取整数6)把把 p1代入代入 n 阶模拟巴特沃斯滤波器传输函数阶模拟巴特沃斯滤波器传输函数H(s)中,并对中,并对 H(s)进行双线性变换得到进行双线性变换得到 n 阶数字传输阶数字传输 函数函数 H(z)。滤波器实现所需的差分方程可直接从。滤波器实现所需的差分方程可直接从 传输函
展开阅读全文