IIR 滤波器设.docVIP

第6章 IIR数字滤波器的设计 数字滤波器是数字信号处理技术的重要内容。和模拟滤波器一样，数字滤波器的主要功能是对数字信号进行处理，最常见的处理是保留数字信号中的有用频率成分，去除信号中的无用频率成分。本章首先介绍数字滤波器的原理和分类，使大家对数字滤波器有个整体的认识，然后介绍由模拟滤波器传递函数转换为数字滤波器传递函数的方法----脉冲响应不变法和双线性变换法。在进行MATLAB设计IIR数字滤波器的讲解之前，第三节介绍滤波器特性分析及使用方法。最后大篇幅举例介绍IIR数字滤波器的经典设计、完全设计函数和直接设计方法。 6.1 概 述 6.1.1 数字滤波器的工作原理 数字滤波器是一个离散时间系统，输入x(n)是一个时间序列，输出y(n)也是一个时间序列。如数字滤波器的系统函数为H(z)，其脉冲响应为h(n)，则在时间域内存在下列的关系。 (6-1) 在z域内，输入和输出存在下列关系： (6-2) 式中，X(z)、Y(z)分别为输入x(n)和输出y(n)的z变换。 同样在频率域内，输入和输出存在下列关系： (6-3) 式中，为数字滤波器的频率特性; 和分别为x(n)和y(n)的频谱。为数字角频率，单位rad。通常设计在某些频段的响应值为1，在某些频段的响应为0。和的乘积在频率响应为1的那些频段的值仍为，即在这些频段的振动可以无阻碍地通过滤波器，这些频带为通带。和的乘积在频率响应为0的那些频段的值不管大小如何均为零，即在这些频段里的振动不能通过滤波器，这些频带称为阻带。 为讨论方便起见，我们在本章中规定，为数字角频率，单位为弧度(rad)，表示模拟角频率，单位弧度/秒（rad/s）。数字角频率在0～范围内。 一个合适的数字滤波器系统函数H(z)可以根据需要改变输入x(n)的频率特性。经数字滤波器处理后的信号y(n)保留信号x(n)中的有用频率成分，去除无用频率成分。 6.1.2 滤波器的分类： 按时间域特性，数字滤波器可以分为无限冲激（脉冲）响应数字滤波器（Infinite impulse response digital filter,简称IIR滤波器）和有限冲激（脉冲）响应数字滤波器（Finite impulse response digital filter,简称FIR滤波器）两类。 IIR滤波器的传递函数为： (6-4) h(n)为滤波器的脉冲响应，n=0~均有值。M和N为分解的分子和分母多项式的系数个数。FIR滤波器的传递函数为： (6-5) 该滤波器的脉冲响应h(n)在n=0,1,…,N-1的有限个点(N个点)上有值。注意（6-4）式中分母全为零时，H(z)具有全零点形式，IIR滤波器退化为FIR滤波器。 按频率特性来讲，数字滤波器和模拟滤波器一样可分为低通、高通、带通和带阻等。数字滤波器是一个离散时间系统，在频率特性中具有周期性，因此我们讨论的频率范围仅在的范围内，相应的归一化频率在0~1之间，和1对应于Nyquist频率。 和模拟滤波器一样，理想数字滤波器的频率特性在通带内必须满足： (6-6) 式中，K,均为常数。 和模拟滤波器一样，数字滤波器的设计目的是使滤波器的频率特性达到所给定的性能指标。其性能指标也包括通带波纹Rp(dB)、阻带衰减Rs(dB)、通带边界频率（Hz）、阻带边界频率（Hz）等。 6.1.3 IIR数字滤波器的设计方法 IIR滤波器的设计就是根据滤波器的性能指标要求，设计滤波器的分子和分母多项式系数。它和FIR滤波器相比优点是在满足相同性能指标要求的前提下，IIR滤波器的阶数要明显低于FIR滤波器；但IIR滤波器的相位是非线性的。 在第5章我们学习了模拟滤波器的设计。在数字滤波器设计中，首先让我们想到的是能否利用模拟滤波器的设计成果进行数字滤波器设计，答案是肯定的。现在我们所讲的IIR经典设计就是将已设计好的模拟滤波器按一定变换原理转换为数字滤波器。该方法先根据滤波器的技术指标设计出相应的模拟滤波器，然后再将设计好的模拟滤波器变换成数字滤波器。在MATLAB中，经典法设计IIR数字滤波器主要采用以下步骤：