2022年课程方案双线性变换法IIR数字带通滤波器方案.docx

《2022年课程方案双线性变换法IIR数字带通滤波器方案.docx》由会员分享，可在线阅读，更多相关《2022年课程方案双线性变换法IIR数字带通滤波器方案.docx（30页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、精选学习资料 - - - - - - - - - 个人资料整理 仅限学习使用郑州轻工业学院课程设计说明书题目：姓名：院系）：专业班级：学号：指导老师：成果：时间：年月日至年月日名师归纳总结 - - - - - - -第 1 页,共 19 页精选学习资料 - - - - - - - - - 个人资料整理 仅限学习使用郑州轻工业学院课程设计任务书题目基于双线性变换法的IIR 数字高通滤波器设计专业、班级电子信息工程 08 级 2 班学号姓名主要内容、基本要求、主要参考资料等：主要内容 ：第一依据给定的性能指标,采纳双线性变换法设计IIR 数字高通滤波器；然后利用MATLAB 软件的 wavread

2、 函数读取 .wav 格式的语音信号,并利用所设计的滤波器对音频信号进行滤波处理,画出滤波前后信号的时域波形及频谱；最终回放语音信号,分析滤波前后的语音变化； 基本要求 ： 1、滤波器技术指标为：fp=3200Hz； Ap=1dB；fs=3000Hz；As=100dB 2、采纳双线性变换法设计 IIR 数字低通滤波器；3、把握利用 wavread函数读取 .wav 格式语音信号的方法；4、对语音信号进行滤波,并画出滤波前后信号的时域波形及频谱；5、回放语音信号,分析滤波前后的语音变化；主要参考资料：1、从玉良 .数字信号处理原理及其MATLAB 实现 M. 北京：电子工业出版社 .2022.7

3、 2、胡广书 .数字信号处理理论、算法与实现 M. 北京：清华高校出版社 .2003,8 完成期限： 2022.6.242022.6.28指导老师签名：课程负责人签名：2022 年 6 月 24日名师归纳总结 - - - - - - -第 2 页,共 19 页精选学习资料 - - - - - - - - - 个人资料整理 仅限学习使用摘要随着信息时代和数字世界的到来,数字信号处理已成为当今一门极其重要的学科和技术领域；在数字信号处理中起着重要的作用并已获得广泛应用的是数字滤波器 、稳固性好 仅运行在 0 与 l 两个电平状态 、敏捷性强等优点；时域离散系统的频域特性 :, 其中、分别是数字滤波

4、器的输出序列和输入序列的频域特性 数字滤波器和有限长冲激响应FIR 数字滤波器； IIR 数字滤波器的特点是,具有无限连续时间冲激响应,需要用递归模型来实现,其差分方程为： 1-1 系统函数为：,使其频率响应 Hz满意所期望得到的频域指标,即符合给定的通带截止频率、阻带截止频率、通带衰减系数和阻带衰减系数；名师归纳总结 - - - - - - -第 4 页,共 19 页精选学习资料 - - - - - - - - - 1.2 IIR数字滤波器设计原理个人资料整理仅限学习使用IIR 数字滤波器是一种离散时间系统,其系统函数为的FIR子系统的级联； IIR 数字滤波器的设计实际上是求解滤波器的系数

5、 和 ,它是数学上的一种靠近问题,即在规定意义上滤波器、切比雪夫 Chebyshev滤波器、椭圆 Cauer滤波器、贝塞尔 Bessel 滤波器等,这些典型的滤波器各有特点；这里介绍切比雪夫滤波器；切比雪夫滤波器特点：误差值在规定的频段上等波纹变化；巴特沃兹滤波器在通带内幅度特性是单调下降的,假如阶次肯定,就在靠近截止 处,幅度下降许多,或者说,为了使通带内的衰减足够小,需要的阶次N很高,为了克服这一缺点,采纳切比雪夫多项式来靠近所期望的；切比雪夫滤波器的在通带范畴内是等幅起伏的,所以在同样的通常内衰减要求下,其阶数较巴特沃兹滤波器要小；切比雪夫滤波器的振幅平方函数为 1-4 式中 c 为有效

6、通带截止频率,表示与通带波纹有关的参量,值越大通带不动愈大；VN切比雪夫滤波器的振幅平方特性如下列图：N为偶数, cos 2=1,得到 min,1-6）N为奇数, cos 2,得到 max,1-7）名师归纳总结 - - - - - - -第 6 页,共 19 页精选学习资料 - - - - - - - - - 个人资料整理 仅限学习使用图 1-2 切比雪夫滤波器的振幅平方特性有关参数的确定 : a、通带截止频率：预先给定；b、与通带波纹有关的参数 , 通带波纹表示成 1-10 1-11名师归纳总结 - - - - - - -第 7 页,共 19 页精选学习资料 - - - - - - - -

7、- 个人资料整理 仅限学习使用；1.2.2 双线性变换法为了克服冲激响应法可能产生的频率响应的混叠失真,这是由于从 S 平面到平面是多值的映射关系所造成的；为了克服这一缺点,可以采纳非线性频率压缩方法,将整个频率轴上的频率范畴压缩到 转换到 Z 平面上；也就是说,第一步先将整个- / T / T 之间,再用sT z=eS 平面压缩映射到S1 平面的 - / T / T 一条横带里；其次步再通过标准变换关系z=e s1 T 将此横带变换到整个Z 平面上去；这样就使 S 平面与 Z 平面建立了一一对应的单值关系,排除了多值变换性,也就排除了频谱混叠现象,映射关系如图 1-3j j 1 jImz/

8、Too11o1RezS平面/ TZ平面S1平面图 1-3 双线性变换的映射关系为了将 S 平面的整个虚轴j 压缩到 S1 平面 j 1 轴上的 - / T 到 / T 段上,可以通过以下的正切变换实现 2tan1T1-13）式中 , T 仍是采样间隔；T2当 1 由- / T 经过 0 变化到 / T 时, 由- 经过 0 变化到 +,也即映射了整个 j 轴；将式 Tej1 T/2e1T/2将此关系解读延拓到整个S 平面和 S1 平面,令 j =s,j 1=s1,就得名师归纳总结 - - - - - - -第 8 页,共 19 页精选学习资料 - - - - - - - - - s2e s 1

9、T/2es 1T/22tanh个人资料整理仅限学习使用s 1 T21es 1 T1-15 e s 1T/2s 1T/2s 1 TTeT2T1e再将 S1 平面通过以下标准变换关系映射到Z 平面z=e s1 T从而得到 S 平面和 Z 平面的单值映射关系为：s21z12s T1z1T 1s2 TT 2z1ss2T式1-10）与式 1-11）是 S 平面与 Z 平面之间的单值映射关系,这种变换 都是两个线性函数之比,因此称为双线性变换 式1-9）与式 1-10）的双线性变换符合映射变换应满意的两点要求；第一 , 把zejw可得tan2j1-18）s21ejj2T1ejT即 S 平面的虚轴映射到Z

10、平面的单位圆；j其次,将 s= +j 代入式 1-12）,得 2zT 2jT因此|z|2221-19）T222T由此看出,当 0 时,| z|0 时,| z|1 ；也就是说, S 平面的左 半平面映射到 Z 平面的单位圆内, S 平面的右半平面映射到 Z 平面的单位圆外,S 平面的虚轴映射到 Z 平面的单位圆上；因此,稳固的模拟滤波器经双线性变 换后所得的数字滤波器也肯定是稳固的；双线性变换法优缺点：双线性变换法与脉冲响应不变法相比,其主要的优名师归纳总结 - - - - - - -第 9 页,共 19 页精选学习资料 - - - - - - - - - 个人资料整理 仅限学习使用点是防止了频

11、率响应的混叠现象；这是由于S 平面与 Z 平面是单值的一一对应关系； S 平面整个 j 轴单值地对应于 Z 平面单位圆一周,即频率轴是单值变换关系；这个关系如式 1-12）所示,重写如下：2tan21-20）1-4 所示；T上式说明, S平面上 与 Z 平面的 成非线性的正切关系,如图由图 1-4 看出,在零频率邻近,模拟角频率 与数字频率 之间的变换关系接近于线性关系；但当 进一步增加时, 增长得越来越慢,最终当 时, 终止在折叠频率 = 处,因而双线性变换就不会显现由于高频部分超过 折叠频率而混淆到低频部分去的现象,从而排除了频率混叠现象；2tan2To图 1-4 双线性变换法的频率变换关

12、系但是双线性变换的这个特点是靠频率的严峻非线性关系而得到的,如式 1-12）及图 1-4 所示；由于这种频率之间的非线性变换关系,就产生了新的问题；第一,一个线性相位的模拟滤波器经双线性变换后得到非线性相位的数字滤波器,不再保持原有的线性相位了；其次,这种非线性关系要求模拟滤波器的幅 频响应必需是分段常数型的,即某一频率段的幅频响应近似等于某一常数 确定性能指标 在设计高通滤波器之前 , 第一依据工程实际的需要确定滤波器的技术指标 : 通带截止频率 fp=3200Hz；阻带截止频率 fs=3000Hz；阻带最小衰减 As=100dB和 通带最大衰减 Ap=1dB；2把频率转化为数字角频率 wp

13、=2*pi*fp*T ；ws=2*pi*fs*T ；3频率预畸变以及模拟高通性能指标转换成模拟低通性能指标用 =2/T*tanw/2 对高通数字滤波器Hz的数字边界频率预畸变 , 得到高通模拟滤波器Hs的边界频率主要是通带截止频率Wp；阻带截止频率Ws的转换；抽样频率 FS=8000Hz；通带截止频率 Wp2=2/T*tanwp/2 阻带截止频率 Ws2=2/T*tanws/2 名师归纳总结 - - - - - - -第 11 页,共 19 页精选学习资料 - - - - - - - - - 个人资料整理 仅限学习使用4模拟低通滤波器的构造借助切比雪夫 Chebyshev滤波器得到模拟低通滤波

14、器的传输函数Has；N,Wn=buttordwp2,ws2,Ap,As,s；z,p,k=buttapN；Bap,Aap=zp2tfz,p,k；G1=tfBap,Aap ；5模拟低通滤波器转换成模拟高通滤波器调用 lp2bp 函数将模拟低通滤波器转化为模拟高通滤波器；Bbs,Abs=lp2hpBap,Aap,Wn ；7模拟高通滤波器转换成数字高通滤波器利用双线性变换法将模拟高通滤波器；Has转换成数字高通滤波器Hz；Bbz,Abz=bilinearBbs,Abs,FS2.2 程序流程图第一确定性能指标,把频率转化为数字角频率,进而在进行频率预畸变,用 =2/T*tanw/2 对高通数字滤波器Hz

15、的数字边界频率预畸变 , 得到高通模拟滤波器 Hs的边界频率主要是通带截止频率Wp；阻带截止频率Ws的转换；抽样频率 fs=10KHz；上述预备工作做好之后,就先把模拟高通性能指标转换成模拟低通性能指标,然后设计模拟低通滤波器,借助切比雪夫Chebyshev滤波器得到模拟低通名师归纳总结 滤波器的传输函数Has；然后将模拟低通滤波器转化为模拟高通滤波器；最第 12 页,共 19 页后利用双线性变换法将模拟高通滤波器Has转换成数字高通滤波器Hz；- - - - - - -精选学习资料 - - - - - - - - - 个人资料整理 仅限学习使用开头读入数字滤波器技术指标将指标转换成归一化模拟

16、低通滤波器的指标设计归一化的模拟低通滤波器阶数N和 3db 截止频率模拟低通滤波器变换成模拟高通滤波器模拟高通滤波器转换成数 字高通滤波器终止 图 2-1 程序流程图 2.3 音频信号部分程序；N= length y；Y=ffty,N；yn=y+noise ；名师归纳总结 soundyn,FS,nbits；第 13 页,共 19 页HPASSyn=filterBbz,Abz,ynHPASSYN=fftHPASSyn；- - - - - - -精选学习资料 - - - - - - - - - 个人资料整理 仅限学习使用2. 仿真结果2.4.1 滤波器性能仿真源程序设计了模拟低通滤波器、模拟高通滤

17、波器与数字高通滤波器,对数 字高通滤波器的性能仿真如下：幅 频 特 性 曲 线 db0|Hz|-5000.511.522.53-100wrad5Hz000.10.20.30.40.50.60.70.80.91-5wrad图 2-2 滤波器性能仿真2.4.2 滤波器除噪性能为了实现滤波器的应用,程序中加入了有噪声的音频信号,通过对其滤波 处理,来显示数字高通滤波器的功能,下面显示未加入噪声,加入噪声和滤波 后的频谱波形图；名师归纳总结 - - - - - - -第 14 页,共 19 页精选学习资料 - - - - - - - - - 个人资料整理 仅限学习使用原 始 信 号 波 形 1 0.5

18、 0-0.5-100.511.522.53x 103.5原 始 信 号 频 谱51000050000 00.511.52x 104图 2-3 未加噪声的音频特性20滤 波 前 信 号 波 形x 10410滤 波 后 信 号 波 形x 10411230.512300-1-0.5-2-155滤 波 前 信 号 频 谱滤 波 后 信 号 频 谱80000x 1048000x 1046000600400040012312320002000055图 2-4 滤波前后的音频特性 由图可知,在 1K 到 2.4K 的频率范畴内,滤波器的滤波成效较好,把噪声 基本上全部滤除掉了；名师归纳总结 - - - -

19、- - -第 15 页,共 19 页精选学习资料 - - - - - - - - - 个人资料整理 仅限学习使用3 总结在拿到题目的那一刹,我被吓着了,赶忙去图书馆找材料,开头解剖题目,知道了什么是切比雪夫、什么是双线性变换,什么是巴特沃斯,学了那么久的 MATLAB从没看到这些个东西,在网上瞎搜了一阵,到是搜到了一些相关的源程 序,头都看大了,最终最终写出来了一个抱负的程序；这次课设下来,对设计带通数字滤波器的整个过程有了很好的把握,懂得 了设计滤波器的基本方法,对双线性变换法,切比雪夫滤波器有了肯定明白,同时呢也熟识了 MATLAB的环境,巩固了相关学问；最大的收成是初步明白了数字滤波器的

20、原理及设计方法,加深了对滤波器的熟识,一切从零开头,虽然没 有以肯定完善终止,但在这么短的时间内能够设计胜利已经出乎意料之外了,总之,收成仍是很大的；在做本次课程设计的过程中,我深深地感受到了自己所学到学问的有限,明 白了只学好课本上的学问是不够的,要通过图书馆和互联网等各种渠道来扩充 自己的学问；在试验过程中我们曾经遇到过问题；但是从中我们学习到了如何 对待遇到的困难,进一步培育了我们一丝不苟的科学态度和不厌其烦的耐心；全部的这些心得会对我以后的学习和工作有帮忙作用,忠心感谢学校给我们提 供这次试验机会；名师归纳总结 - - - - - - -第 16 页,共 19 页精选学习资料 - -

21、- - - - - - - 个人资料整理 仅限学习使用参考文献1 程佩青 . 数字信号处理 . 北京：清华高校出版社 ,2007.2 2 景振毅,张泽兵,董霖 . 北京：中国铁道出版社, 2022.12 3 从玉良 . 数字信号处理原理及其 社.2022.7 MATLAB实现 M. 北京：电子工业出版名师归纳总结 4 胡广书 . 数字信号处理理论、算法与实现M. 北京：清华高校出版社 .2003.8 第 17 页,共 19 页- - - - - - -精选学习资料 - - - - - - - - - 个人资料整理 仅限学习使用附录源程序clear %技术指标FS=8000；T=1/FS； %采

22、样频率大于 3200HZ的两倍fp=3200；fs=3000 ；Ap=1；As=20；wp=fp*2*pi*T ；ws=fs*2*pi*T ；wp2=2*tanwp/2/T ；ws2=2*tanws/2/T ；N,Wn=buttordwp2,ws2,Ap,As,s；z,p,k=buttapN；Bap,Aap=zp2tfz,p,k；G1=tfBap,Aap ；figure1；freqsBap,Aap ；title 模拟低通滤波器的频率响应 Bbs,Abs=lp2hpBap,Aap,Wn ；%模拟低通转化为模拟高通figure2；freqsBbs,Abs ；title 模拟高通滤波器的频率响应Bb

23、z,Abz=bilinearBbs,Abs,FSfigure3；freqzBbz,Abz,512,FS；title 数字滤波器的频率响应 figure4； ；%双线性变化把模拟高通转化为数字高通名师归纳总结 y,FS,nbits=wavread C:SHE.wav；第 18 页,共 19 页N= length y；Y=ffty,N；subplot2,1,1；ploty ；title原始信号波形 ；subplot2,1,2；plotabsY ；title原始信号频谱 ； noise=1+0.1*randnN,1；yn=y+noise ； soundyn,FS,nbits；HPASSyn=filterBbz,Abz,yn；HPASSYN=fftHPASSyn；figure5；subplot2,1,1；plotHPASSyn；title滤波后信号波形 ；- - - - - - -精选学习资料 - - - - - - - - - 名师归纳总结 subplot2,1,2；个人资料整理仅限学习使用第 19 页,共 19 页plotabsHPASSYN；title滤波后信号频谱 ；- - - - - - -