任务11：用窗函数法设计偶对称线性相位FIR高通滤波器课件

1、成 绩数字信号处理课程设计题目：任务11：用窗函数法设计偶对称线性相位FIR高通滤波器班级学号姓名任务分工工作量(%)12级电气x班xxxxxxxxxxx工作安排与报告编辑12.312级电气x班xxxxxxxxxxx报告编辑12.412级电气x班xxxxxxxxxxx报告排版与设计12.412级电气x班xxxxxxxxxxx仿真分析12.412级电气x班xxxxxxxxxxxMATLAB程序,仿真,排版13.512级电气x班xxxxxxxxxxx原始资料搜集与筛选12.412级电气x班xxxxxxxxxxx报告审核、编辑补充12.312级电气x班xxxxxxxxxxx报告排版与修改12.320

2、14年12月用窗函数法设计偶对称线性相位FIR高通滤波器摘要：所谓数字滤波器，是指输入、输出均为数字信号，通过数值运算处理改变输入信号所含频率成分的相对比例，或者滤除某些频率成分的数字器件或程序。数字滤波器处理精度高、稳定、体积小、重量轻、灵活、不存在阻抗匹配问题，可以实现模拟滤波器无法实现的特殊滤波功能。常用的数字滤波器课分为无限长滤波器和有限长滤波器。这里我们主要介绍有限长数字滤波器FIR的类型以及其窗函数设计方法的基本原理思路。本设计介绍的是用窗函数设计方法设计数字频带变换的高通滤波器。简单介绍了几种典型窗函数的基本参数。再设计过程中根据给定的数字高通滤波器的参数我们选定的窗函数是汉宁窗

3、。基本思路是：首先将高通滤波器的参数根据公式转换成低通数字滤波器，求出低通数字滤波器的系统函数，再根据低通与高通之间的映射关系转换成高通数字滤波器的系统函数。得到高通数字滤波器的实现结构，并进行结构及有限字长的分析。确定最适合的滤波器的结构及字长。但是在设计过程中由于没有找到相关的数字低通与数字高通在数字频带内的转化关系式，所以采用的是直接设计，没有进行频带变换；而且由于高通数字滤波器参数设置不是很符合理想状态，求出的N为31，得到的数字滤波器阶数太高，得不出滤波器系统函数H（z）的闭合函数，所以该滤波器的结构只能以直接型实现，没有进行结构对比分析。关键字：高通滤波器、窗函数、数字频带变换、海

4、明窗、布莱克曼窗目录1、简述低通、高通、带通、带阻滤波器的性能特点- 1 -1.1 低通滤波器- 1 -1.2 高通滤波器- 1 -1.3 带通滤波器- 1 -1.4 带阻滤波器- 1 -2、比较模拟巴特沃斯滤波器和切比雪夫滤波器- 2 -2.1 巴特沃斯滤波器- 2 -2.2 切比雪夫滤波器- 2 -3、简述巴特沃斯滤波器和切比雪夫模拟低通滤波器的设计方法- 3 -3.1巴特沃斯滤波器的设计- 3 -3.2切比雪夫滤波器设计- 5 -4、简述模拟低通滤波器转换为数字滤波器（低通、高通、带通、带阻）的设计流程- 7 -5、主要设计内容- 9 -5.1一般

5、模拟低通滤波器：- 9 -5.1.1巴特沃斯低通逼近：- 9 -5.1.2切比雪夫低通逼近：- 10 -5.2设计FIR滤波器的原理- 11 -5.3 各种窗函数的具体参数- 12 -6、设计流程- 14 -6.1采用海明窗- 15 -6.2 采用布莱克曼窗- 16 -7、MATLAB仿真- 18 -7.1 第一种滤波信号fn=sin(2*pi*f1*t)+sin(2*pi*f2*t)+sin(2*pi*f3*t)- 18 -7.1.1 采用海明窗函数- 18 -7.1.2 采用布莱克曼窗函数- 21 -7.2 第二种滤波信号fn=square(1:256,50)+randn(1,256)-

6、24 -7.2.1 采用海明窗函数- 24 -7.2.2 采用海明窗函数- 27 -附件：MATLAB程序附录- 31 -1、freqz-m子函数定义- 31 -2、ideal-lp子函数定义- 31 -3、海明窗高通滤波器- 32 -4、布莱克曼窗高通滤波器- 34 -参考文献- 37 -1、简述低通、高通、带通、带阻滤波器的性能特点 滤波器是一种选频装置，可以使信号中特定的频率成分通过，而极大地衰减其它频率成分。在测试装置中，利用滤波器的这种选频作用，可以滤除干扰噪声获得人们想要的频率范围内的信号或进行频谱分析。根据滤波器的选频作用滤波器可分为以下几种： 1.1 低通滤

7、波器   从0f2频率之间，幅频特性平直，它可以使信号中低于f2的频率成分几乎不受衰减地通过，而高于f2的频率成分受到极大地衰减。 1.2 高通滤波器  与低通滤波相反，从频率f1，其幅频特性平直。它使信号中高于f1的频率成分几乎不受衰减地通过，而低于f1的频率成分将受到极大地衰减。 1.3 带通滤波器  它的通频带在f1f2之间。它使信号中高于f1而低于f2的频率成分可以不受衰减地通过，而其它成分受到衰减。 1.4 带阻滤波器 与带通滤波相反，阻带在频率f1f2之间。它使信号中高于f

8、1而低于f2的频率成分受到衰减，其余频率成分的信号几乎不受衰减地通过。其中低通滤波器和高通滤波器是滤波器的两种最基本的形式，其它的滤波器都可以分解为这两种类型的滤波器，例如：低通滤波器与高通滤波器的串联为带通滤波器，低通滤波器与高通滤波器的并联为带阻滤波器。2、比较模拟巴特沃斯滤波器和切比雪夫滤波器2.1 巴特沃斯滤波器 从幅频特性提出要求，而不考虑相频特性。巴特沃斯滤波器具有最大平坦幅度特性，其幅频响应表达式为：其中n为滤波器的阶数；为滤波器的截止角频率。 当w=时，=1/2,所以，对应的是滤波器的-3db点。巴特沃思低通滤波器是以巴特沃思函数作为滤波器的传递函数H(s)，以最高阶

9、泰勒级数的形式逼近滤波器的理想矩形特性。2.2 切比雪夫滤波器  切贝雪夫滤波器也是从幅频特性方面提出逼近要求的，其幅频响应表达式为：   其中是决定通带波纹大小的波动系数，0<<1，波纹的产生是由于实际滤波网络中含有电抗元件； Wc是通带截止频率，Tn是n阶切贝雪夫多项式。与巴特沃斯逼近特性相比较，这种特性虽然在通带内有起伏，但对同样的n值在进入阻带以后衰减更陡峭，更接近理想情况。值越小，通带起伏越小，截止频率点衰减的分贝值也越小，但进入阻带后衰减特性变化缓慢。 巴特沃斯滤波器，通带内没有任何的纹波，因而也被称为最大平坦度滤波器；但是其幅度

10、响应曲线在截止频率处也较为平缓。也就是说巴特沃斯滤波器的矩形系数较低。 而切比雪夫滤波器，最大的优点是滤波器在通带外拥有陡峭的衰减曲线，有较高的矩形系数；缺点是通带内有纹波，并且，通带内纹波越大，通带到阻带的过渡就越陡峭，矩形系数就越高。切贝雪夫滤波器与巴特沃斯滤波器进行比较，切贝雪夫滤波器的通带有波纹，过渡带轻陡直，因此，在不允许通带内有纹波的情况下，巴特沃斯型更可取；从相频响应来看，巴特沃斯型要优于切贝雪夫型，通过上面二图比较可以看出，前者的相频响应更接近于直线。3、简述巴特沃斯滤波器和切比雪夫模拟低通滤波器的设计方法3.1巴特沃斯滤波器的设计 常用设计巴特沃斯低通滤波器指标：

11、：通带截止频率；：通带衰减，单位：dB；：阻带起始频率；：阻带衰减，单位：dB。说明：（1）衰减在这里以分贝（dB）为单位，即（2）当时为通常意义上的截止频率。（3）在滤波器设计中常选用归一化的频率，即根据设计指标要求，确定归一化巴特沃斯低通滤波器幅度平方函数中的待定系数C及滤波器的阶数N；然后再根据幅度平方函数确定巴特沃斯低通滤波器的传递函数H(s)。（1）将实际频率归一化得，再根据已知的，幅度平方函数 确定C和N。（2） 求C和N ：由并代入 ，得： 由于，故 。由 两边取对数得： 其中，这样可以求出C和N。（3）确定巴特沃斯滤波器的传递函数H(p)。由于，由，解得极点为：将p左半平面的极

12、点赋予即其中为了便于设计，工程上已将当时，各阶巴特沃斯低通滤波器系统函数设计成表格供查阅，如表3-1所示。表3-1 归一化巴特沃斯模拟低通滤波器系统函数表阶次归一化系统函数12345（4）去掉归一化影响 上面设计中采用归一化的频率即，而实际中截止频率为，所以要进行如下的变量代换： 即 综上，归纳出设计巴特沃斯低通滤波器的方法如下：（1）计算归一化频率，。（2）根据设计要求按照和其中计算巴特沃斯滤波器的参数C和阶次N；注意当时 C=1。（3）利用N查表获得归一化巴特沃斯低通原型滤波器的系统函数；（4）令中的得到截止频率为的巴特沃斯低通滤波器的系统函数。3.2切比雪夫滤波器设计已知切比雪夫归一化滤

13、波器的幅度平方函数为 由上式可知，要确定切比雪夫滤波器的幅度平方函数，需要确定三个参数：及N。下面研究如何确定这三个参数，具体步骤如下：（1）将实际频率归一化得 ，再根据已知的，幅度平方函数 确定和N。（2）确定和N：定义通带波纹(即通带衰减)（以分贝为单位）为：代入 ，得即 ，所以，则，其中。这样可以求出和N，其中。在已知、的情况下，就可以根据幅度平方函数求出滤波器的零点和极点，从而确定滤波器的系统函数。表3-2 归一化切比雪夫原型滤波器分母多项式设计系数n波纹（）12.862775221.51620261.425624530.71569381.53489541.252913040.3790

14、5061.02545531.71686621.197385650.178923400.75251811.30957471.93736751.1724909波纹（）11.965226721.10251031.097734330.49130671.23840920.988341240.27562760.74261941.45392480.952811450.12282670.58053420.97439611.68881600.9368201波纹（）11.307560320.63676810.803816430.32689011.02219030.737821640.20576510.516798

15、11.25648190.716215050.08172250.45934910.69347701.49954330.7064606为了设计方便，工程上已将截止频率的切比雪夫低通滤波器的系统函数设计为表格供设计时查阅。归一化原型切比雪夫低通滤波器的系统函数如（5-18）式所示，设计表格如表5-3所示。 （5-18） 再次强调，表5-3是归一化的结果，对于具体的，其系统函数可由（5-19）式得到。 （5-19）综上所述，设计切比雪夫低通滤波器的基本步骤如下：（1）计算归一化频率，。（2）根据通带波纹（通带衰减）db，按照式计算。（3）根据阻带起始频率，阻带衰减和。按照其中式计算滤波器的阶数N。（4

16、）根据滤波器阶数N，查表得归一化原型切比雪夫滤波器系统函数；根据的低频特性求出待定系数，注：当N为偶数时，；当N为奇数时，。（5）去掉归一化影响 根据截止频率，按照式计算切比雪夫滤波器的系统函数。4、简述模拟低通滤波器转换为数字滤波器（低通、高通、带通、带阻）的设计流程实际应用中的数字滤波器有低通、高通、带通、带阻等类型。设计各类数字滤波器有以下两种方法：（1）把一个归一化原型模拟低通滤波器经模拟频带变换成所需要类型（包括高通、带通、带阻与另一截止频率的低通）的模拟滤波器，然后再通过冲激响应不变法或双线性变换法数字化为所需类型的数字滤波器，如图3-1所示。实际设计时是把这一方法中的两步合成一步

17、来实现，即把模拟归一化低通原型变换到模拟低通、高通、带通、带阻等滤波器的公式与双线性变换得到相应数字滤波器的公式合并，就可直接从模拟低通归一化原型通过一定的频率变换关系，一步完成各类型数字滤波器的设计，如图3-2所示。这里只采用了双线性变换法。因为冲激响应不变法有频率响应混叠失真效应，只适用于严格限带的数字低通、带通滤波器；而对于数字高通、带阻滤波器，则不能直接应用。（2）先利用冲激响应不变法或双线性变换法，把模拟低通原型数字化成数字低通滤波器，然后利用数字频带变换法，将它变换成所需要的各型数字滤波器（另一截止频率的数字低通、数字高通、数字带通、数字带阻等），如图3-3所示。模拟归一化原型c=

18、1模拟-模拟频带变换数字化模拟归一化原型c=1数字低通、高通、带通、带阻模拟低通、高通、带通、带阻数字-数字频带变换数字化数字低通、高通、带通、带阻数字低通图3-1 先模拟频带交换，再数字化模拟归一化原型c=1数字低通、高通、带通、带阻图3-2把3-1的两步合成一步直接设计图3-3 先数字化，再进行数字频带变换5、主要设计内容利用窗函数法设计FIR滤波器，绘制出滤波器的特性图。利用所设计的滤波器对多个频带叠加的正弦信号进行处理，对比滤波前后的信号时域和频域图，验证滤波器的效果。1、学习用窗函数法设计FIR数字滤波器的原理及其设计步骤；2、用Matlab软件对FIR数字滤波器数学模型的系数进行求

19、解；并用窗函数法完成高通滤波器的设计；3、用海明窗、布莱克曼窗完成FIR数字滤波器的设计，对两种信号进行滤波，并对结果进行比较和分析。5.1一般模拟低通滤波器：5.1.1巴特沃斯低通逼近：=巴特沃斯逼近又称为最平幅度逼近，巴特沃斯低通滤波器幅度平方函数定义为： 式中N为正整数，代表滤波器的阶次，成为截止频率。当时，有， 所以又称为巴特沃斯滤波器的3分贝带宽。巴特沃斯滤波器有如下特点：a. 当=0时，即=0处无衰减。b. 当=c时， 通带最大衰减在止带内的逼近是单调变化的，不管N为多少，所有幅频特性曲线都经过-3dB点，或说衰减3dB,这就是3dB不变性。c.在<c的通带内：有最大平坦的幅

20、度特性，即N阶Butterworth低通滤波器在=0处：前(2N-1)阶导数为零，因而Butterworth又称最平幅度特性滤波器。随着由0变到c ，单调减小，N越大，减小越慢，也就是通带内特性越平坦。d.在>c，即在过渡带及阻带中，也随增加而单调减小，但是/c >1，故比通带内衰减的速度要快得多，N越大，衰减速度越大。当=s，即频率为阻带截止频率时，衰减为： 为阻带最小衰减。 在已知c 、 s和=p的衰减Ap的情况下，求Butterworth DF阶数N： 5.1.2切比雪夫低通逼近：切比雪夫滤波器的幅度特性就在一个频带中（通带或阻带）具有这种等纹波特性。具可分为切比雪夫I型和切

21、比雪夫II型。各自特点分别为：a. Chebyshev I型：在通带中是等波纹的，在阻带内是单调的；b. Chebyshev II型：在通带中是单调的，在阻带内是等波纹的；由应用的要求，决定采用哪种型式的Chebyshev滤波器。Chebyshev I型模拟滤波器的振幅平方函数为： ，他是表示通带纹波大小的一个参数，越大，纹波也越大。为对的归一化频率，为截止频率，也是滤波器的某一衰减分贝出的通带宽度（这一分贝数不一定是3分贝，也就是说，在切比雪夫滤波器中，不一定是3分贝的带宽）。是N阶切比雪夫多项式，定义为 切比雪夫滤波器的幅度响应所具有的特点如下：a. 当=0，N为偶数是，；当N为奇数，。b

22、. 当有 即 时，所有幅度响应曲线都通过点，则该点为切比雪夫滤波器的截止频率。在这个截止频率下，复读函数不一定下降3分贝。c.在通带内，当时，幅度响应在1之间等纹波起伏 d.在通带外部，当时，随着的增大 使幅度响应迅速趋近于零。切比雪夫滤波器参数求解： 求，当时，所得到的即为 求，与通带纹波有关，应该有； 求阶数N，可得最后，取滤波器阶数为大于由上式求得的N的一个整数。 5.2设计FIR滤波器的原理如果所希望的滤波器的理想的频率响应函数为，则其对应的单位脉冲响应为 （5.1）窗函数设计法的基本原理是用有限长单位脉冲响应序列逼近。由于往往是无限长序列，而且是非因果的，所以用窗函数将截断，并进行加

23、权处理，得到： （5.2）就作为实际设计的FIR数字滤波器的单位脉冲响应序列，其频率响应函数为 （5.3）式中，N为所选窗函数的长度。我们知道，用窗函数法设计的滤波器性能取决于窗函数的类型及窗口长度N的取值。设计过程中，要根据对阻带最小衰减和过渡带宽度的要求选择合适的窗函数类型和窗口长度N 。各种类型的窗函数可达到的阻带最小衰减和过渡带宽度见表1。这样选定窗函数类型和长度N之后，求出单位脉冲响应，并按照式（5.3）求出。是否满足要求，要进行演算。一般在尾部加零使长度满足2的整数次幂，以便用FFT计算。如果要观察细节，补零点数增多即可。如果不满足要求，则要重新选择窗函数类型和长度N，再次验算，直

24、至满足要求。如果要求线性相位特性，则还必须满足：根据上式中的正、负号和长度N的奇偶性又将线性相位FIR滤波器分成四类。要根据所设计的滤波特性正确选择其中一类，例如，要设计线性相位低通特性，可以选择这一类，而不能选择这一类。5.3 各种窗函数的具体参数表1 各种窗函数的基本参数窗函数旁瓣峰值幅度/dB过渡带宽阻带最小衰减/dB矩形窗-134/N-12三角形窗-258/N-25汉宁窗-318/N-44海明窗-418/N-53布拉克曼窗-5712/N-74凯泽窗(=7.865)-5710/N-80设计思路：1FIR滤波器的设计方法和IIR滤波器的设计方法有很大不同。FIR DF设计的含义是：根据设计

25、指标，求解所选运算结构要求的h(n)或H(z)：2线性卷积和快速卷积型结构，求FIR DF的h(n)。3级联和频率采样型结构，求FIR DF 的H(z)。(1)先给定所要求设计的理想滤波器的频率响应Hd(ejw)。(2)设计一个可实现的FIR滤波器频率响应H(ejw)。(3)由于设计是在时域中进行，使所设计滤波器的h(n)去逼近理想单位取样响应hd(n)。6、设计流程开始读入窗口长度N计算hd(n)调用窗函数子程序求w(n)调用freqz_m函数计算H(k)调用绘图子程序（函数）绘制H(k)幅度相位曲线结束 设计流程图计算h(n)= hd(n) w(n)由于窗函数法的通带、阻带纹波是一样的，最

26、大值都是，用（6.1）式可由可求出，即： （6.1）=50dB，则。由（6.2）式可得=0.055dB,由于题目中给出的0.1dB的要求，故设计时直接采用所给阻带衰减的指标。否则，若由算出的不满足要求，则要利用通带指标重新计算阻带所需衰减，以此作为设计的出发点。 （6.2）6.1采用海明窗（1） 理想线性相位高通滤波器及分别见于（6.3）、（6.4）两式，其中。 ， = 0 ， （6.3） ， 1- ， （6.4）对于高通DF，可直接利用（6.5）式得理想线性相位高通滤波器截止频率 （6.5）（2） 求窗函数由阻带衰减，查表海明窗阻带衰减为53dB，满足要求，由过渡带宽度为确定滤波器长度点数N

27、，而海明窗过渡带，则有，由于是线性相位高通滤波器，只能取N=奇数，为了满足过渡带的要求，取大一点的N，取N=35, 则，于是海明窗（3） 求线性相位FIR高通滤波器的h(n)。先将参数带入上式（6.4）中的中，可得 0.7 , 则有： ， 0.7， （4） 求，可得系统幅度响应如下图（5） 检验的各项指标，满足要求。6.2 采用布莱克曼窗（1） 理想线性相位高通滤波器及分别见于（6.6）、（6.7）两式，其中。 ， = 0 ， （6.6） ， 1- ， （6.7）对于高通DF，可直接利用（6.8）式得理想线性相位高通滤波器截止频率 （6.8）（2） 求窗函数由阻带衰减，查表布莱克曼窗阻带衰减为

28、74dB，满足要求，由过渡带宽度为确定滤波器长度点数N，而海明窗过渡带，则有，由于是线性相位高通滤波器，只能取N=奇数，为了满足过渡带的要求，取大一点的N，取N=57, 则，于是海明窗为（3） 求线性相位FIR高通滤波器的h(n)。先将参数带入上式（6.7）中的中，可得 0.7 , 则有： ， 0.7 ， （4） 求，可得系统幅度响应（5） 检验的各项指标是否满足要求。7、MATLAB仿真7.1 第一种滤波信号fn=sin(2*pi*f1*t)+sin(2*pi*f2*t)+sin(2*pi*f3*t)7.1.1 采用海明窗函数海明窗函数图形为：(1).滤波器的幅度响应曲线：(2).时域波形和

29、频域波形：(3).滤波后的时域和频域波形：总的效果图：（4）滤波结果分析由MATLAB计算出：海明窗的窗宽N=35采样频率s=100000通带截止频率为p=ps2=0.4×1000002=20000Hz阻带起始频率为st=sts2=0.2×1000002=10000Hz截止频率为c=p+st2=15000Hz实际通带最大衰减Ap =0.0287dB实际最小阻带衰减 As = 55dB滤波效果来看，对比滤波前后信号频谱发现，频率小于截止频率15000Hz的信号基本被滤去，通带截止频率以上的频率幅值稍微有所衰减，但整体效果达到了预期标准。7.1.2 采用布莱克曼窗函数布莱克曼窗

30、函数为; (1).滤波器的幅度响应曲线： (2).f(n)的时域波形和频域波形：(3).滤波后的时域和频域波形：总的效果图： (4)滤波结果分析：由MATLAB计算出：布莱克曼窗的窗宽N=57采样频率s=100000通带截止频率为p=ps2=0.4×1000002=20000Hz阻带起始频率为st=sts2=0.2×1000002=10000Hz截止频率为c=p+st2=15000Hz实际通带最大衰减Ap =0.003dB实际最小阻带衰减 As = 76dB滤波效果来看，对比滤波前后信号频谱发现，频率小于截止频率15000Hz的信号基本被滤去，通带截止频率以上的频率幅值与采

31、用海明窗相比衰减更大，但整体效果达到了预期标准。7.2 第二种滤波信号fn=square(1:256,50)+randn(1,256)7.2.1 采用海明窗函数海明窗函数图形为：(1).滤波器的幅度响应曲线： (2)时域波形和频域波形： (3).滤波后的时域和频域波形：总的效果图：（4）滤波结果分析由MATLAB计算出：海明窗的窗宽N=35采样频率s=100000通带截止频率为p=ps2=0.4×1000002=20000Hz阻带起始频率为st=sts2=0.2×1000002=10000Hz截止频率为c=p+st2=15000Hz实际通带最大衰减Ap =0.0287dB实

32、际最小阻带衰减 As = 55dB滤波效果来看，对比滤波前后信号频谱发现，频率小于截止频率15000Hz的信号基本被滤去，但大于10000Hz小于20000Hz的过渡带宽内仍存在少量的噪声信号，通带截止频率以上的频率幅值稍微有所衰减，但整体效果达到了预期标准。7.2.2 采用海明窗函数布莱克曼窗函数为; (1).滤波器的幅度响应曲线： (2).f(n)的时域波形和频域波形： (3).滤波后的时域和频域波形：总的效果图：(4)滤波结果分析：由MATLAB计算出：布莱克曼窗的窗宽N=57采样频率s=100000通带截止频率为p=ps2=0.4×1000002=20000Hz阻带起始频率为

33、st=sts2=0.2×1000002=10000Hz截止频率为c=p+st2=15000Hz实际通带最大衰减Ap =0.003dB实际最小阻带衰减 As = 76dB滤波效果来看，对比滤波前后信号频谱发现，频率小于截止频率15000Hz的信号频谱比采用海明窗抑制得更厉害，噪声信号明显降低，过渡带宽与海明窗相比变宽，通带截止频率以上的频率幅值与采用海明窗相比衰减更大，但整体效果达到了预期标准。综上所述，对比布拉克曼窗与海明窗的对数幅度响应曲线可以看出，旁瓣衰减加大，但主瓣宽度也相应加宽了。而且最小阻带衰减只由窗形状决定，不受N的影响；而过渡带的宽度则既和窗形状有关，且随窗宽N的增加而

34、减少。对于用窗函数法设计滤波器，可以总结一下几点：1) 窗函数频谱的主瓣应尽可能地窄, 即能量尽可能集中在主瓣内, 以提高谱估计时的频域分辨率和减小泄漏, 在数字滤波器设计中获得较小的过渡带。2) 尽量减少窗函数频谱的最大旁瓣的相对幅度, 以使旁瓣高度随频率尽快衰减，如这两条不能同时得到满足, 往往是增加主瓣宽度以换取对旁瓣的抑制。总之, 在应用窗函数时, 除要考虑窗函数频谱本身的特性外, 还应充分考虑被分析信号的特点及具体处理要求。3)从上面的两个窗函数的仿真，可以明显看出用不同的窗函数设计的滤波器满足的指标相差很大，所以在设计滤波器时要根据具体的指标来选择相应的窗函数，而且要经过多次仿真验

35、证，最后选最佳方案。附件：MATLAB程序附录1、freqz-m子函数定义%freqz_m的定义：function db,mag,pha,grd,w=freqz_m(b,a);%求取系统的绝对幅度响应、相对的db值幅度响应、相位响应和群延时响应的函数%db为相对振幅(dB)%mag为绝对振幅%pha为相位响应%grd为群延时%w为频率样本点向量%b为Ha(z)分子多项式系数(对FIR而言，b=h)%a为Hz(z)分母多项式系数(对FIR而言，a=1)%H,w=freqz(b,a,1000,'whole');%freqz显示数字滤波器频域中的图形%H,W?=?FREQZ(B,A,

36、N,'whole')uses N points around the whole unit circle.H=(H(1:501)'w=(w(1:501)'mag=abs(H);db=20*log10(mag+eps)/max(mag);%dB=-20log10|H(ejw)|/|H(ejw)|max>=0pha=angle(H);grd=grpdelay(b,a,w);2、ideal-lp子函数定义 function hd=ideal_lp(wc,M);%Ideal Lowpass filter computation%hd=ideal_lp(wc,M)%

37、 hd=ideal impulse response between 0 to M-1% wc=cutoff frequency in radians% M=length of the ideal filter%alpha=(M-1)/2;n=0:1:(M-1);m=n-alpha+eps;hd=sin(wc*m)./(pi*m);3、海明窗高通滤波器clc;clear all;clf;%初始条件wp=0.4*pi;deltaw=0.2*pi;ws=wp-deltaw;N0=ceil(6.6*pi/deltaw); %布莱克曼窗11*pi/N；海明窗6.6*pi/NN=N0+mod(N0+1,

38、2)+2n=0:N-1;wd=(hamming(N)' %海明窗wc=(wp+ws)/2;hd=ideal_lp(pi,N)-ideal_lp(wc,N); %理想高通的冲激响应h=hd.*wd; %截取得到的实际单位脉冲响应%求取系统的绝对幅度响应、相对的db值幅度响应、相位响应和群延时响应的函数db,mag,pha,grd,w=freqz_m(h,1);%dw=2*pi/1000;Ap=-min(db(wp/dw+1:501) %计算实际通带最大衰减 As=-round(max(db(1:ws/dw+1) %计算实际最小阻带衰减 %海明窗函数绘图figure(1);dbh,magh

39、,phah,grdh,wh=freqz_m(wd,1);subplot(121);plot(n,wd);grid on;title('海明窗');xlabel('n');ylabel('w(n)');subplot(122);plot(w/pi,dbh);title('海明窗的对数幅度响应曲线（dBb）');xlabel('w/pi');ylabel('dBb');axis(0,1,-130,5);grid on;%滤波器幅度响应曲线等图figure(2);subplot(221);stem(n,h

40、d);grid on; %plot(n,hd)title('理想单位秒冲响应 hd(n) ');subplot(222);stem(n,h);grid on; %plot(n,h)title('实际单位脉冲响应 h(n) ');subplot(223);plot(n,wd);grid on;title('海明窗');xlabel('n');ylabel('w(n)');axis(0,N,0,1);subplot(224);plot(w/pi,db);grid on;title('高通滤波器的幅度响应曲线（d

41、B）');xlabel('w/pi');ylabel('dB');axis(0,1,-130,5);%滤波程序figure(3)f1=10000;%信号频率Hz（低频信号）f2=20000;%信号频率Hzf3=40000;%信号频率Hzfs=100000; %采样频率HzN=256;%采样点数t=(0:N-1)/fs;%采样时间sfn=sin(2*pi*f1*t)+sin(2*pi*f2*t)+sin(2*pi*f3*t);%信号采样值% fn=square(1:256,50)+randn(1,256);%方波信号叠加高斯噪声的信号f(n)% 注意：若要

42、使用题目中所给的f(n)方波叠加噪声信号，只需将上一行代码中的三个%号去掉. .即去除注释就可覆盖前几行所设置的f(n)信号。n=1:256;subplot(221);plot(fn);grid on;title('噪声f(n)时域波形')FN=fft(fn,512);w=(0:255)/256*(fs/2);subplot(222);plot(w,abs(FN(1:256);title('噪声f(n)频域波形')xlabel('频率(Hz)');ylabel('频率响应'); grid on;hfn=filter(h,1,fn)

43、;subplot(223);plot(hfn);grid on;title('滤波后时域波形');SF=fft(hfn,512); w=(0:255)/256*(fs/2);subplot(224);plot(w,abs(SF(1:256);grid on; title('滤波后频域波形');xlabel('频率(Hz)');ylabel('响应'); grid on;4、布莱克曼窗高通滤波器clc;clear all;clf;%初始条件wp=0.4*pi;deltaw=0.2*pi;ws=wp-deltaw;N0=ceil(11

44、*pi/deltaw); %布莱克曼窗11*pi/N；海明窗6.6*pi/NN=N0+mod(N0+1,2)+2n=0:N-1;wd=(blackman(N)' %布莱克曼窗wc=(wp+ws)/2;hd=ideal_lp(pi,N)-ideal_lp(wc,N); %理想高通的冲激响应h=hd.*wd; %截取得到的实际单位脉冲响应%求取系统的绝对幅度响应、相对的db值幅度响应、相位响应和群延时响应的函数db,mag,pha,grd,w=freqz_m(h,1);%dw=2*pi/1000;Ap=-min(db(wp/dw+1:501) %计算实际通带最大衰减 As=-round(max(db(1:ws/dw+1) %计算实际最小阻带衰减 %布莱克曼窗函数绘图figure(1);dbb,magb,phab,grdb,wb=freqz_m(wd,1);subplot(121);plot(n,wd);grid on;title('布莱克曼窗');xlabel('n');ylabel('w(n)');subplot(122);plot(w/pi,dbb);title('布莱克曼窗的