提到有源电力滤波器APF的谐波抑制搞电力电子的同行应该都不陌生。这几年电网里非线性负载越来越多变频器、充电桩、整流设备到处都是谐波污染问题从工厂配电房一路蔓延到楼宇和园区传统的无源滤波器PPF体积大、容易谐振、补偿效果还固定越来越顶不住。于是有源滤波器成了研究热点而怎么把它的电流跟踪控制做好直接决定补偿效果的上限。我之前在Matlab/Simulink里完整做过一版“PI重复控制”并联型有源滤波器的谐波抑制策略模型从主电路搭建到控制算法离散化再到波形分析和参数整定踩了不少坑也积累了一套稳定可复现的调参思路。这篇就围绕这个项目的完整实现过程把控制策略的设计逻辑、Simulink建模细节、波形调试方法还有我自己总结的避坑经验一次说清楚。无论你是刚接触APF控制的研究生还是想在企业项目里快速验证算法可行性的工程师这篇文章应该都能帮你省下不少折腾时间。1. 控制策略选型为什么是“PI重复控制”而不是单一控制器1.1 单一PI控制的局限在哪里并联型有源滤波器的核心任务只有一个让变流器输出一个与负载谐波电流大小相等、方向相反的补偿电流把电网侧电流“修”回正弦。这个任务落到控制层面本质上就是一个高精度的电流跟踪问题。麻烦的是需要跟踪的目标并不是单一的直流量而是包含了5、7、11、13次等一系列特定频率的谐波分量而且这些分量会随着负载工况动态变化。如果只用传统的PI控制器情况会有点尴尬。PI控制器在同步旋转坐标系下对直流指令信号有无静差跟踪能力但APF的补偿指令是交流量直接用PI去跟幅值和相位都会出现稳态误差。要知道谐波补偿最讲究的就是相位精度——哪怕只差几度补偿电流和负载谐波电流叠加之后不仅消不掉还可能放大谐波。我见过不少初学的人在这里翻车仿真波形一看补偿效果还行但THD始终降不到理想水平问题往往就出在PI对交流信号跟踪的固有相位滞后上。1.2 重复控制为什么适合周期性谐波抑制重复控制的理论基础是内模原理——系统若要无静差地跟踪某个参考输入那么这个输入的数学模型必须包含在稳定的闭环控制器内部。谐波电流是周期性的它的模型本质上就是一个周期延迟环节。重复控制器的做法很简单也很巧妙把上一周期的误差记忆下来在本周期提前修正输出等效于将误差逐周期累积、逐周期消除最终实现对周期信号的无静差跟踪。放到APF这个场景里重复控制的优势特别明显。电网基波频率是50Hz谐波在工频周期内严格重复5次谐波周期是20ms的1/57次是1/7但它们都在20ms这个基波周期上完成整数个循环。只要把重复控制器的延迟时间设置为一个基波周期20ms就能同时对所有整数次谐波产生增益这正是APF需要的宽频带谐波抑制能力。而PI控制只能针对某一个频率点做精细调节要覆盖所有谐波频率光靠一个PI是不够的。1.3 复合控制结构各取所长既然PI擅长动态响应、重复控制擅长稳态精度最合理的选择自然是把两者组合起来。我在项目中采用的是并联复合结构PI通道和重复控制通道并行两者的输出相加作为电流环的控制指令。这个结构的好处体现在两个方面。负载突变时PI通道能立刻产生响应保证补偿电流快速跟上指令变化稳态阶段重复控制通道逐步积累误差修正量把各个谐波频率点的跟踪误差压到极小而PI通道主要负责维持系统稳定性和提供基础带宽。两者分工明确兼顾了动态和稳态性能这也是目前工程上APF电流控制比较主流和成熟的方案。提示如果只用重复控制系统在负载突变的前几个周期补偿效果会比较差因为需要时间“学习”新的误差模式。调参时最好先把PI通道的带宽做高再叠加重复控制会省去很多麻烦。2. 谐波检测与指令电流提取决定补偿效果的天花板2.1 基于瞬时无功功率理论的检测方法控制策略再先进如果谐波检测不准一切都是空谈。目前工程上用得最多的谐波检测方法是基于瞬时无功功率理论pq法或ip-iq法它的核心思想是把三相电流变换到两相同步旋转坐标系然后利用低通滤波器分离出基波正序分量再用总电流减去基波分量得到谐波指令。具体到Simulink里的实现我采用的是ip-iq法它比传统的pq法在电网电压畸变情况下具有更好的适应性因为它的鉴相环节只依赖电网电压的相位信息对电压幅值波动不敏感。检测流程是这样的采集三相负载电流 ia、ib、ic用锁相环PLL获取电网电压的相位角 θ并生成 sinθ 和 cosθ通过 Park 变换把三相电流变换到 d-q 坐标系得到 ip 和 iqip 和 iq 分别经过低通滤波器滤除交流谐波分量得到直流分量 ip_dc 和 iq_dc这两个直流分量对应基波正序有功和无功电流用 Park 逆变换把直流分量变换回三相静止坐标系得到基波正序电流 iaf、ibf、icf用实际负载电流减去基波正序电流就是需要补偿的谐波电流指令 iah、ibh、ich。按这个流程最后得到的就是APF需要跟踪的谐波指令电流。需要注意这个检测方法得到的是“谐波无功”的指令如果只想补偿谐波不补无功在d-q坐标系里只需要提取ip通道的直流分量iq通道直接当作零即可。2.2 低通滤波器截止频率怎么选谐波检测环节里关键技术参数是低通滤波器的截止频率。截止频率选高了滤不干净交流分量基波提取结果里会混入谐波导致指令电流不准选低了动态响应变慢负载突变时检测延迟大补偿效果会拖后腿。我在项目中经过多轮对比最终取的是截止频率 50Hz配合二阶巴特沃斯滤波器在检测精度和响应速度之间取得了比较好的平衡。实测下来用一阶低通滤波器的动态响应快但衰减特性不够陡基波提取精度略差二阶巴特沃斯虽然加入一定延迟但幅频特性平坦、相频特性线性对基波分量的还原度更高对谐波分量的衰减更彻底。我建议优先选二阶追求检测实时性再降阶。2.3 锁相环PLL的建模细节PLL是整个检测环节的“眼睛”它的输出相位直接作用于Park变换相位不准则检测结果全盘皆错。Simulink的Simscape Electrical库里有现成的PLL模块可以直接用但需要把它的动态响应调一调。默认参数的PLL锁定时间偏长会影响整个检测环节的动态性能我最终把PLL的带宽调到了大约200Hz左右在快速锁定和抗电压扰动之间折中。自己搭PLL的话结构不复杂三相电压经过Clark变换得到αβ分量然后通过Park变换得到d-q分量q轴分量经过PI调节器输出角频率修正量再积分得到相位角。PI参数按典型二阶系统整定即可阻尼比取0.707自然角频率取1000rad/s左右锁定时间大约在几个毫秒级对于谐波检测来说足够了。参数项推荐值/范围说明PLL带宽150-300 Hz过高会引入电网噪声过低调速慢低通滤波器类型二阶巴特沃斯幅频特性平坦利于基波精确提取低通截止频率40-80 Hz兼顾检测精度与动态响应Park变换角度源PLL输出θ不要直接用理想电源相位否则失去实用意义3. PI重复控制器设计与离散化实现3.1 重复控制器的内模与完整传递函数重复控制器的经典结构我在第1部分已经提到过它的内模形式在离散域里写作[ G_{rc}(z) \frac{z^{-N}}{1 - Q(z) \cdot z^{-N}} ]其中 N 是一个基波周期内的采样点数。假设采样频率 fs10kHz基波频率 f050Hz那么 Nfs/f0200。z^-N 这一项实现了“延迟一个基波周期”的存储功能Q(z) 是为了增强系统稳定性和鲁棒性而加入的衰减因子通常取一个小于1的常数比如0.95~0.98也可以取一个低通滤波器。完整的重复控制系统还需要在重复控制器前面串联补偿器用来补偿被控对象的幅频和相频特性我这里用了超前环节加增益的组合。整个重复控制器包含以下环节周期延迟存储环节 z^-N衰减因子 Q(z)重复控制增益 Kr相位超前补偿环节 z^k被控对象逆模型的近似补偿 S(z)。S(z) 的设计比较关键它的作用是校正电流环被控对象的幅频和相频特性。LCL滤波器的高频谐振峰如果不做校正直接套重复控制很容易在高频段产生不稳定。我在项目里用二阶低通滤波器加陷波器的方式实现 S(z)把谐振峰压掉同时在整个谐波频段保持近似单位增益。3.2 离散化方法与采样频率选择控制器最终是在数字处理器里运行的所以Simulink模型里需要把连续域设计转换成离散域。这里有一个原则性问题Simulink里的连续仿真模块和离散控制模块混用时会引入计算步长问题我一般直接把整个控制回路建成离散系统这样更贴近实际DSP运行的逻辑也能避免仿真速度和精度的矛盾。离散化方法我采用的是零阶保持器ZOH法也就是在控制器的输出端加一个Zero-Order Hold把连续对象近似为离散对象。PI控制器本身可以直接用增量式数字PI公式u(k) u(k-1) Kp*[e(k)-e(k-1)] KiTse(k)在Simulink里用离散积分器Discrete-Time Integrator或者直接写S-Function实现都行。采样频率的选择直接决定重复控制器的周期延迟点数N。我在项目里用的是10kHz采样N200。如果采样频率提高到20kHzN400重复控制的存储点数翻倍对内存的需求也翻倍但动态响应会更快高频段的补偿效果也会更好。需要注意的是采样频率越高数字控制延迟占控制周期的比例越小对系统稳定性越有利但对硬件的要求也越高。仿真阶段建议先用10kHz跑通整个模型再逐步提高采样频率验证效果。3.3 PI与重复控制的复合连接方式PI控制器和重复控制器的复合不是简单把两个控制器并联接上就行这里有一个重要的连接细节重复控制器的输入应该是电流误差而不是直接并联在PI前面或后面。我最终采用的结构是电流误差 e(k) i_ref(k) - i_f(k)指令电流与反馈电流之差e(k) 分别输入PI控制器和重复控制器的误差输入端PI控制器输出 u_pi(k)重复控制器输出 u_rc(k)总控制输出 u(k) u_pi(k) u_rc(k)。这种并联复合结构的好处前面说过了动态靠PI、稳态靠重复。但要注意一点PI控制器如果参数过强会和重复控制器的修正作用打架导致系统在特定频率点出现振荡。解决办法是PI参数不要追求过强的动态特性Kp和Ki取中等偏保守的值把输出波形精细度交给重复控制去弥补。另一种复合方式是串联结构误差先经过重复控制再进入PI控制。这种结构在高频段的抗干扰能力稍强但参数整定复杂度更高两个控制器的参数会互相耦合调起来更费劲。我自己的经验是并联结构更直观调试效率更高推荐初次搭建时使用。3.4 重复控制的关键参数整定方法重复控制器的参数整定是项目的核心难点这里我把自己的整定过程拆成三步分享出来第一步确定N和Q(z)。N由采样频率和基波频率唯一确定Nfs/f0。Q(z)我取的是常数0.98这个值越大稳态误差越小但对模型参数误差的容忍度越低取0.95虽然省心但稳态精度会略受影响。建议先取0.98如果仿真出现不稳定迹象再降。第二步设计补偿器S(z)和相位超前环节。这一步的目的是让重复控制器带宽范围内即各次谐波频率点的环路增益保持较高且相位接近零。我在项目里是先把电流环开环bode图画出来观察被控对象在3.8kHz附近的谐振峰然后用一个陷波器把它削平再用低通滤波器把5kHz以上的高频增益压下去。相位超前环节我取了z^2也就是超前两个采样周期用来补偿数字控制延迟和LC滤波器的相位滞后。第三步调节重复控制增益Kr。Kr的取值一般从0.1开始逐步增加观察仿真波形中误差收敛速度和稳定裕度。Kr太小重复控制的修正强度不够稳态谐波残留大Kr太大系统容易在高频段振荡。我在项目里最终取的是0.5在仿真模型下表现比较理想。如果系统在全负载范围内都稳定可以考虑把Kr加到0.6~0.7进一步压低谐波但要小心重载工况下的稳定性。4. 主电路与控制系统Simulink建模实操4.1 仿真主电路拓扑与参数配置项目里用的并联型APF主电路是最经典的三相两电平电压源变流器VSC通过LCL滤波器和电网并联。非线性负载用三相不控整流桥带阻感负载来模拟这是最常见的谐波源。具体的仿真参数如下参数项数值说明电网线电压有效值380V50Hz三相交流电网频率50Hz基波周期20ms非线性负载三相不控整流桥RL负载整流桥直流侧电阻20Ω、电感10mHAPF直流侧电压800V电压外环给定值直流侧电容5000μF支撑直流电压、抑制电压波动LCL滤波器逆变器侧电感L10.5mH靠近变流器一侧LCL滤波器网侧电感L20.1mH靠近电网一侧LCL滤波器滤波电容Cf10μF与阻尼电阻串联阻尼电阻2Ω抑制LCL谐振峰采样频率10kHz控制周期100μs仿真步长1μs满足LCL高频特性仿真需求这张表的参数配置是按照常见APF工程样机的比例缩放的直接跑仿真效果就挺理想。有一点需要特别提醒LCL滤波器的谐振频率一定要设计在远离开关频率和主要谐波频率范围之外我用上面的参数算下来谐振频率大约是3.8kHz左右落在主要补偿频段2kHz以内和开关频段10kHz之间是比较合适的。4.2 电流内环控制模型搭建步骤Simulink模型里电流内环的控制链路是从谐波检测模块输出指令电流经过PI重复控制器再得到调制信号最后经过PWM产生IGBT驱动脉冲。搭建过程我习惯按以下顺序先搭谐波检测模块从三相负载电流提取谐波指令电流再搭电流误差计算模块指令电流减去APF输出补偿电流得到误差信号搭建PI控制器用离散积分器实现并联一个比例通路搭建重复控制器核心是N200的延迟线Delay模块参数设200配合Q(z)增益0.98、Kr增益0.5、相位超前z^2、补偿器S(z)两者输出相加进入PWM调制模块。PWM调制我采用的是双极性SPWM方式。相比SVPWMSPWM在Simulink里搭建起来更直接调制逻辑清晰对波形调试更友好后期如果想进一步提升直流电压利用率可以换成SVPWM但控制逻辑部分要额外增加扇区判断和时间计算复杂度会上升不少。4.3 直流侧电压外环控制设计APF的直流侧电压不是恒定不变的有功功率双向流动会导致直流电容充放电电压会波动。直流电压太低变流器输出电压指令无法满足补偿需求谐波补偿能力直接下降直流电压太高开关管电压应力增大器件容易损坏。因此需要一个电压外环把直流侧电压稳定在设定值附近。电压外环的设计相对简单直流电压参考值和实际值相减经过一个PI控制器输出一个有功功率指令调整量叠加到d轴电流指令上。这个d轴电流指令经过Park逆变换后会改变变流器从电网吸收的有功功率从而维持直流电压的平衡。整个外环的响应速度要设计得比电流内环慢很多才能保证多环系统稳定。我用的是电压外环PI参数 Kp0.8 Ki50带宽大约在40Hz左右明显比电流内环的带宽低实测直流电压响应平稳负载突变时电压跌落也控制在10V以内。关于电压外环的采样我在实际建模中发现一个细节直流电压纹波里含有大量的6k±1次纹波分量因为三相整流负载产生6脉动功率脉动如果直接把含纹波的电压信号送进PI控制器PI输出会叠加谐波干扰电流指令。解决办法是在电压反馈路径上串联一个截止频率约80Hz的低通滤波器滤掉纹波后再送PI控制质量会明显提升。4.4 S-Function方式实现控制器的注意事项热搜词里体现了S-Function的高关注度确实有一部分原因就是大家想把控制算法写成C代码风格的S-Function模块这样更容易向实际DSP工程移植。我的项目里PI控制器和重复控制器都用的是Simulink标准模块搭建但重复控制器的S-Function实现我也单独验证过一遍效果等价只是可读性稍差。用S-Function实现重复控制器时最核心的挑战是状态变量管理。重复控制器的内部状态是一个长度为N的数组用来存储过去一个周期的误差或输出值。在S-Function里这个数组要作为连续/离散状态变量保存我的建议是在mdlInitializeSizes里声明为DWork向量大小N200在mdlUpdate函数里完成移位更新在mdlOutputs里根据当前误差和历史状态计算输出。配速比整体仿真步长设置一致时仿真结果和标准模块搭建的完全一致。提示S-Function里的Delay数组索引计算方法要格外小心。因为MATLAB数组索引从1开始而数字控制系统里延迟N个周期对应的索引偏移是N建议用一个计数器记录当前周期数然后用mod(周期数, N)获取写入位置和读取位置避免繁琐的数组搬移操作。5. 仿真结果分析与波形解读5.1 补偿前后电网电流对比把模型跑起来之后最先要观察的是电网电流在APF投入前后的对比变化。APF未投入时三相不控整流桥负载导致的电网电流波形象一个带尖顶的扁波谐波含量非常丰富我用FFT分析过主要成分是5次约23%、7次约11%、11次约7%、13次约5%THD能达到28%左右。APF投入并稳定运行之后大概3~5个工频周期后电网电流波形会明显趋向正弦。5次谐波从23%压到1%以下7次压到0.8%以下THD从28%降到3%以内这个水平已经满足IEEE-519的谐波标准要求了。重复控制在稳态阶段的补偿精度优势特别明显如果只用PI控制单独跑THD大概只能降到8%左右5次谐波还有4%~5%的残留叠加重复控制后那些残留谐波会被逐周期“啃掉”最终压得很干净。5.2 直流侧电压波形与动态响应分析直流侧电压波形是判断控制系统是否正常的另一个关键指标。理想情况下APF稳定运行时直流电压应该稳定在800V附近纹波电压在±5V以内。如果纹波偏大说明电压外环带宽太高或者电流内环补偿效果不佳导致有功波动耦合到直流侧。负载突变试验也是仿真里非常重要的验证环节。我在仿真中设置了0.15s时负载电阻从20Ω切换到15Ω观察系统动态过程。负载突变瞬间指令电流突然增大PI通道立即响应补偿电流在1~2个周期内跟上指令直流电压经过一个短暂的跌落幅度约8V然后快速恢复整个过程大约经过80ms重新回到稳态。重复控制通道则在动态过程的几个周期内重新“学习”新的谐波模式大约经过5个工频周期后补偿精度恢复到稳态水平。这个动态响应速度对大多数工业场景是足够的。5.3 THD与补偿效果的量化对比为了直观展示控制策略的改进效果我把相同工况下“无补偿”、“仅PI控制”、“PI重复控制”三种情况做了量化对比工况5次谐波含量7次谐波含量THD无补偿23.1%11.2%28.4%仅PI控制4.6%2.8%8.2%PI重复控制0.7%0.5%2.6%这个表格清楚地说明了重复控制的“最后一公里”价值。PI控制把THD从28%压到8%左右距离标准要求仍然遥远而叠加重复控制后THD直接压到2.6%补偿精度提升了一个数量级。在很多实际项目中能否达到3%以下的THD直接决定方案能不能通过验收所以重复控制这种“稳态精度利器”是非常值得投入时间去实现的。6. 常见问题与调试验证经验6.1 Simulink仿真报错与排查思路Simulink仿真APF模型时有一些高频报错我整理几个最常见的情况和对应的排查思路。代数环报错Algebraic Loop这是搭建电流内环时最容易踩的坑。当控制器的输出直接反馈到输入端而且中间没有存储器或传递延迟时Simulink会检测到代数环仿真速度会急剧下降甚至报错。解决办法是在反馈回路里串联一个单位延迟Unit Delay或者在控制器输出端加一个Memory模块打断代数环。这样做的代价是引入一个采样周期的延迟但只要系统有足够的相位裕度影响可以忽略。Simulink提示“Caused by: LAPACK加载错误: mllapack.dll”这个问题在仿真较大模型时偶尔会出现尤其是模型里包含多个快速开关器件的时候。这通常和MATLAB版本、本机BLAS库配置有关并不是模型本身的问题。我遇到过几次处理方法是重新启动MATLAB并清理工作区必要时用rehash toolboxcache命令刷新工具箱缓存。如果还不行检查一下模型里是不是用了过大的矩阵运算有些自定义函数会触发LAPACK的边界情况换成基础的Simulink模块表达可以规避。仿真步长太大导致波形异常APF模型里含有LCL滤波器和IGBT开关开关频率10kHz的情况下如果仿真步长设成10μsLCL谐振峰和高频开关行为都无法精确还原。我的经验是最小仿真步长设在1μs固定步长不要用可变步长否则在开关动作时刻容易触发步长骤减仿真速度反而更慢。6.2 重复控制不稳定的常见原因重复控制器在某些参数组合下会出现高频振荡波形上表现为每个周期末端的纹波越来越大最终发散。我总结了三个最常见的原因第一补偿器S(z)没有把谐振峰压住。前面说过LCL滤波器在3.8kHz附近有一个谐振峰如果S(z)设计时没把这段增益压下来重复控制在这个频率点的增益就会很高闭环极点越过单位圆必然振荡。解决办法是在S(z)里加入与谐振频率对应的陷波器把谐振峰增益压到0dB以下。第二相位超前量不足。重复控制在每个谐波频率点的相位裕度是保证稳定性的关键。数字控制延迟、PWM零阶保持效应、LCL滤波器的相位滞后都要靠z^k去补偿。如果k取得太小高频段的相位滞后累积会造成振荡如果k取得太大中频段的相位会过冲同样会不稳定。我在项目里k2是按照10kHz采样频率下大约40μs的数字控制延迟补偿来设计的直接套用多数情况问题不大。第三Kr取值过大。这个前面已经提过Kr调大会加速误差收敛但稳定裕度会下降。调试时应该从Kr0.1起逐步加观察波形不要一步到位。如果发现在某个人为整定值上振荡可以先退回上一步的值再继续找边界参数。6.3 控制器参数调试顺序的建议从项目开始到最终调通我经历了不少参数组合总结出一套比较高效的调试顺序供参考在只接PI控制器的情况下先把电流环调稳用阶跃响应看跟踪速度和超调接入重复控制器但把Kr调到0.1的小值确认系统稳定、无异常振荡逐步增大Kr观察FFT频谱中高次谐波的幅值是否持续下降若高频段出现抬升趋势说明相位补偿不足或S(z)没压住回头调整S(z)和z^k重复控制稳态波形满意后再回头微调PI参数以不引发系统振荡为上限最后做负载突变、电压波动等鲁棒性测试确保在非标称工况下依然稳定。按这个顺序整个仿真模型的调试时间一般可以控制在几天到一周左右前提是模型搭建本身没有结构性错误。6.4 仿真与实验的差异预判跑仿真的时候容易产生一个错觉算法效果这么好拿下实验肯定也没问题。实际从仿真到实验的鸿沟比想象中大得多。仿真中默认三相电网电压完全对称、无畸变、无频率波动而实际电网里这些理想条件几乎不存在。重复控制依赖精确的基波周期一旦电网频率偏离50Hz比如49.8Hz固定周期延迟N200的重复控制器就会出现性能下降。工程上通常需要设计一个频率自适应机制根据PLL测出的实际基波频率动态调整延迟点数或者内插权重这正是仿真模型可以直接扩展的方向。另外仿真里IGBT开关管的损耗、死区效应和管压降都是理想化的实际运行时这些非线性因素会产生额外的低次谐波。要更接近实际情况可以在Simulink里加入IGBT导通压降和死区模型观察控制策略对这些扰动的抑制能力。这些升级方向都能在现有模型基础上逐渐迭代不必推翻重来。这套基于PI重复控制的APF谐波抑制模型我在多个负载工况下反复验证过稳态补偿精度、动态响应能力都表现稳定。它的价值不只是把THD压到3%以下更在于控制结构清晰、参数整定有章法无论是学习研究还是工程验证都是一个扎实的起点。我在实现过程中最大的感触是谐波抑制这块控制策略的先进性是显式的但真正决定最终指标的往往是被很多人忽略的细节——低通滤波器阶数、LCL谐振峰的处理、数字控制延迟的补偿、离散化方法的取舍。这些细节一个个抠完模型才能真正做到“能看又能用”。如果你也在折腾APF仿真按我这条路线搭一遍应该能更早看到一条平滑的电网电流波形出现在Scope里。
