数字信号处理-第六章.ppt
《数字信号处理-第六章.ppt》由会员分享,可在线阅读,更多相关《数字信号处理-第六章.ppt(119页珍藏版)》请在三一办公上搜索。
1、第 6 章IIR数字滤波器设计,Design of IIR Digital Filters,本章主要内容,数字滤波器的基本概念模拟滤波器设计巴特沃思、切比雪夫模拟滤波器设计模拟滤波器的频带变换数字滤波器(Digital Filter,DF)设计用脉冲响应不变法设计IIR数字滤波器用双线性变换法设计IIR数字滤波器IIR数字滤波器的频率变换,滤波的目的为了压制输入信号的某些频率成分,从而改变信号频谱中各频率分量的相对比例滤波技术包括:滤波器设计:根据给定滤波器的频率特性,求得满足该特性的传输函数,滤波过程的实现:获得传输函数后,以何种方式达到对输入信号的进行滤波的目的,6.1 数字滤波的基本概念
2、,图6.1.1 用经典滤波器从噪声中提取信号,数字滤波器输入、输出均为数字信号,通过一定的运算关系,滤除输入信号中某些频率成分或改变输入信号频谱中各频率分量的相对比例的器件具有某种特定频率特性的线性时不变系统,数字滤波器的特点精度高、稳定、灵活 不要求阻抗匹配、实现特殊的滤波器设计数字滤波器的任务寻求一个因果稳定的线性时不变系统,使其系统函数H(z)具有指定的频率特性,1.数字滤波器的分类,经典滤波器(一般滤波器):信号和干扰的频带互不重叠时采用现代滤波器:信号和干扰的频带相互重叠时采用(例如:维纳滤波器、卡尔曼滤波器、自适应滤波器等)经典滤波器的基本功能高通、低通、带通、带阻,图6.1.2
3、理想低通、高通、带通和带阻滤波器幅度特性,数字滤波器分为:无限长单位脉冲响应滤波器(IIR DF)递归系统 有限长单位脉冲响应滤波器(FIR DF)非递归系统,2 数字滤波器的技术指标数字滤波器的频率响应函数H(ej)用下式表示:H(ej)=|H(ej)|ej()|H(ej)|称为幅频特性函数,幅频特性表示信号通过该滤波器后各频率成分振幅衰减情况;()称为相频特性函数;而相频特性反映各频率成分通过滤波器后在时间上的延时情况。两个滤波器幅频特性相同,而相频特性不同,对相同的输入,滤波器输出的信号波形也是不一样的。,理想滤波器不可实现,只能以实际滤波器逼近,通带最大衰减:,阻带最小衰减:,3.滤波
4、器的设计方法,直接设计法在时域或频域直接设计数字滤波器间接设计法先根据指标要求设计对应的模拟滤波器再将模拟滤波器转换为数字滤波器本章介绍间接设计法,6.2 模拟滤波器设计,模拟滤波器设计,模拟滤波器的设计步骤确定指标选择滤波器的类型计算滤波器的阶数查表或计算滤波器的参数,确定其系统函数综合实现及调试滤波器的传输函数(频率响应特性):选频滤波器一般只考虑幅频特性,对相频特性不作要求。幅频特性体现了各频率成分在幅度上的衰减,而相频特性体现的是不同成分在时间上的延时。对输出波形有明确要求时,则需考虑线性相位问题。,6.2.1 模拟滤波器的设计指标,通带边界(截止)频率,阻带边界(截止)频率,3db通
5、带截止频率系统通带和阻带的误差要求通带特性要求:通带最大衰减,或通带峰值波纹,或通带波纹幅度,用分贝表示:,阻带特性要求:阻带最小衰减,或 阻带峰值波纹,或 阻带波纹幅度,用分贝表示:,损耗函数(衰减函数):用来描述滤波器幅频响应特性3dB 通带截止频率c:损耗函数相对放大了小幅度特性:,一般幅频特性,损耗函数,两个附加参数,过渡比或选择性参数反应过渡带的性能,过渡带越窄,k 值趋近于 1对“低通滤波器”:偏离参数 越小、越大,越小,通带、阻带的纹波越小一般有,模拟滤波器的设计模拟滤波器的理论和设计方法已发展得相当成熟,且有若干典型的模拟滤波器供选择。这些滤波器都有严格的设计公式、现成的曲线和
6、图表供设计人员使用。典型滤波器巴特沃斯(Butterworth)滤波器:具有单调下降的幅频特性;切比雪夫(Chebyshev)滤波器:幅频特性在通带或阻带内有波动起伏,这可提高选择性;椭圆(Ellipse)滤波器:幅频特性在通带和阻带都有波动,其选择性最好;贝塞尔(Bessel)滤波器:具有较好的线性相位。,6.2.2 巴特沃斯(Butterworth)模拟低通滤波器设计,巴特沃斯低通滤波器N 阶巴特沃斯低通滤波器的幅度平方函数为:其中N 为滤波器的阶次,为 3dB通带截止频率。,巴特沃斯低通滤波器特点:在 点,的 n(n 2N)阶导数等于零,因此滤波器在 点具有最大平坦幅度;滤波器幅频响应随
7、 的增大而单调下降,因为幅度平方函数的导数小于零;损耗函数,滤波器的特性由 3dB 截止频率 和阶数 N 确定滤波器的给定指标为通带边界频率(或3dB通带截止频率)阻带边界频率通带最小幅度(或通带最大衰减)阻带最大波纹(或阻带最小衰减),3dB截止频率与阶数如何确定?,滤波器幅频响应随频率的增大而单调下降于是取整数 N:,若取,则满足通带指标,阻带指标有富裕若取,则满足阻带指标,通带指标有富裕,若滤波器的给定指标为通带最大衰减阻带最小衰减则先求然后确定阶数 与3dB通带截止频率,或,模拟滤波器系统函数的求取,已知,则巴特沃斯低通滤波器幅度平方函数确定。定义幅度平方函数 Ha(s)模拟滤波器 系
8、统函数Ha(j)滤波器的频率响应|Ha(j)|滤波器的幅频响应,极点和零点总是“成对出现”,共轭对称;对称于S平面的实轴和虚轴,选用 的对称极、零点的任一半作为 的极、零点,则可得到。为了保证 的稳定性,应选用 在S平面左半平面的极点作为 的极点。零点的分布则无此限制,只和滤波器的相位特性有关,如果要求是最小相位延迟特性,则 应取左半面零点,若无特殊要求,则可将对称零点的任一半(为共轭对)取为 的零点。,确定的方法如下:由 得到象限对称的 S 平面函数;将 因式分解,得到各零极点,将左半平面极点归于。零点(或极点)都为偶次,应取一半(应为共轭对)零点归于 的零点。按照 与 的低频或高频特性的对
9、比就可以确定出增益常数。由求出的零点,极点及增益常数,则可完全确定系统函数,例 根据以下幅度平方函数 确定系统函数,解:,Step1:,Step2:,Step3:,将幅度平方函数转换成s域函数,(令 代入)有令分母为零,得 式中,k0,1,2,(2N1),巴特沃斯低通滤波器的系统函数,可见,Butterworth滤波器 的振幅平方函数有2N个极点,它们均匀对称地分布在|s|=c的圆周上。(没有零点)例:N=3 阶BF振幅平方函数的极点分布,如图。,考虑到系统的稳定性,系统函数是由 S 平面左半部分的极点(SP3,SP4,SP5)组成的,它们分别为:系统函数为:或 令,则有,归一化的三阶滤波器的
10、系统函数,巴特沃斯低通滤波器的系统函数,三种形式归一化函数(,):,去归一化:,极点,分母多项式系数,分母因式,低通巴特沃斯滤波器设计步骤:,由 求滤波器阶次 N 和;由 N 查表,求出归一化系统函数 G(p)(或归一化极点)(pp 157)将 G(p)中的 p 代换为(去归一化),得实际滤波器传输函数 Ha(s)。,滤波器的设计步骤:,根据技术指标求出滤波器阶数N:,确定技术指标:,由,得:,同理:,令,则:,求出归一化系统函数:,或者由N,直接查表得,其中技术指标Wc给出或由下式求出:,其中极点:,去归一化,阻带指标有富裕,或,通带指标有富裕,例1 已知通带截止频率fp=5kHz,通带最大
11、衰减p=2dB,阻带截止频率fs=12kHz,阻带最小衰减s=30dB,按照以上技术指标设计巴特沃斯低通滤波器。解:(1)确定阶数N。,(2)其极点为,归一化传输函数为,上式分母可以展开成为五阶多项式,或者将共轭极点放在一起,形成因式分解形式。这里不如直接查表简单,由N=5,直接查表得到:极点:-0.3090j0.9511;-0.8090j0.5878;-1.0000,式中 b0=1.0000,b1=3.2361,b2=5.2361,b3=5.2361,b4=3.2361,(3)为将Ha(p)去归一化,先求3dB截止频率c。,得到:,【将c代入,得到:,这里算得的s比给定的指标要小一些,可见实
12、际设计所得的阻带指标有富裕量。】将p=s/c代入G(p)中得到:,例2:已知通带边界频率,通带最大衰减,阻带截止频率,阻带最小衰减,设计巴特沃斯低通滤波器。解:1)确定阶次 N:可得 N按照以上技术指标设计巴特沃斯低通滤波器。,查表求分母因式:系统函数,例:设计Butterworth低通滤波器,要求在频率低于0.2p rad/s的通带内幅度特性下降小于1dB。在频率0.3p rad/s之后的阻带内,衰减大于15dB。,1)模拟滤波器的技术指标,a)确定参数,用通带技术指标,使阻带特性较好,改善混迭失真,2)设计Butterworth模拟低通滤波器,b)求出极点(左半平面),c)构造系统函数,c
13、)去归一化,6.2.3 切比雪夫(chebyshev)滤波器,切比雪夫型通带内等波纹,阻带内单调下降切比雪夫型通带内单调下降,阻带内等波纹巴特沃兹滤波器在通带内幅度特性是单调下降的。巴特沃兹滤波器如果阶次一定,则在靠近截止频率处,幅度下降很多,或者说,为了使通带内的衰减足够小,需要的阶次(N)很高。切比雪夫型滤波器的波纹在通带范围内是等幅起伏的,同样的通带衰减,其阶数较巴特沃兹滤波器要小。,切比雪夫型的幅度平方函数:通带边界频率通带波纹参数(为小于1的正数)N 阶切比雪夫多项式:切比雪夫多项式的递推公式:,切比雪夫滤波器的振幅平方特性,通带内 变化范围 p,随,当=0时,N为偶数,min,N为
14、奇数,max,,阻带边界频率,N 大于等于上式的最小整数,滤波器的系统函数,实际应用中可以利于查表的方法求得系统函数:,1.按 求出滤波器阶次N;2.由N和(或p)查表得到归一化滤波器分母多项式:3.归一化切比雪夫低通滤波器函数为:4.去归一化,得到通带截止频率为p的低通滤波器函数:,切比雪夫型通带内单调下降,阻带内等波纹幅度平方函数:阶数 N 的计算,系统函数,椭圆滤波器(考尔滤波器),幅值响应在通带和阻带内都是等波纹的对于给定的阶数和给定的波纹要求,椭圆滤波器能获得较其它滤波器更窄的过渡带宽,就这点而言,椭圆滤波器是最优的。通带和阻带内波纹固定时,阶数越高,过渡带越窄;阶数固定,通带和阻带
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数字信号 处理 第六
链接地址:https://www.31ppt.com/p-6294337.html