matlab的带通带阻滤波器的设计 精品.doc

1、 MATLAB简介1、 MATLAB介绍MATLAB是由美国公司发布的主要面对科学计算、可视化以及交互式的高科技计算环境。它将、科学数据可视化以及非动态系统的和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效的众多科学提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式（如C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。 MATLAB的基本数据单位是矩阵，它的与、工程中常用的形式十分相似，故用MATLAB来解算问题要比用C，等语言完成相同的事情简捷得多，并且MATLAB也吸收了像Maple等软件的优点，使MATLAB成为一个强大的。在新的版本中也加入了对，的支持。可以直接调用,用户也可以将自己编写的实用程序导入到MATLAB库中方便自己以后调用，此外许多的MATLAB爱好者都编写了一些经典的程序，用户可以直接进行下载就可以用。 2、 功能特性主要功能 1.数值分析2.数值和符号计算3.工程与科学绘图4.控制系统的设计与仿真5.数字图像处理6.数字信号处理7.通讯系统设计与仿真8.财务与金融工程3、 应用MATLAB 产品族可以用来进行以下各种工作： 1、 数值分析2、 数值和符号计算3、 工程与科学绘图4、控制系统的设计与仿真5、技术6、技术7、通讯系统设计与仿真8、与9、管理与调度优化计算（运筹学）MATLAB 的应用范围非常广，包括信号和、通讯、控制系统设计、测试和测量、和分析以及计算学等众多应用。附加的工具箱（单独提供的专用MATLAB集）扩展了MATLAB 环境，以解决这些应用内特定类型的问题。1、 设计原理 模拟带通、带阻滤波器的设计过程是：先将希望设计的各种滤波器转换为低通滤波器技术指标，然后设计相应的低通滤波器，最后采用频率转换法将低通滤波器转换成所希望的各种滤波器。1、带阻滤波器的设计原理 本设计中分别用ha(t)、Ha(s)、Ha(j)表示模拟滤波器的单位脉冲相应、系统函数、频率响应函数，三者的关系如下Ha(s)=LTha(t)=ha（t）exp(-st)dtHa(j)=FTha(t)=ha（t）exp(-jt)dt 可以用ha(t)、Ha(s)、Ha(j)中任一个描述模拟滤波器，也可以用线性常系数微分方程描述模拟滤波器。但是设计模拟滤波器时，设计指标一般由幅频相应函数|Ha(j)|给出，而模拟滤波器设计就是根据设计指标，求系统函数Ha(s)。 工程实际中通常用所谓的损耗函数即衰减函数A（）来描述滤波器的幅频响应特性，对归一化幅频响应函数，A（）定义如下：A（）=-20lg|Ha(j)|=-10lg|Ha(j)|2dB模拟带阻滤波器的设计指标参数有1p、1s、0、su、pu。p1和pu分别表示带阻滤波器的通带下边界频率和通带上边界频率；1s和su分别表示带阻滤波器的阻带下边界频率和阻带上边界频率。1和2分别称为通带和阻带波纹幅度 图1 带阻滤波器带阻滤波器的技术指标要求必须是几何对称的。如果带阻滤波器不是几何对称的，为了使对称的带阻滤波器满足最小阻带衰减要求，在由非对称的带阻滤波器变换为对称的带阻滤波器时，保留非对称带阻滤波器的阻带截止频率调整两个通带截止频率中的一个，调整过程步骤如下:（1） 计算(0)2=1ssu（2）计算 p1=(0)2/pu ，如果p11p,用p1代替1p；（3）如果p11p，计算pu=(0)2/pu，并用pu代替pu；（4） 如果A1pA2p，选择Ap=minA1p,A2p。 所有类型的变换，包括各自与低通原型的频率变换如表1所示。用频率变换法设计一个模拟滤波器的步骤如下：1、 确定低通、高通、带通和带阻模拟滤波器的技术要求。当要求的带通或带阻滤波器是非几何对称时，根据上面描述的过程使这些技术要求成几何对称的；2、 2、根据表1中的变换关系，确定归一化低通滤波器的技术要求：通带截止频率为，阻带截止频率为s，阻带衰减为Ap(dB)，阻带衰减为As(dB)；3、根据p、s、Ap和As，用巴特沃斯设计归一化低通滤波器；4、根据表1最右栏的变换关系，得到要求的非归一化模拟滤波器。 表1模拟滤波器的频率变换2、带通滤波器的设计原理 设计巴特沃斯带通滤波器时，首先应根据参数要求设计出相应的模拟滤波器，其步骤如下：1）由模拟滤波器的设计指标wp，ws，Ap，As和式（1）确定滤波器阶数N。 （1）2）由式（2）确定wc。 （2）3）若N是奇数，则有式（3）和式（4）确定滤波器的系统函数；若N是偶数，则式（3）和式（5）确定滤波器的系统函数。然后把设计好的模拟Butterworth滤波器变换成数字滤波器，常采用的方法有脉冲响应不变法和双线性变换法。 （3） （4） （5）巴特沃什滤波器是根据幅频特性在通频内具有最平坦特性而定义的滤波器。一维巴特沃什滤波器的平方幅频特性函数为 （6）可以看出，滤波器的幅频特性随着滤波器阶次N的增加而变得越来越好。在截止频率c处的函数值始终为的情况下，在通带内更多的频带区的值接近1；在阻带内函数更迅速的趋近0。巴特沃什滤波器的主要特征：（1）对于所有N，（2）对于所有N，（3）是的单调下降函数。（4）随着阶次N的增大而更加接近于理想滤波器。二、软件实现1、带阻滤波器设计带阻滤波器指标： 阻带上边界频率：5Kz； 阻带下边界频率：7Kz； 通带上边界频率：2Kz； 通带下边界频率：9Kz； 通带最大衰减：1dB； 阻带最小衰减：20dB；设计程序如下：wp=2*pi*2000,9000;ws=2*pi*5000,7000;Rp=1;As=20;Nb,wc=buttord(wp,ws,Rp,As,s);BSB,ASB=butter(Nb,wc,stop,s);hk,w=freqs(BSB,ASB);subplot(2,2,1)plot(w,20*log(abs(hk),b)xlabel(频率/Hz)ylabel(幅度/dB)title(幅频特性)subplot(2,2,2)plot(angle(hk)xlabel(频率/Hz)ylabel(相位)title(相频特性)程序运行结果：巴特沃斯模拟带阻滤波器阶数：Nb=4巴特沃斯模拟带阻滤波器系统函数分子多项式系数向量：巴特沃斯模拟带阻滤波器系统函数分母多项式系数向量： 图2 带阻滤波器的幅频特性和相频特性2、 带通滤波器的设计 带通滤波器指标： 通带上边界频率：4Kz； 通带下边界频率：7Kz； 阻带上边界频率：2Kz； 阻带下边界频率：9Kz； 通带最大衰减：1dB； 阻带最小衰减：20dB；设计程序如下：wp=2*pi*4000,7000;ws=2*pi*2000,9000;Rp=1;As=20;N,wc=buttord(wp,ws,Rp,As,s);BB,AB=butter(N,wc,s);hk,w=freqs(BB,AB);subplot(2,2,1)plot(w,abs(hk),b)xlabel(频率/Hz)ylabel(幅度/dB)title(幅频特性)subplot(2,2,2)plot(angle(hk)xlabel(频率/Hz)ylabel(相位)title(相频特性)程序运行结果：巴特沃斯模拟带通滤波器阶数：N=5巴特沃斯模拟带通滤波器系统函数分子多项式系数向量：巴特沃斯模拟带通滤波器系统函数分母多项式系数向量： 图4 带通滤波器的幅频特性和相频特性 三、仿真设计1、 带阻滤波器仿真设计 利用matlab中的simulink来进行带阻滤波器的仿真 图4 带阻滤波器的仿真设计图 图5 带阻滤波器参数设计下图所示为示波器scope输出波形 图6 输入带阻滤波器的波形下图所示为示波器scope1输出波形，经过带阻滤波器的滤波后只输出不在阻带范围内的波 图7 带阻滤波器输出波形2、 带通滤波器的仿真设计利用matlab中的simulink来进行带通滤波器的仿真 图8 带通滤波器的仿真设计图 图9 带通滤波器参数设计下图所示为示波器scope输出波形 图10 输入带通滤波器的波形下图所示为示波器scope1输出波形，经过带通滤波器的滤波后只输出在通带范围内的波 图11 带通滤波器输出波形致谢在做计算机课设的这段时间内，通过查找各种资料和自己的努力终于通过软件将带通和带阻做出来了。能够把其做出来除了自身努力外，与各位老师、同学和朋友的关心、支持和鼓励是分不开的。非常感谢陈常婷老师在这段时间给自己的指导，从最初的定题，到资料收集，到写作、修改，到定稿，她给了我