基于matlab的潮流计算PQ分解法4节点

clc;
y=0;                     %输入原始数据，求节点导纳矩阵
y（12）=1/（0.17+0.33i）;
y（14）=1/（0.04+0.45i）;
y（23）=1/（1.1*0.3i）;
y（24）=1/（0.25+0.55i）;
y（15）=1/（-34i）;
y（25）=1/（-29i）+0.1/（1.1*0.3i）;
y（35）=-0.1/（1.1*1.1*0.3i）;
y（45）=1/（-30i）;
for i=1:4
for j=i:4
y（ji）=y（ij）;
end
end
Y=0;
                          %求互导纳
for i=1:4
for j=1:4
if i~=j
Y（ij）=-y（ij）;
end
end
end
                          %求自导纳
for i=1:4
Y（ii）=sum（y（i:））;
end
Y                         %Y为导纳矩阵
G=real（Y）;
B=imag（Y）;
for i=1:3                 %B‘
for j=1:3   
B1（ij）=B（ij）;
end
end

for i=1:2                    %B‘‘
for j=1:2   
B2（ij）=B（ij）;
end
end
 B1
 B2
 %原始节点功率
S（1）=-0.32-0.2i;         %12PQ  3PV   4平衡
S（2）=-0.55-0.3i;
S（3）=0.50;
S（4）=0;
P=real（S）;
Q=imag（S）;
%赋初值
U=