电力系统分析报告

1、山东交通学院电力系统分析课程设计报告书院（部）别信息科学与电气工程学院班 级学 号姓 名指导教师时 间 2014069-20130613课程设计任务书题目电力系统分析课程设计学院信息科学与电气工程学院专业电气工程及其自动化班级学生姓名学号月日至月厘日共周指导教师（签字）院长（签字）2014年6月 13日一、设计内容及要求复杂网络牛顿一拉夫逊法潮流分析与计算的设计电力系统潮流计算是电力系统中一项最基本的计算，设计内容为复杂网络潮 流计算的计算机算法 牛顿-拉夫逊法。首先，根据给定的电力系统简图，通过手算完成计算机算法的两次迭代过程， 从而加深对牛顿-拉夫逊法的理解，有助于计算机编程的应用。其次，

2、利用计算机编程对电力系统稳态运行的各参数进行解析和计算；编程 完成复杂网络的节点导纳矩阵的形成；电力系统支路改变、节点增减的程序变化; 编程完成各元件的功率损耗、各段网络的电压损耗、各点电压、功率大小和方向 的计算。二、设计原始资料给出一个46节点、环网、两电源和多引出的电力系统；参数给定，可以选用直角坐标表示的牛拉公式计算，也可以选用极坐标表示的牛 拉公式计算。具体题目详见附录题单三、设计完成后提交的文件和图表计算说明书部分设计报告和手算潮流的步骤及结果图纸部分：电气接线图及等值电路；潮流计算的计算机算法，即程序；运算结果等以图片的形式附在设计报告中。四、进程安排第一天上午：选题，查资料，制

3、定设计方案；第一天下午一一第三天下午：手算完成潮流计算的要求；第四天上午一一第五天上午：编程完成潮流计算，并对照手算结果，分析误差 第五天下午：答辩，交设计报告。五、主要参考资料电力系统分析（第三版）于永源主编，中国电力出版社，2007年电力系统分析，何仰赞温增银编著，华中科技大学出版社，2002年版;电力系统分析，韩桢祥主编，浙江大学出版社，2001年版；电力系统稳态分析，陈珩编，水利电力出版社；成绩评定表指导教师成绩答辩小组成绩总评成绩潮流计算，电力学名词，指在给定电力系统网络拓扑、元件参数和发电、负 荷参量条件下，计算有功功率、无功功率及电压在电力网中的分布。潮流计算是 根据给定的电网结

4、构、参数和发电机、负荷等元件的运行条件，确定电力系统各 部分稳态运行状态参数的计算。通常给定的运行条件有系统中各电源和负荷点的 功率、枢纽点电压、平衡点的电压和相位角。待求的运行状态参量包括电网各母 线节点的电压幅值和相角，以及各支路的功率分布、网络的功率损耗等。传统的 潮流计算程序缺乏图形用户界面，结果显示不直观，难于与其他分析功能集成， 网络原始数据输入工作量大且易于出错。结合电力系统的特点，对于复杂电力系 统，根据定条件，应用牛顿-拉夫逊法进行计算，在手工计算中，由于涉及大量 变量、微分方程、矩阵计算，求解很繁琐，计算不同系统时需要重新计算。运用 MATLAB软件进行仿真潮流计算，图形界

5、面直观，运行稳定，计算准确，提高 了计算速度，各个类的有效封装又使程序具有很好的模块性.可维护性和可重用 性。关键字：潮流计算牛拉法Matlab目录 TOC o 1-5 h z HYPERLINK l bookmark24 o Current Document 1潮流计算1 HYPERLINK l bookmark27 o Current Document 1.1潮流计算概述1 HYPERLINK l bookmark34 o Current Document 1.2潮流计算的要求2 HYPERLINK l bookmark40 o Current Document 1.3潮流计算的优势2 H

6、YPERLINK l bookmark43 o Current Document 1.4潮流计算的用途3 HYPERLINK l bookmark46 o Current Document 2 MATLAB简介4 HYPERLINK l bookmark49 o Current Document 2.1MATLAB 概述43牛顿-拉夫逊法概述53.1牛顿-拉夫逊基本原理5 HYPERLINK l bookmark59 o Current Document 3.2直角牛顿-拉夫逊法潮流计算求解过程6 HYPERLINK l bookmark205 o Current Document 4设计过程

7、124.1 题目 D124.2牛顿拉夫逊法程序流程图2114.3设计程序22 HYPERLINK l bookmark242 o Current Document 心得体会25 HYPERLINK l bookmark245 o Current Document 参考文献261潮流计算1.1潮流计算概述潮流计算是研究电力系统稳态运行情况的一种基本电气计算，常规潮流计算 的任务是根据给定的运行条件和网路结构确定整个系统的运行状态，如各母线上 的电压(幅值及相角)、网络中的功率分布以及功率损耗等。潮流计算的结果是 电力系统稳定计算和故障分析的基础。通过潮流计算可以判断电网母线电压、支 路电流和功率

8、是否越限，如果有越限，就应采取措施，调整运行方式。对于正在 规划的电力系统，通过潮流计算，可以为选择电网供电方案和电气设备提供依据。 潮流计算还可以为继电保护和自动装置定整计算、电力系统故障计算和稳定计算 等提供原始数据。具体表现在以下方面：在电网规划阶段，通过潮流计算，合理规划电源容量及接入点，合理规划 网架，选择无功补偿方案，满足规划水平的大、小方式下潮流交换控制、调峰、调 相、调压的要求。在编制年运行方式时，在预计负荷增长及新设备投运基础上，选择典型方 式进行潮流计算，发现电网中薄弱环节,供调度员日常调度控制参考，并对规划、 基建部门提出改进网架结构，加快基建进度的建议。正常检修及特殊运

9、行方式下的潮流计算，用于日运行方式的编制，指导发 电厂开机方式，有功、无功调整方案及负荷调整方案，满足线路、变压器热稳定要 求及电压质量要求。预想事故、设备退出运行对静态安全的影响分析及作出预想的运行方式 调整方案。总结为在电力系统运行方式和规划方案的研究中，都需要进行潮流计算以比 较运行方式或规划供电方案的可行性、可靠性和经济性。同时，为了实时监控电 力系统的运行状态，也需要进行大量而快速的潮流计算。因此，潮流计算是电力 系统中应用最广泛、最基本和最重要的一种电气运算。在系统规划设计和安排系 统的运行方式时，采用离线潮流计算；在电力系统运行状态的实时监控中，则采 用在线潮流计算。此外，电力系

10、统潮流计算也是计算系统动态稳定和静态稳定的基础。所以潮流计算是研究电力系统的一种很重要和基础的计算。1.2潮流计算的要求电力系统运行必须满足一定技术和经济上的要求。这些要求够成了潮流问题 中某些变量的约束条件，常用的约束条件如下：节点电压应满足U U .i minGiGi maxPQ节点的有功功率和无功功率，以及PU节点的有功功率，在给定是就必须 满足上述条件，因此，对平衡节点的P和Q以及PU节点的Q应按上述条件进行 检验。节点之间电压的相位差应满足I。1=1。-0 .Ivl。-0.1建茴为了保证系统运行的稳定性，要求某些输电线路两端的电压相位不超过一定 的数值。这一约束的主要意义就在于此。因

11、此，潮流计算可以归结为求解一组非线性方程组，并使其解答满足一定 的约束条件。常用的方法是迭代法和牛顿法，在计算过程中，或得出结果之后用 约束条件进行检验。如果不能满足要求，则应修改某些变量的给定值，甚至修改 系统的运行方式，重新进行计算。1.3潮流计算的优势电力系统潮流计算是电力系统分析中的一种最基本的计算，是对复杂电力系 统正常和故障条件下稳态运行状态的计算。潮流计算的目标是求取电力系统在给 定运行状态的计算。即节点电压和功率分布，用以检查系统各元件是否过负荷。 各点电压是否满足要求，功率的分布和分配是否合理以及功率损耗等。对现有电 力系统的运行和扩建，对新的电力系统进行规划设计以及对电力系

12、统进行静态和 暂态稳定分析都是以潮流计算为基础。潮流计算结果可用于电力系统稳态研究， 安全估计或最优潮流等，实际电力系统的潮流计算主要采用牛顿-拉夫逊法。借 由MATLAB可以轻松实现计算复杂的电力系统潮流分布。1.4潮流计算的用途流计算是电力系统非常重要的分析计算，用以研究系统规划和运行中提出的 各种问题。对规划中的电力系统，通过潮流计算可以检验所提出的电力系统规划 方案能否满足各种运行方式的要求；对运行中的电力系统，通过潮流计算可以预 知各种负荷变化和网络结构的改变会不会危及系统的安全，系统中所有母线的电 压是否在允许的范围以内，系统中各种元件（线路、变压器等）是否会出现过负荷， 以及可能

13、出现过负荷时应事先采取哪些预防措施等。潮流计算是电力系统分析最基本的计算。除它自身的重要作用之外， 在电力系统分析综合程序（PSASP）中，潮流计算还是网损计算、静态安 全分析、暂态稳定计算、小干扰静态稳定计算、短路计算、静态和动态等 值计算的基础2MATLAB 简介2.1MATLAB 概述MATLAB的含义是矩阵实验室(Matrix Laboratory)，是美国MathWork公司 于1982推出的一套高性能的数值计算可视化软件，包括MATLAB主程序、 SIMULINK动态系统仿真包和各种专业工具箱它集数值分析，矩阵计算，信号处 理和图形显示于一体，构成一个方便的，界面友好的用户环境，具

14、有极强大的计 算功能和极高的编程效率，特别适合于科学计算、数值分析、系统仿真和信号处 理等任务。MATLAB程学设计语言结构完整，且具有优良的移植性，它的基本数 据元素是不需要定义的数组。它可以高效率的解决工业计算问题，特别是关于矩 阵和矢量的计算。通过MATLAB语言，可以用类似数学公式的方式来编写算法， 大大降低了程序需要的难度别难并节省了时间，从而可把主要的经历集中在算法 的构思而不是编程上。学习运用MATLAB计算电力系统潮流分布是本次课程设计 的重点，可以说，作为工科学生，会运用MATLAB来解决工程问题已成为必须。 到目前为止，MATLAB 已发展成为国际上最优秀的科技应用软件之一

15、。它的功能 十分强大，不仅仅可以实现计算潮流分布，还可以模拟仿真各式各样的数值系统， 工程。这里将借助MATLAB来完成用直角牛顿-拉夫逊法计算电力系统潮流分布。3牛顿-拉夫逊法概述3.1牛顿-拉夫逊基本原理潮流计算的目标是求取电力系统在给定运行状态的计算。即节点电压和功率 分布，用以检查系统各元件是否过负荷。各点电压是否满足要求，功率的分布和 分配是否合理以及功率损耗等。对现有电力系统的运行和扩建，对新的电力系统 进行规划设计以及对电力系统进行静态和暂态稳定分析都是以潮流计算为基础。 潮流计算结果可用如电力系统稳态研究，安全估计或最优潮流等对潮流计算的模 型和方法有直接影响。实际电力系统的潮

16、流技术那主要采用牛顿-拉夫逊法。牛顿一拉夫逊法(简称牛顿法)在数学上是求解非线性代数方程式的有效方 法。其要点是把非线性方程式的求解过程变成反复地对相应的线性方程式进行求 解的过程。即通常所称的逐次线性化过程。对于非线性代数方程组：f (x) = 0艮口 fi(x ,x ,.,x ) = 0(i = 1,2,n)(3-1-1)在待求量x的某一个初始估计值x(0)附近，将上式展开成泰勒级数并略去二 阶及以上的高阶项，得到如下的经线性化的方程组：f (x(0) + f(x(0) Ax(0)= 0(3-1-2)上式称之为牛顿法的修正方程式。由此可以求得第一次迭代的修正量Ax (0) =- f(x (

17、0) )-i f (x (0)(3-1-3)将Ax(0)和x(0)相加，得到变量的第一次改进值x(i)。接着就从x(i)出发，重复上述计算过程。因此从一定的初值x(0)出发，应用牛顿法求解的迭代格式为：f( x (k) Ax (k) = 一 f ( x (k)(3-1-4)x (k+1) = x (k) + Ax (k)(3-1-5)上两式中：f(x)是函数f (x)对于变量x的一阶偏导数矩阵，即雅可比矩阵J;k为迭代次数。有上式可见，牛顿法的核心便是反复形式并求解修正方程式。牛顿法当初始 估计值双0）和方程的精确解足够接近时，收敛速度非常快，具有平方收敛特性。牛顿潮流算法突出的优点是收敛速度

18、快，若选择到一个较好的初值，算法将 具有平方收敛特性，一般迭代45次便可以收敛到一个非常精确的解。而且其迭 代次数与所计算网络的规模基本无关。牛顿法也具有良好的收敛可靠性，对于对 以节点导纳矩阵为基础的高斯法呈病态的系统，牛顿法也能可靠收敛。牛顿法所 需的内存量及每次迭代所需时间均较高斯法多。牛顿法的可靠收敛取决于有一个良好的启动初值。如果初值选择不当，算法 有可能根本不收敛或收敛到一个无法运行的节点上。对于正常运行的系统，各节 点电压一般均在额定值附近，偏移不会太大，并且各节点间的相位角差也不大， 所以对各节点可以采用统一的电压初值（也称为平直电压），如假定：U（0）= 1 0（0）= 0

19、或 e（0）= 1 f（0）= 0 （i = q1,2, .,n;i 丰 s）（3-1-6）iiii这样一般能得到满意的结果。但若系统因无功紧张或其它原因导致电压质量 很差或有重载线路而节点间角差很大时，仍用上述初始电压就有可能出现问题。 解决这个问题的办法可以用高斯法迭代12次，以此迭代结果作为牛顿法的初 值。也可以先用直流法潮流求解一次以求得一个较好的角度初值，然后转入牛顿 法迭代。3.2直角牛顿-拉夫逊法潮流计算求解过程以下讨论的是用直角坐标形式的牛顿一拉夫逊法潮流的求解过程。当采用直 角坐标时，潮流问题的待求量为各节点电压的实部和虚部两个分量 e,f ,e ,f .e ,f由于平衡节点

20、的电压向量是给定的，因此待求两共2（n-1） 1122 n n需要2（n-1）个方程式。事实上，除了平衡节点的功率方程式在迭代过程中没有 约束作用以外，其余每个节点都可以列出两个方程式。对PQ节点来说，尸和Q是给定的，因而可以写出a P = P - e (G e - B f)- f (G f + Be)= 0(3-2-1)(3-2-2)i is i . .，jj j j j 一，jj j jaQ = Q - f (G e - B f)+e (G f + B e)= 0i is iij j ij j jij j ij j对PV节点来说，给定量是P和V.,因此可以列出a P = P - e (Ge

21、 - B f)- f (G f +Be)= oiij j ij ji is iij j ij jj iaV 2=V 2-(e2+f 2)= 0iisii求解过程大致可以分为以下步骤：（1）形成节点导纳矩阵Y（2）将各节点电压设初值U，（3）将节点初值代入式（2-2-1）或式（2-2-2），求出修正方程式的常数项向 量（4）将节点电压初值代入求式，求出雅可比矩阵元素（5）求解修正方程，求修正向量（6）求取节点电压的新值（7）检查是否收敛，如不收敛，则以各节点电压的新值作为初值自第3步 重新开始进行狭义次迭代，否则转入下一步（8）计算支路功率分布，PV节点无功功率和平衡节点柱入功率。以直角坐标系形

22、式表示.迭代推算式采用直角坐标时，节点电压相量及复数导纳可表示为:(3-2-3)将以上二关系式代入上式中，展开并分开实部和虚部；假定系统中的第1,2, .,号为PQ节点，第m+1,m+2, . . ,n-1为PV节点，根据节点性质的 不同，得到如下迭代推算式：对于PQ节点AP= P - e (G e - B f ) - f (G f + B e ) TOC o 1-5 h z i i i ij j ij j i j j j jj=1j=1 HYPERLINK l bookmark76 o Current Document AQ = Q - f (G e - B f ) + e (G f + B

23、 e ) i i i ij j ij j i ij j ij j j=1j=1i = 1,2，m对于PV节点AP = P - e (G e - B f ) - f (G f + B e )i i i ij j ij j i ij j ij j j=1j=1AV 2 = V 2 -(e 2 + f 2)Ii i i(3-2-4)(3-2-5)i = m +1, m + 2,，n 1对于平衡节点平衡节点只设一个，电压为已知，不参见迭代，其电压为:匕=匕+ jfn(3-2-6),修正方程式迭代式共包括2(n-1)个方程.选定电压初值及变量修正量符号，代入方程并按泰勒级数展开，略去Ae,驾二次方程及以

24、后各项，得到一组线性方程组或线性化了的方程组，常称修正方程组:AW = - J AU(3-2-7)AW =AP1AQ1Ae1Af1APmAQmAPm+1AU 2m+1AU =AemAmAem+1Afm+1APn-1AU 2n-1 JAen-1An-1 二SAPSAPrrSeSfSAQSAQ HYPERLINK l bookmark127 o Current Document r rSeSf.1.1SAPSAPmSe1SAQmSe1SPm+1Se1SAU 2mrrSe.1sAQmSf1SPmrrSfrSAU 2mrrSf.1SAPnrSe1SAU 2 nrSe1SAP nr fSAU 2nr气S

25、APrSemSAQrSe.mSAPmSemSAQmSemSPm十1SemSAU 2m十iSe.mSAPnTSemSAU 2nrSemSAPrSfmSAQrSf.mSAPmSAQmSfmSPmrrSfmSAU 2m十rSf.mSAPnrSfmSAU 2nrSfmSAPSAPr rSeSfm+rm+rSAQSAQr rSeSfm+rm+rSAPSAPmmSeSfm+1m+1SAQSAQmmSeSfm+1m+1SPSPm+rm+rSeSfm+1m+1SAU 2SAU 2 HYPERLINK l bookmark143 o Current Document m+rm+rSeSfm+rm+rSAPSAP

26、nrnrSeSfm+1m+1SAU 2SAU 2h-1n-rSe Sfm+1m+1SAPSAPr rSeSfu-1u-1SAQSAQr rSeSfn-1n-1SAPSAPmmSeSfn-1n-1SAQSAQmmSeSfn-1n-1SPSPm+1m+1SeSfn-1n-1SAU 2SAU 2m+rm+rSe Sfn-1n-1SAPSAPnrnrSeSfn-1n-1SAU 2SAU 2n-1nrSe Sfn-1n-1(3-2-8).雅可比矩阵各元素的算式式(3-2-8)中，雅可比矩阵中的各元素可通过对式(3-2-4)和(3-2-5)进行偏导而求得.当j。，时，雅可比矩阵中非对角元素为SAPSAQ=

27、 -(G e + B / )SAfij 1 ij 1jrSAejg _ SAQSAfSAeSAU 2 SAU 2 八=0SeSf(3-2-9)当j = i时,雅可比矩阵中对角元素为:r (G e B f ) Ge B fdeij j ij j ii i i i TOC o 1-5 h z ij1dAPV衬-乙(G f + B e ) G f + B edfij j ij j ii i ii i HYPERLINK l bookmark133 o Current Document jj=1(3-2-10)dAQ - X (G f + B e ) G f + B e deij j ij j ii

28、i ii i HYPERLINK l bookmark161 o Current Document ij=1dAQ -X (G Ae B f ) + G e + B f dfij j ij j ii i ii ij j=1dAU 2卜dejdAU 2 idf i由式(2-2-9)和(2-2-10)看出，雅可比矩阵的特点:.矩阵中各元素是节点电压的函数，在迭代过程中，这些元素随着节点电压 的变化而变化；.导纳矩阵中的某些非对角元素为零时，雅可比矩阵中对应的元素也是为零.若七=0,则必有J = 0 ;3.雅可比矩阵不是对称矩阵;(i = q1,2,.,n;i主s)雅可比矩阵各元素的表示如下式(2-

29、2-11):(Ge + B f)dAPij i ij i(j 丰 i)H =i- sijdeX (Ge B f ) Ge B f (j = i)jij j ij j ii i ii iI jeiN -竺ijdfjLXjeiB e G f)ij iiji(j 丰 i)(G f + Be ) + Be G f (j = i)ij j ij j ii i ii idAQM iij dej=IXB e G f)ij iiji(j 丰 i)(G f + Be ) + Be G f (j = i)ij j ij j ii i ii i侦jei TOC o 1-5 h z G e +B f).dAQ Iij

30、 i ij i(ji)ij = -E-B f) + Ge +B f (j = i)jij j ，j j，,6AC/2 _ I 0 (j 药)=衬=j de 2e (j = i)ji5AE/2o (E)ij-2/ (j = i)4设计过程4.1题目D网络接线如图，各支路导纳和各节点功率已标么值标于图中。其中节点1 是按给定功率发电的发电厂。设节点5电压保持为定值1.06不变，试分析该网 络的潮流分布。方法不限，求解精度为10e-5。上 45+jO, 150.431jO 05系统接线图0.45+j0.150.08+j0.2430.4+j0.050.01+j0.0342jpo0.6+j0.1 V-(

31、0.2+j0.2)G等值阻抗电路图解:节点1为平衡节点,U1=1.06+J0为一值，其它四个节点都是PQ节点给定的注入 功率分别为：2 =0.20+J0.20,S3=-0.45-J0.15,S4=-0.40-J0.05,S5=-0.60-J0.10.形成节点导纳矩阵YBY =6.250-18.750i-5.000+15.000i -1.250+3.750i00-5.000+15.000i10.834-32.500i-1.667+5.000i-1.667+5.000i-1.250-3.750i 0 -2.500+7.500i 0 -1.250+3.750i3.750-11.250i;计算各节点功

32、率的不平衡量：取 U1=1.06+J0U2=1.00+j0U3=1.00+j0U4=1.00+j0U5=1.00+j0AP = P - e Z (G e - B f ) - f Z (G f + B e ) i i i ij ji.j j ii.j j ij jj=ij=iAQ = Q - f Z (Ge - B f ) + e Z (G f + Be ) ii iij j ij j iij j ij jj =ij =iAP2(0) = 0.5000, AQ 2(0) = 1.1000AP3(0) = -0.3750, AQ3(0) = 0.0750Ap4(0) = -0.4000, AQ4(

33、0) = -0.0500AP5(0) = -0.6.000, AQ4(0) = -0.0500计算雅可比矩阵中各元素当j丰时，雅可比矩阵中非对角元素为6A PSA Q广1 (G e + B f )SA eSAfij i a isA PSA Qi =卜B e G f TOC o 1-5 h z SAfSA eij iij ijjSA U 2 SA U 2 八0Se.Sf当j=时，雅可比矩阵中对角元素为:ijB_fj)iiii(GijijiiiiaA qae(GijijiiiiaA qaf(G AeijijiiiiaA uaeaA uafi把数据代入上边公式，可得:雅克比矩阵各元素H22=33.4

34、J22=-11.134N22=10.534L22=31.6H23=-5J23=1.667N23=-1.667L23=-5H24=-5J24=1.667N24=-1.667L24=-5H25=-7.5J25=2.5N25=-2.5L25=-7.5相似可得雅可比矩阵中其它元素33.40010.534一 5.000一 1.667一 5.000一 1.667- 7.5一 2.500一11.13431 .6001.667一 5.0001.667一 5.0002.500一 7.500- 5.000一 1.66738.97512.842- 30.000一 10.0000.0000.000J 0 =1.667

35、一 5.00012.99238.52510.000- 30.0000.0000.000一 5.000一 1.667一 30.000一 10.00038.75012.917一 3.750-1.2501.667一 5.00010.000一 30.000一 12.91738.7501.250一 3.750- 7.500一 2.5000.0000.000一 3.750-1.25011.2503.750_ 2.500一 7.5000.0000.0001.250一 3.750一 3.75011.250列写修正方程式求各节点电压值采用矩阵求各节点电压的修正量求得雅可比矩阵的逆阵节点功率不平衡量 节点电压修正

36、量，从而求得电压新值修正方程的解DY:-0.04880.0358-0.08730.0058-0.09290.0038-0.1076-0.0032节点电压的第1次近似值：1.0358 - 0.0488i1.0058 - 0.0873i1.0038 - 0.0929i0.9968 - 0.1076i1.0600各点的电压实部ei（单位:V）为（节点号从小到大排列）：1.03581.00581.00380.99681.0600各点的电压虚部fi（单位:V）为（节点号从小到大排列）：-0.0488-0.0873-0.0929-0.10760平衡节点的功率：S1 = 1.29816 + 0.24447各

37、支路功率：00.2469 + 0.0815i0.2793 + 0.0806i0.5489 + 0.1333i-0.8751 - 0.0954i-0.2431 - 0.0701i00.1891 - 0.0121i0-0.3960 - 0.0677i0.0633 + 0.0033i-0.2746 - 0.0664i-0.1887 + 0.0132i-0.5370 - 0.0977i0-0.0630 - 0.0023i000.8895 + 0.1387i0.4087 + 0.1058i000Shcrtculs bow io AdJdiscs New描总数片:翰节点为:ms =1PQ节点为:JD =2

38、356.2500 -18. 75001-5. DOOO+15.0D001-1.2500 + 3.7500100-5.0000 -+15. 0000110.成-32.5D001-1.6670 + 5.00001-l.6T0 + 5JCO01-2.5000+ 7.5D001-1.2500 + 3. 75001-LB670+ 5.0J00112.9170 -38.75001 -10. 0000 +30.0000100-1.6670+ 5.0J001-10. 0000 +30.0000112.9H0 -38.75001- 2500+ 3.750010-2.5COO+ 7.5D0010-1.2500 +

39、 3.750013.7500-11.25001导帼车为圣6,2500 -18. 75001-5. DOOO4-15.0D001-1.2500 + 3.7500100-5.0000 +15. 0000i10.制。-32.5)001-1.6570 + 5. OOOOi-1.6610 + 5. OOOOi-2.5000 + 7.5J00i-1.2500 + 3. 75001-1.6670+ 5.0J00112.970 -38.?500i -10. 0000 +30.0000100-1.6670+ 5.0D001-10. 0000 +30.00001：2,r0 -38.rEOOi-：.2500+ 3.

40、7:001涯真比矩陂35. 400010. 5340-5.QOQO-1. 60FQ-S- DOQQ-1. 667D-7_ 5000-2.5000-11. 134031.60DD1.66ro-5. DDDD1- 6670-5. DDDD2-5000-T.5DD0-5. nnnn-la 66703Sa 9750I2a sasn-30- Dnnn-1CL ODODLL 6670-5a OOOD-12a 992038.525010.DnnnTCL ODODL-5. nnnn-la 6670-30a ongn-ina nnnn焚-75nnIM 9170- 7500-la 25001. 6670-5.00

41、0010.0000-3LL OLIIJLI-12- 917038. 75001.2500-3.7500-7. 60IJU-2.6000IIIU-3. ?6ULI-1. 250DL1.2EOO3.7B002. EOOO-7.5000UU】.25 ULI-3 7500-3,760011.2500第M.欠仓iz万程的鄙iDT：-0. 04730. 0430-0. 08630. 0154-0. 09220. 0141-0. 107C0. 0093节点电压的第c (k)次近侃值:|. innni. he - |. 11 nj - n.nF-fir. ii. th i - n.np?;iI _ 】ng -

42、 n. | rrri答点的电压实部药（单位：粉为（节点号从小到大排列）:1. 06001. 04301.Ulb41. 01411. 00?3各点的电压虚部fi （单位:粉为（节点号从小到大排列）: 0-0. 0473-u. UBtiJ-0. 0922-n. inzA迭代次数：1雅克比拒F$JJ：33. 159413. 0920-E. 1360-1. 9751-5.1360-1. 9751-7. 7040-2.9621-12. RR1 133. 60831 一 gl-5. 1 3001 一 9751-5. 13602. 9621-己-4.9331-2. 124138. 384816. 0302-

43、29.5988-12. 7427002. 1241-4. 9331-1C.889038. 078812.7427-29. 598800-4.9168-2. 1516-29.5009-12.907838.155316. 2729-3. 6876-1. 61352. 15164. 916812. 9078-29.5009-17.072938. 05531. 6135-3.6876-7.3011-3. 330400-3.C50C-1. GG521 1. 05 1G4.39GC3.3304-7. 3011001.6652-3. 65055. 595610. 8516Shortcuts 直How to

44、Add 田 Whats MewEl七叭敬：2雅克比矩阵工I：32.933412. 9537-5. 0977-1. 9648-5. 0977-1. 9648-7. 6466-2. 9466-12.585033. 33711. 9648-5. 09771. 9648-5. 09772. 9466-7. 6466-4.8857-2. 098638. 049415.8301-29. 3142-12. 5897002.0986-4.8857-16.694137.679012.5897-29.314200-4.8664-2. 1233-29. 1988-12.737937.799716. 0624-3.

45、6499-1. 59222. 1233-4.866412. 7379-29.1988-16. 844637. 63051. 5922-3. 6499-7. 2110-3. 273800-3. 6055-1. 63C910.97744.32533. 2738-7. 2110001. 6369-3. 6055-5.496010.6556第M次修正方程的解DY:1. 0e-003 *-0. 0025-0. 08920. 0098-0. 12310. 0128-0. 13270. 0241-0. 1731节点电压的第C（k）次近似值：1. 06001. 0354 - 0. 0477i1. 0052 -

46、 0. 0845i1.0032 - 0.0901i0. 9961 - 0. 10441各点的电压实部疝（单位：舛为（节点号从小到大排列）：冬点的电压实部日】（单位:们为（节点号M小到大排列）：1.06001.03541.00521.Q0320.9961各点的电压虚部（单位;硝为（节点号从小到大排列）;0-0.04F7-0.0845-0.0901-0.1044迭代次教：3迭代:次食n3雅克比姑阵JJ：32. 930712. 95245.09 73-1. 964S-5,0973-1. 3646-F. 645 9一巳 9464-1ZB 5S4653. 3340E9646-5. 09751-9646-

47、5. 1973Z. 9464T 6459-4a S851-2a 098438a04 5115. 8281-29a 3 106-12a 5882UU2. U9S4-4. 8851- 16652237.674012.5082-29.3106UU-4.S65S-2. 1230-29.1950-12.73623f.795116.0602-3. 6494-1. 59202a 1230-4a 865812.7362-29a1950-IS. 84243?.S252la 592 0-3a 6494-7. 2098-3. 27310Li-3,6049-1. G36i31 0. 97644. 32433. 273

48、1-7. 209S0Li1.6566-3. 0049-5.4951ILL 6529第上次修正方程的:解D Y:E Oe-LlLir *-0. ULH3-0. 1352U. 0054-LL 18760. 0087-0. 20380. 0232-0. 2889节点电压的第C近似谊：1. 06001. U354 - 0. 04 771K DD5Z - 0. 0S451la 0032 - LL 09011U.9961 - 0. lU44iinoncuis 凶 now to ago |_fj vvnsi s new次惟正方程的航m-U.0004-0.00750. 0017-LL 01010- 0021-

49、Ua 01080.0032-0.0131节点电压的第/k）次诉位值=la 0600L 03B6 - 0-04771L 0053 - 0.0846i0033 - a.0901iD 9963 - 0- 10441各点的电压实部E r单位：仍为节点号M小到大排列）%1.06001.03551. 0053la 00330-9963各诗的电压压都以单位二V）为节点号职小到大排列）:0-Q477-D 0B46-LL 0901一 . 1044Shortcuts 回 How io Add 回 Whafs New U I -I .-4平衡节点的物率SI =1. 29S2 十 0. 2445ic -c =0H O

50、UF2 + LL U2161c 二0. 0136 + 0. 0407ic 0. 0136 + 0. 0407ic =U. U136 + 0. 0407i . 04S2 + 0, 145lLL U182 + LL 14I5l督夷暗;TLn功率00-8895+0.1387iD. 40874-0-1058iU0srei - cl 09641：|D. 24694-0-uai5iQ. 2793+0.DSOSi0. B4S9+0. 1333i-0.3960 - d 06771-0.2431一, Liroii0Q. L891一i. D121iU0-LL 2T46一|：|B 0664i-0. 18S74-0-

51、0132i0LL 0633+0. ai：i33i一 5370一 . 09771n-0.0630一0.O023i0各玄路玖-it功挛-ssgs十 . 13BF1 . clDB?上 - 1D5S851 - 0. 0954l00. 2694-a.oaisi0. 2793+0.080610. E489+D. 13331-LL 3960 - LL 06?Zl-Lia 2431一LL IJUliLi0. L891一0. 01211UU-LL 2746-l：la 0B61i-0. 1887+oa 0132iU口0633+CL O033i0-LL 537D0. 007?i0-a. 0S3ULL D023i0咎

52、支蹈由率顷耗uU-UL44 +0.0433iD. 01274-0- 038 Liu00. 014 + 0, 04331|D. 00384-0. i：U13iaH 0047 + ijH 0142iLU 01 19 + 0. 0356i0, 0127 + 0, 03S1L0-0036 4.OL131n0. 0OD4 4 0. OOl li0-004 4.0L42i0. 00044-0-OOl li0. 0003 + 0. OOlOi .DL19 4- . 0351 . DDD3 4 . DDlDi总损耗 . 0482 + CL 140.00001for m=2:5pp=0;qq=0;for n=1

53、:5pp=pp+ei(m)*(G(m,n)*ei(n)-B(m,n)*fi(n)+fi(m)*(G(m,n)*fi(n)+B(m,n)*ei(n);% Pi (0)qq=qq+fi(m)*(G(m,n)*ei(n)-B(m,n)*fi(n)-ei(m)*(G(m,n)*fi(n)+B(m,n)*ei(n); I(m)= (pp-qq*i)/conj(ei(m)+fi(m)*i);%节点注入电流 Ii(0) end DP(m)=P(m)-pp; %有功功率的不平衡量 DQ(m)=Q(m)-qq; %无功功率的不平衡量 H(m,m)=-B(m,m)*ei(m)+G(m,m)*fi(m)+imag(I(m); N(m,m)=G(m,m)*ei(m)+B(m,m)*fi(m)+real(I(m); J(m,m)=-G(m,m)*ei(m)-B(m,m)*fi(m)+real(I(m); L(m,m)=-B(m,m)*ei(m)+G(m,m)*fi(m)-imag(I(m);JJ(2*m-3,2*m-3)=H(m,m);JJ(2*m-3,2*m-2)=N(m,m);JJ(2*m-2,2*m-3)=J(m,m);JJ