巴特沃斯高通数字滤波器设计_第1页
巴特沃斯高通数字滤波器设计_第2页
巴特沃斯高通数字滤波器设计_第3页
巴特沃斯高通数字滤波器设计_第4页
巴特沃斯高通数字滤波器设计_第5页
已阅读5页,还剩5页未读 继续免费阅读

下载本文档

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

文档简介

1、巴特沃斯高通数字滤波器设计要求:3dB数字截止频率为,阻带下边频rad,阻带衰减为。一、课程设计目的:数字信号处理(Digital Signal Processing DSP)是20世纪60年代以来,随着信息学科和计算机学科的高速发展而迅速发展起来的一门新兴学科。数字信号处理是利用计算机或其他专用处理设备,以数值计算的方法对信号进行采集、变换、滤波、压缩、传输、估值与识别等加工处理,借以达到提取信息和便于应用的目的的一种技术。数字信号处理随着计算机技术信息技术的进步获得了飞速的发展。数字信号处理已广泛应用于科学研究和工程技术的各个领域,是新一代IT工程师必须掌握的信息处理技术。它在越来越多的应

2、用领域中迅速替代传统的模拟信号处理技术,并且开辟出许多新的领域。数字信号处理有很多深奥的数学概念,理论也相对抽象,而且是一门理论与实践密切结合的课程。我们通过课程设计深入掌握课程内容,深入理解与消化关于巴特沃斯滤波器的基本理论,锻炼我们独立解决问题的能力,培养我们的创新意识,加强我们的实践学习。二、设计原理:1、数字滤波器所谓数字滤波器,是指输入输出均为数字信号,通过数字运算处理改变输入信号所含频率成分的相对比例,或者滤除某些频率成分的数字器件或程序。正因为数字滤波通过数值运算实现滤波,所以数字滤波器处理精度高、稳定、体积小、重量轻、灵活、不存在阻抗匹配问题,可以实现模拟滤波器无法实现的特殊滤

3、波功能。按照不同的分类方式,数字滤波器可以有很多种类型,但总起来可以分为两大类:经典滤波器和现代滤波器。经典滤波器的特点是其输入信号中有用的频率成分和希望滤除的成分分别占有不同的频带,通过一个合适的选频滤波器滤除干扰,得到纯净信号,达到滤波目的。但是,如果信号和干扰的频谱相互重叠,则经典滤波器无法有效滤除干扰,最大限度恢复信号,这就需要现代滤波器。现代滤波器是根据随机信号的一些统计特性,在某种最佳准则下,最大限度抑制干扰,同时最大限度恢复信号,达到最佳的滤波效果的目的。2、模拟滤波器模拟滤波器的理论和设计方法已经发展的相当成熟,且有多种典型的模拟滤波器供我们选择,如巴特沃斯滤波器、切比雪夫滤波

4、器、椭圆滤波器、贝塞尔滤波器等。这些滤波器有严格的设计公式、现成的曲线和图表供设计人员使用,而且所涉及的系统函数都满足电路实现条件。这些典型的滤波器各有特点:巴特沃斯滤波器具有单调下降的幅频特性;切比雪夫滤波器的幅频特性在通带或阻带有等波纹特性,可以提高选择性;贝塞尔滤波器通带内有较好的线性相位特性;椭圆滤波器的特性相对前三者是最好的,但通带和阻带内均呈现等波纹特性,相位特性的非线性也稍严重。设计时,根据具体要求选择滤波器的类型。现在,我们分别使用、H、H表示模拟滤波器的单位冲激响应、频率响应函数、系统函数,三者的关系为:可以用、中任意一个描述模拟滤波器,也可以用线性常系数微分方程描述模拟滤波

5、器。但在设计滤波器时,设计指标一般由幅频响应函数给出,而模拟滤波器设计就是根据设计指标,求系统函数。工程上通常用所谓的损耗函数(也称为衰减函数)来描述滤波器的幅频响应特性,对归一化幅频响应函数,定义如下(其单位是分贝,用dB表示): dB应当注意,损耗函数和幅频特性函数只是滤波器幅频响应特性的两种描述方法。损耗函数的优点是对幅频响应的取值非线性压缩,放大了小的幅度,从而可以同时观察通带和阻带的频响特性的变化情况。模拟低通滤波器的设计指标参数有、和。其中和分别称为通带边界频率和阻带截止频率,称为通带最大衰减(即通带中允许的最大值),称为阻带最小衰减(即阻带上允许的最小值),和的单位是dB。对于单

6、调下降的幅度特性,和可表示为:滤波器的技术指标给定后,需要设计一个系统函数,希望其幅度平方函数满足给定的指标。一般的滤波器单位冲激响应为实函数,因此如果能由、和求出,那么就可以求出,由此可求出所需要的。必须是因果稳定的,因此极点必须落在s平面的左半平面,相应的的极点必然落在右半平面。这就是由求所需要的的具体原则。3、巴特沃斯滤波器巴特沃斯滤波器是电子滤波器的一种。巴特沃斯滤波器的特点是通频带的频率响应曲线最平滑。这种滤波器最先由英国工程师斯替芬·巴特沃斯(Stephen Butterworth)在1930年发表在英国无线电工程期刊的一篇论文中提出的。巴特沃斯滤波器的特点是通频带内的频

7、率响应曲线最大限度平坦,没有起伏,而在阻频带则逐渐下降为零。 在振幅的对数对角频率的波特图上,从某一边界角频率开始,振幅随着角频率的增加而逐步减少,趋向负无穷大。 一阶巴特沃斯滤波器的衰减率为每倍频6分贝,每十倍频20分贝。二阶巴特沃斯滤波器的衰减率为每倍频12分贝、 三阶巴特沃斯滤波器的衰减率为每倍频18分贝、如此类推。巴特沃斯滤波器的振幅对角频率单调下降,并且也是唯一的无论阶数,振幅对角频率曲线都保持同样的形状的滤波器。只不过滤波器阶数越高,在阻频带振幅衰减速度越快。其他滤波器高阶的振幅对角频率图和低级数的振幅对角频率有不同的形状。4、设计IIR数字滤波器的频率变换法(1)、归一化模拟低通

8、原型到数字高通滤波器的变换设归一化模拟低通原型滤波器的系统函数是,p为模拟域内的拉普拉斯变量。在模拟域内从低通到高通变换,即以代替p,有进行反归一化后,即以代入上式,得到反归一化后的高通滤波器的传输函数;再进双线性变换,得数字高通滤波器的系统函数为;直接由归一化低通原型变换为数字高通滤波器时,变换关系为,其中,为模拟滤波器3dB截止频率;为数字滤波器3dB截止频率。(2)归一化模拟低通原型到数字带通滤波器的频率变换直接寻求从模拟低通到数字带通之间的映射关系:,其中,为模拟滤波器3dB截止频率;分别是数字带通滤波器的中心频率、通带上边频、通带下边频和带通的上阻带下边频。(3)归一化模拟低通原型到

9、数字带阻滤波器的变换直接寻求从模拟低通到数字带阻之间的映射关系:其中,为模拟滤波器3dB截止频率;分别为数字带阻滤波器的中心频率、阻带上边频、阻带下边频和带阻的上通带下边频。(4)数字低通滤波器到数字低通的变换由截频为的数字低通滤波器系统函数得到截频为的数字低通滤波器系统函数,变换关系为(5)数字低通滤波器到数字高通变换由截频的数字低通滤波器系统函数得到截频为的数字高通滤波器系统函数,变换关系为(6)数字低通滤波器到数字带通变换由截频为的数字低通滤波器系统函数得到中心频率为,带通上、下截频分别为的数字带通滤波器系统函数,变换关系为5、MATLAB模拟频率变换函数MATLAB信号处理工具箱函数l

10、p2lp、lp2bp、lp2hp、lp2bs是模拟频率变换函数。(1)低通模拟滤波器到低通模拟滤波器变换函数lp2lp调用格式说明如下。bs,as=lp2lp(b,a,Wl)该函数将截止频率为1rad/s(归一化截止频率)的模拟低通原型滤波器变换成截止频率为Wl的低通滤波器。(2)低通模拟滤波器到带通模拟滤波器的变换函数lp2bp调用格式如下。bs,as=lp2bp(b,a,W0,Bw)该函数将截止频率为1rad/s(归一化截止频率)的模拟低通原型滤波器变换成具有指定带宽和中心频率为的带通滤波器。如果要求滤波器的低端截止频率为,高端截止频率为,则按下式可计算带宽和中心频率为的带通滤波器。Bw=

11、W2-W1W0=sqrt(W1*W2)(3)低通模拟滤波器到高通模拟滤波器变换函数lp2hp调用格式如下。bs,as=lp2hp(b,a,Wh)该函数将截止频率为1rad/s(归一化截止频率)的模拟低通原型滤波器变换成截止频率为的高通滤波器。(4)低通模拟滤波器到带阻模拟滤波器变换函数lp2bs调用格式如下。bs,as=lp2bs(b,a,Wh)该函数将截止频率为1rad/s(归一化截止频率)的模拟低通原型滤波器变换成具有指定带宽和中心频率为的带阻滤波器。三、设计步骤与流程框图:求出模拟低通滤波器指标设计模拟低通滤波器将模拟滤波器转换为数字滤波器利用MATLAB编程实现滤波器得到结果开始四、设

12、计过程:1、设计要求:设计一巴特沃斯高通数字滤波器,3dB数字截止频率为,阻带下边频rad,阻带衰减为。2、设计过程:(1) 求模拟低通滤波器指标。(2) 设计模拟低通滤波器,取N=4,查表得(3) 将模拟滤波器转换为数字滤波器(4) 利用MATLAB设计巴特沃斯高通滤波器源程序:频率变换函数编程:wp=0.2*pi;ws=0.05*pi;Ap=3;As=48;Wp=tan(wp/2);Ws=tan(ws/2);N,wn=buttord(Wp,Ws,Ap,As,'s');fprintf('滤波器阶数N=%.0fn',N)b,a=butter(N,1,'s

13、');numa,dena=lp2hp(b,a,Wp);numd,dend=bilinear(numa,dena,0.5);disp('分子系数b');fprintf('%.4e ',numd);fprintf('n');disp('分母系数a');fprintf('%.4e ',dend);fprintf('n');w=linspace(0,pi,1024);h=freqz(numd,dend,w);plot(w/pi,20*log10(abs(h);axis(0 1 -50 0);grid

14、;xlabel('归一化频率');ylabel('幅度/dB');程序运行结果如下:滤波器阶数N=4分子系数b9.9999e-001 -4.0000e+000 6.0000e+000 -4.0000e+000 9.9999e-001分母系数a2.3103e+000 -5.4742e+000 5.3459e+000 -2.4366e+000 4.3290e-001其幅度响应函数如图所示:五、结果分析:由最终所得到的运行结果以及该滤波器的幅度响应函数图像可以看出,我们所设计的滤波器完全符合最初的要求,是一个巴特沃斯高通滤波器,3dB数字截止频率为,阻带下边频 rad,阻带衰减为。六、总结:通过这次课程设计,使我对数字信号处理又有了更深一层的了解。巴特沃斯滤波器看似简单,但是在

温馨提示

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

评论

0/150

提交评论