1、宀完 成 内 容收集消化资料、学习 MATLAB软件， 进行相关参数计算编写仿真程序、调试指导教师签字基层教学单位主任签字第一章 前言 4第二章数字信号处理基础知识部分 52.1巴特沃斯滤波器的设计 52.2数字滤波器一一IIR 9第三章 matlab 部分基础知识 143.1MATLAB 介绍 143.2巴特沃斯 matlab函数介绍 15第四章 仿真部分及仿真图 164.1数字滤波器的设计步骤 164.2MATLAB 程序 174.3仿真结果 18第五章 结论 19第六章参考文献 20第一章 前言数字信号处理课程设计是在学生完成数字信号处理和 MATLAB的结合后的基本实验以后开设的。本课
2、程设计的目的是 为了让学生综合数字信号处理和 MATLAB并实现一个较为完整 的小型滤波系统。滤波器是自动控制、信号处理和通信领域的重要组成部分 ,广泛地应用于各种系统中。MATLAB语言是一种简单、高效的高 级语言,是一种内容丰富、功能强大的分析工具，其应用范围几乎 覆盖了所有的科学和工程计算领域。MATLAB中提供了丰富的用 于模拟滤波器设计的函数，通过编程可以很容易实现低通、高通、 带通、带阻滤波器，并能画出滤波器的幅频特性曲线，大大简化了 模拟滤波器设计。本文将通过利用 MATLAB滤波器设计函数直 接实现巴特沃斯低通滤波器的设计，找到应用 MATLAB来设计低通滤波器的方法。介绍了用
3、脉冲响应不变法设计巴特沃斯低通 滤波器的基本理论和设计思想，给出了基于 MATLAB设计巴特沃斯低通滤波器的具体步骤，文中还给出了利用 MATLAB产生的滤波器的幅频及相频特性曲线。详细介绍了在基于 MATLAB设计低通滤波器过程中常用到的工具和命令。第二章数字信号处理基础知识部分2.1巴特沃斯滤波器的设计1 .巴特沃斯滤波器的幅度平方函数及其特点巴特沃斯模拟滤波器幅度平方函数的形式是12N1 - （ / jc）中N为整数，是滤波器的阶次。0 =0时，| H a（ j。）=1时；当0 =0c 时，Ha（jCc）=1/p2，所以称0 c为半功率点。此时通带内衰减 6=-20|g Ha（j0c）
4、=3dB ， Qc又被称为 3dB截止频率。巴特沃斯滤波器的振幅特性如图2.1.1所示。由图可见，当门c 时，随门加大，幅度迅速下降。下降的速度与阶次 N有关，N越 大，幅度下降速度越快，过渡带越窄。图2.1.1 巴特沃斯幅度特性和N的关系2幅度平方函数的极点分布及H a（S）的构成2将幅度平方函数H a （j0）写成S的函数Ha（S）Ha（-S）二 2N1 +（ j0/ j0c）“ *1（S/jc）2N（ 2-1-2）1 2k 1j （亠 ）i.e 2 2N k=0,1,2 ,此式表明幅度平方函数有2N个极点，极点Sk用下式表示2 k 仃1/2N 小 加 2N） j2cSk -（ -1） （
5、ji ）=e *e 2i 黒2N-1 （2-1-3 ）这2N个极点分布在s平面半径为11的圆（称为巴特沃斯圆）上,角度间隔是n /N弧度。N=3时，极点间隔为n /3弧度或60度，如图2.1.2所示。图2.1.2 三阶巴特沃斯滤波器极点分布 极点对虚轴是对称的，且不会落在虚轴上。当 N是奇数时，实轴 上有极点；当N为偶数时，则实轴上没有极点。巴特沃斯滤波器 的N个极点为1 2k JL j ：】（ ）Sk -2 2N ；k =0,1,“ ：，N -1 （2-1-4 ）则H a （S）的表达式即滤波器的系统函数为NHa（s）二 一（2-1-5 ）I （S - Sk）k _03.频率归一化问题式（2
6、-1-5 ）即为所求滤波器的系统函数，可看出 Ha（s）与门c有 关，即使滤波器的幅度衰减特性相同，只要不同，Ha（S）就不 一样。为使设计统一，可将所以的频率归一化。这里采用对 3dB 截止频率“C见归一化，归一化后的Ha（S）表示为Ha（S）= （ 2-1-6 ），令九=乡幺 ，称人为归一化S Sk 八i【（.一.）0 匚： c - - c频率，则S门=j门门=j ，记为p称归一化复变量，这样归一化巴特沃斯的传输函数为 Ma（p）二 一1 （2-1-7 ）式中Pk为丨丨（P - Pk）k -0归一化极点，再按照式（2-1-7）得到归一化的传输函数:a（p）4阶数N的确定阶数N的大小主要影响
7、幅度特性下降速度，它应该由： p,s和0s确定。令扎。=%p，则N由下式表示:k$pN （2-1-9）关于3dB截止频率I%,如果技术指标未给出，则按下式求出:0.1 - c p（io p -1） 2N （2-1-10 ）。5.设计步骤总结以上讨论，低通巴特沃斯滤波器的设计步骤如下:1） 根据技术指标:ps/S和、，用式（2-1-9 ）求出No2） 按照（2-1-8 ），求出归一化极点pk，代入（2-1-7 ），得到归 一化传输函数Ha（p）。也可以直接查表.3） 将Ha（P）去归一化。将P = S/l.代入Ha（P），得到实际的滤波 器传输函数H a （S）。如果技术指标没有给出3dB截止频
8、率-c，可 以按照式（2-1-10 ）求出。2.2数字滤波器一一IIR1.IIR滤波器基础IIR滤波器的设计就是在给定了滤波器的技术指标后，确定滤波 器的阶数N和系数；a,b二在设计IIR滤波器时，常用的方法是 利用模拟滤波器来设计数字滤波器，设计方法是首先将数字滤波 器技术指标转换为对应的模拟滤波器技术指标，然后设计满足技 术指标的模拟滤波器H a S ,最后将设计出的模拟滤波器H a（S）转 换为满足技术指标的数字滤波器 Hz。将系统函数H（z）从S平 面转换到Z平面的方法很多，主要有两种：一种是使数字滤波器 的h n近似于模拟滤波器的ha t ，可导出脉冲响应不变法；另一 种是使数字滤波
9、器的差分方程近似于滤波器的微分方程，由此可 导出双线性变换法。为保证转换后的H z稳定且满足技术要求，对转换关系提出以下 两点要求：（1）因果稳定的模拟滤波器转换成数字滤波器，仍是因果稳定的。 模拟滤波器因果稳定要求其系统函数 H a S的极点全部位于S平 面的左半平面；数字滤波器因果稳定则要求H z的极点全部在单 位圆内。因此，转换关系应是 S平面的左半平面映射到Z平面的 单位圆内。（2） 数字滤波器的频率响应模仿模拟滤波器频率响应， S平面的 虚轴映射z平面单位圆。2脉冲响应不变法脉冲响应不变法的核心是通过对连续函数 ha （t）等间隔采样得 到离散序列ha nT，使h（n）二ha（nT）
10、（其中T为采样间隔），因此 脉冲响应不变法是一种时域上的转换方法，转换步骤如下设模拟滤波器Ha（S）只有单阶极点，且分母多项式的阶次高于分子多项式的阶次，将Ha（S）用部分分式表示，则Ai（2-2-1）式中，Si为Ha（S）的单阶极点。将Ha（S）进行拉氏逆变换得到s tha（t） Aeu（t） （ 2-2-2）i式中，u（t）单阶阶跃函数。对ha（t）进行等间隔采样，采样间隔为T，得到l s n Th（n）=ha（nT ）八 Ai e u（nt） （2-2-3）i仝对上式进行Z变换，得到数字滤波器的系统函数N aH （z） sT_7 （ 2-2-4）i 二 1 e z由这一转换过程看出，它对
11、部分分式表达的模拟系统函数更为方便，对任一极点Si，Ha（S）到H（z）得转换可直接用下式来完成（2-2-5）Ai A sT _Ls - Sj 1 - e z将式（2-2-1）的Ha S和式（2-2-4）的H z加以比较，可以看出:1）S平面的单极点s 7变换到z平面上z =esT处的单极点。2）Ha（s）与H（z）的部分分式的系数是相同的，都是 Ai3）如果模拟滤波器是稳定的，所有极点Si位于s平面的左半平面， 及极点的实部小于零，则变换后的数字滤波器的全部极点在单位 圆内，即模小于1，因此数字滤波器也是稳定的。3. H （ej- ）与 H a（）的关系下面分析从模拟滤波器转换到数字滤波器，
12、 S平面和z平面之间的映射关系，从而找到这种转换方法的优缺点。这里以采样 信号ha（t ）作为桥梁，推导其映射关系。将ha（t）的采样信号用ha （t）表示，toh?a（t） = 、 ha（t）：.（t- nT ）n对ha（t）进行拉式逆变换，得到H?a（S）st s n Tdt 八 ha （nT ）e - n式中，ha（ nT）是ha （t）在采样点t=nT时的幅度值，它与序列的幅度相等，即h（n）=ha（nT ），因此得到Ha（s）=送 h（n）eT =送 h（n）zzr =H （z） zt （2-2-6）n n上式说明采样信号的拉氏变换与相应的 Z变换之间的映射关系为sT z = e （
13、2-2-7）已知模拟信号ha（t）的傅立叶变换Ha（V）和其采样信号ha（t）的傅 立叶变换Ha（r）之间的关系满足1旳Ha（j） Ha（j jks） （ 2-2-8）T k n :：其中，嘔=2二/T，将S=B 1带入上式，得H a（s） H a（S jks） （2-2-9 ）T kH （z） Ha（s jks） （2-2-10 ）上式表明将模拟信号ha（t）的拉氏变换在s平面上沿虚轴按照周期.-s=2：-JT延拓后，再按上式的映射关系映射到 Z平面上，就得到H（z）。 Z二esT可称为标准映射关系。Ret图2.2.1冲击响应不变法的映射关系-C - j-.1jw二 re则得到:（2-2-1
14、1）=1 0, r那么 虚轴映射z平面的单位圆上，s平面右半平面映射z平面单位圆 外。这说明：第一，如果Ha（s）因果稳定，转换后得到H（z）仍因 果稳定；第二，数字滤波器频率响应H（eJ模仿模拟滤波器的频 率响应H a（j），满足转换关系的两点要求。H ej-与Ha（j）二者关系可由式（2-2-10）导出，考虑s =厂，且 T ， 则 得 到1 旳 1 旳 閒_2nkH e二 Ha - jk）二 a H a（j ） （2-2-12）。这T T T就是说，数字滤波器的频率响应是模拟滤波器频率响应以 s为周期的周期延拓。因而正如采样定理所讨论的，只有当模拟滤波器的频率响应是有限的，且带限于折叠频
15、率以内时，即H a（j0 ）=0, C兰兀，才能使数字滤波器的频率响应在折叠T 2频率以内重现模拟滤波器的频率响应，而不产生混叠失真，即沁 1 H （e ） =H a（j ）, 兀。T T4频率混叠现象实际上，任何一个模拟滤波器的频率响应都不是严格带限的，变 换后就会产生周期延拓分量的频谱交叠，即产生频率响应的混叠 失真。如果原模拟信号ha t的频带不是限于+/- n /T之间，则会在 他的奇数倍附近产生频率混叠，从而映射到 z平面上，在3 =+/-n附近产生频率混叠。脉冲响应不变法的频率混叠现象如下图所 示5.脉冲响应不变法的优点脉冲响应不变法的优点是频率坐标变换是线性的， 即 T。另一个优
16、点是数字滤波器的单位脉冲响应完全模仿模拟滤波器的单 位脉冲响应，时域特性逼近好。第三章 matlab部分基础知识3.1 MATLAB 介绍设计平台MATLAB7.1 MATLAB是美国MathWorks公司生产 的一个为科学和工程计算专门设计的 交互式大型软件，是一个可以完成各种精确计算和数据处理的、可视化的、强大 的计算工具。 它集图示和精确计算于一身，在应用数学、物理、化工、机 电工程、医学、金融和其他需要进行复杂数值计算的领域得到了 广泛应用。它不仅是一个 在各类工程设计中便于使用的计算工 具，而且也是一个在数学、数值分析和工程 计算等课程教学中的优秀的数学工具，在世界各地的高和大型计算
17、机上运行， 适 用于Windows、UNIX 等多种系统平台。 MATLAB 作为一种科学计算的高级语言之所以受欢迎， 就是因为它有丰富的 函数资源和工具箱资源，编程人员可以根据自己的需要选择函数，而无需再 去编 写大量繁琐的程序代码，从而减轻了编程人员的工作负担， 被称为第四代编程语言的MATLAB最大的特点就是简洁开放 的程序代码和直观实用的开发环境3.2巴特沃斯matlab函数介绍1.buttordN,wc=buttord（wp，ws， a p， a s）用于计算巴特沃斯数字滤波 器的阶数N和3dB截止频率wc。调用参数wp，ws分别为数字 滤波器的通带、阻带截止频率的归一化值，要求：
18、OWwp 1，0 ws 1且对于低通滤波器 wpvws。1表示数字频率pi。2.butterb， a=butter （N，wc， ftype）计算N阶巴特沃斯数字滤波器系统函数分子、分母多项式的系数 向量b、a。调用参数N和wc分别为巴特沃斯数字滤波器的阶数和 3dB截止 频率的归一化值（关于pi归一化），一般是调用N,wc=buttord（wp， ws， a p， a s）格式计算N和wc。系数b、a是按照z-1的升幕排 列。3.freqzH , 1= freqz B, A, N, Fs计算离散系统在0二范围内的N个频率等分点的频率响应的值。4.absAbs （H）为取H得幅度。5.an g
19、leAngle （ H）为取H的相位角，单位为弧度，用pha显示计算的相 位。6.axisaxis（xmin xmax ymin ymax）中分别给出x轴和y轴的最大值、 最小值7.subplotsubplot（m,n,p）将多个图画到一个平面上的工具。 其中，m表示是 图排成m行，n表示图排成n列，p是指你现在要把曲线画到figure 中哪个图上。8.plotplot（x,y）是绘制二维图形的最基本函数， 以x元素为横坐标值，y元素为纵坐标值绘制曲线。9.Xlabel x 轴注解Ylabel y轴注解Title 标题注解Grid on 开启网格格式第四章 仿真部分及仿真图4.1数字滤波器的设
20、计步骤设计流程图4.2MATLAB 程序Fs=1000; fp=100; fs=300; Ap=3; As=20;Wp=fp/（Fs/2）;Ws=fs/（Fs/2）; N,Wc=buttord（Wp,Ws,Ap,As）; b,a=butter（N,Wc）;H,F=freqz（b,a,500,Fs）;subplot（2,2,2） plot（F,20*log10（abs（H） xlabel（频率（Hz）;ylabel（幅度（dB） axis（0 500 -40 3）;title（ 数字滤波器分频曲线grid onsubplot（2,2,1）plot（F,abs（H）;xlabel（幅度 数字滤波器
21、幅频曲线grid on;subplot（2,2,3）pha=a ngle（H）;plot（F,pha）;相位（rad） 数字低通滤波器相频曲线4.3仿真结果第五章结论本次课程设计的数字低通滤波器首先要设计出相应的巴特沃 斯低通滤波器，由我们课上学过的知识，根据所给的设计参数确 定滤波器的阶数，归一化等。再由通过matlab得出的模拟滤波器， 得到相应的数字滤波器。经过一周的课程设计，让我再次熟悉了所学过的关于数字信 号处理中设计滤波器部分的知识，同时，由于设计需要，我通过 查阅资料，自学了用于本次课设的 matlab程序仿真部分的知识。 MATLAB在以前的学习中并没有怎么接触，但是通过网上查
22、找资 料，我很快的理解了其基本原理。 总的来说，通过这次的课程设 计我对MATLAB有了全面的认识，对数字信号处理的知识又有 了深刻的理解，让我感受到只有在充分理解课本知识的前提下， 才能更好的应用这个工具；并且熟练的应用 MATLAB也可以很好的加深我对课程的理解，方便我的思维。通过本次数字信号来处理课程设计，我通过系统设计，软件 仿真，程序安排与调试，初步掌握了工程设计的具体步骤和方法， 提高了我的分析问题，解决问题的能力，提高实际应用水平。感谢学院给了我们这一次课程设计的机会，让我们受益匪浅， 感谢老师们在我们遇到问题时耐心的指导，希望今后能有更多的 这种亲身实践的机会。机械工业出版社，2009 清华大学出版社，2003 电子工艺出版社，2005第六章 参考文献谢平，王娜. 数字信号处理.薛年喜.MATLAB信号处理中的应用MATLAB7辅助信号处理与应用燕山大学课程设计评审意见表指导教师评语：成绩：年 月日答辩小组评语：评阅人：课程设计总成绩：答辩小组成员签字：年月 日
