1. 项目背景与核心目标
在电力电子系统中，谐波污染一直是影响电能质量的关键问题。非线性负载（如整流器、变频器等）的广泛应用导致电网电流波形严重畸变，传统无源滤波器难以满足动态补偿需求。有源电力滤波器(APF)因其灵活性和高效性成为谐波治理的主流方案，而控制策略的设计直接决定了APF的性能上限。
本项目针对这一技术痛点，提出了一种基于PI控制与重复控制复合策略的谐波抑制方案。通过Simulink建模仿真验证，该方案实现了以下突破性指标：
- 总谐波畸变率(THD)从初始的5%以上降至1%以内
- 动态响应时间控制在10ms以内
- 兼容Matlab 2015/2017双版本环境
关键创新点：将PI控制的快速响应特性与重复控制的周期精确跟踪能力相结合，通过并联复合控制架构实现动态与稳态性能的平衡。
2. 控制策略原理深度解析
2.1 重复控制的内模原理
重复控制的核心思想源于内模原理(Internal Model Principle)——要使控制系统实现对特定信号的零稳态误差跟踪，控制器必须包含该信号的动力学模型。对于周期性谐波抑制场景，其数学模型可表述为：
连续域传递函数： $$ G_{rc}(s) = \frac{e^{-Ts}}{1-e^{-Ts}} $$
离散化后（Z变换）： $$ G_{rc}(z) = \frac{z^{-N}}{1-z^{-N}} $$ 其中N=fs/f1表示一个基波周期对应的采样点数。
这个结构本质上是一个时滞正反馈系统，其工作原理可以类比"学习记忆"过程：
- 误差累积：每个周期将新的误差信号叠加到历史误差上
- 周期修正：利用记忆的误差模式生成补偿信号
- 相位补偿：通过超前环节zk抵消数字控制延迟
2.2 PI控制的关键作用
虽然重复控制在稳态精度上表现优异，但其固有缺陷需要通过PI控制来弥补：
| 性能指标 | 重复控制缺陷 | PI控制补偿方案 |
|---|---|---|
| 动态响应 | 滞后1个基波周期(20ms@50Hz) | 比例环节提供瞬时响应 |
| 非周期扰动 | 仅对周期性信号有效 | 积分环节抑制随机扰动 |
| 参数敏感性 | Q(z)需精确整定 | 提供基础稳定裕度 |
实测数据表明，单独使用重复控制时，负载突变下的恢复时间长达30ms，而PI+重复复合控制可将此时间缩短至8ms。
3. Simulink实现细节
3.1 系统整体架构
模型采用模块化设计，主要包含以下子系统：
APF_System.slx ├── Power_Circuit/ # 主电路 │ ├── LCL_Filter/ # 三阶LCL滤波器 │ └── IGBT_Inverter/ # 三相全桥逆变器 ├── Control_Algorithm/ # 控制算法 │ ├── PLL/ # 锁相环 │ ├── PI_Controller/ # 电流环PI控制 │ └── Repetitive_Control/ # 重复控制器 └── Measurement/ # 测量分析 ├── THD_Analysis/ # 谐波分析 └── Scope_Manager/ # 示波器管理3.2 关键模块参数配置
LCL滤波器设计
% 设计参数 Prated = 10e3; % 额定功率10kW Vgrid = 220; % 电网电压(相电压) fgrid = 50; % 电网频率 fsw = 10e3; % 开关频率 % 计算公式 L1 = 0.1*Vgrid^2/(2*pi*fgrid*Prated); % 网侧电感3mH Cf = 0.05*Prated/(2*pi*fgrid*Vgrid^2); % 滤波电容15uF L2 = L1/5; % 逆变侧电感0.6mH重复控制器实现
function [u_rc] = Repetitive_Control(e, N, Kr, Q, k) % e: 误差输入 % N: 延迟点数 = fs/f1 % Kr: 增益系数(0.95) % Q: 低通滤波器系数(0.95) % k: 相位超前拍数(4) persistent buffer; if isempty(buffer) buffer = zeros(N+1,1); end % 更新延迟线 buffer = [e; buffer(1:end-1)]; % 相位补偿 if k > 0 u_rc = Kr * (buffer(k) + Q*buffer(end)); else u_rc = Kr * buffer(1); end end3.3 调试技巧与注意事项
稳定性保障：
- 先单独调试PI环，确保系统基本稳定
- 引入重复控制时，初始设置Kr=0.5逐步增加
- 用Bode图验证相位裕度>45°
参数整定流程：
graph TD A[确定基波频率f1] --> B[设置N=fs/f1] B --> C[调整PI参数Kp,Ki] C --> D[设置Q=0.9~0.95] D --> E[逐步增加Kr至0.95] E --> F[优化超前拍数k]常见问题处理：
- 振荡现象：降低Kr或Q值，检查PLL同步状态
- 补偿不足：确认N值是否正确，检查ADC采样同步
- 高频噪声：在Q(z)中加入二阶低通滤波器
4. 仿真结果分析
4.1 谐波抑制效果对比
负载条件：三相不控整流带阻感负载(R=10Ω, L=10mH)
| 谐波次数 | 补偿前含量(%) | 单独PI控制(%) | PI+重复控制(%) |
|---|---|---|---|
| 5th | 23.5 | 4.2 | 0.8 |
| 7th | 12.1 | 2.7 | 0.5 |
| 11th | 8.3 | 1.9 | 0.3 |
| THD | 28.7 | 5.4 | 0.9 |
注：测试数据基于Matlab 2017b，采样频率10kHz
4.2 动态响应波形解读
图1展示了负载突加时的电流跟踪过程：
- t=0.1s时负载投入，电流THD瞬时升至25%
- PI控制器在3ms内实现初步抑制(THD→8%)
- 重复控制器在后续两个周期(40ms)内将THD降至1%以下
5. 工程实践建议
5.1 硬件实现考量
处理器选型：
- 最小要求：150MHz主频DSP(TMS320F28335级别)
- 推荐方案：Xilinx Zynq-7020 SoC(FPGA+ARM架构)
关键时序约束：
// 中断服务程序时间分配 void ISR_10kHz(void) { ADC_Read(); // 50us PLL_Update(); // 20us Control_Calc(); // 80us (PI+重复控制) PWM_Update(); // 10us // 总耗时<160us (<80% CPU负载) }
5.2 进阶优化方向
参数自适应：
% BP神经网络参数在线调整示例 net = feedforwardnet([10 5]); net = train(net, [THD; dTHD]', [Kp; Ki]');多采样率控制：
- 电流环：10kHz高速采样
- 重复控制：1kHz更新周期
- 电压环：500Hz带宽
新型拓扑扩展：
- T型三电平APF
- 模块化多电平APF(MMC-APF)
6. 模型获取与使用说明
本模型提供两个版本适配不同Matlab环境：
- 2015版：使用传统Simulink库元件
- 2017版：采用新型电力系统元件库
使用前请确认：
- 安装SimPowerSystems工具箱
- 运行前执行初始化脚本Init_APF.m
- 建议步长设置为1e-6s
典型问题排查：
- 仿真不收敛：尝试减小步长或改用ode23tb求解器
- 波形异常：检查接地连接和信号单位(标幺值/实际值)
- THD分析失败：确认FFT窗口包含整数个周期
通过实际工程验证，该方案在工业现场应用中可将THD稳定控制在2%以内（符合GB/T 14549-93标准），特别适用于轧钢机、电弧炉等严苛工况。
