配电网的潮流计算

1、 毕业设计（论文）题，目 配电网潮流计算与程序设计 学生姓名 石昊晨 学号 2010151107 专业 发电厂及电力系统 班级 20109091 指导教师 刘会家 评阅教师 完成日期 年月日 目 录 摘要一配电网潮流概述41.1配电网潮流计算的目的与意义51.2潮流计算方法概述51.2.1 牛顿拉夫逊法61.2.2 快速解耦法71.2.3 回路阻抗法91.2.4 前推回代法111.3 本文工作11二配电网网络模型112.1元件模型112.1.1 电力线路的数学模型112.1.2 变压器的等值电路132.2网络模型15三：基于matlab的配电网潮流计算算法153.1配电网潮流计算算法原理153

2、.2 matlab的概述193.3程序设计203.3.1 牛顿-拉夫逊法潮流求解过程203.3.2牛顿拉夫逊法的程序框图24四：算例26参考文献27致 谢28配电网潮流计算与程序设计 学 生：石昊晨指导教师：刘会家 （三峡大学 国际文化交流学院）摘 要：本文首先分析了配电网的特点及对算法的要求，然后建立配电网潮流计算模型。针对配电网潮流计算的现状进行了全面分析，深入讨论了目前各方法的特点，并从收敛性及其他性能指标进行了比较分析；详细研究用的比较广泛的牛顿拉夫逊法，并以广度优先顺序搜索策略作为理论基础。针对某地区配电网的具体情况，选取IOKV的配电网子系统进行潮流计算。利用MATLAB 2009

3、a 进行了基于牛顿拉夫逊法的配电网的潮流计算程序。由计算结果可知，该算法具有一定的优越性，软件的开发具有一定的实用性。关键词：电力系统，配电网潮流，牛顿拉夫逊法，MATLAB程序设计Abstract：In this paper, ungrounded system, the characteristics of non-zero sequence path, a three-phase decoupled power flow calculation method. This method ignores the influence of zero sequence components, m

4、aking the three-phase asymmetrical load caused by phase coupling decoupling to be achieved by the phase flow calculation. The algorithm flow algorithm to the existing distribution network in the three-phase node voltage equation 3n-order decomposition of the node voltage equation of three n-order, s

5、o no matter what kind of algorithm can greatly save memory and computation for the distribution network to achieve by phase analysis provides a good way. In this paper, a system of 36 nodes to verify the results show that the method can fully into account the impact of unbalanced three-phase loads,

6、a better computational speed and accuracy.Keywords: power systems, phase decoupling, power flow, back/forward sweep algorithm一 电力系统潮流概述1.1配电网潮流计算的目的与意义电力系统潮流计算是研究电力系统稳态运行情况的一种基本电气计算，是电力系统规划和运营中不可缺少的一个重要组成部分。可以说，潮流计算是电力系统分析中最基本、最重要的计算，是电力系统安全、经济分析和实时控制与调度的基础。潮流计算的任务是根据给定的运行条件和网络结构确定整个系统的运行状态，如各母线上的电压

7、（幅值及相角）、网络中的功率分布及功率损耗等。即潮流计算是对电力系统的功率分布和电压分布的计算, 其具体任务就是编制系统的调度计划和电气设备检修计划, 确定电力系统中变压器分接头位置和系统中枢点与电压控制点的电压曲线， 进行事故运行方式的分析, 为电力系统短路和稳定计算提供数据, 为继电保护及自动装置整定与电力系统设计和规划提供依据等。潮流计算的目的是对现有电力系统的正常运行状态进行分析, 以提示必要的改进措施, 同时为新建系统或扩建系统的有关分析、计算打下基础。配电网潮流计算是配电网经济运行、系统分析等的重要基础，但由于配电网与输电网有着明显的差异:配电网具有环形结构, 而通常以开环方式运行

8、。通常呈辐射状，支路比值较大，分支线较多；配电线路中的R/X 比值偏大使输电网中常用的潮流计算算法如传统的牛顿法和快速分解法在应用于配电网潮流计算时容易形成病态而无法收敛，因此，研究适合于配电网的潮流算法也是至关重要的。目前，输电系统潮流计算方法已较为成熟，而且获得了广泛的实际应用。但随着电力系统规模的不断扩大，潮流方程的阶数越来越高，对这种规模的方程并不是采用任何数学方法都能保证给出正确答案的，因此，这也成为促使电力系统研究人员不断寻求新的、更可靠的潮流计算方法的动力。随着现代电力系统大系统、强非线性与多元件的特点日益突出, 其计算量与计算复杂度急剧增加。旧的计算机软件在处理潮流计算时, 其

9、速度已无法满足大电网模拟和实时控制的仿真要求, 而高效的潮流问题的相关软件的研究已成为大规模电力系统仿真计算的关键。 1.2潮流计算方法概述与输电网相比，配电网的网络结构有着明显的差异：配电网的网络呈现辐射状，在正常运行是开环的，只有在倒换负荷或发生故障时才有可能出现短时环网运行或多电源运行的情况；配电线路的总长度较输电网络要长且分支较多，配电线的线径比输电线细，导致配电网的 R/ X较大，无法满足 << 的 PQ 解耦条件，所以在输电网中常用的快速解耦算法在配电网中难以收敛；由于配电网络直接面向用户，所以网络节点众多。八十年代中期到九十年代中期，随着国际国内电力企业对配电网管理的

10、重视程度的不断加深，对配电网潮流的研究也广泛开展起来，这期间出现了众多结合配电网特殊结构而开发的简单迭代算法。从模型求解过程上可分为牛顿拉夫逊法、快速解耦法、回路阻抗法和前推回代法。1.2.1 牛顿拉夫逊法电力系统潮流计算是电力系统分析中的一种最基本的计算，是对复杂电力系统正常和故障条件下稳态运行状态的计算。潮流计算的目标是求取电力系统在给定运行状态的计算。即节点电压和功率分布，用以检查系统各元件是否过负荷。各点电压是否满足要求，功率的分布和分配是否合理以及功率损耗等。对现有电力系统的运行和扩建，对新的电力系统进行规划设计以及对电力系统进行静态和暂态稳定分析都是以潮流计算为基础。潮流计算结果可

11、用如电力系统稳态研究，安全估计或最优潮流等对潮流计算的模型和方法有直接影响。实际电力系统的潮流技术那主要采用牛顿-拉夫逊法。牛顿-拉夫逊法(简称牛顿法)在数学上是求解非线性代数方程式的有效方法。其要点是把非线性方程式的求解过程变成反复地对相应的线性方程式进行求解的过程。即通常所称的逐次线性化过程。对于非线性代数方程组： 即 (3-1)在待求量x的某一个初始估计值附近，将上式展开成泰勒级数并略去二阶及以上的高阶项，得到如下的经线性化的方程组： (3-2)上式称之为牛顿法的修正方程式。由此可以求得第一次迭代的修正量 (3-3)将和相加，得到变量的第一次改进值。接着就从出发，重复上述计算过程。因此从

12、一定的初值出发，应用牛顿法求解的迭代格式为： (3-4) (3-5)上两式中：是函数对于变量x的一阶偏导数矩阵，即雅可比矩阵J；k为迭代次数。有上式可见，牛顿法的核心便是反复形式并求解修正方程式。牛顿法当初始估计值和方程的精确解足够接近时，收敛速度非常快，具有平方收敛特性。牛顿潮流算法突出的优点是收敛速度快，若选择到一个较好的初值，算法将具有平方收敛特性，一般迭代45次便可以收敛到一个非常精确的解。而且其迭代次数与所计算网络的规模基本无关。牛顿法也具有良好的收敛可靠性，对于对以节点导纳矩阵为基础的高斯法呈病态的系统，牛顿法也能可靠收敛。牛顿法所需的内存量及每次迭代所需时间均较高斯法多。牛顿法的

13、可靠收敛取决于有一个良好的启动初值。如果初值选择不当，算法有可能根本不收敛或收敛到一个无法运行的节点上。对于正常运行的系统，各节点电压一般均在额定值附近，偏移不会太大，并且各节点间的相位角差也不大，所以对各节点可以采用统一的电压初值(也称为平直电压)，如假定： 或 (3-6) 这样一般能得到满意的结果。但若系统因无功紧张或其它原因导致电压质量很差或有重载线路而节点间角差很大时，仍用上述初始电压就有可能出现问题。解决这个问题的办法可以用高斯法迭代12次，以此迭代结果作为牛顿法的初值。也可以先用直流法潮流求解一次以求得一个较好的角度初值，然后转入牛顿法迭代。1.2.2 快速解耦法为了改进牛顿法在内

14、容占用量及计算速度方面的不足，早在 1974 年有人提出的快速解耦法(对称 P-Q 分解法)是较成功的一种算法；它是密切结合高压电力系统固有特点，对牛顿法改进后得到的一种方法。原理是根据系统有功决定于电压相角的变化，而无功主要决定于电压模值的变化这一特性，并进行合理假设：(1)线路两端的相角差不大，且，即认为 ；(2)与节点无功功率对应的导纳远小于节点的自导纳，即。最后得修正方程式： （215）式中：、是由节点导纳矩阵的虚部构成的常数对称矩阵，可有 XB、BX 等方案。这种方法具有简单、快速、内存节省且收敛可靠的优点，是广泛应用于高压网在线处理计算的方法。该方法存在的问题是 R/ X比值敏感，

15、用于配电网可能迭代次数过多或不收敛。针对这一问题，提出了一种改进的快速解耦法。该方法的特点是，它根据配电网的辐射型特点，从一种新概念上构造出潮流方程，即前一节点的电压电流用含后一节点的电压和电流的关系式表示，即 （216）其中，为第 k 个节点的电压和对应的支路电流矩阵，为前后两个节点的关系方程。根据边界条件 ，可建立潮流方程如下所示： （217）其中，为按(2-16)从末端递推到始端形成的以末端电压为变量的方程，的雅可比矩阵可以表示为从馈线末端到始端所有支路雅可比矩阵的乘积，即 （218）其中， （219）这样，一方面可以减少方程的数目，使之等于支路数；另一方面能够充分利用配电网的辐射型结构

16、导致的数值特性，将雅可比矩阵简化为一个三角矩阵，使其求解的实质变为一种前推回推算法，从而简化了运算，并极大提高了其收敛性能。文献6通过以下假设将式(2-19)中的化为单位矩阵：(1)节点K电压的微小变化，将引起前一节点几乎相同的变化，因此左上角项的所有元素近似为1；(2)电流 的微小变化对影响很小，因此右上角项的所有元素近似为0；(3)的微小变化对影响很小，因此左下角项的所有元素近似为0；(4)电流的微小变化时，将引起几乎相同的变化，因此左下角项的所有元素近似为1。1.2.3 回路阻抗法在一般电力系统(发、输电网络)中，各节点和大地间有发电机、负荷、线路电容等对地支路，节点和节点间也有输电线路

17、和变压器支路，使得系统的节点方程式数小于回路方程式数。因而，一般电力系统的分析计算采用节点电压方程为宜。但对于低电压配电网络，由于一般不计配电线路对地充电电容的影响，并忽略变压器的对地导纳，网络中树支数将总大于连支数，因而适合采用回路电流方程进行分析。因此提出了一种基于回路方程的潮流算法，并称之为直接解(Direct SolutionMethod)。由于它基于回路阻抗方程，称之为回路阻抗法。该方法将各节点的负荷用恒定阻抗表示，从馈线节点到每一个负荷节点形成一条回路，以回路电流为变量，根据基尔霍夫电压定律，可列出回路电流方程式组： (220)式中，为根节点电压，为第 i 条回路上的回路电流(等于

18、负荷节点 i 的负荷电流)，为第 i 条回路的自阻抗(等于节点 i 与根节点 s 之间的支路阻抗和，加上节点 i 的负荷阻抗)，为第 i 条回路和第 j 条回路的互阻抗(等于节点 i 与节点 j 到根节点 s 的共同支路阻抗和)。设负荷节点数为 L，则回路阻抗矩阵 Z 是一个 L×L 维的不含零元素的方阵。采用 LU 分解方法对方程式(2-20)进行分解，可求出回路电流，也就得到各个负荷节点的负荷电流。然后可求出各条支路上的电压降，进而可求得各节点的电压和负荷节点的功率，反复迭代，直到求得的负荷节点功率与给定负荷的差值满足一定的精度要求为止。在回路阻抗阵中有许多相同的元素，实际上只有

19、网络支路数目个不同元素。但是在一般的编号方式下，这些不同的元素交叉混杂，无规律性可言。为了减少占用计算机的存储容量，文献8采用了一种特别的节点和支路编号方案，在这种编号方案下，回路阻抗矩阵 Z 和它三角分解得到的上三角矩阵 U 中的元素能够有规律地排列，即许多相同的元素集中排列在一起，因而可以借用“稀疏存储”技术，只存储其中不同的元素，只是这种编号方案太复杂而不易实现。在求 U 矩阵的元素时，文献6也通过采用一些求解技巧，提高了计算速度。但这些技巧不适用于在 U 矩阵中占很大比例的对角元素和同一行与它紧相邻的元素，因而限制了求解速度的提高。特别地，回路阻抗法处理网孔的能力较强，它对增加一条环路

20、后的处理方法比较简单：假定连接节点和（<）形成一条环路，则回路阻抗阵中将只有下面有限几个元素发生变化：(1)节点的自阻抗和节点的自阻抗；(2)和节点的互阻抗。因此只需对回路阻抗阵中的这几个元素进行修改即可。只是由于的改变，将可能在 U 阵的第列的第到第-1行产生-个“注入元素”，使系统的存储容量稍有增加。回路阻抗法中对已有环路的处理方法是，将环路在环路上 i 节点(设 i 节点的负荷为，电压为 )处分解为和节点，使节点和各连有值为的负荷阻抗，这样形成一个等值辐射网。求得这一辐射网的回路阻抗阵，并对矩阵元素进行修正，只需休整元素和即可，设其修正值分别为和。则 （221）由此可见，回路阻抗法

21、处理环路非常简单，处理弱环网的能力较强，因而有特别的应用价值。但是，由上已知回路阻抗法尚存在下述缺点，即编号方案比较麻烦，网络拓扑描述比较复杂，且由于它只对负荷节点进行编号，无法计算确定中间节点的状态(电压幅值和相角)，计算速度也有待提高等，因此有必要对它进行有效的改进，以促进它的应用。1.2.4 前推回代法基于前推回代法思想的算法很多。一般给定配电网络的始端电压和末端负荷，以馈线为计算基本单位。开始时由末端向始端推算，设全网电压都为额定电压，根据负荷功率由末端向始端逐段推导，仅计算各元件中的功率损耗而不计算电压，求得各条支路上的电流和功率损耗，并据此获得始端功率，这是回代过程；再根据给顶的始

22、端电压和求得的始端功率向末端逐段算电压降落，求得各节点电压，这是前推过程；如此重复上述过程，直至各个节点的功率偏差满足收敛条件为止。这种算法对于纯辐射型网络或单环网络编程简单，求解速度快，但处理网孔能力差，随着网孔数量的增加，算法的收敛性变差，甚至发散。1.3 本文工作本课题目的在于通过对配电网的潮流计算，在满足各种给定的条件下以确定配电网在整个电力系统中稳态运行的情况，要求设计配电网潮流计算算法，并用matlab编程实现。二配电网网络模型2.1元件模型2.1.1 电力线路的数学模型文献4指出：电力线路按照结构可分为架空线路和电力电缆两大类，但是它们可以等效为相同的等值电路。在本文中，用单相等

23、值电路代替三相，一方面由于本文中讨论的是三相对称运行方式，另一方面也假设架空线路都已经整循环换位。以单相等值电路代表三相，虽已简化了不少的运算，但由于电力线路的长度长短不一，例如将每公里的电阻、电抗、电纳、电导都一一绘于图上，所得的等值电路仍十分复杂。何况，严格说来，电力线路的参数是均匀分布的，即使是极短的一段线路，都有相应大小的电阻、电抗、电纳、电导。换言之，即使是如此复杂的等值电路，也不能认为精确。但好在电力线路一般不长，需分析的又往往只是它们的端点状况两端电压、电流、功率，通常不考虑线路的这种分布参数特性，只是在个别情况下才要用双曲函数研究具有均匀分布参数的线路。以下，讨论一般线路的等值

24、电路：所谓一般线路，指中等及中等以下长度线路。对架空线路，这长度大约为 300km对电缆线路，大约为 100km。线路不超过这些数值时，可不考虑它们的分布参数特性，而只用将线路参数简单的集中起来的电路来表示。在以下的讨论中，以 R()、X()、G()、B()分别表示全线路每相总电阻、电抗、电导、电纳。显然，线路长度为 L (km)时R =×L ；X =×L G =×L ；B =×L （21）通常，由于线路导线截面积的选择，如前所述，以晴朗天气不发生电晕为前提，而沿绝缘子的泄漏又很少，可设 G =0。一般线路中，又有短线路和中等长度线路之分。所谓短线路，指

25、长度不超过 100km 的架空线路。线路电压不高时，这种线路电纳 B 的影响不大，可略去。从而，这种线路的等值电路最简单，只有一种串联的总阻抗 Z = R+jX，如图 2-1 所示。 图 2-1 短线路的等值电路显然，如电缆线路不长，电纳的影响不大时，也可以采用这种等值电路。所谓中等长度线路，是指长度在 100-300km 之间的架空线路和不超过 100km 的电力电缆线路。这种线路的电纳B一般不能略去。这种线路的等值电路有型等值电路和 T 型等值电路，如图 2-2、图 2-3 所示。其中，常用的是型等值电路 图 2-2 型等值电路 图 2-3 T 型等值电路在型等值电路中，除串联的线路总阻抗

26、 Z = R+jX外，还将线路的总导纳Y = jB分成两半，分别并联在线路的是末端。在 T 型等值电路中，线路的总导纳集中在中间，而线路的总阻抗则分成两半，分别串联在它的两侧。因此，这两种电路都是近似的等值电路，而且，相互间并不等值，即它们不能用-Y 变换公式相互变换。2.1.2 变压器的等值电路文献4指出：当配电网中存在配电变压器时，通常采用型等值电路和 T 型等值电路两种等值电路，分别如图 2-4、图 2-5 所示（这里只画出双绕组的等值电路），其中，图 2-4 中各参数的计算公式如下：; ; （22）式中：变压器高低压绕组的总电阻()；变压器高低压绕组的总电抗()；变压器的电导(S)；变

27、压器的电纳(S)；变压器的短路损耗(kW)；变压器的额定容量(MVA)；变压器的额定电压(kV)；变压器的短路电压百分值；变压器的空载电流百分值；图 2-4 双绕组变压器的 T 型等值电路型等值电路也就是等值变压器模型：图 2-5 双绕组变压器型等值电路不论采用有名制或标幺制，凡涉及多电压级网络的计算，都必须将网络中所有参数和变量归算至同一电压级。这是因为形或形等值电路做变压器模型时，这些等值电路模型并不能体现变压器实际具有的电压变换功能。但是等值变压器模型则具有这种电压变换功能，它也是运用计算机进行电力系统分析时采用的变压器模型，虽然运用这种模型时并不排斥手算。既然这种模型体现电压变换，在多

28、电压等级网络计算中采用这种变压器模型后，就可以不必进行参数和变量的归算，这正是这种变压器模型的主要特点之一。以下，即介绍这种变压器模型。首先，从一个未作电压归算的简单网络入手。设图 2-6、图 2-7 中变压器的导纳或励磁支路和线路的导纳支路都可略去；设变压器两侧线路的阻抗都未经归算，即分别为高低压侧或、侧线路的实际阻抗，变压器本身的阻抗归在低压侧；设变压器的变比为 k，其值为高、低压绕组电压之比。图 2-6 变压器模型(1)图 2-7 变压器模型(2)显然，在这些假设条件下，如在变压器阻抗左侧串联一变比为K 的理想变压器如图 2-8：图 2-8 变压器模型(3)其效果就如同将变压器及其低压侧

29、线路的阻抗都归算至高压侧，或将高压侧线路的阻抗归算至低压侧，从而实际上获得将所有参数和变量都归算到同一侧的等值网络，只要变压器的变比取的是实际变比，这一等值网络无疑是严格的。因此很容易知道图 2-5 中的参数： （23）附带指出，可以证明，变压器不仅有改变电压大小而且有移相功能时，其变比k 将为复数，这时，仍将得到上面所示的、，但其中、不相等，无源电路的互易特性不复存在，不能用形等值电路表示这种变压器模型，虽然这样不影响运用这种模型进行计算。2.2网络模型有名制：所有参数和变量都以有名单位，如、S、kV(V)、kA(A)、MVA(VA)等表示。标幺制：所有参数和变量都以与他们同名基准值相对的标

30、幺值表示，因此都没有单位。对多电压级网络，变压器模型：采用等值变压器模型时，所有参数和变量可不进行归算；采用有名制或标幺制取决于习惯。在我国，电力工程界使用标幺值已有多年；但在国外，有名制的使用也很普遍。至于变压器模型的使用范围，则泾渭分明。手算时，都是用形或T型等值电路模型；计算机计算时，都是用等值变压器或型等值电路模型。 此外，在制定电力网络等值电路模型时，有时还同时作某些简化，常见的有：线路的电导通常都被略去；变压器的电导有时以具有定值的有功功率损耗的形式出现在电路中；100km以下架空线路的电纳被略去；100300km架空线路或变压器的电纳有时以具有定值的容性或感性无功功率损耗的形式出

31、现在电路中。有时，整个元件，甚至部分系统都可能不包括在等值电路中。例如，将某些发电厂的高压母线看作为可维持给定电压、输出给定功率的等值电源时，这些发电厂内部的元件就不再包括在等值电路中。三：基于matlab的配电网潮流计算算法3.1配电网潮流计算算法原理牛顿拉夫逊迭代法是常用的解非线性方程组的方法，也是当前广泛采用的计算潮流的方法，其标准模式如下。设有非线性方程组 式（3.1）其近似解为。设近似解与精确解分别相差，则如下的关式应该成立： 式（3.2）上式中任何一式都可按泰勒级数展开，由此可得： 式（3.3） 以第一式为例， ，式子中：，分别表示以带入这些偏导数表示式时的计算所得，则是一包含，的

32、高次方与的高阶偏导数相乘的函数。如果与精确解相差不大，则的高次方可以略去，从而也可以略去。由此可得： 式（3.4）或简写为： 式（3.5）式中：称函数的雅克比矩阵，为由组成的列向量，则称不平衡向量的列向量。将带入，可得中的各个元素。然后运用任何一组解线性代数方程的方法，可求得，从而球的经第一次迭代后的新值。再将求得的代入，又可以求得的新值，从而解得以及。如此循环而已，最后可获得足够精确的解。运用这种方法计算时，的初值要选择比较接近他们的精确解，否则迭代过程可能不收敛。将这种情况简单说明如下。设函数的图像如图所示，运用这种方法解算时的修正方程式为按着修正方程式迭代求解过程就如图3-1中由求，的过

33、程。由图可见，如的初值选择的接近其精确解，迭代过程将循序收敛；反之，将不收敛。正因为这样，某些运用牛顿拉夫逊计算潮流的程序中，第一、第二次迭代采用高斯赛德尔法，这是因为后者对的初值的选择没有严格要求。f(x)f(x(0)x(1)x(0)x(2)xx(0)x(1)XY图3-1牛顿拉夫逊发的收敛过程与运用高斯赛德尔法时不同，运用牛顿法拉夫逊法时，可以直接用以求解功率方程。 式（3.6）而为此需将代入 式（3.7）并将实数部分和虚数部分分列 式（3.8a） 式（3.8b）此外，由于系统中还有电压大小给定的PV节点，还应补充一组方程式 式（3.8c）和分别表示迭代过程中求得的节点电压实部与虚部，为PQ

34、节点和PV节点的注入有功功率，为PQ节点的注入有功功率，为PV节点的电压大小。对照式(3.8)、式(3.1)可见，式(3-8)的右端项、分别是给定的注入功率和节点电压大小的平方值，他们就对应于式(3.1)右端项；式(3.2)的左端函数分别是由迭代过程求得的节点电压确定的注入功率和节点电压大小的平方值，它就对应于式(3.1)中的左端函数；于是，式(3.8)中的和就对应于式(3.1)中的。至于修正方程式(3.4)中雅可比矩阵的各个元素，显然就是迭代过程中求得的注入功率各个节点电压大小的平方值相对应的的偏导数。牛顿法的核心便是反复形成并求解修正方程。牛顿法当初始估计值和方程的精确解足够接近时，收敛速

35、度非常快。3.2 matlab的概述目前电子计算机已广泛应用于电力系统的分析计算，潮流计算是其基本应用软件之一。现有很多潮流计算方法。对潮流计算方法有五方面的要求：（1）计算速度快（2）内存需要少（3）计算结果有良好的可靠性和可信性（4）适应性好，亦即能处理变压器变比调整、系统元件的不同描述和与其它程序配合的能力强（5）简单。 MATLAB是一种交互式、面向对象的程序设计语言，广泛应用于工业界与学术界，主要用于矩阵运算，同时在数值分析、自动控制模拟、数字信号处理、动态分析、绘图等方面也具有强大的功能。MATLAB程序设计语言结构完整，且具有优良的移植性，它的基本数据元素是不需要定义的数组。它可

36、以高效率地解决工业计算问题，特别是关于矩阵和矢量的计算。MATLAB与C语言和FORTRAN语言相比更容易被掌握。通过M语言，可以用类似数学公式的方式来编写算法，大大降低了程序所需的难度并节省了时间，从而可把主要的精力集中在算法的构思而不是编程上。另外，MATLAB提供了一种特殊的工具：工具箱（TOOLBOXES）.这些工具箱主要包括：信号处理（SIGNAL PROCESSING）、控制系统（CONTROL SYSTEMS）、神经网络（NEURAL NETWORKS）、模糊逻辑(FUZZY LOGIC)、小波(WAVELETS)和模拟（SIMULATION）等等。不同领域、不同层次的用户通过相

37、应工具的学习和应用，可以方便地进行计算、分析及设计工作。MATLAB设计中，原始数据的填写格式是很关键的一个环节，它与程序使用的方便性和灵活性有着直接的关系。原始数据输入格式的设计，主要应从使用的角度出发，原则是简单明了，便于修改。矩阵是MATLAB数据存储的基本单元，而矩阵的运算是MATLAB语言的核心，在MATLAB语言系统中几乎一切运算均是以对矩阵的操作为基础的。矩阵的基本数学运算包括矩阵的四则运算、与常数的运算、逆运算、行列式运算、秩运算、特征值运算等基本函数运算，这里进行简单介绍。四则运算矩阵的加、减、乘运算符分别为“+，*” ，用法与数字运算几乎相同，但计算时要满足其数学要求 在M

38、ATLAB中矩阵的除法有两种形式：左除“”和右除“/”。在传统的MATLAB算法中，右除是先计算矩阵的逆再相乘，而左除则不需要计算逆矩阵直接进行除运算。通常右除要快一点，但左除可避免被除矩阵的奇异性所带来的麻烦。在MATLAB6中两者的区别不太大。与常数的运算 常数与矩阵的运算即是同该矩阵的每一元素进行运算。但需注意进行数除时，常数通常只能做除数。基本函数运算矩阵的函数运算是矩阵运算中最实用的部分，常用的主要有以下几个：det(a) 求矩阵a的行列式eig(a) 求矩阵a的特征值inv(a)或a (-1) 求矩阵a的逆矩阵rank(a) 求矩阵a的秩trace(a) 求矩阵a的迹（对角线元素之

39、和）我们在进行工程计算时常常遇到矩阵对应元素之间的运算。这种运算不同于前面讲的数学运算，为有所区别，我们称之为数组运算。基本数学运算数组的加、减与矩阵的加、减运算完全相同。而乘除法运算有相当大的区别，数组的乘除法是指两同维数组对应元素之间的乘除法，它们的运算符为“.*”和“./”或“.”。前面讲过常数与矩阵的除法运算中常数只能做除数。在数组运算中有了“对应关系”的规定，数组与常数之间的除法运算没有任何限制。另外，矩阵的数组运算中还有幂运算（运算符为 . ）、指数运算（exp）、对数运算（log）、和开方运算（sqrt）等。有了“对应元素”的规定，数组的运算实质上就是针对数组内部的每个元素进行的

40、。矩阵的幂运算与数组的幂运算有很大的区别。逻辑关系运算 逻辑运算是MATLAB中数组运算所特有的一种运算形式，也是几乎所有的高级语言普遍适用的一种运算。3.3程序设计3.3.1 牛顿-拉夫逊法潮流求解过程以下讨论的是用直角坐标形式的牛顿拉夫逊法潮流的求解过程。当采用直角坐标时，潮流问题的待求量为各节点电压的实部和虚部两个分量由于平衡节点的电压向量是给定的，因此待求两共2(n-1)需要2(n-1)个方程式。事实上，除了平衡节点的功率方程式在迭代过程中没有约束作用以外，其余每个节点都可以列出两个方程式。对PQ节点来说，是给定的，因而可以写出 （3-7）对PV节点来说，给定量是，因此可以列出 （3-

41、8）求解过程大致可以分为以下步骤：（1）形成节点导纳矩阵；（2）将各节点电压设初值U（3）将节点初值代入相关求式，求出修正方程式的常数项向量；（4）将节点电压初值代入求式，求出雅可比矩阵元素；（5）求解修正方程，求修正向量；（6）求取节点电压的新值；（7）检查是否收敛，如不收敛，则以各节点电压的新值作为初值自第3步重新开始进行狭义次迭代，否则转入下一步；（8）计算支路功率分布，PV节点无功功率和平衡节点注入功率。以直角坐标系形式表示：迭代推算式采用直角坐标时,节点电压相量及复数导纳可表示为: (3-9)将以上二关系式代入上式中,展开并分开实部和虚部;假定系统中的第1,2,m号为PQ节点,第m+

42、1,m+2,n-1为PV节点,根据节点性质的不同,得到如下迭代推算式: 对于PQ节点 (3-10)对于PV节点 (3-11)对于平衡节点平衡节点只设一个,电压为已知,不参见迭代,其电压为: (3-12)修正方程两组迭代式中包括2(n-1)个方程.选定电压初值及变量修正量符号之后代入,并将其按泰勒级数展开,略去二次方程及以后各项,得到修正方程如下: (3-13)其中，； 。 (3-14)雅可比矩阵各元素的算式式(3-14)中, 雅可比矩阵中的各元素可通过对式(3-10)和(3-11)进行偏导而求得.当时, 雅可比矩阵中非对角元素为 (3-15)当时,雅可比矩阵中对角元素为: (3-16)由式(3

43、-15)和(3-16)看出,雅可比矩阵的特点：矩阵中各元素是节点电压的函数,在迭代过程中,这些元素随着节点电压的变化而变化；导纳矩阵中的某些非对角元素为零时,雅可比矩阵中对应的元素也是为零.若,则必有；雅可比矩阵不是对称矩阵;雅可比矩阵各元素的表示如下： （3-17） （3-18） （3-19） （3-19） （3-20）3.3.2牛顿拉夫逊法的程序框图输入原始数据启动形成导纳矩阵 给定电压初值、置对于PQ节点，按式3-10计算、对于PU节点，按式3-11计算、是否用式3-15,3-16计算雅克比矩阵各元素解修正方程式，求用，修正节点电压以按系统的潮流分布计算平衡节点功率及线路功率输出以四：算

44、例原始资料及主要参数：22334574567891011121314151617181716151413121110098181919202021212222232324242525262627272828292930303131323233原始资料及主要参数：辐射配电网的线路参数和负载参数支路起点终点电阻(欧姆)电抗(欧姆)有功(KW)无功(KVA)1 120.00580.002910013182 230.03080.0157907503340.0228 0.01151204794450.0238 0.0121606605560.0511 0.0441606306670.0117 0.038

45、620010497780.1068 0.077120014258890.0643 0.0462160110699100.0651 0.0462606021010110.0123 0.0041452981111120.0234 0.0077608951212130.0916 0.0721609611313140.0338 0.04451205991414150.0369 0.0328606211515160.0466 0.0340609471616170.0804 0. 1074609361717180.0457 0.035890371182190.0102 0.0998909701919200

46、.0939 0.08469010992020210.0255 0.0298906582121220.0442 0.052590431223230.0282 0.0192907762323240.0560 0.044242010532424250.0559 0.0437420673256260.0127 0.0065604952626270.0177 0.0090607722727280.0661 0.0583609632828290.0502 0.04371206842929300.0317 0.01611007323030310.0608 0.06011509093131320.0194 0

47、.02262104083232330.0213 0.0331100261参考文献1王淳.一种实用的辐射网潮流算法. 继电器. 2008年12期  2王明岗.基于前推后代法解决环网的配电网潮流计算. 南京工程学院学报(自然科学版). 2009年01期3余加喜，郭志忠.一种面向对象的辐射状配电网潮流计算方法.哈尔滨理工大学学报. 2007年06期4谢开贵，周平. 配电网络潮流计算的递推算法. 电力系统自动化. 2007年04期5张荣，王秀和.改进的带二阶项配电网快速潮流算法. 电工技术学报. 2007年07期6刘耀年，岂小梅.基于回路阻抗法的配电网潮流计算.继电器.2007年08期7卢本

48、初. 陈红坤.辐射配电网潮流计算实用方法. 福建电力与电工.2007年04期  8马瑾. 徐建源.辐射状配电网的潮流计算. 沈阳工业大学学报.2008年01期9昝加强. 骆平.配电网的叠加法潮流计算. 现代电力.2009年06期10颜伟. 刘方.辐射型网络潮流的分层前推回代算法. 中国电机工程学报.2009年08期11汪卫华. 杨期余.配电网潮流计算中简单编号方法. 电力自动化设备.2010年01期 12曹亮. 孔峰.一种配电网的实用潮流算法. 继电器.2012年05期13戴雯霞. 吴捷.基于支路电流的配网潮流前推后代法. 继电器.2012年05期14许梦黔 . 赵武智.基于图形化界

49、面配电网潮流计算的开发. 电力标准化与计量.2010年03期15苏宏田. 王秀玲.二叉树理论在配电网潮流计算中的应用. 内蒙古电力技术.2010年02期16 Ulas Eminoglu and M. Hakan Hocaoglu A new power flow method for radial distribution systems including voltage dependent load models Electric power system 2005.5 、致 谢衷心感谢导师刘会家老师在学术上的悉心指导、深深教诲，在生活上的热情关怀，人生路上给予我的指引。本人大学

50、学习期间在刘老师教导下求学，做人，刘老师渊博的学识、敏锐的思维、民主而严谨的作风使学生受益匪浅、永生难忘；追求真理、献身科学、学而不厌、诲人不倦的崇高品质对学生将是永远的鞭策，他不但教会了我许多理论上的知识、实践中的经验，同时我还从他那里学到了人生的哲理、指引我今后的人生方向。感谢我的父母多年来对我的支持和鼓励，比起他们的付出，我所做的一切是那么的微不足道。最后向鼓励和帮助的各位老师、同学致以最诚挚的谢意和最衷心的祝福。附 录：源程序n=input('请输入节点数：n=');n1=input('请输入支路数：n1=');isb=input('请输入平衡母线节点号：isb=');pr=input('请输入误差精度：pr=');B1=input('请输入由支路参数形成的矩阵：B1=');B2=input(&#