第六章IIR数字滤波器设计3详解.ppt

  1. 1、本文档共22页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
第六章IIR数字滤波器设计3详解

6.2.2 用双线性变换法设计IIR数字滤波器 上节介绍的脉冲响应不变法的最大缺点是存在频谱混叠失真。如果 时, ,则数字滤波器产生频谱混叠失真。 双线性变换法采用非线性频率压缩方法,从原理上彻底消除了频谱混叠,所以双线性变换法在IIR数字滤波器的设计中得到更广泛的应用。 * 双线性变换法的基本设计思想 脉冲响应不变法是使模拟滤波器与数字滤波器的冲激响应互相模仿,从而达到两者频响之间互相模仿。 双线性变换法的基本思想是让两种滤波器在输入输出上互相模仿,从而达到频响的互相模仿。 用线性常系数差分方程逼近线性常系数微分方程。 * * * * 设 * 系数关系表 * * 由于数字频率与模拟频率之间的非线性映射关系,典型模拟滤波器幅频响应曲线经过双线性变换后,所得数字滤波器幅频响应曲线有较大的失真。但是,将数字滤波器指标转换成相应的过渡模拟滤波器指标时,只要按照非线性关系式计算模拟滤波器边界频率 (6.2.23) 则过渡模拟滤波器再经过双线性变换后,所得数字滤波器的边界频率就一定满足所要求的数字边界频率指标,幅频响应特性必然满足所给的数字域片断常数幅频响应特性指标。称式(6.2.23)为“预畸变校正”。 * * * * 脉冲响应不变法的优点是频率坐标变换是线性的,如果不考虑频率混叠现象,用这种方法设计的数字滤波器会很好地重现原模拟滤波器的频率特性。一个线性相位模拟滤波器可以映射成一个线性相位的数字滤波器。 另外一个优点是数字滤波器的单位脉冲响应完全模仿模拟滤波器的单位冲激响应,时域特性逼近性好。 缺点是会产生频率混叠现象,适合低通、带通滤波器的设计,不适合高通、带阻滤波器的设计。 双线性变换法避免了频率响应的混叠现象。可以设计高通、带阻滤波器。 模拟频率与数字频率不再是线性关系,所以一个线性相位模拟滤波器经双线性变换后所得到的数字滤波器不再保持原有的线性相位了。 * 下面我们总结利用模拟滤波器设计IIR数字低通滤波器的步骤。 (1)确定数字低通滤波器的技术指标:通带截止频率ωp、通带衰减αp、阻带截止频率ωs、阻带衰减αs。 (2)将数字低通滤波器的技术指标转换成模拟低通滤波器的技术指标。如果使用脉冲响应不变法,则 如果采用双线性变换法,边界频率的转换关系为 (3)按照模拟低通滤波器的技术指标设计模拟低通滤波器。 (4)将模拟滤波器Ha(s),从s平面转换到z平面,得到数字低通滤波器系统函数H(z)。 * 例6.2.3用双线性变换法设计数字低通滤波器,指标要求 要求通带和阻带具有单调下降特性。 * %例6.2.3 close all; clear; wp=0.2; ws=0.35; rp=1; rs=10; [N,wc]=buttord(wp,ws,rp,rs);%[N,wc]=cheb2ord(wp,ws,rp,rs); [B,A]=butter(N,wc) %[B,A]=cheby2(N,rs,wc) [H,W]=freqz(B,A); figure(1) plot(W/pi,20*log10(abs(H))); grid figure(2) plot(W/pi,unwrap(angle(H))); grid * (1) 确定所需类型数字滤波器的技术指标。 (2) 将所需类型数字滤波器的技术指标转换成所需类型模拟滤波器的技术指标。转换公式为 (3)将所需类型模拟滤波器技术指标转换成模拟低通滤波器技术指标。 (4)设计模拟低通滤波器。 (5)将模拟低通通过频率变换,转换成所需类型的模拟滤波器。(6)采用脉冲响应不变法或双线性变换法,将所需类型的模拟滤波器转换成所需类型的数字滤波器。 6.2.3 高通、带通、带阻IIR数字滤波器设计 * 补充例题 用双线性变换法设计数字高通滤波器,指标要求 例6.2.3 改成高通滤波器 close all; clear; ws=0.2; %注意此处修改 wp=0.35; %注意此处修改 rp=1; rs=10; [N,wc]=buttord(wp,ws,rp,rs); [B,A]=butter(N,wc,high) %注意此处修改 freqz(B,A); * close all; clear; fpl=560; fpu=780; fsl=375; fsu=1000; Fs=10000; wp=[2*fpl/Fs,2*fpu/Fs]; ws=[2*fsl/Fs,2*fsu/Fs]; rp=0.5; rs=50; 损耗函数 * [N,wpo]=ellipord(

您可能关注的文档

文档评论(0)

ee88870 + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档