电力系统稳态分析课程设计——清华大学.doc

清华大学电力系统稳态分析课程设计说明书题 目：两机五节点网络潮流计算牛拉法学生姓名：*学 号：*专 业：电气工程*班 级：电气1班指导教师：*34清华电力系统稳态分析课程设计说明书基于 C语言的电力系统潮流计算摘 要潮流计算是电力系统的一项重要分析功能, 是进行故障计算, 继电保护整定, 安全分析的必要工具。另外，传统的程序设计方法是结构化程序设计方法，该方法基于功能分解，把整个软件工程看作是一个个对象的组合，由于对某个特定问题域来说，该对象组成基本不变，因此，这种基于对象分解方法设计的软件结构上比较稳定，易于维护和扩充。本设计主要阐述了电力系统潮流分布，设计及运行的理论和计算方法。内容包括：电力系统的功率平衡，电力网的接线方案的技术论证及经济比较，又包括电压等级的选择，电力网接线方案的初步选择，电力网接线方案的最优选择，发电厂及变电所主接线的选择，发电厂电气主接线的选择及电气设备的选择校验，又包括主变压器的选择，各级电压接线方式的设计，短路电流计算，电气设备的选择及校验，电力网的潮流与电压计算，包括最大最小负荷时的运算负荷及电源运算功率，变压器分接头的选择。关键词：电力系统；潮流分析；C语言清华大学电力系统稳态分析课程设计说明书Power Flow Calculation Based on Language CAbstractFlow calculation is an important analysis function of power system and is the necessary facility of fault analysis, relay protection setting and security analysis. In addition, the traditional design method is a structured program design method based on functional decomposition, the entire software engineering as a combination of objects, as the domain of a particular issue, the composition of the object will remain basically unchanged Therefore, this decomposition method based on object design software structure relatively stable, easy to maintain and expand. Describes the design of the main electrical power plant and a substation system structure, design and operation of the theory and calculation methods. Include: power system power balance, power grid connection schemes of technical feasibility and economic comparison, but also the choice of voltage, power grid connection schemes of the initial choice of grid connection schemes the optimal choice, Power plants and substations main connection options, the choice of power plant electrical main wiring and electrical equipment selection verification, including the main transformer and the choice of the design at all levels of voltage wiring, short circuit current calculation, electrical equipment selection and check, the trend of power network and voltage calculations, including computing the maximum and minimum load power load and power operation, transformer tap selection.Key words: Electric power system; flow calculation; LanguageC清华大学电力系统稳态分析课程设计说明书目 录摘 要IAbstractII第一章 简介和原理说明21.1 研究背景2 1.1.1 C语言简介2 1.1.2电力系统概述 . 5 1.1.3电力系统的组成61.2 牛顿拉法逊法概述7 1.2.1牛顿-拉夫逊基本原理.7 1.2.2求解过程大致可以分为以下步骤81.2.3以直角坐标系形式表示计算过 9 第二章 程序及其说明142.1设计资料及参数14 2.1.1牛顿-拉夫逊法潮流计算程序流程图172.1.2C语言编写潮流计算程序.18 2.1.3程序运行结果322.2结论和总结33总结 33 参考文献.34简介和原理说明1.1 研究背景电力系统潮流计算是研究电力系统稳态运行中最为常用的基本计算，即根据给定电力系统的网络结构和参数，在满足电力系统运行状态参数的边界条件情况下，确定电力系统稳态运行状态的最基本的方法。潮流计算的研究是从20世纪50年代随着电网的产生而开始的，涌现出各种算法都是围绕着算法的可靠性或收敛性；对计算机内存量的要求和计算速度；计算的方便性和灵活性等。 C语言是一种结构化语言。它层次清晰，便于按模块化方式组织程序，易于调试和维护。C语言的表现能力和处理能力极强。它不仅具有丰富的运算符和数据类型，便于实现各类复杂的数据结构。它还可以直接访问内存的物理地址，进行位(bit)一级的操作。由于C语言实现了对硬件的编程操作，因此C语言集高级语言和低级语言的功能于一体。既可用于系统软件的开发，也适合于应用软件的开发。此外，C语言还具有效率高，可移植性强等特点。因此广泛地移植到了各类各型计算机上，从而形成了多种版本的C语言。1.1.1C语言简介C语言是在70年代初问世的。一九七八年由美国电话电报公司(AT&T)贝尔实验室正式发表了C语言。同时由B.W.Kernighan和D.M.Ritchit合著了著名的“THE C PROGRAMMING LANGUAGE”一书。通常简称为K&R，也有人称之为K&R标准。但是，在K&R中并没有定义一个完整的标准C语言，后来由美国国家标准学会在此基础上制定了一个C 语言标准，于一九八三年发表。通常称之为ANSI C。当代最优秀的程序设计语言。早期的C语言主要是用于UNIX系统。由于C语言的强大功能和各方面的优点逐渐为人们认识，到了八十年代，C开始进入其它操作系统，并很快在各类大、中、小和微型计算机上得到了广泛的使用。成为当代最优秀的程序设计语言之一。include称为文件包含命令扩展名为.h的文件也称为头文件或首部文件定义两个实数变量，以被后面程序使用显示提示信息从键盘获得一个实数x求x的正弦，并把它赋给变量s显示程序运算结果main函数结束程序的功能是从键盘输入一个数x，求x的正弦值，然后输出结果。在main()之前的两行称为预处理命令(详见后面)。预处理命令还有其它几种，这里的include 称为文件包含命令，其意义是把尖括号或引号内指定的文件包含到本程序来，成为本程序的一部分。被包含的文件通常是由系统提供的，其扩展名为.h。因此也称为头文件或首部文件。C语言的头文件中包括了各个标准库函数的函数原型。因此，凡是在程序中调用一个库函数时，都必须包含该函数原型所在的头文件。在本例中，使用了三个库函数：输入函数scanf，正弦函数sin,输出函数printf。sin函数是数学函数，其头文件为math.h文件，因此在程序的主函数前用include命令包含了math.h。scanf和printf是标准输入输出函数，其头文件为stdio.h，在主函数前也用include命令包含了stdio.h文件。 需要说明的是，C语言规定对scanf和printf这两个函数可以省去对其头文件的包含命令。所以在本例中也可以删去第二行的包含命令#include。同样，在例1.1中使用了printf函数，也省略了包含命令。在例题中的主函数体中又分为两部分，一部分为说明部分，另一部分执行部分。说明是指变量的类型说明。例题中未使用任何变量，因此无说明部分。C语言规定，源程序中所有用到的变量都必须先说明，后使用，否则将会出错。这一点是编译型高级程序设计语言的一个特点，与解释型的BASIC语言是不同的。说明部分是C源程序结构中很重要的组成部分。本例中使用了两个变量x，s，用来表示输入的自变量和sin函数值。由于sin函数要求这两个量必须是双精度浮点型，故用类型说明符double来说明这两个变量。说明部分后的四行为执行部分或称为执行语句部分，用以完成程序的功能。执行部分的第一行是输出语句，调用printf函数在显示器上输出提示字符串，请操作人员输入自变量x的值。第二行为输入语句，调用scanf函数，接受键盘上输入的数并存入变量x中。第三行是调用sin函数并把函数值送到变量s中。第四行是用printf 函数输出变量s的值，即x的正弦值。程序结束。输入和输出函数在前两个例子中用到了输入和输出函数scanf和 printf，在第三章中我们要详细介绍。这里我们先简单介绍一下它们的格式，以便下面使用。scanf和 printf这两个函数分别称为格式输入函数和格式输出函数。其意义是按指定的格式输入输出值。因此，这两个函数在括号中的参数表都由以下两部分组成： “格式控制串”，参数表 格式控制串是一个字符串，必须用双引号括起来，它表示了输入输出量的数据类型。各种类型的格式表示法可参阅第三章。在printf函数中还可以在格式控制串内出现非格式控制字符，这时在显示屏幕上将原文照印。参数表中给出了输入或输出的量。当有多个量时，用逗号间隔。例如：printf(sine of %lf is %lfn,x,s);其中%lf为格式字符，表示按双精度浮点数处理。它在格式串中两次现，对应了x和s两个变量。其余字符为非格式字符则照原样输出在屏幕上1.1.2电力系统概述1831年法拉第发现了电磁感应定律。在此基础上，很快出现了原始的交流发电机、直流发电机和直流电动机。由于当时发电机发出的电能仅用于电化学工业和电弧灯，而电动机所需的电能又来自蓄电池，电机制造和电力输送技术的发展最初集中于直流电。原始的 电力线路使用的就是100400V低压直流电。由于输电电压低，输送的距离不可能远，输送的功率也不可能大。与100余年前电力系统的雏形相比，近代电力系统不仅在输电电压、输送距离、输送功率等方面有了千百倍的增长，而且电源构成、负荷成分等方面也有了很大变化。电力系统加上发电机的原动机(如汽轮机、水轮机)，原动机的力能部分(如热力锅炉、水库、原子能电站的反应堆)、供热和用热设备，则称为动力系统。电力工业发展初期，电能是直接在用户附近的发电站(或称发电厂)中生产的，各发电站孤立运行。随着工农业生产和城市的发展，电能的需要量迅速增加，而热能资源(如煤田)和水能资源丰富的地区又往往远离用电比较集中的城市和工矿区，为了解决这个矛盾，就需要在动力资源丰富的地区建立大型发电站，然后将电能远距离输送给电力用户。同时，为了提高供电可靠性以及资源利用的综合经济性，又把许多分散的各种形式的发电站，通过送电线路和变电所联系起来。这种由发电机、升压和降压变电所，送电线路以及用电设备有机连接起来的整体，即称为电力系统。现代电力系统提出了“灵活交流输电与新型直流输电”的概念。灵活交流输电技术是指运用固态电子器件与现代自动控制技术对交流电网的电压、相位角、阻抗、功率以及电路的通断进行实时闭环控制，从而提高高压输电线路的输送能力和电力系统的稳定水平。新型直流输电技术是指应用现电力电子技术的最新成果，改善和简化变流站的造价等。 运行方式管理中，潮流是确定电网运行方式的基本出发点；在规划领域，需要进行潮流分析验证规划方案的合理性；在实时运行环境，调度员潮流提供了电网在预想操作情况下电网的潮流分布以校验运行可靠性。在电力系统调度运行的多个领域都涉及到电网潮流计算。潮流是确定电力网络运行状态的基本因素，潮流问题是研究电力系统稳态问题的基础和前经验主义的研究方法1.1.3电力系统的组成电力系统中生产、变换、输送、消费电能的四大部分发电机组、变压器、电力线路、负荷特性。描述一个电力系统需要了解它的产生过程，以及一些基本概念、参数等。1、原始能源火电：锅炉汽轮机发电机 水电：水库水轮机发电机核电：核反应堆汽轮机发电机其它：如风能、地热能、太阳能、潮汐等2、基本概念电力网络是由变压器、电力线路等变换、输送、分配电能设备所组成的部分。动力系统在电力系统的基础上，把发电厂的动力部分（例如火力发电厂的锅炉、汽轮机和水力发电厂的水库、水轮机以及核动力发电厂的反应堆等）包含在内的系统。总装机容量指该系统中实际安装的发电机组额定有功功率的总和，以千瓦（KW）、兆瓦（MW）、吉瓦（GW）为单位计。年发电量指该系统中所有发电机组全年实际发出电能的总和，以千瓦时（KWh）、兆瓦时（MWh）、吉瓦时（GWh）为单位计。 最大负荷指规定时间内，电力系统总有功功率负荷的最大值，以千瓦（KW）、兆瓦（MW）、吉瓦（GW）为单位计。 额定频率按国家标准规定，我国所有交流电力系统的额定功率为50Hz。 最高电压等级是指该系统中最高的电压等级电力线路的额定电压。1.2 牛顿拉法逊法概述牛顿拉夫逊潮流计算的核心是修正方程的建立和求解。潮流计算是电力系统最基本，最常用的计算。在数学上是多元非线性方程组的求解问题，求解的方法有很多种。Newton-Raphson法是数学上解非线性方程式的有效方法，有较好的收敛性。将N-R法用于潮流计算是以导纳矩阵为基础的，由于利用了导纳矩阵的对称性、稀疏性及节点编号顺序优化等技巧，使N-R法在收敛性、占用内存、计算速度等方面都达到了一定的要求。 电力系统潮流计算也分为离线计算和在线计算两种，前者主要用于系统规划设计和安排系统的运行方式，后者则用于正在运行系统的经常监视及实时控制。 利用电子数字计算机进行电力系统潮流计算从50年代中期就已经开始。在这20年内，潮流计算曾采用了各种不同的方法，这些方法的发展主要围绕着对潮流计算的一些基本要求进行的。对潮流计算的要求可以归纳为下面几点：（1）计算方法的可靠性或收敛性；（2）对计算机内存量的要求；（3）计算速度；（4）计算的方便性和灵活性。电力系统潮流计算问题在数学上是一组多元非线性方程式求解问题，其解法都离不开迭代。因此，对潮流计算方法，首先要求它能可靠地收敛，并给出正确答案。由于电力系统结构及参数的一些特点，并且随着电力系统不断扩大，潮流计算的方程式阶数也越来越高，对这样的方程式并不是任何数学方法都能保证给出正确答案的。这种情况成为促使电力系统计算人员不断寻求新的更可靠方法的重要因素1.2.1牛顿-拉夫逊基本原理以下讨论的是用直角坐标形式的牛顿拉夫逊法潮流的求解过程。当采用直角坐标时，潮流问题的待求量为各节点电压的实部和虚部两个分量由于平衡节点的电压向量是给定的，因此待求两共需要2(n-1)个方程式。事实上，除了平衡节点的功率方程式在迭代过程中没有约束作用以外，其余每个节点都可以列出两个方程式。对PQ节点来说，是给定的，因而可以写出 （2-2-1）对PV节点来说，给定量是，因此可以列出 （2-2-2）1.2.2求解过程大致可以分为以下步骤（1）形成节点导纳矩阵（2）将各节点电压设初值U，（3）将节点初值代入相关求式，求出修正方程式的常数项向量（4）将节点电压初值代入求式，求出雅可比矩阵元素（5）求解修正方程，求修正向量（6）求取节点电压的新值（7）检查是否收敛，如不收敛，则以各节点电压的新值作为初值自第3步重新开始进行狭义次迭代，否则转入下一步（8）计算支路功率分布，PV节点无功功率和平衡节点柱入功率。1.2.3以直角坐标系形式表示计算过程. 迭代推算式 采用直角坐标时,节点电压相量及复数导纳可表示为: (2-2-3)将以上二关系式代入上式中,展开并分开实部和虚部;假定系统中的第1,2,m号为PQ节点,第m+1,m+2,n-1为PV节点,根据节点性质的不同,得到如下迭代推算式: 对于PQ节点 (2-2-4)对于PV节点 (2-2-5)对于平衡节点平衡节点只设一个,电压为已知,不参见迭代,其电压为: (2-2-6). 修正方程式迭代式共包括2(n-1)个方程.选定电压初值及变量修正量符号,代入方程并按泰勒级数展开,略去二次方程及以后各项,得到一组线性方程组或线性化了的方程组，常称修正方程组: (2-2-7) (2-2-8).雅可比矩阵各元素的算式式(2-2-8)中, 雅可比矩阵中的各元素可通过对式(2-2-4)和(2-2-5)进行偏导而求得.当时, 雅可比矩阵中非对角元素为 (2-2-9)当时,雅可比矩阵中对角元素为: (2-2-10)由式(2-2-9和(2-2-10)看出,雅可比矩阵的特点:矩阵中各元素是节点电压的函数,在迭代过程中,这些元素随着节点电压的变化而变化;导纳矩阵中的某些非对角元素为零时,雅可比矩阵中对应的元素也是为零.若,则必有;雅可比矩阵不是对称矩阵;雅可比矩阵各元素的表示如下: 潮流计算的意义及其发展电力系统潮流计算是电力系统分析中的一种最基本的计算，是对复杂电力系统正常和故障条件下稳态运行状态的计算。潮流计算的目标是求取电力系统在给定运行状态的计算。即节点电压和功率分布，用以检查系统各元件是否过负荷。各点电压是否满足要求，功率的分布和分配是否合理以及功率损耗等。对现有电力系统的运行和扩建，对新的电力系统进行规划设计以及对电力系统进行静态和暂态稳定分析都是以潮流计算为基础。潮流计算结果可用如电力系统稳态研究，安全估计或最优潮流等对潮流计算的模型和方法有直接影响。实际电力系统的潮流技术那主要采用牛顿-拉夫逊法。在运行方式管理中，潮流是确定电网运行方式的基本出发点；在规划领域，需要进行潮流分析验证规划方案的合理性；在实时运行环境，调度员潮流提供了多个在预想操作情况下电网的潮流分布以校验运行可靠性。在电力系统调度运行的多个领域都涉及到电网潮流计算。潮流是确定电力网络运行状态的基本因素，潮流问题是研究电力系统稳态问题的基础和前提。在用数字计算机解电力系统潮流问题的开始阶段，普遍采取以节点导纳矩阵为基础的逐次代入法。在解决电力系统潮流计算问题时，是以导纳矩阵为基础的，因此，只要我们能在迭代过程中尽可能保持方程式系数矩阵的稀疏性，就可以大大提高牛顿法潮流程序的效率。自从60年代中期，在牛顿法中利用了最佳顺序消去法以后，牛顿法在收敛性。内存要求。速度方面都超过了阻抗法，成为60年代末期以后广泛采用的优秀方法。 第二章 程序及其说明2.1 设计资料及参数清华大学课程设计任务书课程名称电力系统稳态分析课程设计设计题目两机五节点网络潮流计算牛拉法指导教师时间1周一、教学要求电力系统稳态分析课程设计以设计和优化电力系统的潮流分析为重点，提高学生综合能力为目标，尽可能结合实际工程进行。设计内容的安排要充分考虑学校现有的设备，设计时间及工程实际需要，并使学生初步学会运用所学知识解决工程中的实际问题。二、设计资料及参数（一）设计原始资料1、待设计电气设备系统图2、电力系统网络各元件参数3、电力系统电气元件的使用规范4、电力工程电气设计手册（二）设计参考资料1、电力系统稳态分析，陈珩，中国电力出版社，2007，第三版2、电力系统分析，韩祯祥，浙江大学出版社，2005，第三版3、电力系统分析课程实际设计与综合实验，祝书萍，中国电力出版社，2007，第一版三、设计要求及成果1.根据给定的参数或工程具体要求，收集和查阅资料；学习相关软件（软件自选）。2.在给定的电力网络上画出等值电路图。3.运用计算机进行潮流计算。4.编写设计说明书。基本要求：1.编写潮流计算程序；2.在计算机上调试通过（？）；3.运行程序并计算出正确结果（？）；4.写出课程设计报告（包括以下内容）（1份）（1）程序框图；（2）源程序；（3）符号说明表；（4）算例及计算结果5.编写计算说明书（1份）。四、进度安排根据给定的参数或工程具体要求，收集和查阅资料（半天）学习软件（MATLAB或C语言等）（一天半）编程计算复杂系统潮流计算（三天）编写计算设计书（一天）五、评分标准课程设计成绩评定依据包括以下几点：1） 工作态度（占10%）；2） 基本技能的掌握程度（占20%）；3） 程序编写是否合理是否有运行结果(40%)；4） 课程设计说明书编写水平(占30%)。5） 分为优、良、中、合格、不合格五个等级。考核方式：设计期间教师现场检查；评阅设计报告系统接线图其中节点1为平衡节点，节点2、3、4、5为PQ节点。2.1.1牛顿 拉夫逊法的程序框图2.1.2C语言编写潮流计算程序#include#include#define N 5 /节点数#define n_PQ 4 /PQ节点数#define n_PV 0 /PV节点数#define n_br 7 /串联支路数void main()void disp_matrix(float *disp_p,int disp_m,int disp_n); /矩阵显示函数double Us2*N=1.06,0,1.00,0,1.00,0,1.00,0,1.00,0; /电压初值double PsN=0,0.2,-0.45,-0.4,-0.6; /有功初值 double QsN=0,0.2,-0.15,-0.05,-0.1; /无功初值 float GNN,BNN; /各几点电导电纳struct /阻抗参数 int nl; /左节点 int nr; /右节点 float R; /串联电阻值 float X; /串联电抗值 float Bl; /左节点并联电导 float Br; /右节点并联电纳 ydatan_br=1, 2, 0.02 ,0.06 ,0 , 0 ,1, 3, 0.08 ,0.24 ,0 , 0 ,3, 4, 0.01 ,0.03 ,0 , 0 ,4, 5, 0.08 ,0.24 ,0 , 0 ,2, 3, 0.06 ,0.18 ,0 , 0 ,2, 5, 0.04 ,0.12 ,0 , 0 ,2, 4, 0.06 ,0.18 ,0 , 0 ;float Z2; /Z2=R2+X2 各串联阻抗值的平方float eN,fN,dfe2*N; /e,f存储电压的x轴分量和y轴分量,dfe存储电压修正值float mid1N,mid2N,dS2*N; /mid1、mid2存储计算雅克比行列式对角线元素的中间值,dS存储PQU的不平衡量float Jacob2*N2*N,inv_J2*N2*N; /雅克比行列式float dPQU=1.0; /PQU不平衡量最大值int kk=0; /迭代次数int i,j,k; float t; float Pijn_br; /存储线路i-j的有功float Qijn_br; /存储线路i-j的无功float Pjin_br; /存储线路j-i的有功float Qjin_br; /存储线路j-i的无功float dPijn_br; /存储线路i-j的有功损耗 float dQijn_br; /存储线路i-j的无功损耗float AA,BB,CC,DD; /存储线路潮流计算时的中间值float ppp=0,qqq=0;/形成导纳矩阵-for(i=0;iN;i+) for(j=0;jN;j+) Gij=0; Bij=0; for(i=0;in_br;i+) if(ydatai.nl!=ydatai.nr) Z2=(ydatai.R)*(ydatai.R)+(ydatai.X)*(ydatai.X); /串联阻抗等效导纳值 /非对角元素 Gydatai.nl-1ydatai.nr-1=(-ydatai.R)/Z2; Bydatai.nl-1ydatai.nr-1=ydatai.X/Z2; Gydatai.nr-1ydatai.nl-1=(-ydatai.R)/Z2; Bydatai.nr-1ydatai.nl-1=ydatai.X/Z2; /对角元素 Gydatai.nl-1ydatai.nl-1+=ydatai.R/Z2; Gydatai.nr-1ydatai.nr-1+=ydatai.R/Z2; Bydatai.nl-1ydatai.nl-1+=(-ydatai.X/Z2); Bydatai.nr-1ydatai.nr-1+=(-ydatai.X/Z2); /并联导纳等效导纳值 Bydatai.nl-1ydatai.nl-1+=ydatai.Bl; Bydatai.nr-1ydatai.nr-1+=ydatai.Br; else Gydatai.nl-1ydatai.nr-1+=ydatai.R; Bydatai.nl-1ydatai.nr-1+=ydatai.X; printf(G=n);disp_matrix(*G,N,N);printf(B=n);disp_matrix(*B,N,N);/分离e,ffor(i=0;i0.00001) /计算功率不平衡量 for(i=1;iN;i+) mid1i=0; mid2i=0; for(j=0;jN;j+) mid1i=mid1i+Gij*ej-Bij*fj; mid2i=mid2i+Gij*fj+Bij*ej; dS2*i=Psi-(ei*mid1i+fi*mid2i); if(in_PQ) dS2*i+1=Qsi-(fi*mid1i-ei*mid2i); else dS2*i+1=Us2*i*Us2*i-(ei*ei+fi*fi); dPQU=0; for(i=2;i2*N;i+) if(dSi0&dPQU0&dPQU0.00001) kk+; /形成雅克比行列式 for(i=2;i2*N;i+) for(j=2;j2*N;j+) Jacobij=0; for(j=1;jN;j+) /求H,N for(i=1;iN;i+) if(i!=j) Jacob2*i2*j=Bij*ei-Gij*fi; Jacob2*i2*j+1=-Gij*ei-Bij*fi; else Jacob2*i2*i=Bii*ei-Gii*fi-mid2i; Jacob2*i2*i+1=-Gij*ei-Bij*fi-mid1i; /求J,L for(i=1;i=n_PQ;i+) if(i!=j) Jacob2*i+12*j=Gij*ei+Bij*fi; Jacob2*i+12*j+1=Bij*ei-Gij*fi; else Jacob2*i+12*i=Gij*ei+Bij*fi-mid1i; Jacob2*i+12*i+1=Bij*ei-Gij*fi+mid2i; /求R,S for(i=n_PQ;iN-1;i+) if(i=j) Jacob2*i+12*i=-2*fi; Jacob2*i+12*i+1=-2*ei; /雅克比行列式求逆 for(i=2;i2*N;i+) for(j=2;j2*N;j+) if(i!=j) inv_Jij=0; else inv_Jij=1; for(i=2;i2*N;i+) for(j=2;j2*N;j+) if(i!=j) t=Jacobji/Jacobii; for(k=2;k2*N;k+) Jacobjk-=Jacobik*t; inv_Jjk-=inv_Jik*t; for(i=2;i2*N;i+) if(Jacobii!=1) t=Jacobii; for(j=2;j2*N;j+) inv_Jij=inv_Jij/t; /求电压修正值- for(i=2;i2*N;i+) dfei=0; for(j=2;j2*N;j+) dfei-=inv_Jij*dSj; for(i=1;iN;i+) ei+=dfe2*i+1; fi+=dfe2*i; else break;/循环结束-/求平衡节点功率-mid10=0;mid20=0;for(j=0;jN;j+) mid10=mid10+G0j*ej-B0j*fj; mid20=mid20+G0j*fj+B0j*ej;Ps0=e0*mid10+f0*mid20;Qs0=f0*mid10-e0*mid20;/ 显示输出结果printf(kk=%dn,kk);printf(P=);for(i=0;iN;i+) printf(%9.4f,Psi);printf(nQ=);for(i=0;iN;i+) printf(%9.4f,Qsi);for(i=0;iN;i+)ppp+=Psi;qqq+=Qsi;printf(nSS=);printf(%9.4f+j%9.4f,ppp,qqq);/输出线路总损耗功率printf(ne=);for(i=0;iN;i+) printf(%9.4f,ei);printf(nf=);for(i=0;iN;i+) printf(%9.4f,fi);printf(n);/求线路上的潮流/计算Sijfor(i=0;in_br;i+) if(ydatai.nl!=ydatai.nr) Z2=(ydatai.R)*(ydatai.R)+(ydatai.X)*(ydatai.X); AA=-fydatai.nl-1*ydatai.Bl+(eydatai.nl-1-eydatai.nr-1)*ydatai.R/Z2+(fydatai.nl-1-fydatai.nr-1)*ydatai.X/Z2; BB=-eydatai.nl-1*ydatai.Bl-(fydatai.nl-1-fydatai.nr-1)*ydatai.R/Z2+(eydatai.nl-1-eydatai.nr-1)*ydatai.X/Z2; Piji=eydatai.nl-1*AA-fydatai.nl-1*BB; Qiji=eydatai.nl-1*BB+fydatai.nl-1*AA; printf(S%d%d=%9.4f+j%9.4fn,ydatai.nl,ydatai.nr,Piji,Qiji); dPiji=Piji+Pjii; dQiji=Qiji+Qjii; p