iir和fir数字滤波器的设计及其结构研究.docVIP

IIR和FIR数字滤波器的设计及其结构研究 1. 摘要 数字滤波器是一个离散系统，其系统函数一般可以表示为的有利多项式形式，即 （1-1) 当{;i=1,2,...,N}都为零时，由式1-1描述的系统为有限脉冲响应数字滤波器，简称FIR数字滤波器。当系数{;i=1,2,...,N}中至少有一个为非零时，式1-1描述的系统称为无限脉冲响应数字滤波器，简称IIR数字滤波器。对于数字滤波器，一般满足，这是系统称为N阶IIR数字滤波器。对于FIR数字滤波器，系统函数中的有理多项式的最高次幂就是其阶数. 2. 数字滤波器设计的意义 随着数字集成电路，设备和系统技术的快速进步，通过数字方法进行信号处 理已变得越来越有吸引力。目前主要有两类滤波器，模拟滤波器和数字滤波 器，它们在物理组成和工作方式上完全不同，而模拟滤波器的技术发展已相当成 熟，与模拟滤波器相比，数字滤波器是DSP(数字信号处理)系统独特而又重要的一类，是通过计算算法将输入数字序列转换为不同输出序列的离散时间系统，具有更高的精确度和可靠性，使用灵活、方便，已经成为数字信号处理技术中的重要手段。如频谱分析，数字图像处理和语音处理等等。目前随着计算机技术和数字信号处理器芯片的发展，研究不同数字滤波器的设计原理和稳定性分析对于满足军事、航空、民营等等各个领域的信号处理要求具有十分重要的意义。 3. IIR数字滤波器设计的基本过程 IIR滤波器的设计就是根据给定的数字滤波器指标确定滤波器的阶数N和系数﹛﹜。在满足技术指标的条件下，滤波器的阶数应尽可能的低，从而降低成本。我们在设计IIR数字滤波器时一般是通过模拟滤波器来设计数字滤波器。其中IIR数字滤波器的设计过程如图3.1。 待设计数字滤波器指标 待设计数字滤波器指标 模拟滤波器指标数字滤波器H(z)模拟滤波器H(s) 频率 设计模拟 s到z域 模拟 滤波器指标 数字滤波器H(z) 模拟滤波器H(s) 转换 滤波器 转换 图3.1 IIR数字滤波器的设计过程 3.1 模拟滤波器的设计及方法 常用的IIR滤波器设计是以模拟滤波器设计为基础。在设计模拟滤波器时，先将待设计的模拟滤波器技术指标转换为模拟低通滤波器技术指标，然后设计模拟低通滤波器，再通过频率变换将模拟低通滤波器转换为所需要的滤波器。因此模拟低通滤波器的设计是模拟滤波器设计的基础。 3.1.1 butterworth模拟低通滤波器 butterworth模拟低通滤波器简称BW型低通滤波器，其幅度相应的模仿为： 式中N为滤波器阶数，为滤波器的3截止频率， 当时，butterworth模拟低通滤波器称为归一化的butterworth模拟低通滤波器。设计步骤为： （1）由滤波器的设计指标和公式确定滤波器的阶数N。 （2）由公式确定。 （3）由公式 ,计算S左半平面的N个极点。 （4）由公式确定滤波器的系统函数H(s)。 3.1.2 chebyshev模拟低通滤波器 butterworth模拟低通滤波器通带和阻带内都存在许多裕量。chebyshev模拟低通滤波器的幅度响应在一个频带中具有等波纹特性可以降低滤波器的阶数。Chebyshev I型模拟低通滤波器（简称CB I）的幅度响应在通带是等波纹的，在阻带单调下降。Chebyshev = 2 \* ROMAN II型模拟低通滤波器（简称CB = 2 \* ROMAN II）的幅度响应在通带单调下降，在阻带是等波纹的。 CB I 型模拟低通滤波器其幅度响应为：，其中N是滤波器的阶数，是滤波器的参数，是N阶Chebyshev多项式。设计步骤如下： （1）由通带频率确定： =。 （2）由通带衰减确定： 。 （3）确定后，根据阻带衰减由幅度响应可求出N的阶数为： (4)由模方求滤波器的极点。 ， 其中： (5)由极点确定系统函数H(s)： CB = 2 \* ROMAN II型低通滤波器其幅度相应模方为： 其中N是滤波器阶数，是滤波器的参数。设计的具体步骤如下： （1）由通带频率确定： =。 （2）由通带衰减确定： 。 （3）确定后，根据阻带衰减由幅度响应可求出N的阶数为： （4）由表查获CB = 2 \* ROMAN II型低通滤波器的系统函数H(s) CB I 型和CB = 2 \* ROMAN II型低通滤波器最大的区别是CB I 型模拟低通滤波器在通带等波纹波动，参数控制通带波动，而CB = 2 \* ROMAN II型低通滤波器在阻带等波纹波动，参数控制阻带波动