在电能质量治理这个行当干了这些年，三相整流设备带来的谐波问题一直是块硬骨头。前阵子把基于PI+重复控制的有源滤波器（APF）谐波抑制策略在Simulink里完整搭了一遍，从检测环节到复合控制器再到主电路联动，整个过程踩了不少坑，但也把每个环节的逻辑理清了。这篇就把整个建模思路、参数计算和调试过程掰开揉碎讲清楚，给打算在Simulink里复现APF控制策略的同学一条能直接走通的路。
1. 方案选型与整体设计思路
1.1 有源滤波器为什么能治谐波
先聊清楚APF的本质。电网里的谐波源大多是非线性负载，比如三相不控整流桥带阻感负载，它从电网吸取的电流不是正弦波，而是含有大量5、7、11、13次等特征谐波的畸变电流。传统的无源滤波器（PPF）是靠LC谐振回路把特定次谐波短路掉，缺点很突出：只能滤固定次数谐波，电网阻抗变化时容易失谐，甚至可能与系统发生谐振放大谐波。
有源滤波器的思路完全不同——它不是"堵"，而是"顶"。APF实时检测负载电流中的谐波分量，生成一个幅值相等、相位相反的补偿电流注入电网，让电网侧电流只剩下基波。这就好比你在一个嘈杂的房间里，其他人同时说话，你放出一个"反向噪声"把特定频段抵消掉。关键在于"实时"和"反向"，所以谐波检测速度和控制精度直接决定补偿效果。
用Simulink做APF仿真，最大的价值在于可以把整个链条拆开看：电网、非线性负载、APF主电路、检测环节、控制环节、PWM调制，每一个模块的信号都能拉出来观察。实际装置里你很难单独看谐波检测的误差，但仿真里随便一根线都能挂示波器，这对理解控制策略的内部机理帮助非常大。
1.2 PI和重复控制，为什么要凑在一起用
单用PI控制做电流内环，跟踪正弦指令是存在稳态误差的。原因不复杂：PI控制器对直流信号可以做到无静差，但对某个频率的正弦信号，开环增益是有限的，而且频率越高增益越低。谐波电流包含多个频率分量，5次是250Hz，7次是350Hz，想要把这么多频率的正弦信号全部精准跟踪，单靠一个PI环根本忙不过来。实测下来，纯PI控制的APF补偿后电网电流THD通常只能压到8%到12%，远达不到国标5%的要求。
重复控制的引入就是冲着"周期性波形跟踪"来的。它基于内模原理——如果控制器内部包含一个与外部信号同周期的信号发生器模型，那么对这个周期信号的跟踪误差可以趋近于零。换句话说，重复控制器相当于把谐波电流的波形"记住"了，第一个周期动态响应差一点，但从第二个周期开始，它能精确预测并补偿出同样的畸变成分。
那为什么不全用重复控制？因为它先天动态响应慢——它的补偿动作要延迟一个基波周期（20ms）才生效。负载突变时，这20ms里谐波几乎完全靠PI硬扛，效果很差。所以工程上最常见的做法是并联复合：PI负责动态响应，快速跟上负载变化；重复控制负责稳态精度，把周期性谐波压到极致。这个"PI打天下、重复控精度的"思路，正是这个项目核心的选型逻辑。
1.3 整体仿真架构怎么搭
整个Simulink模型我分成了四大部分：主电路部分（三相电源、非线性负载、APF逆变器、滤波电感）、谐波检测部分（基于瞬时无功功率理论的dq检测法）、复合控制器部分（PI环节、重复控制环节、并联合成）、PWM调制部分（三角载波比较或SVPWM）。
架构上最需要注意的是控制信号与主电路信号的隔离。Simulink里电力电子仿真容易出数值发散问题，如果不加隔离直接连，经常出现仿真一步就崩的情况。我实际的做法是：主电路用Simscape Electrical（以前叫SimPowerSystems）模块搭建，控制部分全部用普通Simulink离散模块，两者之间通过电压/电流测量模块转接，信号地注意区分。这样既保留了主电路的物理模型精度，控制部分又能按数字控制器的逻辑走离散采样周期。
2. 谐波检测与指令电流生成
2.1 dq坐标系下的检测原理
谐波检测是APF的"眼睛"。电流采歪了，后面控制再精确也没用。目前工程上最主流的是基于瞬时无功功率理论的检测法，核心逻辑就是坐标变换滤波。
三相静止abc坐标系下的负载电流，经过Clarke变换变成两相静止αβ坐标系下的分量，再经过Park变换转到与电网电压同步旋转的dq坐标系下。关键在于：在同步旋转坐标系里，基波正序分量变成了直流分量，而5、7、11、13次等谐波分量变成了不同频率的交流分量。这时候只要用低通滤波器把直流分量提取出来，剩下被滤掉的就是谐波分量。再用反变换把谐波分量换回abc坐标系，就是APF需要的补偿指令电流。
用Simulink实现时，Park变换的角度θ必须来自锁相环（PLL），这个角度一旦抖动，检测结果直接变差。我用的是三相锁相环模块，仿真里PLL调参要注意：环路带宽设得太宽，锁相结果里会混入谐波噪声，角度曲线毛刺多；太窄则动态跟踪慢，电网频率偏移时跟不上。实测下来，PLL的自然频率设在20Hz到30Hz之间，阻尼比取0.707，效果比较均衡。
2.2 低通滤波器选型与延时取舍
检测环节里LPF的参数直接影响指令电流的精度和延迟。LPF的截止频率选高了，会有一部分低次谐波漏过去当成"基波"滤掉，补偿指令里就缺少这部分分量；选低了，动态响应变慢，负载突变时基波提取跟不上，补偿指令会产生偏差。
我试过二阶Butterworth低通滤波器，截止频率放在150Hz到200Hz之间。为什么不放更低？因为基波是50Hz，谐波最低是250Hz（5次），中间这段频率就是天然的隔离带，150Hz到200Hz在两者之间留了足够余量。但注意，滤波器阶数越高或截止频率越低，相位延迟越大，这会让检测出的谐波指令电流滞后于实际谐波，补偿效果打折。实际装置里通常还要加相位补偿环节，仿真阶段可以先忽略，但要把这个影响记在心里。
还有一点容易被忽略：LPF会让负载突变后的指令电流建立过程变慢。我在模型里实测过，二阶Butterworth、截止频率200Hz时，负载阶跃后大约经过3到5ms指令电流才能趋于稳定。这个延迟靠PI的快速响应去弥补，系统整体动态指标基本不受影响。
2.3 直流侧电压控制与指令修正
APF要输出补偿电流，逆变器直流侧必须维持稳定电压。直流侧电容电压跌了，补偿电流的幅值就上不去，波形还会畸变。直流电压控制通常单独用一个PI环，采样的电容电压与给定值（我用的750V）做差，经PI调节后输出一个有功指令分量，叠加到谐波指令电流的d轴上。
这里有个关键操作：谐波检测环节里LPF滤掉的是基波分量，反变换出来的纯粹是谐波指令，它的d轴分量平均值是零。如果把直流电压环的输出直接叠加到指令电流的d轴上，就能在不影响谐波补偿的前提下，让APF从电网吸收少量有功来维持直流母线电压。这个修正量必须加在坐标变换之后的d/q轴上，不能在abc坐标系里直接加，我在模型里踩过这个坑，一开始加错了位置，直流电压完全稳不住。
3. 复合控制器设计与参数整定
3.1 电流环PI参数计算
电流内环PI参数的计算要依据主电路的滤波电感值。APF的等效电路里，逆变器输出电压加在滤波电感上，产生补偿电流，所以电流环的被控对象近似为一阶惯性环节，其增益和电感值强相关。
我用的主电路参数是：电网线电压380V、直流侧电压750V、滤波电感3mH、开关频率10kHz。电流环PI参数按照零极点对消法来整定，先定电感L=3mH，考虑开关频率对应的控制周期Ts=0.0001s。Kp的初始估值按L/(Ts×2)算出来在15左右，再根据仿真波形微调。Ki的取值要保证稳态无静差，但积分作用太强会造成超调甚至振荡，我最后定在Kp=12、Ki=1000附近，这个组合在阶跃响应测试中表现得比较稳。
这里提醒一下：把Kp调大确实能加快动态响应，但过大的Kp会把开关纹波放大，电感电流上叠加明显的高频毛刺，这对PWM调制非常不友好，电流环输出噪声太大会导致调制波抖动。整定原则是先保稳定性，再追响应速度。
3.2 重复控制器参数配置
重复控制器的离散化传递函数结构是：
Grc(z) = Kr × z^(-N) × S(z) / (1 - Q(z) × z^(-N))
N是一个基波周期内的采样点数。我开关频率和采样频率都是10kHz，基波50Hz，那么N = 10000/50 = 200。这个N不能设错，设成180或者220，重复控制器的谐振峰就会移到错误频率上，补偿效果直接崩掉。
Q(z)一般取常数0.95到0.98，它决定内模的"记忆衰减程度"。取1时内模是一个理想周期信号发生器，控制精度最高但鲁棒性差，一旦系统参数有偏差容易不稳定。取0.95时鲁棒性好一些，但稳态误差会略微增大。我在仿真里取的是0.97，稳定性和精度都比较满意。
S(z)是相位幅值补偿器，通常用一个超前环节或滤波器实现。它补偿被控对象在中低频段的相位滞后，让重复控制的补偿指令更准。这里我用的S(z)直接把电流环闭环模型逆过来近似，再串一个低通滤波器滤掉高频增益。这个补偿器的设计直接决定重复控制的收敛速度和稳态误差，是整个重复控制器调试中最需要耐心的一环。
Kr是重复控制增益，取1时理论上收敛最快，但容易引发振荡。实际测试中Kr取0.8到0.9区间，系统稳定性明显更好，收敛速度也没有牺牲太多。
3.3 PI与重复控制并联结构的融合
并联复合控制的核心问题是两个控制器输出如何合成。我用的结构是：PI控制器和重复控制器各算各的，输出直接相加得到总控制量，送入PWM调制。
但直接相加有一个隐患：PI环带宽内的低频段和重复控制的基波倍频处增益重叠，两个控制器"抢活干"，可能造成合成控制量波动增大。解决办法是在PI通道上串一个高通滤波器，让PI主要管高频动态；或者在重复控制通道上加一个基波周期的延迟匹配，让两者合成的相位一致。我实际用的是后一种思路，重复控制器天然有一个周期的延迟，PI那边加上同样的延迟补偿，确保两个控制器的输出在时间上对齐，合成波形明显更光滑。
还有一点经验：启动阶段重复控制器还没建立记忆，输出是零，整个系统相当于纯PI控制。要等一个周期之后重复控制才逐渐发挥作用并在几个周期后收敛。这个收敛过程可以在仿真里看得非常清楚——补偿后的THD是逐周期下降的，前两三个周期甚至能看到明显的波动，之后才稳定下来。
4. Simulink仿真建模实操
4.1 主电路建模
主电路用Simscape Electrical模块搭建：三相交流电压源做电网，线电压380V、频率50Hz，内阻设一个小值模拟线路阻抗。非线性负载用三相不控整流桥带RL负载来实现，整流桥交流侧电流就是典型的尖顶波谐波源，能够提供充足的5、7次谐波。
APF主电路部分包括：IGBT逆变桥、直流侧电容、滤波电感（L=3mH）。直流侧电容的取值也是个门道，太小则直流电压纹波大，太大则建压过程太慢。我用的电容是4700uF，母线电压750V，建压到稳定的时间大约在0.2s到0.3s。滤波电感的选择直接影响纹波电流大小，我直接取3mH，与电流环计算参数保持统一。注意电感值改小会显著降低谐波补偿能力，这是主电路参数里最敏感的一个量。
4.2 控制模块建模
控制部分全部用离散模块实现，采样周期Ts=1e-4s。整个链路是：电流采样→Clarke变换→Park变换→LPF提取基波→反变换得谐波指令→叠加直流电压环的有功修正量→电流误差进入复合控制器→输出控制量→PWM调制。
锁相环输出角度θ供Park变换使用，角度信号需要做unwrap处理，否则跨过2π边界时角度跳变会让坐标变换结果突变，补偿电流会出现莫名其妙的大尖峰。Simulink里直接用PLL模块的输出信号就可以，但要注意数据类型的匹配。
PWM调制我用的是三角载波比较，载波频率10kHz。把复合控制器的输出与三角波比较生成IGBT的门极信号。APF主电路里调制需要带死区吗？仿真里如果用的是理想开关模型，死区可以设得非常小或者不设，但如果用热模型IGBT，就一定要加死区时间，否则桥臂直通。我的做法是先用理想开关调通控制算法，最后再换成带死区的模型验证工程可行性——先保证控制部分正确，再逐步逼近实际装置。
4.3 仿真参数设置与调参顺序
仿真参数这块我得说重点。Simulink做电力电子仿真，步长设置不对是大部分发散问题的根源。我这里用的是离散定步长仿真，步长Ts_Sim=1e-6s，这比控制周期1e-4s小两个数量级，保证PWM高频开关过程能够被准确捕获。如果你想加快仿真速度，可以把PWM调制也离散化到采样周期里，变成"采样-计算-输出"的零阶保持模式，这样仿真步长可以放大到2e-5s，速度提升明显，但波形细节会损失一些。
调参顺序上，我的经验是：先只开PI环，把电流内环调稳，观察负载阶跃时补偿电流能否快速跟随指令；再接入谐波检测，看指令电流是否干净；最后再叠加重复控制，逐周期看THD下降趋势。每步都确认之前没动过的模块输出正常，再动下一步，否则出了问题根本不知道往哪里查。
5. 典型问题定位与调试心得
5.1 仿真发散或持续振荡怎么办
我在调模型时遇到过两次典型的振荡问题。第一次是刚把复合控制器闭环接好，仿真直接发散到NaN，查了半天发现是Park变换的角度没有unwrap，角度在2π处跳变导致坐标变换结果突变，控制量跟着爆炸。第二次是电流环Kp调太大了（试到30），高频毛刺被放大，PWM调制输出乱跳，电感电流出现持续振荡。把Kp降回12附近，加上滤波器之后波形立刻恢复正常。
最让我意外的问题出在重复控制的N值上。有一版模型我把采样频率看错了，以为采样是20kHz，N设成了400，结果重复控制补偿的频率偏了一倍，补偿后THD不减反增。这件事给我的教训很深刻：重复控制器频率匹配不准，效果可能比不加还差，而且很难通过调其他参数救回来。
5.2 谐波补偿效果不达标的排查
如果接上APF后电网电流THD仍然居高不下，一般按这个优先级排查：先看检测环节的LPF截止频率是否合适，指令电流是否包含了主要的5、7次谐波分量；再看电流内环的跟踪带宽是否足够，补偿电流与指令之间是否存在明显幅值相位误差；最后看重复控制的收敛情况，波形应该是逐周期改善的。
还有个容易被忽略的环节是直流侧电压。如果直流母线电压跌到600V以下，逆变器输出能力受限，高幅值的谐波补偿指令根本执行不出来，THD自然下不来。我测试过直流电压从750V降到650V，补偿效果肉眼可见地变差。这种情况不是控制策略的问题，而是主电路容量约束的问题。
5.3 参数整定的几条实操经验
关于参数整定，我把最终确定的一套参数列在下面，这套参数在模型里验证过补偿效果是稳定的，负载从空载到满载突变也能在几个周期内恢复：
| 参数项 | 取值 | 说明 |
|---|---|---|
| 电网电压 | 380V / 50Hz | 三相交流源 |
| 直流侧电压给定 | 750V | PI闭环控制 |
| 直流侧电容 | 4700uF | 电压纹波与建压速度折中 |
| 滤波电感 | 3mH | 电流环对象参数 |
| 采样/开关频率 | 10kHz | Ts=1e-4s |
| 电流环PI | Kp=12，Ki=1000 | 初始零极点对消，波形微调 |
| 直流电压环PI | Kp=0.5，Ki=20 | 慢环，避免与电流环耦合 |
| 重复控制N | 200 | 10kHz/50Hz |
| Q(z) | 0.97 | 稳定与精度折中 |
| Kr | 0.85 | 收敛速度与稳定折中 |
| LPF截止频率 | 200Hz | 二阶Butterworth |
几条我实际摸索出来的经验：重复控制的Kr不要一上来就取1，先用0.7验证稳定性再逐步加；LPF截止频率宁可高一点，动态响应会好很多，稳态精度靠重复控制去兜底；直流电压环的带宽必须远低于电流环，两个环路之间一旦发生耦合振荡，现象非常难排查。
做完整套模型下来，我最深的体会是：复合控制策略的难度不在理论推导，而在数字实现细节。坐标变换的角度处理、离散化步长的匹配、补偿器的相位对齐，每一项偏差都会在最终THD指标上体现出来。Simulink的优势在于你可以把每个中间信号单独捞出来看，我在调重复控制收敛过程的时候，就是盯着指令电流误差信号，看着它一个周期一个周期地衰减，那种感觉比看最终THD数据直观得多。
最后再分享一个小技巧：调试重复控制时，把补偿后的电网电流波形和负载电流波形放到同一个示波器里对比，你会发现前几个周期补偿效果很差，但波形正在逐渐"变干净"。如果看不到这种逐周期改善的趋势，说明重复控制参数有问题，趁早回头检查N值匹配，别在那里瞎调增益浪费时间。这组模型调通之后，我把它改成LCL滤波器结构做进一步验证，效果更好，但参数整定的复杂度又上了一个台阶——那就是另一个故事了。
