电力系统分析课程设计

1、精品资料推荐16电力系统分析课程设计学 院电气工程学院班 级学 号姓 名指导教师 时 间.、八、一前言电力系统潮流计算是研究电力系统稳态运行情况的一种计算， 它根据给定的 运行条件及系统接线情况确定整个电力系统各部分的运行状态：各母线的电压， 各元件中流过的功率， 系统的功率损耗等等。 在电力系统规划的设计和现有电力 系统运行方式的研究中， 都需要利用潮流计算来定量地分析比较供电方案或运行 方式的合理性。 可靠性和经济性。 此外，电力系统潮流计算也是计算系统动态稳 定和静态稳定的基础。 所以潮流计算是研究电力系统的一种很重要和很基础的计 算。随着科学技术的发展， 电力系统变得越来越复杂， 电气

2、工程师掌握一种好的 能对电力系统进行仿真的软件是学习和研究的需要。 与众多专门的电力系统仿真 软件相比，MATLAB软件具有易学、功能强大和开放性好，是电力系统仿真研究 的有力工具。目录1设计题目 31.1系统图的确定 31.2各节点的初值及阻抗参数 42潮流计算 52.1潮流计算概述与发展 52.2复杂电力系统潮流计算 52.3 MATLAB 概述 62.4牛顿-拉夫逊法原理 72.5牛顿-拉夫逊法解决潮流计算问题 82.6计算机潮流计算的步骤 92.7计算机潮流计算流程图 113手算潮流计算 123.1确定节点类型 123.2求节点导纳矩阵Yb 123.3计算各节点功率的修正方程的初始值（

3、不平衡量） 134计算机算法潮流计算 154.1 计算机法潮流计算过程 154.2计算机法潮流计算结果 234.3系统功率分布图 25个人心得 29附录：源程序 30参考文献 401设计题目1.1系统图的确定选择六节点、环网、两电源和多引出的电力系统，简化电力系统图如图 1-1 所示，等值导纳图如图1-2所示。运用以直角坐标表示的牛顿-拉夫逊计算如图 1-1所示系统中的潮流分布。计算精度要求各节点电压的误差或修正量不大于103.7+j1.3j.25HI卜j0.252+j11.8+j4.01.6+j0.8I1_Ij0.25j0.25 0 二 j0.25C=H0.04+j0.251.05:1rr&

4、gt;rb二寸0.25- 1.05:13.7+j1.3j 0. 03U11 . 050。图1-2电力系统等值导纳图1.2各节点的初值及阻抗参数该系统中，节点为平衡节点，保持U1=1.05+j0为定值，节点为PV节点, 其他四个节点都是PQ节点。给定的注入电压标幺值如表 1-1、线路阻抗标幺值 如表1-2、输出功率标幺值如表1-3。表1-1各节点电压标幺值参数UIU2U3U4U5U61.051.001.001.001.001.05表1-2线路、变压器阻抗标幺值线路L2L3L4L5T1T2Y/2阻抗0.06+j0.250.04+j0.250.08+j0.300.1+j0.35j0.03j0.015

5、j0.25表1-3节点输出功率节点功率2+j11.8+j0.41.6+j0.83.7+j1.35注：各PQ节点的电压取1是为了方便计算和最后验证程序的正确性。2潮流计算2.1潮流计算概述与发展电力系统潮流计算也分为离线计算和在线计算两种，前者主要用于系统规划 设计和安排系统的运行方式，后者则用于正在运行系统的经常监视及实时控制。利用电子数字计算机进行电力系统潮流计算从50年代中期就已经开始。在这20年内，潮流计算曾采用了各种不同的方法，这些方法的发展主要围绕着对 潮流计算的一些基本要求进行的。牛顿-拉夫逊法作为一种实用的，有竞争力的电力系统潮流计算方法，是在 应用了稀疏矩阵技巧和高斯消去法求修

6、正方程后。牛顿-拉夫逊法是求解非线性 代数方程有效的迭代计算。2.2复杂电力系统潮流计算电力系统潮流计算是对复杂电力系统正常和故障条件下稳态运行状态的计算。潮流计算的目标是求取电力系统在给定运行方式下的节点电压和功率分布， 用以检查系统各元件是否过负荷、 各点电压是否满足要求、功率的分布和分配是 否合理以及功率损耗等。对现有电力系统的运行和扩建，对新的电力系统进行规 划设计都是以潮流计算为基础。潮流计算结果的用途，例如用于电力系统稳定研究、安全估计或最优潮流等 也对潮流计算的模型和方法有直接影响。节点类型：（1）PV节点：柱入有功功率P为给定值，电压也保持在给定数值。（2）PQ节点：诸如有功功

7、率和无功功率是给定的。（3）平衡节点：用来平衡全电网的功率。选一容量足够大的发电机担任平 衡全电网功率的职责。 平衡节点的电压大小与相位是给定的， 通常以它的相角为 参考量，即取其电压相角为 0。一个独立的电力网中只设一个平衡点。基本步骤：（1）形成节点导纳矩阵；（2）将各节点电压设初值 U；（3）将节点初值代入相关求式，求出修正方程式的常数项向量；（4）将节点电压初值代入求式，求出雅可比矩阵元素；（5）求解修正方程，求修正向量；（6）求取节点电压的新值；（7）检查是否收敛，如不收敛，则以各节点电压的新值作为初值自第3 步重新开始进行狭义次迭代，否则转入下一步；（8）计算支路功率分布，PV节点

8、无功功率和平衡节点柱入功率。2.3 MATLAB 概述目前电子计算机已广泛应用于电力系统的分析计算， 潮流计算是其基本应用 之一。现有很多潮流计算方法。对潮流计算方法有五方面的要求：（ 1）计算速度快；（ 2 ）内存需要少；（ 3）计算结果有良好的可靠性和可信性；（ 4）适应性好，亦即能处理变压器变比调整、系统元件的不同描述和与其 它程序配合的能力强；（ 5 ）简单。MATLAB 是一种交互式、面向对象的程序设计语言，广泛应用于工业界与学 术界，主要用于矩阵运算，同时在数值分析、自动控制模拟、数字信号处理、动 态分析、绘图等方面也具有强大的功能。MATLAB 程序设计语言结构完整，且具有优良的

9、移植性，它的基本数据元素 是不需要定义的数组。 它可以高效率地解决工业计算问题， 特别是关于矩阵和矢 量的计算。MATLA与C语言和FoRTRA语言相比更容易被掌握。通过 M语言,可以用类似数学公式的方式来编写算法，大大降低了程序所需的难度并节省了时 间，从而可把主要的精力集中在算法的构思而不是编程上。2.4牛顿-拉夫逊法原理£仏(O3X2XnO)Xi),2O).nO)假设有n个联立的非线性代数方程：fi(Xi,X2Xn)0f2(xi,x2Xn)0fn(xi,x2Xn)0NZ(O)假设以给出各变量的初值Xi，(0)X2x(°)，Xn ，令其分别为个变量的修正,使满足以上方程

10、，所以：将上式中的n个多元函数在初始值附近分别展开成泰勒级数，并略去含有(OXn的二次及以上阶次的各项，便得:fi(Xi(O),x2O),XnO)f2(Xi(O),x2O),XnO)fn(x(O),2O),XnO)XiMXi(0)Xi(0)Xi )iX2f2X2(0)X2(O)x2 )iX3f2X3(0)Xi )nX2(0)2 )nX3方程可写成:l1O -l1O.l1Of (X(O) X(O) f1 (X1 ,X2 ,.XnO)1<2Xn(O)X1f (X(O) X(O).XnO)f 2O -f2O.f2*O(O)X2fn(,2O),.XnO)、f<1nX2 fnXn fnX(O

11、)入nO -O.OX1<2Xnf1 f1 f1v(0) v(0)以上方程是对于修正量X1 ,x2 ,(O)Xn的线性方程组，称为牛顿法的修正方程，可解出XF，X)，x(°)Xn。对初始近似解进行修正:(1)Xi(0)Xi(0)Xi(i=1,2, ，n)反复迭代，在进行k+1次迭代时，从求解修正方程式:f(Ik),x2k),Xnk) f2(Ik),x2k),Xnk) fn(Xl(k),x2k),Xnk)f1f11kkk<1<2Xnf2f2f2kkk<1<2Xnfnfnfnk k"kX1X2Xn(k)(k)得到修正量 X1 ,X2,Xr)(k)(k

12、 1)(k)Xn ,对各量进行修正XiXi(i=1,2, ， n )迭代过程一直进行到满足收敛判据ma fi (XIk) ,x2k),，Xnk) |12.5牛顿-拉夫逊法解决潮流计算问题节点总数为n； PQ节点有m; PV节点有n-m-1 ,平衡节点有1个，节点编 号按照先PQ节点，再PV节点，最后平衡节点的顺序进行编号，即:1，2, m 为PQ节点；m+1, m+2，n-1为PV节点；n为平衡节点。可形成结点导纳矩 阵。导纳矩阵元素可表示为Yj gij jBij ,本文中节点电压以直角坐标形式表示，即Vieijfi。由此下列公式可求出Pi , QinnPie (Gijej Bijfj) f

13、i (Gijfj Bijej)j 1j 1nnQi fi (Gijej Bijfj) ei (Gij f j Bij ej )j 1j 1假设系统中的第1,2,，m号节点为PQ节点，第i个节点的给定功率为PiS和 Qis ，对该节点可列方程：Pi PisPiQiQisQinPis ei(GijejBi fj )j1nQisfi(Gij ejBi f j )j1nfi(Gi f jj1nei(Gi f jj1Bij ej )(i Bijej),m)假设系统中的第m+1, m+2n-1号节点为PV节点，则对其中每一个节点可列方程 :PiPisPiPisVi2nei(GijejBi f j)j1Vi

14、s2Vi2Vis2(ei2ni(Gi f jj1fi2)Bij ej ) (i m 1,m 2, ,n 1)第 n 号节点为平衡节点，其电压为是给定的Vnjf n ，故不参加迭代。修正方程 W J V 可写成分块矩阵的形式：W1J11J12J1,n 1V1W2J21J22J2,n 1V2W3J31J32J3,n 1*V3Wn 1Jn 1,1Jn 1,2Jn 1,n 1Vn 1通过反复求解修正方程，解出各节点的未知量， 再通过收敛判据判定是否已为真值。从而求得PQ节点的电压V及相角的真值，PV节点的Q 真值，平 衡节点的P、Q真值，以上即为牛顿-拉夫逊迭代法的潮流计算过程，其优点为计算精确，运行

15、速度快。其中的各个环节都可通过MATLA程序来实现2.6 计算机潮流计算的步骤1) 对电力网络的所有参数设初值 , 包括电压、相角、有功、无功等。2) 处理非标准变比支路 , 使其变成标准变比为 1 的变压器支路。3) 形成节点导纳矩阵 Y。(4)计算有功功率的不平衡量 Pi,从而求出R/UiQ 123, n,i S)5)根据节点的类型形成 J。(6) 解修正方程式,求各节点的电压的变化量 ei(i=1,2,3n,i S)(7) 求各节点相角的新值 ei=ei+ ei (i=1,2,3. n,is)(8) 计算无功功率的不平衡量 Qi,从而求出 Q/Ui (i=l,2,3.n,iS)(9) 解

16、修正方程式 , 求各节点的电压大小的变化量 Ui (i=1,2,3.,n,i s) 。(10) 求各节点的电压大小的新值 Ui Ui Ui (i=1,2,3.,n,i s)( 1 1 )运用个节点的电压的新值自第四步开始下一次迭代。计算平衡节点的功率和线路功率。其中平衡节点的功率的计算公式为. i n . .SS U S Y Ui Pi JQi i 1 si线路上的功率为 :SijUi Iij UiUi yi0(Ui Uj)yij PijJQijSjiU jI ji UjUi yj0 (Uj Ui)yji PjiJQji从而线路上的损耗的功率为SijSij S jiPijJ Qij2.7 计算

17、机潮流计算流程图选用牛顿 -拉夫逊法进行计算机潮流计算。 牛顿 - 拉夫逊法潮流计算流程图 如图 2-1 。图2-1牛顿-拉夫逊法潮流计算流程图3手算潮流计算3.1确定节点类型节点为平衡节点节点为PV节点，其他四个节点都是PQ节点3.2求节点导纳矩阵Yb由图1可知，该系统以串联支路的导纳标幺值和对地并联导纳标幺值得等 值电路如图2所示。可得图2的等值节点导纳矩阵。对角线上的元素为：y-YtiKK1YtiKj 33.33y22YtiK1 KYK2 YTIy 20y30y23y251.53-j37.42y33y,30y,40 y23y341.74j6.39丫44y40y,50 y34y451.58

18、j5.25y55丄YKYTII-KYK2 YT2y'20y50y 25y451.37j 66.51y66丄YKYTIK-IYK YT2j66.67非对角线上的元素为:y12y21j31.75y13y310；14410；15510；1661°y23y320.91j3.78y24y420；26y62 0y25y520.62j3.90y34y430.83j3.11y35y530；36y630y45y540.75j2.6446640y56y65j63.49所以节点导纳矩阵为:Y=j3333j31.750000j31.751.53 j37.420.91 j3.7800.62 j3.90

19、000.91 j3.781.74 j6.390.83 j3.1100000.83 j3.111.58 j5.250.75 j2.64000.62 j3.9000.75 j2.641.37 j6651j63490000j6349j66673.3计算各节点功率的修正方程的初始值（不平衡量）设各节点电压初始值为：?U 1 1.05 j0>-A* = .oo+rZ） = IOO+ U-=i=IOO ÷= IO5 *o根据上述N-R法的求解过程，将各节点电压的初始值代入（1）式和（2）式,进行潮流计算，得节点功率和节点电压:-1.53-41.020.913.780-33.811.533.

20、78-0.9100.913.78-1.746.890.833.780.915.891.743.11000.833.111.58J= 003.110.834.750.623.90000.753.900.62002.640000000000雅可比矩阵：修正各节点电压：e2(0)0.095f2(O)0.123e3(0)0.061f3O)0.604e4(0)0.131f4O)0.638e5(0)0.096f 5O)0.159e6(0)0f (60)0.880 0.623.900003.90-0.64003.1100000.8300005.75 0.752.64001.58 2.640.75002.64

21、 1.3873.21063.490.75 59.81 1.38 63.4900066.67063.490002.10e2(1)1.095f2(1)0.123e(1)30.939f (1)30.640e(1)40.869f (1)f40.638e(1) e51.096f (1)f50.159e(1)61f (1)60.88P22; Q22.61;P31.8; Q30.10;P41.6; Q40.30;P53.7; Q55.40;P65.0; U620.00;4 计算机算法潮流计算4.1 计算机法潮流计算过程导纳矩阵 Y=0 -30.2343i 0 +31.7460i 0 0 00 +31.746

22、0i 14.8252-42.6506i -14.2012+ 5.9172i 0 -0.62403.9002i 0 0 -0.8299+3.1120i 0 0 0 -0.8299 00 -0.6240 + 3.9002i0 +63.4921i0 0 0 0-60.4686i-14.2012+5.9172i 15.0311-8.5292i0+ 3.1120i 1.5846 - 5.5035i -0.7547-0.7547 + 2.6415i1.3787 -72.9583i0 +63.4921i 0+ 2.6415i功率方程第 (1) 次差值：Columns 1 through 1200 4.261

23、9-2.10000.1000-1.80000.5500-1.6000 8.4738 -3.7000 05.0000形成的第 (1) 次 Jacobi 矩阵：Columns 1 through 110000000000000000 00000000-37.388714.82525.9172-14.2012003.9002-0.6240000-14.8252-47.912614.20125.9172000.62403.90020005.9172-14.2012-8.029215.03113.1120-0.82990000014.20125.9172-15.0311-9.02920.82993.11

24、2000000003.1120-0.8299-5.25351.58462.6415 -0.7547 000000.82993.1120-1.5846-5.75350.7547 2.6415 0003.9002-0.624000精品资料推荐02.6415 -0.7547 -63.1845 1.3787 63.49210 0 0.6240 3.9002 00.7547 2.6415 -1.3787 -82.7321 00000000000-2.400000000000076.1905 0Columns 12 through 13000004.26190-2.100000.10000-1.80000

25、-0.55000-1.600008.473863.4921 -3.700000-63.4921 5.0000Jacobi 矩阵第 (1) 次回代运算Columns 1 through 120 -0.1490 0.1103-0.00240.16570.1015 0.3342 -0.1435 0.02970 -0.0431各个节点电压模0 1.1543 1.01600.9586 1.14391.2008功率方程第 (2) 次差值：Columns 1 through 1200-0.9196-0.3298-0.28470.5652 0.0384-1.51850.1960 -0.0019-0.3928形

26、成的第 (2) 次 Jacobi 矩阵：Columns 1 through 11-0.01270 0 0 00 0 0 0 00 0 0 0 0 0 0 0 017精品资料推荐00 -47.156723.25685.2333-16.97000 04.4126-1.1470000 -20.2170-47.588316.97005.23330 01.14704.4126000 3.5786-15.2150-5.660518.19672.9819-1.347400 000 15.21503.5786-14.7623-6.45821.34742.981900 000002.5188-1.7857-3.

27、58994.77952.1211-1.5609000001.78572.5188-1.7466-5.24071.56092.12110004.4413 -0.8294002.9982-0.9415-83.48957.152972.6032000.82944.4413000.94152.9982-0.3333 -83.28441.8854000 0000 000 -2.4000000 0000 0-2.737576.1905-1.8854Columns 12 through 1300000-0.91960-0.32980-0.28470-0.01270-0.565200.0384-1.8854

28、-1.518572.6032 0.1960-0.0862 -0.0019-72.6032 -0.3928Jacobi 矩阵第 (2) 次回代运算Columns 1 through 120 0.04410.0155 0.0700 -0.02100190.1790 -0.0275 0.0288 0.0422 -0.001 0.0487各个节点电压模精品资料推荐0 1.1120 0.9435 0.78211.1170 1.2010功率方程第 (3) 次差值：Columns 1 through 1200 -0.04170.0017-0.0185-0.1340 -0.0217-0.0367 0.0646

29、 -0.0024-0.0916-0.012647形成的第 (3) 次 Jacobi 矩阵：Columns 1 through 11000 00000000000 00 00000000-44.191723.52304.7526-16.4350004.2309-1.1798000-19.9621-46.331516.43504.7526001.17984.23090003.4621-14.0963-5.087217.05792.7814-1.22400000014.09633.4621-13.4379-6.46741.22402.781400000001.9845-1.5515-1.89914.

30、35041.6691-1.3531 000001.55151.9845-1.3054-5.04841.35311.6691 0004.3026-0.9762002.8902-1.031379.882410.0760 70.7751000.97624.3026001.03132.89020-3.4952 -82.57374.56770000 00000-2.40190000 00000.353576.2524 -4.5677Columns 12 through 130 00 00 -0.04170 0.00170 -0.01850 -0.01260 -0.13400 -0.0217 -4.567

31、7 -0.0367 70.7751 0.0646 0.0111 -0.0024 -70.7751 -0.0916Jacobi 矩阵第 (3) 次回代运算Columns 1 through 12000.00780.00430.0191-0.00630.0768 -0.00550.0056 0.01560.00110.0181各个节点电压模0 1.10480.9237 0.70981.11261.2001功率方程第 (4) 次差值：Columns 1 through 1200-0.00040.0002-0.0002-0.0012-0.0276-0.00580.0006 0.0016 -0.0003

32、-0.0053形成的第 (4) 次 Jacobi矩阵：Columns 1 through 1100000 00 000000000 00 000000-43.744923.59274.6455-16.35040 04.1980-1.1918000-20.0302-45.989216.35044.64550 01.19184.19800003.4389-13.7881-4.989916.76802.7273-178813.4389-13.0467-6.42921.18852.727300000001.7502-1.4706-1.12174.24801.4705 -1.2

33、806 000001.47061.7502-1.1039-4.99841.2806 1.4705 0004.2712 -1.0337002.8637-1.0684-79.3737 11.1438 70.4222001.03374.2712001.0684 2.8637 -4.6945 -82.2282 5.56080 0 0 0 0 0 0 0 0 0 -2.39980 0 0 0 0 00 1.5031 76.1844 -5.5608Columns 12 through 130 00 00 -0.00040 0.00020 -0.00020 -0.00120 -0.02760 -0.0058

34、 -5.5608 0.0006 70.4222 0.0016 0.0473 -0.0003-70.4222 -0.0053Jacobi 矩阵第 (4) 次回代运算Columns 1 through 1200.0259 -0.00180.00160 0.00210.0048 0.00020.00140.00530.0056-0.0020各个节点电压模0 1.10290.91780.6857 1.11141.2000功率方程第 (5) 次差值：Columns 1 through 120 0 -0.00000.00000.0000-0.0001-0.0032 -0.0007 -0.0000 0.00

35、01 -0.0000-0.0005形成的第 (5) 次 Jacobi 矩阵：Columns 1 through 1100000 0000000000000000000-43.630023.62104.6137 -16.3286004.1889-1.1958000 -20.0561-45.884116.32864.61370 01.19584.1889000 3.4337-13.6963-4.968816.68032.7114-1.177600 000 13.69633.4337-12.9312-6.41371.17762.711400 000001.6711-1.4437-0.85644.21

36、621.4034-1.2564000001.44371.6711-1.0344-4.98431.25641.40340004.2621-1.0517002.8559-1.0801-79.237411.491170.3226001.05174.2621001.08012.8559-5.0505 -82.12235.8687000 0000 000 -2.3993000 0000 01.841976.1690-5.8687Columns 12 through 130 00 00 -0.00000 0.00000 0.00000 -0.00010 -0.00320 -0.0007-5.8687 -0

37、.000070.3226 0.00010.0580 -0.0000-70.3226 -0.0005Jacobi 矩阵第 (5) 次回代运算Columns 1 through 120 0.00030.0002 0.0008 -0.00030.0007 0.0000 0.00080.6820 1.1113 1.20000.0039 -0.0003 0.0002 各个节点电压模0 1.1026 0.9169功率方程第 (6) 次差值：1.0e-004 *Columns 1 through 120 0 -0.00260.00070.0032-0.0182-0.7419 -0.1647 -0.0016

38、0.0184 -0.0063-0.1036形成的第 (6) 次 Jacobi 矩阵：Columns 1 through 1100000 00000000000 00000000-43.612823.62534.6089004.1876-1.1964000-20.0599-45.868316.3254001.19644.18760003.4330-13.6825-4.96582.7090-168253.4330-12.91391.17602.709000000001.6591-0.81614.21141.3932 -1.2528 000001.4396-1.0238-4

39、.98221.2528 1.3932 0004.26071.054302.8548-1.081879.2170 11.5431 70.3078001.05434.260701.08182.8548-5.1036 -82.1067 5.914800000 00000-2.399300000 0001.892476.1670 -5.9148-16.32544.608916.6670-6.4113-1.43961.659100Columns 12 through 1300000-0.000000.000000.00000-0.00000-0.00010-0.0000-5.9148 -0.000070

40、.3078 0.00000.0596 -0.0000-70.3078 -0.0000Jacobi 矩阵第 （6） 次回代运算1.0e-004 *Columns 1 through 1200 0.07640.04960.2048-0.07280.9503 -0.06410.05640.1753 0.00740.1923各个节点电压模0 1.10260.91690.6819 1.11131.2000功率方程第 （7） 次差值：1.0e-007 *Columns 1 through 120 -0.0015 0.00040.0019 -0.0106-0.4343 -0.0965 -0.0010 0.0

41、108 -0.0037 -0.06064.2 计算机法潮流计算结果迭代次数：6 没有达到精度要求的个数：9 10 10 10 8 3 0各节点的电压复数值 E 为 （节点号从小到大排列 ）：1.2000 1.0947 - 0.1316i 0.9068 - 0.1361i 0.6128 - 0.2991i1.1073 - 0.0932i 1.1996 - 0.0298i各节点的电压模值大小 V 为（节点号从小到大排列 ）：1.2000 1.1026 0.9169 0.6819 1.1113 1.2000各节点的电压相角 sida 为（节点号从小到大排列 ）：0 -6.8549 -8.5347 -

42、26.0187 -4.8097 -1.4241各节点的功率 S 为（节点号从小到大排列 ）：5.0135 + 1.8332i -2.1000 - 1.0000i -1.8000 - 0.4000i -1.6000 - 0.8000i-3.7000 - 1.3000i 5.0000 + 2.5555i各条支路的首端功率 Si 为（顺序同您输入 B1 时一致 ）：S(1,2)=5.0135+1.8332iS(2,3)=3.0895+0.48941iS(2,5)=-0.176-0.31075iS(3,4)=0.78738+0.39432iS(4,5)=-0.90639-0.43087iS(6,5)=

43、5+2.5555i各条支路的末端功率 Sj 为(顺序同您输入 B1 时一致 )：S(2,1)=-5.0135-1.1787iS(3,2)=-2.5874-0.79432iS(5,2)=0.17703-0.29553iS(4,3)=-0.69361-0.36913iS(5,4)=1.123+1.1889iS(5,6)=-5-2.1934i各条支路的功率损耗DS为(顺序同您输入B1时一致):DS(1,2)=-8.8818e-016+0.65451iDS(2,3)=0.50212-0.3049iDS(2,5)=0.0010207-0.60628iDS(3,4)=0.093766+0.025185iDS(4,5)=0.21659+0.75805iDS(6,5