风电控制中陷波器技术的原理与实现

1. 风电控制中的陷波器技术概述

在风力发电机组控制系统中，陷波器(Notch Filter)是一种专门用于抑制特定频率干扰的关键滤波器。这种滤波器在风电领域有着不可替代的作用，因为它能精准地消除传动链扭振、叶片通过频率等周期性干扰。与普通滤波器不同，陷波器的独特之处在于它能在保持其他频率信号基本不变的情况下，对特定窄带频率信号进行深度衰减。

陷波器的实现方式主要分为模拟电路和数字算法两种。随着现代风力发电机组普遍采用全数字化控制系统，基于差分方程的数字陷波器已成为行业主流。这种数字实现方式不仅参数调整灵活，还能方便地与其他控制算法集成。在时域中，数字陷波器通过差分方程来描述其输入输出关系，而差分方程的系数直接决定了滤波器的中心频率、带宽和衰减深度等关键特性。

2. 陷波器差分方程的基本形式

2.1 二阶IIR陷波器的标准方程

典型的数字陷波器通常采用二阶无限脉冲响应(IIR)结构，其差分方程的一般形式为：

y[n] = b0·x[n] + b1·x[n-1] + b2·x[n-2] - a1·y[n-1] - a2·y[n-2]

其中x[n]和y[n]分别表示当前时刻的输入和输出信号，系数b0、b1、b2、a1、a2决定了滤波器的频率响应特性。这种结构计算量小、实现简单，非常适合在风电控制器的实时系统中运行。

2.2 系数与频率参数的映射关系

陷波器系数与其频率特性之间存在精确的数学关系。设采样周期为Ts，需要抑制的频率为f0，带宽为Δf，则系数可以通过以下步骤计算：

1. 计算数字角频率：
   ω0 = 2π·f0·Ts
   Δω = 2π·Δf·Ts

2. 计算极点半径r（决定陷波深度）：
   r ≈ 1 - Δω/2 （通常取0.9~0.99）

3. 最终系数表达式：
   b0 = (1 + r²)/2
   b1 = -2cos(ω0)·(1 + r²)/2
   b2 = (1 + r²)/2
   a1 = -2r·cos(ω0)
   a2 = r²

3. 时域差分方程系数求解实例

3.1 示例参数设定

假设某2MW风力发电机组需要抑制传动链的扭振频率，已知：

• 采样频率：1kHz (Ts=0.001s)
• 目标抑制频率：8Hz
• 陷波带宽：1Hz
• 陷波深度：-