00000
7
10250
92800
16300
00665
8
10250
46648
08500
01086
9
10400
00000
07164
02705
同课本上给出的潮流相比较，结果完全一致，证明计算过程与程序编写正确。
最后得出迭代收敛后各支路的功率和功率损耗：
ijPijQijIijPjiQjiIji
PL
QL
120409402289045720406803869056390002601579
130307000103029950305401654034300001701551
240843201131085450866200838084840023001969
360594601346060200608201807061460013503153
450763800080074470759001070075460004801150
560241002430033680241800312023620000902118
910716402705073630716402392073630000000312
741630000665159161630000918159160000001583
860850001086083600850001496083600000000410
三、源程序及注释
由于计算流程比较简单，所以编写程序过程中没有采用模块化的形式，直接按顺序一步步进行。disp【节点数】
1xlsreadi
putxlsA3A3节点数disp【支路数】
xlsreadi
putxlsB3B3支路数disp【精度】Accuracyxlsreadi
putxlsB4B4精度bra
chxlsreadi
putxlsE4K12
odexlsreadi
putxlsM4S12Data_B1bra
ch支路参数Data_B2
ode节点参数T1zeros
2
fT2zeros
13
isqrt1
formatshort
forj1

T1j1Data_B1j3Data_B1j41i
T1j2Data_B1j51i
e
d
forj1
1
T2j1Data_B2j1Data_B2j21i
T2j2Data_B2j3Data_B2j41i
e
d
B1zeros
6
B2zeros
15
forj1

B1j1Data_B1j1
B1j2Data_B1j2
B1j3T1j1
B1j4T1j2
B1j5Data_B1j6
B1j6Data_B1j7
e
d
forj1
1
B2j1T2j1
B2j2T2j2
B2j3Data_B2j5
B2j4Data_B2j6
B2j5Data_B2j7
e
d
disp【支路参数矩阵】B1显示支路参数矩阵disp【节点参数矩阵】B2显示节点参数矩阵以上为从excel中导入初值的程序
Yzeros
1
fori1

ifB1i60
不含变压器的支路
pB1i1
qB1i2
YpqYpq1B1i3
YqpYpq
YppYpp1B1i305B1i4
YqqYqq1B1i305B1i4
else
含有变压器的支路
pB1i1
fqB1i2YpqYpq1B1i3B1i5YqpYpqYppYpp1B1i3YqqYqq1B1i52B1i3e
de
ddisp【导纳矩阵】Y显示导纳矩阵m0fori1
1ifB2i52mm1e
de
dmPQ节点个数l0fori1
1ifB2i51ll1e
de
dlPV节点个数Mismatch_powerzeroslm21fori1
11Pj0forj1
1
PjPjB2i3B2j3realYijcosB2i4B2j4imagYijsi
B2i4B2j4
e
dMismatch_poweri1realB2i1realB2i2Pje
dfork
1lm2
Qj0forj1
1
QjQjB2k
113B2j3realYk
11jsi
B2k
114B2j4r