电力系统课程设计潮流计算

《电力系统课程设计潮流计算》由会员分享，可在线阅读，更多相关《电力系统课程设计潮流计算（29页珍藏版）》请在装配图网上搜索。

1、word摘 要潮流计算是电力系统非常重要的分析计算，用以研究系统规划和运行中提出的各种问题。对规划中的电力系统，通过潮流计算可以检验所提出的电力系统规划方案能否满足各种运行方式的要求；对运行中的电力系统，通过潮流计算可以预知各种负荷变化和网络结构的改变会不会危与系统的安全，系统中所有母线的电压是否在允许的X围以内，系统中各种元件(线路、变压器等)是否会出现过负荷，以与可能出现过负荷时应事先采取哪些预防措施等。潮流计算是电力系统分析最根本的计算。除它自身的重要作用之外，潮流计算还是网损计算、静态安全分析、暂态稳定计算、小干扰静态稳定计算、短路计算、静态和动态等值计算的根底。实际电力系统的潮流计算

2、主要采用牛顿-拉夫逊法。按电压的不同表示方法，牛顿-拉夫逊潮流计算分为直角坐标形式和极坐标形式两种。本次计算采用直角坐标形式下的牛顿-拉夫逊法，牛顿-拉夫逊法有很好的收敛性，但要求有适宜的初值。传统的潮流计算程序缺乏图形用户界面，结果显示不直接难与其他分析功能集成。网络原始数据输入工作大量且易于出错。本文采用MATLAB语言运行WINDOWS操作系统的潮流计算软件。目前MATLAB已成为国际控制界最流行、使用最广泛的语言了。它的强大的矩阵处理功能给电力系统的分析、计算带来很多方便，而且采用MATLAB界面直观，运行稳定，计算准确。所以本次课程设计程序设计采用MATLAB计算。关键词： 电力系统

3、潮流计算 牛顿拉夫逊法潮流计算 MATLAB目 录一、 概述1.1 设计目的与要求.31.1.1 设计目的.31.1.2 设计要求.31.2 设计题目.31.3 设计内容.3二 电力系统潮流计算概述.42.1 电力系统简介.42.2 潮流计算简介.42.3 潮流计算的意义与其开展. .5三 潮流计算设计题目.63.1 潮流计算题目. 63.2 对课题的分析与求解思路. 7四 潮流计算算法与手工计算.74.1 极坐标下P-Q法的算法.74.2 节点电压方程.84.3节点导纳矩阵.94.4 导纳矩阵在潮流计算.104.5 潮流计算的手工计算.12五 Matlab概述.135.1 Matlab简介.

4、 145.2 Matlab的应用.145.3 矩阵的运算. 145.3.1 与常数的运算. 145.3.2 根本数学运算.145.3.3 逻辑关系运算.145.4 Matlab中的一些命令. 15六 潮流计算流程图与源程序.186.1 潮流计算流程图.186.2 潮流计算源程序图.196.3 运行计算结果.27七 总结.29八 参考文献.29第一章 系统概述1.1 设计目的与要求1.1.1 设计目的1. 掌握电力系统潮流计算的根本原理；2. 掌握并能熟练运用一门计算机语言MATLAB语言或C语言或C+语言；3. 采用计算机语言对潮流计算进展计算机编程。1.1.2 设计要求1. 程序源代码；2.

5、 给定题目的输入，输出文件；3. 程序说明；4. 给定系统的程序计算过程；5. 给定系统的手算过程至少迭代2次。1.2 设计题目电力系统潮流计算牛顿-拉夫逊法、P-Q分解法1.3 设计内容1. 根据电力系统网络推导电力网络数学模型，写出节点导纳矩阵；2. 赋予各节点电压变量直角坐标系形式初值后，求解不平衡量；3. 形成雅可比矩阵；4. 求解修正量后，重新修改初值，从2开始重新循环计算；5. 求解的电压变量达到所要求的精度时，再计算各支路功率分布、功率损耗和平衡节点功率；6. 上机编程调试；7. 计算分析给定系统潮流分析并与手工计算结果做比拟分析；8. 书写课程设计说明书。第二章 电力系统潮流计

6、算概述电力工业开展初期，电能是直接在用户附近的发电站或称发电厂中生产的，各发电站孤立运行。随着工农业生产和城市的开展，电能的需要量迅速增加，而热能资源和水能资源丰富的地区又往往远离用电比拟集中的城市和工矿区，为了解决这个矛盾，就需要在动力资源丰富的地区建立大型发电站，然后将电能远距离输送给电力用户。同时，为了提高供电的可靠性以与资源利用的综合经济性，又把许多分散的各种形式的发电站，通过送电线路和变电所联系起来。这种由发电机、升压和降压变电所，送电线路以与用电设备有机连接起来的整体，即称为电力系统。现代电力系统提出了“灵活交流输电和新型直流输电的概念。灵活交流输电技术是指运用固态电子器件与现代自

7、动控制技术对交流电网的电压、相位角、阻抗、功率以与电路的通断进展实时闭环控制，从而提高高压输电线路的诉讼能力和电力系统的稳态水平。新型直流输电技术是指应用现电力电子技术的最新成果，改善和简化变流站的造价等。运营方式管理中，潮流是确定电网运行方式的根本出发点：在规划领域，需要进展潮流分析验证规划方案的合理性；在实时运行环境，调度员潮流提供了电网在预想操作预想下的电网的潮流分布以与校验运行的可靠性。在电力系统调度运行的多个领域都涉与到电网潮流计算。潮流是确定电力网咯运行状态的根本因素，潮流问题是研究电力系统稳态问题的根底和前提。电力系统潮流计算是研究电力系统稳态运行情况的一种计算，它根据给定的运行

8、条件与系统接线情况确定整个电力系统各局部的运行状态：各母线的电压。各元件中流过的功率，系统的功率损耗等等。在电力系统规划的设计和现有电力系统运行方式的研究中，都需要利用潮流计算来定量的分析比拟供电方案或运行方式的合理性。可靠性和经济性。此外，电力系统的潮流计算也是计算机系统动态稳定和静态稳定的根底，所以潮流计算是研究电力系统的一种和重要和根底的计算。电力系统潮流计算也分为离线计算和在线计算两种，前者主要用于系统规划设计和安排系统的运行方式，后者如此用于正在运行系统的经常监视与实时控制。利用电子数字计算机进展潮流计算从50年代中期就已经开始了。在这20年内，潮流计算曾采用了各种不同的方法，这些方

9、法的开展主要围绕着对潮流计算的一些根本要求进展的，对潮流计算的要求可以归纳为以下几点：（1） 计算方法的可靠性或收敛性；（2） 对计算机内存量的要求；（3） 计算速度；（4） 计算的方便性和灵活性。电力系统潮流计算是电力系统分析中的一种最根本的计算，是对复杂电力系统正常和故障条件下稳态运行状态的计算。潮流计算的目标是求取电力系统在给定运行状态的计算，即节点电压和功率分布，用以检查系统各元件是否过负荷。各点电压是否满足要求，功率的分布和分配是否合理以与功率损耗等。对现有的电力系统的运行和扩建，对新的电力系统进展规划设计以与对电力系统进展静态和稳态分析都是以潮流计算为根底。潮流计算结果可用如电力系

10、统稳态研究，安全估计或最优潮流等对潮流计算的模型和方法有直接影响。实际电力系统的潮流技术那主要采用牛顿拉夫逊法。运行方式管理中，潮流是确定电网运行方式的根本出发点；在规划领域，需要进展潮流分析验证规划方案的合理性；在实时运行环境，调度员潮流提供了多个在预想操作情况下电网的潮流分布以与校验运行可靠性。在电力系统调度运行的多个领域问题是研究电力系统稳态问题的根底和前提。在用数字解算计算机解电力系统潮流问题的开始阶段，普遍采取以节点导纳矩阵为根底的逐次代入法。这个方法的原理比拟简单，要求的数字计算机内存量比拟差下，适应50年代电子计算机制造水平和当时电力系统理论水平，但它的收敛性较差，当系统规模变大

11、时，迭代次数急剧上升，在计算中往往出现迭代不收敛的情况。这就迫使电力系统的计算人员转向以阻抗矩阵为根底的逐次代入法。阻抗法改善了系统潮流计算问题的收敛性，解决了导纳无法求解的一些系统的潮流计算，在60年代获得了广泛的应用，阻抗法德主要缺点是占用计算机内存大，每次迭代的计算量大。当系统不断扩大时，这些缺点就更加突出，为了克制这些缺点，60年代中期开展了以阻抗矩阵为根底的分块阻抗法。这个方法把一个大系统分割为几个小的地区系统，在计算机内只需要存储各个地区系统的阻抗矩阵与它们之间联络的阻抗，这样不仅大幅度的节省了内存容量，同时也提高了计算速度。克制阻抗法缺点是另一个途径是采用牛顿-拉夫逊法。这是数学

12、中解决非线性方程式的典型方法，有较好的收敛性。在解决电力系统潮流计算问题时，是以导纳矩阵为根底的，因此，只要我们能在迭代过程中尽可能保持方程式系数矩阵的稀疏性，就可以大大提高牛顿法潮流程序的效率。自从60年代中期，牛顿法中利用了最优顺序消去法以后，牛顿法在收敛性。内存要求。速度方面都超过了阻抗法，成为了60年代末期以后广泛采用的优秀方法。 第三章 潮流计算设计题目3.1 潮流计算课题题目：在图1所示的简单电力系统中，系统中节点1、2为节点，节点3为节点，节点4为平衡节点，已给定，网络各元件参数的标幺值如表2所示，给定电压的初始值如表2所示，收敛系数。试求：图1 简单电力系统表1 网络各元件参数

13、的标幺值支路电阻电抗输电线路变压器变比k1213232434表2 各节点电压初值标幺值参数节点i12343.2 对课题的分析与求解思路此电力系统是一个4节点，5支路的电力网络。综合比拟牛顿拉夫逊法直角坐标、极坐标、PQ分解法等多种求解方法的特点，最后确定采用牛顿拉夫逊法极坐标。因为此方法所需解的方程组最少。第四章 潮流计算算法与手工计算4.1 极坐标下P-Q法的算法4.1.1 节点导纳矩阵Y根据题目提供的各节点的参数，求得节点导纳矩阵=4.1.2 简化雅可比矩阵B/和B/通过上一步的导纳矩阵，形成有功迭代和无功迭代的简化雅可比矩阵B/和B/对雅可比矩阵进展三角分解，形成因子表，为后面进展修正方

14、程计算作好准备。4.1.3 修正和迭代第一步，给定PQ节点初值和各节点电压相角初值。 第二步，作第一次有功迭代，按公式计算节点有功功率不平衡量。第三步，做第一次无功迭代，按公式计算无功功率不平衡量，计算时电压相角最新的修正值。解修正方程式，可得各节点电压幅值的修正量。第四步，第一轮有功迭代和无功迭代便做完了。第五步，按公式计算平衡节点功率。直到节点不平衡功率下降到10-5以下，迭代便可以完毕。4.2 潮流计算算法此题采用了题目要求的牛顿拉夫逊潮流计算的方法。牛顿-拉夫逊法潮流计算的公式。把牛顿法用于潮流计算，采用极坐标形式表示的如式1-3所示的形式。其中电压和支路导纳可表示为： 将上述表示式1

15、-2代入1-1式的右端，展开并分出实部和虚部，便得： 1-3按照以上的分类，PQ节点的输出有功功率和无功功率是给定的，如此第i节点的给定功率设为和称为注入功率。 假定系统中的第1、2、m节点为PQ节点，对其中每一个节点的N-R法表达式F(x)=0如、形式有些如下方程：1-4=1、2、m PV节点的有功功率和节点电压幅值是给定的。假定系统中的第m+1、m+2、n-1节点为PV节点，如此对其中每一PV节点可以列写方程： 1-5=m+1、m+2、n-1(6)形成雅可比矩阵。N-R法的思想是；本例；对F(x)求偏导的式1-6、式1-7，即式1-4、式1-5中的、是多维变量的函数，对多维变量求偏导、，并

16、以矩阵的形式表达称为雅可比矩阵。当j=i时，对角元素为 1-6当时,矩阵非对角元素为： 1-7 由上式不难看出，雅可比矩阵有以下特点。 雅可比矩阵中的诸元素都是节点电压的函数，因此在迭代过程中，它们将随着节点电压的变化而不断的变化。 雅可比矩阵具有结构对称性，数据不对称。如非对角，。 由式1-7可以看出，当导纳矩阵中非对角元素为零时，。雅可比矩阵中相应的元素也为零，即矩阵是非常稀疏的。因此，修正方程的求解同样可以应用稀疏矩阵的求解技巧。正是由于这一点才使N-R法获得广泛的应用。4.3 手工计算4.3.1 节点导纳矩阵求得节点导纳矩阵Y=各节点的导纳值如下：;Y11= Y13=-6.8966 +

17、17.2414i Y14=0Y41=04.3.2 简化雅可比矩阵 迭代中的雅克比矩阵:4.3.3 修正、迭代给定PQ节点初值和各节点电压相角初值V1=1.0，V2(0)=V3(0)=1.0，V42(0)=3(0)=0， 4(0)=01 作第一次有功迭代，按公式计算节点有功功率不平衡量 迭代中的P：2做第一次无功迭代，按公式计算无功功率不平衡量，计算时电压相角最新的修正值。 迭代中的Q： 解修正方程式，可得各节点电压幅值的修正量为 迭代中电压的模: 到这里为止，第一轮有功迭代和无功迭代便做完了。3 按公式计算平衡节点功率，得： P1+jQ1经过四轮迭代，节点不平衡功率也下降到10-5以下，迭代到

18、此完毕。 全线路各个点的功率分配如下： 0 -0.0582 - 0.3378i -0.3418 + 0.0172i 0 0 0.2654 + 0.9332i 0.0505 + 1.2558i 0 第五章 Matlab概述5.1 Matlab简介 MATLAB是由美国mathworks公司发布的主要面对科学计算、可视化以与交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以与非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以与必须进展有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言如C

19、、Fortran的编辑模式，代表了当今国际科学计算软件的先进水平MATLAB和Mathematica、Maple并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进展矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。5.2 Matlab的应用MATLAB的根本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用MATLAB来解算问题要比用C，FORTRAN等语言完成一样的事情简捷得多，并且mathwork也吸收了像Maple等软

20、件的优点,使MATLAB成为一个强大的数学软件。在新的版本中也参加了对C，FORTRAN，C+ ，JAVA的支持。可以直接调用,用户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用，此外许多的MATLAB爱好者都编写了一些经典的程序，用户可以直接进展下载就可以用。MALAB 产品族可以用来进展以下各种工作： 数值分析 数值和符号计算 工程与科学绘图 控制系统的设计与仿真 数字图像处理 技术 数字信号处理 技术 通讯系统设计与仿真 财务与金融工程 MATLAB 的应用X围非常广，包括信号和图像处理、通讯、控制系统设计、测试和测量、财务建模和分析以与计算生物学等众多应用领域。附

21、加的工具箱单独提供的专用 MATLAB 函数集扩展了 MATLAB 环境，以解决这些应用领域内特定类型的问题。 常数与矩阵的运算即是同该矩阵的每一元素进展计算。但需注意进展数除时，常数通常只能做除量。 根本函数运算中，矩阵的函数运算是矩阵预算中最实用的局部，常用的主要有以下几个： det(a) 求矩阵a的行列式 eig(a) 求矩阵a的特征值 inv(a)或a(-1) 求矩阵a的逆矩阵 rank(a) 求矩阵a的秩 trace(a) 求矩阵a的迹对角线元素之和 我们进展工程计算时常常遇到矩阵对应元素之间的运算。这种运算不同于前面讲的数学运算，为有所区别，我们称之为数组运算。 数组的加、减与矩阵

22、的加、减运算完全一样。而乘除法运算有相当大的区别，数组的乘除法是指两同维数组对应元素之间的乘除法，它们的运算符为“.*和“./或“.。前面讲过常数与矩阵的除法运算中常数只能做除数。在数组运算中有了“对应关系的规定，数组与常数之间的除法运算没有任何限制。 另外，矩阵的数组运算中还有幂运算运算符.、指数运算exp、对数运算(log)、和开方运算sqrt、等，有了“对应元素的规定，数组的运算实质上就是针对数组内部的每个元素进展的。矩阵的幂运算与数组的幂运算有很大的区别。 5.5 逻辑关系运算 逻辑运算是MATLAB中数组运算所特有的一种运算形式，也是几乎所有的高级语言普遍适用的一种运算。5.6 Ma

23、tlab中的一些命令1一般MATLAB命令格式为 输出参数1，输出参数2，=命令名输入参数1，输入参数2， 输出参数用方括号，输入参数用圆括号如果输出参数只有一个可不使用括号。 2可用、键来重现已输入的数据或命令。用、键来移动光标进展修改。3所有MATLAB命令都用小写字母。大写字母和小写字母分别表示不同的变量。4常用数有特定的名字，如pi=3.141596 、Inf=、NaN如此表示不定型求得的结果如0/0。5矩阵的输入要一行一行的进展，每行各元素用空格或，分开，每行用；分开。 6MATLAB书写格式为A=1 2 3 ；4 5 6 ；7 8 9 在MATLAB中运行如下程序可得到A矩阵 a=

24、1 2 3;4 5 6;7 8 9 a =1 2 3 4 5 6 7 8 97需要显示命令的计算结果时，如此语句后面不加“；号，否如此要加“；号。 运行下面两种格式可以看出他们的区别 a=1 2 3;4 5 6;7 8 9 a=1 2 3;4 5 6;7 8 9; a = 不显示计算结果 1 2 3 4 5 68当输入语句过长需要换行时，应加上“后再回车，如此可连续输入。9diary命令使用该命令可以在窗口中以ASCII码形式记录所有的输入和输出。但这个命令不是存储数据，而是存储输入与屏幕上输出的内容。它可以记录下工作的过程。在每个工作过程之前使用该命令，工作完毕后使用diary off 如此

25、能将整个工作过程记录下来。 格式diary(文件名)扩展名. diary off一般来说扩展名可取，m这样就可在MATLABBIN目录下存入该文件。 10save 命令该命令存储定义的变量或演算结果，也可以用来存储指定的变量。 命令格式为 save 文件名 . 扩展名 11what 命令 该命令可以在当目录下显示MATLAB文件和MAT数据文件12dir 命令 显示当前目录下的所有文件.13clear命令 14d1,d2,d3,.=size(a) 求矩阵的大小，对m*n二维矩阵，第一个为行数m,第二个为列数n。如果输入 calear a b c ，如此表示去除工作空间中指定变量a,b,c；如果

26、仅仅输入calear命令，如此去除整个工作空间。 与此同时，MATLAB具有强大的矩阵运算功能, 但由于我们在求节点导纳矩阵时用的不多, 因此这里我们只作简单介绍。1) 在MATLAB中表示一个矢量要用方括号， 而列矢量的输入只需在行矢量输入格式根底上加转置符()即可。 如 x=1 2 3;4 5 6 x = 1 2 3 4 5 6 而 x=1 2 3;4 5 6 (加转置符) x = 1 4 2 5 3 6注意上面两式的区别。2) 下面三条命令可以产生一个行矢量1 a=linspace(x,y,n) 2 a=logspace(x,y,n) 3 a=x:n:y第一条命令可以在线性空间产生一个值

27、在10x至10y之间间隔点数为n的行矢量(一组数据)。第二条命令可以在对数空间产生一个值在x至y之间等间隔的行矢量(一组数据)。其行矢量的起始值是x, 终值为y, 点数为n。第三条命令产生X至y步长为n的行矢量。但是, 三个命令之间存在差异，下面的例子可以说明这一点。例一 x=logspace(0,5,6) x = 1 10 100 1000 10000 100000例二 x=linspace(0,10,11) x = 0 1 2 3 4 5 6 7 8 9 10例三 x=0:1:10 x = 0 1 2 3 4 5 6 7 8 9 10 通过上面三个例子可以看出例一, 例二中n代表选取的点数

28、。而在例三 中n如此表示步长. 我们应当注意它们的区别。3) 矩阵的加, 减, 乘, 除等, 和其它语言书写一样。但要注意的是在运算符 前面加有(.)如此表示是元素对元素的操作. 4)以下是常用的运算命令运算命令名 功能Angle求复数的角Min求最小值Max求最大值Sum求和Roots求多项式的根Poly由多项式的根求多项式的系数Polyval求给定点多项式的值Polyder多项式求导在进展潮流分布计算时，实际上是由多个简单系统构成的复杂系统，在求节点导纳矩阵时的要用到反应的指令，所以在MATLAB中有下面几种命令可以解决两个系统间的连接问题。1系统的并联parallel 命令可以实现两个系

29、统的并联。示意图如下：系统1系统2 u1 y1u + y u2 y2并联后的系统传递函数表示式为：其中n1、d1和n2、d2分别为g1s、g2s的传递函数分子、分母系数行矢量。命令格式： n，d = paralltln1，d1，n2，d2a，b，c，d = paralltla1，b1，c1，d1，a2，b2，c2，d22系统的串联series命令实现两个系统的串联，示意图如下：g2sg1s u1 y u2 y2串联后系统的传递函数为 命令格式： n，d = seriesn1，d1，n2，d2 a，b，c，d = seriesa1，b1，c1，d1，a2，b2，c2，d23系统的反应feedba

30、ck命令实现两个系统的反应连接，示意图如下：u1 + y1 g2sg1s y2 u2连接后系统的传递函数表示为： 命令格式： n，d = feedbackn1，d1，n2，d2或： n，d = feedbackn1，d1，n2，d2，sign a，b，c，d = feedbacka1，b1，c1，d1，a2，b2，c2，d2，sign其中sign是指示y2到u1连接的符号，缺省时默认为负即sign = -1。4系统的闭环cloop命令可以将系统的输出反应到系统的输入构成闭环系统，示意图如下：g1su y正、负反应后闭环系统为： 命令格式： n，d = cloop n1，d1，sign ac，b

31、c，cc，dc = cloopa，b，c，d，sign通过以上对MATLAB根本指令的了解，我们就可以对所求的电力系统网络的节点导纳矩阵进展画编程框架图。第六章 潮流计算流程图与源程序本次课程设计采用极坐标下的牛顿-拉夫逊计算网络的潮流计算。其牛顿-拉夫逊潮流计算程序框图如下所示。输入原始数据形成节点导纳矩阵给定节点电压初值ei(0),fi(0)用公式计算不平衡功率Piki QikiVi2kimax|PiKiQiiiVi2ki|pr for m=1:n2 UD(m)=U1(m); end for m=1:n1+n2 for R=1:n PT(R)=U1(m)*U1(R)*(G(m,R)*cos

32、(O(m)-O(R)+B(m,R)*sin(O(m)-O(R); end PT1(m)=sum(PT); PP(m)=P(m)-PT1(m); PP1(k+1,m)=PP(m); end for m=1:n2 for R=1:n QT(R)=U1(m)*U1(R)*(G(m,R)*sin(O(m)-O(R)-B(m,R)*cos(O(m)-O(R); end QT1(m)=sum(QT); QQ(m)=Q(m)-QT1(m); QQ1(k+1,m)=QQ(m); endPR1=max(abs(PP); PR2=max(abs(QQ); PR=max(PR1,PR2); for m=1:n1+n

33、2 for R=1:n1+n2 if m=R H(m,m)=U1(m)2*B(m,m)+QT1(m); else H(m,R)=-U1(m)*U1(R)*(G(m,R)*sin(O(m)-O(R)-B(m,R)*cos(O(m)-O(R); end end end for m=1:n1+n2 for R=1:n2 if m=R N(m,m)=-U1(m)2*G(m,m)-PT1(m); else N(m,R)=-U1(m)*U1(R)*(G(m,R)*cos(O(m)-O(R)+B(m,R)*sin(O(m)-O(R); end end end for m=1:n2 for R=1:n1+n2

34、 if m=R J(m,m)=U1(m)2*G(m,m)-PT1(m); else J(m,R)=U1(m)*U1(R)*(G(m,R)*cos(O(m)-O(R)+B(m,R)*sin(O(m)-O(R); end end end for m=1:n2 for R=1:n2 if m=R L(m,m)=U1(m)2*B(m,m)-QT1(m); else L(m,R)=-U1(m)*U1(R)*(G(m,R)*sin(O(m)-O(R)-B(m,R)*cos(O(m)-O(R); end end end JJ=H N;J L; disp(JJ); PQ=PP;QQ; DA=-inv(JJ)*

35、PQ; DA1=DA; for m=1:n1+n2 OO(m)=DA1(m); end for m=n:n1+n2+n2 UU1(m-n1-n2)=DA1(m); end UD2=diag(UD); UU=UU1*UD2; for m=1:n1+n2 O(m)=O(m)+OO(m); end for m=1:n2 U1(m)=U1(m)+UU(m); end for m=1:n1+n2 o(k+1,m)=180/pi*O(m); end for m=1:n2 u(k+1,m)=U1(m); end k=k+1;endfor m=1:n b(m)=U1(m)*cos(O(m); c(m)=U1(

36、m)*sin(O(m);endU=b+i*c;for R=1:n PH1(R)=U(isb)*conj(Y(isb,R)*conj(U(R);endPH=sum(PH1);for m=1:n for R=1:n if m=R C1(m,R)=1/C(m,R); else C1(m,m)=C(m,m); end endendfor m=1:n for R=1:n if (C(m,R)=inf)&(m=R) SS(m,R)=U1(m)2*conj(C1(m,m)+U(m)*(conj(U(m)-conj(U(R)*conj(C1(m,R); end endenddisp(迭代中的P：);disp(

37、PP1);disp(迭代中的Q：);disp(QQ1);disp(迭代中相角:);disp(o);disp(迭代中电压的模:);disp(u);disp(平衡结点的功率:);disp(PH);disp(全部线路功率分布:);disp(SS);6.2.2 程序结果请输入结点数：n=4请输入PV结点数：n1=1请输入PQ结点数：n2=2请输入平衡结点：isb=4请输入变比矩阵看:K=0 0 0 0;0 0 0 0.9625;0 0 0 0;0 0 0 0请输入支路阻抗矩阵：C=0 0.03+0.09i 0.02+0.05i inf; 0.03+0.09i 0 0.04+0.08i 0.0+0.05

38、i;0.02+0.05i 0.04+0.08i 0 0.03+0.07i;inf 0.0+0.05i 0.03+0.07i 0请输入支路导纳矩阵：y=0 0.01i 0.01i 0;0.01i 0 0 0;0.01i 0 0 0;0 0 0 0请输入结点电压矩阵：U=1+0i 1+0i 1.02+0i 1.05+0i请输入各结点的功率：S=-0.4-0.3i -0.3-0.2i 0.4 0结点导纳矩阵: 10.2299 -27.2214i -3.3333 +10.0000i -6.8966 +17.2414i 0 迭代中的雅克比矩阵:迭代中的P：迭代中的Q：迭代中相角:迭代中电压的模:平衡结点

39、的功率:全部线路功率分布: 0 -0.0582 - 0.3378i -0.3418 + 0.0172i 0 0 0.2654 + 0.9332i 0.0505 + 1.2558i 0 第七章 总 结 此次课程设计首先让我明白了要使电力系统运行的稳定，必须经过精细的设计和计算。在进展课题设计的过程中，加深了我对潮流计算的认识，尤其是对牛顿拉夫逊潮流计算的求解思路有了比拟透彻的理解。同时由于求解过程中用到求节点导钠矩阵，求矩阵的逆等等，又使我对以前所学的知识有了一次很好的温习。同时也看到了研究性学习的效果，从研究中去学习，理论结合实际，将理论运用到实际，同时在实践中发现问题，然后解决问题， 通过这次课程设计我在MATLAB编程，潮流计算，WORD文档的编辑方面均有所提高，但也暴露出了一些问题；我发现了自己的根底知识有很多的不足，对MATLAB的性能和特点还不能有一个全面的把握，这些根底的缺乏给我的设计计划造成了不小的障碍，相信通过以后的学习能弥补这些不足，从而达到一个新的层次。 第八章 参考文献1 何仰赞等.电力系统分析M. 2 某某交通大学等.电力系统计算M.3 陈衍.电力系统稳态分析M4 李光琦.电力系统暂态分析M.5 于永源，杨绮雯. 电力系统分析第二版M