并联型有源电力滤波器APF三相三线模型都包括simulink仿真 利用基于瞬时无功功率理论的ip-iq谐波检测算法对三相三线制并联型APF控制系统进行建模与Matlab仿真最近在搞三相三线制并联型APF的仿真发现基于ip-iq谐波检测的方案确实挺有意思。这玩意儿主要靠实时分解电流的谐波分量来生成补偿信号咱们今天就来扒一扒这个实现过程。先说说坐标系变换这个基本功。在Simulink里搭模型的时候三相电流先得走个Clark变换把ABC坐标系转成αβ静止坐标系。不过要注意相位对齐问题我一般用这个公式来写Matlab函数块function [i_alpha, i_beta] clark_transform(ia, ib, ic) % 三相转两相变换 i_alpha (2/3)*(ia - 0.5*ib - 0.5*ic); i_beta (2/3)*(sqrt(3)/2*ib - sqrt(3)/2*ic); end这里有个坑实际硬件采样可能会有相位延迟仿真时最好在电流输入端加个零阶保持器模拟采样保持。我之前没注意这个结果谐波检测总带着奇怪的毛刺。接下来是重头戏——ip-iq运算环节。需要构造与电网电压同步的正余弦信号这里推荐用锁相环(PLL)模块。不过新手容易犯的错是直接拿标准正弦信号做同步实际电网电压可能有畸变得用增强型PLL才行。并联型有源电力滤波器APF三相三线模型都包括simulink仿真 利用基于瞬时无功功率理论的ip-iq谐波检测算法对三相三线制并联型APF控制系统进行建模与Matlab仿真谐波分离的关键在于低通滤波器设计这个直接决定补偿效果。用二阶Butterworth滤波器的话代码大概长这样function [ip_filt, iq_filt] harmonic_filter(ip, iq) persistent filter_coeff if isempty(filter_coeff) [b,a] butter(2, 100*2*pi, low, s); % 截止频率100Hz filter_coeff {b, a}; end ip_filt filter(filter_coeff{1}, filter_coeff{2}, ip); iq_filt filter(filter_coeff{2}, filter_coeff{1}, iq); end注意数字滤波会有相位延迟仿真时要对比原始信号和滤波后的波形延迟超过1/4周期就得调参。有次我把截止频率设到300Hz结果5次谐波都没滤干净补偿后THD反而更高了。补偿电流生成环节建议用电流跟踪控制这里PWM调制策略选空间矢量还是三角波比较实测发现三电平SVPWM在开关损耗和补偿精度上平衡得更好。分享个简单的占空比计算逻辑function duty svpwm_calc(v_ref, v_dc) % 归一化参考电压 v_norm v_ref / (v_dc/2); % 扇区判断 theta atan2(v_norm(2), v_norm(1)); sector floor(theta/(pi/3)) 1; % 占空比计算此处简化 duty [0.5 0.5*v_norm(1), 0.5 0.5*v_norm(2), 0.5 - 0.5*(v_norm(1)v_norm(2))]; end调试的时候记得看补偿电流的跟踪波形如果出现明显滞后可能是控制器的PI参数没调好。有个小技巧先把积分项设为零调比例系数到临界震荡状态再慢慢加积分。最后说下整体仿真架构。主电路用Universal Bridge模块模拟IGBT直流侧电容初始电压设到650V左右比较合理。负载建议用三相不控整流带阻感负载这样能产生典型的6k±1次特征谐波。跑完仿真别忘了做FFT分析看补偿前后的电流频谱对比。正常应该能把THD从30%压到5%以内。不过要注意仿真步长别设太大建议用ode23tb求解器特别是开关频率超过10kHz的时候固定步长容易数值震荡。