风电控制中二阶巴特沃斯滤波器的设计与应用

1. 风电控制中的滤波需求解析

在风力发电机组控制系统中，转速信号和功率信号的测量噪声是影响控制性能的关键因素。风速的随机波动会导致发电机转速信号中包含高频噪声成分，这些噪声若直接进入控制系统，可能引发控制器的高频振荡。我在某2MW直驱永磁机组调试过程中就遇到过这样的案例：原始转速信号中的2kHz开关噪声导致变桨系统出现5°的高频抖动。

巴特沃斯滤波器因其在通带内具有最大平坦的幅度特性，成为风电信号处理的首选方案。与切比雪夫滤波器相比，它没有通带纹波；与贝塞尔滤波器相比，它具有更陡峭的过渡带。特别是在变桨控制这类对相位线性度要求不苛刻的场合，二阶巴特沃斯滤波器能有效滤除500Hz以上的高频噪声，同时将关键频段（0-10Hz）的幅值衰减控制在±3dB以内。

2. 二阶巴特沃斯滤波器的数学本质

2.1 归一化低通原型

标准二阶巴特沃斯滤波器的传递函数为：

code复制H(s) = 1 / (s² + √2 s + 1)

这个看似简单的表达式蕴含着三个重要特性：

1. 极点位置在单位圆上等角度分布（本例为±45°）
2. -3dB截止频率严格等于1 rad/s
3. 分子常数项确保直流增益为1

在实际风电应用中，我们需要将其转换为截止频率ωc可调的通用形式。通过频率缩放变换s → s/ωc，得到：

code复制H(s) = ωc² / (s² + √2 ωc s + ωc²)

这个变换相当于将滤波器特性曲线沿频率轴拉伸或压缩。例如当ωc=100π rad/s（对应50Hz）时，滤波器-3dB点就移动到50Hz处。

2.2 极点的几何意义

将分母多项式进行因式分解：

code复制s² + √2 ωc s + ωc² = (s - p₁)(s - p₂)

解得极点位置：

code复制p₁, p₂ = (-√2/2 ± j√2/2)ωc

在复平面上，这两个极点与虚轴的夹角为45°，距原点距离为ωc。这种对称分布正是巴特沃斯滤波器相位响应非线性的根源。我在MATLAB中仿真发现，当输入信号频率达到ωc时，单级二阶滤波器会产生约90°的相位滞后，这对需要快速响应的变距控制是个挑战。

3. 从模拟到数字的离散化实现

3.1 双线性变换法

在风电PLC中，我们需要将模拟滤波器转换为数字形式。采用双线性变换：

code复制s = (2/T)(z-1)/(z+1)

其中T为采样周期。代入后得到数字滤波器传递函数：

code复制H(z) = (b₀ + b₁ z⁻¹ + b₂ z⁻²) / (1 + a₁ z⁻¹ + a₂ z⁻²)

系数计算公式：

code复制b₀ = b₂ = (ωc T)² / D
b₁ = 2(ωc T)² / D 
a₁ = [2(ωc T)² - 8] / D
a₂ = [4 - 2√2 ωc T + (ωc T)²] / D
D = 4 + 2√2 ωc T + (ωc T)²

以一个实际案例说明：当ωc=100π rad/s（50Hz截止），T=0.001s（1kHz采样）时，计算得：

code复制b₀=b₂=0.0675, b₁=0.1350
a₁=-1.4193, a₂=0.5485

3.2 频率预畸变补偿

双线性变换会导致频率轴非线性扭曲。对于风电这类需要精确截止频率的场景，必须进行预畸变校正：

code复制ωc = (2/T) tan(ωd T/2)

其中ωd是期望的数字截止频率。例如当需要数字域50Hz截止频率时，实际应设计的模拟截止频率为：

code复制ωc = 2000π tan(π/20) ≈ 100.53π rad/s

这个约0.5%的偏差在MW级风机控制中足以引起明显的性能差异。

4. 风电场景下的参数整定要点

4.1 截止频率选择原则

根据不同类型信号的特征建议：

• 变桨系统转速反馈：15-30Hz（兼顾响应速度与噪声抑制）
• 功率调节回路：5-10Hz（平滑电网波动）
• 振动监测信号：100-200Hz（需保留故障特征）

一个实用的经验公式：

code复制ωc = 2π × (0.2~0.5) × 叶片通过频率

例如对于转速18rpm的三叶片风机，叶片通过频率为0.9Hz，则ωc可选2π×0.4Hz≈2.5rad/s。

4.2 量化误差预防

在16位定点DSP实现时，需特别注意：

1. 系数归一化：确保所有系数绝对值小于1
2. 采用直接Ⅱ型结构：减少延迟单元数量
3. 增加溢出保护：在每级输出设置饱和限幅

某项目中的教训：未做归一化的系数a₁=-1.4193直接存入Q15格式会导致溢出，正确做法是整体除以2，然后补偿增益：

code复制b₀'=b₂'=0.0337, b₁'=0.0675
a₁'=-0.7096, a₂'=0.2742

最终输出需乘以2补偿增益。

5. 实际工程问题排查指南

5.1 相位滞后补偿方案

当滤波器导致系统相位裕度不足时，可采用：

1. 前馈补偿：在控制指令通道并联相同滤波器
2. 相位超前校正：在滤波器后追加(τs+1)/(ατs+1)环节
3. 降低滤波器阶数：改用两个一阶滤波器串联

5.2 常见异常现象分析

在某海上风电项目调试中，发现变桨响应出现5Hz振荡，最终定位是滤波器ωc设置过高（50Hz）放大了测量噪声。调整至25Hz后，系统抖动消失且响应速度仍满足要求。