用matlab信号处理工具箱进行fir滤波器设计三种方法.pdfVIP

摘 要 介绍了利用 MATLAB 信号处理工具箱进行 FIR 滤波器设计的三种方法：程序设计法 FDATool 设计法和 SPTool 设计法，给出了详细的设计步骤，并将设计的滤波器应用到一个 混和正弦波信号，以验证滤波器的性能。 关键词 MATLAB ，数字滤波器，有限冲激响应，窗函数，仿真 1 前言 数字滤波器是一种用来过滤时间离散信号的数字系统，通过对抽样数据进行数学处理来 达到频域滤波的目的。根据其单位冲激响应函数的时域特性可分为两类：无限冲激响应 （IIR）滤波器和有限冲激响应（FIR）滤波器。与IIR 滤波器相比，FIR 的实现是非递归的， 总是稳定的；更重要的是，FIR 滤波器在满足幅频响应要求的同时，可以获得严格的线性相 位特性。因此，它在高保真的信号处理，如数字音频、图像处理、数据传输、生物医学等领 域得到广泛应用。 2 FIR 滤波器的窗函数设计法 FIR 滤波器的设计方法有许多种，如窗函数设计法、频率采样设计法和最优化设计法等。 窗函数设计法的基本原理是用一定宽度窗函数截取无限脉冲响应序列获得有限长的脉冲响 应序列，主要设计步骤为： (1) 通过傅里叶逆变换获得理想滤波器的单位脉冲响应 hd(n)。 (2) 由性能指标确定窗函数W(n)和窗口长度 N。 (3) 求得实际滤波器的单位脉冲响应 h(n)， h(n) 即为所设计 FIR 滤波器系数向量 b(n)。 (4) 检验滤波器性能。 本文将针对一个含有 5Hz 15Hz 和 30Hz 的混和正弦波信号，设计一个 FIR 带通滤波器， 给出利用 MATLAB 实现的三种方法：程序设计法 FDATool 设计法和 SPTool 设计法。参 数要求：采样频率 fs=100Hz ，通带下限截止频率 fc1=10 Hz ，通带上限截止频率 fc2=20 Hz，过渡带宽 6 Hz，通阻带波动 0.01，采用凯塞窗设计。 2 程序设计法 MATLAB 信号处理工具箱提供了各种窗函数、滤波器设计函数和滤波器实现函数。本文 的带通滤波器设计及滤波程序如下： [n,Wn,beta,ftype]=kaiserord([7 13 17 23],[0 1 0],[0.01 0.01 0.01],100); %得出滤波器的阶数 n=38，beta=3.4 w1=2*fc1/fs; w2=2*fc2/fs;%将模拟滤波器的技术指标转换为数字滤波器的技术指标 window=kaiser(n+1,beta);%使用 kaiser 窗函数 b=fir1(n,[w1 w2],window);使用标准频率响应的加窗设计函数 fir1 freqz(b,1,512);%数字滤波器频率响应 t = (0:100)/Fs; s = sin(2*pi*t*5)+sin(2*pi*t*15)+sin(2*pi*t*30);%混和正弦波信号 sf = filter(b,1,s); ％对信号s 进行滤波 程序执行的结果如图 1 所示： （1）滤波器幅频特性和相频特性 （2）滤波前后的波形 图 1 滤波器特性和滤波效果图 4 FDATool 设计法 FDATool(Filter Design Analysis Tool)是 MATLAB 信号处理工具箱专用的滤波器设计分 析工具，操作简单、灵活，可以采用多种方法设计 FIR 和 IIR 滤波器。在 MATLAB 命令窗口 输入 FDATool 后回车就会弹出 FDATool 界面。 4.1 带通滤波器设计 已知滤波器的阶数 n=38，beta=3.4。本例中，首先在 Filter Type 中选择 Bandpass；在 Design Method 选项中选择 FIR Window，接着在 Window 选项中选取 Kaiser，Beta 值为 3.4；指定 Filter Order 项中的 Specify order 为 38；采样频率 Fs=100Hz,截止频率 Fc1=10Hz,Fc2=20Hz。设置完以后点击窗口下方的 Design Filter，在窗口上方就会看到所设 计滤波器的幅频响应，通过菜单选项 Analysis 还可以看到滤波器的相频响应、组延迟、脉冲 响应、阶跃响应、零极点配置等。设计完成后将结果保存为 kaiser15.fda 文件。 4.2