MATLAB数字滤波器设计和仿真.docVIP

MATLAB的数字滤波器设计与仿真 摘要：根据一实例提出基于MATLAB的三种设计IIR数字带通滤波器的方法，给出了MATLAB完整程序与SPTOOL、FDATOOL工具界面设计的详细步骤，并用SIMULINK进行仿真。最终分析了二种不同设计方法的区别与优势。 关键词： MATLAB；DF；IIR；SPTOOL；FDATOOL；SIMULINK2、数字滤波器 数字滤波器是一种用来过滤时间离散信号的数字系统，通过对抽样数据进行数学处理来达到频域滤波的目的。数字滤波器也是具有一定传输选择特性的数字信号处理装置，其输入、输出均为数字信号，实质上是一个由有限精度算法实现的线性时不变离散系统。它的基木工作原理是利用离散系统特性对系统输入信号进行加工和变换，改变输入序列的频谱或信号波形，让有用频率的信号分量通过，抑制无用的信号分量输出。数字滤波器除了具有数字信号处理的固有优点外，还有滤波精度高、稳定性好、灵活性强等优点。 数字滤波器有多种分类, 根据数字滤波器冲激响应的时域特征, 可将数字滤波器分为两种,即无限长冲激响应(IIR)滤波器和有限长冲激响应(FIR)滤波器；根据数字滤波器的实现方法和形式可将数字滤波器分为递归型、非递归型、快速卷积型；根据数字滤波器的频率特性可以分为低通、高通、带通、带阻数字滤波器。 3、IIR数字滤波器的设计 IIR滤波器系统函数的极点可以在单位圆内的任何位置，实现IIR滤波器的阶次较低，所用的存储单元较少，效率高，又由于IIR数字滤波器能够保留一些模拟滤波器的优良特性，因此应用很广。 一个 N阶 IIR 数字滤波器的系统函数可以表示为以下的有理分式: 图1 从上式可以看出设计IIR数字滤波器也即确定上式的系数, 使滤波器的频率特性满足给定的技术指标要求。 目前设计数字滤波器的方法可以分为两大类：人工设计、计算机辅助设计。在人工设计中，由于模拟滤波器设计理论已比较成熟，因此数字滤波器设计的关键就是将H(s)转化成为相应的H(z)即把H(s)平面映射到H(z)平面。H(s)平面映射到H(z)平面的方法主要有基于冲激响应不变法的IIR数字滤波器设计，基于双线性Z变换法的IIR数字滤波器设计。在实际工程应用中，由于手工设计计算复杂、效率低下，设计人员更多地采用计算机辅助设计。 下面以一个具体的例子来说明不同的方法。如产生有三个正弦成分（0Hz，15 Hz和30Hz）的信号，设计一滤波器来去除0Hz和30Hz的正弦信号，保留15 Hz的信号。 3.1 用FDATOOL+SIMULINK设计 IIR 数字滤波器 Matlab中提供了功能强大的SIMULINK仿真软件。SIMULINK环境下可以对所设的滤波器进行仿真。 FDATOOL（Filter Design Analysis Tool）是 MATLAB 信号处理工具箱专用的滤波器设计分析工具，操作简单、灵活，可以采用多种方法设计 FIR和IIR滤波器。在 MATLAB 的命令窗口输入FDATOOL, 按回车进入FDATOOL设计界面进行设计。在实际设计中可以采用FDATOOL+SIMULINK的设计方法，即在SIMULINK环境中调入Digital Filter Design模块，双击该模块即可打开如FDATOOL的界面。然后再调入一些信号模块和测量模块就可以仿真了。如图2所示即为仿真模型。 图2 采用FDATOOL+SIMULINK法 在图2中，双击Digital Filter Design模块进入如图3所示的界面，具体配置参数如下图所示。 图3 Digital Filter Design模块参数配置 滤波前后的波形图如图4所示，可以看出滤波后的波形含有较少的谐波，达到了滤波效果。 图4 滤波前后的波形图 对滤波前后的波形图进行频谱分析可以得出如下图5、图6所示的频谱图。 图5 滤波前的信号频谱图 图6 滤波后的信号频谱图 3.2 用SPTOOL设计 IIR 数字滤波器 SPTOOL是信号处理工具箱中一个具有交互式图形用户界面的信号处理工具，专门用于完成常用的数字信号处理任务。SPTOOL可用于信号浏览、滤波器的设计分析以及频谱分析。它主要由三个交互式信号处理模块构成：Signal Browser，主要用于数字信号的显示、分析；Filter模块，主要用于数字滤波器的幅频响应及相位响应、阶跃及脉冲响应等的查看，各种滤波器的分析及设计；Spectra模块，用于信号的频谱分析。 下面用SPTOOL设计上例的滤波器，具体方法如下： （1）创建并导入信号源。此过程中采用从工作空间中导入的方法，需预先在命令窗口中输入命令：clear all;Fs=100;t=(0:100)/Fs; s= sin(2*