第4章-滤波器设计方法课件.ppt
- 【下载声明】
1. 本站全部试题类文档,若标题没写含答案,则无答案;标题注明含答案的文档,主观题也可能无答案。请谨慎下单,一旦售出,不予退换。
2. 本站全部PPT文档均不含视频和音频,PPT中出现的音频或视频标识(或文字)仅表示流程,实际无音频或视频文件。请谨慎下单,一旦售出,不予退换。
3. 本页资料《第4章-滤波器设计方法课件.ppt》由用户(晟晟文业)主动上传,其收益全归该用户。163文库仅提供信息存储空间,仅对该用户上传内容的表现方式做保护处理,对上传内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知163文库(点击联系客服),我们立即给予删除!
4. 请根据预览情况,自愿下载本文。本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
5. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007及以上版本和PDF阅读器,压缩文件请下载最新的WinRAR软件解压。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 滤波器 设计 方法 课件
- 资源描述:
-
1、4.1 引言引言数字滤波器的总体设计步骤包括:(1)按照设计要求,确定滤波器的性能指标;(2)采用一个因果、稳定的离散LTI系统的系统函数逼近该性能指标;(3)利用有限精度算法来实现该系统函数;(4)采用通用计算机软件、专用或通用的数字信号处理器(DSP)来实现。本章主要讨论其中的第(2)项和第(3)项中的运算结构等内容 4.2 模拟滤波器的设计模拟滤波器的设计 按照频率特性,模拟滤波器可分为低通、高通、带通、带阻和全通等类型。模拟低通滤波器的传输函数经频率变换可以转换成模拟高通、带通和带阻滤波器。设计模拟低通滤波器是根据一组设计规范求得模拟系统函数,使其逼近理想低通滤波器的特性。通常可根据滤
2、波器的幅度平方函数进行逼近,也可根据相位特性或群延迟特性进行逼近。模拟低通滤波器的幅度平方逼近函数有巴特沃斯型、切比雪夫型和椭圆型函数等。4.2.1 低通模拟滤波器特性及幅度 平方逼近 1理想滤波器的幅频特性 上图给出了各种理想滤波器的幅频特性,这些 特性包括:(1)通带内传输函数的幅度是常数,相移特性 是线性的;(2)阻带内传输函数的幅度为零,对相移特性 没有要求;(3)从通带到阻带有一个过渡带,对于理想滤 波器为零。2.低通滤波器的技术指标及其幅频特性对于模拟低通滤波器,技术指标包括:(1)通带截止频率 及通带内的最大衰减(或通带边缘增益 )或通带波纹 ;(2)阻带截止频率 及阻带内最小衰
3、减(或阻带边缘增益 )或阻带波纹 。ppapa-pssasa-s 对于单调下降的幅度特性,如果 且 ,则 称为3dB截止频率或半功率截止频率。通带波纹 用来描述通带内最大和最小增益之差,即 。衰减为 。可得 与 的关系为:2)1(1log10ppacdB3pacp)j(1appH)1log(20ppap/2H1j2ca 阻带波纹 用来描述阻带内的最大增益,即 。衰减为 dB,可得 与 的关系为:21log10ssa)j(assH sslog20sas 给定技术指标后,模拟滤波器的设计任务就是求得一个系统函数 ,使其满足该技术指标。)(sHa3由幅度平方函数确定系统函数模拟滤波器的幅度响应常用幅
4、度平方函数来表示,即)j()j()j(aa2aHHH 由于滤波器冲激响应是 实信号,有 ,所以)j()j(aaHH)t(ahjsaaaa2a)s()s()j()j()j(HHHHH 若已知 ,如何求得 呢?2a)j(H)s(aH 由于 为实信号,其系统函数中的极点(或零点)必以共轭对形式出现。又由于实际可实现的滤波器都是稳定的,其系统函数 的极点一定位于s左半平面,而 的极点一定位于s右半平面,即左半平面的极点属于 ,而右半平面的极点属于 。)t(ah)s(aH)s(aH)s(aH)s(aH4.2.2 巴特沃斯滤波器设计 1巴特沃斯滤波器的特性巴特沃斯低通滤波器的幅度平方函数定义为:N22a1
5、1)j(Hc式中,N为正整数,表示滤波器的阶数,为3dB截止频率。c 应用巴特沃斯幅度平方函数描述的低通滤波器具有如下特点。(1)当 时,。01)0 j(2aH(2)当 时,或 。c21)j(2acH707.021)j(caHdB3)j(lg20caH(4)当 时,即过渡带和阻带内,也随着 的增大单调减小,但比通带内衰减的速度快得多,且N越大,衰减速度越快,过渡带也就越陡峭,如图所示。c2a)j(Hc(3)当 时,即通带内,随着 的增加单调减小,且N越大,实际越慢,通带也就越平坦。c02)(jHac巴特沃斯滤波器幅度特性及其与N的关系2巴特沃斯滤波器的系统函数 N2caa2jsajs11)s(
6、)s()j(HHH将 代入式中,可得:js 1N2,2,1,0k,ej1sN21k221jccN21k 可见,巴特沃斯滤波器是一个全极点型滤波器。共有2N个极点,每个极点表示如下:由上可看出,的极点分布特点为:)s(H)s(Haa (1)极点在s平面是象限对称的,分布在半径为 的 圆(巴特沃斯圆)上,共有2N个极点。极点间的角度间 隔为 。cradN (2)极点绝不能落在虚轴上,这样滤波器才有可能稳 定。(3)N为奇数时,实轴上有极点;N为偶数时,实轴上 无极点。选取在左半平面上的极点作为 的极点,因此有:1N0kkNca)ss()s(H sHa例例4-1 导出三阶巴特沃斯型低通滤波器 的系统
7、函数,已知s/rad2c解解 该巴特沃斯型低通滤波器的幅度平方函数为:62a211)j(H66aa2s11)s()s(HH即即 5,2,1,0k,e2s61k221jk各极点为:3j1e2s32j02e2s j13j1e2s32j23j1e2s31j32e2s2j43j1e2s31j5即 极点的分布如图所示。取s平面左半平面的极点 、和 组成 :0s1s2s)s(aH8s8s4s8)ss)(ss)(ss()s(232103caH3巴特沃斯模拟滤波器的设计方法 确定技术指标,包括:通带截止频率 及通带内的最大衰减 ,阻带截止频率 及阻带内最小衰减 。ppassa根据式得到:N2cp2pa11)j
8、(HN2cs2sa11)j(H从而得到:1101101010Nspspaa则 sp1010lg110110lgNspaa 当 时,简化得:cpps2slg211lgN例例4-24-2 已知一个模拟低通滤波器通带截止频率 通带最大衰减 阻带起始频率 ,阻带最小衰减 ,试根据上述要求设计巴特沃斯低通滤波器。kHz5pfdB3pakHz10sfdB30sa解解 (1)确定滤波器的阶数 因为 ,由 及 可求得 cp2sas)j(1lg10Ha)j(sasH 0316.01010203020ssas则 9829.46021.00002.32lg210316.01lglg211lgN2ps2s取阶数N5。
9、因为 所以dB3pas/rad10523pc(2)计算s左半平面上的N个极点53jc0es 54jc1es jc2es 56jc3es 57jc4es(3)写出滤波器的系统函数2218214394552240kk5ca1006.3s10152.3s10624.1s10168.5s10017.1s1006.3)ss()s(H4.2.3 切比雪夫切比雪夫I型滤波器设计型滤波器设计 1切比雪夫滤波器的特性 若切比雪夫滤波器的幅度特性在通带内是等波纹的,而在阻带内是单调的,则称之为切比雪夫I型滤波器,如图所示 若幅度特性在通带内是单调的,而在阻带内是等波纹的,则称之为切比雪夫II型滤波器。这里主要介绍
10、切比雪夫I型滤波器的设计。N阶切比雪夫I型滤波器形状N阶模拟切比雪夫I型滤波器,其幅度平方函数定义为:p2N22aC11)j(H 式中,参数 为小于1的正数,为低通滤波器的截止频率(不一定是3dB频率)。是N阶切比雪夫多项式,定义为:p)x(NC1x)x(harccosNcosh1x)xarccos(Ncos)x(NC可知,幅度平方函数由阶数N和参数 确定。2切比雪夫模拟滤波器的设计方法 模拟切比雪夫I型滤波器的阶数N可由下式确定:psharccosharccosN式中,由阻带边缘增益 (如图所示)计算得到:112s参数 由通带边缘增益计算得到:1)1(12p则系统函数为:线性增益频谱图 ii
11、sssH11a2 例例4-3 一个模拟低通滤波器,要求通带截止频率为 3kHz,通带最大衰减为0.1dB,阻带截止频率为12kHz,阻带最小衰减为60dB。应用切比雪夫I型滤波器的设计 方法求该滤波器的幅度函数。解解由 ,求得 ,则60lg20s001.0102060s10001001.011122s1523.019886.011)1(122p1.01lg20p 由 ,求得 。则9886.0101201.0p所需滤波器的阶数为:5957.40634.24828.9312harccos1523.01000harccosharccosharccosNps取N=5,则该模拟低通滤波器的幅度表示为:2
12、1523.01111j252p22aCCHN4.3 无限冲激响应(IIR)数字滤波器设计 无限冲激响应滤波器也称为递归滤波器,其所对应的差分方程如下:00kMkkkknxbknya一般IIR滤波器的传输函数为:NN2211MM22110zzz1zzz)z(Haaabbbb具有M个零点和N个极点,可见递归滤波器不能保证其稳定性。所以,稳定性检验是许多IIR滤波器设计的重要组成部分。递归滤波器很难实现非递归滤波器所具有的线性相位;递归滤波器的优势在于,在实现类似的性能要求时,递归滤波器比非递归滤波器所需要的系数(或阶数)要小得多。设计IIR数字滤波器一般有以下两种方法。递归滤波器与非递归滤波器之间
13、的区别还包括:(1)先按照技术指标设计一个模拟滤波器 ,然后再按照一定的映射关系将 转换成数字滤波器的 。(2)计算机辅助设计法。这是一种最优化的设计方法。这种设计需要进行大量的迭代运算,故需要计算机的支持。)s(aH)s(aH)(zH 我们主要介绍第一种方法,即利用模拟滤波器设计数字滤波器。那么,如何将模拟滤波器的传输函数 转换成数字滤波器的传输函数 呢?)s(aH)(zH 模拟滤波器到数字滤波器的转换可以在时域实现,也可以在频域实现。时域转换法是使数字滤波器的时域响应与模拟滤波器的时域采样值相等,具体方法有:冲激响应不变法、阶跃响应不变法和匹配z变换法。频域变换法是使数字滤波器在范围内的幅
14、度特性与模拟滤波器在范围内的幅度特性一致,具体方法有:双线性变换法、微分映射法。上述方法中得到广泛使用的有冲激响应不变法和双线性变换法。无论采用哪种方法,为了保证由 转换得到的 仍然满足技术指标并具有稳定性,复变量s到复变量z之间的映射关系必须满足以下两个基本要求:)s(aH)(zH(1)的频率响应要能够模仿 的频率响应。)(zH)s(aH)s(aH)(zH(2)因果稳定的 能够映射成因果稳定的 。4.3.1 冲激响应不变法 1 变换原理 冲激响应不变法是依据数字滤波器的冲激响应 与模拟滤波器的冲激响应 在采样点上值 相等,即 )(nh)t(ahnTta)t(h)nT()t()n(anTtah
15、hh 如果令 是 的拉普拉斯变换,是 的z变换,则由模拟系统的系统函数 求拉普拉斯逆变换得到模拟的冲激响应 ,然后采样得到 ,再取z变换得 。s平面和z平面之间的对应关系:)s(aH)t(ah)(zH)n(h)s(aH)t(ah)nT()n(ahh)(zHTesz 下图显示了冲激响应不变法的s域到z域的映射关系。若将 和 代入 ,可得 jrezjsTer T 分析结果表明:(1)当 ,即 时,即s平面的虚轴映射到z平面的单位圆上;(2)当 时,即s平面的左半平面映射到z平面的单位圆内;(3)当 时,即s平面的右半平面映射到z平面的单位圆外。js 01z01z01zTesz 分析图不难发现,从s
16、平面到z平面是一种多值映射关系,因此,要求 必须在 内严格带限,且带限于折叠频率以内时,即 2 混叠失真混叠失真)j(HTT时才能使数字滤波器的频率响应在折叠频率以内重现模拟滤波器的频率响应而不产生混叠失真,即,)Tj(T1)e(ajHH 但是,任何严格设计的模拟滤波器的频率响应都不是严格带限的,变换后会产生周期延拓分量的频谱交叠。0jaH,2s 由于 是实数,因而 的极点必成共轭对存在,则 变成 的映射关系为:,)t(ah)s(aH)s(aH)(zH3 模拟滤波器数字化kskkkzAssAe11ks1e1zAssAkkk例例4-4 已知模拟滤波器的系统函数为试用冲激响应不变法求出相应的数字滤
17、波器的系统函数。解22a)s(s)s(baaHbabaHjs21js21)s(a根据冲激响应不变法中z平面和s平面的映射关系可知,相应的数字滤波器在 处有一对极点,其系统函数为:这是一个二阶递归滤波器,在 处有一对共轭极点,在 和 处有两个零点。T)j(ebaz1TjT1TjT1T1TjT1TjTzee1zee1z)Tcos(e1zee121zee121)z(babaabababH()ajb Tze0z)Tcos(ezaTb将进行部分分式展开得4 冲激响应不变法的优缺点 冲激响应不变法的最大优点是保持了模拟滤波器的时域瞬态特性,具有良好的时域逼近,而且模拟频率和数字频率之间呈线性映射关系。冲激
18、响应不变法的主要缺点是有频域混叠,所以它只适用于带限的模拟滤波器。4.3.2 双线性变换法 冲激响应不变法是使数字滤波器在时域上模仿模拟滤波器,但是由于从s平面到z平面是多值映射关系,所以产生了频率响应的混叠失真。为了克服这一缺点,可以采用双线性变换法。1变换原理 )e(jH)j(HT2TTTs1ez 双线性变换法是使数字滤波器的频率响应 与模拟滤波器的频率响应 相近似的一种变换方法。它通过将整个s平面压缩变换到某一中介的s1平面的一条横带里(宽度为 ,即 ),然后再通过标准变换关系 将此横带变换到整个z平面上去,如图所示,4.3.2 双线性变换法双线性变换法的映射关系将s平面整个 轴压缩到s
19、1平面轴 上的 一段上,可以通过以下的正切映射来实现:j1jTT11jTjT1e1e1c)2Ttan(jcj从而得到s平面和z平面的单值映射关系为:1z1zcz1z1cs11scscz则这两式就是s平面与z平面之间的单值映射关系,这种变换方法称为双线性变换法。(1)(2)将 和 代入式得到模拟角频率 和数字角频率 之间的变换关系:js jez 2tanc 上式的逆双线性变换式为:ctan21下面讨论双线性变换法对变换条件的可满足性:(1)将 代入式(1),可得即s平面的虚轴映射成z平面的单位圆。(2)将 代入式(2),可得 因此有 j2tanjce1e1csjjj)c(j)c(scscz222
20、2)c()c(zjezjs 可见,当 时,;当 时,;当 时,。因此,稳定的模拟滤波器经过双线性变换后,所得的数字滤波器也一定是稳定的。01z01z01z2.频率预畸变 2tan c 双线性变换的最大优点是避免了频率响应的混叠现象。但是,由于 变换关系是非线性的,所以 和 之间存在着严重的非线性关系,这种非线性关系要求模拟滤波器的幅度响应必须是分段恒定的。对应分段恒定的滤波器,经过双线性变换后,仍得到幅频特性为分段恒定的滤波器,但是各个分段边缘的临界频率点产生了畸变。这种现象可以通过频率的预畸变来加以校正,也就是将临界频率事先加以畸变(利用公式 ),然后再经过双线性变换后映射到所需要的频率上。
21、2tan c 3.模拟滤波器数字化模拟滤波器数字化 )(asH)(zH 模拟滤波器的数字化可以由模拟滤波器的系统函数 通过式(1),得到数字滤波器的系统函数 ,即 11a11a1111zzcHsHzHzzcs 也可以先将模拟系统函数分解成并联的子系统或级联的子系统,使每个子系统函数都变成低阶的(如一、二阶),然后再对每个子系统函数分别采用双线性变换法。例例4-5 已知模拟滤波器的系统函数为 ,试利用双线性变换法将其转换成数字滤波器 。ssHa11)()(zH解解1111111112121111211112)()(11zfzzzzfzzfHsHzHsssazzfsas1131)(zzzH当T=1
22、或 时,则1sf1 冲激响应不变法设计数字低通滤波器 步骤如下:4.3.3 4.3.3 低通数字滤波器设计低通数字滤波器设计()确定待求通带边缘频率 、待求阻带边缘频率 和待求阻带边缘增益 ,通带边缘频率一般对应3dB增益。()用 把以Hz为单位表示的待求边缘频率转换成以rad为单位的数字频率,得到 和 。()模拟频率 和数字频率 之间的转换关系为线性关系,即:,或 ,求得 和 ,单位是rad/s。()由已知的阻带边缘增益 dB,确定阻带边缘增益 。1pfs1fslg20s2ff psT sfTpsslg20s()用下式计算所需滤波器的阶数。()将代入N阶巴特沃斯滤波器的传输函数中,求出的N个
23、极点。为整数NNpss,lg211lg2 2aa2jaj11jssHsHHs(7)用下面公式求得所需数字滤波器的传输函数。NkTskzeAzHk111)(1akkkssAsH2 双线性变换法设计数字低通滤波器 步骤如下:()确定待求通带边缘频率 、待求阻带边缘频率 、待求的通带边缘增益 dB和待求阻带边缘增益 dB。(注意:这时,通带边缘频率不必一定对应-3dB增益。)p1fs1fp1lg20slg20()用公式 把以Hz为单位的待求边缘频率转换成以rad为单位的数字频率,得到 和 。s2ff ps()计算预扭曲模拟频率,以避免双线性变换带来的失真。由 求得 和 ,单位是rad/s。2tan2
24、fsps(4)由指定的通带边缘增益 dB,确定通带边缘增益 ,并由下式计算参数 :(5)由指定的阻带边缘增益 dB,确定阻带边缘增益 ,并用下式确定 。p1lg20p11)1(12pslg20s112s(6)用下式确定所需滤波器的阶数N。(7)把 和 代入N阶切比雪夫I型滤波器的传输函数 中(由于系统函数通常非常复杂所以一般借助于滤波器设计软件来完成),并对 进行双线性变换得到N阶数字系统函数 。psN11coshcoshp)(sHa)(sHa)(zH(8)把 代入下式即可得到数字滤波器幅频特性 。2tan2fs 112azzfsssHzHjeH22a11jCH4.3.4 高通、带通和带阻高通
25、、带通和带阻IIR数字滤波数字滤波 器设计器设计 实际应用中的高通、带通和带阻数字滤波器的设计主要有如图所示的三种方法。数 字 化 模 拟 域 频 率 变 换 模 拟 归 一 化 原 型 c 1 模 拟 高 通、带 通、带 阻 数 字 低 通、高 通、带 通、带 阻 (a)方法一:先进行模拟域频率变换,再数字化 数字域频率变换 数字化 模拟归一化原型 c1 数字低通 数字低通、高通、带通、带阻(b)方法二:先数字化,再进行数字域频率变换 数字化 模拟归一化原型 c1 数字低通、高通、带通、带阻(c)方法三:直接对模拟原型进行数字化1.方法一:先进行模拟域频率变换,再数字化 设计步骤如下:(1)
展开阅读全文