基于ＭＡＴＬＡＢ的ＩＩＲ滤波器设计与仿真.docVIP

基于ＭＡＴＬＡＢ的ＩＩＲ滤波器设计与仿真 　　摘要：介绍了利用MATLAB快速有效地设计无限长冲激响应数字滤波器的方法。给出了使用MATLAB 语言进行程序设计和FDATool工具进行界面设计的详细步骤。介绍了如何利用MATLAB环境下的仿真软件Simulink对所设计的滤波器进行模拟仿真。 　　关键词：数字滤波器；MATLAB；无限长冲激响应 　　中图分类号：TN713 TN702文献标识码： A文章编号：1009-3044(2008)23-1072-03 　　Design and Simulation of IIR Filter Based onMATLAB 　　WANG Zhan-zhong 　　(Department of Computer Science and Information Engineering, Anyang Institute of Technology, Anyang 455000, China) 　　Abstract: MATLAB-based methods are introduced to design the infinite impulse response (IIR) filter effectively. To design the IIR filter , it can be programmed by MATLAB language or FDATool tool. Detailed steps of theseways are given. With the Simulink of MATLAB, the simulation for designated filter is also introduced. 　　Key words: digital filter;MATLAB;IIR 　　 　　1 引言 　　 　　数字滤波器(Digital Filter,DF)是数字信号处理的重要内容，在对信号的过滤、检测与参数的估计等信号处理中，数字滤波器是使用最为广泛的装置，无论在工业、农业和其它行业均有应用。数字滤波器实质上是一个有限精度算法实现的线性时不变离散系统，它的基本工作原理是利用离散系统特性对系统输入信号进行加工和变换，改变输入序列的频谱或信号波形，让有用的频率信号分量通过，抑制无用的信号分量输出[1]。数字滤波器按照其冲激响应函数的时域特性，可分为无限长冲激响应(Infinite Impulse Response,IIR)滤波器和有限长冲激响应(Finite Impulse Response,FIR)滤波器。在满足相同指标下，IIR滤波器的阶数明显小于FIR,硬件实现容易且大大减少了运算量，遗憾的是IIR是非线性相位,在不要求严格线性相位的情况下，IIR滤波器的应用相当广泛。 　　本文基于 MATLAB7.1 讨论 IIR 数字滤波器的设计问题。在 MATLAB 里提供了很多设计数字滤波器的方法，这里只介绍两种：通过程序设计和通过 FDATool 工具设计。通过一个具体的实例，说明详细的操作，最后对实现的滤波器进行仿真。 　　 　　2 常规IIR滤波器的基本设计法 　　 　　IIR滤波器设计的最基本的方法就是模拟原型法。数字滤波器和模拟滤波器有千丝万缕的联系，它们之间的转换是s平面和z平面的转换，转换的基本方式就是冲激响应不变法和双线性变换法。根据设计指标用合适的传递函数去实现之，对于要求特殊的滤波器，自我构建传递函数是一个非常复杂的工作，不过经过努力可以得到极佳的效果（当然MATLAB也提供从指定的频幅特性的直接设计，在此不作讨论）。对于常规的应用而言，MATLAB 提供有许多经典的模拟原型函数候选[2-3]。设计人员要做的工作是选取适宜的函数，指定合适的参数，最后对得到滤波器进行验证。 　　选择原型函数时的考虑如下：巴特沃斯滤波器的频率特性在通带和阻带内都是随着频率单调的变化，显然，如果在通带的边缘能满足指标，在通带的内部肯定超过设计指标要求，造成滤波器的阶数比较高，在各频带内没有幅度的波动，如果要求幅值稳定性极好时，巴特沃斯滤波器是一个不错的选择。在许多情况下，降低滤波器的阶数至为重要，因为降低阶数意味着实现简单、成本低，所以一般设计人员都比较重视降低阶数。如果将指标的精度要求均匀地分布在整个通带内，或者均匀地分布在整个阻带内，更有效的方法是同时均匀地分布在通带和阻带内，可以设计出满足设计要求的、阶数又比较低的的滤波器。切比雪夫I型滤波器特性是频幅在通带内是等波纹的，在阻带内是单调的；切比雪夫II型则相反，在通带内是单调的，在阻带内是等波纹的。所以一般地讲，切比雪夫滤波器的阶数比巴特沃斯滤波器要低些。还有一种值得关注的滤波器是椭圆