电力系统潮流程序开发(共59页).doc

《电力系统潮流程序开发(共59页).doc》由会员分享，可在线阅读，更多相关《电力系统潮流程序开发(共59页).doc（59页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、精选优质文档-倾情为你奉上毕业设计（论文）任务书一、设计题目：1、题目名称 电力系统潮流计算程序开发（直角坐标形式） 2、题目来源 （自拟） 二、目的和意义电力系统潮流计算是电力系统运行分析、规划设计、稳定分析等工作的基础，研究电力系统潮流计算方法和计算手段有一定的实用价值。电力系统的结构越来越复杂，其节点规模已达数千个或上万个，应用编制计算机程序的方式来进行电力系统分析成为唯一的选择。作为从事电力系统运行、规划、设计、管理的毕业生掌握一定的应用计算机程序解决实际问题的技巧和能力是非常必要的。三、原始资料潮流方程：直角坐标形式方程求解方法：牛顿拉夫孙法四、设计说明书应包括的内容1、电网结构的表

2、示方法和读入程序；2、节点优化方法和程序3、导纳矩阵形成方法和程序；4、雅克比矩阵的形成方法和程序；5、修正方程的解法和程序；6、潮流计算结构输出程序。五、设计应完成的图纸六、主要参考资料1. 何仰赞，温增银. 电力系统分析(上、下册)(第三版). 武汉：华中科技大学出版社，20022. 杨健霑. C语言程序设计. 武汉：武汉大学出版社，20063. 苏小林，阎晓霞.电力系统分析.北京：中国电力出版社，20074. 陈亚民.电力系统计算程序及其实现. 北京：中国电力出版社，1995.11（1988重印）5. 谭浩强.C语言程序设计. 北京：清华大学出版社，2000.1（2007.5重印）七、进

3、度要求1、实习阶段 第 1周（ 月 日）至第 6 周（ 月 日）共 周2、设计阶段 第 7 周（ 月 日）至第 周（ 月 日）共 14 周3、答辩日期 第 14 周（ 年 月 日） 八、其它要求专心-专注-专业电力系统潮流程序开发摘要潮流计算是电力系统的各种计算的基础，同时它又是研究电力系统的一项重要分析功能，是进行故障计算，继电保护鉴定，安全分析的基础。电力系统潮流计算是计算电力系统动态稳定和静态稳定的基础。在电力系统规划设计和现有电力系统运行方式的研究中，都需要利用电力系统潮流计算来定量的比较供电方案或运行方式的合理性、可靠性和经济性。因此，电力系统潮流计算是电力系统中一项最基本的计算，既

4、具有一定的独立性，又是研究其他问题的基础。本文以潮流计算软件的开发设计为重点，在数学模型与计算方法的基础上，利用C+语言进行软件编写，和进行了数据测试工作，结果较为准确，收敛效果较好。设计主要采用牛顿-拉夫逊算法为算法背景。关键词：电力系统；潮流计算；牛顿-拉夫逊算法；节点优化；导纳矩阵；C+Power Flow program developmentABSTRACTThe flow calculation is calculated on the basis of a variety of power system , but it also is an important power sy

5、stem analysis function , fault calculations , relay identification, the basis of security analysis . Power flow calculation is the basis for computing power system dynamic stability and static stability . In power system planning and design and research of the existing power system operation , power

6、 flow calculation to quantify the power supply plan or run the way rationality , reliability and economy . Therefore, the power flow calculation is a basic calculation of the power system , not only has a certain independence , but also to study the basis of other issues . Trend calculation software

7、, the development and design , on the basis of the mathematical model and calculation method of using C + + language for software development, and conduct a data test work results more accurately , the convergence effect is better. Design using the Newton - Raphson algorithm algorithm background .Ke

8、ywords: Power system ; power flow calculation ;Newton - Raphson algorithm ; node optimization; the admittance matrix ; C + + 引言电力系统的传统潮流算法一般适合输电网，由于配电网的特殊结构：配电网具有环形结构, 而通常以开环方式运行。通常呈辐射状，支路比值较大，分支线较多；配电线路中的R/X 比值偏大使输电网中常用的潮流计算算法如传统的牛顿法和快速分解法在应用于配电网潮流计算时容易形成病态而无法收敛，所以需要根据传统的潮流算法进行改进以适应低压配电网，如前推回代法、并行算

9、法、分层前推回代法等。目的和意义：电力系统潮流计算是研究电力系统稳态运行情况的一种基本电气计算，是电力系统规划和运营中不可缺少的一个重要组成部分。可以说，潮流计算是电力系统分析中最基本、最重要的计算，是电力系统安全、经济分析和实时控制与调度的基础。潮流计算的任务是根据给定的运行条件和网络结构确定整个系统的运行状态，如各母线上的电压（幅值及相角）、网络中的功率分布及功率损耗等。即潮流计算是对电力系统的功率分布和电压分布的计算, 其具体任务就是编制系统的调度计划和电气设备检修计划, 确定电力系统中变压器分接头位置和系统中枢点与电压控制点的电压曲线， 进行事故运行方式的分析, 为电力系统短路和稳定计

10、算提供数据, 为继电保护及自动装置整定与电力系统设计和规划提供依据等。潮流计算的目的是对现有电力系统的正常运行状态进行分析, 以提示必要的改进措施, 同时为新建系统或扩建系统的有关分析、计算打下基础。配电网潮流计算是配电网经济运行、系统分析等的重要基础，但由于配电网与输电网有着明显的差异:配电网具有环形结构, 而通常以开环方式运行。通常呈辐射状，支路比值较大，分支线较多；配电线路中的R/X 比值偏大使输电网中常用的潮流计算算法如传统的牛顿法和快速分解法在应用于配电网潮流计算时容易形成病态而无法收敛，因此，研究适合于配电网的潮流算法也是至关重要的。目前，输电系统潮流计算方法已较为成熟，而且获得了

11、广泛的实际应用。但随着电力系统规模的不断扩大，潮流方程的阶数越来越高，对这种规模的方程并不是采用任何数学方法都能保证给出正确答案的，因此，这也成为促使电力系统研究人员不断寻求新的、更可靠的潮流计算方法的动力。 随着现代电力系统大系统、强非线性与多元件的特点日益突出, 其计算量与计算复杂度急剧增加。旧的计算机软件在处理潮流计算时, 其速度已无法满足大电网模拟和实时控制的仿真要求, 而高效的潮流问题的相关软件的研究已成为大规模电力系统仿真计算的关键。第1章 题目的名称及意义设计题目为电力系统潮流程序开发。根据所给的电力系统，编制潮流计算程序，通过计算机进行调试，最后完成一个切实可行的电力系统计算应

12、用程序。通过自己设计电力系统计算程序使同学对电力系统分析有进一步理解，同时加强计算机实际应用能力的训练。1.1潮流计算简介应用计算机进行电力系统计算，首先要掌握电力系统相应计算的数学模型；其次是运用合理的计算方法；第三则是选择合适的计算机语言编制计算程序。建立电力系统计算的相关数学模型，就是建立用于描述电力系统相应计算的有关参数间的相互关系的数学方程式。该数学模型的建立往往要突出问题的主要方面，即考虑影响问题的主要因素，而忽略一些次要因素，使数学模型既能正确地反映实际问题，又使计算不过于复杂。运用合理的计算方法，就是要求所选用的计算方法能快速准确地得出正确结果，同时还应要求在解算过程中占用内存

13、少，以利提高计算机的解题规模。选择合适的语言编写程序，就是首先确定用什么计算机语言来编制程序；其次是作出计算的流程图；第三根据流程图用选择的语言编写计算程序。然后上机调试，直到语法上无错误。本程序采用C+进行编程。所编制的程序难免存在逻辑错误，因此先用一个已知结果的系统作为例题进行计算。用程序计算的结果和已知结果相比较，如果结果相差甚远就要逐步分析程序的计算步骤，查出问题的出处；如果结果比较接近，则逐步分析误差来源；直到结果正确为止。1.2潮流计算的意义其意义在于：1) 在电网规划阶段,通过潮流计算,合理规划电源容量及接入点,合理规划网架,选择无功补偿方案,满足规划水平的大、小方式下潮流交换控

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

15、系统的运行状态，也需要进行大量而快速的潮流计算。因此，潮流计算是电力系统中应用最广泛、最基本和最重要的一种电气运算。在系统规划设计和安排系统的运行方式时，采用离线潮流计算；在电力系统运行状态的实时监控中，则采用在线潮流计算。 第2章 潮流计算的方法和步骤潮流计算是电力系统分析中的一种最基本的计算，它的任务是对给定的运行条件确定系统的运行状态，如母线上的电压（幅值及相角）、网络中的功率分布及功率损耗等。目前计算机潮流计算的方法主要有牛顿-拉夫逊算法和PQ分解法。牛顿-拉夫逊算法是数学上求解非线形方程组的有效方法，具有较好的收敛性，是潮流计算中应用比较普遍的方法。 牛顿-拉夫逊法计算电力系统潮流的

16、基本步骤：1) 形成节点导纳矩阵；2) 给各节点电压设初值；3) 将节点电压初值代入，求出修正方程式是常数项向量；4) 将节点电压初值代入，求出雅可比矩阵元素；5) 求解修正方程式，求出变量的修正向量；6) 求出节点电压的新值；7) 如有PV节点，则检查该类节点的无功功率是否越限；8) 检查是否收敛，如不收敛，则以各节点电压的新值作为初值自第3步重新开始下一次迭代，否则转入下一步。1) 计算支路功率分布，PV节点无功功率和平衡节点注入功率，最后输出结果，并结束其流程图如下图所示： 牛顿拉夫逊法求解潮流程序流程图（直角坐标法）开 始读入网络参数及系统运行参数数据文件形成节点导纳矩阵初 始 化 设

17、定电压初值V（0）=1，（0）=0，最大迭代数Kmax=50，迭代记数器k=0计算节点功率不平衡量P，Qmax|Pi，Qi|=e，i=1,2,n 打印不收敛Y判断kKmax 或eKmaxNN 计算支路功率，形成雅可比矩阵，修正方程，求出V， 输出潮流结果计算，k=k+1 结束计算节点功率不平衡量P，Qmax|Pi，Qi|=e，i=1,2,第3章 原始数据读入和存储程序开发电力系统原始数据是电力系统计算的基础。电力系统每个计算程序都要求输入一定的原始数据，这些数据可以反映电力网络结构、电力系统正常运行条件、电力系统各元件参数和特性曲线。不同的计算程序需要不用的原始数据。3.1电力网络的描述电力网

18、络是由输电线路、电力变压器、电容器和电抗器等元件组成。这些元件用集中参数的电阻、电抗和电容表示。为了表示电力网络中各元件是怎样互相连接的，要对网络节点进行编号。电力网络的结构和参数由电力网络中各支路的特性来描述。3.1.1 线路参数在本程序设计中，线路参数采用线路的型数学模型，即线路用节点间的阻抗和节点对地容性电纳来表示，由于线路的对地电导很小，可忽略不计。其等价回路如下： r+jx i j-jb -jb对于线路参数的数据文件格式可写为：线路参数（1，节点i，节点j，r，x，b）3.1.2 变压器参数在本程序设计中，变压器参数采用型等值变压器模型，这是一种可等值地体现变压器电压变换功能的模型。

19、在多电压级网络计算中采用这种变压器模型后，就可不必进行参数和变量的归算。双绕组变压器的等值回路如下： k ZT k:1 ZT (a)接入理想变压器后的等值电路 (b) 等值电路以导纳表示 YT/k (c) 等值电路以导纳表示因此，对于变压器参数的数据文件格式可写为：变压器参数（2，节点i，节点j，r，x，k）其中，k表示变压器变比。3.1.3 相关程序按上述线路、变压器参数形式读入原始数据到datafile文本文档中，然后按所给参数生成支路实例,并将生成的支路对象指针保存到数组。if(sr0!=/ & sr0!=!) /判断是否为注释行tok=strtok_s(sr,&ntok);int t=

20、atoi(tok);/获得数据类型if(t=1 | t=2) /线路或变压器数据tok=strtok_s(NULL,&ntok);int n1=atoi(tok);/第一节点号tok=strtok_s(NULL,&ntok);int n2=atoi(tok);/第二节点号tok=strtok_s(NULL,&ntok);double r=atof(tok);/电阻tok=strtok_s(NULL,&ntok);double x=atof(tok);/电抗tok=strtok_s(NULL,&ntok);double b=atof(tok);/电纳或变比PBranch *pB=new PBra

21、nch(n1,n2);/生成支路实例pB-SetR(r);/赋值RpB-SetX(x);/赋值XpB-SetB(b);/赋值BpB-SetType(t-1);/支路类型BranchArray.Add(pB);/将生成的支路对象指针保存到数组3.2电力系统运行条件数据电力系统运行条件数据包括发电机（含调相机）所连接的节点号、有功与无功功率；负荷所连接的节点号、有功与无功功率；PV节点与给定电压值；平衡节点的节点号与给定电压值。3.2.1节点功率参数电力系统中有流入流出功率的称为功率节点，有流入功率的称发电节点，一般为各发电站、枢纽变电站等节点；有流出功率的称负荷节点。对于电力系统稳态计算来说，功

22、率节点都用有功功率P和无功功率Q来简单表示。其等价回路如下：QGPG PLQL3.2.2 节点类型参数根据给定节点变量的不同，可以有以下三种类型的节点：1. PV节点（电压控制母线） 这种节点的注入有功功率Pi为给定值，电压Ui也保持在给定数值。这种类型节点相当于发电机母线节点，其注入的有功功率由汽轮机调速器设定，而电压则大小由装在发电机上的励磁调节器控制；或者相应于一个装有调相机或静止补偿器的变电所母线，其电压由可调无功功率的控制器设定。要求有连续可调的无功设备，调无功来调电压值。2. PQ节点 这种节点的注入有功和无功功率是给定的，相应于实际电力系统中的一个负荷节点，或有功和无功功率给定的

23、发电机母线。3. 平衡节点这种节点用来平衡全电网的功率，一般选用一容量足够大的发电厂（通常是承担系统调频任务的发电厂）来担任。平衡节点的电压和相位大小是给定的，通常以它的相角为参考量，即取其电压相角为0。一个独立的电力网络只设一个平衡节点。三类节点的划分并不是绝对不变的。PV节点之所以能控制其节点的电压为某一设定值，重要原因在于它具有可调节的无功功率出力。一旦它的无功功率出力达到可调节的上限或下限，就不能使电压保持在设定值，PV节点将转化成PQ节点。对于这三种类型的节点参数可如下表示平衡节点：给出节点编号，节点电压。节点功率数据（3，节点i，3, Ui）PQ节点：在节点功率参数中就可表示。节点

24、功率数据（3，节点i，1, PLi，QLi）PV节点：需单列，其数据文件格式可写为：节点功率数据（3，节点i，2, PLi，Ui）3.2.3 相关程序首先对节点号进行类的定义，然后按节点参数形式读入原始数据到datafile文本文档中，然后按所给参数生成支路实例,并将生成的支路对象指针保存到数组。相关程序：else if(t=3)/节点数据tok=strtok_s(NULL,&ntok);int n=atoi(tok);/节点号PNode *Node=NodeExist(n);if(Node)tok=strtok_s(NULL,&ntok);int nt=atoi(tok);/节点类型if(n

25、t=1)/PQ节点tok=strtok_s(NULL,&ntok);double s1=atof(tok);/有功设定tok=strtok_s(NULL,&ntok);double s2=atof(tok);/无功设定Node-SetS1(s1);Node-SetS2(s2);else if(nt=2)/PV节点tok=strtok_s(NULL,&ntok);double s1=atof(tok);/有功设定tok=strtok_s(NULL,&ntok);double s2=atof(tok);/电压设定Node-SetS1(s1);Node-SetS2(s2);Node-SetVolta

26、ge(s2,0.0);Node-SetNodeType(1);else if(nt=3)/平衡节点tok=strtok_s(NULL,&ntok);double s1=atof(tok);/电压设定Node-SetS1(s1);Node-SetVoltage(s1,0.0);Node-SetNodeType(2);3.3各类数据文件格式潮流数据文件格式1) 节点数，平衡节点，平衡节点电压，计算精度2) 线路参数（1，节点i，节点j，r，x，b）3) 变压器参数（2，节点i，节点j，r，x，k）4) 节点功率数据（序号，节点i，nt，S1，S2）第4章 节点优化程序开发4.1 节点优化标号及方法

27、节点的编号顺序对于计算效力的影响至关重要，特别是采用了稀疏技术后，它直接影响到矩阵A的因子表矩阵的稀疏度。严格的说，最优编号是一个组合优化问题，求其最优解是困难的，但在实际工程中，有许多实用的次优的编号方法得到了广泛的应用。根据节点优化编号实现的复杂程度和最终的编号效果的不同，可有如下分类：1. Tinney-1编号方法这种方法也称静态节点优化编号方法。这种方法在有向图上统计每一个节点的出线度，即该节点和其他节点相连结的支路树，然后按节点出线度由小到大按顺序进行编号。对于出线度相同的节点，哪个排在前边是任意的。这种编号方法的出发点是认为在图上因子分解的过程中出线度小的节点消去时产生新因子的可能

28、性也小。这种编号方法简单，但编号效果较差。2. Tinney-2编号方法这种方法也称最小度算法，或半动态节点优化编号方法。这种方法首先统计所有节点的出线度，然后选择出线度最小的节点进行编号。编号过程中，按图上因子分解的方法消去该节点，只进行网络结构变化的处理，而不进行边权计算。然后消去已编号的节点和其相关支路，在剩下的子图上重复上述编号过程。这种方法也比较简单，图上因子分解产生新支路以及处理过的支路这些变化可用在原来的图上修正来实现。这种编号方法可使有向因子图上新增加的支路数大大减少，而程序复杂性和计算量又增加不多，是一种使用十分广泛的编号方法。其步骤如下：1) 网络节点进行随意的人工编号。2

29、) 统计原始网络各节点所连接的支路数，并记存各节点所连接支路对端的节点号。3) 令新的节点号I=1。4) 在尚未编号网络中，查找连接支路数最少的节点J，将其编号取为节点号I。5) 消去J节点。其效应有二：去掉与J节点相连接的所有支路，也就是使与J节点相连接的所有节点连接的支路数各自减一，并去掉支路对端的节点号J。使原来与J节点相连接的所有节点每两个之间如果原来没有连接支路，则增加一条新的支路，同时，新支路两端的节点各自记存对端的节点号。6) I=I+1；7) 判别I是否对于N（网络节点数）。若对于N，则节点编号优化结束；否则转到步骤）。3. Tinney-3编号方法这种方法也称动态节点优化编号

30、方法。它和上面的Tinney-2编号方法的不同之处是对所有待编号的节点，统计消去该节点时产生的新支路的数目，并以该数目最小为优先编号的准则。某一节点编号完成之后，要立即修改因子图。其优化步骤为：1) 将n个节点网络的每个节点轮流进行一次消去运算，统计各节点消去后各自增加新的支路数，将增加新支路数最少的节点编为第号，随后消去该节点。2) 将n-1个节点网络的每个节点轮流进行一次消去运算，统计各节点消去后各自增加新的支路数，将增加新支路数最少的节点编为第2号，随后消去该节点。依次类推，进行n步操作，完成了节点编号优化。从理论上说，这种方法效果最好，但在每步编号前后要对所有待编号节点统计消去后产生的

31、新支路数，程序复杂程度和编号时的计算量都很大，所以不常用。4.2 相关程序1）节点中只有一个平衡节点，将它放到最后的位置，即先找到它，与最后一个节点进行交换。OpNodeArray.Append(NodeArray);for(int i=0;iGetNodeType()=2)/如果节点为平衡节点，则与最后一个节点对换位置OpNodeArrayi=OpNodeArrayOpNodeArray.GetSize()-1;OpNodeArrayOpNodeArray.GetSize()-1=Nodei;break;2）找出关联节点最少的节点，放在第i个位置上。上述步骤4）for(int i=0;iOp

32、NodeArray.GetSize()-2;i+)/找出关联节点最少的节点，放在第i个位置上PNode *Nodei=OpNodeArrayi;for(int j=i+1;jGetAssNodeNum()GetAssNodeNum()OpNodeArrayi=Nodej;OpNodeArrayj=Nodei;Nodei=OpNodeArrayi;3）消去J节点/模拟i节点的消去过程，修改剩余节点关联的节点int ND=Nodei-GetAssNodeNum();for(int j=0;jGetAssNode(j);for(int k=j+1;kGetAssNode(k);Node1-AddAs

33、sNode(Node2);Node2-AddAssNode(Node1);for(int i=0;iSetNumNew(i+1);第5章 导纳矩阵形成程序开发5.1节点导纳矩阵其中：对角元Yii称为自导纳，数值上等于该节点直接连接的所有支路导纳的总和；非对角元Yij称为互导纳，数值上等于连接节点i，j支路导纳的负值。N个节点的电力网络的节点导纳矩阵的特点：1） nn阶方阵；2） 对称；3） 复数矩阵；4） 每一非对角元素Yij是节点i和j间支路导纳的负值，当i和j间没有直接相连的支路时，为0。根据一般电力系统的特点，每一节点平均与3-5个相邻节点有直接联系，所以导纳矩阵是一高度稀疏矩阵。互导纳

34、，不包括对地支路；5） 对角元素Yii为所有联结于节点i的支路的导纳之和。5.2 形成节点导纳矩阵的方法节点导纳矩阵的自导纳（对角元素）等于节点所连支路的导纳之和，互导纳（非对角元素）等于两个节点之间连接支路导纳取负号。因此，采用支路追加法来形成导纳矩阵。上章进行节点优化后对生成支路进行排序。void BranchOrder()/支路排序/支路节点互换for(int i=0;iGetNodeAddr1();PNode *Node2=Branch-GetNodeAddr2();if(Node1 & Node2)/两个节点号都不为if(Node1-GetNumNew()Node2-GetNumNe

35、w()Branch-SetNode1(Node2);Branch-SetNode2(Node1);if(Branch-GetType()=1)/变压器支路,变比符号变换double k=Branch-GetK();Branch-SetK(-k);else if(!Node1 | !Node2)/有一个节点号为if(!Node2)Branch-SetNode1(Node2);Branch-SetNode2(Node1);/支路按首节点号排序for(int i=0;iGetNodeAddr1();for(int j=i+1;jGetNodeAddr1();if(Nodei-GetNumNew()N

36、odej-GetNumNew()BranchArrayi=Branchj;BranchArrayj=Branchi;Branchi=BranchArrayi;Nodei=Branchi-GetNodeAddr1();/首节点号相同。按末节点号排序for(int i=0;iGetNodeAddr1();PNode *Nodei2=Branchi-GetNodeAddr2();for(int j=i+1;jGetNodeAddr1();PNode *Nodej2=Branchj-GetNodeAddr2();if(Nodei1-GetNumNew()!=Nodej1-GetNumNew()brea

37、k;if(Nodei2-GetNumNew()Nodej2-GetNumNew()BranchArrayi=Branchj;BranchArrayj=Branchi;Branchi=BranchArrayi;Nodei1=Branchi-GetNodeAddr1();Nodei2=Branchi-GetNodeAddr2();5.33.1.1节点导纳矩阵的存储其为高度稀疏的N阶复数对称方阵。因此记录矩阵的下三角即可。 数组表示法：数组1：记录矩阵对角元素的数值；数组2：记录矩阵非对角元素的数值（按列存储） ；数组3：记录矩阵非对角元素的行号；数组4：记录矩阵非对角元素的按行排的位置数；数组5：

38、记录矩阵非对角元素的按行存储对应按列存储的位置数。void FormY()/形成导纳矩阵int nd=OpNodeArray.GetSize();int nb=BranchArray.GetSize();YDArray.SetSize(nd);YIArray.SetSize(nd);YUArray.SetSize(nb);YJArray.SetSize(nb);int count=0;for(int i=0;ind;i+)YDArrayi.g=0.0;YDArrayi.b=0.0;YIArrayi=0;double R,X,Z,K;for(int i=0;iGetNodeAddr1();PNo

39、de *Node2=Branch-GetNodeAddr2();if(Branch-GetType()=0)/普通支路if(Node1=NULL)/接地支路int n=Node2-GetNumNew();R=Branch-GetR();X=Branch-GetX();Z=R*R+X*X;YDArrayn-1.g+=R/Z;YDArrayn-1.b+=-X/Z;elseint n1=Node1-GetNumNew();int n2=Node2-GetNumNew();R=Branch-GetR();X=Branch-GetX();Z=R*R+X*X;/互导纳YUArraycount.g=-R/Z;YUArraycount.b=X/Z;YJArraycount=n2-1;YIArrayn1-1+;count+;/自导纳YDArrayn1-1.g