计算机潮流计算

1、学号工厂供电课程设计（2013级本科）题 目： 某工厂变电所一次系统设计 系（部）院： 物理与机电工程学院 专 业： 电气工程及其自动化 作者姓名： 杨兴海 指导教师： 田娜 职称： 助教 完成日期： 2016 年 1 月 5 日河西学院本科生课程设计任务书设 计 题 目计算机潮流计算作 者 姓 名杨兴海学院、专业、年级物电学院电气工程专业13级指导教师姓名、职称田娜 助教任务下达日期2015年6月1日一、设计内容：结合电力系统分析课程，进一步学习并理解电力系统潮流计算的物理及数学模型，利用PQ分解方法对某电力系统进行潮流计算，根据潮流计算结果进行电压调整分析。1.网络接线图及计算数据如下：2

2、.线路数据： L1、L2和L3型号均为LGJQ4400，长度分别为340km、280km和300km；L4和L5型号分别为LGJ240和LGJ185，长度分别为80km和8km。变压器数据： 自行选择变压器型号。母线2、10、6、11、8的电压等级分别为：500kV、35kV、220kV、110kV和10kV。变压器标幺值变比取1.0。节点类型： 节点1、4分别为平衡节点和PV节点，其余为PQ节点。发电机组数据： 等值发电机G2有功出力700MW，等值发电机G1有功最大出力700MW。负荷数据： 节点2：450+j100MVA；节点3：300+j120MVA； 节点5：350+j130MVA；

3、节点7：120+j80MVA； 节点9：2.5+j0.9MVA；节点10：30+j10MVA。二、任务要求：（1） 查设计手册，确定变压器、线路等的电气参数；（2） 选择线路、变压器等模型，进行有名值和标幺值参数计算，并归算形成等值网络；（3） 手算Jaccobi矩阵的部分元素；（4） 应用程序进行潮流计算（机端母线电压取为1.0p.u.）；（5） 分析不同节点电压水平的差别及其原因；（6） 选择适当调压方式进行电压调整；（7） 撰写设计报告。三、论文(设计)进度安排阶段论文（设计）各阶段名称起止日期1熟悉设计任务书、设计题目及设计背景资料6.16.52查阅有关资料6.66.73阅读设计要求必

4、读的参考资料6.86.104书写设计说明书6.116.195小组答辩质疑6.216.226上交设计成果6.30四、需收集和阅读的资料及参考文献（指导教师指定）1: 陈珩.电力系统稳态分析（第三版）M，北京，中国电力出版社，20072：何仰赞.温增银.电力系统分析第三版M，武汉,华中科技大学出版社，20023：陈悦.电气工程毕业设计指南电力系统分册M，北京，中国水利水电出版社，2008教 研 室 意 见 负责人签名： 年 月 日 目录摘要：1第一章 电力系统潮流计算概述21.1 潮流计算简介21.2 潮流计算的意义及其发展2第二章 电气参数的确定32.1变压器、线路参数确定32.2 电力线路、变

5、压器模型的建立4第三章 PQ分解方法潮流计算分析53.1 PQ分解法的极坐标表示及简化算法53.1.1潮流计算的定义53.1.2潮流计算的约束条件63.1.3节点电压用极坐标表示时的牛顿-拉夫逊潮流计算63.1.4对牛顿拉夫逊法潮流计算的数学模型进行简化修正83.2 PQ分解法潮流计算的简化算法103.3 PQ 分解法潮流计算的基本步骤13第四章 程序编写及结果分析134.1程序编程134.2执行结果20第五章 课程设计心得22参考文献22摘要潮流计算是电力系统最基本最常用的计算。根据系统给定的运行条件，网络接线及元件参数，通过潮流计算可以确定各母线的电压（幅值和相角），各支路流过的功率，整个

6、系统的功率损耗。潮流计算是实现电力系统安全经济发供电的必要手段和重要工作环节。因此，潮流计算在电力系统的规划计算，生产运行，调度管理及科学计算中都有着广泛的应用。潮流计算在数学上是多元非线性方程组的求解问题，PQ分解法是数学上解非线性方程组的有效方法，有较好的收敛性。运用电子计算机计算一般要完成以下几个步骤：建立数学模型，确定解算方法，制订计算流程，编制计算程序。关键字：PQ分解法 计算机潮流计算 MATLAB 第一章 电力系统潮流计算1.1 潮流计算简介 电力系统潮流计算是研究电力系统稳态运行情况的一种计算，它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态：各母线的电压，各

7、元件中流过的功率，系统的功率损耗等等。在电力系统规划的设计和现有电力系统运行方式的研究中，都需要利用潮流计算来定量地分析比较供电方案或运行方式的合理性。可靠性和经济性。此外，电力系统潮流计算也是计算系统动态稳定和静态稳定的基础。所以潮流计算是研究电力系统的一种很重要和基础的计算。 电力系统潮流计算也分为离线计算和在线计算两种，前者主要用于系统规划设计和安排系统的运行方式，后者则用于正在运行系统的经常监视及实时控制。 利用电子数字计算机进行电力系统潮流计算从50年代中期就已经开始。在这20年内，潮流计算曾采用了各种不同的方法，这些方法的发展主要围绕着对潮流计算的一些基本要求进行的。对潮流计算的要

8、求可以归纳为下面几点： （1）计算方法的可靠性或收敛性； （2）对计算机内存量的要求； （3）计算速度； （4）计算的方便性和灵活性。 电力系统潮流计算问题在数学上是一组多元非线性方程式求解问题，其解法都离不开迭代。因此，对潮流计算方法，首先要求它能可靠地收敛，并给出正确答案。由于电力系统结构及参数的一些特点，并且随着电力系统不断扩大，潮流计算的方程式阶数也越来越高，对这样的方程式并不是任何数学方法都能保证给出正确答案的。这种情况成为促使电力系统计算人员不断寻求新的更可靠方法的重要因素。 1.2 潮流计算的意义及其发展 (1)在电网规划阶段,通过潮流计算,合理规划电源容量及接入点,合理规划 网

9、架,选择无功补偿方案,满足规划水平的大、小方式下潮流交换控制、调峰、调相、调压的要求。 (2)在编制年运行方式时,在预计负荷增长及新设备投运基础上,选择典型方式进行潮流计算,发现电网中薄弱环节,供调度员日常调度控制参考,并对规划、基建部门提出改进网架结构,加快基建进度的建议。 (3)正常检修及特殊运行方式下的潮流计算,用于日运行方式的编制,指导发电厂开机方式,有功、无功调整方案及负荷调整方案,满足线路、变压器热稳定要求及电压质量要求。 (4)预想事故、设备退出运行对静态安全的影响分析及作出预想的运行方式调整方案。 总结为在电力系统运行方式和规划方案的研究中，都需要进行潮流计算以比较运行方式或规

10、划供电方案的可行性、可靠性和经济性。同时，为了实时监控电力系统的运行状态，也需要进行大量而快速的潮流计算。因此，潮流计算是电力系统中应用最广泛、最基本和最重要的一种电气运算。在系统规划设计和安排系统的运行方式时，采用离线潮流计算；在电力系统运行状态的实时监控中则采用在线潮流计算。 近20多年来，潮流算法的研究仍然非常活跃，但是大多数研究都是围绕改进牛顿法和P-Q分解法进行的。此外，随着人工智能理论的发展遗传算法、人工神经网络、模糊算法也逐渐被引入潮流计算。但是，到目前为止这些新的模型和算法还不能取代牛顿法和P-Q分解法的地位。由于电力系统规模的不断扩大，对计算速度的要求不断提高，计算机的并行计

11、算技术也将在潮流计算中得到广泛的 应用，成为重要的研究领域。第二章 电气参数的确定2.1变压器、线路参数确定查设计手册：LGJQ4400型导线:=0.020/km,=0.276/km,=4.02S/km; LGJ240型导线:=0.131/km,=0.432/km; LGJ185型导线:=0.170/km,=0.440/km。变压器型号:T1,T2为SF-/220；T3部分参数为:额定容量240/120/120MVA,额定电压500/242/38.5kv。T4部分参数为:额定容量120/120/60MVA, ,额定电压242/121/10.5kv,等值电路中所有参数都归算到高压侧。线路L1,L

12、2的电导,L4,L5的导纳都可略去。表2-1短路电压（未经归算）绕组高压-中压高压-低压中压-低压短路电压百分值(%)1030202.1.1电力变压器及线路参数的计算1.采用有名值： 2.采用标幺值：取基准功率为1000MVA，基准电压为500kv 2.2.2电力变压器及线路参数的计算结果 电力线路参数见表2-1，电力变压器参数见表2-2表2-1电力线路参数 线 路参 数LGJQ4400LGJ240L4LGJ185L5L1L2L36.85.66.010.481.360.02720.02240.02400.04190.005493.8477.2882.834.563.520.37540.3091

13、0.33120.13820.01416.8+j93.845.6+j77.286.0+j82.810.48+j34.561.36+j3.520.0272+j0.37540.0224+j0.30910.024+j0.33120.0419+j0.13820.0054+j0.0141S1.36681.12561.206/ 3.4172.8143.015/表2-2电力变压器参数 变压器参数T1T2T3T425.5425.540.2430.1420.4860.1420.4860.284 0.10220.10221.921.92104.1748.800208.3397.60.00768 0.007680.4

14、1670.1952000.83330.3904 25.54+j1.9225.54+j1.920.243+j104.170.142+j48.80.4860.1420.486+j208.330.284+j97.60.1022+j0.007680.1022+j0.00768+j0.4167+j0.1952+j0.8333+j0.3904 S17.7717.779.565.460.004440.004440.002390. S20.8320.8316.3912.290.0.0.0. 17.77+j20.8317.77+j20.839.56+j16.395.46+12.294.44+j5.2084.44

15、+j5.2082.39+j4.0981.365+j3.0732.2 电力线路、变压器模型的建立线路采用型等值电路，变压器采用型等值电路，有名值和标幺值等值电路分别见：图2-1、图2-2。图2-1 有名值等值网络图2-2标幺值等值网络第三章 PQ分解方法潮流计算分析3.1 PQ分解法的极坐标表示及简化算法3.1.1潮流计算的定义潮流计算是根据给定的电网结构、参数和发电机、负荷等元件的运行条件，确定电力系统各部分稳态运行状态参数的计算。通常给定的运行条件有系统中各电源和负荷点的功率、枢纽点电压、平衡点的电压和相位角。待求的运行状态参量包括电网各母线节点的电压幅值和相角，以及各支路的功率分布、网络的

16、功率损耗等。PQ分解法的极坐标表示是派生于以极坐标表示的牛顿拉夫逊法。3.1.2潮流计算的约束条件为了保证电力系统的正常运行潮流问题中某些变量应满足一定的约束条件，常用的约束条件有：(1)所有节点电压必须满足 (i=1,2，n) 从保证电能质量和供电安全的要求看，电力系统的所有电气设备都必须运行在额定电压附近。PV节点的电压幅值必须按上述条件给定。因此，这一约束主要是对PQ节点而言。(2)所有电源节点的有功功率和无功功率必须满足 PQ节点的有功功率和无功功率以及PV节点的有功功率，在给定时就必须满足上述条件。因此，对平衡节点的P和Q以及PV节点的Q应按上述条件进行检验。(1)某些节点之间电压的

17、相位差应满足 为了保证系统运行的稳定性，要求某些输电线路两端的电压相位差不超过一定的数值。因此，潮流计算可以归纳为求解一组非线性方程组，并使其解答满足一定的约束条件。如果不满足，则应修改某些变量的给定值，甚至修改系统的运行方式，重新进行计算。3.1.3节点电压用极坐标表示时的牛顿-拉夫逊潮流计算采用极坐标时，节点电压表示为节点功率方程（11-25）将写成 公式 1式中，是i，j两节点电压的相差角。方程式（公式1）把节点功率表示为节点电压的幅值和相角的函数。在有n节点的系统中，假定第1m号节点PQ节点，第m+1n-1号节点的PV节点，第n号节点为平衡节点。和是给定的，PV节点的电压幅值Vm+1V

18、n-1也是给定的。因此，只剩下n-1个节点的电压相角和m个节点的电压幅值,是未知量。实际上，对于每一个PQ节点或每一个PV节点都可以列写一个有功功率不平衡量方程式=-=- （i=1,2，n-1） 公式 2而对于每一个PQ节点还可以再列写一个无功功率不平衡量方程式=-=- （i=1,2，m） 公式 3式2和式3一共包含了n-1+m个方程式，正好同未知量的数目相等，而比直角坐标形式的方程式少了n-1-m个。对于方程式2和式3可以写出修正方程式如下： 公式 4式中； ； ； 公式 5H是(n-1)(n-1)阶方阵，其元素为；N是(n-1)m阶矩阵，其元素为；K是m (n-1)阶矩阵，其元素为；L是m

19、 m阶方阵，其元素为。3.1.4对牛顿拉夫逊法潮流计算的数学模型进行简化修正在交流高压电网中，输电线路的电抗要比电阻大得多，系统中母线有功功率的变化主要受电压相位的影响，无功功率的变化则主要受母线电压幅值变化的影响。在修正方程式的系数矩阵中，偏导数和的数值相对于偏导数和是相当小的。作为简化的第一步，可以将上式子块N和K略去不计，即认为它们的元素都等于零。这样，n-1+m阶的方程式便分解为一个n-1阶和一个m阶的方程。这一简化大大地节省了机器内存和解题时间。以上方程式表明，节点的有功功率不平衡量只用于修正电压的相位，节点的无功功率不平衡量只用于修正电压的幅值。这两组方程分别轮流进行迭代，这就是所

20、谓有功-无功功率分解法。但矩阵H和L的元素都是节点电压幅值和相角差的函数，其数值在迭代过程中是不断变化的。因此，最关键的一步简化就在于，把系数矩阵H和L简化成常数矩阵。它的根据是什么呢?在一般情况下，线路两端电压的相角差是不大的(不超过100200)，因此可以认为 此外，与系统各节点无功功率相适应的导纳必远小于该节点自导纳的 或 考虑到以上的关系，矩阵H和L的元素的表达式便被简化成 （i，j=1，2，n-1） 公式 6 （i，j=1，2，m） 公式 7而系数矩阵H和L则可以分别写成H= 公式 8 L= = = 公式 9代入得: 公式 10 公式 11这就是简化了的修正方程，也可展开为=- 公式

21、 12=- 公式 13在这两个修正方程式中，系数矩阵都由节点导纳矩阵的虚部构成，只是阶次不同，矩阵为n-1阶，不含平衡节点对应的行和列，矩阵为m阶，不含平衡节点和PV节点所对应的行和列。由于修正方程的系数矩阵为常数矩阵，只要作一次三角分解，即可反复使用，结合采用稀疏技巧，还可进一步的节省机器内存和计算时间。3.2 PQ分解法潮流计算的简化算法PQ分解法潮流计算修正方程的简化根据极坐标表示的牛顿-拉夫逊法可得：每一个PQ或PV节点的有功功率不平衡方程：=-=- （i=1,2，n-1） 公式 14每一个PQ节点的无功功率不平衡方程：=-=- （i=1,2，m） 公式 15由此可以写出修正方程如下：

22、 公式 16对式24和式25求偏导数，可以得到雅克比矩阵元素的表达如下：a) 当ij时b) 当ij时对修正方程的第一个简化是：将式26中的N、K子阵略去而将其简化为 公式 17上式可分别写成以下两式 公式 18 公式 19对修正方程式的第二个简化式：在一般情况下，线路两端电压的相角差是不大的(不超过100200)，因此可以认为 因此可得： （i，j=1，2，n-1） （i，j=1，2，m） 对修正方程的第三个简化是：当采用标幺值计算是，可近似地认为各节点电压的大小，且当忽略接地支路时，导纳元素的对角元为非对角元之和，这样可得： 且此时式29也可以简化为 公式 20 这是雅克比矩阵的两个子阵H、

23、L具有相同的表达式，只是其阶数不同，前者为（n-1）阶，后者为m阶，两个子阵的展开式如下：将式28代入式30可得：= 公式 21= 公式 22式31和式32可简写为 公式 23式33为简化PQ分解法的修正方程式，原PQ分解法的修正方程的简化形式 通过两者的比较可知，简化PQ分解法的修正方程式比原PQ分解法的修正方程更为简单。3.3 PQ 分解法潮流计算的基本步骤（1） 形成系数矩阵 ,并求其逆矩阵。（2） 设各节点电压的初值为 (i=1,2，n，is)和(i=1,2，m， is)。 （3） 通 过 有 功 功 率 的 不 平 衡 方 程 计 算 有 功 功 率 的 不 平 衡 量 ,从 而 求

24、 出 (i=1,2，n，is)。 （4）解修正方程式，求各节点电压相位角的变量 (i=1,2，n，is)。（5） 求各节点电压相位角的新值 (i =1,2，n，is)。 （6） 通 过 无 功 功 率 的 不 平 衡 方 程 计 算 无 功 功 率 的 不 平 衡 量 ， 从 而 求 出 (i=1,2，m，is)。（7）解修正方程式，求各节点电压大小的变量 (i=1,2，m，is)。（8）求各节点电压大小的新值 (i=1,2，m，is)。 （9）运用各节点电压的新值自第三步开始进入下一次迭代。 （10）计算平衡节点功率和线路功率。第四章 程序编写及结果分析4.1程序编程n=11;nl=11;i

25、sb=1;pr=0.00001; B1=25.54+j1.92 25.54+j1.92 -j0.7316 0 -j0.0884 0 0 0 0 0 0 ;0 -j0.7316 32.3466+j94.8853 25.54+j1.92 -j0.8292 0.0266-j0.0445 0 0 0 0.0266-j0.0445 0; 0 -j0.0884 -j0.8292 -j0.9176 0 0 0 0 0 0 0 ;0 0 0.729+j104.17 0 0 12.181+j347.06 10.48+j34.56 0 0 0.972+j208.33 0;0 0 0 0 0 10.48+j34.5

26、6 11.19+j229.76 0.426+j146.4 0 0 0.284+j48.8;0 0 0 0 0 0 0.426+j146.4 2.212+j247.52 1.36+j3.52 0 0.426+j97.6;0 0 0 0 0 0 0 1.36+j3.52 1.36+j3.52 0 0;0 0 0.729+j312.5 0 0 0.972+j208.33 0 0 0 1.701+j520.83 0;0 0 0 0 0 0 0.284+j48.8 0.426+j97.6 0 0 0.71+j146.4 ;B2=25.54+j1.92 25.54+j1.92 -j0.7316 0 -j0

27、.0884 0 0 0 0 0 0 ;0 -j0.7316 32.3466+j94.8853 25.54+j1.92 -j0.8292 0.0266-j0.0445 0 0 0 0.0266-j0.0445 0; 0 0 0.0237-j0.0278 0.0237-j0.0278 0 0 0 0 0 0 0 ;0 -j0.0884 -j0.8292 -j0.9176 0 0 0 0 0 0 0 ;0 0 0.729+j104.17 0 0 12.181+j347.06 10.48+j34.56 0 0 0.972+j208.33 0;0 0 0 0 0 10.48+j34.56 11.19+j

28、229.76 0.426+j146.4 0 0 0.284+j48.8;0 0 0 0 0 0 0.426+j146.4 2.212+j247.52 1.36+j3.52 0 0.426+j97.6;0 0 0 0 0 0 0 1.36+j3.52 1.36+j3.52 0 0;0 0 0.729+j312.5 0 0 0.972+j208.33 0 0 0 1.701+j520.83 0;0 0 0 0 0 0 0.284+j48.8 0.426+j97.6 0 0 0.71+j146.4X=1 0;2 0;3 0;4 0;5 0;6 0;7 0;8 0;9 0;10 0;11 0;na=3

29、;Y=zeros(n);YI=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);O=zeros(1,n);for i=1:n if X(i,2)=0; p=X(i,1); Y(p,p)=1./X(i,2); endendfor i=1:nl if B1(i,6)=0 p=B1(i,1);q=B1(i,2); else p=B1(i,2);q=B1(i,1); end Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5); YI(p,q)=YI(p,q)-1./B1(i,3); Y(q,p)=Y(p,q); YI(q,p)=YI(p,q

30、); Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)2)+B1(i,4)./2; YI(q,q)=YI(q,q)+1./B1(i,3); Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2; YI(p,p)=YI(p,p)+1./B1(i,3);endG=real(Y);B=imag(YI);BI=imag(Y);for i=1:n S(i)=B2(i,1)-B2(i,2); BI(i,i)=BI(i,i)+B2(i,5);endP=real(S);Q=imag(S);for i=1:n e(i)=real(B2(i,3); f(i)=imag(B2(i

31、,3); V(i)=B2(i,4);endfor i=1:n if B2(i,6)=2 V(i)=sqrt(e(i)2+f(i)2); O(i)=atan(f(i)./e(i); endendfor i=2:n if i=n B(i,i)=1./B(i,i); else IC1=i+1; for j1=IC1:n B(i,j1)=B(i,j1)./B(i,i); end B(i,i)=1./B(i,i); for k=i+1:n for j1=i+1:n B(k,j1)=B(k,j1)-B(k,i)*B(i,j1); end end endendp=0;q=0;for i=1:n if B2(

32、i,6)=2 p=p+1;k=0; for j1=1:n if B2(j1,6)=2 k=k+1; A(p,k)=BI(i,j1); end end endendfor i=1:na if i=na A(i,i)=1./A(i,i); else k=i+1; for j1=k:na A(i,j1)=A(i,j1)./A(i,i); end A(i,i)=1./A(i,i); for k=i+1:na for j1=i+1:na A(k,j1)=A(k,j1)-A(k,i)*A(i,j1); end end endendICT2=1;ICT1=0;kp=1;kq=1;K=1;DET=0;ICT3

33、=1;while ICT2=0|ICT3=0ICT2=0;ICT3=0;for i=1:n if i=isb C(i)=0; for k=1:n C(i)=C(i)+V(k)*(G(i,k)*cos(O(i)-O(k)+BI(i,k)*sin(O(i)-O(k); end DP1(i)=P(i)-V(i)*C(i); DP(i)=DP1(i)./V(i); DET=abs(DP1(i); if DET=pr ICT2=ICT2+1; end endendNp(K)=ICT2; if ICT2=0 for i=2:n DP(i)=B(i,i)*DP(i); if i=n IC1=i+1; for

34、 k=IC1:n DP(k)=DP(k)-B(k,i)*DP(i); end else for LZ=3:i L=i+3-LZ; IC4=L-1; for MZ=2:IC4 I=IC4+2-MZ; DP(I)=DP(I)-B(I,L)*DP(L); end end end end for i=2:n O(i)=O(i)-DP(i); end kq=1;L=0; for i=1:n if B2(i,6)=2 C(i)=0;L=L+1; for k=1:n C(i)=C(i)+V(k)*(G(i,k)*sin(O(i)-O(k)-BI(i,k)*cos(O(i)-O(k); end DQ1(i)=

35、Q(i)-V(i)*C(i); DQ(L)=DQ1(i)./V(i); DET=abs(DQ1(i); if DET=pr ICT3=ICT3+1; end end end else kp=0; if kq=0; L=0; for i=1:n if B2(i,6)=2 C(i)=0;L=L+1; for k=1:n C(i)=C(i)+V(k)*(G(i,k)*sin(O(i)-O(k)-BI(i,k)*cos(O(i)-O(k); end DQ1(i)=Q(i)-V(i)*C(i); DQ(L)=DQ1(i)./V(i); DET=abs(DQ1(i); end end end end Nq

36、(K)=ICT3; if ICT3=0 L=0; for i=1:na DQ(i)=A(i,i)*DQ(i); if i=na for LZ=2:i L=i+2-LZ; IC4=L-1; for MZ=1:IC4 I=IC4+1-MZ; DQ(I)=DQ(I)-A(I,L)*DQ(L); end end else IC1=i+1; for k=IC1:na DQ(k)=DQ(k)-A(k,i)*DQ(i); end end end L=0; for i=1:n if B2(i,6)=2 L=L+1; V(i)=V(i)-DQ(L); end end kp=1; K=K+1; else kq=0

37、; if kp=0 K=K+1; end end for i=1:n Dy(K-1,i)=V(i); endenddisp(迭代次数)disp(K);disp(每次没有达到精度要求的有功功率个数为);disp(Np);disp(每次没有达到精度要求的无功功率个数为);disp(Nq);for k=1:n E(k)=V(k)*cos(O(k)+V(k)*sin(O(k)*j; O(k)=O(k)*180./pi;enddisp(各节点的电压标么值E为);disp(E);disp(各节点的电压V大小);disp(V);disp(各节点的电压相角O);disp(O);for p=1:n C(p)=0

38、; for q=1:n C(p)=C(p)+conj(Y(p,q)*conj(E(q); end S(p)=E(p)*C(p);enddisp(各节点的功率为);disp(S);disp(各条支路的首端功率为);for i=1:nl if B1(i,6)=0 p=B1(i,1);q=B1(i,2); else p=B1(i,2);q=B1(i,1); end Si(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); disp(Si(p,q);enddisp(各条

39、支路的末端功率为);for i=1:nl if B1(i,6)=0 p=B1(i,1);q=B1(i,2); else p=B1(i,2);q=B1(i,1); end Sj(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); disp(Sj(q,p);enddisp(各条支路的功率损耗为);for i=1:nl if B1(i,6)=0 p=B1(i,1);q=B1(i,2); else p=B1(i,2);q=B1(i,1); end DS(i)=Si(p,q)+Sj(q,p); disp(DS(i);endfor i