%%IEEE33節點潮流計算程序%%
countnum=0;
BranchData=[
????1?2?0.0922?0.0470;
????2?3?0.4930?0.2511;
????3?4?0.3660?0.1864;
????4?5?0.3811?0.1941;
????5?6?0.8190?0.7070;
????6?7?0.1872?0.6188;
????7?8?0.7114?0.2351;
????8?9?1.0300?0.7400;
????9?10?1.0440?0.7400;
????10?11?0.1966?0.0650;
????11?12?0.3744?0.1238;
????12?13?1.4680?1.1550;
????13?14?0.5416?0.7129;
????14?15?0.5910?0.5260;
????15?16?0.7463?0.5450;
????16?17?1.2890?1.7210;
????17?18?0.7320?0.5
????2?19?0.1640?0.1565;
????19?20?1.5042?1.3554;
????20?21?0.4095?0.4784;
????21?22?0.7089?0.9373;
????3?23?0.4512?0.3083;
????23?24?0.8980?0.7091;
????24?25?0.8960?0.7011;
????6?26?0.2030?0.1034;
????26?27?0.2842?0.1447;
????27?28?1.0590?0.9337;
????28?29?0.8042?0.7006;
????29?30?0.5075?0.2585;
????30?31?0.9744?0.9630;
????31?32?0.3105?0.3619;
????32?33?0.3410?0.5302;];%支路阻抗
NodeData=[
????2?100.00?60.00;
????3?90.00?40.00;
????4?120.00?80.00;
????5?60.00?30.00;
????6?60.00?20.00;
????7?200.00?100.00;
????8?200.00?100.00;
????9?60.00?20.00;
????10?60.00?20.00;
????11?45.00?30.00;
????12?60.00?35.00;
????13?60.00?35.00;
????14?120.00?80.00;
????15?300.00?200.00;
????16?60.00?20.00;
????17?60.00?20.00;
????18?90.00?40.00;
????19?90.00?40.00;
????20?90.00?40.00;
????21?90.00?40.00;
????22?90.00?40.00;
????23?90.00?50.00;
????24?420.00?200.00;
????25?300.00?200.00;
????26?60.00?25.00;
????27?60.00?25.00;
????28?60.00?20.00;
????29?120.00?70.00;
????30?200.00?600.00;
????31?150.00?70.00;
????32?210.00?100.00;
????33?300.00?40.00;];%節點，負荷
pv=[15;25];pvg=2;hk=0;Qx=[-100?-500;-100?-500];pvv=[1;1];%PV節點
pii=29;piig=1;piii=0.0633;Qk=0;?????????????????????????%PI節點
pqv=33;pqvg=1;pqx=[0.8?0.5?50];????????????????????????%PQ（V）節點
dv=0;
UB=12.66;%電壓基準kV
SB=10;%功率基準MVA
ZB=UB^2/SB;%阻抗基準0hm
BranchData（:[34]）=BranchData（:[34]）/ZB;%阻抗標幺化
NodeData（:[23]）=NodeData（:[23]）/SB/1000;%功率標幺化
Qx=Qx/SB/1000;
NN=33;??%節點數
A0=zeros（NN）;
for?n=1:NN-1
????A0（BranchData（n1）BranchData（n2））=1;
end?%形成A0

X=zeros（pvg）;
for?n=1:pvg??%形成PV節點的電抗矩陣
????s=pv（n）;
????while?s>1
????????X（nn）=X（nn）+BranchData（s-14）;%求自導納
????????s=BranchData（s-11）;
????end;
end;
for?m=1:pvg
????for?