1. 为什么要做这个模型APF谐波抑制的背景与方案选型1.1 谐波问题的现实场景先聊几句题外话。这几年做并网变流器、储能PCS、充电桩电源的朋友应该都有一个明显的感受电网的谐波环境越来越差。大量非线性负载接入整流器、变频器、开关电源、LED驱动每个都是谐波源。我见过工厂配电房里的电流波形THD能到30%以上中性线电流比相线还大变压器嗡嗡响、发热严重。这种环境下单纯靠无源滤波器已经很难兼顾全频段补偿体积大、容易谐振、对电网阻抗变化敏感于是有源电力滤波器APF成了主流选择。APF的核心思想就一句话你产生多少谐波我反向注入多少用主动抵消的方式把电网侧电流掰回正弦。这个反向注入听起来简单实际做起来是典型的强耦合、非线性、多频段跟踪控制问题。谐波次数从2次到50次都有幅值有大有小相位随时在变控制器的跟踪带宽和稳态精度必须同时拉满。我这次搭的模型就是一套基于PI重复控制的APF谐波抑制策略仿真。整个模型在Simulink里实现覆盖了谐波提取、电流内环控制、PWM调制、主电路、负载扰动全链路。它解决的核心问题是单纯PI控制器补偿多频次谐波时存在稳态误差而单纯重复控制动态响应又偏慢二者结合恰好取长补短。适合正在做APF课题的学生、刚入门的电源工程师以及想快速验证控制算法效果的研发人员。1.2 核心矛盾在哪里APF从电路拓扑上看就是一个并网变流器但和常规并网逆变器有一个非常本质的区别并网逆变器只输出基波电流用一个PI就能把基波稳态误差压得很低APF输出的电流里面同时包含几十个频率分量而且每个频率的幅值相位都被谐波提取环节实时刷新。这时候如果只用PI控制情况会变成什么PI控制器在旋转坐标系下对直流信号有无限增益所以基波跟踪没问题。但是谐波在旋转坐标系下是交流信号比如5次谐波转到dq系是300Hz扰动PI的增益在这个频率上已经衰减得很厉害稳态压缩比根本不够。你加大PI增益系统带宽上去了但采样延迟和PWM调制延迟带来的相位滞后也跟着放大很容易引发谐振整个系统可能直接啸叫。那为什么不能只用重复控制重复控制基于内模原理在谐波频率处提供很高的增益理论上能实现多频次谐波的无静差跟踪。但重复控制本质上是逐周期的误差修正它天生带一个周期的延迟动态性能差。负载突变的时候电流要过好多个基波周期才能重新收敛这在很多工业现场是没法接受的。所以PI重复控制是一个很自然的互补结构。PI管动态把基波和低频段快速稳住重复控制管稳态把周期性谐波逐个压下去。两者并联工作的时候PI负责快速逼近重复控制负责精确消除配合得当的话既能把THD降下来又能在负载切换时保持相对快的响应速度。这也是这个模型最值得拆开细看的地方。1.3 方案选型背后的逻辑我在搭建这套模型之前其实是先做了一轮方案对比的。备选方案包括单纯PI、PI谐振控制器组合PR、无差拍控制、重复控制。逐一说下取舍。PR控制器的思路是在特定频率点叠加谐振增益比如5、7、11、13次各放一个谐振器稳态效果确实好但有两个麻烦。第一你需要明确知道要补偿哪些频次电网频率波动或者负载特性变化导致谐波频次偏移时谐振点的选择性反而成了短板。第二谐振器数量多参数整定维度爆炸实际调试工作量很大。无差拍控制动态响应快但对模型参数极敏感电感漂移、开关管压降、死区效应都会导致预测误差在实验平台上往往表现为高频噪声明显、鲁棒性差。做仿真很理想上硬件就露馅。重复控制不需要逐个频点设计它对所有基波整数倍频率都提供高增益结构统一、参数少特别适合谐波频谱不确定的场景。唯一的问题就是动态响应慢。那为什么不直接堆一个PIPR重复控制的豪华组合因为控制器的复杂度上去了系统稳定性分析、参数鲁棒性验证、离散化实现的难度都会指数上升。做工程不是把能加的都加上而是在满足指标的前提下选一个够用、好调、可维护的结构。PI重复控制在我个人的实践中是性价比最高的组合。2. 谐波提取与控制框架拆解2.1 瞬时无功功率理论谐波电流怎么拿出来的任何APF都绕不开一个问题你需要补偿的参考电流从哪来早期做法是频域FFT分解逐次提取谐波分量再合成算法复杂度高、实时性差而且对电网频率波动敏感。现在工程上几乎都用瞬时无功功率理论也就是基于dq坐标变换或p-q变换的方法。我在这套模型里用的是基于dq旋转坐标变换的提取方式。三相电流经Clark变换从abc系得到αβ静止坐标系分量再经过Park变换转到dq旋转坐标系。这里有一个关键细节d轴和q轴的旋转角度必须锁相于电网A相电压一般用SRF-PLL来保证也就是说dq变换的旋转速度严格跟随电网基波频率。进了dq坐标系之后基波正序分量变成直流量谐波分量变成交流量。接下来把直流量滤掉剩下的全是谐波。最简单的做法是接一个低通滤波器但低通滤波器的截止频率和阶数会导致基波提取有相位延迟进而影响谐波补偿精度。所以更稳的做法是用陷波器或者直接用高通通路也就是全通减去低通把基波成分减掉得到谐波指令。模型里我用的是陷波器加延迟补偿的组合实测下来比单纯一阶低通的效果好很多。还有一点容易被忽略dq坐标系下同时包含了正序谐波和负序谐波。正序分量的次数在dq系里表现为6k±1次变换后的6k次分量负序分量同样会表现成高频扰动。所以提取环节出来的是一个综合的谐波dq信号后面控制器要的是把这个dq信号跟踪好而不是把某个特定频次单独抠出来。2.2 电流内环的控制对象APF的电流控制对象是从逆变器输出电压到并网电流的传递函数。忽略电网电压扰动和开关死区后在L型滤波器的情况下传递函数就是1/(LsR)。这里L是并网电感R是寄生电阻加线路阻抗。但要注意实际模型里我并没有直接把这个理想传递函数丢给控制器设计而是加了几个现实的约束。第一PWM调制有一个等效延迟一般等效为1.5个开关周期的纯延迟这会显著降低相位裕度。第二数字控制至少再加入一个采样周期的计算延迟。这两个延迟叠加起来在20kHz开关频率下意味着约75微秒的相位滞后在600Hz的25次谐波处就已经有超过16度的相位衰耗了如果不做补偿重复控制的稳定性会非常紧张。所以我在电流内环里加了有源阻尼本质是在控制输出里注入一个与电感电流微分项相关的阻尼项。更工程化的实现是用电容电流反馈或者电感电流比例反馈来等效增加系统阻尼。这里我用的是在控制器输出前叠加高通滤波阻尼简单说就是给谐振峰压一压。没有这一步只是单纯调PI增益去压稳态误差系统很容易在某些谐波频段出现增益尖峰然后整个电流波形开始振。2.3 整体控制框图整条链路的信号走向我用文字描述一下电网三相电压经PLL得到θ负载电流经Clark/Park得到dq电流dq电流经高通/陷波提取得到谐波dq参考值APF并网电流同样做Clark/Park得到dq反馈谐波dq参考值与dq反馈的偏差送入PI重复控制并联结构控制器输出经反Park变换得到αβ电压指令αβ电压指令送入PWM调制模块生成开关信号开关信号驱动三相逆变桥逆变桥通过L滤波器并网完成谐波注入。还有一层是直流母线电压控制。APF要和电网交换有功功率必须有直流母线电压支撑母线电压不能掉。常规做法是取dq电流的d轴基波分量做电压外环把母线电压误差经PI控制后叠加到d轴电流参考上。但在APF里要注意这个d轴参考是基波有功分量不能跟谐波提取出来的谐波dq参考混在一起。我在模型里是把这两个通道分开的母线电压环输出加在谐波dq参考之前经过坐标变换后自然和基波对应而谐波通道只负责无功和谐波电流这样相互干扰最小。3. 控制器设计从理论到参数3.1 内模原理与重复控制的本质先讲明白重复控制为什么能对周期性误差做到无静差。内模原理是控制理论的一个经典结论如果控制系统开环传递函数中包含参考信号或扰动信号的动态模型那么这个系统对该信号就能实现无静差跟踪。比如PI控制器里面其实内嵌了一个积分器1/s它恰好就是阶跃信号的模型所以PI对阶跃输入没有稳态误差。谐波扰动是周期性的它的模型在连续域是一个周期延迟环节。重复控制把这个周期延迟环节放进控制器内部对基波周期的延时取正反馈等效于在每一个基波整数倍频率上都放了一个谐振器。这个思想非常优雅你不需要知道扰动包含哪几次谐波只要它是周期性的重复控制就自动覆盖所有频点。传递函数上可以这样理解。重复控制器中有一个延时代Z^-NN为一个基波周期内的采样点数配合一个正反馈环就构成了周期信号发生器。误差信号进入这条支路后周期部分会被记住并在下一个周期继续输出补偿信号非周期部分无法通过所以它相当于一个梳状滤波器形式的积分器。3.2 数字化实现要点离散化是重复控制的落地关键。我在模型里采用了一拍延迟型重复控制结构带一个低通滤波器Q(z)。这个Q(z)是稳定性的核心。理论上Q(z)可以直接取1这样谐波频段增益无穷大稳态性能最好但实际不行。原因在于对象模型在高频段存在未建模动态和谐振峰1的Q会让系统在高频段失去增益裕度极易振荡。工程上Q(z)一般取一个低截止频率的低通滤波器比如0.95或者一阶低通牺牲少量高频谐波补偿能力换整个回路的稳定性。第二个关键是相位补偿。重复控制环里必须设计一个相位超前补偿环节来抵消从控制输出到电流反馈之间的采样延迟、PWM延迟、滤波器相位滞后。这个补偿环节一般用超前校正或者直接移相配合一个增益Krc来调节补偿强度。我见过的很多模型跑出来波形很差八成就是这个相位补偿没做对导致每个谐波频率处的补偿相位对不上误差不但没消掉反而被励起来。第三个关键是周期延迟N的计算。N等于基波周期除以采样周期必须取整。50Hz基波、20kHz采样下N400。这里有个容易踩的坑如果Simulink里仿真步长不固定或者PLL输出的频率有小幅波动N对应的延迟会偏离真实基波周期重复控制的效果立刻劣化。所以这个模型里重复控制这块的触发脉冲最好也同步到PLL周期脉冲上而不是简单地用固定采样步数。3.3 PI与重复控制的并联结构PI和重复控制并联不是简单把两个控制器输出相加就完事。要考虑两个控制器的频带分工和动态重叠问题。我的结构是这样的误差信号e同时进入PI支路和重复控制支路两支路输出相加得到总控制量。PI支路带宽相对低主要针对基波和谐波包络的快速变化重复控制支路主要针对整周期重复的稳态误差。但在PI带宽覆盖到的频段内两个支路的相位特性不同可能在某几个谐波频点出现输出相互削弱的情况。要避免这个问题有两个常用手段。一是把重复控制的误差通路加一个高通或者带通滤波器让重复控制主要补偿高于PI带宽的频段。二是在重复控制支路输出前加一个与PI支路增益匹配的比例系数让两支路的幅值合理分配。我在模型里选了方案二并在仿真里观察了支路输出波形确保重复控制支路的补偿电压没有和PI支路明显掐架。还有一个细节启动阶段。重复控制有周期延迟启动瞬间它输出的补偿信号是零此时全靠PI压误差如果负载谐波很大PI可能在大误差下出现饱和。所以我在模型里给重复控制支路的参考输入加了斜坡使能前两个基波周期逐渐切入。这个做法在实验结果里非常有用能明显降低启动冲击。4. Simulink仿真模型搭建步骤4.1 主电路与系统参数设定模型整体在Simulink里用电气模块搭建我用的是Simscape Electrical专业库同时保留了一个纯数学建模的简化版本作对比验证。主电路参数我按一个10kW级别的APF来设计。电网线电压有效值380V频率50Hz电网阻抗设为每相0.1mH加0.05Ω。APF并网电感L取2mH寄生电阻R取0.1Ω。直流母线电容取4mF母线电压参考值设750V。逆变桥用三相两电平桥开关频率20kHz直流侧用一个恒定直流源模拟已整定的母线方便先验证控制算法再在后面的模型里接入母线电容和电压环。负载侧我搭了一个三相不可控整流桥带阻感负载输出端还并了一个整流型负载产生的谐波特征以5、7、11、13次为主。之所以模拟这个负载是因为工业现场最常见的谐波源就是这种整流负荷做出来的仿真结果有代表性。所有电气模块的采样时间统一设置控制部分用离散控制器采样频率20kHz。整个模型的主电路仿真步长我设成固定步长1e-6秒控制器的采样在零阶保持器里以2e-5秒同步触发。仿真时间设0.4秒前0.2秒让系统进入稳态0.2秒时做一个负载切换方便观察动态响应。4.2 谐波提取模块搭建谐波提取模块我用的是abc到dq的同步旋转坐标变换。锁相环我用的是Simulink里自带的SRF-PLL但做了两处改造。一处是把PI参数按捕获带和锁定带分别整定避免启动阶段锁相震荡另一处是PLL输出角度经过一个wrap到0到2π的限幅再做周期脉冲生成这个周期脉冲直接提供给重复控制器做周期对齐用。dq变换之后基波分量在d轴和q轴上表现为直流。为了把谐波分离出来我用了陷波器。陷波器的中心频率设在100Hz也就是dq坐标系下的2倍基波频率。之所以用100Hz而不是直接用低通是因为dq坐标系下最低的谐波扰动频率就是100Hz附近的2次谐波分量陷波器可以把这一带滤得干净又不影响更高频段的谐波。陷波器的带宽我设得比较窄品质因数在2左右这样滤基波的同时不会把相邻谐波成分过度衰减。提取出来的谐波dq参考值会经过一个限幅和低通平滑防止负载电流快速变化时参考电流出现尖峰避免过调制。这个平滑滤波器的截止频率我设成2kHz对5次300Hz和7次200Hz谐波几乎没有影响但可以把PWM谐波和采样噪声挡掉。4.3 电流控制器与PWM调制实现电流控制器分三块PI控制器、重复控制器、合成限幅模块。PI控制器在dq坐标系下运行分d轴和q轴各一个参数我先按典型工程法初定然后在仿真里微调。需要强调的是这个PI控制器的带宽不需要太高我把它压到大约800Hz左右也就是只覆盖到13次谐波附近。带宽更低的好处是相位裕度充足不容易震荡同时给重复控制留出足够的高频工作空间。重复控制器的实现方式比较关键。我在Simulink里用离散延时节实现Z^-NN400配一个正反馈环前向通路里串Q(z)低通滤波器。Q(z)我用的是离散一阶低通截止频率取开关频率的五分之一也就是4kHz。在Q(z)之后是相位补偿环节我用的是超前校正相位超前量设25度对应补偿两拍延迟。增益Krc设在0.6左右这个值不能太激进大了高频会振小了谐波压不住。最后合成控制量经过反Park变换得到αβ轴电压指令。PWM调制我用的是规则采样法设定载波频率20kHz。逆变桥开关信号的生成不直接用Simulink自带PWM模块而是自己做三角载波比较好处是可以很方便地加入死区时间并且能精确控制调制波更新时刻这对验证延迟补偿算法有帮助。4.4 参数整定与仿真配置要点参数整定我强烈建议分步走不要一上来就全参数同时跑。我的顺序是这样第一步先把重复控制支路旁路掉只用PI控制跑通整个系统。调PI的时候观察两个指标一个是启动瞬间的动态超调一个是基波电流的稳态误差。调试中我发现PI的Kp主要决定动态响应速度Ki决定稳态收敛但INR穿过某值之后电流波形开始出现60Hz附近的振荡那就是带宽拉太高了要退回来。第二步接上重复控制支路但Krc先取0.2观察稳态波形。此时THD应该已经开始改善如果波形反而变差先查相位补偿方向是不是反了再查Q(z)截止频率是不是太高。我遇到过一次选4kHz正好选10kHz就开始振的现象后面细看就是高频相位损耗过大。第三步逐步加大Krc到0.6同时观察电网电流FFT。正常情况应该是5次、7次谐波先降下来然后11次、13次跟着降。如果高频段谐波不降反升说明相位补偿在高频不够需要增加超前角度或调整Q(z)。仿真配置上有几个细节要留意。第一Simscape电气模块的求解器必须配合固定步长否则控制器的离散延迟和实际物理模型之间会产生不可控的时序错位。第二PLL的初始角度要设置合理不然启动前半个周期会出现很大的冲击电流。第三数据采集的采样率要高于控制频率我在测量THD时用的是1MHz的示波器采样率在谐波分析前先做数据对齐避免FFT频谱泄漏。5. 仿真中常见的坑与排查方法5.1 稳态效果差重复控制似乎没起作用这是最常见的现象。现象是PI控制跑出来的THD是15%接了重复控制之后变成14.5%几乎没变化。这时候先不要怀疑控制结构先检查两个地方。第一N值对不对。如果N400但周期脉冲和PLL之间存在微小偏差重复控制器的补偿信号会出现周期性相位漂移积攒起来就是误差无法收敛。验证方法很简单在重复控制器输出端接一个Scope看补偿波形是否是严格的周期函数有没有逐周期漂移。我排查过一台类似的模型发现是PLL的dq角度在跨周期接缝处有跳变导致重复控制对齐失败。第二误差通路的方向对不对。重复控制补偿的是电流误差但如果谐波提取模块的输出极性接反了控制器会朝着错误方向补偿结果谐波被放大而不是被抑制。这时候看FFT会发现某几个频点反而升高。排查方法更直接把参考电流设成零注入一个固定谐波扰动看控制器输出的相位是否和误差反向。5.2 重复控制引发系统振荡另一个典型问题是重复控制接入后系统在某个高频段出现持续振荡电流波形上叠加着一团密集的毛刺。这个毛刺频率往往在2到5kHz之间很容易被误认为是开关噪声其实它是控制回路自身的不稳定。根本原因通常是重复控制回路在未建模的高频谐振峰处增益过高。解决方法优先顺序是第一步降低Krc第二步降低Q(z)截止频率第三步增加相位补偿角度。我把这三步分别试过降低Krc是最快的止血手段但会牺牲谐波补偿力度降低Q(z)截止频率不伤中低频段但对25次以上谐波压制减弱增加相位补偿是最理想的手段但要注意超前校正在高频段会放大噪声过量了反而坏事。还有一个此前困扰过我的细节仿真模型里的电感如果带了非线性饱和模型或者电气模型里包含杂散电容高频谐振特性会和理想电感完全不同。这种时候不能只看控制器的理论设计要看实际仿真的频率响应我后来给电流环做了扫频分析从20Hz扫到10kHz直接在Bode图上看到谐振峰然后针对性地设计阻尼参数振荡问题就解决了。5.3 仿真发散与计算速度问题Simulink仿真发散很多时候不是控制算法的问题而是电气系统模型本身数值刚性太强。现象包括仿真步长减到纳秒级才能跑稳定或者某几个节点上电压出现几百千伏的尖峰。排查思路我总结为三步。第一步检查是否有代数环尤其是电气模块和控制模块直接相连而没有串入状态空间或延迟模块的位置代数环会严重恶化数值稳定性。第二步检查物理模型参数有没有数量级错配比如电感2mH碰上电容4mF特征频率很低在1e-6步长下会引入不必要的数值阻尼。第三步检查PWM调制信号是否频繁触发开关动作如果载波频率和仿真步长的倍数关系不合适会出现伪开关事件导致状态跳变。计算速度方面Simscape模型本身偏重。我的做法是把电路模型和控制模型拆分运行先用简化数学模型做控制参数的粗调再用Simscape模型做最终验证。这样可以节省大量等待时间。如果一定要跑全模型建议把无关的测量模块关掉示波器采样率不要设太高仿真数据存储用双缓冲处理完一批再存一批。5.4 THD分析优化经验THD的测量和计算也有讲究。Simulink里用FFT工具分析电流波形时要注意先去掉启动暂态段取稳态后至少5个完整基波周期并且窗函数要选汉宁窗避免频谱泄漏。我实测下来的优化经验是影响最终THD的不只是控制器的谐波补偿能力还有几个被忽视的因素。第一个是死区时间。逆变桥死区会在电流中引入3次和5次附近的分量如果不做死区补偿THD可能卡在3%左右就下不去了。我在仿真里设了2微秒死区然后加了基于电流极性的死区补偿THD从3.0%降到1.8%效果非常明显。第二个是直流母线电压纹波。母线电容不够大时母线电压波动会通过调制比进入输出电流产生额外谐波。我在4mF母线下THD是1.5%换成2mF后升到2.3%这个趋势在理论上是可预期的。第三个是谐波提取环节的相位延迟。如果陷波器滤基波的过程引入了延迟参考电流和实际谐波之间就存在相位差重复控制只能部分补偿残留误差较大。我在提取环节做了延迟补偿后THD进一步降到了1.2%左右。这套模型最终在0.2秒负载切换工况下稳态THD能做到1.4%动态恢复时间在40毫秒以内效果还是令人满意的。6. 调试经验与优化建议6.1 调参顺序与测试清单如果让我给出一个可以直接照抄的调试流程我会这样列先空载启动APF观察母线电压闭环是否稳定母线电压波动在±10V以内。投入负载但APF补偿通道关闭记录负载电流的FFT数据确认谐波特征。打开PI电流环观察dq电流跟踪效果确认基波部分无误。打开谐波提取通路接上参考电流但控制器增益先调低确认补偿方向正确。投入重复控制逐步加Krc每次加完都记录THD和波形形态。做负载切换测试确认动态恢复时间满足指标。做电网电压跌落测试确认系统不失稳。这个顺序看起来琐碎但每步都在缩小故障排查范围。我见过很多人直接全功能打开出了问题根本不知道是哪个环节引起的回头逐段排查反而更浪费时间。6.2 离散化与延迟处理的几点心得离散化对重复控制的影响我写了三年理论和仿真之后才真正重视起来。一个容易被忽略的点是重复控制的周期延迟N是整数但PLL输出周期和采样周期之间不是严格整数倍关系时会引入周期偏差。工程上可以做一个分数延迟滤波器来实现非整数延迟补偿但在仿真阶段先用周期脉冲同步的方式基本够用。另外控制器的离散化方式也要注意。PI控制器我用的是前向欧拉离散重复控制的延迟直接用整数拍延迟相位补偿则用离散超前校正。大可以放心地把这些模块都按离散时间设定因为后续如果要去硬件平台移植这些离散结构几乎可以无缝改写为C代码这是建模型时就要考虑到的。电流内环的PI参数在离散域和连续域之间转换时积分系数要乘以采样周期比例系数不需要变。很多人在这里搞混导致仿真里的参数一上硬件就发散因为离散化增益没有换算对。我在模型里特别标注了这个换算关系避免后续移植时踩坑。6.3 后续可以怎么扩展这套模型做完之后扩展空间其实很大。我列几个我自己验证过或者正在推进的方向供参考。第一个方向是把并网电感换成LCL滤波器。LCL在高频段有更好的衰减特性但谐振峰问题更突出。目前在模型里已经测过只要加上电容电流反馈有源阻尼再套用这套PI重复控制依然能稳定工作而且高频谐波压制效果比L型更好。第二个方向是自适应重复控制。电网频率在49.5到50.5Hz波动时固定N值的效果会下降。可以通过实时检测PLL频率动态调整延迟拍数N实现频率自适应的重复控制。这个方向我测试过初步版本频率偏移0.5Hz时THD劣化幅度能控制在0.5%以内。第三个方向是把谐波提取改为基于自适应陷波器的方法在电网频率波动和背景谐波变化环境下提取精度更稳定。这个方向适合后续做深入研究配合神经网络或在线参数辨识做自适应滤波也是当前学术界比较活跃的分支。一些个人经验结尾最后说点实际的。这套PI重复控制的APF模型我从最早的开环仿真模型一路调到能稳定复现预期结果中间踩过的坑都在上面的章节里。如果说有什么最想提醒同行的话那就是电路模型和控制模型之间的时序关系永远是仿真的第一优先级。Simulink里看似简单的离散控制器一旦和电气模型的计算步长没有配合好结果完全可能得出错误的结论。一定要花时间把采样时序、PWM更新时刻、PLL角度对齐这三件事搞清楚再去调控制参数。另外就是调试要留证据。每次调参前先把当前版本的波形和FFT数据存下来完成后对比这样才能知道每一步改动到底起到了什么作用而不是拍脑袋感觉好像好了点。这套习惯是做好任何控制仿真的基础。今天先聊到这里希望能给正在做APF项目的朋友省下几周摸索时间。