窗函数法设计FIR数字低通滤波器设计毕业设计doc.docx

1、窗函数法设计FIR数字低通滤波器设计毕业设计doc基于窗函数法的FIR数字低通滤波器设计摘 要数字滤波器是一种用来过滤时刻离散信号的数字系统，通过对抽样数据进行数学处置来达到频域滤波的目的。按照其单位冲激响应函数的时域特性可分为两类：无穷冲激响应（IIR）滤波器和有限冲激响应（FIR）滤波器。与IIR滤波器相较，FIR的实现是非递归的，老是稳固的；更重要的是，FIR滤波器在知足幅频响应要求的同时，能够取得严格的线性相位特性。因此，它在高保真的信号处置，如数字音频、图像处置、数据传输、生物医学等领域取得普遍应用。滤波器的设计是信号处置的核心问题之一。按照FIR滤波器的原理，提出了FIR滤波器的窗

2、函数设计法，给出了在MATLAB环境下，用窗函数法设计FIR滤波器的进程和设计实例。通过利用不同的窗函数方式设计FIR滤波器，对所设计的滤波器进行分析比较，得出各类方式设计的滤波器的优缺点及其不同的利用处合，从而能够在设计滤波器时能够正确的选择FIR数字滤波器的窗函数的选取及设计方式。关键词：FIR滤波器，MATLAB，窗函数1 概述 FIR滤波器简介 FIR数字滤波器设计最简单的方式是窗函数法，通常也称为傅立叶级数法。它是在时域进行的，因此必需由理想滤波器的频率响应推导出其单位冲激响应，在设计一个FIR数字滤波器的单位冲激响应去逼近。按照冲激响应的时域特性，数字滤波器可分为无穷长冲激响应（I

3、IR）和有限长冲激响应滤波器（FIR），FIR的突出长处是：系统老是稳固的、易于实现线性相位、允许设计多通带（或多阻带）滤波器，但与IIR相较，在知足一样阻带衰减的情形下需要的阶数较高，滤波器的阶数越高，占用的运算时刻越多，因此在知足指标要求的情形下应尽可能减少滤波器的阶数。 窗函数设计法窗函数设计法是一种通过截短和计权的方式使无穷长非因果序列成为有限长脉冲应响应序列的设计方式，通常在设计滤波器之前，应该先按照具体的工程应用肯定滤波器的技术指标，在大多数实际应用中，数字滤波器常常被用来实现选频操作，所以指标的形式一般为在频域中以分贝值给出的相对幅度响应和相位响应。 用窗函数法设计FIR滤波器的

4、步骤如下： （1）根据过渡带宽及阻带衰减要求，选择窗函数的类型并估计窗口长度N（或阶数M=N-1），窗函数类型可根据最小阻带衰减As独立选择，因为窗口长度N对最小阻带衰减As没有影响，在确定窗函数类型以后，可根据过渡带宽小于给定指标确定所拟用的窗函数的窗口长度N，设待求滤波器的过渡带宽为，它与窗口长度N近似成反比，窗函数类型确定后，其计算公式也确定了，不过这些公式是近似的，得出的窗口长度还要在计算中逐步修正，原则是在保证阻带衰减满足要求的情况下，尽量选择较小的N，在N和窗函数类型确定后，即可调用MATLAB中的窗函数求出窗函数W（n）。 （2）根据待求滤波器的理想频率响应求出理想单位脉冲响应h

5、d（n），若是给出待求滤波器频率应为Hd（ejw），则理想的单位脉冲响应能够用下面的傅里叶反变换式求出： （1）在一般情形下，hd（n）是不能用封锁公式表示的，需要采用数值方式表示；从w=0到w=2采样N点，采用离散傅里叶反变换（IDFT）即可求出。 （3）计算滤波器的单位脉冲响应h（n），它是理想单位脉冲响应和窗函数的乘积。 (4)算技术指标是不是知足要求，为了计算数字滤波器在频域中的特性，可挪用freqz子程序，若是不知足要求，可按照具体情形，调整窗函数类型或长度，直到知足要求为止。使用窗函数法设计时要满足以下两个条件：窗谱主瓣尽可能地窄，以取得较陡的过渡带；尽可能减少窗谱的最大旁瓣的相对

6、幅度，也就是使能量尽可能集中于主瓣，减小峰肩和纹波，进行增加阻带的衰减。窗函数的选择原则是：1具有较低的旁瓣幅度，尤其是第一旁瓣的幅度；2旁瓣的幅度下降的速度要快，以利于增加阻带的衰减；3主瓣的宽度要窄，如此能够取得比较窄的过渡带。通常上述的几点难以同时知足。实际中设计FIR数字滤波器往往要求是线性相位的，因此要求知足线性相位的条件，即要求w(n)知足： （2）所以，窗函数不仅有截短的作用，而且能够起到光滑的作用，在很多领域取得了应用。2 设计原理 大体原理设计低通FIR数字滤波器，寻求一系统函数，使其频率响应逼近滤波器要求的理想频率响应，其对应的单位脉冲响应 （3）其中 若是所希望的滤波器的

7、理想的频率响应函数为，则其对应的单位脉冲响应为 （4）窗函数设计法的大体原理是用有限长单位脉冲响应序列逼近。由于往往是无穷长序列，而且是非因果的，所以用窗函数将截断，并进行加权处置，取得： （5）就作为实际设计的FIR数字滤波器的单位脉冲响应序列，其频率响应函数为 （6）式中，N为所选窗函数的长度。 典型的窗函数（1）矩形窗(Rectangle Window) （7） 其频率响应和幅度响应别离为：， （8）（2）三角形窗(Bartlett Window) （9） 其频率响应为： （3）汉宁(Hanning)窗，又称升余弦窗 （10）其频率响应和幅度响应别离为：（4）汉明(Hamming)窗，又

8、称改良的升余弦窗 （11）其幅度响应为： （5）布莱克曼(Blankman)窗，又称二阶升余弦窗 （12）其幅度响应 ： （6）凯塞(Kaiser)窗 （13） 其中：是一个可选参数，用来选择主瓣宽度和旁瓣衰减之间的互换关系，一般说来，越大，过渡带越宽，阻带越小衰减也越大。I0()是第一类修正零阶贝塞尔函数。 若阻带最小衰减表示为，的肯定可采用下述经验公式： 若滤波器通带和阻带波纹相等即p=s时，滤波器节数可通过下式肯定： 式中：咱们明白，用窗函数法设计的滤波器性能取决于窗函数的类型及窗口长度N的取值。设计进程中，要按照对阻带最小衰减和过渡带宽度的要求选择适合的窗函数类型和窗口长度N 。各类类

9、型的窗函数可达到的阻带最小衰减和过渡带宽度见下表1。表2-1 各类窗函数的大体参数窗函数旁瓣峰值幅度/dB过渡带宽阻带最小衰减/dB矩形窗-134/N-12汉明窗-418/N-53布莱克曼窗-5712/N-74如此选定窗函数类型和长度N以后，求出单位脉冲响应，并依照式（6）求出。是不是知足要求，要进行演算。一般在尾部加零使长度知足2的整数次幂，以便用FFT计算。若是要观察细节，补零点数增多即可。若是不知足要求，则要从头选择窗函数类型和长度N ，再次验算，直至知足要求。若是要求线性相位特性，则还必需知足 (14)按照上式中的正、负号和长度N的奇偶性又将线性相位FIR滤波器分成四类。要按照所设计的

10、滤波特性正确选择其中一类，例如，要设计线性相位低通特性，能够选择这一类，而不能选择这一类。咱们在设计滤波器时，希望窗谱主瓣尽可能窄，以取得较陡的过渡带，同时尽可能减小最大旁瓣的相对幅度，咱们发觉这是不可能的，因为这本身就是一个矛盾体，所以在设计滤波器时只是按如实际情形来选择适合的窗函数。3 几种数字低通滤波器的窗函数设计 采用矩形窗设计FIR数字低通滤波器图3-1 采用矩形窗设计FIR数字低通滤波器的仿真图 采用汉明窗设计FIR数字低通滤波器 图3-2 采用汉明窗设计FIR数字低通滤波器的仿真图 采用布莱克曼窗设计FIR数字低通滤波器图3-3 采用布莱克曼窗设计FIR数字低通滤波器的仿真图图3

11、-4采用布莱克曼窗冲激响应仿真图图3-5采用布莱克曼窗对含噪语音信号及滤波后的语音信号仿真图参考文献1 从玉良.数字信号处置原理及其MATLAB实现M.北京：电子工业出版社.2 胡广书.数字信号处置理论、算法与实现M.北京：清华大学出版社.2003,8 3 万永革.数字信号处置的MATLAB实现M.北京: 科学出版社.2007. 4 薛山.MATLAB基础教程M.北京: 清华大学出版社.2011. 5 陈怀琛.数字信号处置教程MATLAB释义与实现M.北京: 电子工业出版社.2002 6 缪家鼎，徐文娟 光电技术M 杭州：浙江大学出版社，1994：112-156 7 李晓东，张庆红，叶瑾琳气候

12、学研究的若干理论问题J北京大学学报：自然科学版，1999，35(1)：101-106 8 郑开青通信系统模拟及软件D北京：清华大学无线电系，1987 9 Online Computer Library Center, Inc. History of OCLCEB/OL2000-01-08 org/about/ history/default. htm. 附 录矩形窗设计程序passrad=*pi;w2=boxcar(16); n=1:1:16;hd=sin(passrad*(n-8)./(pi*(n-8);hd(8)=passrad/pi;h2=hd.*rot90(w2);title(desi

13、gned by Hanning window);mag2,rad=freqz(h2);subplot(2,1,1);plot(rad,20*log10(abs(mag2);grid on;h2,w2=freqz(h2,1,100,2);subplot(2,1,2);plot(w2,unwrap(angle(h2);grid on;汉明窗设计程序passrad=*pi;w2=hamming(32);n=1:1:32;hd=sin(passrad*(n-16)./(pi*(n-16);hd(16)=passrad/pi;h2=hd.*rot90(w2);title(designed by Hann

14、ing window);mag2,rad=freqz(h2);subplot(2,1,1);plot(rad,20*log10(abs(mag2);grid on;h2,w2=freqz(h2,1,100,2);subplot(2,1,2);plot(w2,unwrap(angle(h2);grid on;布莱克曼窗函数设计程序passrad=*pi;w2=blackman(64);n=1:1:64;hd=sin(passrad*(n-32)./(pi*(n-32);hd(32)=passrad/pi;h2=hd.*rot90(w2);title(designed by Hanning win

15、dow);mag2,rad=freqz(h2);subplot(2,1,1);plot(rad,20*log10(abs(mag2);grid on; h2,w2=freqz(h2,1,100,2);subplot(2,1,2);plot(w2,unwrap(angle(h2);grid on;用布莱克曼语音信号进行滤波，并得滤波前后信号的时域波形及频谱程序clear;clc;x,fs=wavread(D:);sound(x,fs);if x/2=0 x=x 0;endfnoise =10000;T=length(x)/fs;t=T/length(x):T/length(x):T;for j=

16、1:length(x)noise(j)=sin(2*pi*fnoise*t(j);endx=x+noise;fp=1000;fst=1200;delta_w=2*pi*(fst-fp)/fs;wc=pi*(fst+fp)/fs;N=31;if N/2=0N=N+1;endtau=(N-1)/2;for n=1:Nh(n)=sin(wc*(n-tau)/(pi*(n-tau);endh(N-1)/2)=(fst+fp)/fs;wn=(blackman(N);h1=h(n).*wn;subplot(2,1,1);stem(h1);F=abs(fft(h1);subplot(2,1,2);plot(F(1:(N-1)/2);X=abs(fft(x);figure;subplot(2,1,1);plot(X(1:(length(x)+1)/2);y=conv(x,h1);Y=abs(fft(y);subplot(2,1,2);plot(Y(1:(length(x)+1)/2);sound(x,fs);sound(y,fs);