1. 项目背景与核心目标
在工业电力系统中，谐波污染一直是影响电能质量的关键问题。随着电力电子设备的广泛应用，非线性负载产生的谐波电流会导致电网电压畸变、设备过热甚至保护误动作。传统无源滤波器虽然结构简单，但存在谐振风险且无法动态适应负载变化。有源电力滤波器(APF)因其动态补偿能力成为解决这一问题的有效方案。
本项目针对APF的核心挑战——如何实现高精度谐波抑制，提出了一种基于PI控制与重复控制相结合的复合控制策略。通过Simulink仿真验证，该方案能将总谐波畸变率(THD)控制在1%以下，显著优于单一控制方案。特别值得关注的是，模型提供了对Matlab 2015和2017版本的双重支持，确保了工程实践的兼容性。
2. 控制策略原理深度解析
2.1 重复控制的内模原理本质
重复控制的核心思想源于内模原理(Internal Model Principle)——要使控制系统实现对特定信号的零稳态误差跟踪，控制器必须包含该信号的动力学模型。对于周期性谐波抑制这一典型场景：
-
基础内模结构：传统方法需要对每个谐波频率设计独立的内模，例如对5次谐波需构建250Hz(50Hz×5)的正弦模型。这种方案在工程实现上存在明显缺陷：随着谐波次数增加，控制器复杂度将呈指数级增长。
-
周期信号统一建模：观察到电力谐波具有严格的周期性特征（每个基波周期波形重复），我们采用改进的内模结构：
math复制G_{IM}(s) = \frac{e^{-Ts}}{1 - e^{-Ts}}其中T为基波周期(20ms@50Hz)。该模型通过时延环节实现了对所有谐波成分的统一描述，避免了多模型并联的复杂性。
2.2 PI与重复控制的协同机制
动态响应与稳态精度的平衡
-
PI控制优势：比例环节提供快速误差响应，积分环节消除稳态偏差。在APF中，PI控制器能在10ms内响应负载突变，但其对周期性谐波的抑制存在固有局限——积分作用仅能保证直流误差为零，对交流谐波的抑制效果有限。
-
重复控制特性：通过周期延迟正反馈构建"学习"机制，能渐进消除所有周期性误差。但其固有延迟（至少一个基波周期）导致暂态响应迟缓。实测表明，单独使用重复控制时动态响应时间超过20ms。
并联复合控制结构
采用如图1所示的并联架构，实现优势互补：
code复制 +-------+
e(t) --->| PI |---> u_PI(t)
+-------+ +
===> u_total(t) --> APF
+-------+ +
e(t) --->| 重复 |---> u_RC(t)
+-------+
- 信号合成：
u_total(t) = u_PI(t) + u_RC(t) - 分工协作：
- PI控制器快速补偿突变量（如负载阶跃变化）
- 重复控制器逐步修正周期性残余误差
- 通过增益调节（如Kr=0.95）避免两者输出冲突
2.3 离散化实现关键步骤
为适应数字控制需求，需将连续域模型转换为离散形式：
-
时延环节离散化：
math复制e^{-Ts} \xrightarrow{Z变换} z^{-N}, \quad N=T/T_s=f_s/f_1例如当采样频率fs=10kHz、基频f1=50Hz时，N=200。
-
稳定性增强设计：
- 添加低通滤波器Q(z)（常取0.95）抑制高频噪声
- 相位补偿器S(z)=z^k补偿计算延迟（典型k=4）
-
最终离散传递函数：
math复制G_{RC}(z) = \frac{K_r z^k Q(z)}{1 - z^{-N} Q(z)}
3. Simulink建模与参数整定
3.1 主电路建模要点
LCL滤波器设计
APF输出侧采用LCL滤波器（图2），其参数选择遵循：
-
谐振频率应避开主要谐波频段：
math复制f_{res} = \frac{1}{2π}\sqrt{\frac{L_1+L_2}{L_1 L_2 C}}通常设计在开关频率(如10kHz)与最高次谐波(如25次=1.25kHz)之间。
-
阻尼设计：为避免谐振峰，可采用：
- 无源阻尼：串联电阻（增加损耗）
- 有源阻尼：电容电流反馈（推荐方案）
坐标变换实现
通过abc-dq0变换解耦三相电流：
matlab复制% Clarke变换
Iα = (2/3)*(Ia - 0.5*Ib - 0.5*Ic);
Iβ = (1/sqrt(3))*(Ib - Ic);
% Park变换
Id = Iα*cosθ + Iβ*sinθ;
Iq = -Iα*sinθ + Iβ*cosθ;
在dq坐标系中，谐波表现为交流量，便于控制器设计。
3.2 控制器参数整定指南
PI参数设计
采用带宽法确定Kp、Ki：
- 确定目标带宽fc（通常取1kHz）
- 计算比例增益：
math复制其中Leq为等效电感K_p ≈ 2π f_c L_{eq} - 积分时间常数：
math复制Tsw为开关周期T_i = \frac{K_p}{K_i} ≈ \frac{1}{10} T_{sw}
重复控制参数
| 参数 | 设计准则 | 典型值 |
|---|---|---|
| Kr | 接近1但需留稳定裕度 | 0.95 |
| Q(z) | 决定高频衰减特性 | 0.95 |
| k | 补偿计算延迟（1-2个采样周期） | 4 |
| N | 严格匹配基波周期 | fs/f1 |
调试技巧：先单独调试PI至稳定，再逐步增加Kr观察THD改善
4. 仿真结果与性能分析
4.1 稳态性能对比
在整流负载条件下测试（表1）：
| 控制策略 | THD(%) | 主要谐波残余 |
|---|---|---|
| 传统PI | 4.03 | 5次(3.2%)，7次(2.1%) |
| 单独重复控制 | 2.11 | 各次均<1%但存在相位滞后 |
| PI+重复复合 | 1.97 | 所有谐波<0.5% |
关键波形对比如图3所示：
- 复合控制下网侧电流正弦度显著改善
- 负载电流突变后10ms内恢复稳定
4.2 动态响应测试
设计阶跃负载变化场景：
- t=0.1s时负载从50%突增至100%
- 观测指标：
- 调节时间：<15ms
- 超调量：<20%
复合控制展现出优于单独重复控制的动态特性，同时保持PI控制的快速性。
5. 工程实践中的挑战与解决方案
5.1 非理想工况应对
电网频率波动
当电网频率偏移±2Hz时：
- 问题：固定N值导致内模失配
- 解决方案：
matlab复制N = round(fs/f_actual); % 实时更新延迟节拍
非周期干扰
对随机脉冲干扰：
- 增加滑动平均滤波预处理
- 设置重复控制死区阈值
5.2 硬件实现考量
计算资源优化
重复控制的200拍延迟（fs=10kHz时）需大量存储：
- 采用环形缓冲区实现
- 使用DSP的延迟指令（如TI C2000系列的DELAY_US）
抗混叠设计
- ADC采样前需至少2阶抗混叠滤波器
- 推荐截止频率：
math复制f_c = \frac{f_s}{2.5}
6. 模型使用指南
6.1 版本适配说明
模型包包含两个版本：
- 2015版：使用R2015a语法
- 2017版：支持新型电力系统模块集
6.2 关键模块配置
-
重复控制器：
- 设置正确的N值（在InitFcn中计算）
- 调试阶段可先设Q=1观察稳定性
-
THD分析模块：
- 建议测量窗口取10个基波周期
- 谐波次数设置到50次
6.3 扩展接口
模型预留了以下扩展点：
- 神经网络参数整定接口
- 第三方硬件在环(HIL)测试接口
实际部署时建议先进行控制器在环(CIL)测试验证实时性