基于MATLAB环境下FIR滤波器的设计和仿真技术.docVIP

基于MATLAB环境下FIR滤波器的设计和仿真技术 摘要：数字滤波技术在数字信号处理中占有重要地位，信号处理的核心就是滤波器。针对FIR滤波器的原理，有人提出了FIR滤波器的窗函数设计法，本文对运用最为频繁的几种窗函数而进行了几点比较。根据仿真结果来看，FIR滤波器的各方面指标都可以达到设计要求，而且设计过程非常简便。这种方法也为快速而高效地设计FIR滤波器多提供了一个渠道。 关键词：窗函数；FIR滤波器；MATLAB 在信息数据采集系统中，滤波所占有的地位非常重要，而数字滤波也因为它本身具有的设计灵活、实现简单的特点而广泛被人们所关注和接受，它也因此而成为了数据处理和数据分析的基本方法。数字滤波器的主要组成部分包括乘法器、加法器和延时单元，最终组合成为一种算法或者装置。数字滤波器的主要功能是对输入离散信号的数字代码进行运算和处理，以此来完成改变信号频谱和符合技术要求的目标。 一、FIR数字滤波器和硬件的实现方式 在数字通信、语音信号处理、图像处理和自适应处理等多个领域一般都会对信号的传输过程有要求，即不能出现明显的相位失真情况。虽然IIR的运用是存在一定的频率色散问题，但是FIR滤波器对IIR滤波器的使用中可以显示出很多比较独特的优越性，也就是在它确保能够满足滤波器幅频响应要求的同时还可以获得相对严格的线性相位特性，也能够因此而保持稳定。所以FIR滤波器也能够得到更加广泛的使用。 1、FIR数字滤波器的原理 数字滤波器的函数通常是表示为 H(Z)== （1） 根据（1）中的H(Z)可以得出表示输入和输出关系的一个常系数线性差分方程： y(n)=+ (2) 根据上面两个方程可以看出来，数字滤波器主要是把输入序列经过一些运算，如式2，而得出来的一个输出序列。但是多数一般的数字滤波器是线性而非时变的滤波器，所以在对应FIR系统中它的函数是只有零点（除了z=0的极点外），而系数全部都是0，所以公式（2）就能简化成为 y(n)= (3) 公式（3）可以理解为x(n)与单位脉冲响应h(n)的直接卷积。 2、根据滤波器设计和分析工具而来的设计技术简介 MATLAB在信号处理工具箱中有一个先进的可视化滤波器集成设计环境，也就是滤波器的设计和分析工具（Filter Design and Analysis Tool ,FDATool）。要设计滤波器的性能指标和直接赋值滤波器或者是从MATLAB的工作空间内导入滤波器参数都可以通过FDATOOL的交互式形式用户界面的多种方式来完成滤波器设计。同时FDATOOL还能够给出一组滤波器分析工具用来频率响应、零—极点图、冲激响应、群延迟等分析数据。在安装了信号处理的工具箱之后就可以通过FDATOOl而取代滤波器去设计函数，从而进行高效的滤波器分析和设计。 FDATOOL的功能是多方面的，但是主要表现在以下几个方面：1）滤波器在选择类型的时候，可以选择低通、高通、带通、带阻或者其它等多种选择类型；2）滤波器主要可以选择IIR和FIR这两种设计类型；3）能够灵活设定滤波器的指标；4）对滤波器的系数能进行快速的计算和分析；5）它能够转换滤波器的结构；6）可以导入或者倒出滤波器设计。 二、窗函数的设计 窗函数的设计一般都是根据窗函数的具体用途来确定设计方案的，窗函数主要是通过截断或者集权的方式而使得无限长非因果序列便成为有限长脉冲响应序列，在设计之前都会有具体工程用途的参数给出，以给出的标准再来具体设计滤波器。滤波器在实际中是一种选频手段，所以给出的指标也都是在频域中以分贝值形式出现的相对幅度相应和相位响应。 如何使用窗函数设计滤波器： 1、在选择窗函数的类型的时候一般都是要依据过渡带宽和阻带衰减估算窗口的长度N或者阶数M=N-1。在选择窗函数类型的时候可以以最小阻带衰减数值单独选择，主要是因为窗口的长度对最小阻带衰减值As没有影响。 窗函数类型明确后，根据之前拟定的窗函数数值来假设，设要求的滤波器过渡带宽是△，△与窗口长度成反比，确定窗口函数后，计算公式也随着确定。在计算窗口长度的过程中要保证在阻带衰减满足要求的情况下，选择较小的窗口长度N值。N和窗函数类型都确定，就可以通过MATLAB而求出窗函数值。 2、滤波器一般在设计情况下都是根据理想的频率响应而得出来的单位脉冲响应值，在这里假设要求的频率响应值是，由此可以得出理想状态下脉冲响应值的公式： = 一般值都是用数值来表示，因为通常都不采用闭合的公式去表示，从=0到=2π所采用的N点，都可以运用以上公式求出来。 3、计算滤波器中的单位脉冲响应值。滤波器最理想的单位脉冲响应与窗函数乘积：，在MATLAB运算环境中用点乘命令来表示就是。 3、这一步主要是要求出滤波器单位脉冲响应值。设计滤波器通常情况最好的值就是单位脉冲响应和窗函数之