【风电控制】二阶IIR滤波器的数学原理与负序处理应用：从理论到工程实践

引言

在电力系统、通信工程等领域，滤波器是信号处理的核心工具，用于分离特定频率分量、抑制噪声或消除干扰。其中，二阶无限脉冲响应（IIR）滤波器因结构简单、频率选择性灵活，被广泛应用于负序分量提取等场景。本文结合滤波器模型，从数学原理出发，解析二阶IIR滤波器的工作机制，并说明其在负序处理中的具体应用。

一、二阶IIR滤波器的数学基础

二阶IIR滤波器的核心是差分方程，它描述了输入信号与输出信号的时域递推关系；通过Z变换可将其转化为频域的传递函数，便于分析频率响应。

1. 差分方程（时域模型）

二阶IIR滤波器的输入输出关系由以下差分方程定义：
[ y(n) = b_0 x(n) + b_1 x(n-1) + b_2 x(n-2) - a_1 y(n-1) - a_2 y(n-2) ]

• ( x(n) )：当前时刻输入信号；
• ( y(n) )：当前时刻输出信号；
• ( x(n-1), x(n-2) )：输入的过去值（延迟1、2个采样周期）；
• ( y(n-1), y(n-2) )：输出的过去值（反馈，延迟1、2个采样周期）；
• ( b_0, b_1, b_2 )：前馈系数（决定输入对输出的贡献）；
• ( a_1, a_2 )：反馈系数（决定输出对自身的反馈，影响极点位置）。
  差分方程的物理意义是：当前输出由当前输入、过去输入和过去输出共同决定，体现了“前馈+反馈”的结构。

2. 传递函数（频域模型）

对差分方程做Z变换（离散时间域的拉普拉斯变换），得到传递函数：
[ H(z) = \frac{Y(z)}{X(z)} = \frac{b_0 + b_1 z^{-1} + b_2 z^{-2}}{1 + a_1 z^{-1} + a_2 z^{-2}} ]

• 分子：前馈路径的Z变换（对应输入的过去值，决定零点位置）；
• 分母：反馈路径的Z变换（对应输出的过去值，决定极点位置）；
• ( z^{-1} )：延迟单元（表示信号延迟1个采样周期）。
  传递函数是频域分析的核心，通过它可计算滤波器的频率响应（幅频特性和相频特性）。

二、频率响应与滤波特性

滤波器的频率选择性由**极点（分母的根）和零点（分子的根）**的位置决定。将 ( z = e^{j\omega T_s} )（( T_s ) 为采样周期，( \omega ) 为角频率）代入传递函数，得到频率响应 ( H(e^{j\omega T_s}) )，其幅频特性决定了滤波器的类型：

1. 极点零点与滤波类型

• 低通滤波器：极点靠近 ( z=1 )（直流，( \omega=0 )），允许低频信号通过，抑制高频噪声；
• 高通滤波器：极点靠近 ( z=-1 )（奈奎斯特频率，( \omega=\pi/T_s )），允许高频信号通过，抑制低频信号；
• 带通滤波器：极点在复平面的共轭位置（如 ( z = re^{\pm j\theta} )，( r < 1 )），允许某一频段（如中心频率 ( \omega_0 )）通过，抑制其他频段。

2. 二阶带通滤波器的设计

在负序处理中，二阶滤波器通常设计为带通滤波器，其中心频率设为2倍电网频率（2ω）（如电网频率50Hz时，2ω=100Hz）。此时：

• 正序分量（直流，ω=0）：被带通滤波器抑制（不在通带内）；
• 负序分量（2ω交流）：在带通滤波器的通带内，被保留并输出；
• 其他频率（如噪声、谐波）：被抑制（不在通带内）。

三、负序处理的实际应用：从Park变换到滤波器提取

在电力系统中，负序分量是三相不平衡的产物，其特性决定了滤波器的设计目标。

1. 负序分量的Park变换特性

通过Park变换（以电网同步角 ( \theta ) 为旋转角度），三相电压的分量在d-q旋转坐标系中表现为：

• 正序分量：直流（与旋转坐标系同速，相对静止）；
• 负序分量：2倍电网频率（2ω）的交流分量（与旋转坐标系反向旋转，相对速度为2ω）。

2. 二阶滤波器的提取作用

二阶带通滤波器的设计目标是提取负序分量，具体逻辑如下：

• 输入：d-q坐标系下的电压信号（包含正序直流和负序2ω交流）；
• 滤波器：中心频率设为2ω，带宽足够窄以抑制正序直流；
• 输出：负序2ω交流分量（即负序分量）。

四、设计与实现要点

1. 系数选择（稳定性与带宽）

• 稳定性：极点必须在单位圆内（( |z| < 1 )），否则滤波器不稳定；
• 带宽：通过调整 ( a_1, a_2 ) 控制极点位置，从而调整带通滤波器的带宽（如窄带宽可提高负序提取的精度）。

2. 延迟单元的作用

图中“( z^{-1} )”“( z^{-2} )”模块表示延迟单元，其作用是存储信号的过去值，实现差分方程中的“( x(n-1) )”“( y(n-1) )”等项。延迟单元的数量决定了滤波器的阶数（二阶滤波器需要2个延迟单元）。

五、总结

二阶IIR滤波器的数学原理是差分方程+Z变换，其频率选择性由极点零点位置决定。在负序处理中，通过设计带通滤波器（中心频率2ω），可提取d-q坐标系下的负序分量，解决三相不平衡问题。这种从理论到工程的应用，体现了数学工具在电力系统中的核心价值。
封面图：展示了二阶IIR滤波器的结构、差分方程、Z变换公式、频率响应曲线和Park变换示意图，直观呈现了滤波器的数学原理与应用场景。
（注：本文中的公式和模型均基于离散时间信号处理理论，实际工程中需结合采样频率、电网参数等调整系数。）