三相并联型有源电力滤波器APF仿真探索
三相并联型有源电力滤波器APF仿真（电压外环电流内环均为PI控制），id-iq谐波检测方法，SVPWM调制方法。
在电力系统领域，谐波问题一直是影响电能质量的关键因素。三相并联型有源电力滤波器（APF）作为一种有效的谐波治理手段，备受关注。今天咱们就来唠唠基于电压外环电流内环均为PI控制，搭配id - iq谐波检测方法以及SVPWM调制方法的APF仿真。
id - iq谐波检测方法
id - iq谐波检测法是一种常用的检测三相电路谐波和无功电流的方法。它基于三相电路的瞬时功率理论，通过坐标变换将三相电流从abc坐标系变换到dq坐标系。
咱们来看段简单代码示例（以Python为例，实际工程可能用C、Matlab等语言）：
import numpy as np # 假设已经获取到三相电流ia, ib, ic ia = np.array([1, 2, 3, 4, 5]) ib = np.array([2, 3, 4, 5, 6]) ic = np.array([3, 4, 5, 6, 7]) # 坐标变换矩阵C32 C32 = np.array([[2 / 3, -1 / 3, -1 / 3], [1 / np.sqrt(3), 1 / np.sqrt(3), -2 / np.sqrt(3)], [0.5, 0.5, 0.5]]) # 将三相电流变换到αβ坐标系 i_alpha_beta = np.dot(C32, np.array([ia, ib, ic])) # 再变换到dq坐标系 theta = np.pi / 4 # 假设角度 P_dq = np.array([[np.cos(theta), np.sin(theta)], [-np.sin(theta), np.cos(theta)]]) i_dq = np.dot(P_dq, i_alpha_beta[:2, :]) # 提取出直流分量，假设这里就是基波分量 i_dc_d = np.mean(i_dq[0, :]) i_dc_q = np.mean(i_dq[1, :]) # 反变换回去得到基波电流 i_dq_fundamental = np.array([[i_dc_d], [i_dc_q]]) i_alpha_beta_fundamental = np.dot(np.linalg.inv(P_dq), i_dq_fundamental) i_abc_fundamental = np.dot(np.linalg.inv(C32), np.vstack((i_alpha_beta_fundamental, np.zeros(len(ia))))) # 谐波电流就是原电流减去基波电流 i_abc_harmonic = np.array([ia, ib, ic]) - i_abc_fundamental代码分析：首先我们假设获取了三相电流ia、ib、ic。然后通过C32矩阵将三相电流变换到αβ坐标系，接着再利用角度theta通过P_dq矩阵变换到dq坐标系。在dq坐标系中提取直流分量作为基波分量，再通过反变换得到三相基波电流，最后原电流减去基波电流就得到了谐波电流。
电压外环电流内环均为PI控制
PI控制是一种经典的控制策略。在APF中，电压外环主要负责维持直流侧电压的稳定，电流内环则用于跟踪补偿电流的指令信号。
下面用伪代码简单示意下PI控制：
# 初始化参数 kp_voltage = 0.5 ki_voltage = 0.1 kp_current = 0.3 ki_current = 0.05 error_voltage_last = 0 error_current_last = 0 integral_voltage = 0 integral_current = 0 # 假设获取到直流侧电压Udc和参考电压Udc_ref Udc = 300 Udc_ref = 350 # 假设获取到电流指令值i_ref和实际电流i i_ref = 10 i = 8 # 电压外环PI控制 error_voltage = Udc_ref - Udc integral_voltage = integral_voltage + error_voltage output_voltage_pi = kp_voltage * error_voltage + ki_voltage * integral_voltage # 电流内环PI控制 error_current = i_ref - i integral_current = integral_current + error_current output_current_pi = kp_current * error_current + ki_current * integral_current代码分析：在这段伪代码里，我们先初始化了电压环和电流环的比例系数kp与积分系数ki，还有上一时刻的误差和积分值。然后获取直流侧电压和参考电压计算电压误差，通过PI控制公式得出电压环输出。同样的方法，对电流指令值和实际电流计算误差并得出电流环输出。
SVPWM调制方法
SVPWM（空间矢量脉宽调制）调制方法的目的是使逆变器输出的电压尽量接近正弦波。它通过控制逆变器开关的通断时间，在空间上合成所需的电压矢量。
% 假设已经有参考电压矢量Vref Vref = [1 + 1j; 2 + 2j; 3 + 3j]; % 假设为三相参考电压矢量 T = 0.001; % 采样周期 theta = 0; % 初始角度 for n = 1:length(Vref) % 计算参考电压矢量的幅值和角度 Vmag = abs(Vref(n)); theta = angle(Vref(n)); % 计算扇区 sector = floor(theta / (pi / 3)) + 1; % 计算时间 T1 = (sqrt(3) * Vmag * sin(pi / 3 - mod(theta, pi / 3))) / (2 * Vdc); T2 = (sqrt(3) * Vmag * sin(mod(theta, pi / 3))) / (2 * Vdc); T0 = T - T1 - T2; % 计算占空比 D1 = T1 / T; D2 = T2 / T; D0 = T0 / T; % 根据扇区确定开关状态 switch sector case 1 % 这里确定具体的开关状态逻辑，比如S1 = 1, S2 = 1, S3 = 0 等 case 2 % 类似确定开关状态 % 其他扇区同理 end end代码分析：在Matlab代码中，我们先假设已经得到了三相参考电压矢量Vref。在循环中，每次根据参考电压矢量计算其幅值、角度，进而确定所在扇区。然后根据扇区计算各个基本电压矢量的作用时间T1、T2、T0，从而得到占空比D1、D2、D0。最后根据扇区确定逆变器的开关状态，以实现SVPWM调制。
通过这一套组合拳，即id - iq谐波检测方法检测谐波，电压外环电流内环PI控制调节补偿电流，SVPWM调制方法输出期望电压，三相并联型有源电力滤波器就能有效地对电力系统中的谐波进行治理啦。这其中每一个环节都有其精妙之处，希望本文能让你对APF仿真有更深入的理解。