用MATLAB语言实现IIR滤波器的设计.doc

  1. 1、本文档共9页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
第17卷第4期 2001年12月齐齐哈尔大学学报Journal of Q i q ihar Universit y Vol. 17,No. 4Dec. ,2001 用MATL AB 语言实现IIR 滤波器的设计 柳春锋 (齐齐哈尔大学自动化系, 齐齐哈尔161006 摘要本文通过对双线性变换法的分析, 着重探讨了用MA TLAB 语言设计IIR 滤波器的软件方法及应注意的问题。 关键词:IIR 滤波器;MA TLAB 语言; 设计; 方法 中图分类号:TN 713+. 1; TP312文献标识码:A文章编号:1007-984X (2001 04-0051-04 在信号处理中, 软件实现大多使用BASIC 、FOR TRAN 和C 语言, 这在处理如FF T 等复杂的数学问题和输出图形时有许多不便。MA TAL AB 是一种可视化的功能强大的系统分析和仿真工具, 其工具箱中丰富的函数和作图功能正是其他语言所缺乏的。本文在分析IIR 滤波器设计法的基础上, 用MA TAL AB 语言实现了IIR 滤波器的设计, 通过修改某些参数和比较结果, 提出了IIR 滤波器设计中应注意的问题。1设计思想 IIR 滤波器是无限脉冲响应滤波器, 不同。本文使用双线性变换法设计IIR ; 其次, 设计模拟滤波器( a (, 将模拟滤波器的H a (s , 从s 平面转换到z 平面, H ; 最后通过对H (z 的处理, 输出幅频特性等曲线图。 实例:Butterwort h 滤波器, 设计指标是: π, R p =1dB , w s =0. 3π, A s =15dB w p =0. 2 设计分析: Ωs 和Ωc 及阶数N 1 求模拟低通的技术指标(各截止频率Ωp 、 ΩP =(2/T 3tan (w p /2 , 采样周期T =1(原则上T 值可任取 , Ωs =2/T 3tan (w s /2 , Ωc =Ωs (100.13A s -1 -(1/2N , N =-(l g k s p /l g λs p , 其中k s p =((10R p /10-1 /(10A s/10-1 λ, s p =Ωs /Ωp 。2 设计Butterwort h 低通滤波器, 得到归一化函数H a (p , 去归一化后得到实际的传输函数H a (s 。3 用双线性变换法将H a (s 变换为数字低通滤波器的系统函数H (z : H (z =H a (s |-1 /s =2*(1-z -1(1+z 2几个重要的函数 1 butta p (函数。利用该函数, 在给定阶数N 后, 可运算出Butterwort h 模拟滤波器原型, 其传输函数为: H a (s =k/((s -p (1 (s -p (2 …(s -p (n 2 bilinear (函数。利用它实现双线性变换, 将S 域变换为Z 域, 其映射关系满足 H (z =H a (s |s =2*(1-z -1 /-1(1+z 3 sub p lot (。该函数用来分割图形窗口, 可在一个图形窗口中输出多幅图。其格式为sub p lot (m ,n , I , 表示图形窗口被分割成m 行n 列个子窗口, 本子窗口是第I 个。适合于有多个结果图的输出。 3MA TL AB 语言程序 1 主程序main. m 收稿日期:2001-06-11 作者简介:柳春峰, 男,1965年生, 大学本科, 讲师, 主要从事数字信号处理研究。 ?52?齐齐哈尔大学学报2001年%Imp ulse Invariance Transformation %Butterworth Low p ass Filter Desi g n w p =0. 23p i ; ws =0. 33p i ; R p =1; As =15; T =1; Fs =1/T ; Ome g a P =(2/T 3tan (w p /2 ; Ome g aS =(2/T 3tan (ws/2 ; e p =s q rt (10^(R p /10 -1 ; Ri pp le =s q rt (1/(1+e p 3e p ; Att n =1/(10^(As/20 ; [cs ,ds ]=af d-butt (Ome g a P ,Ome g aS ,R p ,As ; [b ,a ]=bilinear (cs ,ds , T ; [C ,B ,A ]=dir2cas (b ,a %p lot fi g ure (1 ; [db ,ma g , p ha , g rd ,w ]=f re q z-m (b ,a ; sub p lot (2,2,1 ; p lot (w/p i ,ma g ;title (′Ma g nit ude Rres p ;

您可能关注的文档

文档评论(0)

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

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

1亿VIP精品文档

相关文档