- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
实验四 数字低通巴特沃斯滤波器的设计
实验四实验目的数字低通巴特沃斯滤波器的设计2?掌握HR数字 滤波器的设计方法
实验原理
IIR数字滤波器的设计,主要采用间接法,即:首先设计出低通 模拟滤波器H⑸;进行频率变换,将其转换为高通、带通、带阻滤波 器;再用脉冲响应不变法或双线性变换法从模拟滤波器转换为数字滤 波器。
对单极点的N阶H⑸用部分展开式:
H(S)??AK
K?1 (S-Sk) N
冲激响应不变法取H(S)的单阶极点Sk的指数函数e
H(z)??NSKFs 作为 H (Z)的极点 AK
(l-eskFsK?lz?l)
2?1?Z?1?双线性变换法是用S??代换H (S)中的S得到H(Z),双线 性变换?1?T?1?Z?
法可完全消除频率混叠失真但存在非线性频率失真,而冲激响应 不变法存在混叠失真。
在不同的设计阶段MATLAB的信号处理工具箱都给出了相应的 滤波器设计函数,这些函数代表了不同类型的逼近函数的滤波器,常 用的有巴特沃斯滤波器、切比雪夫滤波器、椭圆滤波器类型。
设计流程图
下图所示为IIR滤波器设计过程及所用到的MATLAB工具箱中的 函数。
合成一步的设计函数
(butter, chebl, cheb2, ellip, besself)
模拟原型低通滤波器设计
(buttap, cheblap,
cheb2ap, bessllap, ellipap
函数)
求最小阶数N (buttord,
cheblord、 cheb2ord,
ellipord函数)频率转换(低通转低通Ip2lp,低通转高通Ip2hp, 低通转带通Ip2bp,低通转带阻函数Ip2bs)模拟滤波器转数字滤波器
(bilinear, impinvar 滤波器系数 A,B
实验内容
?设计一个低通巴特沃斯模拟滤波器,绘制滤波器的的幅频响应 及零极点图。指标如下:通带截止频率:WP = 1000HZ/通带最大衰减: RP=3dB
阻带截止频率:Ws = 2000HZ,阻带最小衰减:Rs=40dB
参考程序butterl.m
?用冲激响应不变法和双线性变换法将一模拟低通滤波器转换
为数字低通滤波器 并图释H⑸和H(Z),采样频率Fs = 1000Hz
参考程序butter2.m
实验报告
1?列出本次实验编写的所有文件及各项实验结果的曲线,加注必 要的说明频率特性应包括幅频、相频特性
2?理论计算模拟低通原型滤波器的阶数和极点分布,说明实验中
所用的冲激响应不变法和双线性变换法的原理及其编程根据
?总结实验体会及实验中存在的问题
思考:
1?双线性变换法和冲激响应不变法比较有何优点?
IIR设计函数:
巴特沃斯滤波器为例介绍设计函数的功能和用法,其它设计函数 的用法类似。
求最小阶数N的函数buttord
调用格式 1: [N,Wn]=BUTTORD(Wp,Ws,Rp,Rs,#39;s#39;)
?功能:求出巴特沃斯模拟滤波器的最小阶数N和频率参数Wn
?说明:Wp,Ws,为通带截止频率和阻带截止频率如式所示;
Ws?2?fstop
fs, Wp?2?fpassfs,其中fs为时域采样频率
Rp,Rs为通带最大衰减Apass(dB)和阻带最小衰减Astop(dB)如式
Ap?20lgH(ej0)
H(ej2?fpass)(dB)As?20lg
)H(ej0)H(ej2?fstop)(dB))??20lgH(ej2?fpass??20lgH(ej2?fstop
#39;s#39;对应模拟滤波器。
模拟低通滤波器设计函数buttap
调用格式:[ZFK]二BUTTAP(N)
?功能:按给定的阶数N设计出巴特沃斯模拟低通滤波器
?说明:乙P,K为返回的N阶模拟滤波器的零点、极点和增益系数。 如:
H(Z) = k?(Z?zj?l
i?lNj)[式 3.5.4]?(Z?p)i
其中k为增益常数,zj和pi为系统函数H(Z)的N个零点和N个 极点。
模拟数字变换函数一双线性变换函数bilinear或脉冲响应不变 法函数 impinvar 调用格式:[NUMd’DENd]二BILINEAR(NUMQEN’Fs)
?功能:把模拟滤波器系数为NUMQEN变为近似等价的数字滤波 器系数为
NUMd,DENd
调用格式:[BZ,AZ]=IMPINVAR(BAFs)
设计出数字滤波器其单位抽样响应h(n)为模拟滤波器的单位冲 激响应h(t)?功能:
的采样值,采样频率为Fs
B,A为模拟滤波器的系统函数H(S)的分子和分母多项式的系数,
如下所示:?说明:
blsm?b2sm?l??bms?bm?lB(S)H ⑸??l?2???n?n?l
B=[bl/b2,...,bm,bm+l],A=[al,a2,...,an/an+l]
合为一步的数字滤波器设计函数butt
文档评论(0)