1. IIR滤波器设计核心思路解析

在数字信号处理领域，IIR（无限脉冲响应）滤波器因其高效的频率选择特性被广泛应用。与FIR滤波器相比，IIR滤波器在相同性能指标下通常需要更低的阶数，这使得它在实时处理场景中具有明显优势。但非线性相位特性也带来了设计上的特殊考量。

1.1 IIR与FIR的本质差异

IIR滤波器的递归结构决定了其脉冲响应在理论上是无限长的，这源于系统函数分母多项式的存在。从实现角度看，这种结构带来两个关键特征：

• 存储单元需求少（n阶滤波器只需n个延迟单元）
• 能够用较低阶数实现陡峭的过渡带

但递归结构也引入了稳定性问题，所有极点必须位于单位圆内。这点在MATLAB实现时尤其需要注意，使用zp2sos转换为二阶节形式能有效提升数值稳定性。

经验提示：实际工程中建议优先使用[z,p,k]或sos格式，避免直接使用传递函数系数[b,a]，可减少量化误差带来的影响

2. 经典设计方法全解

2.1 模拟原型法设计流程

最传统的IIR设计方法基于模拟滤波器转换，核心步骤包括：

1. 参数规格化：

   • 将数字频率要求转换为模拟角频率
   • 对低通原型进行频率预畸变

2. 原型选择：

   matlab复制% 巴特沃斯原型示例
   [z,p,k] = buttap(n);
   [b,a] = zp2tf(z,p,k);
   

3. 频率变换：

   matlab复制% 低通到带通变换
   [bt,at] = lp2bp(b,a,w0,bw);
   

4. 离散化处理：

   matlab复制% 双线性变换
   [bz,az] = bilinear(bt,at,fs);
   

2.2 直接设计法实践

对于特殊响应需求，yulewalk提供了直接拟合频响曲线的方案：

matlab复制f = [0 0.3 0.4 0.6 0.7 1]; % 归一化频率
m = [1 1 0 0 1 1];         % 期望增益
[b,a] = yulewalk(10,f,m);   % 10阶设计

这种方法特别适合：

• 多通带/阻带需求
• 非标准形状的频响曲线
• 快速原型验证