窗函数FIR滤波器设计方案与SIUMLINK仿真

1、个人资料整理仅限学习使用电子技术综合设计姓 名：学号：专 业： 电子信息工程题 目：基于MATLAB勺FIR滤波器的设计与仿真课程名称： 专业创新设计设计地点：电工电子实验中心设计日期：成绩：指导教师:年月电子技术综合设计任务书设计日期：设计专题：电子技术综合设计设计题目：基于MATLAB勺FIR滤波器的设计与仿真设计内容和要求:1 .主要内容：了解FIR滤波器的原理。了解使用MATLA酷言设计FIR滤波器的方法。了解使用MATLAB勺仿真软件SIUMLINK的仿真方法。2 .设计要求：用窗函数法设计一个数字低通滤波器。通带截止频率：0.2兀，阻带起始频率：0.4兀，阻带最小衰减：-50dB。

2、采用窗函数法设计，并且采用不同的窗做比较。使用FDATOOLS计滤波器并分析。 使用SIMULINKM设计的滤波器进行仿真。b5E2RGbCAP3 .设计滤波器的MATLAB程序汉宁窗，哈明窗，矩形窗，三角窗设计滤波器的程序计算机打印）4 .程序运行结果的图汉宁窗，哈明窗，矩形窗，三角窗程序运行结果图计算机打印）5 .仿真模块原理图使用SIUMLINK设计仿真模块原理图 计算机打印）指导教师年月日摘要传统的FIR滤波器的设计过程复杂，计算工作量大，滤波特性调整困难，影响了它的应用。本次设计利用MATLAB言号处理工具箱快速有效的设计由软件组成的常规FIR滤波器。给出了使用MATLAB言进行程序

3、设计和利用信号处理工具箱的 FDATool工具进行滤波器设计的步骤。利用MATLA豉计滤波器，可以随时对比设计要求和滤波器特性调整参数，直观简 便，极大的减轻了工作量，有利于滤波器设计的最优化。同时 还利用MATLA酬境下的仿真软件Simulink对所设计的滤波器 进行模拟仿真。plEanqFDPw【关键字】FIR滤波器MATLAB FIR IIR使用目录第一章 绪论6DXDiTa9E3d1.1 FIR滤波器的设计背景 61.2 FIR滤波器的优缺点 6第二章FIR数字滤波器的介绍 62.1 FIR数字滤波器的特点 62.1.1 线性相位FIR数字滤波器的特点 72.2 FIR 数字滤波器的设

4、计原理 82.3 数字滤波器的性能指标 9第三章窗函数设计法 103.1 .窗函数设计原理分析 113.2 窗函数的一些通用性质 113.3 设计方法113.4 窗函数介绍 143.5 窗函数法设计步骤 17第四章设计实例 184.1 滤波器设计指标要求 184.2 设计步骤 184.3 仿真结果 194.4 仿真结果分析 21第五章 基于siumLInk的FIR滤波器的仿真 22RTCrpUDGiT5.1 使用FDATOOL计滤波器 225.2 滤波器的分析 245.3 设计仿真模块 26使用5.4 设置输入信号，并导入滤波器 275.5对设计好的滤波器进行仿真与分析 27第六章总结 29参

5、考文献 29附录一 32第一章绪论1.1 FIR滤波器的设计背景随着电子工业的发展，对滤波器的性能要求越来越高，功能也越来越多， 并且要求它们向集成方向发展。我国滤波器研制和生产与上述要求相差甚 远，为缩短这个差距，电子工程和当代年轻人负有重大的责任。5PCzVD7HxA根据数字滤波器冲激响应函数的时域特性。可将数字滤波器分为两种，即无限长冲激响应（IIR 滤波器和有限长冲激响应（FIR滤波器。IIR滤波 器的特征是具有无限持续时间的冲激响应。 FIR滤波器冲激响应只能延 续一定时间。其中FIR滤波器很容易实现严格的线性相位，使信号经过处 理后不产生相位失真，舍入误差小，稳定等优点。能够设计具

6、有优良特性的 多带通滤波器、微分器和希尔伯特变换器，所以在数字系统、多媒体系统中获得极其广泛的应用。FIR数字滤波器的设计方法有多种，如窗函数设 计法、最优化设计和频率取样法等等。而随着MATLA瞅件尤其是MATLAB的信号处理工具箱和 Simulink仿真工具的不断完善，不仅数字滤波器的 计算机辅助设计有了可能，而且还可以使设计达到最优化。jLBHrnAILg1.2 FIR滤波器的优缺点FIR滤波器具有以下优点：1、具有严格的线性相位，同时又具有任意的幅度特性，因此满足要求信 道具有线性相位特性的图像处理以及数据传输。XHAQX74J0X2 、FIR滤波器的单位抽样响应是有限长的，因而滤波器

7、一定是稳定 的。3 、FIR滤波器由于单位冲激响应是有限长的，因而可以用快速傅立 叶变换（FFT算法来实现过滤信号，从而可大大提高运算效率。 LDAYtRyKfE使用FIR 滤波器的缺点：对于相同的设计指标，FIR 滤波器所要求的阶数比IIR滤波器要高，信号的延迟偏大。Zzz6ZB2Ltk第二章FIR数字滤波器的介绍2.1 FIR数字滤波器的特点有限长单位冲激响应(FIR数字滤波器可以做成具有严格的线性相位，同时又可以具有任意的幅度特性。此外， FIR滤波器的单位抽样响应 是有限长的，因而滤波器一定是稳定的。再有，只要经过一定的延时，任 何非因果有限长序列都能变成因果的有限长序列，因而总能用因

8、果系统来 实现。最后，FIR滤波器由于单位冲激响应是有限长的，可以用快速傅立 叶变换(FFT算法来实现过滤信号，从而可大大提高运算效率。但是，要 取得很好的衰减特性，FIR滤波器H(z的阶次比IIR滤波器的要高。dvzfvkwMII2.1.1 线性相位FIR数字滤波器的特点一、单位冲激响应h(n的特点FIR滤波器的单位冲激响应h(n是有限长(0 & n& N-1,其Z变换为：N 1H (z) =、'、h(n)z4mm3式 2-1)在有限Z平面有(N-1个零点，而它的(N-1个极点均位于原点z=0处。二、线性相位的条件如果FIR滤波器的单位抽样响应h(n为实数而且满足以下

9、任一条件：偶对称：h(n=h(N-1-n奇对称：h(n=-h(N-1-n其对称中心在n=(N-1/2处，则滤波器具有准确的线性相位。三、线性相位特点和幅度函数的特点四种线性相位FIR滤波器的特性可以总结如下：第一种情况，偶对称、奇数点，四种滤波器都可设计；第二种情况，偶对称、偶数点，可设计低、带通滤波器，不能设计高通和 带阻；第三种情况，奇对称、奇数点，只能设计带通滤波器，其它滤波器都不能 设计；第四种情况，奇对称、偶数点，可设计高、带通滤波器，不能设计低通和 带阻。2.2 FIR 数字滤波器的设计原理一个截止频率为0c (rad/s>的理想数字低通滤波器，其传递函数的表 达式是：Hd(

10、e=<式 <2-2>由式2-2可以看出，这个滤波器在物理上是不可实现的，因为冲激响 应具有无限性和因果性。为了产生有限长度的冲激响应函数，我们取样响应为h(n),长度为N,其系数函数为H (z) ： rqyni4ZNXIN 1式 <2-3>H (z)八 h(n)z”n 0用h(n)表示截取hd(n)后冲激响应，即h(n)=®(n)hd(n),式子中o(n) 为窗函数，长度为 No当r=(N-1>/2时，截取的一段h(n)对(N-1>/2对称，可保证所设计的滤波器具有线性相位。EmxvxOtOco一般来说，FIR数字滤波器输出y(n)的Z变换形

11、式Y(z>与输入x(n) 的Z变换形式之间的关系如下：_ _ _1. n _ _Y(z) =H(z)X(z) =(h(0)+h(1)z + +h(n)z )X(z)式<2-4>从上面的Z变换和结构图可以很容易得出FIR滤波器的差分方程表示形式。对式2-4进行反Z变换，可得：式 <2-5>y(n) =h(1)x(n) h(2)x(n -1)h(n)x(1)使用图2-1卷积型滤波器式2-5为FIR数字滤波器的时域表示方法，其中x(n)是在时间n的滤波器的输入抽样值。根据式 2-5即可对滤波器进行设计。从上面的公式 我们可以看出，在对滤波器实际设计时，整个过程的运算量很

12、大。SixE2yXPq52.3 数字滤波器的性能指标我们在进行滤波器设计时，需要确定其性能指标。一般来说，滤波器 的性能要求往往以频率响应的幅度特性的允许误差来表征。以低通滤波器特性为例，频率响应有通带、镇度带及阻带三个范围g在通带内：1- AP < H(e ) <10(e巧 Ja6ewMyirQFL 0< c在阻带中:其中At为阻带误差。st<CO6 c为通带截止频率，3 st为阻带截止频率，Ap为通带误差，图2-2低通滤波器的幅度特性与模拟滤波器类似，数字滤波器按频率特性划分为低通、高通、带 通、带阻、全通等类型，由于数字滤波器的频率响应是周期性的，周期为 2 九

13、o kavU42VRUs由于频率响应的周期性，频率变量以数字频率0来表示，所以数字滤波器设计中必须给出抽样频率。1、一个高通滤波器相当于一个全通滤波器减去一个低通滤波器。2、一个带通滤波器相当于两个低通滤波器相减。使用3 、一个带阻滤波器相当于一个低通滤波器加上一个高通滤波器。第三章窗函数设计法3.1.窗函数设计原理分析设数字滤波器的传输函数为 应，H (z)为系统函数。H (ej。), hd (n)是与其对应的单位脉冲响NH(ej8)=£ h(n)e,。式<3-1>n =01 二i.儿二，Hd(e )ej%。式 <3-2>N JH(z)=E h(n)z口式&

14、lt;3-3>n=0一般说来，hd(n)是无限长的，需要求对H d(ej巧的一个逼近。采用窗 函数设计法时，可通过对理想滤波器的单位采样响应加窗设计滤波器y6V3ALoS89h(n) =co(n)hd(n)式<3-4>其中，®(n)是一个长度有限的窗，在区间0 关于中间点对称< n < N外值为0 ,且6(n)=0(N 1 n)式 <3-5>频率响应根据式3-5 ,由卷积定理得出H(e )=1Hjd(e ) (e )-式 <3-6>理想的频率响应被窗函数的离散时间傅立叶变换采用窗函数设计法设计出来的滤波器的频率S(j) “平滑”

15、了。响应对理想响应Hd(ej*的逼近程度，由两个因素决定：8(6")主瓣的宽度。切仁侬)旁瓣的幅度大小。M2ub6vSTnP理想的情况是切仁0)主瓣的宽度窄，旁瓣的幅度小。但对于一个长度使用固定的窗函数来说，这些不能独立地达到最小。3.2 窗函数的一些通用性质(1> 、窗函数的长度N增加，主瓣的宽度减小，使得过渡带变小。 关系为:NB = C其中:B是过渡带的宽度。C是取决于窗函数的一个参数。 如矩形窗为4冗。调整N可以有效地控制过渡带的宽度，但N的改变不改 变主瓣和旁瓣的相对比例。随着 N值增加,过渡带变窄，波动频率也随着增 加，虽然总的幅度有所减少，但截止频率附近的肩峰并不

16、减少，而只是随着 N值的增加,肩峰被抑制在愈来愈小的范围内,使肩峰宽度变窄。OYujCfmUCw(2> 、窗函数的旁瓣的幅度大小取决于窗函数的选择。选择恰当的窗 函数使主瓣包含更多的能量，相应旁瓣的幅度就减小。旁瓣幅度的减小，可 以减少通带和阻带的波动，使通带尽可能趋近水平，阻带尽可能达到最大衰 减。但通常此时过渡带会变宽。eUts8ZQVRd(3> 、取不同的窗函数对幅度特性的整形效果比单纯的增加窗口长度 要强得多。3.3 设计方法这种方法也叫傅里叶级数法。一般是先给出所要求的理想的滤波器的 频率响% Hd(e，要求设计一个FIR滤波器频率响应 H (ej'nfhmNj5

17、1来逼近Hd(eje)。设计是在时域进行的，因而先由 H d(e*的而里叶反变换导出hd(n),即sQsAEJkW5Thd (n)= H d (ej°)ejO1d0 式<3-7)2 二一二由于Hd(ejO)是矩形频率响应特性，故hd(n) 一定是无限长序列，且是非因果的，而FIR滤波器的h(n)必然是有限长的，所以要用有限长的h(n)来逼近无限长的hd (n),最有效的方法是截断hd(n)或者说用一个有限长度 的窗口函数序列钏n)来截取hd(n),即 GMsIasNXkAh(n) =(n)hd(n)式<3-8)因而窗函数序列的形状及长度的选择就是关键。我们以一个截止频率为

18、6c的线性相位的理想矩形幅度特性的低通滤 波器为例来讨论。设低通特性的群延时为 a ,即e j 0<6 <6Hd(e/=Je, °c -0c式<3-9)0,Oc < co < n ,-n < o W 0c这表明，在通带网范围内，Hdleja的幅度是均匀的，其值为1,相位是一03 。hd (n)=,ej%3d82 M f, -c sin c(n - - )J二 c(n-：)式 <3-10)hd(n)是中心点在支的偶对称无限长非因果序列，要得到有限长的h(n),一种最简单的方法就是取矩形窗RN(n),即m(n) = RN(n)式 <3-11

19、)但是按照线形相位滤波器的约束，h(n)必须是偶对称的，对称中心应为长度的一半(N-1>/2 ,因而必须=(N-1>/2 ,所以有TIrRGchYzgh(n) = hd(n)0 (n)='hd(n),0 M n M N -1、0,n为其他式 <3-12)N -12将式3-9代入式3-10 ,可得h(n)=0,n为其他值此时，一定满足h(n) = h(N 一1一可这一线性相位的条件。下面求h(n)的傅里叶变换，也就是找出待求FIR滤波器的频率特7EqZcWLZNX性，以便能看出加窗处理后究竟对频率响应有何影响o按照复卷积公式，在时域是相乘、频域上是周期性卷积关系，即H

20、(e心)=工广 H d (ej9)ej(0>9d6 式<3-14)2 二-二因而H(e0逼近Hd(ej国的好坏，完全取决于窗函数的频率特性W(e%窗函数s(n)的频率特性W(e巧为N 1W(e心)=£ co(n)ejn 式<3-15)n =0对矩形窗RN (n),则有NWR(ej)=、e-j n n=0.N-1-j1-N sin() 2式 <3-16) /N、sin(w)也可表示成幅度函数与相位函数/N、 i-l( ) WN(ej°) =WR©)e2式 <3-17)其中sin(-N)Wr(6)=N2式 <3-18)sin(-)W

21、R(e9就是频域抽样内插函数，其幅度函数WR0)在0 =±2兀/N之内为一个主瓣，两侧形成许多衰减振荡的旁瓣，如果将理想频率响应也写成 lzq7IGf02E/N f j( ) Hd(ejC0) = Hd®)e 2式 <3-19)使用则其幅度函数为H d (。)= «式 <3-20)3.4 窗函数介绍实际应用的窗函数，可分为以下主要类型和特点:<1 ）矩形窗矩形窗属于时间变量的零次幕窗，函数形式为切=仃0, t <T<T式 <3-21)相应的窗谱为:2si 、 2sin TW( )二eT 式 <3-22)矩形窗使用最多，习惯

22、上不加窗就是使信号通过了矩形窗。这种窗的 优点是主瓣比较集中，缺点是旁瓣较高，并有负旁瓣，导致变换中带进了 局频干扰和泄漏，甚至出现负谱现象。zvpgeqJIhk使用图3-1矩形窗的时域及频域波形<2）三角窗亦称费杰（Fejer>窗，是幕窗的一次方形式，其函数形式 是：(t)t(1t MT0,t力式 <3-23)三角窗与矩形窗比较，主瓣宽约等于矩形窗的两倍，但旁瓣小，而且 无负旁瓣，如图3-2所示。图3-2三角窗的时域波形图<3）汉宁（Hanning> 窗汉宁窗又称开余弦窗，其时域表达式为：,i i ± i( cos=T 220,t MT；)，t MT式

23、 <3-24)相应的窗谱为:/ 、 sin T 1 sin( T 二) sin( T -二) W()：式 <3-25)T 2 IL :?T .二 T - -：由此式可以看出，汉宁窗可以看作是 3个矩形时间窗的频谱之和，或者是 3个sin（t>型函数之和，而括号中的两项相对于第一个谱窗向左、右各 移动了九/T,从而使旁瓣互相抵消，消去高频干扰和漏能。可以看出，汉宁窗主瓣加宽并降低，旁瓣则显著减小，从减小泄漏观点出发，汉宁窗 优于矩形窗。但汉宁窗主瓣加宽，相当于分析带宽加宽，频率分辨率下 降o NrpoJac3v1<4）海明（Hamming物海明窗也是余弦窗的一种，又称改进

24、的开余弦窗，其时间函数表达式为:1二 tt(0.54 0.4cosy),t <T0,t >T式<3-26)其窗谱为:sin TW( ) =1.08T0.46 "T二)"一)IL ' :T ，” -二式<3-27)海明窗与汉宁窗都是余弦窗，只是加权系数不同。海明窗加权的系数能使旁瓣达到更小。分析表明，海明窗的第一旁瓣衰减为 -42dB。海明窗的频谱 也是由3个矩形窗的频谱合成，但其旁瓣衰减速度为 20dB/（10oct> ,这比 汉宁窗衰减速度慢。海明窗与汉宁窗都是很有用的窗函数。1nowfTG4KI<5）高斯窗高斯窗是一种指数窗。

25、其时域函数为：0, t >T式 <3-28)式中a为常数，决定了函数曲线衰减的快慢。a值如果选取适当，可以使 截断点（T为有限值>处的函数值比较小，则截断造成的影响就比较小。高 斯窗谱无负的旁瓣，第一旁瓣衰减达一 55 dBo高斯窗的主瓣较宽，故而 频率分辨率低。高斯窗函数常被用来截断一些非周期信号，如指数衰减信使用号等。不同的窗函数对信号频谱的影响是不一样的，这主要是因为不同的 窗函数，产生泄漏的大小不一样，频率分辨能力也不一样。信号的截断产 生了能量泄漏，而用FFT算法计算频谱又产生了栅栏效应，从原理上讲这 两种误差都是不能消除的，但是我们可以通过选择不同的窗函数对它们的

26、 影响进行抑制。图3-3是几种常用的窗函数的时域和频域波形，其中矩形 窗主瓣窄，旁瓣大，频率识别精度最高，幅值识别精度最低；布莱克曼窗 主瓣宽，旁瓣小，频率识别精度最低，但幅值识别精度最高。fjnFLDa5Zo图3-3几种常用的窗函数的时域和频域波形对于窗函数，还有一些要求：(13dB带宽B,它是主瓣归一化的幅度下降到-3dB时的带宽。当数据长 度为N时，矩形窗主瓣两个过零点之间的宽度为4冗/N。tfnNhnE6e5(2最大边瓣峰值A(dB。3边瓣谱峰渐进衰减速度 D(dB/cot。所以，理想的窗函数应当具有最小的 B和A,和最大的D3.5窗函数法设计步骤(1)、首先是给定所要求的频率响应函数

27、Hd(ejE)；hd (n)(2)、其次，求单位冲激响应1 TT. Hd(ej )ej nd .2n(3、再次，有过渡带宽及阻带最小衰减的要求，查表选定窗函数及N的大小，一般N的大小要通过几次试探而后确定；HbmVN777sL(4、求得所设计的FIR滤波器的单位冲激响应；h(n)=6(n)hd，n=0, 1,N-1 ;N 1(5求H(ej")=£ h(n)e,5 ,检验是否满足设计要求，如不满 n为足，则需要重新设计。第四章设计实例4.1 滤波器设计指标要求用窗函数法设计一个数字低通滤波器。要求：通带截止频率:0.2冗,阻带起始频率：0.4兀,阻带最小衰减：-50dB。V7

28、l4jRB8Hs4.2 设计步骤(1设H (eja为理想线性相位滤波器H(ej )=式 <4-1)由所需低通滤波器的过渡带求出理想低通滤波器的截止数字频率=0.3冗,得出：使用儿ejndJ：""T = N -1,为缎出和位由需的移os 15sin c(n- ),n冗n二，n 二式<4-2)w2(n) =0.54 -0.4 cos n R 33n(2>由阻带衰减好用嗜R (n中过渡带会离短N值。阻带最小衰减50dB,比对6种窗函浪把座数呢壮窗函数为汉宁窗，哈明窗，汉宁窗，83lcPA59W92n.0 £ n £ 17小 33 一一W4(n

29、)=所要求的过渡带宽:9 2n2 - 17 _ n _ 3333A =0.4 兀-0.2 兀=0.2 几N=6.6tt /0.2 兀 +1=34, r =(N-1>/2=16设® 1(n)为汉宁窗，仍2”)为哈明窗，(n)为矩形窗切式可为三角 窗，以下为其满足题意的窗函数。式<4-3)式 <4-4)(3>加窗得到设计结果h(n) =hd(n) (n)<4 )进行 matlab仿真使用4.3 仿真结果(1>汉宁窗仿真结果理想单位脉冲响应hd(n)汉宁窗 w(n)实际单位脉冲响应hd(n)0-0.1010203040图4-1汉宁窗函数

30、仿真结果ap =0.0714 , as =44, tr_wdith =0.6283<2)哈明窗仿真结果理想单位脉冲响应hd(n)0-0.10-0.1010203040010203040哈明窗 w(n)幅度响应(dB)00.51图4-2哈明窗函数仿真结果ap =0.0477 ,仅可.as =52, tr_wdith =0.6283<3)矩形窗仿真结果图4-3矩形窗函数仿真结果ap =1.3722 , as=29, tr_wdith =0.6283<4)三角窗仿真结果0.3脉冲响应hd(n)bartlett(n)0.20.10010203040

31、0301020理想单位-0.140仅可.4-4三角窗函而真结果ap =0.1999, as =27,tr_wdith =0.6283 ;4.4仿真结果分析<1）汉宁窗的通带的波纹 ap =0.0714 ,阻带最小衰减 as =44,通带波纹 小，滤波器在通带中的平稳性好，并且阻带衰减很接近指标要求的50db,旁瓣衰减逐渐减小，能量更加集中在主瓣中，通带、阻带波动减小，从而 增大了阻带的衰减，能更好的达到预期目标。mZkklkzaaP<2 ）哈明窗的通带的波纹 ap =0.0477 ,阻带最小衰减as =52 ,与汉 宁窗相比通带波纹更小，滤波器在通带中的平稳性更好好，并且阻带衰减

32、更接近指标要求的50db,旁瓣衰减逐渐减小，能量更加集中在主瓣中，通 带、阻带波动减小，能更好的达到预期目标。AVktR43bpw<3）矩形窗通带波纹 ap =1.3722 ,阻带最小衰减 as =29,通带波纹太 大，影响滤波器在同代的稳定性，而且阻带衰减太小，与预期的相差很 大，旁瓣波动太大，不符合滤波器的设计指标。ORjBnOwcEd< 4 ）三角窗通带波纹 ap =0.1999 ,阻带最小衰减 as =27,通带波纹 相对较大，阻带衰减太小，不符合设计要求。2MiJTy0dTT< 5 ）对于同一性能指标下不同窗函数生成的过渡带宽和阻带最小衰减 是不同的，随着旁瓣的减小

33、，主瓣宽度相应增加了.gIiSpiue7A< 6 ）从五个窗比较三角窗，矩形窗，汉宁窗，汉明窗，窗设计的滤波器的仿真结果旁瓣衰减逐步增加，主瓣相应增宽。uEh0U1Yfmh第五章基于SIUMLINK的FIR滤波器的仿真前面我们已经介绍了 FIR滤波器的性能以及设计方法，并且实现了不 同窗函数对同一个滤波器的设计方法，matlab的siumlink工具箱中提供 了强大的功能，下面我们用 siumlink对以上所设计的滤波器进行仿真。IAg9qLsgBX使用5.1 使用FDATOOL计滤波器<1）首先根据前面用窗函数设计的滤波器确定下滤波器的指标:前面已经算过，滤波器的阶数为 34;低

34、通滤波器；采用窗函数法设计；fs取220kHZ,fc取33kH乙 取fs的归一化频率为1,则fc的归一化频率就是 1.5 o WwghWvVhPE<2 ）打开MATLAB勺FDATOOL计滤波器并保存单击 MATLAB主窗口下方的“start ”按钮，选择“ ToolBox ” 一"FilterDesign" 一"Filter Design&Analysis Tool"命令，打开FDATool,如图 5-1 所示。asfpsfpi4k图5-1设计FIR滤波器按照以上指标设计滤波器，单击“ Design Filter ”图标，指标要求的滤波器

35、就设计好了。设计好的滤波器如图6-1-2所示。ooeyYZTjj1Msgnftuide Response (dB)-205.2 滤波器的分析如果要了解滤波器是否满足要求，分析步骤如下:<1）滤波器的幅频响应Freqiut&ifiey (Hz)使用图5-2设计好的滤波器将设计好的滤波器以untitled.fda名字保存。选择 FDATool 菜单的 “Analysis ” - “Magnitude Response ”，启动幅频响应，分析如图5-3所示，X为频率，Y为幅度g<dB）。BkeGuInkxI弃H昌国i声即应X 口同IO困困给表田l噩扇国目K?-Magnitude

36、Response (dB)Current Filter Information!一 1 口口FJormalized Fr&iue-ncy a* radL«,snmple)FreciuencvMe 口nitgluUnits:r .llriimuiTi fThe atteniuation at cutatfTr duenctes is nxed 司 6 elfWindow:HannthsIT the passband gain)Butterworthq FHWindow口 口登$2内产唱仃；Structure.Order: hEble: 与aura：Dirftet-Form Fl

37、F?33VeaDesignedFilter o Specify o. 34-prrfiinr-：0 Scale Pass.File Edit Analysis Targets View Window HelpHiahoassBandslopDifferentiatorDesbgn CJ Block Parameters Digital Filter Design可城金至IftM. 一Filt-er Ml&nsiq炉rnrwww使用图5-3汉宁窗的幅频响应<2）滤波器的相频响应选择“Analysis ” - “ Phasc Rcsponsc”，启动相频响应分析，如图5-4所示。图5

38、-4布莱克曼窗设计的滤波器的相频响应通过FDATool,我们可以吧幅频响应图和相频响应的图放在一起来比较。<3）幅频响应与相频响应比较“Analysis " 一"Magnitude and Phase Response"来启动分析，如图 5-5所示 0 PgdO0sRlMo§名皂匚晏-100C.153S-6.9335-14.0208-21 10S151015Frequency (Hz)20_ Magnitude (dB) and! Phase Responseso O21 -图5-5幅频响应与相频响由图可以看出该滤波器是线性相位的FIR滤波器冲击响

39、应和阶跃响应使用“Analysis " -"Impulse Response "如图 5-6 所示“Analysis " 一"Step Response"如图 5-7 所示D1002QD 300400500 &QDTime (ms-econds)700_ pulse Response图5-6冲击响应_ Step Response _01002003g 4M iin& rmscrnG'SCO600700总,64.S D o o Q-图5-7阶跃响应由图5-6可知，滤波器的实际脉冲响应是偶对称即h(n>=h(

40、m-1-n>,故该滤波器满足 FIR线性相位的条件，该滤波器是线性相位FIR滤波器。3cdXwckm155.3 设计仿真模块打开simulink的模型窗口，构建出如图5-8所示模型.图5-8simulink 仿真图5.4 设置输入信号，并导入滤波器<1)构造以基波为主的原始信号 x<t) =15sin(10冗t >+ 10sin(30 兀 t >+15sin(600 兀 t >.<2)把设计好的滤波器导入到Digital Filter Design 模块，通过波形观察窗口和频谱观查仿真结果。h8c52WOngM导入方法：打开 Digital Filter Design 模块，“file ” 一“session ”然后选择前面保存的untitled.fda 文件。v4bdyGious5.5 对设计好的滤波器进行仿真与分析单击开始按钮开始仿真。仿真结果：使用<1）滤波前后波形图T