1 基本概念
1.1 定义与命名
IIR滤波器是指其输出不仅取决于当前及过去的输入信号，也与过去的输出信号有关的一类数字滤波器。由于这种反馈结构，理论上其单位脉冲响应不会在有限长度内完全结束，因此称为“无限脉冲响应”滤波器。它通常由差分方程描述，并可用递归运算实现。
在工程实践中，IIR滤波器常被用于在较低阶数下逼近所需的频率特性，因此在计算资源有限的场景中具有较高的使用价值。
1.2 与FIR滤波器的区别
与FIR滤波器相比，IIR滤波器最显著的特点是包含反馈项。FIR滤波器只依赖当前和过去输入，不涉及输出反馈，因此结构上更易保证稳定性，且相位特性通常更容易设计为线性相位。
IIR滤波器的优势在于同等幅频指标下往往所需阶数更低，计算量也较小；但其设计更复杂，对稳定性、舍入误差以及实现结构更为敏感。此外，IIR滤波器的相位一般不呈线性变化，某些对波形保真度要求较高的应用会因此受到限制。
1.3 典型应用场景
IIR滤波器广泛应用于音频均衡、通信接收机、传感器信号平滑、振动分析以及生物信号预处理等领域。在这些场景中，它常被用来实现低通、高通、带通、带阻等基本滤波功能。
由于其计算效率较高，IIR滤波器也常见于嵌入式系统和实时处理平台。例如，在心电信号去噪、无线信道选择性抑制以及控制系统中的信号整形任务中，都可以见到其应用。
2 数学基础
2.1 差分方程表示
IIR滤波器通常可表示为线性常系数差分方程。其输出由输入序列与历史输出序列共同决定，典型形式为
y[n] = Σ b_k x[n-k] - Σ a_m y[n-m]
其中，x[n]为输入，y[n]为输出，b_k和a_m为滤波器系数。该表达式体现了递归关系，也是IIR滤波器最基本的数学描述。
2.2 Z变换与传递函数
对差分方程进行Z变换后，可得到滤波器的传递函数。IIR滤波器的系统函数通常写作
H(z) = Y(z) / X(z)
其分子与分母均为z的多项式，因此可形成具有极点和零点的有理函数形式。这种表示方法便于分析频率特性、稳定性以及系统结构。
2.2.1 极点与零点
零点决定某些频率分量的抑制能力，极点则与系统的共振、增益峰值及稳定性密切相关。对于IIR滤波器而言，极点位置尤为关键，因为极点若接近单位圆，会使系统在对应频率附近具有更强的响应。
若极点落在单位圆外，则系统通常不稳定；若全部极点位于单位圆内，则系统在常见条件下可视为稳定。
2.2.2 频率响应
将传递函数中的z替换为单位圆上的复变量e^{jω}，即可得到频率响应H(e^{jω})。其幅度特性反映滤波器对不同频率成分的衰减或放大能力，相位特性则描述各频率分量的相移情况。
IIR滤波器常借助极点和零点的分布来塑造频响，因此能够用较少的阶数获得较陡峭的过渡带。
2.3 系统函数的有理表达式
IIR滤波器的系统函数通常是输入输出比值的有理表达式，分子对应零点，分母对应极点。其一般形式可以写成
H(z) = (b_0 + b_1 z^-1 + ... + b_N z^-N) / (1 + a_1 z^-1 + ... + a_M z^-M)
这种结构便于通过极点和零点的配置实现目标滤波特性，也是后续各种设计方法的基础。
3 工作原理
3.1 反馈机制
IIR滤波器的核心在于反馈。系统当前输出的一部分会重新作为后续计算的输入，从而形成递归过程。正是这种反馈，使得滤波器可以在较少系数下产生较复杂的频率响应。
反馈机制虽然提高了效率，但也使系统更容易受到参数误差影响，因此在实现时需要谨慎处理。
3.2 脉冲响应特性
当输入为单位脉冲时，IIR滤波器的输出通常不会在有限时间后完全归零，而是以指数衰减或其他递归方式延续下去。这一特征构成其“无限脉冲响应”名称的来源。
在理想稳定条件下，脉冲响应会逐渐衰减到接近零；若系统不稳定，则响应可能持续增大，甚至发散。
3.3 稳定性与因果性
对于IIR滤波器，稳定性通常与极点位置直接相关。若系统全部极点位于单位圆内，则响应对有界输入保持有界输出，通常可认为是稳定的。
因果性则表示输出仅依赖当前与过去的输入、输出，不依赖未来值。工程中常用的IIR滤波器一般都设计为因果系统，以便实时实现。
4 结构形式
4.1 直接型结构
直接型结构是IIR滤波器最常见的实现方式之一，通常按照差分方程直接搭建延迟单元、乘法器和加法器。它形式直观，便于理解与编程实现。
4.1.1 直接型I
直接型I结构将输入通道与反馈通道分开实现，分别使用对应的延迟链。其结构清晰，便于从差分方程直接映射，但所需存储单元相对较多。
4.1.2 直接型II
直接型II结构将输入与输出的延迟单元合并，通常比直接型I占用更少的存储资源。由于内部状态更紧凑，它在很多数字信号处理系统中更受欢迎，但也可能更容易受到数值误差影响。
4.2 级联型结构
级联型结构把高阶IIR滤波器分解为若干个二阶节串联而成。每个二阶节只负责一部分极点和零点的实现，有利于改善数值稳定性，并降低高阶多项式直接实现带来的误差积累。
这种结构特别适合高阶滤波器设计，也便于分段调试和参数调整。
4.3 并联型结构
并联型结构将整体传递函数分解为若干个简单子系统的和，再将各支路输出叠加。它适合处理可分解为部分分式的系统函数，并有助于分析不同极点簇对整体响应的贡献。
在某些情况下，并联结构还能提高实现灵活性，便于针对不同频段分别设计子滤波器。
4.4 格型结构
格型结构以递推关系和反射系数为基础，常见于自适应滤波和参数估计等场景。其特点是结构规整、数值性质较好，且对系数量化有一定鲁棒性。
在需要稳定递推和便于在线更新参数的应用中，格型结构具有较高实用价值。
5 设计方法
5.1 模拟滤波器原型法
模拟滤波器原型法是IIR设计中常用的一类思路。其基本做法是先从经典模拟滤波器原型出发，再通过离散化变换得到数字滤波器。
这种方法能够继承成熟的模拟滤波理论，便于获得平滑、稳定且指标明确的数字滤波器。
5.1.1 巴特沃斯设计
巴特沃斯滤波器以幅频响应单调平滑著称，通带内没有波纹，整体过渡较为自然。它常用于对通带平坦性要求较高的场合。
相较于其他原型，巴特沃斯滤波器的滚降速度通常不算最陡，但设计和解释都较为直观。
5.1.2 切比雪夫设计
切比雪夫滤波器通过允许通带或阻带出现波纹，换取更陡的过渡带。根据波纹位置不同，可分为不同类型，适用于对选择性要求较强的场景。
这类设计在相同阶数下往往比巴特沃斯获得更强的频率分离能力。
5.1.3 椭圆设计
椭圆滤波器在通带和阻带都允许波纹，因此通常能在给定阶数下实现最陡峭的幅频过渡。它的频率选择性很强，适合对滤波边界要求严格的应用。
不过，椭圆滤波器的相位特性通常也更复杂，工程实现时需权衡其优缺点。
5.2 双线性变换法
双线性变换法是将模拟域s平面映射到数字域z平面的一种常用方法。它可以将模拟滤波器原型无失真地映射为数字滤波器，并避免频率混叠问题。
该方法的一个重要特点是频率预畸变，即在设计阶段对关键频率进行修正，以补偿频率映射带来的非线性失真。
5.3 冲激响应不变法
冲激响应不变法通过使数字滤波器的冲激响应在采样点上与模拟原型一致来进行设计。它在低频段常能较好保留模拟滤波器特性。
但该方法可能引入频谱重叠，因此更适合带宽较窄或对高频折叠不敏感的情形。
5.4 数字频率变换
数字频率变换是在数字域内直接将一种基本滤波器变换为另一种类型的方法，例如由低通原型推导高通、带通或带阻滤波器。它可在保持结构规律的基础上扩展滤波器用途。
这种方法使得已设计好的滤波器原型能够更方便地服务于多种频段需求。
6 性能分析
6.1 幅频特性
幅频特性是评价滤波器最直观的指标之一，反映其对不同频率成分的增益变化。对于IIR滤波器而言，幅频曲线通常能较快从通带过渡到阻带，这也是其高效性的体现。
实际分析时，常关注通带纹波、阻带衰减和过渡带宽度等参数。
6.2 相频特性
相频特性描述滤波器对不同频率分量所施加的相移。IIR滤波器一般不具备线性相位，因此不同频率的相移往往不一致。
在音频、测量和波形重建等应用中，相位畸变可能影响信号形态，需要结合具体需求进行评估。
6.3 群时延
群时延表示信号包络或窄带分量通过滤波器时所经历的延迟变化。它通常由相位对频率的导数决定，用于衡量滤波器对时间结构的影响。
对于IIR滤波器，群时延常在过渡带附近变化较大，这可能导致波形展宽或失真。
6.4 稳定裕度
稳定裕度用于描述系统距离失稳状态的安全范围。对于IIR滤波器而言，极点若远离单位圆，则通常具有更好的稳定裕度；若极点靠近单位圆，则虽能获得更尖锐的频响，但系统对误差更敏感。
在工程设计中，稳定裕度与频率选择性之间往往需要折中。
7 实现与优化
7.1 定点与浮点实现
IIR滤波器可以在定点或浮点平台上实现。浮点实现精度高、范围大，适合一般软件处理；定点实现则更节省硬件资源，常见于嵌入式和专用数字电路。
在定点环境下，系数量化和溢出处理尤为重要，需要通过缩放和结构优化加以控制。
7.2 量化误差
量化误差来源于系数、状态变量或中间计算结果的有限精度表示。对于IIR滤波器，这类误差可能引起极点偏移，进而改变实际频率响应，严重时甚至导致不稳定。
因此，在高精度要求的应用中，系数量化策略和字长选择都是关键问题。
7.3 数值稳定性
数值稳定性不仅关乎理论稳定，还涉及有限精度运算下的实际表现。某些在数学上稳定的IIR结构，在定点实现中也可能因舍入误差放大而出现异常。
采用二阶节级联、适当排序极点和零点、避免内部状态过大等措施，通常有助于提升数值稳定性。
7.4 运算效率优化
IIR滤波器的效率优化主要包括减少乘法次数、合并延迟单元、利用二阶节结构以及针对硬件特性进行并行化处理等。对于实时系统，运算效率直接影响延迟和功耗。
在实际部署中，通常需要在精度、复杂度和资源占用之间取得平衡。
8 应用领域
8.1 音频信号处理
在音频处理中，IIR滤波器常用于均衡器、低音增强、噪声抑制和音色修饰等任务。由于其结构紧凑，适合实时音频设备和嵌入式音频芯片。
不过，在需要严格线性相位的场合，设计者往往会更谨慎地选择滤波方案。
8.2 通信系统
通信系统中，IIR滤波器可用于信道选频、载波抑制、基带整形以及接收端噪声过滤等环节。其较低计算量有助于在高速数据处理链路中降低负担。
在调制解调、同步与链路前端处理中，它常作为基础模块之一。
8.3 生物医学信号处理
在心电、脑电和肌电等信号处理中，IIR滤波器常用于去除工频干扰、基线漂移和高频噪声。其计算效率适合长期监测和便携设备。
但由于生物信号对波形保真较敏感，设计时通常需要控制相位失真和瞬态响应。
8.4 控制与自动化系统
在控制与自动化领域，IIR滤波器可用于传感器数据平滑、反馈信号预处理以及执行器驱动信号整形。它能够改善测量噪声对控制回路的影响。
在实时控制场景中，滤波器的稳定性与延迟性能直接关系到系统响应质量。
9 相关概念
9.1 滤波器阶数
滤波器阶数通常指系统中延迟项或极点的数量，反映其复杂程度。阶数越高，通常可获得更陡的频率过渡，但计算量和数值风险也随之增加。
9.2 截止频率
截止频率是通带与阻带之间的重要分界参数，用于描述滤波器开始明显衰减信号的频率位置。它是设计和评价滤波器时的核心指标之一。
9.3 通带与阻带
通带是希望信号尽量保留的频段，阻带则是希望被抑制的频段。两者之间的过渡区域构成滤波器选择性的体现，也是幅频特性分析的重点。
9.4 采样定理
采样定理说明了连续信号离散化时所需的最低采样条件，是数字滤波器设计的前提之一。若采样率不足，可能产生混叠，影响IIR滤波器的实际效果。