潮流计算程序

合集下载

电力系统潮流计算完整c语言程序(含网损计算的最终版)

for(j=0;j<Bus_Num;j++)
{
ia[i]=ia[i]+gY_G[n][j]*ge[j]-gY_B[n][j]*gf[j];
ib[i]=ib[i]+gY_G[n][j]*gf[j]+gY_B[n][j]*ge[j];
}
}
for(i=0,n=1;i<Bus_Num-1;i++,n++)
{
gDelta_PQ[2*i]=gDelta_P[i];
gDelta_PQ[2*i+1]=gDelta_Q[i];
}
if((fp=fopen("C:\\Documents and Settings\\Zorro\\桌面\\1\\data\\unbalance.txt","w"))==NULL)
if(gBus[n].Type==1)
gDelta_Q[i]=gDelta_Q[i]-gf[n]*(gY_G[n][j]*ge[j]-gY_B[n][j]*gf[j])+ge[n]*(gY_G[n][j]*gf[j]+gY_B[n][j]*ge[j]);
}
}
for(i=0;i<Bus_Num-1;i++)
{
gY_G[i][j]=0.0;
gY_B[i][j]=0.0;
}
for(l=0;l<Line_Num;l++)
{
i=gLine[l].No_I-1;
j=gLine[l].No_J-1;
r=gLine[l].R;
x=gLine[l].X;

附录_潮流计算程序

附录牛顿-拉夫逊法潮流计算程序本附录所介绍的牛顿-拉夫逊法潮流计算程序采用极坐标形式，其中所涉及的计算公式与第四章中的式(4-42)-(4-59)完全相同，计算程序框图与图4-6基本一致。

程序采用C语言编制。

为了便于初学者阅读，节点导纳矩阵和雅可比矩阵都用满阵存储而未采用稀疏技巧；但为了适当照顾使用的方便性，在输入数据中对节点编号次序不作任何要求。

下面先介绍输入文件的格式和要求，然后列出程序，最后说明潮流结果的输出。

建议读者先从原始数据的输入中了解和熟悉它们在程序中对应的变量、结构体数组及其成员的名称，然后对照图4-6和式(4-42)-(4-59)仔细和耐心地阅读程序，最好能在计算机上亲自实现和进行调试。

在调试时，可以用下面给出的对应于[例4-3]系统的输入数据，以及程序运行中得出的中间结果，逐步与[例4-3]中所给出的中间结果进行对比，从而查找错误所在并进行改正。

一、原始数据的输入程序通过“输入数据.txt”文件输入以下5个数据段。

1.信息(共6个)(1)总节点数(变量num_node)(2)线路和并联电容器总数(变量num_line)(3)变压器支路总数(变量num_tran)(4)发电机节点总数(变量num_gene)(5)负荷节点总数(变量num_load)(6)节点功率不平衡量的容许误差(变量error)2.线路和并联电容器数据(结构体数组line)：每一线路或并联电容器包括5个数据I侧节点编号和J侧节点编号可以对换；线路和并联电容器之间的次序可以任意，而且允许多条线路或多个电容器并联。

3.变压器支路数据(结构体数组tran)：每一变压器支路包括5个数据变压器电阻、电抗和非标准变比与两侧节点编号之间的关系服从图2-29,即电阻和电抗在1侧而非标准变比在2侧；三绕组变压器需按图2-26化成3个变压器支路，其中3侧变压器支路的非标准变比为1。

变压器支路之间的次序可以任意。

4.发电机节点数据(结构体数组gene)：每一发电机节点包括5个数据对于PQ节点，节点种类为1，电压可给任意值；对于PV节点，节点种类为-1，发出无功功率可给任意值；对于平衡节点，节点种类为0，发出有功功率和发出无功功率都可给任意值。

潮流计算程序

潮流程序说明书1.程序原始数据输入格式如下：B1是支路参数矩阵，第一列和第二列为节点编号，按照从小到大编写；第三列为支路串联阻抗参数；第四列为对地导纳参数；第五列为变压器变比；第六列为是否含有变压器，“1”为含有变压器，“2”为不含变压器。

B2是节点参数矩阵，第一列为节点注入发电功率参数；第二列为负荷功率参数；第三列为节点电压参数；其中“1”为平衡节点，“2”为PQ节点，“3”为PV节点。

X是节点号和对地参数矩阵，第一列为节点号；第二列为对地情况，“0”表示没有对地。

2.基本原理：采用牛顿-—拉夫逊潮流算法（直角坐标系）主要变量定义：B1—支路参数矩阵；B2—节点参数矩阵；X—节点号和对地参数；S—节点功率数值，序号为奇数的储存有功功率，偶数的储存无功功率；DS—序号为奇数的储存有功功率不平衡量，对PQ节点，序号为偶数的储存无功功率不平衡量，对PV节点，序号为偶数的储存电压不平衡量；I—节点电流数值；DetaU—节点电压修正值；Number—迭代次数；dS—线路上损耗的功率；3.程序流程框图如下：启动输入原始数据形成节点导纳矩阵Y设置节点电压初始值B2（i,3）设置迭代次数Number=1根据给定初始值计算功率数值PV 节点计算有功功率功率不平衡量DS 和电压不平衡量DetaUPQ 节点计算有功功率和无功功率不平衡量DSJacobi 矩阵是否全部形成，i>n计算Jacobi 矩阵的对角元素与非对角元素增大节点号i=i+1置节点号i=1解修正方程求出功率增量和各节点电压增量求出电压增量最大值max(DetaU)迭代是否收敛，最大值是否小于ac求出各节点功率S、各条支路首末端功率Si 、Sj,节点电压U 、节点电流I停止增大迭代次数Number=Number+1计算各节点电压值是否是否4.计算结果如下（完整计算结果见文本文档）其中：3节点运行结果如下：请输入节点数:n=3请输入支路数:n1=3请输入平衡节点号:bl=3请输入误差精度:ac=0.000001请输入支路参数：B1=[1 2 0.15i 0 1 1;2 3 0.93i 0 1 0;2 3 0.5i 0 1 0]请输入节点参数：B2=[0.9 0 1 0 0 3;0 0 1 0 0 2;0 0 0.9001 1.05 0 1]节点号和对地参数:X=[1 0;2 0;3 0]Y =0 + 6.6667i0 - 6.6667i0 + 3.0753i0 + 6.6667i0 - 9.7419i0 - 3.0753i0 0 + 3.0753iU =0.96941.00000.9001各节点的功率S为（按节点号排列）：0 - 0.1976i0 + 0.5110i0 - 0.2765i各条支路的首端功率Si为（同B1顺序一样）：0 + 0.2038i0 + 0.1074i0 + 0.1998i各条支路的末端功率Sj为（同B1顺序一样）：0 - 0.1976i0 - 0.0967i0 - 0.1798i各条支路的功率损耗dS为（同B1顺序一样）：0 + 0.0062i0 + 0.0200i0 + 0.0200i9节点运行结果如下：请输入节点数:n=9请输入支路数:n1=9请输入平衡节点号:bl=1请输入误差精度:ac=0.001请输入支路参数：B1=[1 4 0+0.0576i 0 1 1;2 7 0+0.0625i 0 1 1;3 9 0+0.0586i 0 1 1;4 5 0.01+0.085i 0.176i 1 0;4 6 0.017+0.092i 0.158i 1 0;5 7 0.032+0.161i 0.306i 1 0;6 9 0.039+0.17i 0.358i 1 0;7 0.0085+0.072i 0.149i 1 0;8 9 0.0119+0.1008i 0.209i 1 0]请输入节点参数：B2=[0 0 1.04 1.05 0 1;1.63+0i 0 1.025 1.05 0 3;0.85+0i 0 1.025 1.05 0 3;0 0 1 0 0 2;0 -1.25-0.50i 1 0 0 2;0 -0.90-0.30i 1 0 0 2;0 0 1 0 0 2;0 -1-0.35i 1 0 0 2;0 0 1 0 0 2]节点号和对地参数:X=[1 0;2 0;3 0;4 0;5 0;6 0;7 0;8 0;9 0]Y =Columns 1 through 50 -17.3611i 0 0 0 +17.3611i0 0 -16.0000i 0 0 00 0 0 -17.0648i 0 00 +17.3611i 0 0 3.3074 -39.3089i-1.3652 +11.6041i0 0 0 -1.3652 +11.6041i2.5528 -17.3382i0 0 0 -1.9422 +10.5107i0 0 +16.0000i 0 0-1.1876 + 5.9751i0 0 0 0 00 0 0 +17.0648i 0 0Columns 6 through 90 0 0 00 0 +16.0000i 0 00 0 0 0 +17.0648i0 0 0-1.9422 +10.5107i0 00 -1.1876 + 5.9751i0 0 -1.2820 + 5.5882i3.2242 -15.8409i-1.6171 +13.6980i0 2.8047 -35.4456i-1.1551 + 9.7843i2.7722 -23.3032i0 -1.6171 +13.6980i-1.2820 + 5.5882i 0 -1.1551 + 9.7843i 2.4371 -32.1539i -0.0168 -0.1328iU =1.04003.2481 - 0.9555i3.2076 - 0.8593i3.2428 - 0.5596i3.2235 - 0.7453i1.8842 - 0.3587i2.3853 - 0.6321i2.4410 - 0.8344i各节点的功率S为（按节点号排列）：1.0e+002 *0.1010 - 0.3977i-0.0361 + 0.4978i-0.0571 + 0.0252i0.0687 + 1.7360i-0.0214 + 0.1542i-0.0808 - 0.4329i0.0531 - 0.5394i-0.1579 - 0.1256i0.2074 + 0.4176i-0.00180.0071-0.00710.00410.0044各条支路的首端功率Si为（同B1顺序一样）： -1.0105e+001 +1.2945e+002i3.6128 -36.2000i5.7087 - 2.3123i7.0539 - 2.2715i9.9160 +46.4206i4.8627 +14.6652i-1.8285 -15.5465i5.7060 - 4.7477i-10.1578 -16.9569i各条支路的末端功率Sj为（同B1顺序一样）： 10.1046 -39.7723i-3.6128 +49.7840i-5.7087 + 2.5192i-7.0064 + 0.7595i-6.2529 -27.7428i-4.0131 -12.9970i4.2137 +23.3616i-5.6348 + 4.4012i10.8152 +20.7070i各条支路的功率损耗dS为（同B1顺序一样）： -0.0000 +89.6771i0 +13.5840i0 + 0.2069i3.6631 +18.6778i0.8496 + 1.6681i2.3852 + 7.8152i0.0712 - 0.3465i0.6574 + 3.7501i39节点运行结果如下：请输入节点数:n=39请输入支路数:n1=46请输入平衡节点号:bl=31请输入误差精度:ac=0.1请输入支路参数：B1=[1 2 0.00350+0.04110i 0.6987i 1 0;1 39 0.00100+0.02500i 0.7500i 1 0;2 3 0.00130+0.01510i 0.2572i 1 0;2 25 0.00700+0.00860i 0.1460i 1 0;30 2 0.00000+0.01810i 0.00 102.5 1;3 4 0.00130+0.02130i 0.2214i 1 0;3 18 0.00110+0.01330i 0.2138i 1 0;4 5 0.00080+0.01280i 0.1342i 1 0;4 14 0.00080+0.01290i 0.1382i 1 0;5 6 0.00020+0.00260i -0.02170 1 0;5 8 0.00080+0.01120i 0.1476i 1 0;6 7 0.00060+0.00920i 0.1130i 1 0;6 11 0.00070+0.00820i 0.1389i 1 0;31 6 0.00000+0.02500i 0.00 107.000 1;7 8 0.00040+0.00460i 0.0780i 1 0;8 9 0.00230+0.03630i 0.3804i 1 0;9 39 0.00100+0.02500i 1.2000i 1 0;10 11 0.00040+0.00430i 0.0729i 1 0;10 13 0.00040+0.00430i 0.0729i 1 0;32 10 0.00000+0.02000i 0.00 107.000 1;11 12 0.00160+0.04350i 0.00 100.600 1;13 12 0.00160+0.04350i 0.00 100.600 1;13 14 0.00090+0.01010i 0.1723i 1 0;14 15 0.00180+0.02170i 0.3660i 1 0;15 16 0.00090+0.00940i 0.3420i 1 0;16 17 0.00070+0.00890i 0.1342i 1 0;16 19 0.00160+0.01950i 0.3040i 1 0;16 21 0.00080+0.01350i -0.12740 1 0;16 24 0.00030+0.00590i 0.0680i 1 0;17 18 0.00070+0.00820i 0.1319i 1 0;17 27 0.00130+0.01730i 0.3216i 1 0;20 19 0.00070+0.01380i 0.00 106.000 1;33 19 0.00070+0.01420i 0.00 107.000 1;34 20 0.00090+0.01800i 0.00 100.900 1;21 22 0.00080+0.01400i 0.2565i 1 0;22 23 0.00060+0.00960i 0.1846i 1 0;35 22 0.00000+0.01430i 0.00 102.500 1;23 24 0.00220+0.03500i 0.3610i 1 0;36 23 0.00050+0.02720i 0.00 100.000 1;25 26 0.00320+0.03230i 0.5130i 1 0;37 25 0.00060+0.02320i 0.00 102.500 1;26 27 0.00140+0.01470i 0.2396i 1 0;26 28 0.00430+0.04740i 0.7802i 1 0;26 29 0.00570+0.06250i 1.0290i 1 0;28 29 0.00140+0.01510i 0.2490i 1 0;38 29 0.00080+0.01560 0.00 102.500 1;] 请输入节点参数：B2=[0 0 1 0 0 2;0 0 1 0 0 2;1 0 0 2;0 -3.22-0.024i1 0 0 2;0 -5.00-1.84i0 0 1 0 0 2;0 0 1 0 0 2;1 0 0 2;0 -2.338-0.84i1 0 0 2;0 -5.22-1.76i0 0 1 0 0 2;0 0 1 0 0 2;0 0 1 0 0 2;1 0 0 2;0 -0.085-0.88i0 0 1 0 0 2;0 0 1 0 0 2;1 0 0 2;0 -3.20-1.53i1 0 0 2;0 -3.29-0.323i0 0 1 0 0 2;0 -1.58-0.30i1 0 0 2;0 0 1 0 0 2;1 0 0 2;0 -6.80-1.03i1 0 0 2;0 -2.74-1.15i0 0 1 0 0 2;1 0 0 2;0 -2.475-0.846i0 -3.086+0.922i 1 0 0 2;1 0 0 2;0 -2.24-0.472i1 0 0 2;0 -1.39-0.17i1 0 0 2;0 -2.81-0.755i1 0 0 2;0 -2.06-0.276i0 -2.835-1.269i1 0 0 2;2.50+1.3621i 0 1.04750 0 0 3;0 0 1.050 1.05 0 1;6.50+1.7590i 0 0.98310 1.05 0 3;6.32+1.0335i 0 0.99720 1.05 0 3;5.08+1.6440i 0 1.01230 1.05 0 3;6.50+2.0484i 0 1.04930 1.05 0 3;5.60+0.9688i 0 1.06350 1.05 0 3;5.40-0.0444i 0 1.02780 1.05 0 3;8.30+0.1939i 0 1.02650 1.05 0 3;1.03 1.05 0 3]10.0+0.68460i -11.04-2.50i节点号和对地参数:X=[1 0;2 0;3 0;4 0;5 0;6 0;7 0;8 0;9 0;10 0;11 0;12 0;13 0;14 0;15 0;16 0;17 0;18 0;19 0;20 0;21 0;22 0;23 0;24 0;25 0;26 0;27 0;28 0;29 0;30 0;31 0;32 0;33 0;34 0;35 0;36 0;37 0;38 0;39 0]Y =1.0e+002 *Columns 1 through 50.0365 - 0.6337i -0.0206 + 0.2416i 0 0 00 0-0.0206 + 0.2416i-0.0566 + 0.6574i0.6465 - 1.5929i-0.0285 + 0.4677i0.1469 - 1.8684i0 -0.0566 + 0.6574i-0.04860 0 -0.0285 + 0.4677i0.1251 - 2.0157i+ 0.7782i0 0 0 -0.0486 + 0.7782i0.4061 - 5.4887i0 0 0 0-0.2941 + 3.8235i0 0 0 0 00 0 0 0-0.0635 + 0.8883i0 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 -0.0479 + 0.7722i0 0 0 0 00 0 0 0 00 0 -0.0618 + 0.7468i 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 -0.5693 + 0.6994i 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 + 0.0054i 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 0-0.0160 + 0.3994i 0 0 0 0Columns 6 through 100 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 0-0.2941 + 3.8235i 0 -0.0635 + 0.8883i 0 00.4679 - 6.1153i -0.0706 + 1.0824i 0 0 00 0-0.0706 + 1.0824i-0.1876 + 2.1576i0.2582 - 3.2390i-0.0174 + 0.2744i0.2685 - 3.3173i0 -0.1876 + 2.1576i0 0 -0.0174 + 0.2744i 0.0334 - 0.6658i0 0 0 00.4290 - 4.6106i-0.1034 + 1.2107i 0 0 0-0.2145 + 2.3056i0 0 0 0 00 0 0 0-0.2145 + 2.3056i0 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 + 0.0037i 0 0 0 00 0 0 0 0 + 0.0047i0 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 -0.0160 + 0.3994iColumns 11 through 150 0 0 0 00 0 0 0 00 0 0 -0.0479 + 0.7722i0 0 0 0 0-0.1034 + 1.2107i 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 0-0.2145 + 2.3056i 0 -0.2145 + 2.3056i 0 00.3263 - 3.7448i -0.0001 + 0.0023i 0 0 00 0-0.0001 + 0.0023i-0.0001 + 0.0023i0.0000 - 0.0000i-0.0875 + 0.9823i0 -0.0001 + 0.0023i0.3105 - 3.5163i-0.03800.1734 - 2.2088i0 0 -0.0875 + 0.9823i+ 0.4577i0 0 0 -0.0380 + 0.4577i0.1389 - 1.5083i0 0 0 0-0.1009 + 1.0542i0 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 0Columns 16 through 200 0 0 0 00 0 0 0 00 0 -0.0618 + 0.7468i 0 00 0 0 0 00 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 0-0.1009 + 1.0542i 0 0 0 00.3596 - 5.1047i0 -0.0418 + 0.5094i-0.0878 + 1.1167i-0.1034 + 1.2107i0 00.2344 - 2.8992i-0.0878 + 1.1167i0 -0.1034 + 1.2107i 0.1651 - 1.9557i 0 00 0 0.0418 - 0.5080i-0.0003-0.0418 + 0.5094i+ 0.0068i0 0 0 -0.0003 + 0.0068i0.0367 - 0.7228i-0.0437 + 0.7381i 0 0 0 00 0 0 0 00 0 0 0 0-0.0860 + 1.6905i 0 0 0 00 0 0 0 00 0 0 0 00 -0.0432 + 0.5748i 0 0 00 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 -0.0003 + 0.0066i0 0 0 0-0.0003 + 0.0055i0 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 0Columns 21 through 250 0 0 0 00 0 0 0-0.5693 + 0.6994i0 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 0-0.0437 + 0.7381i 0 0 -0.0860 + 1.6905i0 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00.0838 - 1.4488i -0.0407 + 0.7120i 0 0 0-0.0649 + 1.0376i0 00.1055 - 1.7474i-0.0407 + 0.7120i-0.0179 + 0.2846i0.0827 - 1.3195i0 -0.0649 + 1.0376i0 0 -0.0179 + 0.2846i 0.1038 - 1.9730i0 0 0 00.5997 - 1.0027i0 0 0 0-0.0304 + 0.3066i0 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 + 0.0068i 0 0 00 0 -0.0001 + 0.0037i 0 00 0 0 0-0.0001 + 0.0042i0 0 0 0 00 0 0 0 0Columns 26 through 300 0 0 0 00 0 0 0 0 + 0.0054i0 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 00 0 0 0 00 0 0 0 00 -0.0432 + 0.5748i 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 0-0.0304 + 0.3066i 0 0 0 0-0.0190 + 0.2092i-0.0145 + 0.1587i0.1280 - 1.3359i-0.0642 + 0.6742i-0.0642 + 0.6742i 0.1074 - 1.2461i 0 0 0-0.0609 + 0.6566i0 0.0799 - 0.8607i-0.0190 + 0.2092i0.0754 - 0.8089i0 -0.0609 + 0.6566i-0.0145 + 0.1587i0 0 0 0 0 - 0.5525i0 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 -0.00590 0 0 0 0Columns 31 through 350 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 + 0.0037i 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 + 0.0047i 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 -0.0003 + 0.0066i 00 0 0 -0.0003 + 0.0055i0 0 0 0 00 0 0 0 0 + 0.0068i0 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 - 0.4000i 0 0 0 00 0 - 0.5000i 0 0 00 0 0.0346 - 0.7025i 0 00 0 0 0.0277 - 0.5542i0 0 0 0 0 - 0.6993i0 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 0Columns 36 through 390 0 0 -0.0160 + 0.3994i 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0160 + 0.3994i 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 00 0 0-0.0001 + 0.0037i0 0 0 00 00 -0.0001 + 0.0042i0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0059 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 00 0 00.0068 - 0.3675i0 00 0.0111 - 0.4307i0 0 0.6098 0 0 0 0 0.0319 - 0.7890i U =1.0e+003 *0.0915 - 0.1339i0.0712 - 0.1171i0.0618 - 0.0724i0.0680 - 0.0540i0.0666 - 0.0494i0.0745 - 0.0550i0.0792 - 0.0574i0.1659 - 0.0776i0.0535 - 0.0308i0.0555 - 0.0296i3.8741 +4.1652i0.0514 - 0.0316i0.0546 - 0.0614i0.0497 - 0.1055i0.0496 - 0.1237i0.0547 - 0.1336i0.0605 - 0.1277i0.0504 - 0.1256i-0.0002 + 0.0015i0.0466 - 0.1265i0.0461 - 0.1286i0.0455 - 0.1290i0.0484 - 0.1255i0.0744 - 0.1436i0.0516 - 0.1896i0.0520 - 0.1647i0.0334 - 0.2574i0.0283 - 0.2765i0.0016 - 0.0779i0.00110.0014 - 0.0726i-0.0072 - 0.0790i-0.0062 - 0.0964i0.0011 - 0.0800i-0.0023 - 0.0823i-0.0031 - 0.0825i-4.3537 + 0.6080i0.2219 - 0.0873i各节点的功率S为（按节点号排列）： 1.0e+009 *-0.0000 + 0.0000i-0.0000 + 0.0000i -0.0000 - 0.0000i -0.0000 + 0.0000i -0.0000 + 0.0000i -0.0000 - 0.0000i -0.0000 - 0.0000i -0.0000 + 0.0000i -0.0000 + 0.0000i 0.0000 - 0.0001i 0.0002 + 0.0001i 0.0000 - 0.0001i -0.0000 + 0.0000i -0.0000 - 0.0000i -0.0000 + 0.0000i 0.0000 + 0.0000i -0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i -0.0000 - 0.0000i 0.0000 + 0.0000i -0.0000 + 0.0000i -0.0000 + 0.0000i -0.0000 + 0.0000i -0.0000 + 0.0000i -0.0000 + 0.0000i -0.0000 + 0.0000i -0.0000 - 0.0000i -0.0000 + 0.0003i 0.0000 + 0.0000i -0.0000 + 0.0003i 0.0000 + 0.0004i 0.0000 + 0.0005i -0.0000 + 0.0004i 0.0000 + 0.0002i0.0000 + 0.0003i1.1785 + 0.0007i 0.0005 + 0.0008i 0.00000.0000-0.00000.0000-0.00000.00000.0000-0.00000.0000-0.00000.0000-0.00000.0000-0.00000.0000-0.00000.0000-0.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.0000各条支路的首端功率Si为（同B1顺序一样）： 3.5400e+005 +2.4736e+005i-3.5581e+005 -2.3628e+005i1.0071e+005 +2.6013e+005i2.3682e+005 -1.6375e+005i3.7264e+003 +1.4471e+006i-8.0159e+004 +2.7957e+005i1.4664e+005 -5.0894e+004i-1.1893e+005 +8.1179e+004i-6.8733e+003 +9.6236e+004i-8.1580e+004 +1.3584e+005i-3.9332e+004 -8.2485e+004i-7.5059e+004 +2.1608e+005i-1.9406e+001 +2.7536e+005i-2.5689e+004 -1.0282e+005i-1.0682e+005 -2.1610e+005i-1.2535e+005 -4.1708e+005i-2.9483e+004 -1.4331e+004i2.6881e+004 +1.7017e+004i1.7943e+003 +1.8913e+005i2.7401e+007 +7.4282e+008i2.7399e+007 +7.4282e+008i1.3064e+005 -1.2139e+005i1.1489e+005 -1.2327e+005i7.7164e+004 -2.1257e+005i-2.8359e+004 -1.6599e+005i-1.7425e+003 -1.7302e+004i3.5700e+004 -1.7516e+004i3.8835e+004 -3.2323e+004i-1.2909e+005 +6.6903e+004i1.0067e+005 -2.4294e+005i6.7162e+004 +1.3247e+006i6.6365e+004 +1.2810e+006i-7.0330e+000 +2.1637e+002i9.7641e+003 -2.0069e+004i1.0245e+004 -4.5881e+003i2.4215e+003 +1.2980e+006i-1.4470e+004 +6.4486e+003i1.4051e+004 +6.8360e+005i1.9061e+005 -1.7744e+005i3.1770e+004 +1.1214e+006i-6.2319e+004 +3.1983e+005i1.2274e+005 -2.7770e+005i1.1898e+005 -2.7519e+005i9.8516e+004 -3.3183e+005i4.8839e+006 -7.0589e+005i各条支路的末端功率Sj为（同B1顺序一样）： -3.3804e+005 -8.4052e+004i3.5993e+005 +3.0184e+005i-9.6779e+004 -2.2023e+005i-2.1492e+005 +1.8682e+005i-3.7264e+003 -5.6637e+003i8.6098e+004 -1.8535e+005i-1.4524e+005 +6.3681e+004i1.2077e+005 -5.2863e+004i8.2090e+004 -1.2718e+005i4.0208e+004 +9.3498e+004i8.3567e+003 +9.6235e+004i8.0400e+004 -1.5426e+005i19.4058 -26.1530i2.6210e+004 +1.0809e+005i1.2059e+005 +4.2529e+005i1.3051e+005 +4.9198e+005i2.9595e+004 +1.5259e+004i-2.6774e+004 -1.6139e+004i-1.7943e+003 -1.0552e+003i-7.9694e+004 -1.7952e+004i-7.7356e+004 -1.2573e+004i-1.2279e+005 +2.0853e+005i-1.0739e+005 +2.0990e+005i-7.3846e+004 +2.4186e+005i2.9462e+004 +1.7743e+005i1.7620e+003 +1.2057e+004i-3.7913e+004 +1.8782e+004i-3.8792e+004 +3.1940e+004i1.2980e+005 -6.1207e+004i-9.6457e+004 +2.9084e+005i3.9717e+001 +1.3853e+002i-3.5171e+003 -6.0815e+003i2.4385e+001 +1.3068e+002i-9.7460e+003 +1.5660e+004i-1.0242e+004 +1.1971e+003i-2.4215e+003 -7.0055e+003i1.4506e+004 -1.2520e+004i-1.5552e+003 -3.8100e+003i-1.8260e+005 +2.4172e+005i-2.8930e+003 -4.7791e+003i6.6279e+004 -2.8646e+005i-1.1337e+005 +3.3958e+005i-1.0726e+005 +3.4408e+005i-9.6140e+004 +3.3945e+005i2.8547e+005 +7.0589e+005i各条支路的功率损耗dS为（同B1顺序一样）： 1.5956e+004 +1.6331e+005i4.1140e+003 +6.5558e+004i3.9344e+003 +3.9904e+004i2.1901e+004 +2.3079e+004i-7.8671e-011 +1.4414e+006i1.4002e+003 +1.2787e+004i 1.8394e+003 +2.8317e+004i 8.3252e+002 +1.2332e+004i 5.0953e+002 +8.6579e+003i 8.7681e+002 +1.1013e+004i 6.5898e+002 +9.2305e+003i 5.3409e+003 +6.1812e+004i -1.3866e-011 +2.7534e+005i 5.2034e+002 +5.2760e+003i 1.3774e+004 +2.0919e+005i 5.1660e+003 +7.4900e+004i 1.1258e+002 +9.2738e+002i 1.0692e+002 +8.7808e+002i 7.5033e-012 +1.8808e+005i 2.7321e+007 +7.4280e+008i 2.7322e+007 +7.4281e+008i 7.8440e+003 +8.7133e+004i 7.4944e+003 +8.6625e+004i3.3182e+003 +2.9296e+004i 1.1030e+003 +1.1435e+004i 1.9504e+001 -5.2453e+003i -2.2133e+003 +1.2653e+003i4.2504e+001 -3.8293e+002i 7.1612e+002 +5.6968e+003i 4.2119e+003 +4.7899e+004i6.7202e+004 +1.3248e+006i 6.2848e+004 +1.2749e+006i 1.7352e+001 +3.4704e+002i 1.8043e+001 -4.4090e+003i3.6386e+000 -3.3909e+003i4.4565e-011 +1.2910e+006i 3.5978e+001 -6.0715e+003i 1.2496e+004 +6.7979e+005i 8.0137e+003 +6.4281e+004i2.8877e+004 +1.1166e+006i3.9594e+003 +3.3374e+004i 9.3637e+003 +6.1880e+004i1.1719e+004 +6.8892e+004i2.3756e+003 +7.6167e+003i 5.1694e+006。

潮流计算程序及计算结果

附表1：计算机计算潮流程序：%本程序的功能是用牛顿——拉夫逊法进行潮流计算% B1矩阵：1、支路首端号；2、末端号；3、支路阻抗；4、支路对地电纳% 5、支路的变比；6、支路首端处于K侧为1，1侧为0% B2矩阵：1、该节点发电机功率；2、该节点负荷功率；3、节点电压初始值% 4、PV节点电压V的给定值；5、节点所接的无功补偿设备的容量% 6、节点分类标号clear;n=13;%input('请输入节点数：n=');nl=13;%input('请输入支路数：nl=');isb=1;%input('请输入平衡母线节点号：isb=');pr=0.00001;%input('请输入误差精度：pr=');B1=[];%input('请输入由支路参数形成的矩阵：B1=');B2=[];%input('请输入各节点参数形成的矩阵：B2=');Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1=zeros(nl); %－－－－－－－修改部分－－－－－－－－－－－－ym=0;SB=100;UB=220;%ym=input('您输入的参数是标么值？（若不是则输入一个不为零的数值）'); if ym~=0%SB=input('请输入功率基准值:SB=');%UB=input('请输入电压基准值:UB=');YB=SB./UB./UB;BB1=B1;BB2=B2;for i=1:nlB1(i,3)=B1(i,3)*YB;B1(i,4)=B1(i,4)./YB;enddisp('B1矩阵B1=');disp(B1)for i=1:nB2(i,1)=B2(i,1)./SB;B2(i,2)=B2(i,2)./SB;B2(i,3)=B2(i,3)./UB;B2(i,4)=B2(i,4)./UB;B2(i,5)=B2(i,5)./SB;enddisp('B2矩阵B2=');disp(B2)end% % %---------------------------------------------------for i=1:nl %支路数if B1(i,6)==0 %左节点处于低压侧p=B1(i,1);q=B1(i,2);elsep=B1(i,2);q=B1(i,1);endY(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5)); %非对角元Y(q,p)=Y(p,q);Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2; %对角元K侧Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2; %对角元1侧end%求导纳矩阵disp('导纳矩阵Y=');disp(Y)%----------------------------------------------------------G=real(Y);B=imag(Y); %分解出导纳阵的实部和虚部for i=1:n %给定各节点初始电压的实部和虚部e(i)=real(B2(i,3));f(i)=imag(B2(i,3));V(i)=B2(i,4); %PV节点电压给定模值endfor i=1:n %给定各节点注入功率S(i)=B2(i,1)-B2(i,2); %i节点注入功率SG-SLB(i,i)=B(i,i)+B2(i,5); %i节点无功补偿量end%=========================================================== ========P=real(S);Q=imag(S);ICT1=0;IT2=1;N0=2*n;N=N0+1;a=0;while IT2~=0IT2=0;a=a+1;for i=1:nif i~=isb %非平衡节点C(i)=0;D(i)=0;for j1=1:nC(i)=C(i)+G(i,j1)*e(j1)-B(i,j1)*f(j1);%Σ(Gij*ej-Bij*fj)D(i)=D(i)+G(i,j1)*f(j1)+B(i,j1)*e(j1);%Σ(Gij*fj+Bij*ej)endP1=C(i)*e(i)+f(i)*D(i);%节点功率P计算eiΣ(Gij*ej-Bij*fj)+fiΣ(Gij*fj+Bij*ej)Q1=C(i)*f(i)-e(i)*D(i);%节点功率Q计算fiΣ(Gij*ej-Bij*fj)-eiΣ(Gij*fj+Bij*ej)%求P',Q'V2=e(i)^2+f(i)^2; %电压模平方%========= 以下针对非PV节点来求取功率差及Jacobi矩阵元素=========if B2(i,6)~=3 %非PV节点DP=P(i)-P1; %节点有功功率差DQ=Q(i)-Q1; %节点无功功率差%=============== 以上为除平衡节点外其它节点的功率计算=================%================= 求取Jacobi矩阵===================for j1=1:nif j1~=isb&j1~=i %非平衡节点&非对角元X1=-G(i,j1)*e(i)-B(i,j1)*f(i); % dP/de=-dQ/dfX2=B(i,j1)*e(i)-G(i,j1)*f(i); % dP/df=dQ/deX3=X2; % X2=dp/df X3=dQ/deX4=-X1; % X1=dP/de X4=dQ/dfp=2*i-1;q=2*j1-1;J(p,q)=X3;J(p,N)=DQ;m=p+1;J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X4;J(m,q)=X2;elseif j1==i&j1~=isb %非平衡节点&对角元X1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);% dP/deX2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);% dP/dfX3=D(i)+B(i,i)*e(i)-G(i,i)*f(i); % dQ/deX4=-C(i)+G(i,i)*e(i)+B(i,i)*f(i);% dQ/dfp=2*i-1;q=2*j1-1;J(p,q)=X3;J(p,N)=DQ;%扩展列△Qm=p+1;J(m,q)=X1;q=q+1;J(p,q)=X4;J(m,N)=DP;%扩展列△PJ(m,q)=X2;endendelse%=============== 下面是针对PV节点来求取Jacobi矩阵的元素===========DP=P(i)-P1; % PV节点有功误差DV=V(i)^2-V2; % PV节点电压误差for j1=1:nif j1~=isb&j1~=i %非平衡节点&非对角元X1=-G(i,j1)*e(i)-B(i,j1)*f(i); % dP/deX2=B(i,j1)*e(i)-G(i,j1)*f(i); % dP/dfX5=0;X6=0;p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV;m=p+1;J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6;J(m,q)=X2;elseif j1==i&j1~=isb %非平衡节点&对角元X1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);% dP/deX2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);% dP/dfX5=-2*e(i);X6=-2*f(i);p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV;m=p+1;J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6;J(m,q)=X2;endendendendend%========= 以上为求雅可比矩阵的各个元素===================== for k=3:N0 % N0=2*n （从第三行开始，第一、二行是平衡节点）k1=k+1;N1=N; % N=N0+1 即N=2*n+1扩展列△P、△Qfor k2=k1:N1 % 扩展列△P、△QJ(k,k2)=J(k,k2)./J(k,k); % 非对角元规格化endJ(k,k)=1; % 对角元规格化if k~=3 % 不是第三行%======================================================== ====k4=k-1;for k3=3:k4 % 用k3行从第三行开始到当前行前的k4行消去for k2=k1:N1 % k3行后各行下三角元素J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);%消去运算endJ(k3,k)=0;endif k==N0break;end%==========================================for k3=k1:N0for k2=k1:N1J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);%消去运算endJ(k3,k)=0;endelsefor k3=k1:N0for k2=k1:N1J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);%消去运算endJ(k3,k)=0;endendend%====上面是用线性变换方式将Jacobi矩阵化成单位矩阵=====for k=3:2:N0-1L=(k+1)./2;e(L)=e(L)-J(k,N); %修改节点电压实部k1=k+1;f(L)=f(L)-J(k1,N); %修改节点电压虚部end%------修改节点电压-----------for k=3:N0DET=abs(J(k,N));if DET>=pr %电压偏差量是否满足要求IT2=IT2+1; %不满足要求的节点数加1endendICT2(a)=IT2;ICT1=ICT1+1;end%用高斯消去法解"w=-J*V"disp('迭代次数：');disp(ICT1);disp('没有达到精度要求的个数：');disp(ICT2);for k=1:nV(k)=sqrt(e(k)^2+f(k)^2);sida(k)=atan(f(k)./e(k))*180./pi;E(k)=e(k)+f(k)*j;end%=============== 计算各输出量=========================== disp('各节点的实际电压标幺值E为(节点号从小到大排列)：');disp(E);EE=E*UB;disp(EE);disp('-----------------------------------------------------');disp('各节点的电压大小V为(节点号从小到大排列)：');disp(V);VV=V*UB;disp(VV);disp('-----------------------------------------------------');disp('各节点的电压相角sida为(节点号从小到大排列)：');disp(sida);for p=1:nC(p)=0;for q=1:nC(p)=C(p)+conj(Y(p,q))*conj(E(q));endS(p)=E(p)*C(p);enddisp('各节点的功率S为(节点号从小到大排列)：');disp(S);disp('~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~');SS=S*SB;disp(SS);disp('-----------------------------------------------------');disp('各条支路的首端功率Si为(顺序同您输入B1时一致)：');for i=1:nlp=B1(i,1);q=B1(i,2);if B1(i,6)==0Si(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)*B1(i,5))-conj(E(q)))*conj(1./( B1(i,3)*B1(i,5))));Siz(i)=Si(p,q);elseSi(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)./B1(i,5))-conj(E(q)))*conj(1./( B1(i,3)*B1(i,5))));Siz(i)=Si(p,q);enddisp(Si(p,q));SSi(p,q)=Si(p,q)*SB;ZF=['S(',num2str(p),',',num2str(q),')=',num2str(SSi(p,q))];disp(ZF);%disp(SSi(p,q));disp('-----------------------------------------------------');enddisp('各条支路的末端功率Sj为(顺序同您输入B1时一致)：');for i=1:nlp=B1(i,1);q=B1(i,2);if B1(i,6)==0Sj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)./B1(i,5))-conj(E(p)))*conj(1./( B1(i,3)*B1(i,5))));Sjy(i)=Sj(q,p);elseSj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)*B1(i,5))-conj(E(p)))*conj(1./( B1(i,3)*B1(i,5))));Sjy(i)=Sj(q,p);enddisp(Sj(q,p));SSj(q,p)=Sj(q,p)*SB;ZF=['S(',num2str(q),',',num2str(p),')=',num2str(SSj(q,p))];disp(ZF);%disp(SSj(q,p));disp('-----------------------------------------------------');enddisp('各条支路的功率损耗DS为(顺序同您输入B1时一致)：');for i=1:nlp=B1(i,1);q=B1(i,2);DS(i)=Si(p,q)+Sj(q,p);disp(DS(i));DDS(i)=DS(i)*SB;ZF=['DS(',num2str(p),',',num2str(q),')=',num2str(DDS(i))];disp(ZF);%disp(DDS(i));disp('-----------------------------------------------------');endfigure(1);subplot(2,2,1);plot(V);xlabel('节点号');ylabel('电压标幺值');grid on;subplot(2,2,2);plot(sida);xlabel('节点号');ylabel('电压角度');grid on;subplot(2,2,3);bar(real(S));ylabel('节点注入有功');grid on;subplot(2,2,4);bar(Siz);ylabel('支路首端无功');grid on;1.冬季最大运行方式潮流计算结果:计算机运行的B1,B2阵如下：B1=[ 1 2 0.0318+0.0454*i 0.282*i 1 01 11 0.0114+0.0374*i 0.2332*i 1 011 3 0.001975+0.0695*i 0 1.025:1.1 11 12 0.0087+0.029*i 0.1788*i 1 012 4 0.0031+0.103*i 0 1:1.05 11 5 0.043+0.142*i 0.22*i 1 05 6 0.0031+0.103*i 0 1:1.05 15 7 0.043+0.142*i 0.22*i 1 01 7 0.051+0.168*i 0.26*i 1 07 8 0.00198+0.0695*i 0 1:1.1 17 13 0.01+0.033*i 0.204*i 1 013 9 0.0025+0.083*i 0 0.9956 11 10 0.00239+0.084*i 0 1.048 1]B2=[0 0 1.1 0 0 10 1.43+0.886*i 1.05 1.05 0 30 0.88+0.545*i 1 0 0 20 0.77+0.4772*i 1 0 0 20 0 1 0 0 20 0.77+0.4772*i 1 0 0 20 0 1 0 0 20 0.88+0.545*i 1 0 0 20.642+0.3817*i 0 1.05 1.05 0 31+0.75*i 0.2+0.1549*i 1.05 1.05 0 30 0 1 0 0 20 0 1 0 0 20 0 1 0 0 2]计算机运行结果如下表：节点号电压值相角值支路标号首端功率末端功率支路功率损耗1 242.0000 0 1-2 148.9391+3.327637i -143-27.45476i 5.93913-24.1271i2 231.0000 -3.0341 1-11 89.1956+37.193i -88.19803-61.4687i 0.997519-24.2757i3 227.4154 -4.4383 11-3 88.19803+61.4687i -88-54.5i 0.19803+6.9687i4 226.3304 -4.9004 1-12 77.8364+36.4186i -77.2404-55.7062i 0.596039-19.2876i5 229.8318 -3.9547 12-4 77.2404+55.7062i -77-47.72i 0.24036+7.9862i6 217.7226 -8.2720 1-5 63.5438+12.2204i -61.8773-32.0321i 1.66656-19.8116i7 234.9326 -3.1047 5-6 77.2597+56.3501i -77-47.72i 0.25974+8.6301i8 226.0991 -6.2429 5-7 15.3825-24.3181i 15.5354+0.274108i 0.152961-24.044i9 231.0000 0.5851 7-8 88.20085+61.55009i -88-54.5i 0.20085+7.0501i10 231.0000 3.4950 1-7 40.806-6.05737i -40.0647-22.0555i 0.741264-28.1128i11 236.1931 -1.3348 7-13 -63.6716-39.7687i 64.0965+17.5832i 0.424934-22.1855i12 237.9346 -0.8892 13-9 -64.0965-17.5832i 64.2+21.0187i 0.10348+3.4355i13 238.1868 -2.2028 1-10 79.8613+4.23546i 80+0.641048i 0.13875+4.8765i 计算机计算结果图形：2.冬季最小运行方式潮流计算结果：计算机运行的B1B2矩阵如下：B1=[ 1 2 0.0318+0.0454*i 0.282*i 1 01 11 0.0114+0.0374*i 0.2332*i 1 011 3 0.001975+0.0695*i 0 1.025:1.1 11 12 0.0087+0.029*i 0.1788*i 1 012 4 0.0031+0.103*i 0 1:1.05 11 5 0.043+0.142*i 0.22*i 1 05 6 0.0031+0.103*i 0 1:1.05 15 7 0.043+0.142*i 0.22*i 1 01 7 0.051+0.168*i 0.26*i 1 07 8 0.00198+0.0695*i 0 1:1.1 17 13 0.01+0.033*i 0.204*i 1 013 9 0.0025+0.083*i 0 0.9956 11 10 0.00239+0.084*i 0 1.048 1]B2=[0 0 1.1 0 0 10 1.43+0.886*i 1.05 1.05 0 30 0.616+0.3817*i 1 0 0 20 0.539+0.3817*i 1 0 0 20 0 1 0 0 20 0.539+0.334*i 1 0 0 20 0 1 0 0 20 0.539+0.334*i 1 0 0 20.642+0.3817*i 0 1.05 1.05 0 31+0.75*i 0.1+0.06197*i 1.05 1.05 0 30 0 1 0 0 20 0 1 0 0 20 0 1 0 0 2]电压调整措施：变电所1、4变压器变比：+2.5% 水电厂变压器变比：-2.5%5 234.2241 -2.0051 (12--4) 54.0234+42.2706i -53.9-38.17i 0.12342+4.1006i6 226.2159 -4.8577 (1--5) 34.1669+4.22856 -33.6427-28.276i 0.524165-24.0474i7 235.1128 -0.4737 (5--6) 54.0179+37.3169i -53.9-33.4i 0.11789+3.9169i8 223.9941 -2.4603 (5--7) -20.3752-9.04094i 20.5371-15.4558i 0.161947-24.4968i9 231.0000 3.4429 (1--7) 11.0013+1.45186i -10.8258-31.4513i 0.175442-29.9995i10 231.0000 3.9321 (7--8) 53.9768+36.0957i -53.9-33.4i 0.076797+2.6957i11 238.4187 -0.9561 (7--13) -63.6881+10.8115i 64.0874-32.7763i 0.39932-21.9648i12 239.1742 -0.6185 (13--9) -64.0874+32.7763i 64.2-29.0386i 0.11258+3.7377i13 234.9468 0.6942 (1--10) -89.8244+5.1673i 90+1.0049i 0.17561+6.1722i 计算机运行结果的图形：3.夏季最大运行方式计算机计算结果：计算机运行B1B2阵如下：B1=[ 1 2 0.0318+0.0454*i 0.282*i 1 01 11 0.0114+0.0374*i 0.2332*i 1 011 3 0.001975+0.0695*i 0 1.025:1.1 11 12 0.0087+0.029*i 0.1788*i 1 012 4 0.0031+0.103*i 0 1:1.05 11 5 0.043+0.142*i 0.22*i 1 05 6 0.0031+0.103*i 0 1:1.05 15 7 0.043+0.142*i 0.22*i 1 01 7 0.051+0.168*i 0.26*i 1 07 8 0.00198+0.0695*i 0 1:1.1 17 13 0.02+0.066*i 0.102*i 1 013 9 0.0025+0.083*i 0 0.9956 11 10 0.00239+0.084*i 0 1.048 1]B2=[0 0 1.1 0 0 10 1.26+0.7814*i 1.05 1.05 0 30 0.776+0.481*i 1 0 0 20 0.543+0.3367*i 1 0 0 20 0 1 0 0 20 0.543+0.3367*i 1 0 0 20 0 1 0 0 20 0.776+0.481*i 1 0 0 21.35+0.6538*i 0.2+0.124*i 1.05 1.05 0 31+0.75*i 0.25+0.1549*i 1.05 1.05 0 30 0 1 0 0 20 0 1 0 0 20 0 1 0 0 2]计算机运行结果如下表：10 231.0000 3.4950 (7,8) 77.7585+53.6634i -77.6-48.1i0.1585+5.5634i11 237.0938 -1.1858 (7,13) -113.4984+9.476406i 114.6926-28.38885i 1.19422-18.9124i12 239.1742 -0.6185 (13,9) -114.6926+28.38885 115-18.18369i0.307384+10.2052i13 233.3912 3.2303 (1,10) -79.8613+4.23546i 80+0.641048i 0.13875+4.8765i 计算机计算结果如图：4.夏季最小运行方式：计算机运行B1B2阵如下：B1=[ 1 2 0.0318+0.0454*i 0.282*i 1 01 11 0.0114+0.0374*i 0.2332*i 1 011 3 0.001975+0.0695*i 0 1.025:1.1 11 12 0.0087+0.029*i 0.1788*i 1 012 4 0.0031+0.103*i 0 1:1.05 11 5 0.043+0.142*i 0.22*i 1 05 6 0.0031+0.103*i 0 1:1.05 15 7 0.043+0.142*i 0.22*i 1 01 7 0.051+0.168*i 0.26*i 1 07 8 0.00198+0.0695*i 0 1:1.1 17 13 0.02+0.066*i 0.102*i 1 013 9 0.0025+0.083*i 0 0.9956 11 10 0.00239+0.084*i 0 1.048 1]B2=[0 0 1.1 0 0 10 1.26+0.7814*i 1.05 1.05 0 30 0.543+0.336*i 1 0 0 20 0.4753+0.2947*i 1 0 0 20 0 1 0 0 20 0.4753+0.2947*i 1 0 0 20 0 1 0 0 20 0.543+0.336*i 1 0 0 21.35+0.6538*i 0.2+0.124*i 1.05 1.05 0 31+0.75*i 0.25+0.1549*i 1.05 1.05 0 30 0 1 0 0 20 0 1 0 0 20 0 1 0 0 2]10 231.0000 3.9321 (7,8) 54.3735+36.2509i -54.3-33.67i 0.073528+2.5809i11 239.0099 -0.8518 (7,13) -113.4428+24.39707i 114.6754-43.79301i 1.23255-19.3959i12 239.8726 -0.5689 (13,9) -114.6754+43.79301i 115-33.01613i 0.324605+10.7769i13 235.9565 4.4510 (1,10) -89.8244+5.1673i 90+1.0049i 0.17561+6.1722i 计算机计算结果如图：5.夏季故障运行状态:调压及无功补偿措施如下:变电所3的变压器变比为-2.5%,无功补偿容量为20Mvar。

电力系统潮流计算完整程序及详细理论说明

电力系统潮流计算完整程序及详细理论说明——秦羽风在我刚开始学习潮流程序时,总是找不到一个正确的程序开始模仿学习。

后来经过多方努力，终于自己写出了一个结构清晰、完整的潮流程序。

此程序是一个通用的程序，只需要修改输入数据的子函数（PowerFlowsData_K)里面的母线、支路、发电机、负荷，就能算任意一个网络结构的交流系统潮流。

很适合初学者学习.为了帮助电力系统的同学一起学习,我将我编写的潮流计算程序分享下来给大家；此程序是在基于牛顿拉夫逊算法的基础上，编写的快速解耦算法。

每一个子程序我都有备注说明。

如果有不对的地方，希望大家指正！下文中呈现的顺序为：网络结构、子程序、主程序、运算结果、程序设计理论说明。

一、网络结构：5节点网络如下图。

二、子程序（共有9个子程序）子程序1：（其他系统，只需要修改Bus、Branch、Generator、Load,这四个矩阵就行了）function [Bus,Branch，Generator,Load］=PowerFlowsData_K％％节点数据% 类型:1-平衡节点；2-发电机PV节点；3—负荷PQ节点；4-发电机PQ节点;Bus=［% 类型电压相角1 1。

06 0；2 1 0；3 1 0;3 1 0；3 1 0]；％% 线路数据Branch=[% 发送接收电阻电感(电导电容)并联1 2 0.02 0.06 0 0.06；1 3 0。

08 0。

24 0 0。

05；2 3 0.06 0.18 0 0。

04;2 4 0。

06 0。

18 0 0.04;2 5 0.04 0.12 0 0。

03；3 4 0.01 0.03 0 0。

02;4 5 0.08 0.24 0 0.05］;%% 发电机数据Generator=［％节点定有功定无功(上限下限）无功1 0 0 5 —5;2 0。

4 03 —3];％％负载数据Load=［% 节点定有功定无功2 0.2 0.1；3 0。

Matlab实现潮流计算程序

程序代码如下:111111.%读入数据clcclearfilename='123.txt';a=textread(filename)n=a(1,1);pinghengjd=a(1,2);phjddianya=a(1,3);jingdu=a(1,4);b=zeros(1,9);j1=0;[m1,n1]=size(a);for i1=1:m1if a(i1,1)==0j1=j1+1;b(j1)=i1;endendb;%矩阵分块a1=a(b(1)+1:b(2)-b(1)+1,1:n1);a2=a(b(2)+1:b(3)-1,1:n1);a3=a(b(3)+1:b(4)-1,1:n1);a4=a(b(4)+1:b(5)-1,1:n1);a5=a(b(5)+1:b(6)-1,1:n1);%设置初值vcz=1;dcz=0;kmax=20;k1=0;%求节点导纳矩阵a11=zeros(4,6);for i0=1:3for j0=1:6a11(i0,j0)=a1(i0,j0);a11(4,j0)=a2(1,j0);endenda11;linei=a11(1:4,2);linej=a11(1:4,3);liner=a11(1:4,4);linex=a11(1:4,5);lineb=a11(1:4,6);branchi=0;branchj=0;branchb=0;G=zeros(4,4);B=zeros(4,4);for k=1:4i2=linei(k,1);j2=linej(k,1);r=liner(k,1);x=linex(k,1);b=0;GIJ=r/(r*r+x*x);BIJ=-x/(r*r+x*x);if k>=4 & lineb(k)~=0k0=lineb(k);G(i2,j2)=-GIJ/k0;G(j2,i2)=G(i2,j2);B(i2,j2)=-BIJ/k0;B(j2,i2)=B(i2,j2);G(i2,i2)=G(i2,i2)+GIJ/k0/k0; B(i2,i2)=B(i2,i2)+BIJ/k0/k0;elseG(j2,i2)=-GIJ;G(i2,j2)=G(j2,i2);B(j2,i2)=-BIJ;B(i2,j2)=B(j2,i2);G(i2,i2)=G(i2,i2)+GIJ;b=lineb(k);B(i2,i2)=B(i2,i2)+BIJ+b;endG(j2,j2)=G(j2,j2)+GIJ;B(j2,j2)=B(j2,j2)+BIJ+b;endG;B;B=B.*i;Yf=G+BY=abs(Yf);alf=angle(Yf);%赋Jacobian矩阵参数P=zeros(n,1);Q=zeros(n,1);Pd=zeros(1,n);Qd=zeros(1,n);dP=zeros(1,n);dQ=zeros(1,n);PG=a4(:,3);PD=a4(:,5);QG=a4(:,4);QD=a4(:,6);i8=a4(:,2);for j8=1:length(i8)P(i8(j8))=PG(i8(j8))-PD(i8(j8));Q(i8(j8))=QG(i8(j8))-QD(i8(j8));enddelt=zeros(n,1);V=ones(n,1);V(3)=1.10;V(4)=1.05;ddelt=zeros(n,1);dV=zeros(n,1);A=zeros(2*n,2*n);B=zeros(2*n,1);Jacobian=Jaco(V,delt,n,Y,alf)%求取矩阵功率for j5=1:kmaxdisp(['第' int2str(j5) '次计算结果'])if k>=kmaxbreakendfor i10=1:4Pd(i10)=0;Qd(i10)=0;for j10=1:nPd(i10)=Pd(i10)+V(i10)*Y(i10,j10)*V(j10)*cos(d elt(i10)-delt(j10)-alf(i10,j10));Qd(i10)=Qd(i10)+V(i10)*Y(i10,j10)*V(j10)*sin(d elt(i10)-delt(j10)-alf(i10,j10));endendfor i4=1:3dP(i4)=P(i4)-Pd(i4);endfor j4=1:2dQ(j4)=Q(j4)-Qd(j4);endA=Jaco(V,delt,n,Y,alf)for i14=1:nB(i14*2-1)=-dP(i14);B(i14*2)=-dQ(i14);endif max(abs(B))>jingduX=A\B;for i16=1:nddelt(i16)=X(2*i16-1);dV(i16)=X(2*i16)*V(i16);endV=V+dVdelt=delt+ddeltelsebreakenddisp('----------------')end%流氓算法% for ii=1:2% V(ii)=V(ii)+dV(ii);% end% V222222.function A=Jaco(V,delt,n,Y,alf)%计算Jacobian矩阵for i7=1:nHd1(i7)=0;Jd1(i7)=0;for j7=1:nHd1(i7)=Hd1(i7)+V(i7)*Y(i7,j7)*V(j7)*sin(delt(i7)-delt(j7)-alf(i7,j7));Jd1(i7)=Jd1(i7)+V(i7)*Y(i7,j7)*V(j7)*cos(delt(i7)-delt(j7)-alf(i7,j7));endendfor i6=1:nfor j6=1:nif i6~=j6H(i6,j6)=-V(i6)*Y(i6,j6)*V(j6)*sin(delt(i6)-delt(j6)-alf(i6,j6));N(i6,j6)=-V(i6)*Y(i6,j6)*V(j6)*cos(delt(i6)-delt(j6)-alf(i6,j6));J(i6,j6)=-N(i6,j6);L(i6,j6)=H(i6,j6);elseH(i6,i6)=Hd1(i6)-V(i6)*Y(i6,i6)*V(i6)*sin(delt(i6)-delt(j6)-alf(i6,j6));J(i6,j6)=-Jd1(i6)+V(i6)*Y(i6,j6)*V(j6)*cos(delt(i6)-delt(j6)-alf(i6,j6));N(i6,j6)=-Jd1(i6)-V(i6)*Y(i6,i6)*V(i6)*cos(alf(i6,i6));L(i6,i6)=-Hd1(i6)+V(i6)*Y(i6,i6)*V(i6)*sin(alf(i6,i6));endendend%修正Jacobian矩阵for j9=3for i9=1:nN(i9,j9)=0;L(i9,j9)=0;J(j9,i9)=0;L(j9,i9)=0;endendL(j9,j9)=1;for j9=4for i9=1:nH(i9,j9)=0;N(i9,j9)=0;J(i9,j9)=0;L(i9,j9)=0;H(j9,i9)=0;N(j9,i9)=0;J(j9,i9)=0;L(j9,i9)=0;endendH(j9,j9)=1;L(j9,j9)=1;%Jaco=[H N;J L];%Jaco=zeros(2*n,2*n);for i11=1:nfor j11=1:nJaco(2*i11-1,2*j11-1)=H(i11,j11); Jaco(2*i11-1,2*j11)=N(i11,j11); Jaco(2*i11,2*j11-1)=J(i11,j11);Jaco(2*i11,2*j11)=L(i11,j11);endendA=Jaco;33333.数据:4 4 1.05 0.000011 12 0.1 0.40 0.015282 1 4 0.12 0.50 0.019203 24 0.08 0.40 0.014131 1 3 0 0.3 0.909090911 1 0 0 0.30 0.182 2 0 0 0.55 0.133 3 0.5 0 0 01 3 1.10 0 0。

潮流计算程序（3.0版）说明

潮流计算程序(3.0版)说明1特点∙既可以用于高压输电网的潮流解算，又可以用于低压配电网的潮流解算，还可以同时解算输电网加配电网的混合潮流问题。

∙可以处理多平衡（机）节点问题，使用者只需输入各平衡节点的电压幅值和相位角，计算出的系统不平衡功率部分将自动在各平衡节点间进行分摊。

∙可以同时解算相互解列的几个子系统的潮流问题（只要每个子系统内均含有平衡节点），甚至某些子系统退化成孤立节点也不会影响其它子系统的潮流解算。

∙可以同时解算多配电馈线的潮流问题，而无需一个馈线一个馈线地分别计算。

∙由于该程序能够处理多平衡节点问题，故在解算具有多根节点和环状配电网潮流时，无需象以往方法那样只保留一个平衡节点，而将其余平衡节点全部人为地改成PV节点。

∙该版程序为PQ分解法。

2一般变量说明∙LINEMAX 程序所能处理的最大线路数，可以在#define说明语句中进行修改∙GENERA TORMAX 程序所能处理的最大发电机节点数，可以在#define说明语句中进行修改∙LOADMAX 程序所能处理的最大负荷节点数，可以在#define说明语句中进行修改∙NODEMAX 程序所能处理的最大系统节点总数，可以在#define说明语句中进行修改∙SWINGMAX 程序所能处理的最大平衡节点数，可以在#define说明语句中进行修改∙PVMAX 程序所能处理的最大PV节点数，可以在#define说明语句中进行修改∙NODEFACTOR 导纳矩阵的上三角阵中程序所能处理的最大非零非对角元素的个数相对于最大节点数（NODEMAX）的倍数∙Deg_to_Rad 度到弧度的转换系数，在#define说明语句中定义∙Rad_to_Deg 弧度到度的转换系数，在#define说明语句中定义∙Num_Line 系统的实际总线路数∙Num_Gen 系统的实际总发电机（节点）数∙Num_Load 系统的实际总负荷（节点）数∙Num_Node 系统的实际总节点数∙Num_Swing 系统的实际总平衡节点数∙Num_GPV 系统发电机节点中的PV节点数∙Num_GPQ 系统发电机节点中的PQ节点数∙Eps 节点功率失配量值收敛限值∙Iter_Max 迭代次数限值（在不满足电压收敛误差限值的条件下强行中止收敛的最大迭代次数）∙V olIni_Flag 是否读取电压初值标志：0-不读（以0-1启动）；1-读初值∙V olRes_Flag 是否将电压收敛值保存以备以后计算时当启动初值：0-不保存；1-保存3参数结构定义说明∙线路参数结构定义struct Line{int Node_No[2];int Flag;double RXBK[3];}LLine[LINEMAX]其中，Node_No[2]代表线路两端节点号：0—左节点号；1—右节点号。

潮流计算步骤

潮流计算步骤
潮流计算是电力系统分析中的一种基本计算方法，用于确定电网中的电压分布和功率流动情况。

以下是潮流计算的基本步骤：
1、输入原始数据和信息：包括电网的结构信息、设备参数、负荷和电源的分布及大小等。

2、建立数学模型：根据电路理论和电力系统网络模型，建立描述电力系统中电压、电流和功率关系的数学模型。

3、形成节点导纳矩阵：根据电网结构，形成节点导纳矩阵，用于描述系统中各节点之间的电气联系。

4、确定待求状态变量初值：根据实际情况，为待求的状态变量（如节点电压）设定初值。

5、迭代求解：使用迭代法对数学模型进行求解，逐步更新状态变量的值，直到满足收敛条件为止。

6、计算节点电压：根据迭代求解的结果，计算出各节点的电压值。

7、计算功率分布：根据节点电压和网络参数，计算出各支路的功率流动情况。

8、结果分析：对计算结果进行整理和分析，评估电网的运行状态，为进一步优化和调整提供依据。

需要注意的是，潮流计算的具体步骤可能会因不同的计算方法和电力系统分析软件而有所差异。

在实际应用中，需要根据具体的软件
和要求进行操作。

电力系统潮流计算程序(详细)

潮流计算，电力学名词，指在给定电力系统网络拓扑、元件参数和发电、负荷参量条件下，计算有功功率、无功功率及电压在电力网中的分布。

潮流计算是根据给定的电网结构、参数和发电机、负荷等元件的运行条件，确定电力系统各部分稳态运行状态参数的计算。

通常给定的运行条件有系统中各电源和负荷点的功率、枢纽点电压、平衡点的电压和相位角。

待求的运行状态参量包括电网各母线节点的电压幅值和相角，以及各支路的功率分布、网络的功率损耗等。

潮流计算程序，相关的原始数据数据数据输入格式如下：%B1是支路参数矩阵，第一列和第二列是起始节点编号和终点节点编号%第三列、第四列、第五列、第六列、第七列、第八列分别为：支路电阻、电抗、电导、电纳、变压器变比、是否有变压器（1为有、0为无）。

%B2为节点参数矩阵，其中第一列到第六列为节点编号；为节点类型；注入有功、注入无功、电压幅值、电压相位。

%“1”为平衡节点，“2”为PQ节点，“3”为PV节点参数。

n=input('请输入节点数:n=');n1=input('请输入支路数:n1=');isb=input('请输入平衡节点号:isb=');pr=input('请输入误差精度:pr=');B1=input('请输入支路参数:B1=');B2=input('请输入节点参数:B2=');Y=zeros(n);N=1;%建立节点导纳矩阵for i=1:n1if B1(i,8)==0p=B1(i,1);q=B1(i,2);Y(p,q)=Y(p,q)-1/(B1(i,3)+B1(i,4))-B1(i,5);Y(q,p)=Y(p,q);Y(p,p)=Y(p,p)+1/(B1(i,3)+B1(i,4))+0.5*B1(i,6)+B1(i,5);Y(q,q)=Y(q,q)+1/(B1(i,3)+B1(i,4))+0.5*B1(i,6)+B1(i,5);elsep=B1(i,1);q=B1(i,2);Y(p,q)=Y(p,q)-1/((B1(i,3)+B1(i,4))*B1(i,7)) -B1(i,5);Y(q,p)=Y(p,q);Y(p,p)=Y(p,p)+1/(B1(i,3)+B1(i,4))+B1(i,5);Y(q,q)=Y(q,q)+1/(B1(i,7)^2*(B1(i,3)+B1(i,4)))+B1(i,5);endendYG=real(Y);B=imag(Y);PriS=zeros(2*n-2,1);ImbS=zeros(2*n-2,1);%创建PriS，用于存储初始功率参数h=0;j=0;for i=1:nif i~=isb&B2(i,2)==1h=h+1;for j=1:nPriS(2*h-1,1)=PriS(2*h-1,1)+B2(i,5)*(G(i,j)*B2(j,5)-B(i,j)*B2(j,6))+B2(i,6)*(G(i,j)*B2(j,6)+B(i, j)*B2(j,5));PriS(2*h,1)=PriS(2*h,1)+B2(i,6)*(G(i,j)*B2(j,5)-B(i,j)*B2(j,6))-B2(i,5)*(G(i,j)*B2(j,6)+B(i,j)*B 2(j,5));endendendfor i=1:nif i~=isb&B2(i,2)==2h=h+1;for j=1:nPriS(2*h-1,1)=PriS(2*h-1,1)+B2(i,5)*(G(i,j)*B2(j,5)-B(i,j)*B2(j,6))+B2(i,6)*(G(i,j)*B2(j,6)+B(i,j)*B2(j,5));PriS(2*h,1)=PriS(2*h,1)+B2(i,6)*(G(i,j)*B2(j,5)-B(i,j)*B2(j,6))-B2(i,5)*(G(i,j)*B2(j,6)+B(i,j)*B2(j,5))endendendPriSU3=zeros(n-h-1,2);%U3存储PV节点的初始电压t=0;for i=1:nif B2(i,2)==2t=t+1;U3(t,1)=B2(i,5);U3(t,2)=B2(i,6);endendU3%ImbS于存储有功功率、无功功率和电压幅值的不平衡量h=0;for i=1:nif i~=isb&B2(i,2)==1h=h+1;ImbS(2*h-1,1)=B2(i,3)-PriS(2*h-1,1);ImbS(2*h,1)=B2(i,4)-PriS(2*h,1);endendt=0;for i=1:nif i~=isb&B2(i,2)==2h=h+1;t=t+1;ImbS(2*h-1,1)=B2(i,3)-PriS(2*h-1,1);ImbS(2*h,1)=U3(t,1)^2+U3(t,2)^2-B2(i,5)^2-B2(i,6)^2;endendImbSI=zeros(n-1,1);%I，存储节点电流参数h=0;for i=1:nif i~=isbh=h+1;I(h,1)=(PriS(2*h-1,1)-PriS(2*h,1)*sqrt(-1))/conj(B2(i,5)+B2(i,6)*sqrt(-1));endendIJacbi=zeros(2*n-2);%Jacbi(雅可比矩阵)h=0;k=0;for i=1:nif B2(i,2)==1h=h+1;for j=1:nif j~=isbk=k+1;if i==jJacbi(2*h-1,2*k-1)=-B(i,j)*B2(i,5)+G(i,j)*B2(i,6)+imag(I(h,1));Jacbi(2*h-1,2*k)=G(i,j)*B2(i,5)+B(i,j)*B2(i,6)+real(I(h,1));Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k)+2*real(I(h,1));Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1)-2*imag(I(h,1));elseJacbi(2*h-1,2*k-1)=-B(i,j)*B2(i,5)+G(i,j)*B2(i,6);Jacbi(2*h-1,2*k)=G(i,j)*B2(i,5)+B(i,j)*B2(i,6);Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k);Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1);endif k==(n-1)k=0;endendendendendk=0;for i=1:nif B2(i,2)==2h=h+1;for j=1:nif j~=isbk=k+1;if i==jJacbi(2*h-1,2*k-1)= -B(i,j)*B2(i,5)+G(i,j)*B2(i,6)+imag(I(h,1));Jacbi(2*h-1,2*k)= G(i,j)*B2(i,5)+B(i,j)*B2(i,6)+real(I(h,1));Jacbi(2*h,2*k-1)=2*B2(i,6);Jacbi(2*h,2*k)=2*B2(i,5);elseJacbi(2*h-1,2*k-1)= -B(i,j)*B2(i,5)+G(i,j)*B2(i,6);Jacbi(2*h-1,2*k)= G(i,j)*B2(i,5)+B(i,j)*B2(i,6);Jacbi(2*h,2*k-1)=0;Jacbi(2*h,2*k)=0;endif k==(n-1)k=0;endendendendendJacbi%求解修正方程，获取节点电压的不平衡量ImbU=zeros(2*n-2,1);ImbU=inv(Jacbi)*ImbS;ImbU%修正节点电压j=0;for i=1:nif B2(i,2)==1j=j+1;B2(i,5)=B2(i,5)+ImbU(2*j,1);B2(i,6)=B2(i,6)+ImbU(2*j-1,1);endendfor i=1:nif B2(i,2)==2j=j+1;B2(i,5)=B2(i,5)+ImbU(2*j,1);B2(i,6)=B2(i,6)+ImbU(2*j-1,1);endendB2while abs(max(ImbU))>prPriS=zeros(2*n-2,1);h=0;j=0;for i=1:nif i~=isb&B2(i,2)==1h=h+1;for j=1:nPriS(2*h-1,1)=PriS(2*h-1,1)+B2(i,5)*(G(i,j)*B2(j,5)-B(i,j)*B2(j,6))+B2(i,6)*(G(i,j)*B2(j,6)+B(i, j)*B2(j,5));PriS(2*h,1)=PriS(2*h,1)+B2(i,6)*(G(i,j)*B2(j,5)-B(i,j)*B2(j,6))-B2(i,5)*(G(i,j)*B2(j,6)+B(i,j)*B 2(j,5));endendendfor i=1:nif i~=isb&B2(i,2)==2h=h+1;for j=1:nPriS(2*h-1,1)=PriS(2*h-1,1)+B2(i,5)*(G(i,j)*B2(j,5)-B(i,j)*B2(j,6))+B2(i,6)*(G(i,j)*B2(j,6)+B(i,j)*B2(j,5));PriS(2*h,1)=PriS(2*h,1)+B2(i,6)*(G(i,j)*B2(j,5)-B(i,j)*B2(j,6))-B2(i,5)*(G(i,j)*B2(j,6)+B(i,j)*B2(j,5))endendendPriS%创建ImbSh=0;for i=1:n %对PQ节点的处理if i~=isb&B2(i,2)==1h=h+1;ImbS(2*h-1,1)=B2(i,3)-PriS(2*h-1,1);ImbS(2*h,1)=B2(i,4)-PriS(2*h,1);endendt=0;for i=1:n %对PV节点的处理if i~=isb&B2(i,2)==2h=h+1;t=t+1;ImbS(2*h-1,1)=B2(i,3)-PriS(2*h-1,1);ImbS(2*h,1)=U3(t,1)^2+U3(t,2)^2-B2(i,5)^2-B2(i,6)^2;endendImbS%创建II=zeros(n-1,1);h=0;for i=1:nif i~=isbh=h+1;I(h,1)=(PriS(2*h-1,1)-PriS(2*h,1)*sqrt(-1))/conj(B2(i,5)+B2(i,6)*sqrt(-1));endendI%创建JacbiJacbi=zeros(2*n-2);h=0;k=0;for i=1:nif B2(i,2)==1h=h+1;for j=1:nif j~=isbk=k+1;if i==jJacbi(2*h-1,2*k-1)=-B(i,j)*B2(i,5)+G(i,j)*B2(i,6)+imag(I(h,1));Jacbi(2*h-1,2*k)=G(i,j)*B2(i,5)+B(i,j)*B2(i,6)+real(I(h,1));Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k)+2*real(I(h,1));Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1)-2*imag(I(h,1));elseJacbi(2*h-1,2*k-1)=-B(i,j)*B2(i,5)+G(i,j)*B2(i,6);Jacbi(2*h-1,2*k)=G(i,j)*B2(i,5)+B(i,j)*B2(i,6);Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k);Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1);endif k==(n-1)k=0;endendendendendk=0;for i=1:nif B2(i,2)==2h=h+1;for j=1:nif j~=isbk=k+1;if i==jJacbi(2*h-1,2*k-1)= -B(i,j)*B2(i,5)+G(i,j)*B2(i,6)+imag(I(h,1));Jacbi(2*h-1,2*k)= G(i,j)*B2(i,5)+B(i,j)*B2(i,6)+real(I(h,1));Jacbi(2*h,2*k-1)=2*B2(i,6);Jacbi(2*h,2*k)=2*B2(i,5);elseJacbi(2*h-1,2*k-1)= -B(i,j)*B2(i,5)+G(i,j)*B2(i,6);Jacbi(2*h-1,2*k)= G(i,j)*B2(i,5)+B(i,j)*B2(i,6);Jacbi(2*h,2*k-1)=0;Jacbi(2*h,2*k)=0;endif k==(n-1)k=0;endendendendendJacbiImbU=zeros(2*n-2,1);ImbU=inv(Jacbi)*ImbS;ImbU%修正节点电压j=0;for i=1:nif B2(i,2)==1j=j+1;B2(i,5)=B2(i,5)+ImbU(2*j,1);B2(i,6)=B2(i,6)+ImbU(2*j-1,1);endendfor i=1:nif B2(i,2)==2j=j+1;B2(i,5)=B2(i,5)+ImbU(2*j,1);B2(i,6)=B2(i,6)+ImbU(2*j-1,1);endendB2N=N+1; %迭代次数加1endN。

牛拉法潮流计算程序（附3机9节点结果对比）

摘要电力系统潮流计算是研究电力系统稳态运行的一种重要方法,它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态,包括各母线的电压、线路的功率分布以及功率损耗等等。

潮流计算主要用于电网规划和静态安全分析，它可为扩建电力网络，以达到规划周期内所需要的输电能力提供依据；也可以对预想事故进行模拟和分析，校核预想事故下的电力系统安全性。

本文简单介绍了牛顿-拉夫逊潮流计算的原理、模型与算法，然后用具体的实例，利用MATLAB对牛顿-拉夫逊法的算法进行了验证。

关键词：电力系统潮流计算牛顿-拉夫逊法 MATLAB一、牛拉法的数学模型对一个N 节点的电力网路，列写节点电压方程，即I =Y V(1.1)式中，I 为节点注入电流列相量，Y 为节点导纳矩阵，V 为节点电压列相量。

由于异地测量的两个电流缺少时间同步信息，以注入功率替换注入电流作为已知量。

即***1+niij j ij j i i i Y V V I V Q P ∙∙===∑(1.2)其中,Y ij =G ij +j B ij ，带入上式，得到有功功率和无功功率方程 P i =V i V j G ij cos θij +B ij sin θij n j =1 (1.3)Q i =V i V j G ij sin θij −B ij cos θij n j =1(1.4)大部分情况下，已知PQ ，求解V θ。

考虑到电网的功率平衡，至少选择一台发电机来平衡全网有功功率，即至少有一个平衡节点，常选择调频或出线较多的发电机作为平衡节点。

具有无功补偿的母线能保持电压幅值恒定，这类节点可作为PV 节点。

潮流计算中节点分类总结如下：已知电力系统有m 个PQ 节点，r 个PV 节点和1个平衡节点，则可以提取m+r 个有功功率方程和m 个无功功率方程，从而求解出m+r 个θ和m 个V ，其余节点的有功和无功可通过式（1.3）、（1.4）求得，这样就完成了潮流计算。

潮流计算程序

%% included in power flow solution, not necessarily in input
PF = 14; %% real power injected at "from" bus end (MW) (not in PTI format)
% (assume OPF objective function has units, u)
% 20 MU_ANGMIN Kuhn-Tucker multiplier lower angle difference limit (u/degree)
% 21 MU_ANGMAX Kuhn-Tucker multiplier upper angle difference limit (u/degree)
QF = 15; %% reactive power injected at "from" bus end (MVAr) (not in PTI format)
PT = 16; %% real power injected at "to" bus end (MW) (not in PTI format)
% ANGMIN, ANGMAX, MU_ANGMIN, MU_ANGMAX] = idx_brch
%
% Some examples of usage, after defining the constants using the line above,
% are:
%
% branch(4, BR_STATUS) = 0; % take branch 4 out of service
%IDX_BRCH Defines constants for named column indices to branch matrix.

潮流计算程序

潮流程序报告班级：硕721姓名：胡与非学号：3107161033一、概述本程序是PQ分解法中的BX法，本程序中主要用到的三个主要算法概述如下：1．节点优化编号本程序使用Tinney-2编号方法，也称为最小度算法，或半动态节点优化编号方法。

该方法使用最为广泛，且简单有效，可以大大减少矩阵因子分解过程中产生的注入元的数目。

2．稀疏技术稀疏技术包括稀疏矩阵技术和稀疏矢量技术。

本程序使用了稀疏矩阵技术，利用网络矩阵的稀疏结构，减少了存储量和计算量，提高了求解网络方程的计算速度。

3．潮流计算的快速分解法（FDLF）快速分解法不是求解高维数的修正方程，而是分别交替求解两个低维数且系数矩阵是常数矩阵的修正方程，因此计算速度快，收敛性好，应用广泛。

二、输入数据格式说明第一行：支路数第二行：节点数第三行：接地支路数第四行：发电机数第五行：负荷数第六行：平衡节点号第七行：精度以下是支路数据：1、支路一端节点号支路另一端节点号支路电阻支路电抗接地电容(变压器变比)（节点号为负表明变压器支路,负号的节点侧为变压器的非标准变比侧,此时节点电容换为变压器变比；）2、接地支路号电阻电容以下是发电机数据：（最后一行是平衡机，其余发电机看成PV节点）节点号有功功率无功功率电压（电压为负表示PV节点）以下是负荷数据：节点号有功功率无功功率电压(初始为0)三、数据结构介绍//----------------线路数据结构体--------------//struct Line{int I ;//线路一端所连的节点号，节点号绝对值i<j,i按绝对值从小到大排列;//节点号为负表示为变压器支路，负号的节点侧为变压器的非标准变比侧;此时节点电容换为变压器变比int j; //线路另一端所连的节点号double r;//线路电阻double x;//线路电抗double b;//接地电容，若为变压器节点电容换为变压器变比//这里的b为线路导纳的一半b/2int hh; //标记支路号用于还原输出}*line;//----------------接地线路数据结构体--------------//struct Grand_line{int i; //接地支路号double r;//接地支路电阻double b;//接地支路电容int hh; //标记支路号，用于还原输出}*grand;//----------------发电机数据结构体--------------//struct Generator{int i; //发电机节点号,最后一个发电机为平衡节点double P;//发电机有功功率double Q;//发电机无功功率double V;//发电机机端电压,电压为负表示PV节点}*generator;struct Pinghengji{int i;//平衡机节点号,最后一个发电机为平衡节点double P;//平衡机有功功率double Q;//平衡机无功功率double V;//平衡机机端电压}*pinghengji;//----------------负荷数据结构体--------------//struct Load{int i;//负荷节点号double P;//负荷的有功功率double Q;//负荷的无功功率double V;//负荷处的电压(初始为0)}*load;//----------------优化结构体--------------//struct Y ouhua{int nn;//记录节点的度int nn1;int *jd;//记录节点号int *w;}*youhua;//-------------自导纳结构体--------------// struct Yii{double g,b;}*yii;struct Yii1{double g,b;}*yii1;//------------互导纳结构体----------------// struct Yij{double g,b;int j;}*yij;struct Yij1{double g,b;int j;}*yij1;//----------------U--------------//struct U_Type{double value;int j;}*U1,*U2,*U;//-------------------节点电压数组-----------------// struct NodalV oltage_Type{double V,theta;}*NodalV oltage;三、程序流程介绍四、程序中重要函数体的说明1、读数据函数read();读入各个原始数据；输入5，则为IEEE5节点；输入14，则为IEEE14节点；输入30，则为IEEE30；输入57，则为IEEE57；输入118，则为IEEE118；2、节点优化函数op();在此函数里实现节点优化以及支路、发电机、负荷的重排；3、形成导纳矩阵函数daona();形成导纳矩阵，在这个函数中同时调用函数formb(1)形成因子表1；调用formbb(2)形成因子表2；4、计算节点电压的迭代函数Diedai();迭代过程见流程图；5、打印函数dayin();打印各节点电压的辐值和相位，各个支路的功率和整个网络的网损；五、小结1．节点优化编号：节点优化编号直接影响到矩阵A的因子表矩阵的稀疏度，而稀疏度的大小又直接影响到稀疏技术的使用效率的好坏。

牛顿拉夫逊潮流计算程序

f1 x 2 xn f 2 f 2 x 2 xn f n f n x 2 xn
f i
f1
Δ X 1 Δ X 2 Δ X n
（13）
' ' 式（13）等号右边矩阵的 x 等都是对于 X1 ,X 2 的值，这一矩阵称为雅可比（Jacobi）矩阵。
P—Q分解法潮流计算的计算机算法。
电力系统分析
1概述
类型：
导纳法
阻抗法
牛顿-拉夫逊法（N—R法）快速分解法（ PQ分解法）
电力系统分析
2 潮流计算的基本方程
2.1节点的分类 2.2基本方程式
电力系统分析
2.1节点(bus)的分类
根据电力系统中各节点性质的不同，可把节点分成三种类型。 1.PQ节点事先给定的是节点功率（P、Q），待求的是节点电压向量（U、θ）。通常变电所母线都是PQ节点，当某些发电机的出力P、Q给定时，也可作为PQ节点。PQ节点上的发电机称之为PQ机（或PQ给定型发电机）。在潮流计算中，系统大部分节点属于PQ节点。
2
m 1
Δ P1 e m 1 Δ Q 1 e m 1 Δ Pm e m -1 Δ Q m e m -1 Δ Pm 1 e m -1 Δ U
2
m 1
Δ P1 f m 1 Δ Q 1 f m 1 Δ Pm f m -1 Δ Q m f m -1 Δ Pm 1 f m -1 Δ U
电力系统分析
3.2牛顿－拉夫逊法潮流计算
方程式（15）和（16）具备方程组（12）的形式：
Δ W -J Δ U
（17）
Δe 1 Δf 1 Δe m Δf m ΔU Δe m 1 Δf m 1 Δe n - 1 Δf n - 1

前推回代潮流计算程序

前推回代潮流计算程序潮流计算是电力系统分析的基本方法之一，用于计算电力系统中各节点的电压、功率和电流等参数，并确定各设备的潮流分布情况。

在电力系统规划，运行和故障分析等方面都有广泛应用。

过去，潮流计算主要依赖于人工计算，需要大量的手工作业和时间。

随着计算机技术的发展，前推回代法（FDP）逐渐成为常用的潮流计算方法之一前推回代法是一种解决非线性方程组的数值迭代方法，其基本思想是根据节点电压相等和节点功率平衡的条件，将电力系统的潮流计算问题转化为求解非线性方程组的问题。

其核心是通过迭代计算，不断逼近方程的解。

前推回代法的步骤如下：1.建立节点潮流计算方程组：将电力系统的节点电压和节点功率平衡方程表示为非线性方程组。

其中节点电压方程是根据节点注入功率和导纳矩阵计算得到的，节点功率平衡方程是根据节点注入功率与节点出力功率之间的关系建立的。

2.初始化迭代计算：给定节点电压和相位的初始猜测值。

3.前推计算：从发电节点开始，根据节点电压的计算公式，逐个迭代计算各节点的电压和相位值，直到达到收敛条件。

4.回代计算：根据回代公式，从负荷节点开始，逐个迭代计算各节点的潮流值，直到达到收敛条件。

5.收敛判断：根据设定的收敛判断条件，判断迭代计算是否收敛。

如果未达到收敛条件，则返回第3步进行迭代计算；如果达到收敛条件，则结束计算。

前推回代法的优势在于可以较为准确地计算非线性电力系统的潮流分布情况。

与传统的手工计算相比，前推回代法具有计算速度快、准确度高、可靠性好等优点。

而且，通过计算机编程实现前推回代法，可以大大减少人工计算的工作量，提高计算效率。

然而，前推回代法也存在一些问题和限制。

首先，前推回代法在计算过程中需要不断迭代计算，迭代的次数与收敛的速度成正比。

如果系统存在严重的不平衡或不稳定问题，可能导致计算结果不收敛或收敛速度很慢。

其次，前推回代法在计算大规模电力系统时，可能会面临存储和计算能力的限制。

因此，在实际应用中，需要根据具体情况选择合适的潮流计算方法。

潮流计算程序

%潮流计算fprintf('开始潮流计算\n');fprintf('请输入待求网络的相应参数\n');%参数输入部分n=input('网络中的节点数:n=');L=input('网络中的支路数:L=');ss=input('平衡节点ss=');pr=input('误差精度:pr=');X1=input('支路参数:X1=');X2=input('节点参数:X2=');X=input('节点号和对地参数:X=');fprintf('参数输入部分结束\n\n');Y=zeros(n);%置迭代次数mm=1;%创建节点导纳矩阵for i=1:Lif X1(i,6)==0 %不含变压器的支路p=X1(i,1);q=X1(i,2);Y(p,q)=Y(p,q)-1/X1(i,3);Y(q,p)=Y(p,q);Y(p,p)=Y(p,p)+1/X1(i,3)+0.5*X1(i,4); Y(q,q)=Y(q,q)+1/X1(i,3)+0.5*X1(i,4); else %含有变压器的支路p=X1(i,1);q=X1(i,2);Y(p,q)=Y(p,q)-1/(X1(i,3)*X1(i,5)); Y(q,p)=Y(p,q);Y(p,p)=Y(p,p)+1/X1(i,3);Y(q,q)=Y(q,q)+1/(X1(i,5)^2*X1(i,3)); endendY;OrgS=zeros(2*n-2,1);DetaS=zeros(2*n-2,1); %将OrgS、DetaS初始化%创建OrgS，用于存储初始功率参数h=0;j=0;for i=1:n %对PQ节点的处理if i~=ss&&X2(i,6)==2h=h+1;for j=1:nOrgS(2*h-1,1)=OrgS(2*h-1,1)+real(X2(i,3))*(real(Y(i,j))*real(X2(j,3)) -imag(Y(i,j))*imag(X2(j,3)))+imag(X2(i,3))*(real(Y(i,j))*imag(X2(j,3) )+imag(Y(i,j))*real(X2(j,3)));OrgS(2*h,1)=OrgS(2*h,1)+imag(X2(i,3))*(real(Y(i,j))*real(X2(j,3))-ima g(Y(i,j))*imag(X2(j,3)))-real(X2(i,3))*(real(Y(i,j))*imag(X2(j,3))+im ag(Y(i,j))*real(X2(j,3)));endendendfor i=1:n %对PV节点的处理，注意这时不可再将h初始化为0if i~=ss&&X2(i,6)==0h=h+1;for j=1:nOrgS(2*h-1,1)=OrgS(2*h-1,1)+real(X2(i,3))*(real(Y(i,j))*real(X2(j,3)) -imag(Y(i,j))*imag(X2(j,3)))+imag(X2(i,3))*(real(Y(i,j))*imag(X2(j,3) )+imag(Y(i,j))*real(X2(j,3)));OrgS(2*h,1)=OrgS(2*h,1)+imag(X2(i,3))*(real(Y(i,j))*real(X2(j,3))-ima g(Y(i,j))*imag(X2(j,3)))-real(X2(i,3))*(real(Y(i,j))*imag(X2(j,3))+im ag(Y(i,j))*real(X2(j,3)));endendendOrgS;%创建PVU 用于存储PV节点的初始电压PVU=zeros(n-h-1,1);t=0;for i=1:nif X2(i,6)==0t=t+1;PVU(t,1)=X2(i,3);endendPVU;%创建DetaS，用于存储有功功率、无功功率和电压幅值的不平衡量h=0;for i=1:n %对PQ节点的处理if i~=ss&&X2(i,6)==2h=h+1;DetaS(2*h-1,1)=real(X2(i,2))-OrgS(2*h-1,1);DetaS(2*h,1)=imag(X2(i,2))-OrgS(2*h,1);endendt=0;for i=1:n %对PV节点的处理，注意这时不可再将h初始化为0if i~=ss&&X2(i,6)==0h=h+1;t=t+1;DetaS(2*h-1,1)=real(X2(i,2))-OrgS(2*h-1,1);DetaS(2*h,1)=real(PVU(t,1))^2+imag(PVU(t,1))^2-real(X2(i,3))^2-imag(X 2(i,3))^2;endendDetaS;%创建I，用于存储节点电流参数i=zeros(n-1,1);h=0;for i=1:nif i~=ssh=h+1;I(h,1)=(OrgS(2*h-1,1)-OrgS(2*h,1)*sqrt(-1))/conj(X2(i,3));endendI;%创建Jacbi(雅可比矩阵)Jacbi=zeros(2*n-2);h=0;k=0;for i=1:n %对PQ节点的处理if X2(i,6)==2h=h+1;for j=1:nif j~=ssk=k+1;if i==j %对角元素的处理Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(X2(i,3))+real(Y(i,j))*imag(X2(i ,3))+imag(I(h,1));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(X2(i,3))+imag(Y(i,j))*imag(X2(i,3) )+real(I(h,1));Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k)+2*real(I(h,1));Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1)-2*imag(I(h,1)); else %非对角元素的处理Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(X2(i,3))+real(Y(i,j))*imag(X2(i ,3));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(X2(i,3))+imag(Y(i,j))*imag(X2(i,3) );Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k);Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1);endif k==(n-1) %将用于内循环的指针置于初始值，以确保雅可比矩阵换行k=0;endendendendendk=0;for i=1:n %对PV节点的处理if X2(i,6)==0h=h+1;for j=1:nif j~=ssk=k+1;if i==j %对角元素的处理Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(X2(i,3))+real(Y(i,j))*imag(X2(i ,3))+imag(I(h,1));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(X2(i,3))+imag(Y(i,j))*imag(X2(i,3) )+real(I(h,1));Jacbi(2*h,2*k-1)=2*imag(X2(i,3));Jacbi(2*h,2*k)=2*real(X2(i,3));else %非对角元素的处理Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(X2(i,3))+real(Y(i,j))*imag(X2(i ,3));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(X2(i,3))+imag(Y(i,j))*imag(X2(i,3) );Jacbi(2*h,2*k-1)=0;Jacbi(2*h,2*k)=0;endif k==(n-1) %将用于内循环的指针置于初始值，以确保雅可比矩阵换行k=0;endendendendendJacbi;%求解修正方程，获取节点电压的不平衡量DetaU=zeros(2*n-2,1);DetaU=inv(Jacbi)*DetaS;DetaU;%修正节点电压j=0;for i=1:n %对PQ节点处理if X2(i,6)==2j=j+1;X2(i,3)=X2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);endendfor i=1:n %对PV节点的处理if X2(i,6)==0j=j+1;X2(i,3)=X2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);endendX2;%开始循环********************************************************************* *while abs(max(DetaU))>prOrgS=zeros(2*n-2,1); %!!!初始功率参数在迭代过程中是不累加的，所以在这里必须将其初始化为零矩阵h=0;j=0;for i=1:nif i~=ss&&X2(i,6)==2h=h+1;for j=1:nOrgS(2*h-1,1)=OrgS(2*h-1,1)+real(X2(i,3))*(real(Y(i,j))*real(X2(j,3)) -imag(Y(i,j))*imag(X2(j,3)))+imag(X2(i,3))*(real(Y(i,j))*imag(X2(j,3) )+imag(Y(i,j))*real(X2(j,3)));OrgS(2*h,1)=OrgS(2*h,1)+imag(X2(i,3))*(real(Y(i,j))*real(X2(j,3))-ima g(Y(i,j))*imag(X2(j,3)))-real(X2(i,3))*(real(Y(i,j))*imag(X2(j,3))+im ag(Y(i,j))*real(X2(j,3)));endendendfor i=1:nif i~=ss&&X2(i,6)==0h=h+1;for j=1:nOrgS(2*h-1,1)=OrgS(2*h-1,1)+real(X2(i,3))*(real(Y(i,j))*real(X2(j,3)) -imag(Y(i,j))*imag(X2(j,3)))+imag(X2(i,3))*(real(Y(i,j))*imag(X2(j,3) )+imag(Y(i,j))*real(X2(j,3)));OrgS(2*h,1)=OrgS(2*h,1)+imag(X2(i,3))*(real(Y(i,j))*real(X2(j,3))-ima g(Y(i,j))*imag(X2(j,3)))-real(X2(i,3))*(real(Y(i,j))*imag(X2(j,3))+im ag(Y(i,j))*real(X2(j,3)));endendendOrgS;%创建DetaSh=0;for i=1:nif i~=ss&&X2(i,6)==2h=h+1;DetaS(2*h-1,1)=real(X2(i,2))-OrgS(2*h-1,1);DetaS(2*h,1)=imag(X2(i,2))-OrgS(2*h,1);endendt=0;for i=1:nif i~=ss&&X2(i,6)==0h=h+1;t=t+1;DetaS(2*h-1,1)=real(X2(i,2))-OrgS(2*h-1,1);DetaS(2*h,1)=real(PVU(t,1))^2+imag(PVU(t,1))^2-real(X2(i,3))^2-imag(X 2(i,3))^2;endendDetaS;%创建Ii=zeros(n-1,1);h=0;for i=1:nif i~=ssh=h+1;I(h,1)=(OrgS(2*h-1,1)-OrgS(2*h,1)*sqrt(-1))/conj(X2(i,3));endendI;%创建JacbiJacbi=zeros(2*n-2);h=0;k=0;for i=1:nif X2(i,6)==2h=h+1;for j=1:nif j~=ssk=k+1;if i==jJacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(X2(i,3))+real(Y(i,j))*imag(X2(i ,3))+imag(I(h,1));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(X2(i,3))+imag(Y(i,j))*imag(X2(i,3) )+real(I(h,1));Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k)+2*real(I(h,1));Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1)-2*imag(I(h,1)); elseJacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(X2(i,3))+real(Y(i,j))*imag(X2(i ,3));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(X2(i,3))+imag(Y(i,j))*imag(X2(i,3) );Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k);Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1);endif k==(n-1)k=0;endendendendendk=0;for i=1:nif X2(i,6)==0h=h+1;for j=1:nif j~=ssk=k+1;if i==jJacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(X2(i,3))+real(Y(i,j))*imag(X2(i ,3))+imag(I(h,1));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(X2(i,3))+imag(Y(i,j))*imag(X2(i,3) )+real(I(h,1));Jacbi(2*h,2*k-1)=2*imag(X2(i,3));Jacbi(2*h,2*k)=2*real(X2(i,3));elseJacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(X2(i,3))+real(Y(i,j))*imag(X2(i ,3));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(X2(i,3))+imag(Y(i,j))*imag(X2(i,3) );Jacbi(2*h,2*k-1)=0;Jacbi(2*h,2*k)=0;endif k==(n-1)k=0;endendendendendJacbiDetaU=zeros(2*n-2,1);DetaU=inv(Jacbi)*DetaS;DetaU%修正节点电压j=0;for i=1:nif X2(i,6)==2j=j+1;X2(i,3)=X2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1); endendfor i=1:nif X2(i,6)==0j=j+1;X2(i,3)=X2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1); endendX2mm=mm+1; %迭代次数加1endmm。

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

X是节点号和对地参数矩阵，第一列为节点号；第二列为对地情况，“0”表示没有对地。

2.基本原理：采用牛顿-—拉夫逊潮流算法（直角坐标系）主要变量定义：B1—支路参数矩阵；B2—节点参数矩阵；X—节点号和对地参数；S—节点功率数值，序号为奇数的储存有功功率，偶数的储存无功功率；DS—序号为奇数的储存有功功率不平衡量，对PQ节点，序号为偶数的储存无功功率不平衡量，对PV节点，序号为偶数的储存电压不平衡量；I—节点电流数值；DetaU—节点电压修正值；Number—迭代次数；dS—线路上损耗的功率；3.程序流程框图如下：4.计算结果如下（完整计算结果见文本文档）其中：3节点运行结果如下：请输入节点数:n=3请输入支路数:n1=3请输入平衡节点号:bl=3请输入误差精度:ac=0.000001请输入支路参数：B1=[1 2 0.15i 0 1 1;2 3 0.93i 0 1 0;2 3 0.5i 0 1 0]请输入节点参数：B2=[0.9 0 1 0 0 3;0 0 1 0 0 2;0 0 0.9001 1.05 0 1]节点号和对地参数:X=[1 0;2 0;3 0]Y =0 - 6.6667i 0 + 6.6667i 00 + 6.6667i 0 - 9.7419i 0 + 3.0753i0 0 + 3.0753i 0 - 3.0753iU =0.96941.00000.9001各节点的功率S为（按节点号排列）：0 - 0.1976i0 + 0.5110i0 - 0.2765i各条支路的首端功率Si为（同B1顺序一样）：0 + 0.2038i0 + 0.1074i0 + 0.1998i各条支路的末端功率Sj为（同B1顺序一样）：0 - 0.1976i0 - 0.0967i0 - 0.1798i各条支路的功率损耗dS为（同B1顺序一样）：0 + 0.0062i0 + 0.0200i0 + 0.0200i9节点运行结果如下：请输入节点数:n=9请输入支路数:n1=9请输入平衡节点号:bl=1请输入误差精度:ac=0.001请输入支路参数：B1=[1 4 0+0.0576i 0 1 1;2 7 0+0.0625i 0 1 1;3 9 0+0.0586i 0 1 1;4 5 0.01+0.085i 0.176i 1 0;4 6 0.017+0.092i 0.158i 1 0;5 7 0.032+0.161i 0.306i 1 0;6 9 0.039+0.17i 0.358i 1 0;7 8 0.0085+0.072i 0.149i 1 0;8 9 0.0119+0.1008i 0.209i 1 0]请输入节点参数：B2=[0 0 1.04 1.05 0 1;1.63+0i 0 1.025 1.05 0 3;0.85+0i 0 1.025 1.05 0 3;0 0 1 0 0 2;0 -1.25-0.50i 1 0 0 2;0 -0.90-0.30i 1 0 0 2;0 0 1 0 0 2;0 -1-0.35i 1 0 0 2;0 0 1 0 0 2]节点号和对地参数:X=[1 0;2 0;3 0;4 0;5 0;6 0;7 0;8 0;9 0]Y =Columns 1 through 50 -17.3611i 0 0 0 +17.3611i 00 0 -16.0000i 0 0 00 0 0 -17.0648i 0 00 +17.3611i 0 0 3.3074 -39.3089i -1.3652 +11.6041i0 0 0 -1.3652 +11.6041i 2.5528 -17.3382i0 0 0 -1.9422 +10.5107i 00 0 +16.0000i 0 0 -1.1876 + 5.9751i0 0 0 0 00 0 0 +17.0648i 0 0Columns 6 through 90 0 0 00 0 +16.0000i 0 00 0 0 0 +17.0648i-1.9422 +10.5107i 0 0 00 -1.1876 + 5.9751i 0 03.2242 -15.8409i 0 0 -1.2820 + 5.5882i0 2.8047 -35.4456i -1.6171 +13.6980i 00 -1.6171 +13.6980i 2.7722 -23.3032i -1.1551 + 9.7843i-1.2820 + 5.5882i 0 -1.1551 + 9.7843i 2.4371 -32.1539i -0.0168 - 0.1328iU =1.04003.2481 - 0.9555i3.2076 - 0.8593i3.2428 - 0.5596i3.2235 - 0.7453i1.8842 - 0.3587i2.3853 - 0.6321i2.4410 - 0.8344i3.1907 - 0.7505i各节点的功率S为（按节点号排列）：1.0e+002 *0.1010 - 0.3977i-0.0361 + 0.4978i-0.0571 + 0.0252i0.0687 + 1.7360i-0.0214 + 0.1542i-0.0808 - 0.4329i0.0531 - 0.5394i-0.1579 - 0.1256i0.2074 + 0.4176i-0.00180.0071-0.00710.00410.0044各条支路的首端功率Si为（同B1顺序一样）：-1.0105e+001 +1.2945e+002i3.6128 -36.2000i5.7087 - 2.3123i7.0539 - 2.2715i9.9160 +46.4206i4.8627 +14.6652i-1.8285 -15.5465i5.7060 - 4.7477i-10.1578 -16.9569i各条支路的末端功率Sj为（同B1顺序一样）：10.1046 -39.7723i-3.6128 +49.7840i-5.7087 + 2.5192i-7.0064 + 0.7595i-6.2529 -27.7428i-4.0131 -12.9970i4.2137 +23.3616i-5.6348 + 4.4012i10.8152 +20.7070i各条支路的功率损耗dS为（同B1顺序一样）：-0.0000 +89.6771i0 +13.5840i0 + 0.2069i0.0476 - 1.5120i3.6631 +18.6778i0.8496 + 1.6681i2.3852 + 7.8152i0.0712 - 0.3465i0.6574 + 3.7501i39节点运行结果如下：请输入节点数:n=39请输入支路数:n1=46请输入平衡节点号:bl=31请输入误差精度:ac=0.1请输入支路参数：B1=[1 2 0.00350+0.04110i 0.6987i 1 0;1 39 0.00100+0.02500i 0.7500i 1 0;2 3 0.00130+0.01510i 0.2572i 1 0;2 25 0.00700+0.00860i 0.1460i 1 0;30 2 0.00000+0.01810i 0.00 102.5 1;3 4 0.00130+0.02130i 0.2214i 1 0;3 18 0.00110+0.01330i 0.2138i 1 0;4 5 0.00080+0.01280i 0.1342i 1 0;4 14 0.00080+0.01290i 0.1382i 1 0;5 6 0.00020+0.00260i -0.02170 1 0;5 8 0.00080+0.01120i 0.1476i 1 0;6 7 0.00060+0.00920i 0.1130i 1 0;6 11 0.00070+0.00820i 0.1389i 1 0;31 6 0.00000+0.02500i 0.00 107.000 1;7 8 0.00040+0.00460i 0.0780i 1 0;8 9 0.00230+0.03630i 0.3804i 1 0;9 39 0.00100+0.02500i 1.2000i 1 0;10 11 0.00040+0.00430i 0.0729i 1 0;10 13 0.00040+0.00430i 0.0729i 1 0;32 10 0.00000+0.02000i 0.00 107.000 1;11 12 0.00160+0.04350i 0.00 100.600 1;13 12 0.00160+0.04350i 0.00 100.600 1;13 14 0.00090+0.01010i 0.1723i 1 0;14 15 0.00180+0.02170i 0.3660i 1 0;15 16 0.00090+0.00940i 0.3420i 1 0;16 17 0.00070+0.00890i 0.1342i 1 0;16 19 0.00160+0.01950i 0.3040i 1 0;16 21 0.00080+0.01350i -0.12740 1 0;16 24 0.00030+0.00590i 0.0680i 1 0;17 18 0.00070+0.00820i 0.1319i 1 0;17 27 0.00130+0.01730i 0.3216i 1 0;20 19 0.00070+0.01380i 0.00 106.000 1;33 19 0.00070+0.01420i 0.00 107.000 1;34 20 0.00090+0.01800i 0.00 100.900 1;21 22 0.00080+0.01400i 0.2565i 1 0;22 23 0.00060+0.00960i 0.1846i 1 0;35 22 0.00000+0.01430i 0.00 102.500 1;23 24 0.00220+0.03500i 0.3610i 1 0;36 23 0.00050+0.02720i 0.00 100.000 1;25 26 0.00320+0.03230i 0.5130i 1 0;37 25 0.00060+0.02320i 0.00 102.500 1;26 27 0.00140+0.01470i 0.2396i 1 0;26 28 0.00430+0.04740i 0.7802i 1 0;26 29 0.00570+0.06250i 1.0290i 1 0;28 29 0.00140+0.01510i 0.2490i 1 0;38 29 0.00080+0.01560 0.00 102.500 1;] 请输入节点参数：B2=[0 0 1 0 0 2;0 0 1 0 0 2;0 -3.22-0.024i 1 0 0 2;0 -5.00-1.84i 1 0 0 2;0 0 1 0 0 2;0 0 1 0 0 2;0 -2.338-0.84i 1 0 0 2;0 -5.22-1.76i 1 0 0 2;0 0 1 0 0 2;0 0 1 0 0 2;0 0 1 0 0 2;0 -0.085-0.88i 1 0 0 2;0 0 1 0 0 2;0 0 1 0 0 2;0 -3.20-1.53i 1 0 0 2;0 -3.29-0.323i 1 0 0 2;0 0 1 0 0 2;0 -1.58-0.30i 1 0 0 2;0 0 1 0 0 2;0 -6.80-1.03i 1 0 0 2;0 -2.74-1.15i 1 0 0 2;0 0 1 0 0 2;0 -2.475-0.846i 1 0 0 2;0 -3.086+0.922i 1 0 0 2;0 -2.24-0.472i 1 0 0 2;0 -1.39-0.17i 1 0 0 2;0 -2.81-0.755i 1 0 0 2;0 -2.06-0.276i 1 0 0 2;0 -2.835-1.269i 1 0 0 2;2.50+1.3621i 0 1.04750 0 0 3;0 0 1.050 1.05 0 1;6.50+1.7590i 0 0.98310 1.05 0 3;6.32+1.0335i 0 0.99720 1.05 0 3;5.08+1.6440i 0 1.01230 1.05 0 3;6.50+2.0484i 0 1.04930 1.05 0 3;5.60+0.9688i 0 1.06350 1.05 0 3;5.40-0.0444i 0 1.02780 1.05 0 3;8.30+0.1939i 0 1.02650 1.05 0 3;10.0+0.68460i -11.04-2.50i 1.03 1.05 0 3]节点号和对地参数:X=[1 0;2 0;3 0;4 0;5 0;6 0;7 0;8 0;9 0;10 0;11 0;12 0;13 0;14 0;15 0;16 0;17 0;18 0;19 0;20 0;21 0;22 0;23 0;24 0;25 0;26 0;27 0;28 0;29 0;30 0;31 0;32 0;33 0;34 0;35 0;36 0;37 0;38 0;39 0]Y =1.0e+002 *Columns 1 through 50.0365 - 0.6337i -0.0206 + 0.2416i 0 0 0-0.0206 + 0.2416i 0.6465 - 1.5929i -0.0566 + 0.6574i 0 00 -0.0566 + 0.6574i 0.1469 - 1.8684i -0.0285 + 0.4677i 00 0 -0.0285 + 0.4677i 0.1251 - 2.0157i -0.0486 + 0.7782i0 0 0 -0.0486 + 0.7782i 0.4061 - 5.4887i0 0 0 0 -0.2941 + 3.8235i0 0 0 0 00 0 0 0 -0.0635 + 0.8883i0 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 -0.0479 + 0.7722i 00 0 0 0 00 0 0 0 00 0 0 0 00 0 -0.0618 + 0.7468i 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 -0.5693 + 0.6994i 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 + 0.0054i 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 0-0.0160 + 0.3994i 0 0 0 0Columns 6 through 100 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 0-0.2941 + 3.8235i 0 -0.0635 + 0.8883i 0 00.4679 - 6.1153i -0.0706 + 1.0824i 0 0 0-0.0706 + 1.0824i 0.2582 - 3.2390i -0.1876 + 2.1576i 0 00 -0.1876 + 2.1576i 0.2685 - 3.3173i -0.0174 + 0.2744i 00 0 -0.0174 + 0.2744i 0.0334 - 0.6658i 00 0 0 0 0.4290 - 4.6106i-0.1034 + 1.2107i 0 0 0 -0.2145 + 2.3056i0 0 0 0 00 0 0 0 -0.2145 + 2.3056i0 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 + 0.0037i 0 0 0 00 0 0 0 0 + 0.0047i0 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 -0.0160 + 0.3994i 0Columns 11 through 150 0 0 0 00 0 0 0 00 0 0 -0.0479 + 0.7722i 00 0 0 0 0-0.1034 + 1.2107i 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 0-0.2145 + 2.3056i 0 -0.2145 + 2.3056i 0 00.3263 - 3.7448i -0.0001 + 0.0023i 0 0 0-0.0001 + 0.0023i 0.0000 - 0.0000i -0.0001 + 0.0023i 0 00 -0.0001 + 0.0023i 0.3105 - 3.5163i -0.0875 + 0.9823i 00 0 -0.0875 + 0.9823i 0.1734 - 2.2088i -0.0380 + 0.4577i0 0 0 -0.0380 + 0.4577i 0.1389 - 1.5083i0 0 0 0 -0.1009 + 1.0542i0 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 0Columns 16 through 200 0 0 0 00 0 0 0 00 0 -0.0618 + 0.7468i 0 00 0 0 0 00 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 0-0.1009 + 1.0542i 0 0 0 00.3596 - 5.1047i -0.0878 + 1.1167i 0 -0.0418 + 0.5094i 0-0.0878 + 1.1167i 0.2344 - 2.8992i -0.1034 + 1.2107i 0 00 -0.1034 + 1.2107i 0.1651 - 1.9557i 0 0-0.0418 + 0.5094i 0 0 0.0418 - 0.5080i -0.0003 + 0.0068i0 0 0 -0.0003 + 0.0068i 0.0367 - 0.7228i-0.0437 + 0.7381i 0 0 0 00 0 0 0 00 0 0 0 0-0.0860 + 1.6905i 0 0 0 00 0 0 0 00 0 0 0 00 -0.0432 + 0.5748i 0 0 00 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 -0.0003 + 0.0066i 00 0 0 0 -0.0003 + 0.0055i0 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 0Columns 21 through 250 0 0 0 00 0 0 0 -0.5693 + 0.6994i0 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 0-0.0437 + 0.7381i 0 0 -0.0860 + 1.6905i 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00.0838 - 1.4488i -0.0407 + 0.7120i 0 0 0-0.0407 + 0.7120i 0.1055 - 1.7474i -0.0649 + 1.0376i 0 00 -0.0649 + 1.0376i 0.0827 - 1.3195i -0.0179 + 0.2846i 00 0 -0.0179 + 0.2846i 0.1038 - 1.9730i 00 0 0 0 0.5997 - 1.0027i0 0 0 0 -0.0304 + 0.3066i0 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 + 0.0068i 0 0 00 0 -0.0001 + 0.0037i 0 00 0 0 0 -0.0001 + 0.0042i0 0 0 0 00 0 0 0 0Columns 26 through 300 0 0 0 00 0 0 0 0 + 0.0054i0 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 00 0 0 0 00 0 0 0 00 -0.0432 + 0.5748i 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 0-0.0304 + 0.3066i 0 0 0 00.1280 - 1.3359i -0.0642 + 0.6742i -0.0190 + 0.2092i -0.0145 + 0.1587i 0-0.0642 + 0.6742i 0.1074 - 1.2461i 0 0 0-0.0190 + 0.2092i 0 0.0799 - 0.8607i -0.0609 + 0.6566i 0-0.0145 + 0.1587i 0 -0.0609 + 0.6566i 0.0754 - 0.8089i 00 0 0 0 0 - 0.5525i0 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 -0.0059 00 0 0 0 0Columns 31 through 350 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 + 0.0037i 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 + 0.0047i 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 -0.0003 + 0.0066i 00 0 0 -0.0003 + 0.0055i 00 0 0 0 00 0 0 0 0 + 0.0068i0 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 - 0.4000i 0 0 0 00 0 - 0.5000i 0 0 00 0 0.0346 - 0.7025i 0 00 0 0 0.0277 - 0.5542i 00 0 0 0 0 - 0.6993i0 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 0Columns 36 through 390 0 0 -0.0160 + 0.3994i0 0 0 00 0 0 00 0 0 00 0 0 00 0 0 00 0 0 00 0 0 00 0 0 -0.0160 + 0.3994i0 0 0 00 0 0 00 0 0 00 0 0 00 0 0 00 0 0 00 0 0 00 0 0 00 0 0 00 0 0 00 0 0 00 0 0 00 0 0 0-0.0001 + 0.0037i 0 0 00 0 0 00 -0.0001 + 0.0042i 0 00 0 0 00 0 0 00 0 0 00 0 -0.0059 00 0 0 00 0 0 00 0 0 00 0 0 00 0 0 00 0 0 00.0068 - 0.3675i 0 0 00 0.0111 - 0.4307i 0 00 0 0.6098 00 0 0 0.0319 - 0.7890i U =1.0e+003 *0.1757 - 0.1082i0.0915 - 0.1339i0.0712 - 0.1171i0.0618 - 0.0724i0.0680 - 0.0540i0.0666 - 0.0494i0.0745 - 0.0550i0.0792 - 0.0574i0.1659 - 0.0776i0.0535 - 0.0308i0.0555 - 0.0296i3.8741 +4.1652i0.0514 - 0.0316i0.0546 - 0.0614i0.0497 - 0.1055i0.0496 - 0.1237i0.0547 - 0.1336i0.0605 - 0.1277i0.0504 - 0.1256i-0.0002 + 0.0015i0.0466 - 0.1265i0.0461 - 0.1286i0.0455 - 0.1290i0.0484 - 0.1255i0.0744 - 0.1436i0.0516 - 0.1896i0.0520 - 0.1647i0.0334 - 0.2574i0.0283 - 0.2765i0.0016 - 0.0779i0.00110.0014 - 0.0726i-0.0072 - 0.0790i-0.0062 - 0.0964i0.0011 - 0.0800i-0.0023 - 0.0823i-0.0031 - 0.0825i-4.3537 + 0.6080i0.2219 - 0.0873i各节点的功率S为（按节点号排列）：1.0e+009 *-0.0000 + 0.0000i0.0000 + 0.0000i -0.0000 + 0.0000i -0.0000 - 0.0000i -0.0000 + 0.0000i -0.0000 + 0.0000i -0.0000 - 0.0000i -0.0000 - 0.0000i -0.0000 + 0.0000i -0.0000 + 0.0000i 0.0000 - 0.0001i 0.0002 + 0.0001i 0.0000 - 0.0001i -0.0000 + 0.0000i -0.0000 - 0.0000i -0.0000 + 0.0000i 0.0000 + 0.0000i -0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i -0.0000 - 0.0000i 0.0000 + 0.0000i -0.0000 + 0.0000i -0.0000 + 0.0000i -0.0000 + 0.0000i -0.0000 + 0.0000i -0.0000 + 0.0000i -0.0000 + 0.0000i -0.0000 - 0.0000i -0.0000 + 0.0003i 0.0000 + 0.0000i -0.0000 + 0.0003i 0.0000 + 0.0004i 0.0000 + 0.0005i -0.0000 + 0.0004i 0.0000 + 0.0002i0.0000 + 0.0003i1.1785 + 0.0007i 0.0005 + 0.0008i 0.00000.0000-0.00000.0000-0.00000.0000-0.00000.0000-0.00000.0000-0.00000.0000-0.00000.0000-0.00000.0000-0.00000.0000-0.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.0000各条支路的首端功率Si为（同B1顺序一样）：3.5400e+005 +2.4736e+005i-3.5581e+005 -2.3628e+005i1.0071e+005 +2.6013e+005i2.3682e+005 -1.6375e+005i3.7264e+003 +1.4471e+006i-8.0159e+004 +2.7957e+005i1.4664e+005 -5.0894e+004i-1.1893e+005 +8.1179e+004i-6.8733e+003 +9.6236e+004i-8.1580e+004 +1.3584e+005i-3.9332e+004 -8.2485e+004i-7.6977e+003 -8.7005e+004i-7.5059e+004 +2.1608e+005i-1.9406e+001 +2.7536e+005i-2.5689e+004 -1.0282e+005i-1.0682e+005 -2.1610e+005i-1.2535e+005 -4.1708e+005i-2.9483e+004 -1.4331e+004i2.6881e+004 +1.7017e+004i1.7943e+003 +1.8913e+005i2.7401e+007 +7.4282e+008i2.7399e+007 +7.4282e+008i1.3064e+005 -1.2139e+005i1.1489e+005 -1.2327e+005i7.7164e+004 -2.1257e+005i-2.8359e+004 -1.6599e+005i-1.7425e+003 -1.7302e+004i3.5700e+004 -1.7516e+004i3.8835e+004 -3.2323e+004i-1.2909e+005 +6.6903e+004i1.0067e+005 -2.4294e+005i6.7162e+004 +1.3247e+006i6.6365e+004 +1.2810e+006i-7.0330e+000 +2.1637e+002i9.7641e+003 -2.0069e+004i1.0245e+004 -4.5881e+003i2.4215e+003 +1.2980e+006i-1.4470e+004 +6.4486e+003i1.4051e+004 +6.8360e+005i1.9061e+005 -1.7744e+005i3.1770e+004 +1.1214e+006i-6.2319e+004 +3.1983e+005i1.2274e+005 -2.7770e+005i1.1898e+005 -2.7519e+005i9.8516e+004 -3.3183e+005i4.8839e+006 -7.0589e+005i各条支路的末端功率Sj为（同B1顺序一样）：-3.3804e+005 -8.4052e+004i3.5993e+005 +3.0184e+005i-9.6779e+004 -2.2023e+005i-2.1492e+005 +1.8682e+005i-3.7264e+003 -5.6637e+003i8.6098e+004 -1.8535e+005i-1.4524e+005 +6.3681e+004i1.2077e+005 -5.2863e+004i7.7058e+003 -8.3904e+004i8.2090e+004 -1.2718e+005i4.0208e+004 +9.3498e+004i8.3567e+003 +9.6235e+004i8.0400e+004 -1.5426e+005i19.4058 -26.1530i2.6210e+004 +1.0809e+005i1.2059e+005 +4.2529e+005i1.3051e+005 +4.9198e+005i2.9595e+004 +1.5259e+004i-2.6774e+004 -1.6139e+004i-1.7943e+003 -1.0552e+003i-7.9694e+004 -1.7952e+004i-7.7356e+004 -1.2573e+004i-1.2279e+005 +2.0853e+005i-1.0739e+005 +2.0990e+005i-7.3846e+004 +2.4186e+005i2.9462e+004 +1.7743e+005i1.7620e+003 +1.2057e+004i-3.7913e+004 +1.8782e+004i-3.8792e+004 +3.1940e+004i1.2980e+005 -6.1207e+004i-9.6457e+004 +2.9084e+005i3.9717e+001 +1.3853e+002i-3.5171e+003 -6.0815e+003i2.4385e+001 +1.3068e+002i-9.7460e+003 +1.5660e+004i-1.0242e+004 +1.1971e+003i-2.4215e+003 -7.0055e+003i1.4506e+004 -1.2520e+004i-1.5552e+003 -3.8100e+003i-1.8260e+005 +2.4172e+005i-2.8930e+003 -4.7791e+003i6.6279e+004 -2.8646e+005i-1.1337e+005 +3.3958e+005i-1.0726e+005 +3.4408e+005i-9.6140e+004 +3.3945e+005i2.8547e+005 +7.0589e+005i各条支路的功率损耗dS为（同B1顺序一样）：1.5956e+004 +1.6331e+005i4.1140e+003 +6.5558e+004i3.9344e+003 +3.9904e+004i2.1901e+004 +2.3079e+004i-7.8671e-011 +1.4414e+006i5.9389e+003 +9.4226e+004i 1.4002e+003 +1.2787e+004i 1.8394e+003 +2.8317e+004i 8.3252e+002 +1.2332e+004i 5.0953e+002 +8.6579e+003i 8.7681e+002 +1.1013e+004i6.5898e+002 +9.2305e+003i 5.3409e+003 +6.1812e+004i -1.3866e-011 +2.7534e+005i 5.2034e+002 +5.2760e+003i 1.3774e+004 +2.0919e+005i 5.1660e+003 +7.4900e+004i 1.1258e+002 +9.2738e+002i 1.0692e+002 +8.7808e+002i 7.5033e-012 +1.8808e+005i 2.7321e+007 +7.4280e+008i 2.7322e+007 +7.4281e+008i 7.8440e+003 +8.7133e+004i 7.4944e+003 +8.6625e+004i 3.3182e+003 +2.9296e+004i 1.1030e+003 +1.1435e+004i 1.9504e+001 -5.2453e+003i -2.2133e+003 +1.2653e+003i 4.2504e+001 -3.8293e+002i 7.1612e+002 +5.6968e+003i 4.2119e+003 +4.7899e+004i 6.7202e+004 +1.3248e+006i 6.2848e+004 +1.2749e+006i 1.7352e+001 +3.4704e+002i 1.8043e+001 -4.4090e+003i3.6386e+000 -3.3909e+003i4.4565e-011 +1.2910e+006i 3.5978e+001 -6.0715e+003i 1.2496e+004 +6.7979e+005i 8.0137e+003 +6.4281e+004i2.8877e+004 +1.1166e+006i3.9594e+003 +3.3374e+004i 9.3637e+003 +6.1880e+004i1.1719e+004 +6.8892e+004i2.3756e+003 +7.6167e+003i 5.1694e+006。