资源描述：
《基于matlab的数字带通滤波器》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、西安文理学院机械电子工程系课程设计报告专业班级题目基于MATLAB的数字带通滤波器学号学生姓名指导教师2011年12月西安文理学院机械电子工程系11课程设计任务书学生姓名专业班级学号指导教师职称副教授教研室电子信息工程课程数字信号处理题目基于MATLAB的数字带通滤波器设计任务与要求设计任务：要求设计一个IIR带通滤波器，其中通带的中心频率为，通带的截止频率，，通带最大衰减；阻带最小衰减，阻带截止频率，。设计要求：1.根据设计任务要求给出实现方案及实现过程。2.给出所实现的滤波器幅频特性及相频特性曲线并加以分析
2、。3.论文要求思路清晰，结构合理，语言流畅，书写格式符合要求。开始日期2011.12.19完成日期2011.12.302011年12月18日11一、设计任务设计一数字带通滤波器，用IIR来实现，其主要技术指标：通带边缘频率：wp1=0.4，wp2=0.6通带最大衰减：Ap=3dB阻带边缘频率：ws1=0.3，ws2=0.7阻带最小衰减：As=15dB设计总体要求：用MATLAB语言编程进行设计，给出IIR数字滤波器的参数，给出幅度和相位响应曲线，对IIR实现形式和特点等方面进行讨论。二、设计方法IIR数字滤波器
3、具有无限宽的冲激响应，与模拟滤波器相匹配，所以IIR滤波器的设计可以采取在模拟滤波器设计的基础上进一步变换的方法。比较常用的原型滤波器有巴特沃什滤波器（Butterworth）、切比雪夫滤波器（Chebyshev）、椭圆滤波器（Ellipse）和贝塞尔滤波器（Bessel）等。他们有各自的特点，巴特沃什滤波器具有单调下降的幅频特性；切比雪夫滤波器的幅频特性在通带和阻带里有波动，可以提高选择性；贝塞尔滤波器通带内有较好的线性相位特性；椭圆滤波器的选择性最好。本设计IIR数字滤波器采用巴特沃什滤波器[3]。设计巴特
4、沃什数字滤波器时，首先应根据参数要求设计出相应的模拟滤波器，其步骤如下：（1）由模拟滤波器的设计指标wp，ws，Ap，As和式（1）确定滤波器阶数N。（1）（2）由式（2）确定wc。11（2）（3）若N是奇数，则有式（3）和式（4）确定滤波器的系统函数；若N是偶数，则式（3）和式（5）确定滤波器的系统函数。然后把设计好的模拟Butterworth滤波器变换成数字滤波器，常采用的方法有脉冲响应不变法和双线性变换法[4]。（3）（4）（5）巴特沃什滤波器是根据幅频特性在通频内具有最平坦特性而定义的滤波器。一维巴特沃
5、什滤波器的平方幅频特性函数为（6）可以看出，滤波器的幅频特性随着滤波器阶次N的增加而变得越来越好。在截止频率Ωc处的函数值始终为的情况下，在通带内更多的频带区的值接近1；在阻带内函数更迅速的趋近0。巴特沃什滤波器的主要特征：（1）对于所有N，（2）对于所有N，（3）是Ω的单调下降函数。（4）随着阶次N的增大而更加接近于理想滤波器。11最常用于设计IIR数字带通滤波器的方法基于原型模拟传输函数的双线性变换。模拟传输函数通常是：巴特沃什、切比雪夫1型、切比雪夫2型和椭圆传输函数。巴特沃什低通传输函数在dc处有一个最
6、大平坦的幅度响应，且随着频率增加幅度响应单调地减少。本设计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）的分母
7、多项式和分子多项式。三、数字带通滤波器的具体设计及实现3.1用MATLAB设计的IIR数字带通滤波器的具体程序如下：ws=[0.30.7];%ws为归一化阻带边界频率wp=[0.40.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函数%单位圆上将函数转化
8、为256等份w=0:pi/255:pi;%计算相位h=freqz(num,den,w);%频率计算g=20*log10(abs(h));%计算幅度pha=angle(h);%相角计算%画带通滤波器的幅度曲线图11subplot(1,2,1);%创建子图plot(w/pi,g);%二维直角坐标曲线图grid%画坐标网格图axis([01-605]);%轴的刻度与表现xlabel('om
