FIR滤波器的MATLAB设计和实现.docVIP

FIR滤波器MATLAB设计和实现 1、设计任务 1、用MATLAB软件实现FIR滤波器； 2、设计基于DSPFIR滤波器硬件框图； 3、了解用DSP实现FIR滤波器关键问题； 4、完成必需软件步骤图。 2．序言 在通信和电子信息当中，在对信号作分析和处理时，常会碰到有用信号叠加无用噪声问题。这些噪声信号有是和信号同时产生，有是在传输过程中混入，在接收信号中，必需消除或减弱噪声干扰，这是信号处理中十分关键问题。依据有用信号和噪声不一样特征，消除或减弱噪声，提取有用信号过程就称为滤波。滤波器种类很多，实现方法也多个多样。伴随数字技术飞速发展，数字滤波理论也得到了长足进步。所以，对数字滤波系统硬件实现要求也越来越高，而软件模拟方法不仅能立即地提供系统运行信息，还能够随时改变系统结构从而验证全新系统，所以软件仿真变得很关键。Matlab 是含有很强科学计算和图形显示功效软件系统，能够对数字滤波器进行正确设计，而且方便地进行FFT 频谱分析和频谱图显示，从而对数字滤波器进行快速地检验和分析。本文讨论在MATLAB 平台下FIR 数字滤波器设计和分析。 DSP是一个实时、快速、尤其适合于实现多种数字信号处理运算微处理器。因为它由含有丰富硬件资源、高速数据处理能力和强大指令系统，而在通信、航空、航天、雷达、工业控制、网络及家用电器等各个领域得到广泛应用。DSP分为定点和浮点两种，本文以定点DSP芯片为例，讨论FIR滤波器实现多个关键问题。所讨论这些问题，在DSP系统设计中有实际参考和应用价值。 3．FIR 滤波器原理和MATLAB仿真设计 3.1 滤波器概述 数字滤波在数字信号处理中占相关键地位，是广泛使用一个基础线性处理模块，它能够实现模拟器件极难达成正确线性相位关系特征。数字滤波器分为无限冲激响应滤波器(IIR)和有限冲激响应滤波器(FIR)。因为FIR系统只有零点，所以这类滤波器不像IIR 滤波器那样轻易取得比很好通带和阻带衰减特征。要取得好衰减特征，通常要求H(z)阶次要高，即N 要大。FIR 滤波器有自己突出优点，其一是系统总是稳定，其二是易实现线性相位，其三是只要经过一定时延，任何非因果有现场序列全部能变成因果有限长序列，所以总能用因果系统来实现，其四是FIR滤波器因为单位冲击响应是有限长，所以能够用快速傅里叶变换算法来实现过滤信号，从而能够大大提升运算效率。因为FIR 滤波器在阶数相同条件下运算速度比IIR 滤波器快，同时FIR滤波器含有不含反馈环路、结构简单和能够实现严格线性相位等优点，所以在对相位要求比较严格条件下，常常采取FIR 数字滤波器。现在常见设计方法关键有窗函数法、频率取样法及等波纹迫近法。本文应用窗函数法设计FIR 数字低通滤波器。具体参数见表一。 3.2 FIR 数字滤波器基础原理 3.2.1窗函数法基础思想 先结构一个线性相位理想滤波器频率响应，然后用一个N点窗函数w(n), 去截取理想滤波器单位抽样响应（通常为无限长），从而得到含有线性相位实际滤波器有限长单位抽样响。 3.2.2基础方法 (1)结构线性相位理想滤波器频率响应，为简单起见，若没有90°相移特殊要求，通常选择滤波器含有第一类线性相位。 (2)求理想滤波器单位抽样响应。 (3)依据技术指标要求和4种形式线性相位FIR数字滤波器特点，选择适宜窗函数w(n)及其长度N，然后对加窗函数截取，得到含有线性对应实际FIR数字滤波器单位抽样对应，。 (4)检验实际滤波器频率响应是否满足设计指标要求。 3.2.3用窗函数设计FIR滤波器步骤 (1)依据技术要求确定线性相位理想滤波器频率响应。 (2)求理想滤波器单位抽样响应。 (3)依据对过渡带及阻带衰减要求，选择窗函数形式，并估量窗口长度N，设待求滤波器过渡带用△表示，它近似等于窗函数主瓣宽度。 (4)计算滤波器单位抽样响应 (5)验算技术指标是否满足要求，设计出滤波器频率响应用下式计算 3.2.4用窗函数设计FIR滤波器 设滤波器通带截止频率为，阻带截止频率为，系统采样频率为，则其过渡带宽带为 ， 抽样周期为 则过渡带数字角频率为 多种窗函数过渡带宽可用来表示，其中N 是滤波器阶数，x是一个整数，对应于矩形窗，三角窗，汉宁窗，汉明窗，其值分别为4,8,8,8，所以滤波器阶数N可由下式求得 理想低通滤波器冲击响应为 其中, 是3dB通带截止数字角频率，其值为 加窗后低通滤波器冲击响应为 检验实际滤波器频率响应是否满足设计指标要求。 表1 多个常见窗函数对比 窗函数 旁瓣峰值 主瓣 最小阻带 衰减/dB 宽度/ /N 衰减/dB 矩形窗 13 4 -21 三角窗 -25 8 -