基于MATLAB的数字带通滤波器课程设计报告.doc

西安文理学院机械电子工程系课程设计报告专业班级 08级电子信息工程1班 题 目 基于MATLAB的数字带通滤波器 学 号 学生姓名 指导教师 2011 年 12 月西安文理学院机械电子工程系课程设计任务书学生姓名 _专业班级 _ 学 号 _ 指导教师 _ 职 称 副教授 教研室 电子信息工程 课 程 数 字 信 号 处 理 题目 基于MATLAB的数字带通滤波器设计 任务与要求设计任务：要求设计一个IIR带通滤波器，其中通带的中心频率为，通带的截止频率，通带最大衰减；阻带最小衰减，阻带截止频率，。设计要求： 1. 根据设计任务要求给出实现方案及实现过程。2. 给出所实现的滤波器幅频特性及相频特性曲线并加以分析。3. 论文要求思路清晰，结构合理，语言流畅，书写格式符合要求。开始日期 2011.12.19 完成日期 2011.12.30 2011年 12月 18 日一、设计任务设计一数字带通滤波器，用IIR来实现，其主要技术指标：通带边缘频率：wp1=0.4，wp2=0.6通带最大衰减：Ap=3dB阻带边缘频率：ws1=0.3，ws2=0.7阻带最小衰减：As=15dB设计总体要求：用MATLAB语言编程进行设计，给出IIR数字滤波器的参数，给出幅度和相位响应曲线，对IIR实现形式和特点等方面进行讨论。二、设计方法IIR数字滤波器具有无限宽的冲激响应，与模拟滤波器相匹配，所以IIR滤波器的设计可以采取在模拟滤波器设计的基础上进一步变换的方法。比较常用的原型滤波器有巴特沃什滤波器（Butterworth）、切比雪夫滤波器（Chebyshev）、椭圆滤波器（Ellipse）和贝塞尔滤波器（Bessel）等。他们有各自的特点，巴特沃什滤波器具有单调下降的幅频特性；切比雪夫滤波器的幅频特性在通带和阻带里有波动，可以提高选择性；贝塞尔滤波器通带内有较好的线性相位特性；椭圆滤波器的选择性最好。本设计IIR数字滤波器采用巴特沃什滤波器3。设计巴特沃什数字滤波器时，首先应根据参数要求设计出相应的模拟滤波器，其步骤如下：（1）由模拟滤波器的设计指标wp，ws，Ap，As和式（1）确定滤波器阶数N。 （1）（2）由式（2）确定wc。 （2）（3）若N是奇数，则有式（3）和式（4）确定滤波器的系统函数；若N是偶数，则式（3）和式（5）确定滤波器的系统函数。然后把设计好的模拟Butterworth滤波器变换成数字滤波器，常采用的方法有脉冲响应不变法和双线性变换法4。 （3） （4） （5）巴特沃什滤波器是根据幅频特性在通频内具有最平坦特性而定义的滤波器。一维巴特沃什滤波器的平方幅频特性函数为 （6）可以看出，滤波器的幅频特性随着滤波器阶次N的增加而变得越来越好。在截止频率c处的函数值始终为的情况下，在通带内更多的频带区的值接近1；在阻带内函数更迅速的趋近0。巴特沃什滤波器的主要特征：（1）对于所有N，（2）对于所有N，（3）是的单调下降函数。（4）随着阶次N的增大而更加接近于理想滤波器。最常用于设计IIR数字带通滤波器的方法基于原型模拟传输函数的双线性变换。模拟传输函数通常是：巴特沃什、切比雪夫1型、切比雪夫2型和椭圆传输函数。巴特沃什低通传输函数在dc处有一个最大平坦的幅度响应，且随着频率增加幅度响应单调地减少。本设计IIR数字带通滤波器就采用巴特沃什型。巴特沃什型数字带通滤波器的设计，MATLAB调用形式如下：N,wn=buttord(wp,ws,Ap,As)num,den=butter(N,wn)在buttord中，调用参数wp，ws，是数字带通滤波器的归一化通带和阻带截止频率，Ap是通带最大衰减，As是阻带最小衰减。返回参数N是巴特沃什滤波器的阶数，wn是3dB截止频率，函数butter获得数字（den）和（num）是巴特沃什滤波器系统函数H（z）的分母多项式和分子多项式。三、数字带通滤波器的具体设计及实现3.1用MATLAB设计的IIR数字带通滤波器的具体程序如下：ws=0.3 0.7; %ws为归一化阻带边界频率wp=0.4 0.6; %wp为归一化通带边界频率Ap=3; %Ap单位为dB的通带波形As=15; %As单位为dB的最小阻带衰减n1,wn1=buttord(wp,ws,Ap,As);%估计数字滤波器的阶数 %wn1为归一化截止频率num,den=butter(n1,wn1);%设计数字滤波器%计算单位为dB的gain函数%单位圆上将函数转化为256等份w=0:pi/255:pi; %计算相位h=freqz(num,den,w); %频率计算g=20*log10(abs(h); %计算幅度pha=angle(h); %相角计算%画带通滤波器的幅度曲线图subplot(1,2,1); %创建子图plot(w/pi,g); %二维直角坐标曲线图grid %画坐标网格图axis(0 1 -60 5); %轴的刻度与表现xlabel(omega/pi);ylabel(幅度);%横坐标表示频率%纵坐标表示分贝数title(数字带通滤波器幅频曲线); %写标题subplot(1,2,2);plot(w/pi,pha);grid %画带通滤波器的相位曲线图xlabel(omega/pi);ylabel(相位);title(数字带通滤波器相位曲线);3.2IIR数字带通滤波器的仿真结果及其分析从仿真过程和结果可以得出：IIR数字滤波器的阶数为：n1=3截止频率为：wn1 = 0.3759 0.6241系统函数的分子系数为：num = 0.0311 0 -0.0933 0 0.0933 0 -0.0311系统函数的分母系数为：den = 1.0000 -0.0000 1.4703 -0.0000 0.9195 -0.0000 0.2004 根据运行结果可写出IIR滤波器的系统函数为：根据上述系统传递函数得到仿真数字滤波器的频率特性曲线如下图1 IIR数字带通滤波器的幅频响应和相位响应从图可以看出巴特沃什滤波器具有单调下降的幅频特性，并在截止频率wc处的函数值为0.5的情况下，通带内的值为1；在阻带内函数值更迅速的趋近0。相频特性曲线则呈非线性这表明会使信号产生失真。3.3 其它方法的滤波器设计3.3.1 Chebyshev型滤波器设计条件同上，程序如下：Wp=0.4,0.6;Ws=0.0154,0.7;ap=3;as=15; N,Wn=cheb2ord(Wp,Ws,ap,as); num,den=cheby2(N,as,Wn);figure;freqz(num,den)特性曲线如图所示：图2 Chebyshev型滤波器的幅频特性曲线图3 Chebyshev型滤波器的相频特性曲线由图可知，Chebyshev II型滤波器在0.0154到0.6左右为单调递减趋势的通带，在接近0.6左右时即为阻带Chebyshev II型滤波器其通带内为单调，阻带内等波纹，但其阻带内波纹较大。3.3.2椭圆型滤波器设计条件同上，程序如下：Wp=0.4,0.6;Ws=0.0154,0.7;ap=3;as=15;N,Wn=ellipord(Wp,Ws,ap,as);num,den=ellip(N,ap,as,Wn);figure;freqz(num,den)特性曲线如图所示：图4 椭圆型滤波器的幅频特性曲线图5 椭圆型滤波器的相频特性曲线椭圆滤波器的幅频特性曲线与Butterworth和chebyshev滤波器相比，通带附近波形斜率很大，ellip函数可以得到下降斜度更大的滤波器，得到通带和阻带均为等波纹。一般情况下，椭圆滤波器能以最低的阶实现指定的性能指标。其相频特性曲线为三者中波动最大的，其在通带范围时迅速升高和下降，然后恢复正常。四、结论利用MATLAB的信号处理工具箱强大的信号处理功能，采用编程的方法设计数字滤波器，使设计达到了最优化；并且它可以快速的实现数字滤波器的仿真，使设计达到了最简化。本文设计的数字滤波器采用IIR来实现， IIR采用巴特沃什函数法，按设计指标要求进行设计。利用MATLAB采用编程的方法实现。通过对IIR数字滤波器程序的编写和实现可知IIR数字滤波器幅频特性较好，而相频特性曲线呈非线性，会使信号产生失真，阶数比较低，系统相对不稳定，这个问题有待解决。五、心得体会经过两个星期的课程设计，平时看课本时，有时问题老是弄不懂，做完课程设计，那些问题就迎刃而解了。而且还可以记住很多东西。认识来源于实践，实践是认识的动力和最终目的，实践是检验真理的唯一标准。所以这课程设计对我们的作用是非常大的。 在制作基于MATLAB数字滤波器时，发现细心耐心，恒心一定要有才能做好事情，首先是要正确安装MATLAB软件，学习该软件是基本使用方法，学习MATLAB编程语言，兼顾到方方面面去考虑是很需要的，否则只是一纸空话，在设计过程中，经常会遇到这样那样的情况，就是心里想老想着这样的接法可以行得通，但实际写出来的程序，总是出现异常错误，因此耗费在这上面的时间用去进行调试。 我认为我们的工作是一个团队的工作，团队需要个人，个人也离不开团队，必须发扬团结协作的精神。某个人的离群都可能导致整项工作的失败。设计中只有一个人知道原理是远远不够的，必须让每个人都知道，否则一个人的错误，就有可能导致整个工作失败。团结协作是我们实习成功的一项非常重要的保证。而这次实习也正好锻炼我们这一点，这也是非常宝贵的。这次课程设计,虽然短暂但是让我得到多方面的提高：1、提高了我们的逻辑思维能力，使我们在MATLAB的使用和编程语言分析与设计上有了很大的进步。加深了我们对MATLAB编程语言的认识，进一步增进了对逻辑事务的处理能力。另外，我们还更加充分的认识到，MATLAB教程这门课程在科学发展中的至关重要性2，查阅参考书的独立思考的能力以及培养非常重要，我们在编写程序时，遇到很多不理解的东西，有的我们通过查阅参考书弄明白，有的通过网络查到，但由于时间和资料有限我们更多的还是独立思考。通过这次课程设计，加强了我们动手、思考和解决问题的能力，我们