IIR数字滤波器设计：原理、MATLAB实现与选型指南

1. 数字滤波器设计概述

在信号处理领域，IIR（无限脉冲响应）滤波器因其高效的频率选择特性而广泛应用。与FIR滤波器相比，IIR滤波器能够在满足相同性能指标的情况下使用更低的阶数，这使得它在实时处理系统中具有明显优势。本文将深入解析Butterworth、Chebyshev I型、Chebyshev II型和椭圆滤波器这四种经典IIR滤波器，以及Yulewalk这种基于频率响应拟合的设计方法。

提示：IIR滤波器设计本质上是在频域特性与实现复杂度之间寻找平衡，不同设计方法各有其适用场景。

2. 滤波器理论基础与特性对比

2.1 滤波器核心参数解析

设计数字滤波器时，我们需要明确几个关键指标：

• 通带截止频率（ωp）：信号通过滤波器时幅度衰减不超过允许值的最高频率
• 阻带起始频率（ωs）：信号被显著衰减的最低频率
• 通带波纹（Rp）：通带内允许的最大波动（dB）
• 阻带衰减（Rs）：阻带要求的最小衰减量（dB）

这些参数直接影响滤波器的阶数计算和性能表现。例如，一个低通滤波器的典型频率响应曲线中，通带频率范围是0到ωp，过渡带是ωp到ωs，阻带则是ωs到π。

2.2 五种设计方法原理对比

2.2.1 Butterworth滤波器

Butterworth滤波器在通带和阻带都具有最大平坦特性，没有波纹。其幅度响应满足：
|H(ω)|² = 1 / [1 + (ω/ωc)^(2N)]
其中N是滤波器阶数，ωc是截止频率。这种设计在过渡带相对较宽，需要较高阶数才能达到陡峭的衰减。

2.2.2 Chebyshev I型

Chebyshev I型在通带具有等波纹特性，阻带单调衰减。其幅度响应为：
|H(ω)|² = 1 / [1 + ε²T_N²(ω/ωp)]
其中T_N是N阶Chebyshev多项式，ε决定通带波纹大小。相比Butterworth，它在相同阶数下能获得更陡的过渡带。

2.2.3 Chebyshev II型

Chebyshev II型（又称逆Chebyshev）在阻带呈现等波纹，通带单调。其设计公式为：
|H(ω)|² = 1 / [1 + 1/(ε²T_N²(ωs/ω))]
这种设计避免了通带波纹，适合对通带平坦度要求高的应用。

2.2.4 椭圆滤波器

椭圆滤