- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
* * 由上图 d)可见,这时阻带衰减达到了-43dB。为了进一步增加阻带衰减,可再增加一个过渡采样点,并将采样点数增加到60,两个过渡样点值经优化分别为0.5925和0.1099, 相应的 和 分别为 如下图所示。 过渡带为 ,求得的单位脉冲响应 和幅度响应 * * 由上图 d)可见,这时阻带衰减达到了-63dB。还可以通过进一步增加过渡样点来增加阻带衰减, 显然,在保证过渡带宽不变的情况下,相应的采样 点数也就随着成倍的增加,这样将使滤波器的复杂 度大大增加,在实现滤波时计算量也就随着增加。 频率采样法可以直接从频域出发进行设计,比较简单。这种方法最适合于设计窄带滤波器。 与窗函数法比较 1、从性能上进行比较 6.5 FIR滤波器IIR滤波器的比较 IIR滤波器 传输函数的极点可位于单位圆内的任何地方,因此可 用较低的阶数获得高的选择性,所用的存储单元少,所以 经济而效率高。但是这个高效率是以相位的非线性为代价 的,选择性越好,则相位非线性越严重。 FIR滤波器 可以得到严格的线性相位,然而由于FIR滤波器传输函 数的极点固定在原点,所以只能用较高的阶数达到高的选 择性;对同样的滤波器设计指标,FIR滤波器所要求的阶数 可以比IIR滤波器高5-10倍,结果,成本较高,信号延时也 较大。 2、从结构上进行比较 IIR滤波器 必须采用递归结构,极点位置必须在单位圆内,否则 系统将不稳定。另外,在结构中,由于运算过程中对序列 的舍入处理,这种有限字长效应有时会引入寄生震荡。 FIR滤波器 主要采用非递归结构,不论在理论上还是在实际的有 限精度运算中都不存在稳定性问题,运算误差也较小。此 外,FIR滤波器可以采用快速傅立叶变换算法,在相同阶数 的条件下,运算速度可以快很多。 3、从设计工具进行比较 IIR滤波器 可以借助于模拟滤波器的成果,因此一般都有有效的 封闭形式的设计公式可供准确计算,计算工作量比较小, 对计算工具的要求不高。 FIR滤波器 设计则一般没有封闭形式的设计公式,窗口法虽然仅 仅对窗口函数可以给出计算公式,但计算通带阻带衰减等 仍无显式表达式。一般,FIR滤波器的设计只有计算程序可 循,因此对计算工具要求较高。 4、其他 IIR滤波器 虽然设计简单,但主要是用于设计具有片段常数特性 的滤波器,如低通、高通、带通及带阻等,往往脱离不了 模拟滤波器的格局。 FIR滤波器 设计灵活,尤其它能易于适应某些特殊的应用,如构 成微分器或积分器,或用于巴特沃思、切比雪夫等逼近不 可能达到预定指标的情况,例如,由于某些原因要求三角 形振幅响应或一些更复杂的幅频响应,因而有更大的适应 性和更广阔的天地。 总结 IIR与FIR滤波器各有所长,所以在实际应用时应该从 多方面考虑来加以选择。例如,从使用要求上来看,在对 相位要求不敏感的场合,如语言通讯等,选用IIR较为合 适,这样可以充分发挥其经济高效的特点,而对于图像信 号处理,数据传输等以波形携带信息的系统,则对线性相 位要求较高,如果有条件,采用FIR滤波器较好,当然, 在实际应用中应考虑经济上的要求以及计算工具的条件等 多方面的因素。 * ⑷ 滤波器长度N一般由线性相位的斜率 决定,当 未给定时,N可由过渡带确定,N确定后, 也就确定了。它们关系为 。 此题过渡区宽度要求 汉宁窗设计的滤波器过渡带宽度为 ,则 选 ,则 ⑶ 根据阻带指标,查表可知汉宁窗,海明窗和布莱 克曼窗都满足阻带40dB的衰减 * ⑸ ,N和 确定后, 也就确定了: ⑹ 求 : * * %根据给定的参数指标设计fir滤波器 clear all f=[.3 .46]; a=[ 1 0]; Ap=3; As=40; dp=10^(Ap/(-20)); ds=10^(As/(-20)); dev=[dp ds]; [N,Wn,beta]=kaiserord(f,a,dev); %求开赛窗函数的参数 b1=fir1(N,Wn,kaiser(N+1,beta)); %参数window必须是一个长度为N+1的列向量 wp=0.3;ws=0.46; wn1=(ws+wp)/2; n1=ceil(8*pi/(wn1*pi)); b2=fir1(n1,wn1,hann(n1+1)); %求Hanning窗函数的参数 b3=fir1(n1,wn1,hamming(n1+1)); %
文档评论(0)