PI并联重复控制在APF谐波抑制中的Simulink仿真实践

做有源滤波器谐波抑制仿真的朋友多半都有同一种感受单靠PI电流环做补偿电网电流的畸变率总是压不到理想水平波形上明显还留着一条条毛刺FFT一分析5次和7次谐波还是飘在那里。我在Simulink里对比了多种控制方案最后固定下来的策略是PI并联重复控制这也是目前工程上很主流的一种复合结构。这套模型适合两类人参考一是做APF课题、需要把控制策略落地成仿真的研究生二是做并网逆变器、SVG或者电能质量治理装置想快速验证控制算法的工程师。它解决的核心问题很清楚——在保证动态响应的前提下把稳态谐波补偿精度提上去。1. 方案选择的底层逻辑为什么PI必须配重复控制1.1 谐波抑制的本质是周期信号跟踪APF的本质可以理解成“反向注入电流”。非线性负载从电网抽取的电流是畸变的里面除了50Hz基波还叠着5次、7次、11次、13次这些特征次谐波。APF并联在公共连接点检测出负载电流里的谐波成分然后让逆变器输出一路和它等幅反相的电流把电网测看到的电流拉回正弦。控制目标因此非常具体补偿电流要能实时、精确地跟踪一条周期性变化的谐波指令。这条指令不是单一频率的正弦波而是多个频率分量的叠加。以最典型的三相不可控整流负载为例特征谐波是6k±1次在dq旋转坐标系下会表现为6k次纹波也就是300Hz、600Hz、900Hz这些频率。对控制器来说它面对的不是一个直流量而是一串已知周期、多频复合的周期性信号而且幅值还会随着负载变化。所以“谐波抑制”这个控制问题本质上是周期信号跟踪问题不是普通的直流稳压问题。这也是后续选控制器架构的最根本依据。1.2 PI单独扛不住谐波问题出在哪里如果只给电流环配PI控制器从原理上就有先天短板。被控对象是逆变器输出电感的电流模型传递函数近似为1/(LsR)。PI的传递函数里有一个积分环节它擅长的是把直流误差或缓慢变化误差慢慢消除但对周期性正弦指令它的带宽有限在谐波频率上提供的增益不足跟踪结果必然是幅值衰减加相位滞后。用个生活里的类比。PI控制像是一个人跟着前车尾灯开车前车直行时跟得住前车连续变道时就会慢半拍而且这个慢半拍是固定存在的。谐波指令恰恰是那种频繁变道的路况5次谐波还没跟稳7次谐波又压过来了结果就是误差里永远留着一块周期性的残差。如果为了减小这个残差一味加大Kp麻烦又来了电流环带宽一旦延伸到开关频率附近数字控制的一拍采样延时就会让相位裕度恶化高频开关纹波也被放大THD不降反升。很多人在仿真里会遇到这种“加大KpTHD先降后升”的诡异现象本质就是带宽和噪声在打架。1.3 重复控制的“逐周期修正”思路重复控制的理论基础是内模原理。简单说如果扰动是一个周期已知的信号控制器里也放进一个同样周期的信号发生器这个控制器就能对这个周期信号实现无差跟踪。离散化以后这个“周期信号发生器”可以写成G_rc(z)Q(z)·z^-N / (1 - Q(z)·z^-N)的结构其中N是一个基波周期内的采样点数Q(z)是内模增益或一个低通滤波器。这里不要被公式吓到它的实际行为非常直观。重复控制器内部维护了一段长度为N的误差记忆每到一个采样时刻它先把当前误差写进记忆再把一个周期前对应时刻的误差取出来经过增益Kr和补偿器S(z)后叠加到控制输出上。正因为有这段记忆它可以在下一个周期把上一次没补对的地方修正掉下个周期再修正一点几个周期后稳态误差就趋近于零。这就像一个学生反复刷同一套题第一次错得离谱第二次少错几个第三第四次基本就能拿满分了。它不需要反应快它靠的是“记得住”。代价也很明确重复控制在动态上天生慢一拍负载突变后至少要一个基波周期才能把新增误差逐步修完所以它不适合单独拿来扛动态。1.4 并联结构如何分工既然PI管不了稳态精度重复控制又管不了动态把两者并联起来就成了顺理成章的选择。误差信号同时送进PI控制器和重复控制器两路输出相加后作为电流环控制量。PI通道负责即时响应负载突变瞬间靠它把误差先压住重复控制通道负责精细修正逐周期把PI留下的稳态残差抹平。一个管动态一个管稳态各干各的不抢工作。并联结构还有个额外的好处是调试安全。调试时可以先单独把重复控制输出断开系统就退化成普通PI闭环所有原有分析和参数都不会失效确认PI部分没问题后再把重复控制通道慢慢加上来出了问题很容易定位是哪个通道的锅。串联结构虽然也能用但两者耦合更深参数互相影响一上来就调容易两头打架。所以我强烈建议第一次做这个策略时选并联拓扑。2. Simulink模型搭建的核心环节2.1 整体框架与模块划分我在Simulink R2022b里搭的这套模型主要用Simscape Electrical工具箱和Simulink基础库。整体拆成五块主电路、谐波检测、电流控制、PWM调制、测量与FFT分析。搭之前先花半小时画一张信号流图把“从哪测电流、指令怎么算、控制输出怎么变成占空比、PWM如何驱动开关管”这条链路理清楚后面能省下大量排查时间。主电路负责把电网、非线性负载和APF连起来谐波检测模块算出补偿电流指令电流环跑PI加重复控制调制部分用PWM Generator发驱动信号Powergui里的FFT工具负责看谐波指标。模块划分清晰之后每一个环节都可以单独验证这比把所有逻辑堆在一个大模型里要舒服得多。2.2 主电路搭建要点主电路里电网用Three-Phase Source电压设成380V/50Hz内阻按实际配电网取小值。非线性负载用三相二极管整流桥带阻感负载这是最经典的谐波源整流桥后面接一个R10Ω、L20mH的负载电流畸变率通常在25%到30%足够用来验证APF的补偿能力。APF主电路用Universal Bridge搭IGBT逆变器直流侧用电容支撑电压。这里有个细节值得多说一句直流侧电压必须高于网侧线电压峰值否则逆变器没有足够的电压裕度去注入补偿电流。我习惯先给直流电容预充到700V左右电流环调通以后再考虑换成真正的直流电压外环。补偿电感L取2mH串在逆变器输出和PCC之间这个参数直接影响电流环带宽设计后面整PI参数时要用到。搭主电路最容易踩的坑是补偿电流的方向。APF输出电流的参考方向必须和电流传感器的箭头方向、电流闭环的反馈极性保持一致任何一种约定颠倒仿真一启动就会看见电流爆表。我的习惯是先在模型草图上标好所有电流方向再对着Universal Bridge的端口说明核对一遍确认无误后再连线。2.3 谐波检测环节的实现谐波检测我用的同步旋转坐标法也就是常说的dq法。先用PLL锁相电网电压得到基波正序相位θ。负载电流经Clark变换从abc转到αβ再做Park变换从αβ转到dq。在dq坐标系下基波正序分量变成直流谐波则表现为交流纹波比如5次和7次谐波都会出现在300Hz上。用低通滤波器把直流量挑出来反变换回abc得到的就是基波电流负载电流减去基波电流剩下的就是谐波补偿指令。LPF的选择非常影响最终效果。我试过一阶低通截止频率放到100Hz指令会有明显相位滞后补偿电流跟不上谐波的变化。后来换了二阶Butterworth滞后稍微好一点但仍要配合超前补偿。如果检测通道的相位滞后补不过来后面电流环调得再好整体补偿效果也会打折所以检测环节的相位特性一定要重视。更极端的做法是直接用高通滤波器提取id和iq里的交流纹波但高通截止频率不能选太低否则5次分量也会被削掉。我在最终模型里用的是100Hz二阶低通加一阶超前补偿动态和精度兼顾。2.4 重复控制环节的离散实现重复控制在Simulink里没有现成模块必须自己搭。最直观的方式是用N个Unit Delay串成一条延迟链Nfs/f1比如采样频率10kHz、电网频率50Hz时N200。这个方法逻辑上没错但在模型里堆200个delay模块会非常臃肿跑一次仿真还慢。我第一版就这么干过后来实在受不了全换成MATLAB Function实现。用MATLAB Function实现重复控制的核心是维护一个长度为N的循环数组来当误差记忆。每次调用先把当前误差写入当前索引再输出延迟N步后的误差值并做Q因子回代。伪代码框架可以写成这样function y repControl(e, kr, q, N) persistent buf idx if isempty(buf) buf zeros(N, 1); idx 0; end idx idx 1; if idx N idx 1; end y kr * (e - q * buf(idx)); buf(idx) e;这段代码看起来不长坑全在细节里。buf的数据类型一定要是double如果被Simulink推断成uint负误差会被截断成错误数值重复控制就直接失灵。另外S(z)相位补偿环节我建议单独用离散滤波器模块搭不要塞进同一个MATLAB Function里这样方便单独开波德图检查相位补偿是否到位。3. 参数整定与调试顺序3.1 电流环PI参数的工程估算电流环PI参数我按工程设计法来。被控对象是滤波电感L和线路电阻R开环传函近似为PI乘以1/(LsR)。把PI的零点设置成和对象极点对消即Ki/KpR/L然后根据期望电流环带宽ωc确定KpLωc。采样和PWM逆变器带来的等效延时按1.5个采样周期估算所以电流环带宽上限一般取开关频率的1/10到1/5。按前面给的参数算一下L2mH、R0.05Ω、开关频率10kHz取带宽ωc2π×800 rad/s算出来Kp≈10Ki≈251。这个值可以直接填进离散PID控制器但要注意仿真里Kp不要取到这个理论值的100%。我实测Kp超过15之后补偿电流波形上就开始出现高频抖动这就是带宽延伸到开关频率附近离散控制的一拍延时开始捣乱了。取计算值的60%到80%往往更稳。3.2 重复控制的四个关键参数重复控制的参数可以拆成四块来调。N最没有争议就是采样频率除基波频率10kHz对50Hz就是200。Q(z)的取值直接影响高频稳定性可以取常数0.95也可以用低通滤波器。常数0.95实现简单内模对高频误差还有一定记忆低通滤波器能进一步压低高频增益让系统在高频段更稳但代价是对高次谐波的修正能力变弱。我测试下来Q取0.95常数已经够用。Kr是学习增益决定误差收敛速度。我的做法是0.2起步每增加0.1跑一次仿真看误差波形和THD的变化。Kr太小收敛慢Kr太大高频增益被放大误差波形上会出现基波周期倍数的“鼓包”严重时和开关纹波相互作用产生振荡。S(z)是补偿器用来补偿被控对象和PWM环节的相位滞后。通常做成一个低通滤波器串联若干拍超前超前拍数根据主要谐波频率处的相位滞后量来确定。我习惯在Simulink里搭一个开环测试点从控制输出到补偿电流的传递函数用波德图扫一遍然后选超前拍数让相位在3kHz以下尽量平直。这一步没法一蹴而就只能反复试但效果是立竿见影的。3.3 推荐的调试顺序与观察点调整个模型最忌一上来全开。我固定下来的顺序分四步。第一步断开重复控制只保留PI电流环先给一个直流指令验证主电路和电流环方向确保补偿电流能跟上直流量。第二步给电流环一个叠加了5次、7次分量的谐波指令单独看PI的跟踪误差记录稳态误差的幅值。第三步打开重复控制观察误差波形是不是逐周期收敛如果三个周期内没有明显下降说明N或者Q出了问题。第四步接入真实的谐波检测模块看整体THD。每个阶段都在关键位置设好Scope探针尤其是误差信号它是判断重复控制有没有生效的最直观窗口。我还会在每个阶段用Powergui的FFT工具记录THD和5、7次分量一旦数值突然恶化立刻回退到上一步用二分法定位变量。这套顺序定下来之后整个模型的调试可控性会好很多。4. 仿真结果分析与效果对比4.1 补偿前后典型波形仿真里先不看APF的电网电流波形是明显削顶的一个周期内有两个宽大的“驼峰”被切掉FFT分析出来5次谐波占了十多个百分点7次谐波也不小。接上APF但只用PI的时候波形主峰已经回来了整体接近正弦但波峰根部还是能看到一圈毛刺FFT里7次谐波残余明显5次压下去了7次又冒出来典型的此消彼长。打开PI加重复控制之后波形干净了很多毛刺基本消失FFT里各次特征谐波都被压到很低的水平。用Scope同时看负载电流、谐波指令、补偿电流三条曲线会看到补偿电流和指令几乎完全重合误差曲线在启动后的前几个周期内有一条清晰的衰减过程。这种亲眼看到误差逐周收敛的感觉比看一堆理论推导要有说服力得多。4.2 PI与PI重复控制的定量对比我直接给一组典型仿真值不同参数下具体数字会有差异但趋势是一致的。关键看稳态指标指标仅PIPI重复控制PCC电流THD8.6%2.1%5次谐波占比6.2%1.1%7次谐波占比3.8%0.5%11次谐波占比1.6%0.3%动态响应1ms内响应PI主导响应重复控制周期修正单独PI的控制效果其实也不算差从30%左右的背景畸变压到8.6%在工业上已经有应用价值了。但叠加重复控制以后THD直接掉到2.1%这背后就是重复控制把PI留下来那些“周期性残差”逐个吃掉的结果。对于需要满足IEEE 519这类谐波限值的场合这个差距是决定性的。4.3 动态响应表现动态部分看负载突变。我在仿真里跑到0.2s时把负载电阻从10Ω突降到5Ω谐波指令的幅值立刻跳变。单看补偿电流误差PI通道会在突变时刻马上输出一个大尖峰把新增误差先压下去然后在接下来的一个基波周期里重复控制通道的误差记忆逐渐更新每过一个周期误差就变小一点最终回到稳态零点。整个过程电网电流只是轻微波动一下没有出现大幅畸变。有一个细节值得拿出来说如果只看最终的THD指标你可能以为整套控制从头到尾都是重复控制在起作用实际上负载突变那十几个毫秒里真正救场的是PI。把PI通道和重复控制通道的输出分别接Scope就能看得很清楚一个出尖峰一个平滑修正分工非常明确。5. 常见问题与排查技巧5.1 发散与振荡类问题仿真一启动电流就爆表先别急着怀疑控制参数大概率是主电路的问题。我遇到最多的是补偿电流方向接反或者直流侧电容极性反了。把Universal Bridge的端口极性、电流传感器箭头、闭环反馈符号三者对齐多数发散问题能直接解决。还有PWM死区时间设得过大也会在开关瞬间造成补偿电流毛刺用默认死区值先跑通模型再说。高频振荡的问题优先怀疑重复控制的Kr和Q。Kr过大或者Q太接近1重复控制通道在高频段增益就会膨胀误差波形上出现密集的小锯齿。处理办法很直接Kr降到0.5左右Q改成0.95再跑一次看振荡有没有减弱。还有一种隐蔽情况是仿真步长设置成了变步长开关纹波和重复控制的离散记忆互相干扰改成固定步长能解决大半稳定性问题。5.2 收敛慢与精度不足类问题重复控制不收敛或者收敛特别慢第一件事就是核对N。N必须和实际基频精确匹配采样频率10kHz配50Hz就是200如果模型里把电网频率设成60Hz却仍然用N200重复控制的记忆周期就和实际信号周期对不上误差永远收敛不了。还有S(z)补偿不合适时误差波形上会残留一个固定形状的误差无论如何都不消失这就是补偿相位没对齐。排查这类问题的技巧是单频测试法。把谐波检测模块摘掉直接给电流环一个固定幅值的5次谐波指令单独开重复控制观察误差是不是按周期递减。如果单频指令能收敛但复杂谐波指令不行问题多半在S(z)没有覆盖高次谐波的相位区间。另外注意MATLAB Function里的buf数组如果被推断成整型或者初始状态没处理也会出现看似“重复控制没用到”的假象。5.3 仿真效率和模型优化200级delay链跑起来确实很慢尤其要仿真几十个周期看收敛过程时急性子根本等不了。把重复控制换成MATLAB Function加循环数组之后模型运行速度能快一个量级。还有一个优化技巧把整个仿真步长固定为50微秒PWM载波频率保持10kHz谐波检测的LPF截止频率稍微调整在保证精度的前提下能节省不少时间。Scope不要挂太多尤其不要每路信号都加To WorkspaceMATLAB工作区会被大数组占满。另外Powergui里的FFT工具分析窗口要选对。我习惯分析时取5到10个完整周期窗函数选Hanning这样谐波谱更干净。如果只看一个周期FFT分辨率不够5次和7次分量会糊在一起。调试阶段可以先把FFT工具的刷新频率调低等方案定型以后再精细测量。最后分享一点实操体会。这个模型我前前后后调了两周半最深的感受不是仿真本身有多难而是很多看似不起眼的细节决定成败比如补偿电流的极性、延迟数组的数据类型、Q值取0.95而不是0.9、FFT分析窗选对没有。每差一点THD就天差地别。如果你也在做APF或者并网逆变器的谐波抑制建议先把PI和重复控制各自的角色想清楚再动手搭模型先跑通单PI再加重复控制一步步来。等你能看到误差波形在一个周期一个周期地收敛时这套复合策略的原理就真正变成你自己的了。