[计算机]基于MATLAB的FIR滤波器设计及实现.pdfVIP

i于MATLAB的FlR滤波器 设计及实现 邵 毅 温 艳 (宿州学院 宿州 234000) 条件是其乘法器的系数关于中点对称 ，又因为其结 1 引言 构无反馈的回路 ，所 以是无条件 的稳定系统。 数字滤波是数字信号处理的一种重要算法，广 3 窗函数法设计及 MATLAB实现 泛用于对信号的过滤、检测与参数的估计等信号处 理中。数字滤波器是使用最为广泛的装置 ，在工业 、 FIR滤波器 目前常用的设计方法有窗函数法和 农业和其他行业均有应用 。数字滤波器按其单位脉 频率采样法，窗函数法是从时域进行设计，而频率采 冲响应的长度可分为有限脉冲响应 (vm)滤波器和 样法是从频域进行设计 。窗函数法 由于简单、物理意 无限脉冲响应 (IIR)滤波器两类 。FIR滤波器的线性 义清晰，因而得到了较为广泛 的应用 。窗函数法设 与稳定性使其应用更为广泛。本文主要介绍采用窗 计的基本思想是：首先根据技术指标要求，选取合适 函数法设计 FIR滤波器及其 MATLAB的软件实现 的阶数 Ⅳ和窗函数的类型 (n)，使其幅频特性逼近 方法 。 理想滤波器幅频特性 。其次 ，因为理想滤波器的 h(凡)是无限长的，所以需要对 h (n)进行截断，数学 2 FIR数字滤波器设计原理 上称这种方法为窗函数法。 FIR数字滤波器可 以用下面的差分方程来描述 一 1 。j (n)=FT H【d( )] ，v一1 (n)= h(m)x(n—m)= ()h(n) (n)= (n)w(n) m=0 其中h (n)是希望逼近的滤波器的单位脉冲响 其中 (n)是输入序列，v(n)是输出序列 ，h(n) 是系统的单位脉冲响应。由于其采用卷积的方法得 应，H ( )是希望逼近 的滤波器的幅频特性，h(n) 到输出，故可以采用快速算法FFT进行卷积计算，提 是所设计的滤波器的单位脉冲响应 ，(n)是窗函数。 高信号处理速度 。FIR数字滤波器满足线性相位的 在截断的过程中要保证 FIR滤波器的稳定性与线性 基金支持：宿州学院自然科学研究项 目(2008yzk09)；宿州学院硕士科研启动基金(2008yss14)。 83 中国教通信 20102 CH I1A N EW TELECo M M UN 1CA Tf0 N S 圈IFIR数字滤波系统 相位的特性。由于截断效应使所设计滤波器的幅频 的滤波器设计方法，并填人滤波器设计指标，即可为 特性存在误差，称之为吉布斯效应。最后，验证所设 用户生成所需的滤波器频率响应特性。根据技术指 计的滤波器是否满足要求，若不满足，重新设计。 标要求，采用凯塞窗进设计一个低通的34阶FIR数 常用的窗函数有矩形窗、三角窗、汉宁窗和凯塞 字滤波器。通带截止频率为 40Hz，阻带截止频率为 窗等。窗函数的选择一般要满足以下三个条件：(1) lOOHz，通带最大衰减为 ldB，阻带最小衰减为 具有较低的旁瓣幅度，尤其是第一旁瓣幅度。(2)旁 80dB。 瓣幅度下降速度要快，以利于增加阻带衰减。(3)主