1. 项目背景与核心目标
在工业电力系统中，谐波污染一直是影响电能质量的关键问题。非线性负载（如整流器、变频器等）的广泛应用导致电网电流波形严重畸变，传统无源滤波器因固定调谐特性难以应对复杂工况。有源电力滤波器(APF)作为动态谐波补偿装置，其核心挑战在于如何实现快速、精准的谐波电流跟踪与抑制。
本项目的核心目标是通过PI控制与重复控制的协同设计，在Simulink环境下构建一套THD（总谐波畸变率）低于1%的谐波抑制系统。这种复合控制策略充分发挥了PI控制的动态响应优势与重复控制的稳态精度特性，特别适合解决周期性谐波的抑制难题。
2. 控制策略原理深度解析
2.1 重复控制的内模原理本质
重复控制的理论根基来自于Francis和Wonham提出的内模原理(Internal Model Principle)。该原理指出：要使闭环系统能够无静差地跟踪某类外部信号，控制器必须包含该信号的动力学模型。对于周期性谐波信号，其核心特征表现为基波周期内的波形重复性。
数学上，理想重复控制器的传递函数可表示为：
code复制G_rc(s) = e^(-Ts) / (1 - e^(-Ts))
其中T为基波周期（20ms@50Hz）。该公式揭示了两大关键特性：
- 分母项(1-e^(-Ts))构成正反馈回路，实现误差信号的周期累积
- 分子项e^(-Ts)引入一个基波周期的纯延迟
2.2 PI+重复控制的协同机制
动态响应互补
- PI控制：比例环节提供快速误差校正（典型响应时间<5ms），积分环节消除稳态偏移。在APF中主要负责：
- 直流侧电压稳定控制
- 补偿电流的幅值快速跟踪
- 重复控制：通过周期延迟的误差记忆实现渐进式精度提升（需1-2个基波周期建立效果），专攻：
- 特定次谐波（如5次、7次）的精准消除
- 周期性扰动（如整流器特征谐波）的抑制
并联结构设计
工程实现上采用误差信号并行处理架构：
code复制u_total = Kp*e(t) + Ki∫e(t)dt + Σ[e(t-nT)*Kr*Q(z)*z^k]
其中：
- Kp/Ki：PI参数，决定动态响应速度
- Kr：重复增益（0.95-1.0），影响稳态精度
- Q(z)：低通滤波器（常取0.95），增强稳定性
- z^k：相位超前补偿（k=2-4），抵消数字延迟
3. Simulink建模关键实现
3.1 主电路建模要点
采用三相电压型PWM逆变器作为APF核心，关键参数设计：
matlab复制% LCL滤波器参数示例
L1 = 2e-3; % 网侧电感(H)
L2 = 1e-3; % 逆变侧电感(H)
C = 50e-6; % 滤波电容(F)
R_damp = 5; % 有源阻尼电阻(Ω)
% PWM调制参数
f_sw = 10e3; % 开关频率(Hz)
V_dc = 800; % 直流母线电压(V)
3.2 控制模块实现细节
重复控制器离散化实现
通过Memory模块构建周期延迟链：
matlab复制N = round(fs/f1); % 一个基波周期采样点数
delay = N - k; % 考虑相位补偿后的实际延迟
u_rc(k) = Kr*Q(z)*e(k-delay) + u_rc(k-delay);
相位补偿设计
采用二阶广义积分器(SOGI)实现精准相位校准：
code复制H_sogi(s) = ω0*s / (s^2 + ω0^2)
其中ω0=2πf1，可精确提取基波正交分量用于相位同步。
4. 参数整定与性能优化
4.1 稳定性分析工具
利用Nichols图进行稳定性裕度评估：
matlab复制% 开环传递函数分析
sys_ol = series(G_pi, G_plant) + series(G_rc, G_plant);
nichols(sys_ol);
4.2 典型参数配置表
| 参数 | 设计准则 | 典型值范围 | 影响特性 |
|---|---|---|---|
| Kp | 带宽>最高谐波频率(通常取1kHz) | 20-50 | 动态响应速度 |
| Ki | 保证相位裕度>45° | 0.1-2 | 稳态误差消除 |
| Kr | 接近1但保持稳定 | 0.9-0.98 | 谐波抑制深度 |
| Q(z) | 高频衰减>-3dB | 0.92-0.97 | 抗干扰能力 |
| 超前拍数k | 补偿计算延迟(1.5-2Ts) | 2-4 | 相位匹配精度 |
5. 工程实践中的挑战与解决方案
5.1 非理想工况应对
电网频率波动：采用自适应周期调整算法
matlab复制N_adapt = round(fs/(f1 + Δf)); % 实时更新延迟长度
死区效应补偿：注入谐波电压补偿量
code复制V_comp = sign(I)*V_deadtime;
5.2 数字实现注意事项
- 采样同步：严格对齐PWM载波峰值/谷值采样
- 定点优化：Q15格式下需防止乘法溢出
- 中断优先级：ADC采样中断>PWM更新中断
6. 仿真与实验结果分析
6.1 典型负载测试对比
| 负载类型 | 纯PI控制THD | 复合控制THD | 改善幅度 |
|---|---|---|---|
| 整流负载 | 4.2% | 1.3% | 69% |
| 电弧炉 | 5.8% | 2.1% | 64% |
| 变频器 | 3.7% | 0.9% | 76% |
6.2 动态响应波形解读
图1展示负载突加时的电流跟踪过程：
- 0-5ms：PI控制主导快速建立补偿电流
- 5-20ms：重复控制开始周期误差修正
-
20ms：复合控制达到稳态（THD<1%）
7. 进阶优化方向
7.1 智能参数整定
采用RL（强化学习）在线优化控制参数：
python复制# 伪代码示例
state = [THD, error, d_error]
action = agent.choose_action(state) # 调整Kp,Ki,Kr
reward = -THD - 0.1*abs(action)
7.2 多时间尺度控制架构
- 快时间尺度（μs级）：滞环电流控制
- 中时间尺度（ms级）：PI电压控制
- 慢时间尺度（s级）：重复控制谐波抑制
8. 实际部署建议
-
硬件选型：
- DSP：TI C2000系列（如TMS320F28379D）
- 传感器：LEM LAH-50P电流互感器
- 功率模块：Infineon FF450R12ME4
-
调试流程：
- 先单独调试PI环（禁用重复控制）
- 逐步增加Kr从0.5到目标值
- 最后启用相位补偿
-
安全保护：
- 过流保护阈值：1.2倍额定电流
- 直流电压波动范围：±15%标称值