三相不平衡潮流计算、三相潮流 Matlab，编程。 ①方法：前推回代 ②说明：三相节点负荷、三相线路自阻抗和互阻抗参数。 ③适用范围：配电网，附图为程序在IEEE 33 bus节点系统中的应用。
配电网里的三相不平衡问题就像家里三个熊孩子总闹别扭，谁都觉得自己该多分点零食。传统单相潮流计算在这儿根本镇不住场子，得用三相不平衡潮流计算才能把各相电压电流安排明白。今天咱们用Matlab撸一套前推回代法的三相潮流程序，手把手教线路互阻抗怎么调教熊孩子。
先看三相线路阻抗矩阵，这玩意儿像个俄罗斯套娃：
% 三相线路阻抗矩阵初始化 Z_self = 0.4 + 1i*0.8; % 自阻抗 Z_mutual = 0.2 + 1i*0.3; % 互阻抗 Zabc = [Z_self Z_mutual Z_mutual; Z_mutual Z_self Z_mutual; Z_mutual Z_mutual Z_self];别被这个对称矩阵唬住，其实就是对角线上自阻抗，其他位置填互阻抗。这结构让三相电流计算时会产生耦合效应，跟串羊肉串似的相互拉扯。
三相不平衡潮流计算、三相潮流 Matlab，编程。 ①方法：前推回代 ②说明：三相节点负荷、三相线路自阻抗和互阻抗参数。 ③适用范围：配电网，附图为程序在IEEE 33 bus节点系统中的应用。
前推回代的精髓在于先装孙子后当爷。先假设末端电压都是理想值（比如1.0∠0°），然后从末端倒推回电源：
% 前推计算支路电流 for k = length(branch):-1:1 I_load = conj(S_load(k,:)./V_node(k,:)); % 节点负荷电流 I_branch(k,:) = I_load + sum(I_branch(child_branches,:), 1); end % 回代更新节点电压 V_source = [1.05, 1.05*exp(-1i*120*pi/180), 1.05*exp(1i*120*pi/180)]; % 电源电压 for k = 1:length(branch) downstream_node = branch(k).to; V_node(downstream_node,:) = V_node(branch(k).from,:) - Zabc * I_branch(k,:).'; end这里有个坑要注意：负荷功率S_load得按三相分开存储，每相都是复数数组。回代时用矩阵乘法处理阻抗耦合，比单相计算多出来的计算量全在这儿了。
在IEEE33节点系统实测时，收敛判据得设得宽松点。三相系统就像三个醉汉走直线，总得允许他们稍微晃悠几下：
max_error = max(abs(abs(V_node) - V_previous)); if max_error < 1e-4 % 收敛阈值放宽到1e-4 break; end运行结果会看到某相电压可能掉到0.92pu，而另外两相还坚挺在0.98pu以上——这就是典型的三相不平衡症状。程序输出的相电压热力图里，颜色分布跟迷彩服似的，哪块颜色突变就是需要加无功补偿的位置。
遇到振荡不收敛的情况，八成是互阻抗参数填反了。把Zabc矩阵右上角的互阻抗值调小点，相当于给闹得最凶的熊孩子塞个奶嘴，系统立马老实。记住，配电网计算不需要追求极致精度，能抓住主要矛盾才是正经事。
