[精品论文]毕业论文__基于matlab的数字滤波器的设计

基于MATLAB的数字滤波器的设计摘要数字滤波器是数字信号处理的重要环节，数字滤波器可分为IIR和FIR两大类。本文介绍了IIR和FIR的基本设计原理以及在MATLAB环境下如何利用直接程序设计法、SPTOOL设计法和FDATOOL设计法给出IIR和FIR数字滤波器的设计方法和操作步骤，并给出设计设计实例及运行结果，同时利用MATLAB环境下的仿真软件SIMULINK对所设计的滤波器进行模拟仿真，仿真结果表示设计参数设置合理。关键词数字滤波器，仿真，MATLAB，FIR，IIR。1引言MATLAB是“矩阵实验室”（MATRIXLABORATOY）的缩写，它是由美国MATHWORKS公司于1984年正式推出的，是一种以矩阵运算为基础的交互式程序语言，专门针对科学、工程计算及绘图的需求。MATLAB是功能强大的科学及工程计算软件，它不但具有以矩阵计算为基础的强大数学计算和分析功能，而且还具有丰富的可视化图形表现功能和方便的程序计算能力。MATLAB的应用领域极为广泛，除数学计算和分析外，还被广泛地应用于自动控制、系统仿真、数字信号处理、图形图像分析、数理统计、人工智能、虚拟现实技术、通信工程、金融系统等领域，因此，MATLAB是面向21世纪的计算机程序设计及科学计算语言。数字滤波器是指完成信号滤波处理功能的，用有限精度算法实现的离散线性非时变系统，其输入是一组（由模拟信号取样和量化的数字量，其输出是经过变换或说处理的另一组数字量。因此，它本身既可以是用数字硬件装配成德一台完成给定运算的专用数字计算机，也可将所需的运算编成程序，让通用计算机来执行。数字滤波器具有稳定性高、精度高、灵活性大等突出优点。随着数字技术的发展，用数字技术实现滤波器的功能愈来愈受到人们的注意和广泛使用。这里所说的数字滤波器是指理想带通，低通等的频率选择数字滤波器。数字滤波器设计的一个重要步骤是确定一个可实现的传输函数HZ，这个确定传输函数HZ的过程称为数字滤波器设计。数字滤波器的一般设计过程为（1）按照实际需要，确定滤波器的性能要求（通常在频域内给定数字滤波的性能要求）。（2）寻找一满足预定性能要求的离散时间线性系统。（3）用有限精度的运算实现所设计的系统。（4）通过模拟，验证所设计的系统是否符合给定性能要求。2数字滤波器的设计滤波器分为两种，分别为模拟滤波器和数字滤波器。数字滤波器可以理解为是一个计算程序或算法，将代表输入信号的数字时间序列转化为代表输出信号的数字时间序列，并在转化的过程中，使信号按预定的形式变化。数字滤波器有多种分类，从数字滤波器功能上分可分为低通、高通、带阻、带通滤波器，根据数字滤波器冲激响应的时域特征，可将数字滤波器分为两种，即无限长冲激响应滤波器（IIR）和有限长冲激响应滤波器（FIR）。数字滤波器指标一般来说，滤波器的幅频特性是分段常数的，以低通为例，在通带内逼近于1，阻带内逼近与0，实际设计的滤波器并非是锐截止的通带和阻带两个范围，两者之间总有一个过渡带。在设计滤波器时事先给定幅频特性允许误差，在通带范围内幅度响应以误差逼近于1，在阻带内幅度响应以误差逼近于0。2（1）WEHRJWCJ,2|1|1式中WC和WR分别为通带边界频率和阻带边界频率，WRWC为过渡带。在具体的技术指标中往往用通带波动来表示，用最小阻带衰减AT来表示，其具体的对应公式这里就不详述了。121IIR数字滤波器设计IIRDF的冲激响应HN是无限长的，其输入输出的关系为2INXHNYI系统函数为3NNZHZHKKMRRZAB10设计无限长单位脉冲响应（IIR）数字滤波器一般可有三种方法。第一种方法，先设计一个合适的模拟滤波器，然后将其数字话，即将S平面映射到Z平面得到所需的数字滤波器。模拟滤波器的设计技巧非常成熟，不仅得到的是闭合形式的公式，而且设计系数已经表格化了。因此，由模拟滤波器设计数字滤波器的方法准确，简便，得到普遍采用。对于这种方法，工程上有两种常见得变换法脉冲响应不变法及双线性变换法。第二种方法，在Z平面直接设计IIR数字滤波器，给出闭合形式的公式，或者以所希望的滤波器响应作为依据，直接在Z平面上通过多次选定极点和零点的位置，以逼近该响应。第三种方法，利用最优化技术设计参数，选定极点和零点在Z平面上的合适位置，在某种最优化准则意义上逼近所希望的响应。但一般不能得到滤波器的系数（即零，极点的位置）作为给定响应的闭合形式函数表达式。优化设计需要完成大量的迭代运算，这种设计法实际上也是IIR滤波器的直接设计。本文着重介绍由模拟滤波器设计相应的IIR数字滤波器的方法。（1）脉冲响应不变法脉冲响应不变法是使数字滤波器的单位脉冲响应序列HN逼近模拟滤波器的冲激响应，让AHTHN正好等于的采样值。设已有满足要求的模拟滤波器,则可AHTAHSAHTANHTHZ因为,T的图形的图形公式导出具体转换如下设以一阶极点为例41NKAAHSS作拉氏反变换,得11KNSTAAKHTFTHSAEU采样得1KNSNTAKHNEU作Z变换,得51011KKNNSNTNSTNKKAHZAEUZEZ与极点关系为6HZASKSTKZE一般对应关系STJTJJZREE7,TRE图1S平面到Z平面变换示意图所以,模拟系统稳定因果数字系统稳定因果。按照脉冲响应不变法，从S平面到Z平面的映射不是单值关系，而是先将在S平面沿虚轴AHS作周期严拓，再按照映射关系将映射到Z平面，得到，因此，脉冲响应不变法只适用于带AHSHZ限的滤波器（如低通、带通）。在MATLAB中利用M文件IMPINVAR可以对模拟传输函数实行脉冲响应不变法。（2）双线性变换法脉冲响应不变法不适带阻和高通滤波器的设计，因为高频带为通带，前述方法易引起混频。故希望平面虚轴Z平面单位圆一周,且应有S,00,因为,所以选变换TAN/28TAN2K其中可取任意正常数,后面将导出K/TJO1ROS平面0Z平面TREKSKZ设计思路设计出模拟滤波器转化成数字滤波器TAN2,PSSPSPSK图2数字域频率与模拟域频率的对应关系转化公式推导如下SIN/2TAN2COJJJK/2/2/1JJJJJJEEKK因只关心频率转换,故可设,则有SJJZE,称为双线性变换91K所以模拟滤波器转换成数字滤波器的公式为101ZASKHZ由双线性变换公式,可得,视为两复平面变换,再由1,JSJZRE可得11221KR从上式可得时,上虚轴上单位圆周。01RSZ时,上左半平面上单位圆内。时,上右半平面上单位圆外。故若模拟滤波器稳定，则双线性变换后数字滤波器也稳定。由于双线性变换法是一种单值映射，因此消除了频率混叠的现象。双线性变换法的缺点是模拟频率与数字频率间的非线性，这种非线性关系要求被变换的连续系统的幅度响应是分段常数型的（某一频率范围内幅度响应近似于某一常数），不然所映射出的数字频率响应相对于原来的模拟频率响应会产生变形。为解决双线性变换中的频率非线性关系，我们采用预畸的方法，即，其中K2/T。TAN2K在MATLAB中利用M文件BILINEAR可以对模拟传输函数实行双线性变换法。1050510202W2ATANWPPOMATLAB中IIR数字滤波器的设计过程包括两步。第一步，根据给定指标，确定滤波器的阶数N和频率缩放因子WN。第二步，利用这些参数和给定的波纹参数，确定传输函数的关系。阶数估计利用双线性变换法设计数字滤波器时，首先要对IIR数字滤波器的阶数进行估计，相应的M文件为BUTTORD用于巴特沃斯滤波器，CHEB1ORD用于切比雪夫1型滤波器，CHEB2ORD用于切比雪夫2型滤波器，ELLIPORD用于椭圆滤波器。滤波器的设计对于基于双线性变换法的IIR滤波器设计，对应于四种逼近技术（即巴特沃斯、切比雪夫1型和2型及椭圆逼近），MATLAB工具箱中有相应的函数。特别地可以用到下面的M文件BUTTER用于巴特沃斯滤波器的设计，CHEBY1用于切比雪夫1型滤波器的设计，CHEBY2用于切比雪夫2型滤波器的设计，ELLIP用于椭圆滤波器的设计。这些函数的输出可以是滤波器传输函数分子和分母的系数向量，也可以是滤波器的零极点向量和标量增益因子。同时，利用ZP2TF可以由滤波器的零极点向量和标量增益因子得到传输函数分子和分母的系数向量。相应地，利用函数ZP2SOS可以得到传输函数分子和分母系数向量的二次项因子。在计算出传输函数的系数之后，可以利用M文件FREQZ来计算频率响应。22FIR数字滤波器设计FIRDF的冲激响应HN是有限长的，M阶FIRDF可以表示为1210MIINXHNY其系统函数为1310MNNZHZH与IIR数字滤波器的设计不同，FIR滤波器的设计与模拟滤波器的设计没有任何联系。因此，FIR滤波器的设计基于对指定幅度响应的直接逼近，并通常要求其具有线性相位响应。为了保证滤波器具有线性相位特性，滤波器系数必须满足条件HNHM1N。目前关于FIR滤波器的设计方法主要有三种，即窗函数法，频率取样法和切比雪夫等波纹逼近的最优化设计方法。一般应用较多的是第一种和第三种方法。这是因为窗函数法比较简单，可应用现成的窗函数公式，在技术指标要求不严格的情况下市比较灵活的。最优化设计法必须借助计算机计算，但是它能得到最佳的等波纹的线性相位FIR滤波器。目前切比雪夫等波纹的线性相位FIR滤波器的计算机机助设计程序已经比较完善，由于采用了REMEZ迭代算法，所以设计效率也很高，在应用中越来越占优势。（1）窗函数法一般设计过程总是先给定一理想的滤波器频率响应，然后设计一个FIR滤波器，用它的JWDEH频率响应来逼近理想的。这种逼近中最直接的方法，是在时域中用JWNMNJWEHEH10JDFIR滤波器的单位脉冲响应HN去逼近理想的单位脉冲响应。因而，先由的IDTFT导出NHJWDEHNHD14WEHJNDD21由于是矩形频率特性，故一定是无限长的序列，且是非因果的。然而FIR滤波器是有限JWDEH长的，所以用有限长的HN来逼近无限长的，最简单的方法是截取中最重要的一段，将HDNHD无限长的截取成长度为M的有限长序列，等效于再上施加了一个长度为M的矩形窗口，更NHDNHD为一般的，可以用一个长度为M的窗口函数WN来截取，即15NHWD这一方法通常称为窗函数法，窗口函数的形状及长度M的选择是窗函数法的关键。下面我们一低通为例，了解一下窗函数法的运用提出希望频率响应函数（低通）图3理想低通滤波器的频响线性相位,具有片断特点,即|0|JJCDCEH算出1D2JJNDHNECJJN无限长SINC图4理想低通的单位脉冲响应（无限长的一部分）加窗,长,得WNNDHNWN要线性相位,就要关于偶对称,而关于偶对称,故要求H1/21/2NO025|JDHE1010203000510102030010010203所以要求关于偶对称WN1/2N010203000510102030010010203图5窗函数图6加窗后的单位脉冲响应再回过来检验是否满足精度要求JHE图7图4的脉冲响应的频响图8理想频响与实际频响的对比若基本满足,则依截取的,制硬件,编软件HN为便于选择使用,将5种常见的窗函数基本参数如表1所示。表15种常见的窗函数基本参数类型窗函数的旁瓣峰N过渡带宽度B加窗后滤波器的阻带最小衰减SRECTWIN134/N21BARTLET三角258/N25HANNING318/N44HAMMING418/N53BLACKMAN5712/N74（2）频率取样法窗口设计法事从时域出发，把理想的用一定形状的窗口函数截取成有限长的HN，以此HN来NHD近似理想的，从而频率响应也近似于理想的频率响应。我们知道一个有限长NHDJWEHJWDEH序列可以通过其频谱的相同长度的等间隔采样值准确地恢复原有的序列。频率采样法便是从频域出发，对理想的频率响应加以等间隔采样JWE16|2KDKMWJD然后以此作为实际FIR滤波器的频率特性的离散样本HK，即KHD171,20,|2MEHKWJDD01230051O2|JDHEJE01230051由HK通过IDFT可求出有限长序列HN为181,0,102MNEKHMNHKJ利用M个频率的离散样本HK同样可求出FIR滤波器的系统函数HZ及频率响应。JWEH19MJMKKNNEWZHZHZ210110，其中令可得到滤波器的频率响应。如果设计的是线性相位的FIR数字滤波器，其采样值HKJWEZJWE的相位的幅度一定要满足特定的约束条件，这个设计时一定要注意。（3）最优化设计法最优化设计法事以最佳一致逼近（最大误差最小化）理论为基础，利用雷米兹算法设计的具有等波纹特性的设计方法。具体设计步骤如下对设计指标进行归一化处理。确定REMEZORD函数所需要的参数。包括归一化边界频率、各频带的幅度要求和波纹要求等。归一化边界频率总是从0开始到1结束，故只需递增列出中间的边界频率；频带幅度要求不含过渡区，个数是边界频率个数的一半加1；波纹要求是频带内幅度允许的波动要求，与分贝间的关系是20LOG201LOG20,1LOG2021SPRR利用REMEZORD函数确定REMEZ所需参数。调用REMEZ函数进行设计。利用FREQZ函数验算技术指标是否满足要求。23数字滤波器类型的选择IIR和FIR各有优缺点，在实际运用中如何选择它们，这里做一个简单的比较。表2IIR与FIR的比较IIRFIR设计方法利用AF的设计图表，可简单，有效的完成设计一般无解析的设计公式，要借助计算机程序完成设计结果只能得到幅频特性，相频特性未知（缺点），如需要线性相位，需用全通网络校准，但增加滤波器的阶数和复杂性可得到幅频特性（可以多带）和线性相位（优点）稳定性有稳定性问题极点全部在原点（永远稳定），无稳定性问题因果性总是满足，任何一个非因果的有限长序列，总可以通过一定的延时，转变为因果序列结构递归系统非递归运算误差有反馈，由于运算中的四舍五一般无反馈，运算误差小入会产生极限环快速算法无快速运算方法可用FFT减少运算量从以上简单的比较可以得到，IIR与FIR滤波器各有所长，所以应根据实际应用要求，从多方面考虑加以选择。3数字滤波器的MATLAB设计31直接程序设计法（1）IIR的直接程序设计法例如欲设计一数字（IIR）带阻滤波器，其数字域指标为数字阻带边缘频率分别为04和07，数字通带边缘频率为025和08，通带波动为1DB最小阻带衰减为40DB。此题的MATLAB程序为WS04PI07PI数字阻带边缘频率WP025PI08PI数字通带边缘频率RP1通带波动（DB）AS40阻带衰减N,WNCHEB2ORDWP/PI,WS/PI,RP,AS根据给定指标，确定滤波器的阶数N和频率缩放因子WNB,ACHEBY2N,AS,WS/PI,STOP返回的B,A分别为H（Z的分子、分母。H,WFREQZB,A,512返回的H,W分别为滤波器的频率响应及其频率PLOTW/PI,ABSH画出频率响应（以W/PI为横轴）GRIDXLABELW/PIYLABEL幅值TITLE频率响应程序运行结果为图9所设计的带阻滤波器的频率响应在设计中如果该滤波器的特性不满足要求，原有的参数必须做相应的调整，在程序中只需对参数做新的设定就可以得到所需的滤波器。接下来我们来看看此题所设计的滤波器的滤波效果S为含有3个频率成分的信号（归一化频率（W/2）分别为01、03、045），用所设计的滤波器滤除归一化频率为03的成分。N0100S1SINPI02NS2SINPI06NS3SINPI09NS4S1S3SS1S2S3SFFILTERB,A,SSUBPLOT311STEMN,STITLE滤波前的信号SUBPLOT312STEMN,SFTITLE滤波后的信号SUBPLOT313STEMN,S4TITLE想要保留的信号程序运行的结果为图10采用FILTER函数进行数字滤波前后信号比较示意图由图可以看出，滤波后的信号与想要保留的信号基本一致（相位有些许偏差，但基本一致），所以我们可以说该滤波器基本满足了以上所提出的滤波要求。（2）FIR的直接程序设计法例如欲设计一个线性相位数字（FIR）带通滤波器，其数字域指标为数字通带边界频率为035和065，数字阻带边界频率为02和08，通带波动为1DB，最小阻带衰减为60DB。FIR数字滤波器的窗函数法此题的MATLAB程序为WS102PIWP1035PIWP2065PIWS208PIAS60TRMINWP1WS1,WS2WP2MCEIL11PI/TR1滤波器的阶数，程序运行后M75N01M1RM1/2R为群时延WC1WS1WP1/2WC2WP2WS2/2HDSINWC2NREPS/PINREPSSINWC1NREPS/PINREPSHD为理想滤波器的脉冲响应W_BLABLACKMANM长度为M的BLACKMAN窗HHDW_BLAH为滤波器的实际脉冲响应STEMN,HTITLE滤波器的实际单位脉冲响应FREQZH,1,512TITLE幅度响应和相位响应图11所设计的滤波器的实际单位脉冲响应由上图可知滤波器的实际脉冲响应H是偶对称的，即HNHM1N,故该滤波器满足FIR线性相位的条件，该滤波器是线性相位的FIR滤波器。图12所设计的带通滤波器的幅度和相位响应由滤波器的相位特性也可以看出该滤波器是线性相位的FIR滤波器。接下来我们来看看此题所设计的滤波器的滤波效果S为含有3个频率成分的信号（归一化频率（W/2）分别为005、02、045），用所设计的滤波器滤除归一化频率为005和045的成分。L0100S1SIN01PILS2SIN04PILS3SINPI09LSS1S2S3SFFILTERH,1,SSUBPLOT311STEML,STITLE滤波前的信号SUBPLOT312STEML,SFTITLE滤波后的信号SUBPLOT313STEML,S2TITLE想要保留的信号图13采用FILTER函数进行数字滤波前后信号比较示意图由上图可知滤波后的信号和想要保留的信号的幅度和频率基本一致（滤波后的信号相对于想要保的信号有一个相位延迟，这是线性相位FIR滤波器的群延迟引起的，此滤波器留的群延迟RM1/237），所以我们可以说该滤波器基本满足了以上所提出的滤波要求。FIR数字滤波器的频率采样法此题的MATLAB程序为M40取滤波器的阶数为40ALM1/2群时延N0M1T2059417456T10109021HRSZEROS1,5,T1,T2,ONES1,7,T2,T1,ZEROS1,9,T1,T2,ONES1,7,T2,T1,ZEROS1,4采样值的幅值K10FLOORM1/2K2FLOORM1/21M1ANGHAL2PI/MK1,AL2PI/MMK2采样值的相位HHRSEXPJANGHHREALIFFTH,M长度为M的单位脉冲响应STEMN,HTITLE滤波器的实际单位脉冲响应FREQZH,1,512TITLE幅度响应和相位响应图14所设计的滤波器的实际单位脉冲响应由图14可知滤波器的实际脉冲响应H是偶对称的，即HNHM1N,故该滤波器满足FIR线性相位的条件，该滤波器是线性相位的FIR滤波器。图15所设计的带通滤波器的幅度和相位响应由滤波器的相位特性也可以看出该滤波器是线性相位的FIR滤波器。此滤波器的群延时为ALM1/2195。接下来我们来看看此题所设计的滤波器的滤波效果S为含有3个频率成分的信号（归一化频率（W/2）分别为005、02、045），用所设计的滤波器滤除归一化频率为005和045的成分。L0100S1SIN01PILS2SIN04PILS3SINPI09LSS1S2S3SFFILTERH,1,SSUBPLOT311STEML,STITLE滤波前的信号SUBPLOT312STEML,SFTITLE滤波后的信号SUBPLOT313STEML,S2TITLE想要保留的信号图16采用FILTER函数进行数字滤波前后信号比较示意图同上面分析相似，滤波后的信号和想要保留的信号的幅度和频率基本一致（滤波后的信号相对于想要保的信号有一个相位延迟，这是线性相位FIR滤波器的群延迟引起的，此滤波器留的群延迟（RM1/2195），所以我们可以说该滤波器基本满足了以上所提出的滤波要求。FIR数字滤波器的最优设计法此题的MATLAB程序为设计指标WS102PIWP1035PIWP2065PIWS208PIRP1AS60设置边界频率和幅度要求FWS1/PI,WP1/PI,WP1/PI,WS2/PIA0,1,0设置各频带的波纹要求DEVP10RP/201/10RP/201DEVS10AS/20DEVDEVS,DEVP,DEVS确定REMEZ参数，其中滤波器的阶数为（N1），程序运行后得到N26N,FO,AO,WREMEZORDF,A,DEV调用REMEZ函数进行设计HREMEZN,FO,AO,WN0NSTEMN,HTITLE滤波器的单位冲激响应FREQZH,1,512TITLE幅度响应和相位响应图17所设计的滤波器的实际单位脉冲响应由图17可知滤波器的实际脉冲响应H是偶对称的，即HNHNN,故该滤波器满足FIR线性相位的条件，该滤波器是线性相位的FIR滤波器。图18所设计的带通滤波器的幅度和相位响应由滤波器的相位特性也可以看出该滤波器是线性相位的FIR滤波器。此滤波器的群延时为ALN/213。接下来我们来看看此题所设计的滤波器的滤波效果S为含有3个频率成分的信号（归一化频率（W/2）分别为005、02、045），用所设计的滤波器滤除归一化频率为005和045的成分。L0100S1SIN01PILS2SIN04PILS3SINPI09LSS1S2S3SFFILTERH,1,SSUBPLOT311STEML,STITLE滤波前的信号SUBPLOT312STEML,SFTITLE滤波后的信号SUBPLOT313STEML,S2TITLE想要保留的信号图19采用FILTER函数进行数字滤波前后信号比较示意图同上面分析相似，滤波后的信号和想要保留的信号的幅度和频率基本一致（滤波后的信号相对于想要保的信号有一个相位延迟，这是线性相位FIR滤波器的群延迟引起的，此滤波器留的群延迟RN/213），所以我们可以说该滤波器基本满足了以上所提出的滤波要求。32利用信号处理工具箱SPTOOL设计法SPTOOL是信号处理工具箱中一个具有交互式图形用户界面的信号处理工具，专门用于完成常用的数字信号处理任务。通过这个工具，只需要鼠标简单的操纵鼠标（点击或拖动），就可以完成载入、观察、分析、实现和设计数字滤波器并进行谱分析等甚至十分复杂的数字信号处理任务，而不需要用户对数字滤波器的设计原理非常熟悉。下面是用SPTOOL工具设计的IIR滤波器对上述信号进行滤波。SPTOOL使用步骤导入信号、滤波器、频谱，设计滤波器，对信号滤波，分析输入输出信号的谱。1IIR滤波器的设计创建并导入信号源，在MATLAB命令窗口输入命令FS4000（这里取采样频率为4000HZ）N0100S1SINPI02NS2SINPI06NS3SINPI09NS4S1S3SS1S2S3此时，变量FS、N、S、S1、S2、S3、S4将显示在MATLAB的WORKSPACE列表中。在MTLAB命令窗口中输入SPTOOL，即可弹出SPTOOL的主界面,如图所示图20SPTOOL的主界面这里我们首先来设计前面的IIR带阻滤波器，前面的数字指标转化为模拟指标为通带频率分别为500HZ和1600HZ，阻带频率为800HZ和1400HZ，通带波动为1DB，最小阻带衰减为40DB。以上三个信号转换为模拟域的频率分别为400HZ，1200HZ，1800HZ。1）导入信号。使用【FILE/IMPORT】可以导入信号、滤波器和谱。信号的来源可以是MATLAB工作空间变量或MAT数据文件，也可以是在DATA前面的文本框中直接输入信号的数据。这里信号源仍然利用前面程序设计法中的混频信号，将信号S和采样频率FS4KHZ导入并命名为SIG1，将信号S2导入命名为SIG2,将信号S4导入命名为SIG4。2）滤波器的设计。单击FILTER栏中的NEW按钮打开滤波器设计工具。在界面的最上面制定滤波器的名字、采样频率和设计方法。这里命名滤波器名称为FILT1，采样频率FS4KHZ和设计方法为CHEBYSHEVTYPEIIIIR；左侧的SPECIFICATION面板用来指定设计指标；MINIMUMORDER,TYPE为BANDSTOP，通带频率分别为500HZ和1600HZ，阻带频率为800HZ和1400HZ，通带波动为1DB，最小阻带衰减为40DB。完成以上输入之后单击APPLY按钮，则完成了滤波器的设计；中间的FREQUENCYRESPONSE面板则根据设计指标自动绘出滤波器的幅频响应，而右侧的MEASUREMENTS面板显示了滤波器设计完成后的实测参数。此时在FILTERS栏中，选中FILT1按下VIEW按钮即弹出界面，我们还可以观察到滤波器的相位响应、幅频相位响应、脉冲响应、阶跃响应等等（只须点击界面上的相应快捷按钮）。图21所设计的滤波器的幅值响应3）将滤波器应用到SIG1信号序列。分别在SIGNALS、FILTERS栏中选择SIG2、FILT1单击FILTERS栏列表下的APPLY按钮，在弹出的APPLYFILTER对话框中将输出信号命名为SIG3（滤波后信号）。同时按下ALT与SHIFT键选中SIG3与SIG2点击SIGNALS栏下的VIEW可观察到它们的时域波形其中红色是信号SIG3，蓝色是信号SIG2，如图所示。由图可观察出信号S2通过该带阻滤波器后得到很大的衰减，几乎衰减为0。图22信号S2滤波前后的比较分别在SIGNALS、FILTERS栏中选择SIG1、FILT1单击FILTERS栏列表下的APPLY按钮，在弹出的APPLYFILTER对话框中将输出信号命名为SIG5（滤波后信号）。同时按下ALT与SHIFT键选中SIG4与SIG5点击SIGNALS栏下的VIEW可观察到它们的时域波形其中红色是信号SIG5，蓝色是信号SIG4，如图所示。由图可观察出信号S通过滤波器后所得的信号与S4基本吻合，只是相位上有点偏差，即该带阻滤波器基本滤除了信号S2，滤波效果良好。图23滤波后的信号与S4的比较4）进行频谱分析。在SIGNALS中选择SIG1，单击SPECTRA栏下的CREATE按钮，在弹出的SPECTRAVIEWER界面中选择METHOD为FFT，NFFT512，单击APPLY按钮生成S的频谱SPECT1同理得到SIG4的频谱SPECT4，SIG5的频谱SPECT5,同时按下ALT与SHIFT键选中SPECT1与SPECT5，点击SPECTRA栏下的VIEW可观察到它们的频谱,下图是SIG1和SIG5的频谱，其中红色代表滤波后信号的频谱（SIG5），蓝色代表滤波前信号的频谱（SIG1）。图24滤波前后信号的频谱的比较由频谱可以看出，混频信号S中频率为1200HZ的部分基本被滤除了，频率为400HZ和1800HZ的部分基本被保留下来，未受滤波的影响。同时按下ALT与SHIFT键选中SPECT4与SPECT5，点击SPECTRA栏下的VIEW可观察到它们的频谱,下图是SIG4和SIG5的频谱，其中红色代表滤波后信号的频谱（SIG5），蓝色代表滤波前信号的频谱（SIG4），如下图。图25滤波后的信号与信号S4频谱的比较由上图可以看出，S4与滤波后的信号频谱基本一致，所以我们可以得出结论该带阻滤波器的滤波效果较好。2FIR滤波器的设计与1的设计过程类似，这里就简述过程。创建并导入信号源，在MATLAB命令窗口输入命令FS4000HZL0100S1SIN01PILS2SIN04PILS3SINPI09LSS1S2S3S4S1S3前面的数字指标转化为模拟指标为通带频率分别为700HZ和1300HZ，阻带频率为400HZ和1600HZ，通带波动为1DB，最小阻带衰减为60DB，以上三个信号转换为模拟域的频率分别为200HZ，800HZ，1800HZ。1）导入信号。将信号S和采样频率FS4KHZ导入并命名为SIG1，将信号S2导入命名为SIG2,将信号S4导入命名为SIG4。2）滤波器的设计。图26所设计的滤波器的幅值响应3）将滤波器应用到SIG1信号、SIG4信号，生成的信号分别命名为SIG5和SIG3。图27信号S4滤波前后的比较上图中SIG4信号与SIG3信号比较（其中红色代表SIG3，蓝色代表SIG4）可知信号S4通过该带通滤波器后得到大大衰减，几乎衰减为0。图28滤波后的信号与信号S2的比较上图中SIG2和SIG5比较（其中红色代表SIG5，蓝色代表SIG2），可知信号S通过该带通滤波器后的信号与S2基本具有相同的频率和幅值（只是相位上有一点误差），所以我们认为该滤波器较好的滤掉了S4，滤波性能优良。4）进行频谱分析。分别生成SIG1的频谱SPECT1，SIG2的频谱SPECT2，SIG5的频谱SPECT5。下图是SIG1和SIG5的频谱，其中红色代表滤波后信号的频谱（SIG5），蓝色代表滤波前信号的频谱（SIG1）。由图可看出该滤波器基本滤除了频率为200HZ和1800HZ的信号，保留了频率为800HZ的信号。图29滤波前后信号的频谱的比较下图是SIG2和SIG5的频谱，其中红色代表滤波后信号的频谱（SIG5），蓝色代表滤波前信号的频谱（SIG2）。由图可看出滤波后的信号与S2具有基本一致的频谱。图30滤波后的信号与信号S2频谱的比较33滤波器分析设计工具箱FDATOOL法FDATOOL（FILTERDESIGNS1SINPI02NS2SINPI06NS3SINPI09NS4S1S3SS1S2S3修改正玄信号的参数SINETYPE/TIMEBASED、FREQUENCY分别改为02PI，06PI，09PI，SAMPLETIME改为1。其他模块的参数根据上述要求进行调整。将下图分别是仿真框图和滤波前后信号的波形图34仿真框图图35滤波前后信号的比较由上图可以看出信号S经过滤波后大致分离出了归一化频率为03的信号，即滤波后的信号与S4大致相同（与前面两种方法的