基于汉宁窗的FIR数字滤波器精品资料

1、（1）熟悉用窗函数法设计 FIR 数字滤波器的原理与方法；课（2）了解用等波纹最佳逼近法设计 FIR 数字滤波器的原理与方法；程（3）掌握调用 MATLAB 信号处理工具箱中滤波器设计函数（或滤波器设设计分析工具 fdatool）设计各种 FIR 数字滤波器，学会根据滤波需求确计定滤波器指标参数。目（4）通过观察滤波器输入输出信号的时域波形及其频谱，加深对数字滤的波器的概念理解。（1）简述线性相位 FIR 滤波器的特点及条件；课（2）简述窗函数法 FIR 滤波器的设计步骤和等波纹最佳逼近法的设计思程想；设（3）解释吉布斯效应，比较各种典型窗函数的性能特点；计（4）完成以上设计实验，并对结果进行

2、分析和解释；要（5）打印程序清单和要求画出的信号波形；求（6）写出本次课程设计的收获和体会。课程（1）采样间隔采样间隔过大或过小都会产生较大的误差和频谱混叠现设象，因此采样频率通常选为（3 4） fh(fh为信号的最高频率 ) 。计（2）截断、泄露和窗函数为了减小截断的影响，常采用其他的时窗函注意数来对所截取的时域信号进行加权处理。所选的窗函数应力求其频谱的事主瓣宽度窄些，旁瓣幅度小些。项课用汉宁窗函数法设计一个数字FIR 带通滤波器，要求通带边界频率程 为 400Hz，500Hz，阻带边界频率为 350Hz，550Hz，通带最大衰减 1dB，设 阻带最小衰减 40dB，抽样频率为 2000H

3、z，用 MATLAB 画出幅频特性，计画出并分析滤波器系统函数的零极点；内信 号 x(t ) x1 (t )x2 (t ) sin(2 f1t ) sin(2 f 2容f 1 450Hz， f 2600Hz，滤波器的输出 y(t)t) 经 过 该 滤 波 器 ， 其 中1.选题：我选的课题是基于汉宁窗的数字滤波器设计；课2.了解数字滤波器和汉宁窗的设计原理；程3. 掌握了解课程设计的设计要求和设计方法，开发及设计工具的使设计用方法；简4.根据指标要求拟定设计方案， 对设计方案进行分析选择最优设计方要案进行设计；操5.在 Matlab 软件环境下编写仿真程序，运行观察仿真结果进行分析作得出结论；

4、步骤6.整理清单和目录，书写设计报告。通过两周的基于汉宁窗的数字滤波器设计， 对于 MATLAB 语句有了更加深刻的理解，也注意到了一些运算符号的使用，掌握了课一些课程设计的设计要求和设计方法，开发及设计工具的使用方法，程最重要的是通过这一设计实践的过程，我们不再约束在理论上，而且设计锻炼了动手能力和分析解决问题的能力，积累了经验，培养了按部就心班，一丝不苟的态度和对所学知识的综合应用能力，而且通过同学间得的分组合作课题，锻炼了我们的团队合作能力。体会这次的数字信号处理实训，让我受益匪浅，不只是老师和书本带来的知识与收获，还可以通过我们的自我实践来完成目标设计，是对我们很好的锻炼，既巩固知识，

5、又提高了实践动手能力。课程评设语计评语及成指导教师成绩绩（签 名）年月日目录目录3一摘要4二原理说明52.1 数字滤波技术52.2 FIR 滤波器6三窗函数简介：汉宁窗7四程序设计8五结果与分析115.1 汉宁窗函数仿真结果115.2 汉宁窗函数仿真结果分析13六总结与体会13七 . 致谢14八参考文献14一摘要数字滤波器是一种用来过滤时间离散信号的数字系统， 通过对抽样数据进行数学处理来达到频域滤波的目的。 根据其单位冲激响应函数的时域特性可分为两类：无限冲激响应（ IIR ）滤波器和有限冲激响应（ FIR）滤波器。与 IIR 滤波器相比， FIR 的实现是非递归的，总是稳定的；更重要的是，

6、 FIR 滤波器在满足幅频响应要求的同时， 可以获得严格的线性相位特性。 因此，它在高保真的信号处理，如数字音频、图像处理、数据传输、生物医学等领域得到广泛应用。滤波器的设计是信号处理的核心问题之一。 根据 FIR 滤波器的原理，提出了 FIR 滤波器的窗函数设计法，给出了在 MATLAB 环境下，用窗函数法设计 FIR 滤波器的过程和设计实例。 通过利用汉宁窗函数方法设计 FIR 滤波器，对所设计的滤波器进行分析比较， 得出设计的滤波器的优缺点及其不同的使用场合， 从而可以在设计滤波器时能够正确的选择 FIR 数字滤波器的窗函数的选取及设计方法。关键词： MATLAB汉宁窗二原理说明随着信息

7、时代的到来，数字信号处理已经成为一门极其重要的学科和技术，并且在通信、语音、图像、自动控制等众多领域得到了广泛的应用。在数字信号处理中，数字滤波器占有极其重要的地位，它具有精度高、可靠性好、灵活性大等特点。现代数字滤波器可以用软件或硬件两种方式来实现。软件方式实现的优点是可以通过滤波器参数的改变去调整滤波器的性能。MATLAB是一种面向科学和工程计算的语言，它集数值分析、矩阵运算、信号处理和图形显示于一体，具有编程效率高、 调试手段丰富、 扩充能力强等特点。 MATLAB的信号处理工具箱具有强大的函数功能，它不仅可以用来设计数字滤波器，还可以使设计达到最优化，是数字滤波器设计的强有力工具。2.

8、1数字滤波技术数字滤波，就是通过一定的计算或判断程序减少干扰在有用信号中的比重，故实质上是一种程序滤波。 与此对应的就是模拟滤波， 由于模拟滤波牵扯到的其他知识太多在此不详细介绍了，模拟滤波主要无源绿波（直接用电阻、电容、电感等不外接电源的元件组成的）与有源滤波（如运算放大器等需要外接电源组成的），其目的是将信号中的噪音和干扰滤去或者将希望得到的频率信号滤出为我所用。数字滤波的出现克服了模拟滤波的很多不足，具有以下优点：A 是用程序实现的，不需要增加硬设备，所以可靠性高，稳定性好。B可以对频率很低的信号实现滤波，克服了模拟滤波的缺陷。C可以根据信号的不同，采用不同的滤波方法或参数，具有灵活、方

9、便、功能强的特点。几种常用的滤波方法：1 算术平均值法2 中值滤波法3 滑动平均值法4 限幅滤波法5 惯性滤波法数字滤波技术通过数字滤波器实现， 从实现方法上可以分为FIR 数字滤波器和 IIR 数字滤波器，按功能可分为低通滤波器（ LPF）、高通滤波器（ HPF）、带通滤波器（ BPF）和带阻滤波器（ BSF）。本文主要对 FIR 滤波器加以介绍。2.2 FIR滤波器FIR(Finite Impulse Response)滤波器：有限长单位冲激响应滤波器，是数字信号处理系统中最基本的元件， 它可以在保证任意幅频特性的同时具有严格的线性相频特性， 同时其单位抽样响应是有限长的， 因而滤波器是稳

10、定的系统。 因此， FIR 滤波器在通信、图像处理、模式识别等领域都有着广泛的应用。FIR 滤波器的工作原理：在进入FIR滤波器前，首先要将信号通过A/D 器件进行模数转换， 使之成为 8bit 的数字信号，一般可用速度较高的逐次逼进式A/D 转换器，不论采用乘累加方法还是分布式算法设计 FIR 滤波器，滤波器输出的数据都是一串序列，要使它能直观地反应出来，还需经过数模转换，因此由 FPGA构成的 FIR 滤波器的输出须外接 D/A 模块。FPGA有着规整的内部逻辑阵列和丰富的连线资源， 特别适合于数字信号处理任务， 相对于串行运算为主导的通用 DSP芯片来说，其并行性和可扩展性更好， 利用

11、FPGA 乘累加的快速算法，可以设计出高速的 FIR 数字滤波器。窗函数设计法是FIR 滤波器的一种基本设计方法，它的基本思路是直接从立项滤波器的频率特性入手，通过积分求出对应的单位采样响应表达式，最后通过加窗，得到满足要求的FIR滤波器的单位采样响应， 窗函数在很大程度上决定了FIR 滤波器的性能指标， 因此称作“窗函数设计法” 。三窗函数简介：汉宁窗窗函数法是设计FIR 滤波器的最主要方法之一， 实际中遇到的离散时间信号总是有限长的， 因此不可避免的要遇到数据截短的问题，在信号处理中， 对离散序列的截短是通过序列与窗函数相乘来实现的。在信号处理中， 窗函数是一种除在给定区间之外取值均为0

12、的实函数。 譬如：在给定区间内为常数而在区间外为0 的窗函数被形象地称为矩形窗。 任何函数与窗函数之积仍为窗函数，所以相乘的结果就像透过窗口“看”其他函数一样。窗函数在光谱分析、滤波器设计以及音频数据压缩等方面有广泛的应用。汉宁窗（ Hanning Window）又称升余弦窗，汉宁窗可以看作是3 个矩形时间窗的频谱之和，或者说是 3 个 sint谱窗向左、右各移动了 /T ，从而使旁瓣互相抵消，消去高频干扰和漏能。可以看出，汉宁窗主瓣加宽并降低，旁瓣则显著减小，从减小泄漏观。汉宁窗（ Hanning）函数时域形式可表示为：W n11cos 2 nRN ( n)2N 1利用傅利叶变换的调制特性，

13、由上式可得汉宁窗的平谱函数为：22N 1jW e j0.5WR0.25 WRWRe2N1N1N 1jWe2jN 1式中， WR e jWRe2当 N 远大于 1 时，上式可近似表示为：W0.5WR0.25 WR2WR2NN这三部分之和使旁瓣互相抵消，能量更集中在主瓣， 汉宁窗函数的最大旁瓣值比主瓣值低 31dB，但是主瓣宽度比矩形窗函数的主瓣宽度增加了1倍，为8 /N。hanning 函数：生成汉宁窗调用方式：(1)w = hanning(n) ：输入参数 n 是窗函数的长度；输出参数w 是由窗函数的值组成的 n 阶向量。注意：此函数不返回是零点的窗函数的首尾两个元素。(2) w = hann

14、ing(n,'symmetric')：与上面相类似。(3) w = hanning(n,'periodic') ：此函数返回包括为零点的窗函数的首尾两个元素。四程序设计用汉宁窗函数法设计一个数字 FIR 带通滤波器，要求通带边界频率为 400Hz，500Hz，阻带边界频率为 350Hz，550Hz，通带最大衰减 1dB，阻带最小衰减 40dB，抽样频率为 2000Hz，用 MATLAB画出幅频特性，画出并分析滤波器系统函数的零极点；信 号 ()x1()x2(t)s i n2 ()s i n2 ()经过该滤波器，其中x ttf1 tf 2tf1 450Hz， f2

15、600Hz，滤波器的输出 y(t) 。hanning 窗设计程序：%用汉宁窗函数法设计一个数字FIR 带通滤波器 %clearall; clcFs=2000;fpl=400;fpu=500;fsl=350;fsu=550;wpl=2*pi*fpl/Fs;wpu=2*pi*fpu/Fs;%通带截止频率 %wsl=2*pi*fsl/Fs;wsu=2*pi*fsu/Fs;%阻带频率频率 %Bt=wpl-wsl;%计算过渡带宽度 %N0=ceil(6.2*pi/Bt);%计算所需 h（n ）长度 N0，ceil取大于等于 %wc=(wpl+wsl)/2/pi,(wpu+wsu)/2/pi;%计算理想带

16、通滤波器截止频率 %N=N0+mod(N0+1,2);%确保 h （n）长度 N是奇数 %n=0:N-1;hn=fir1(N-1,wc,hanning(N);%调用 firl计算带通 FIR 数字滤波器的h（ n） %H,w=freqz(hn,1,1024);%计算频率响应函数 h（ n）%figure(1);magH=20*log10(abs(H)/max(abs(H);%计算幅度 %subplot(3,1,1);stem(n,hn,'.');xlabel('n' );ylabel('h(n)' );title(' 汉宁窗 FIR 数字

17、带通滤波器的单位脉冲响应 ' );subplot(3,1,2);plot(w/pi*Fs/2,magH);%绘制幅度特性 %title(' 汉宁窗 FIR 数字带通滤波器的幅度特性 ' );xlabel(' 频率 /Hz' );ylabel('20lg|Hg(ejomega)|/max(|Hg(ejomega)|)');gridon ;subplot(3,1,3);plot(w/pi*Fs/2,unwrap(angle(H);%频率响应相位 %title(' 汉宁窗 FIR 数字带通滤波器的相位特性 ' );xlabel(

18、' 频率 /Hz' );ylabel(' 相位 /rad');gridon ;figure(2);zplane(hn,1);%绘制零极点图 %title(' 汉宁窗 FIR 数字带通滤波器系统函数的零极点图' );legend('零点','极点');gridon ;n=0:1023;dt=1/Fs;t=n*dt;Tp=1024*dt;f1=450;f2=600;x=sin(2*pi*f1*t)+sin(2*pi*f2*t);y=filter(hn,1,x);figure(3);subplot(2,1,1);plot

19、(t,x);title( ' 输入信号 x(t)' xlabel( 't/s' ); ylabel( 'x(t)' ););axis(0,Tp/6,min(x),max(x);subplot(2,1,2);plot(t,y);title( ' 输出信号 y(t)' xlabel( 't/s' ); ylabel( 'y(t)' ););axis(0,Tp/6,min(y),max(y);五结果与分析5.1汉宁窗函数仿真结果图 1图 2图 35.2 汉宁窗函数仿真结果分析由仿真波形图可以得出以下几个方

20、面的分析结果：由图 1，汉宁窗是典型的升余弦窗，基本符合汉宁窗函数的时域表达式，所加的窗口是余弦函数 （类似余弦函数） 第一幅图可以看出初始状态为零时单位脉冲响应的变化，大约在 n=62 处偶对称，这说明 FIR 带通滤波器相位是严格的线性关系；再由第二幅幅度特性看出在频率 400500 之间时，幅度为零，即为通带边界频率；由第三幅相位特性图看出频率在350 550 之间时相位呈衰减趋势，即为阻带边界频率也就是说旁瓣幅度远小于主瓣幅度，能量全部集中在主瓣， 主瓣宽度也有所增加，效果比较明显。由图 2 可以看出，极点为 0，零点分布在 1 和-1 附近，是互为倒数的共轭对，还有几个特殊的零点和极点同在0 上。由图 3 可以看出，汉宁窗设计的输入信号x(t) 的波形在通过滤波器的软件运行后，在通带之外的波形被过滤掉，经过大概0.03s, 在通带频率之外的频率的信号的幅度将受到很大的影响，以致信号不能通过滤波器，最后输出y(t) 波形经过滤波后可以认为是较均匀的余弦波信号，其中的不规则由吉布斯效应产生。六总结与体会本次 MATLAB 课程设计主要任务是完成 FIR 滤波器的设计，对我来说这个题目还是很有挑战性的，因为自己对 MATLAB 中有关滤波器的设计知识了解较少，但正是这种有挑战性的题目才能提高自己