带通切比雪夫带通滤波器的设计.doc

燕山大学课 程 设 计 说 明 书 题目：切比雪夫带通滤波器的设计 学院（系）： 电气工程学院年级专业： 08级检测1班学 号： 学生姓名： 指导教师： 教师职称： 讲师 _燕山大学课程设计（论文）任务书院（系）：电气工程学院 基层教学单位： 仪器科学与工程系 学 号080103020098学生姓名王娟专业（班级）08级检测1班设计题目14 切比雪夫带通滤波器的设计设计技术参数带宽：200Hz；中心频率：1000Hz；通带内衰减：3dB；频率小于830Hz或大于1200Hz处衰减不小于25dB。设计要求求该带通滤波器转移函数H(s)及幅频相频特性。参考资料数字信号处理方面资料MATLAB方面资料周次第一周第二周应完成内容收集消化资料、学习MATLAB软件，进行相关参数计算 编写仿真程序、调试指导教师签 字 王娜基层教学单位主任签字说明：此表一式四份，学生、指导教师、基层教学单位、系部各一份。2011年6月18日目录第一章摘要- 3 -第二章 引言·················································3 第三章 切比雪夫带通滤波器的设计思路························· -3- 3.1 设计的基本思路································· 4 3.2 切比雪夫滤波器的基本知识······················· 4 3.3 带通滤波器的设计步骤··························· 8第四章 基于MATLAB的切比雪夫带通滤波器程序设计··········· -9-4.1 仿真用到的命令介绍·····························- -9-4.2 仿真过程与仿真图······························ - 10-4.3 仿真图结论及误差分析························· - 13-第五章 心得体会········································· -13-参考文献··················································· -14-附录I ···················································· - -15-评审意见表················································- - 16- 第一章 摘要MATLAB是当今最优秀的科技应用软件之一，它简单易用，具有强大的科学计算能力、可视化功能、开放式可扩展环境，所附带的工具箱支持30多个领域。它在信号处理、通信、自动控制及科学计算等领域被广泛应用，被认为是最能够提高工作效率、改善设计手段的工具软件，掌握了MATLAB，就掌握了一把开启科学研究领域大门的钥匙。本次课设是基于MATLAB的切比雪夫模拟带通滤波器的设计，全文首先对带通滤波器的前景作了简明的阐述，通过介绍切比雪夫滤波器设计的基本原理，以及在此基础上带通滤波器的设计思路，最终应用MATLAB进行了编程，通过仿真得到了老师要求的参数及图像，并对仿真结果进行了分析总结。第二章 引言模拟带通滤波器在信号检测和传输过程中起着很重要的作用，但是传统的模拟带通滤波器的设计过程中存在着大量繁琐的运算数值计算，如果手工计算则费事又费力。现在伴随着MATLAB的广泛使用，给切比雪夫带通滤波器的设计带来了方便，用户通过已经编写好的程序，只要输入预期的参数值就可以直接求出最终结果。目前高通、带通、带阻滤波器的设计思路本质上都是将给定的技术指标通过某种频率变换关系先转换成模拟低通滤波器的技术指标，并依据这些技术指标求出低通的转移函数，然后再通过频率变换关系转换成所要涉及的带通转移函数。第三章 切比雪夫带通滤波器的设计思路3.1 设计的基本思路对于切比雪夫模拟带通滤波器的设计思路，就是依据给出的带通技术指标通过某种频率变换，转换成模拟低通滤波器的技术指标，并依据这些技术指标求出低通的转移函数，最后再经频率变换关系转换回去。设计流程如图3.1所示。图3.1 切比雪夫带通滤波器设计流程图3.2 切比雪夫滤波器的基本知识模拟滤波器的理论和设计方法己发展得相当成熟，且有一些典型的模拟滤波器供我们选择，如巴特沃斯(Butterworth)滤波器、切比雪夫(Chebyshev)滤波器、椭圆(Cauer)滤波器、贝塞尔(Bessel)滤波器等，这些典型的滤波器各有特点。巴特沃兹滤波器在通带内幅度特性是单调下降的，如果阶次一定，则在靠近截止 处，幅度下降很多，或者说，为了使通带内的衰减足够小，需要的阶次N很高，为了克服这一缺点，采用切比雪夫多项式来逼近所希望的 。切比雪夫滤波器的 在通带范围内是等幅起伏的，所以在同样的通常内衰减要求下，其阶数较巴特沃兹滤波器要小。切比雪夫滤波器特点：误差值在规定的频段上等波纹变化。切比雪夫滤波器是一类重要的归一化模拟低通滤波器，式（3.1）中是次切比雪夫多项式。次切比雪夫多项式是如下定义的： (3.1) 尽管在上式中，表达式分成了两部分，并且各部分都不是多项式的形式，但它却表示了一个次多项式。例如，若=2，1时 （ 1）（3.2）而1时 （ 1）（3.3）无论1还是 1，的表达式都相同，这是一个二次多项式。对=0和1的简单情形，由（3.1）式容易得到 （3.4） （3.5） 由于(3.1)式中两部分表达式给出的结果是相同的，因此可以从任一式出发，得到切比雪夫多项式的具体表达式 。对1的情形，可以应用如下的递推公式计算各次切比雪夫多项式： （1） （3.6）例如 (3.7) 这与前面由定义式得到的结果完全一致。类似地 (3.8)对切比雪夫多项式作了简单介绍后，下面讨论切比雪夫幅方函数。切比雪夫振幅平方函数为 （3.9） 式中c为有效通带截止频率，表示与通带波纹有关的参量，值愈大通带波动愈大。是N阶切比雪夫多项式，在（3.1）式已介绍过。切比雪夫滤波器的振幅平方特性如图3-2所示：N为偶数，得到min， （3.10） N为奇数，得到max， （3.11） 图3-2 切比雪夫滤波器的振幅平方特性 有关参数的确定:  a、通带截止频率： 预先给定; b、与通带波纹有关的参数,通带波纹表示成 （3.12）由上式可以看出，切比雪夫滤波器的通带最大波纹与常数有关。若给出不同的值，滤波器就可以有不同的通带最大衰减。这种情形，与巴特沃兹滤波器不同。因此，称为切比雪夫滤波器的通带衰减常数。另一方面，切比雪夫滤波器的通带最大衰减与无关，这与巴特沃兹滤波器是一致的。若已知通带波纹值，则可以求得    （3.13） 、阶数N:由阻带的边界条件确定。 、A2为事先给定的边界条件,即在阻带中的频率点处 ,要求滤波器频响衰减到1/A2以上。当时， （3.14） （3.15） 1时， 得 （3.16） 因此，要求阻带边界频率处衰减越大，要求N也越大，参数N， 给定后，查阅有关模拟滤波器手册，就可求得系统函数Ha(s)。3.3 带通滤波器的设计步骤前面已经讨论过如何去设计带通滤波器的基本流程，介绍完切比雪夫滤波器的相关内容，我主要针对本次的课设是带通切比雪夫滤波器的设计来说明一下带通滤波器的设计步骤：1） 确定模拟带通滤波器的技术指标，即带通上限频率，带通下限频率；下阻带上限频率，上阻带下限频率；通带中心频率，通带宽度；与以上边界频率对应的归一化边界频率如下，； （3.17）还需要确定的技术指标有：通带最大衰减，阻带最小衰减。2） 确定归一化低通技术要求 （3.18）与的绝对值可能不相等，一般取绝对值小的，这样保证在较大的处更能满足要求。通带最大衰减仍为，阻带最小衰减亦为。3） 设计归一化低通G(p)。4） 由直接将G(p)转换成带通。第四章 基于MATLAB的切比雪夫带通滤波器的程序设计4.1 仿真用到的命令介绍通过对切比雪夫带通滤波器的设计思路的阐述，由于切比雪夫带通滤波器的设计需要大量的计算，计算如此繁琐，我们通常使用MATLAB软件进行编程，这样只要编写好程序，输入我们设定的参数，我们就可以很快地达到我们想要的图像和参数。下面我就着重介绍一下怎么才能通过程序实现我们要求来求出我们想要的转移函数和幅相曲线。首先我们在设计切比雪夫低通滤波器之前已经给定了一些参数，所以程序的第一行就是要给出初始化参数，即B=200;Wo=1000;Wsl=830;Wsh=1200;Rp=3;Rs=25;此处就与第三章介绍的方法有些区别了，因为我们知道，所有带通设计参数要转化成低通参数，但是我们在些命令时就不用具体转化，因为我们接下来用到的cheb1ord函数就可以自动判断给出的参数是什么类型切比雪夫滤波器的参数。这就是用软件编程的好处！接下来就是以下两句：Wp = Wo-B/ 2 Wo + B/ 2 * 2*pi ; %求目的是出通带上下限频率，注意此处与理论值还是有些许误差的。（理论算法在第三章已经介绍）Ws = Wsl Wsh * 2 * pi ; %求出阻带上下限频率，此处的频率是两个。 n ,Wn = cheb1ord (Wp ,Ws ,Rp ,Rs ,'s') b ,a = cheby1 (n ,Rp ,Wn ,'s')%这两句就是得到滤波器的阶数、转移函数和分子、分母系数。W= linspace (1 ,3* Wo ,1000) * 2* pi ;%等分横坐标轴上点，便于观察得出的图像H = freqs (b ,a ,W) ; %求出幅频响应，此处是用一个求幅频响应的函数，因为之前b、a已经求出，所以我们可以直接得到函数的幅频响应。magH = abs ( H) ;%此函数是用来求得幅值。angH=angle(H);% 求得相角。subplot(1,2,1);%此处是要绘制一行两列的图像，最后一个一代表是第一个图像。plot (W/ (2 * pi) ,20 * log10 (magH) ,'b') ; %绘制出幅频响应特性曲线title ('切比雪夫I 型模拟带通滤波器幅频特性曲线') ;%加注所绘制的滤波器的特性曲线名称xlabel ('频率( Hz)') ; %加注横坐标轴标识ylabel ('幅度(dB)') ; %加注纵坐标轴标识axis(500,1500,-50,0) %限定曲线范围set(gca,'XTickMode','manual','XTick',830,900,1000,1200);%画出栅格线，便于观察是否达到要求的参数。set(gca,'YTickmode','manual','YTick',-30,-25,-3,0); grid %功能同上。subplot(1,2,2);plot (W/ (2 * pi) ,angH,'b') ; %绘制出幅频响应特性曲线title ('切比雪夫1 型模拟带通滤波器相频特性曲线');%加注所绘制的滤波器的特性曲线名称xlabel ('频率( Hz)') ; %加注横坐标轴标识ylabel ('相角') ; %加注纵坐标轴标识4.2 仿真过程与仿真图将上述程序输入MATLAB程序仿真窗口中，如图4-1所示的仿真界面所示，方法是新建一个仿真文件，将程序输入到editor中，最后保存并运行程序，得到图4-2的仿真图图4-1 仿真程序界面图4-2 仿真出的幅频相频特性曲线同时我们也得从command window中得到了如图4-3所示图4-3但是此处的a，b是不完全显示的，图4-4为最终的b，a值。图4-44.3 仿真图结论及误差分析通过以上的仿真图及仿真数据，我很成功地达到了老师要求的仿真参数要求，从仿真图4-2图上可以很明显地观察到切比雪夫I型幅频特性曲线在频率小于830Hz或大于1200Hz处阻带衰减是不小于25dB的，其中在830HZ处是-25dB的，在1200HZ处阻带衰减大于25dB的。而在900Hz处阻带衰减达到最大，正好是-3dB的，这也非常符合设计要求。因此这种结果是非常符合设计要求的，同时我们也很明确地知道我们所用的切比雪夫I型带通滤波器的阶次为3阶，而通常如果用巴特沃兹滤波器设计的话阶次可能会更高些。 通过得出的b，a ，以及阶数n 求得带通滤波器的转移函数H(s)为 （4.1） 误差分析：本次误差主要Wp = Wo-B/ 2 Wo + B/ 2 * 2*pi ;目的是出通带上下限频率，但是在理论上我们的计算公式是，。所以在此处会有误差存在，但是一般不影响我们最终结果。第五章 心得体会通过这次的课程设计，我深刻地体会到：数字信号是一门应用性很强的学科，并且通过结合MATLAB可以实现很多的信号的处理功能。我通过这次的课程设计更好地锻炼了自己对数字信号的理解与MATLAB的运用。课程设计是培养我们综合运用所学知识,发现、提出、分析和解决实际问题，锻炼实践能力的重要环节,是对我们实际工作能力的具体训练和考察过程，为自己今后从事与数字信号有关的工作打下了基础。我从一开始拿到自己的题目,我是很茫然的，因为课本上没有过多的关于切比雪夫滤波器设计的知识，所以我首先借阅了关于切比雪夫滤波器方面的书籍，仔细学习了关于切比雪夫滤波器的基本知识，对切比雪夫滤波器有了个大概的了解。并且知道有两种切比雪夫滤波器形式，最后我选择了切比雪夫I型。还有就是老师给定的好多参数，我发现我对基本概念的理解确实不扎实，起初都不知道说的到底是哪个参数，究竟在图上对应的哪个参数，所以我也很认真地翻阅了课本上有关这些知识的地方，从而进一步理解了带通滤波器的设计思路。带通滤波器最终要将有关参量转化到低通然后再去计算转移函数，最后通过频率变换算处带通的转移函数，虽然有思路，但是如何编写程序还是一个棘手的问题，毕竟自己以前没有学过MATLAB，在短短几天里真的很难学的透彻，所以我只是选择性地学习了如何求切比雪夫滤波器转移函数分子分母以及阶数的算法函数，再加上如何指标作出曲线的程序段，最终实现了课设的要求。最后我发现，课设当中虽然有遇到瓶颈的时候，但是自己一定不要放弃，同组同学的启发也是很重要的，所以我也体会到团队的重要性，还有我认为自己哪不懂也不能一味地去寻求外界的帮助，毕竟自己理解的才深刻，多想自己发问，多跟书本问为什么，而且一定要静下心来认真思考，养成解决问题的能力是最重要的。我也非常感谢课设中老师对我的热心帮助！ 参考文献 冷建华 李萍 王良红 数字信号处理 国防工业出版社 2002年7月 郭仕剑 王宝顺 贺志国 杨可心等 MATLAB7.x数字信号处理 人民邮电出版社 2006年12月 谢平 王娜 林洪彬 信号处理原理及应用 机械工业出版社2008年1月 附录I ： 设计程序B=200;Wo=1000;Wsl=830;Wsh=1200;Rp=3;Rs=25;Wp = Wo-B/ 2 Wo + B/ 2 * 2*pi ; Ws = Wsl Wsh * 2 * pi ; n ,Wn = cheb1ord (Wp ,Ws ,Rp ,Rs ,'s') b ,a = cheby1 (n ,Rp ,Wn ,'s')W= linspace (1 ,3* Wo ,1000) * 2* pi ;H = freqs (b ,a ,W) ; magH = abs ( H) ;angH=angle(H);subplot(1,2,1);plot (W/ (2 * pi) ,20 * log10 (magH) ,'b') ; title ('切比雪夫1 型模拟带通滤波器幅频特性曲线') ;xlabel ('频率( Hz)') ; ylabel ('幅度(dB)') ; axis(500,1500,-50,0)set(gca,'XTickMode','manual','XTick',830,900,1000,1200);set(gca,'YTickmode','manual','YTick',-30,-25,-3,0);gridsubplot(1,2,2);plot (W/ (2 * pi) ,angH,'b') ; title ('切比雪夫1 型模拟带通滤波器相频特性曲线');xlabel ('频率( Hz)') ; ylabel ('相角') ; 燕山大学课程设计评审意见表指导教师评语：该生学习态度 （认真 较认真 不认真） 该生迟到、早退现象 （有 无）该生依赖他人进行设计情况 （有 无）平时成绩： 指导教师签字： 2011 年7 月 9 日图面及其它成绩：答辩小组评语：设计巧妙，实现设计要求，并有所创新。 设计合理，实现设计要求。 实现了大部分设计要求。 没有完成设计要求，或者只实现了一小部分的设计要求。 答辩成绩： 组长签字： 2011 年 7 月9日课程设计综合成绩：答辩小组成员签字： 2011年 7月9日