MATLAB数字滤波器课程设计

1、课程设计任务书20102011学年第一学期专业： 通信工程 学号： 0* 姓名： * 课程设计名称： 数字信号处理课程设计 设计题目： 基于MATLAB的数字带通滤波器的设计 完成期限：自 2011 年 1 月 3 日至 2011 年 1 月 9 日共 1 周一、设计目的1、巩固所学的理论知识；2、提高综合运用所学理论知识独立分析和解决问题的能力；3、更好地将理论与实践相结合；4、掌握信号分析与处理的基本方法与实现；5、熟练使用MATLAB语言进行编程实现。二、设计内容设计巴特沃斯数字带通滤波器，要求通带范围为0.25rad0.45rad，通带最大衰减为3dB，阻带范围为00.15rad和0.

2、55radrad，阻带最小衰减为40dB。三、设计要求调用buttord和butter设计，并显示数字滤波器系统函数H(Z)的系数，绘制数字滤波器的损耗函数和相频特性函数。并分析这种设计对应于脉冲响应不变法还是双线性变换法？四、设计条件计算机、MATLAB语言环境五、参考资料1数字信号处理（第三版），丁玉美，高西全.西安电子科技大学出版社，2000.2MATLAB及在电子信息课程中的应用，陈怀堔，吴大正，高西全.电子工业出版社，2006.3MATLAB 7.0从入门到精通，求是科技.人民邮电出版社，2006.4数字信号处理（第三版）学习指导，高西全，丁玉美.西安科技大学出版社，2001.指导教

3、师（签字）： 教研室主任（签字）： 批准日期： 年 月 日摘 要数字滤波器是指输入、输出均为数字信号，通过数值运算处理改变输入信号所含频率成分的相对比例，或者滤除某些频率成分的数字器件和程序。经典数字滤波器从滤波特性上分类，可以分成低通、高通、带通和带阻等滤波器。数字滤波器的幅频特性表示信号通过该滤波器后各频率成分振幅衰减情况，而相频特性反映各频率成分通过滤波器后在时间上的延时情况。本次课程设计先是对巴特沃斯数字带通滤波器有关理论知识作介绍，在性能指标分析基础上运用MATLAB的buttord和butter函数设计程序，得到损耗函数和相频特性函数图像，写出数字滤波器系统函数H(Z)，并对结果进

4、行分析，最后总结课程设计体会。关键词：数字带通滤波器；MATLAB；系统函数；损耗函数；相频特性目 录1 课题描述12设计原理12.1数字滤波器设计方法概述12.2模拟滤波器频率变换模拟带通滤波器的设计22.2.1带通滤波器的指标要求22.2.2模拟低通滤波器变换成数字带通滤波器33设计内容43.1设计思路43.2设计步骤43.3 MATLAB程序63.4结果分析73.4.1损耗函数与相频特性73.4.2双线性变换8总 结11参考文献12课程设计说明书（论文）1 课题描述在社会飞速发展的今天，数字信号处理已成为当今一门极其重要的学科和手段，数字信号处理以其在处理过程中的实时、快速、稳定，故而在

5、雷达、通信、音视频处理、自动控制、航空航天、医疗和家用电器等技术领域得到了广泛的应用。数字信号处理学科的一项重大发展是关于数字滤波器设计方向的研究，数字滤波是数字信号处理的重要内容之一。与模拟滤波器相比，其具有精度、稳定度、灵活性高，不需要阻抗匹配和便于大规模集成等优点。60年代中期，数字滤波器形成一套完整的正规理论，设计方法趋于成熟，如巴特沃斯滤波器、切比雪夫滤波器、椭圆滤波器，贝塞尔滤波器等。数字滤波器从结构上可分为无限冲击响应数字滤波器和有限冲击响应数字滤波器；从功能上分为低通、带通、高通和带阻滤波器。巴特沃斯滤波器具有单调下降的幅频特性，在阻频带逐渐下降为零，其通频带的频率响应曲线最大

6、限度平坦，没有起伏。本文针对的就是巴特沃斯数字带通滤波器的设计。MATLAB是用于科学计算、可视化的高性能软件环境。它集数值分析、矩阵运算、信号处理和图形显示于一体，具有界面友好的用户环境。它的信号处理工具箱包含了各种数字信号处理技术，是一个非常优秀的算法研究与辅助设计的工具。通过调用MATLAB语言中的buttord和butter函数可用于设计数字带通滤波器。本文即是利用MATLAB工具来设计巴特沃斯数字带通滤波器的。2设计原理 2.1数字滤波器设计方法概述IIR滤波器和FIR滤波器的设计方法完全不同。IIR滤波器设计方法有间接法和直接法，间接法是借助于模拟滤波器的设计方法进行的。其设计步骤

7、是：先设计过渡模拟滤波器得到系统函数，然后将按某种方法转换成数字滤波器的系统函数。这是因为模拟滤波器的设计方法已经很成熟，不仅有完整的设计公式，还有完善的图标和曲线供查阅；另外，还有一些典型的优良滤波器可供我们使用。直接法直接在频域或者时域中设计数字滤波器，由于要解联立方程，设计时需要计算机辅助设计。FIR滤波器不能采用间接法，常用的设计方法有窗函数法、频率采样法和切比雪夫等波纹逼近法。模拟低通滤波器的设计是设计其他滤波器的基础。模拟高通、带通和带阻滤波器的设计过程是：先将希望设计的各种滤波器的技术指标转换为低通滤波器技术指标，然后设计响应的低通滤波器，最后采用频率转换法将低通滤波器转换成所希

8、望的各种滤波器。滤波器设计公式较多，计算繁杂。但是各种设计方法都有现成的设计程序设计函数供我们调用，只要掌握了滤波器基本设计原理，在工程实际中采用计算机辅助设计滤波器是常用的现象。2.2模拟滤波器频率变换模拟带通滤波器的设计 1确定“模拟滤波器”设计指标；2 指标转换为相应归一化低通模拟滤波器指标；3 选择类型，设计归一化低通然后去归一化得到低通模拟滤波器；4 频率变换，将转换为实际模拟滤波器系统函数。：需要设计的实际“AF”频响函数，：拉式复变量：实际AF的归一化频率，：实际AF的归一化拉式复变量：实际AF的归一化频响函数，：低通AF传输函数：低通AF拉式复变量，：低通AF归一化频率：低通A

9、F归一化拉式复变量，：低通AF归一化频响函数2.2.1带通滤波器的指标要求，指标转换公式：边界频率转换成低通的边界频率。带通滤波器频率特性是正负对称的，故这个变换必须是一对二的映射，它应该是的二次函数 p对应u，s对应s2，总结模拟带通的设计步骤：(1)确定模拟带通滤波器的技术指标，即：带通上限频率，带通下限频率下阻带上限频率，上阻带下限频率通带中心频率，通带宽度通带最大衰减为，阻带最小衰减为与以上边界频率对应的归一化边界频率如下：，(2) 确定归一化低通技术要求：，通带最大衰减仍为，阻带最小衰减亦为。(3) 设计归一化低通G(p)。(4) 直接将G(p)转换成带通H(s)。 2.2.2模拟低

10、通滤波器变换成数字带通滤波器由模拟低通滤波器变换成数字带通滤波器得变换关系为根据双线性变换，模拟带通与数字带通滤波器之间的S平面与Z平面的关系仍为，3设计内容要求通带范围为0.25rad0.45rad，通带最大衰减为3dB，阻带范围为00.15rad和0.55radrad，阻带最小衰减为40dB。3.1 设计思路根据题目要求，按照一定规则把给定的滤波器性能指标转化为模拟低通滤波器的技术指标，根据模拟滤波器技术指标设计为相应的模拟低通滤波器，根据双线性不变法把模拟滤波器转换为数字滤波器，即先设计过渡模拟滤波器得到系统函数，然后将按双线性变化法转换成数字滤波器的系统函数。3.2设计步骤1数字带通滤

11、波器技术指标通带上截止频率u=0.45rad，通带下截止频率l=0.25rad阻带上截止频率s2=0.55rad，阻带下截止频率s1=0.15rad通带内最大衰减p=3dB，阻带内最小衰减s=40dB。 2 模拟带通滤波器技术指标为了计算简单，可设T=1，则有，（通带中心频率）（带宽）将以上边界频率对带宽B归一化，得到，3 模拟归一化低通滤波器技术指标归一化阻带截止频率：归一化通带截止频率：p=1，p=3dB,s=40dB4 设计模拟低通滤波器，取N=7查表6.2.1，得到归一化低通传输函数G(p),5 模拟低通转换成模拟带通将归一化模拟低通转换成模拟带通6 数字带通滤波器通过双线性变换法将H

12、a(s)转换成数字带通滤波器H(z)。下面将(5)、(6)两步合成一步计算：将上式代入上面的转换公式，得将上面的p等式代入G(p)中，得3.3 MATLAB程序MATLAB程序如下：clear;close allt=1; fs=4000;wpu=0.45*pi,wpl=0.25*pi;wsu=0.55*pi,wsl=0.15*pi;wpz=0.25,0.45;wsz=0.15,0.55;wp=2/t*tan(wpz/2);ws=2/t*tan(wsz/2);rp=3;as=40;n,wc=buttord(wp,ws,rp,as,'s'); %计算带通滤波器阶数N和3dB截止频率

13、wcb,a=butter(n,wc,'s'); %计算带通滤波器系统函数分子分母多项式系数向量b,abz,az=bilinear(b,a,fs);nd,wdc=buttord(wpz,wsz,rp,as);bd,adz=butter(nd,wdc);hk=freqz(bd,adz)subplot(2,2,1)plot(abs(hk)xlabel('Hz'); ylabel(' 幅度'); title('幅度函数曲线')subplot(2,2,2)plot(angle(hk)xlabel('omega/pi'); y

14、label('phi(omega)'); title('相频特性曲线')hp,w=freqz(bd,adz,4000);subplot(2,2,3)plot(w/pi,20*log10(hp);axis(0,1,-100,1.5);xlabel('f/kHz'); ylabel(' -A(f)/dB'); title('损耗函数曲线')程序运行结果如下所示：3.4结果分析3.4.1损耗函数与相频特性常用的数字滤波器一般属于选频滤波器。其频率响应函数用幅频特性函数和相频特性函数表示。幅频特性表示信号通过该滤波器后各频

15、率成分振幅衰减情况，相频特性反映各频率成分通过滤波器后在时间上的延时情况。工程实际中通常用损耗函数（也称为衰减函数） 来描述滤波器的幅频响应特性，对归一化幅频响应函数，定义如下（其单位是分贝，用dB表示）：。应当注意，损耗函数和幅频特性函数只是滤波器幅频响应特性的两种描述方法。损耗函数的优点是对幅频响应的取值非线性压缩，放大了小的幅度，从而可以同时观察通带和阻带频响特性的变化情况。由程序运行结果可知，此损耗函数曲线符合带通滤波器的特性。滤波的过程中，信号通过滤波器，在一定的条件下，可以控制原信号的幅值不发生改变，但是会发生延迟，也就是输出信号的相角与输入信号之间会有一定的差值，在整个过程中，就

16、表现为相频特性。本课题研究的内容中，输出信号与输入信号的相角差值，表现出来的相频特性符合系统的特征。3.4.2双线性变换脉冲响应不变法的主要缺点是会产生频谱混叠现象，使数字滤波器的频响偏离模拟滤波器的频响特性，为了克服这一缺点，可将整个S平面压缩到S平面的一条横带里，通过标准变换关系将此横带变换到整个Z平面上去。由此建立S平面与Z平面一一对应的单值关系，消除多值性，也就消除了混淆现象。S平面与Z平面的单值映射关系：和，即为双线性变换。双线性变换法的主要优点是不存在频谱混迭。由于S平面与Z平面一一单值对应，S平面的虚轴(整个j)对应于Z平面单位圆的一周，S平面的=0对应于Z平面的=0， 对应 ，

17、即数字滤波器的频率响应终止于折迭频率处，所以双线性变换不存在频谱混迭效应。双线性变换具有将s的左半平面映射到z平面单位圆内;j轴映射到单位圆上的基本性质。当=0时，z=1，说明s平面j轴映射到z平面单位圆上；当<0时，上式中的分母大于分子z<1，说明s左半平面映射到z平面单位圆内。因而它们一一对应，有着单值关系，是一种保角交换。因此一个稳定的模拟滤波器，通过双线性变换只能得到一个（唯一）稳定的数字滤波器。由此可见，双线性变换法将频带严格限制在s/2范围内，从根本上消除了频谱混叠。双线性变换法不仅解决了直接z变换中存在的混叠误差，而且它还是个代数变换，只要将s与z之间存在的简单代数关

18、系直接代入模拟滤波器系统函数式中，就可以求得相应的数字滤波器的系统函数和频响，模拟滤波器所具有的优良特性就得以保存。但是双线性变换法由于从ss的变换，频率关系不是线性的，它只有在低频段接近线性，频率越高被压缩得越厉害，因而出现频率非线性畸变（扭曲），原模拟滤波器的转折点频率的位置因非线性而产生了变化，相频特性也因非线性而出现了畸变。此次巴特沃斯数字带通滤波器的设计时对应于双线性变化法，如上所述，双线性变化是一种非线性转换。双线性变换法的缺点：与的非线性关系，导致数字滤波器的幅频响应相对于模拟滤波器的幅频响应有畸变，另外，一个线性相位的模拟滤波器经双线性变换后，滤波器就不再有线性相位特性。虽然双线性变换有这样的缺点，但它目前仍是使用得最普遍、最有成效的一种设计工具。双线性变换比脉冲响应法的设计计算直接和简单。由