matlab滤波器的设计与实现,fir滤波器的matlab设计与实现.doc

fir滤波器的matlab设计与实现.doc

1FIR滤波器的MATLAB设计与实现1、设计任务1、用MATLAB软件实现FIR滤波器；2、设计基于DSP的FIR滤波器硬件框图；3、了解用DSP实现FIR滤波器的关键问题；4、完成必要的软件流程图。2．前言在通信与电子信息当中，在对信号作分析与处理时，常会遇到有用信号叠加无用噪声的问题。这些噪声信号有的是与信号同时产生的，有的是在传输过程中混入的，在接收的信号中，必须消除或减弱噪声干扰，这是信号处理中十分重要的问题。根据有用信号与噪声的不同特性，消除或减弱噪声，提取有用信号的过程就称为滤波。滤波器的种类很多，实现方法也多种多样。随着数字技术的飞速发展，数字滤波理论也得到了长足的进步。因此，对数字滤波系统硬件实现的要求也越来越高，而软件模拟的方法不仅能及时地提供系统运行的信息，还可以随时改变系统结构从而验证全新的系统，所以软件仿真变得非常重要。MATLAB是具有很强的科学计算和图形显示功能的软件系统，可以对数字滤波器进行精确设计，并且方便地进行FFT频谱分析与频谱图显示，从而对数字滤波器进行快速地检验和分析。本文讨论在MATLAB平台下的FIR数字滤波器设计与分析。DSP是一种实时、快速、特别适合于实现各种数字信号处理运算的微处理器。由于它由具有丰富的硬件资源、高速数据处理能力和强大的指令系统，而在通信、航空、航天、雷达、工业控制、网络及家用电器等各个领域得到广泛应用。DSP分为定点和浮点两种，本文以定点DSP芯片为例，讨论FIR滤波器实现的几个关键问题。所讨论的这些问题，在DSP系统设计中有实际的参考和应用价值。3．FIR滤波器的原理与MATLAB仿真设计31滤波器概述数字滤波在数字信号处理中占有重要的地位，是广泛使用的一种基本线性处理模块，它可以实现模拟器件很难达到的准确线性相位关系特性。数字滤波器分为无限冲激响应滤波器IIR和有限冲激响应滤波器FIR。由于FIR系统只有零点，因此这类滤波器不像IIR滤波器那样容易取得比较好的通带与阻带衰减特性。要2取得好的衰减特性，一般要求HZ的阶次要高，即N要大。FIR滤波器有自己突出的优点，其一是系统总是稳定的，其二是易实现线性相位，其三是只要经过一定的时延，任何非因果有现场序列都能变成因果有限长序列，因而总能用因果系统来实现，其四是FIR滤波器由于单位冲击响应是有限长的，因而可以用快速傅里叶变换算法来实现过滤信号，从而可以大大提高运算效率。由于FIR滤波器在阶数相同的条件下运算速度比IIR滤波器快，同时FIR滤波器具有不含反馈环路、结构简单以及可以实现的严格线性相位等优点，因而在对相位要求比较严格的条件下，常常采用FIR数字滤波器。目前常用的设计方法主要有窗函数法、频率取样法及等波纹逼近法。本文应用窗函数法设计FIR数字低通滤波器。具体参数见表一。32FIR数字滤波器基本原理321窗函数法的基本思想先构造一个线性相位理想滤波器的频率响应，然后用一个N点的窗函JWDHE数WN,去截取理想滤波器的单位抽样响应(通常为无限长)，01NNDHN从而得到具有线性相位的实际滤波器的有限长单位抽样响。322基本方法1构造线性相位理想滤波器的频率响应，为简单起见，若没有90°JWDHE相移的特殊要求，一般选择滤波器具有第一类线性相位。2求理想滤波器的单位抽样响应。DHN3根据技术指标要求和4种形式的线性相位FIR数字滤波器的特点，选择合适的窗函数WN及其长度N，然后对加窗函数截取，得到具有线性相应的实D际FIR数字滤波器的单位抽样相应，。HNW01NN4检验实际滤波器的频率响应是否满足设计指标要求。JHEDTFH323用窗函数设计FIR滤波器的步骤1根据技术要求确定线性相位理想滤波器的频率响应。JWDHE2求理想滤波器的单位抽样响应。DHN3根据对过渡带及阻带衰减的要求，选择窗函数的形式，并估计窗口长度3N，设待求滤波器的过渡带用△表示，它近似等于窗函数主瓣的宽度。W4计算滤波器的单位抽样响应DHN5验算技术指标是否满足要求，设计出的滤波器频率响应用下式计算10NJWJWNNHEE324用窗函数设计FIR滤波器设滤波器的通带截止频率为，阻带截止频率为，系统采样频率为，PFSFSAMPF则其过渡带宽带为，TWSPFF抽样周期为2SSAMPT则过渡带数字角频率为2SPTWTSAMPAMFF各种窗函数的过渡带宽可用来表示，其中N是滤波器阶数，X是一个整数，对X应于矩形窗，三角窗，汉宁窗，汉明窗，其值分别为4,8,8,8，所以滤波器的阶数N可由下式求得2SAMPTWXFN理想低通滤波器的冲击响应为SINCDH其中,是3DB通带截止数字角频率，其值为12NCW2PSSPSCAMFTF加窗后的低通滤波器的冲击响应为NDNHW检验实际滤波器的频率响应是否满足设计指标要求。JHEDTFH表1几种常见的窗函数对比4旁瓣峰值主瓣最小阻带窗函数衰减/DB宽度//N衰减/DB矩形窗13421三角窗25825汉宁窗31844海明窗41853例用窗函数设计一个线性相位FIR低通滤波器，并满足性能要求通带边界的归一化频率WP05,阻带边界的归一化频率WS066，阻带衰减不小于30DB,通带波纹不大于3DB，假设一个信号，其中F15HZ,F220HZ。信号的采样频率为50HZ。并将原信号与通过滤波器的信号进行比较。由题意值，阻带衰减不小于30DB,根据表1，选取汉宁窗，因为汉宁窗的第一旁瓣相对主瓣衰减为31DB，满足滤波要求。WP05PIWS066PI滤波器的边界频率WDELTAWSWP过渡带宽度NCEIL8PI/WDELTA根据过渡带宽等于表中汉宁窗函数的主瓣宽度求得滤波器所用常函数的最小带宽WN0506PI/2截止频率取通带和阻带边界频率的中点BFIR1N,WN/PI,HANNINGN1设计FIR滤波器H,FFREQZB,1,512,50采用50HZ的采样频率绘出该滤波器的幅频和相频响应SUBPLOT2,1,1,PLOTF,20LOG10ABSHXLABEL 频率 YLABEL 振幅 GRIDONSUBPLOT2,1,2,PLOTF,180/PIUNWRAPANGLEHXLABEL 频率 YLABEL 相位 GRIDONF13F220DT002T0DT3采样间隔和检测信号的时间序列XSIN2PIF1TCOS2PIF2T检测信号YFFTFILTB,X给出滤波器的输出FIGURE2SUBPLOT2,1,1,PLOTT,X,TITLE 输入信号 绘出输入信号5SUBPLOT2,1,2,PLOTT,Y绘出输出信号HOLDONPLOT11N1/2DT,YLIM, R 绘出延迟到的时刻XLABEL 时间 ,TITLE 输出信号 图1所设计滤波器的幅频响应(上图)和相频响应(下图)6图2所设计滤波器的输入和输出信号程序运行结果如图1,2该例对应于50HZ的采样频率通带边界频率为FP50/205125HZ,FS50/2066165HZ。有图1上图得，在小于125HZ的频段上，几乎看不到下降，即满足通带波纹不大于3DB的要求。在大于165HZ的频段上，阻带衰减大于30DB,满足题目要求。由图1下图得，在通带范围内，相位频率响应为一条直线，表面该