《基于某汉宁窗的FIR数字滤波器》由会员分享，可在线阅读，更多相关《基于某汉宁窗的FIR数字滤波器（14页珍藏版）》请在金锄头文库上搜索。
1、word课 程 设 计 目 的1熟悉用窗函数法设计FIR数字滤波器的原理与方法；2了解用等波纹最优逼近法设计FIR数字滤波器的原理与方法；3掌握调用MATLAB信号处理工具箱中滤波器设计函数或滤波器设计分析工具fdatool设计各种FIR数字滤波器，学会根据滤波需求确定滤波器指标参数。4通过观察滤波器输入输出信号的时域波形与其频谱，加深对数字滤波器的概念理解。课 程 设 计 要 求1简述线性相位FIR滤波器的特点与条件；2简述窗函数法FIR滤波器的设计步骤和等波纹最优逼近法的设计思想；3解释吉布斯效应，比拟各种典型窗函数的性能特点；4完成以上设计实验，并对结果进展分析和解释；5打印程序清单和要
2、求画出的信号波形；6写出本次课程设计的收获和体会。课 程 设 计 注 意 事 项1采样间隔 采样间隔过大或过小都会产生较大的误差和频谱混叠现象，因此采样频率通常选为34fh(fh为信号的最高频率)。2截断、泄露和窗函数 为了减小截断的影响，常采用其他的时窗函数来对所截取的时域信号进展加权处理。所选的窗函数应力求其频谱的主瓣宽度窄些，旁瓣幅度小些。课 程 设 计 内 容用汉宁窗函数法设计一个数字FIR带通滤波器，要求通带边界频率为400Hz，500Hz，阻带边界频率为350Hz，550Hz，通带最大衰减1dB，阻带最小衰减40dB，抽样频率为2000Hz，用MATLAB画出幅频特性，画出并分析滤
3、波器系统函数的零极点；信号经过该滤波器，其中450Hz，600Hz，滤波器的输出课 程 设 计 简 要 操 作 步 骤1.选题：我选的课题是基于汉宁窗的数字滤波器设计；2.了解数字滤波器和汉宁窗的设计原理；3. 掌握了解课程设计的设计要求和设计方法，开发与设计工具的使用方法； 4.根据指标要求拟定设计方案，对设计方案进展分析选择最优设计方案进展设计； 5.在Matlab软件环境下编写仿真程序，运行观察仿真结果进展分析得出结论； 6.整理清单和目录，书写设计报告。课 程 设 计 心 得 体 会 通过两周的基于汉宁窗的数字滤波器设计，对于MATLAB语句有了更加深刻的理解，也注意到了一些运算符号的
4、使用，掌握了一些课程设计的设计要求和设计方法，开发与设计工具的使用方法，最重要的是通过这一设计实践的过程，我们不再约束在理论上，而且锻炼了动手能力和分析解决问题的能力，积累了经验，培养了按部就班，一丝不苟的态度和对所学知识的综合应用能力，而且通过同学间的分组合作课题，锻炼了我们的团队合作能力。这次的数字信号处理实训，让我受益匪浅，不只是教师和书本带来的知识与收获，还可以通过我们的自我实践来完成目标设计，是对我们很好的锻炼，既巩固知识，又提高了实践动手能力。课 程 设 计 评 语 与 成 绩评 语成 绩指导教师签 名年 月 日目录目录4一摘要5二原理说明62.1 数字滤波技术62.2 FIR滤波
5、器7三窗函数简介：汉宁窗8四程序设计9五结果与分析125.1 汉宁窗函数仿真结果1214六总结与体会14七 . 致谢15八参考文献15一摘要数字滤波器是一种用来过滤时间离散信号的数字系统，通过对抽样数据进展数学处理来达到频域滤波的目的。根据其单位冲激响应函数的时域特性可分为两类：无限冲激响应IIR滤波器和有限冲激响应FIR滤波器。与IIR滤波器相比，FIR的实现是非递归的，总是稳定的；更重要的是，FIR滤波器在满足幅频响应要求的同时，可以获得严格的线性相位特性。因此，它在高保真的信号处理，如数字音频、图像处理、数据传输、生物医学等领域得到广泛应用。滤波器的设计是信号处理的核心问题之一。根据FI
6、R滤波器的原理，提出了FIR滤波器的窗函数设计法，给出了在MATLAB环境下，用窗函数法设计FIR滤波器的过程和设计实例。通过利用汉宁窗函数方法设计FIR滤波器，对所设计的滤波器进展分析比拟，得出设计的滤波器的优缺点与其不同的使用场合，从而可以在设计滤波器时能够正确的选择FIR数字滤波器的窗函数的选取与设计方法。关键词：MATLAB 汉宁窗二原理说明随着信息时代的到来，数字信号处理已经成为一门极其重要的学科和技术，并且在通信、语音、图像、自动控制等众多领域得到了广泛的应用。在数字信号处理中，数字滤波器占有极其重要的地位，它具有精度高、可靠性好、灵活性大等特点。现代数字滤波器可以用软件或硬件两种
7、方式来实现。软件方式实现的优点是可以通过滤波器参数的改变去调整滤波器的性能。MATLAB是一种面向科学和工程计算的语言，它集数值分析、矩阵运算、信号处理和图形显示于一体，具有编程效率高、调试手段丰富、扩大能力强等特点。MATLAB的信号处理工具箱具有强大的函数功能，它不仅可以用来设计数字滤波器，还可以使设计达到最优化，是数字滤波器设计的强有力工具。2.1 数字滤波技术数字滤波，就是通过一定的计算或判断程序减少干扰在有用信号中的比重，故实质上是一种程序滤波。与此对应的就是模拟滤波，由于模拟滤波牵扯到的其他知识太多在此不详细介绍了，模拟滤波主要无源绿波直接用电阻、电容、电感等不外接电源的元件组成的
8、与有源滤波如运算放大器等需要外接电源组成的，其目的是将信号中的噪音和干扰滤去或者将希望得到的频率信号滤出为我所用。数字滤波的出现克制了模拟滤波的很多不足，具有以下优点：A是用程序实现的，不需要增加硬设备，所以可靠性高，稳定性好。B可以对频率很低的信号实现滤波，克制了模拟滤波的缺陷。C可以根据信号的不同，采用不同的滤波方法或参数，具有灵活、方便、功能强的特点。几种常用的滤波方法：1 算术平均值法2 中值滤波法3 滑动平均值法4 限幅滤波法5 惯性滤波法数字滤波技术通过数字滤波器实现，从实现方法上可以分为FIR数字滤波器和IIR数字滤波器，按功能可分为低通滤波器LPF、高通滤波器HPF、带通滤波器
9、BPF和带阻滤波器BSF。本文主要对FIR滤波器加以介绍。2.2 FIR滤波器FIR(Finite Impulse Response)滤波器：有限长单位冲激响应滤波器，是数字信号处理系统中最根本的元件，它可以在保证任意幅频特性的同时具有严格的线性相频特性，同时其单位抽样响应是有限长的，因而滤波器是稳定的系统。因此，FIR滤波器在通信、图像处理、模式识别等领域都有着广泛的应用。 FIR滤波器的工作原理：在进入FIR滤波器前，首先要将信号通过A/D器件进展模数转换，使之成为8bit的数字信号，一般可用速度较高的逐次逼进式A/D转换器，不论采用乘累加方法还是分布式算法设计FIR滤波器，滤波器输出的数
10、据都是一串序列，要使它能直观地反响出来，还需经过数模转换，因此由FPGA构成的FIR滤波器的输出须外接D/A模块。FPGA有着规整的内部逻辑阵列和丰富的连线资源，特别适合于数字信号处理任务，相对于串行运算为主导的通用DSP芯片来说，其并行性和可扩展性更好，利用FPGA乘累加的快速算法，可以设计出高速的FIR数字滤波器。窗函数设计法是FIR滤波器的一种根本设计方法，它的根本思路是直接从立项滤波器的频率特性入手，通过积分求出对应的单位采样响应表达式，最后通过加窗，得到满足要求的FIR滤波器的单位采样响应，窗函数在很大程度上决定了FIR滤波器的性能指标，因此称作“窗函数设计法。三窗函数简介：汉宁窗
11、窗函数法是设计FIR滤波器的最主要方法之一，实际中遇到的离散时间信号总是有限长的，因此不可防止的要遇到数据截短的问题，在信号处理中，对离散序列的截短是通过序列与窗函数相乘来实现的。在信号处理中，窗函数是一种除在给定区间之外取值均为0的实函数。譬如：在给定区间内为常数而在区间外为0的窗函数被形象地称为矩形窗。任何函数与窗函数之积仍为窗函数，所以相乘的结果就像透过窗口“看其他函数一样。窗函数在光谱分析、滤波器设计以与音频数据压缩等方面有广泛的应用。汉宁窗Hanning Window又称升余弦窗，汉宁窗可以看作是3个矩形时间窗的频谱之和，或者说是3个sin(t)型函数之和，而括号中的两项相对于第一个
12、谱窗向左、右各移动了/T，从而使旁瓣互相抵消，消去高频干扰和漏能。可以看出，汉宁窗主瓣加宽并降低，旁瓣如此显著减小，从减小泄漏观。 汉宁窗Hanning函数时域形式可表示为：利用傅利叶变换的调制特性，由上式可得汉宁窗的平谱函数为：式中，当N远大于1时，上式可近似表示为：这三局部之和使旁瓣互相抵消，能量更集中在主瓣，汉宁窗函数的最大旁瓣值比主瓣值低31dB，但是主瓣宽度比矩形窗函数的主瓣宽度增加了1倍，为8/N。hanning函数：生成汉宁窗调用方式：(1) w = hanning(n)：输入参数n是窗函数的长度；输出参数w是由窗函数的值组成的n阶向量。注意：此函数不返回是零点的窗函数的首尾两个
13、元素。 (2) w = hanning(n,symmetric)：与上面相类似。(3) w = hanning(n,periodic)：此函数返回包括为零点的窗函数的首尾两个元素。四程序设计用汉宁窗函数法设计一个数字FIR带通滤波器，要求通带边界频率为400Hz，500Hz，阻带边界频率为350Hz，550Hz，通带最大衰减1dB，阻带最小衰减40dB，抽样频率为2000Hz，用MATLAB画出幅频特性，画出并分析滤波器系统函数的零极点；信号经过该滤波器，其中450Hz，600Hz，滤波器的输出。hanning 窗设计程序：%用汉宁窗函数法设计一个数字FIR带通滤波器%clear all; c
14、lcFs=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);%计算所需hn长度N0，ceil取大于等于%wc=(wpl+wsl)/2/pi,(wpu+wsu)/2/pi;%计算理想带通滤波器截止频率%N=N0+mod(N0+1,2);%确保hn长度N是奇数%n=0:N-1;hn=fir1(N-1,wc,hanning(N);%调用fi
15、rl计算带通FIR数字滤波器的hn%H,w=freqz(hn,1,1024);%计算频率响应函数hn%figure(1);magH=20*log10(abs(H)/max(abs(H);%计算幅度%subplot(3,1,1);stem(n,hn,.);xlabel(n);ylabel(h(n);title(汉宁窗FIR数字带通滤波器的单位脉冲响应);subplot(3,1,2);plot(w/pi*Fs/2,magH);%绘制幅度特性%title(汉宁窗FIR数字带通滤波器的幅度特性);xlabel(频率/Hz);ylabel(20lg|Hg(ejomega)|/max(|Hg(ejomega)|);grid on;subplot(3,
