Matlab7.1环境下IIR滤波器设计和仿真.docVIP

Matlab7.1环境下IIR滤波器设计和仿真 　　摘要：在Matlab7.1环境下IIR数字滤波器的设计有三种方法，典型设计法，完全设计法，FDATOOL设计法。由于典型设计法与完全设计法具有复杂的设计过程、大量的计算工作、调整滤波器的特性困难的特点，限制了这两种方法的应用。该文重点介绍第三种快速有效的设计数字滤波器的方法——FDATOOL设计法。并详细说明了运用Matlab语言进行完全设计法和利用FDATOOL进行IIR数字滤波器界面设计的步骤。FDATOOL界面设计法试我们可以直观简便的随时对比设计要求和调整滤波器特性参数，对于减轻设计工作量和最优化滤波器设计起到极大的作用。同时还说明了在Matlab环境下如何利用simulink对所设计的IIR数字滤波器仿真。 关键词：Matlab7.1；IIR数字滤波器；数字信号处理；无限冲激响应 中图分类号：TP391.9 文献标识码：A 文章编号：1009-3044（2013）19-4499-04 IIR数字滤波器属于经典数字滤波器中的一种，应用非常广泛。然而，在实际使用中，数字滤波器的设计要进行复杂的推导和计算， 所以运算量非常庞大。现在，常通过计算机，利用Matlab来实现数字滤波器的仿真与设计。 Matlab7.1是美国MATHWORK公司推出的一套面向科学和数值计算的可视化语言，它集数值分析、矩阵运算、信号处理和图形显示于一体，是一个高度集成系统，具有友好的用户界面和良好的帮助功能。该文是在 Matlab7.1环境，利用其自带的信号处理工具箱（ Signal Processing Toolbox），其具有强大的信号处理和分析功能，其中自带了许多函数， 可以快速有效地实现数字滤波器的设计与仿真[1]。有高效 、简便及可视化程度高的特点 ，并具有存储滤波器参数 向量的功能，可分析显示滤波器系统的幅频 、相频特性，以及对给定信号经所设计的滤波器进行滤波并分析表示了滤波前和滤波后的幅频 、相频特性。 1 IIR数字滤波器及典型设计方法 1.1 数字滤波器 数字滤波器是数字信号处理的重要基础，在对信号的过滤、检测与参数估计等处理过程中，它是使用最为广泛的一种系统。数字滤波器是指完成信号滤波处理功能的、用有限精度算法实现的离散时间线性非时变系统，其输入是一组由模拟信号取样和量化的数字量，其输出是经过数字变换的另一组数字量。数字滤波器具有稳定性高、精度高、灵活性大等突出优点。可将数字滤波器分为两种，即无限长冲激响应（IIR）滤波器和有限长冲激响应（FIR）滤波器。IIR 滤波器是一类LTI 系统，这类系统有无限持续时间的冲击响应[2]。 1.2 IIR数字滤波器设计的基本原理 滤波器的设计实质上就是去寻找一个既能物理实现，又能满足给定频率特性指标要求的系统传输函数。反馈系统中，这一时刻的输出，将会反馈作为下一时刻输入的一部分，由于输出抽样值与之前的输出值有关，从计算的观点出发，这样的滤波器也称为递归滤波器。对于N阶IIR数字滤波器，其系统函数可表示为： 1.3 IIR数字滤波器典型设计法 IIR数字滤波器具有无限宽的冲击响应， 与模拟滤波器相匹配。根据IIR滤波器特性，可通过设计模拟滤波器，再在此基础上进行变换的方法来进行设计，即为IIR数字滤波器典型设计法。其设计步骤为： 1）转换技术指标。即根据一定规则将已知的数字滤波器的技术指标转换为模拟低通滤波器的技术指标； 2）确定最小阶数及固有频率。即根据第一步的结果，通过滤波器阶数选择函数，将固有频率Wn以及最小阶数N确定。 3）生成模拟滤波器系统函数。通过上一步所得到的最小阶数，得到H（s）（模拟滤波器系统函数）。 4）将模拟滤波器的系统函数H（s），从S平面转到z平面 ，得到数字低通滤波器的系统函数H（z）。 5）将模拟滤波器转换成数字滤波器。可通过双线性变换法、冲击响应不变法来实现。 2 IIR数字滤波器完全设计法 Matlab提供了几个直接设计IIR 数字滤波器的专用函数如：butter， chebyl， cheby2， ellip和 bessel。设计时可以直接调用这些函数就能设计出低通、高通、带通和带阻滤波器[3]。 完全设计法实例：设计一个滤波器滤除指定正弦信号中的噪声，还原正弦信号。该正弦信号为sin（2*pi*100*t），噪声信号是sin（2*pi*50*t）和sin（2*pi*150*t）。信号采样频率为2KHz。滤波器的性能指标：巴特沃斯IIR带通滤波器Fs=2000Hz， Wp1=90Hz，Wp2=110Hz，Ws1=80Hz，Ws2=120Hz，Rp=1，Rs=20。 以下为实例的完全设计法程