带阻数字滤波器设计-任务六_第1页
带阻数字滤波器设计-任务六_第2页
带阻数字滤波器设计-任务六_第3页
带阻数字滤波器设计-任务六_第4页
带阻数字滤波器设计-任务六_第5页
已阅读5页,还剩10页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

1、成 绩数字信号处理课程设计题目:任务6:设计带阻数字滤波器班级学号姓名任务分工工作量(%)MATLAB仿真,课程设计总结,汇总14.3数字滤波器设计14.3MATLAB仿真,汇总14.3第一题14.2第二题14.3第三题14.3第四题14.32014年 12 月目录1、简述四种滤波器的性能特点和幅度响应图21.1 低通滤波器21.2 高通滤波器21.3 带通滤波器21.4 带阻滤波器22、比较模拟巴特沃斯滤波器和切比雪夫滤波器的频域特性22.1 巴特沃斯滤波器22.2 切比雪夫滤波器23、简述巴特沃斯滤波器和切比雪夫模拟低通滤波器的设计方法23.1 巴特沃斯低通滤波器低通逼近23.1.1 逼近

2、原理23.1.2 巴特沃斯滤波器的设计23.2 切比雪夫滤波器23.2.1设计原理23.2.2 Chebyshev滤波器的设计24、简述模拟低通滤波器转换为数字滤波器的设计流程24.1 三种主要的方法24.2 设计流程25、任务6,设计带阻数字滤波器25.1说明所选择的设计方法的原理及设计流程25.2设计一个切贝雪夫I型带阻滤波器26、MATLAB滤波器仿真及结果分析26.1 画出滤波器的幅度响应曲线26.2 观察随机噪声通过滤波器效果26.3噪声信号f(n)通过设计的滤波器,观察时域波形和频域波形26.4 产生特定频率正弦波叠加检测数字带阻滤波器性能,滤波结果的分析27、课程设计总结28、参

3、考文献2附件21至3由于版权问题无法上传,但这部分都是数字信号书上抄的4、简述模拟低通滤波器转换为数字滤波器的设计流程4.1 三种主要的方法先模拟频带变换,再数字化:模拟低通、高通、带通、带阻模拟频带变换数字化数字低通、高通、带通、带阻模拟归一化原型c=1直接设计:数字低通、高通、带通、带阻模拟归一化原型c=1数字频带转换先数字化,在频带转换:数字化数字低通数字低通、高通、带通、带阻模拟归一化原型c=14.2 设计流程(1)将模拟低通滤波转换成其他相应模拟滤波器:模拟低通滤波器转换成模拟带阻滤波器的方法:首先定出低通滤波系统HLP(s),而低通的截止频率c与带阻滤波器阻带带宽B成反比,在HLP

4、(s)中代入低通到带阻的变换式,得到带阻滤波器的系统函数HBRp=HLPs,模拟低通滤波器转换成模拟低通滤波器的方法: 首先把数字滤波器的性能要求转换为与之对应的作为“样本”的模拟滤波器的性能要求,根据此性能要求设计模拟滤波器,既可以用查表得方法,也可以用解析的方法。然后通过冲激响应不变法或双线性变换法,将此“样本”模拟低通滤波器Ha(s)数字化为所需的低通滤波器H(z)。 模拟低通滤波器转换成模拟带通滤波器的方法:首先,定出低通系统函数HLP(s),其截止频率c应等于所需模拟带通滤波器的带宽c=B=2-1。在HLP(s)中代入低通到带通的变换式,得到带通滤波器的系统函数HBPp=HLPs,其

5、中s=p2+02p。模拟低通滤波器转换成模拟高通滤波器的方法:首先,确定低通系统函数HLP(s),其确定角频率c由高通参考角频率c选定。为了方便起见,一般都选c=1的归一化原型低通滤波器。在所得的HLP(s)中代入变换关系式,得到高通系统函数HHPp=HLPs,s=cpc。(2)有相应的模拟滤波器到数字滤波器的变换,使用双线性变换。(3)把模拟模拟的变换式与模拟数字的变换式结合起来,直接得到z平面的表达式,在经过化简就可得到相应的系统函数。以下我们将以模拟低通到数字带阻滤波器为例来说明具体步骤:由模拟低通到模拟带阻的变换这一模拟低通到带阻的变换关系为s=02pp2+02式中s为模拟低通原型拉普

6、拉斯变量(s=+j),p为模拟带阻的拉普拉斯变量(p=+j),0是模拟带阻滤波器的几何中心频率。令p=j,可得s=j0202-2故p平面的虚轴与s平面的虚轴相对应,代入s=j,消去j,可得=0202-2画出的0时的图形如图6-29所示,由此图可看出对应频率的关系为 =0 =0,= =c =1 =0 =-c =2其中c为低通滤波器的通带截止频率,故低通的阻带映像到带阻的阻带。将上面的映射结果代入式,可得c=02102-12-c=02202-22将相加、相减并化简,可得以下两个有用的关系式:0=12=2-1=02c=12c2,1分别为带阻滤波器上、下通带的截止频率,B为带阻滤波器的阻带带宽。从()

7、式看出,B与低通原型的带宽成反比。因此,有模拟低通滤波器确定模拟带阻滤波器的方法如下:(1) 定出低通系统函数HLP(s),而低通的截止频率c与带阻滤波器阻带带宽B成反比;(2) 在HLP(s)中的变换关系,可得到带阻滤波器系统函数HBRp=HLPs|s=02p/(p2+02)由模拟带阻到数字带阻的变换仍利用双线性变换p=c1-z-11+z-1把变换式结合起来,可得到直接从模拟低通原型滤波器的s平面变换成数字带阻滤波器的z平面的表达式s=c02pp2+02=02c1-z-11+z-1c2(1-z-1)2(1+z-1)2+02经推导后得s=c02(1-z-1)(c2-02)1-2c2-02c2+

8、02z-1+z-2下面来化简。利用双线性变换的频率间的关系=ctan2是数字带阻滤波器的数字频率。可得tan202=tan12tan22tan22-tan12=c2tan202c考虑到模拟带阻滤波器到数字带阻滤波器是带阻中心频率相对应的映射关系,设数字带阻滤波器的阻带中心频率为0,则有0=ctan02因而0c=0ccot02令D1=c02c2+02可得D1=ctan(2-12)令E1=2c2-02c2+02与带同滤波器一样,可以导出E1=2cos0=cos(2+12)cos(2-12)可得s=D1(1-z-1)1-E1z-1+z-2因此,从模拟低通系统函数HLP(s),进一步转换,就可以变成数

9、字带阻系统函数HBRp=HLPs|s=D1(1-z-1)1-E1z-1+z-2可以看出,数字带阻滤波器的极点数(阶数)等于模拟低通滤波器极点(阶数)的两倍。下面讨论模拟低通滤波器和数字带阻滤波器频率之间的关系。令s=j,z=ej,代入上式,经推导后可得=D1sincos-cos0其映射关系为=0 =0,= = = 0也就是说,确定低通滤波器的通带(=0 附近)映射到带阻滤波器的带阻范围之外(=0,附近及=附近),低通滤波器的带阻(=附近)映射到带阻滤波器上(= 0附近)。5、任务6,设计带阻数字滤波器5.1说明所选择的设计方法的原理及设计流程数字带阻滤波器的设计方法:(1) 是先将要设计的滤波

10、器的技术指标(主要是WP, W s),通过频率转变关系转换成模拟低通滤波器技术指标。(2) 依据这些技术指标设计出低通滤波器的转移函数.(3) 再依据频率转换关系变成所要设计的滤波器转移函数.以上只进行简要的阐述,具体设计看目录3.2.2 Chebyshev滤波器的设计5.2设计一个切贝雪夫I型带阻滤波器设计的切贝雪夫I型带阻滤波器满足,抽样频率fs=20kHz,技术指标为下通带截止频率为fp1=2kHz,上通带截止频率为fp2=7kHz,Rp=2dB,阻带下截止频率fst1=3kHz,阻带上截止频率为fst2=4kHz,As=20dB要求采用从低通滤波器直接变换到带阻滤波器的办法设计(1)

11、各对应的数字域临界频率(不必预畸变)为:wp1=2*fp1fs=2*200020000=0.2wp2=2*fp2fs=2*700020000=0.7wst1=2*fst1fs=2*300020000=0.3wst2=2*fst2fs=2*400020000=0.4(2) 求此带阻DF的中心频率w0cosw0=cos(wp2+wp12)cos(wp2-wp12)=cos(0.45)cos(0.25)=0.2212w0=1.3477515=0.4299(3) 求样本低通AF的p , st .利用预畸变换式可得:p=sinwp1coswp1-cosw0=sin(0.2)cos(0.2)-0.2212

12、=0.9999st1=sinwst1coswst1-cosw0=sin(0.3)cos(0.3)-0.2212=2.2069st2=sinwst2coswst2-cosw0=sin(0.4)cos(0.4)-0.2212=10.8300取st=minst1,st2=2.2069 在较小频率处满足衰减要求,则较大频率处衰减一定会更大,超过要求。(4) 求样本低通AF的阶次:因为:g=100.1As-1100.1Rp-1=13.0101 且s=stp=2.2071所以: Nlgg+g2-1lgs+s2-1=1.41470.6206=2.2794 所以取N=3(5) 查表7.5,N=3 ,Rp=2d

13、B 且 =100.1Rp-1=0.7648得归一化切比雪夫型样本低通AF的Han(s), Han(s)分子的系数为 12N-1=14 Hans=14s3+0.7378216s2+1.0221903s+0.3268901(6) 由于是切比雪夫型滤波器,故p=c ,可以直接利用p去归一化得到样本低通AF的Ha(s)Has=Hansp=p34s3+0.7378216ps2+1.0221903p2s+0.3268901p3=ds3+aps2+bp2s+cp3其中:,p=0.9999 ,=0.7648, d=p34=0.3268 , a=0.7378216 b=1.0221903,c=0.3268901

14、 (7)求所求的带阻DF的系统函数:s=1-z-21-2z-1cosw0+z-2=1-z-21-2ez-1+z-2 其中 e=cosw0=0.2212Hbz=Has|s=1-z-21-2ez-1+z-2=ds3+aps2+bp2s+cp3 |s=1-z-21-2ez-1+z-2=d(1-2ez-1+z-2)3(1-z-2)3+ap1-z-22(1-2ez-1+z-2)+bp2(1-z-2)(1-2ez-1+z-2)2+cp3(1-2ez-1+z-2)3=0.1058(1-1.3271z-1+ 3.587z-2 - 2.7409z-3+ 3.587z-4- 1.3271z-5+ z-6)1-0.

15、53924z-1- 0.43528z-2- 0.078721z-3+ 0.71683z-4 + 0.046705z-5- 0.3102z-6 6、MATLAB滤波器仿真及结果分析6.1 画出滤波器的幅度响应曲线图11 切比雪夫1型幅频带阻滤波器幅度响应6.2 观察随机噪声通过滤波器效果自行产生信号f(n)( 例如f(n)=square(1:256,50)+randn(1,256),这是一个长度为256点的信号,由占空比50%的方波叠加了高斯噪声而形成),画出f(n)的时域波形和频域波形。图12 原噪声的时域波形图13 噪声的频域波形6.3噪声信号f(n)通过设计的滤波器,观察时域波形和频域波形

16、 图14 滤波后的时域波形图15 滤波前后的频域波形滤波结果:可以看到在2000至7000Hz频带内赋值明显衰减。6.4 产生特定频率正弦波叠加检测数字带阻滤波器性能,滤波结果的分析(1) 将一个有1kHz,4kHz,8kHz正弦波进行叠加信号通过采样后通过滤波器观察滤波器性能。yt=sin2*1000*t+sin2*4000*t+sin(2*8000*t)对信号进行采样,采样频率fs=20kHz,采样点数N=256,MATLAB程序段为:f1=1000;%信号频率1kHzf2=4000;%信号频率4kHzf3=8000;%信号频率8kHzfs=20000;%采样频率20kHzN=256;%采

17、样点数t=(0:N-1)/fs;%采样时间ss=sin(2*pi*f1*t)+sin(2*pi*f2*t)+sin(2*pi*f3*t);%信号采样值(2) 将噪声信号通过滤波器后观察的时域与频域波形图16 滤波前后噪声信号时域波形图17 滤波前后的频域波形(3) 分析仿真结果通过滤波器后我们可以看到在4kH信号幅值衰减基本上到0,通过测试:频率/kHz原信号幅值滤波后信号衰减分贝陪数/dB1120105.61.11032125.798.792.09244120.50.242853.9148712697.982.18478125,9120.60.3735由所仿真测试数据得到滤波结果:1) 通带

18、内信号幅值衰减小于2dB(fst=1kHz,8kHz);2)阻带下限截止频率(fst=2kHz)处信号幅值衰减2.0924dB,阻带上限截止频率(fst=7kHz)处信号幅值衰减2.1847dB,基本满足设计要求;3)阻带上截止频率(fst=4kHz)处信号幅值衰减53.9148dB。较设计参数比较衰减倍数增加。滤波结果的分析:通过MATLAB仿真测试,所设计带阻滤波器满足设计要求,通带截止频率2dB,阻带截止频率20dB。7、课程设计总结此次带阻滤波器的课程设计,我们是先设计一个切比雪夫型低通模拟滤波器,然后通过双线性变换法转换等最后得到所需的数字带阻滤波器。后面有经过两种测试信号进行滤波测

19、试。在设计过程中我们小组成员有很多感悟。1、理论结合实际。由课本上知识得知,一个理想的滤波器是物理不可实现的,肯定会有一些误差,我们应尽量减小误差,理想逼近。此次课程设计中,我们一遍遍地改变设计参数,继而调试运行,查看调试出的图形结果,使它能尽量的逼近理想滤波器。 期间我们还发现书上的例题有错误,经过与老师的商量并获得老师的同意我们修改了所设计课程题目的参数。2、优秀的设计工具能大大减轻设计工作。这段时间内,我最大的收获就是对MATLAB有了更深刻的认识,体会到MATLAB强大的功能,它不仅具有高效的计算能力、灵活的图形处理能力、简单易懂的编程语言,更重要的是它对图形有超强的逼近模仿能力,应用

20、起来非常方便。 3、合作创造双赢。本次课程设计,我们体会到了大家一起相互交流、合作、齐心协力的重要性。不管是组内的交流还是组间的交流,我们遇到不懂的问题就进行相互讨论,最后都不会的就汇总起来一起问老师。同时,在老师的悉心指导和严格要求下,我们极大地提高了自己设计与分析的能力,同时也感受到理论与实践之间的差距,使我们对数字信号处理的应用有了更加深刻的认识。8、参考文献1程佩清.数字信号处理教程M.北京:清华大学出版社,2001.2刘树棠.数字信号处理使用MATLAB.西安:西安交通大学出版社,20023刘卫国.MATLAB程序设计与应用(第二版).北京:高等教育出版社,20064罗军辉等编著.M

21、ATLAB7.0在数字信号处理中的应用.北京:机械工业出版社,2005.5高西全.丁玉美.数字信号处理M西安电子科技大学出版社,2008.6桂志国.数字信号处理北京:科学出版社20097薛年喜.MATLAB在数字信号处理中的应用M.北京:清华大学出版社,2003.附件1、数字滤波器的幅度响应曲线代码clearclcb=0.1058 0.14040718 0.3795046 0.28998722 0.3795046 0.14040718 0.1058;a=1 0.53924 0.4352 0.078721 0.078721 0.046705 0.3102;H,W=freqz(b,a,20000);Hr=abs(H)fs=20000;figure(1);pl

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论