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=