牛顿拉夫逊法潮流计算

牛顿拉夫逊法潮流计算

摘要

此外，在进行电力系统静态及暂态稳定计算时，要利用潮流计算的结果作为其计算的根底；一些故障分析以及优化计算也需要有相应的潮流计算作配合；潮流计算往往成为上述计算程序的一个重要组成局部。以上这些，主要是在系统规划设计及运行方式安排中的应用，属于离线计算范畴。

牛顿－拉夫逊法在电力系统潮流计算的常用算法之一，它收敛性好，迭代次数少。

关键词：电力系统潮流计算，牛顿－拉夫逊法，MATLAB

I

ABSTRACT

ThisarticlefirstintroducestheflowcalculationbasedontheprincipleofMALABBankofChina,meaning,andthenusespecificexamples,abriefintroduction,howtouseMALABtotheflowcalculationinpowersystems.

Asweallknow,isthestudyofpowerflowcalculationofpowersystemsteady-stateoperationofacalculation,whichaccordingtothegivenoperatingconditionsandsystemwiringtheentirepowersystemtodeterminetheoperationalstatusofeachpart:thebusvoltageflowingthroughthecomponentspower,systempowerlossandsoon.Inpowersystemplanningpowersystemdesignandoperationmodeofthecurrentstudy,arerequiredtoquantitativelycalculatedusingthetrendanalysisandcomparisonoftheprogramorrunmodepowersupplyreasonable,reliabilityandeconomy.

Inaddition,duringthepowersystemstaticandtransientstabilitycalculation,theresultsofcalculationtotakeadvantageofthetrendasitsbasisofcalculation;numberoffaultanalysisandoptimizationalsorequiresacorrespondingflowcalculationforcooperation;powerflowcalculationprogramoftenbecometheanimportantpart.These,mainlyinthewayofsystemdesignandoperationarrangementsintheapplicationareasareoff-linecalculation.

Newton-Raphsonpowerflowcalculationinpowersystemisonecommonlyusedmethod,itisgoodconvergenceoftheiterationnumberofsmall,introducethetrendofcomputer-aidedpowersystemanalysisofthebasicknowledgeandpowerflowNewton-Raphsonmethod,introducedbythelastmatlabrunresults.

Keywords：powersystemflowcalculation,Newton–Raphsonmethod,matlab

II

目录

1绪论.............................................................1

1.1课题背景......................................................1

1.2电力系统潮流计算的意义........................................2

1.3电力系统潮流计算的开展........................................2

1.4潮流计算的开展趋势............................................4

2潮流计算的数学模型................................................5

2.1电力线路的数学模型及其应用....................................5

2.2等值双绕组变压器模型及其应用..................................6

2.3电力网络的数学模型............................................8

2.4节点导纳矩阵..................................................9

2.4.1节点导纳矩阵的形成.........................................9

2.4.2节点导纳矩阵的修改........................................10

2.5潮流计算节点的类型...........................................11

2.6节点功率方程.................................................122·7潮流计算的约束条件...........................................13

3牛顿－拉夫逊法潮流计算根本原理...................................14

3.1牛顿-拉夫逊法的根本原理......................................14

3.2牛顿-拉夫逊法潮流计算的修正方程..............................17

3.3潮流计算的根本特点...........................................20

3.4节点功率方程.................................................21

4牛顿－拉夫逊法分解潮流程序........................................224·1牛顿－拉夫逊法分解潮流程序原理总框图.........................22

4.2形成节点导纳矩阵程序框图及代码...............................23

4.2。1形成节点导纳矩阵程序框图..................................23

4.2.2形成节点导纳矩阵的程序代码................................244·3雅克比矩阵求取的程序框图及代码...............................25

4·3·1形成节点导纳矩阵程序框图..................................25

4·3·2形成雅克比矩阵程序的代码..................................25

III

4·4求取Df、De的程序框图及代码..................................28

4·4·1求取Df、De的程序框图.....................................28

4·4·2求取Df、De的程序代码.....................................28

5实例与分析.......................................................29

5.1一个6节点算例..............................................295·2根据算例输入相应节点的线路参数..............................31

5.3算例运行....................................................33

5.3.1原始数据的输入............................................33

5.3.2原始数据输入程序段........................................34

5.3.3导纳矩阵的形成............................................35

5.3.4导纳矩阵Y为..............................................36

5.3.5雅克比矩阵J〔k=0〕.......................................36

5.3.6算例运行结果输出程序段....................................36

5.3.7算例运行结果输出..........................................38

5.4潮流计算GUI界面.............................................40

全文总结...........................................................45

参考文献...........................................................46

致谢...............................................................47

附录A潮流计算程序.................................................48

附录B“去除数据〞按钮回调函数......................................57

附录C“关闭〞按钮回调函数..........................................58

IV

1绪论

1.1课题背景

潮流计算是研究电力系统的一种最根本和最重要的计算，最初，电力系统潮流计算是通过人工手算的，后来为了适应电力系统日益开展的需要，采用了交流计算台。随着电子数字计算机的出现，1956年Ward等人编制了实际可行的计算机潮流计算程序。这样，就为日趋复杂的大规模电力系统提供了极其有利的计算手段。经过几十年的时间，电力系统潮流计算已经开展的十分成熟。潮流计算是研究电力系统稳态运行情况的一种计算，是根据给定的运行的条件及系统接线情况确定整个电力系统各个局部运行的状态，如各母线的电压、各元件中流过的电流、系统的功率损耗等等。电力系统潮流计算是计算系统动态稳定和静态稳定的根底。在电力系统规划设计和现有电力系统运行的方式研究中，都需要利用电力系统潮流计算来定量的比拟供电方案或运行方式的合理性、可靠性和经济性。

电力系统潮流计算分为离线计算和在线计算，离线计算主要用于系统规划设计、安排系统的运行方式，在线计算那么用于运算中系统的实时监测和实施控制。两种计算的原理在本质上是相同的。

实际电力系统的潮流技术主要采用牛顿－拉夫逊法。牛顿－拉夫逊法早在50年代末就已应用于求解电力系统潮流问题，但作为一种适用的、有竞争力的电力系统潮流计算方法，那么是在应用了稀疏矩阵技巧和高斯消元法求修正方程式以后。牛顿－拉夫逊法是求解非线性代数方程有效的迭代计算，本设计就是采用牛顿－拉夫逊法计算电力系统潮流的。

1

1.2电力系统潮流计算的意义

(1)在电网规划阶段，通过潮流计算，合理规划电源容量及接入点，合理规划网架，选择无功补偿方案，满足规划水平的大、小方式下潮流交换控制、调峰、调相、调压的要求。

(2)在编制年运行方式时，在预计负荷增长及新设备投运根底上，选择典型方式进行潮流计算，发现电网中薄弱环节，供调度员日常调度控制参考，并对规划、基建部门提出改良网架结构,加快基建进度的建议。

(3)正常检修及特殊运行方式下的潮流计算，用于日运行方式的编制，指导发电厂开机方式，有功、无功调整方案及负荷调整方案，满足线路、变压器热稳定要求及电压质量要求。

(4)预想事故、设备退出运行对静态平安的影响分析及作出预想的运行方式调整方案。

总之在电力系统运行方式和规划方案的研究中，都需要进行潮流计算以比拟运行方式或规划供电方案的可行性、可靠性和经济性。同时，为了实时监控电力系统的运行状态，也需要进行大量而快速的潮流计算。因此，潮流计算是电力系统中应用最广泛、最根本和最重要的一种电气运算。在系统规划设计和安排系统的运行方式时，采用离线潮流计算；在电力系统运行状态的实时监控中，那么采用在线潮流计算。

1.3电力系统潮流计算的开展

利用电子计算机进行潮流计算从20世纪50年代中期就已经开始。此后，潮流计算曾采用了各种不同的方法，这些方法的开展主要是围绕着对潮流计算的一些根本要求进行的。对潮流计算的要求可以归纳为下面几点：

〔1〕算法的可靠性或收敛性

〔2〕计算速度和内存占用量

〔3〕计算的方便性和灵活性

电力系统潮流计算属于稳态分析范畴，不涉及系统元件的动态特性和过渡过程。因此其数学模型不包含微分方程，是一组高阶非线性方程。非线性代数方程组的解法离不开迭代，因此，潮流计算方法首先要求它是能可靠的收敛，并给出正确答案。随着电力系统规模的不断扩大，潮流问题的方程式阶数越来越高，目前已到达几千阶甚至上万阶，对这样规模的方程式并不是采用任何数学方法都能保证给出正确答案的，这种情况促使电力系统的研究人员不断寻求新的更可靠的计算方法。

在用数字计算机求解电力系统潮流问题的开始阶段，人们普遍采用以节点导

2

纳矩阵为根底的高斯-赛德尔迭代法〔以下简称导纳法〕。这个方法的原理比拟简单，要求的数字计算机的内存量也比拟小，适应当时的电子数字计算机制作水平和电力系统理论水平，于是电力系统计算人员转向以阻抗矩阵为主的逐次代入法〔简称阻抗法〕。

20世纪60年代初，数字计算机已经开展到第二代，计算机的内存和计算速度发生了很大的飞跃，从而为阻抗法的采用创造了条件。阻抗矩阵是满矩阵，阻抗法要求计算机储存表征系统接线和参数的阻抗矩阵。这就需要较大的内存量，而且阻抗法每迭代一次都要求顺次取阻抗矩阵中的每一个元素进行计算，因此，每次迭代的计算量很大。

阻抗法改善了电力系统潮流计算问题的收敛性，解决了导纳法无法解决的一些系统的潮流计算，在当时获得了广泛的应用，曾为我国电力系统设计、运行和研究作出了很大的奉献。但是，阻抗法的主要缺点就是占用计算机的内存很大，每次迭代的计算量很大。当系统不断扩大时，这些缺点就更加突出。为了克服阻抗法在内存和速度方面的缺点，后来开展了以阻抗矩阵为根底的分块阻抗法。这个方法把一个大系统分割为几个小的地区系统，在计算机内只需存储各个地区系统的阻抗矩阵及它们之间的联络线的阻抗，这样不仅大幅度的节省了内存容量，同时也提高了节省速度。

克服阻抗法缺点的另一途径是采用牛顿-拉夫逊法〔以下简称牛顿法〕。牛顿法是数学中求解非线性方程式的典型方法，有较好的收敛性。解决电力系统潮流计算问题是以导纳矩阵为根底的，因此，只要在迭代过程中尽可能保持方程式系数矩阵的稀疏性，就可以大大提高牛顿潮流程序的计算效率。自从20世纪60年代中期采用了最正确顺序消去法以后，牛顿法在收敛性、内存要求、计算速度方面都超过了阻抗法，成为直到目前仍被广泛采用的方法。

在牛顿法的根底上，根据电力系统的特点，抓住主要矛盾，对纯数学的牛顿法进行了改造，得到了P-Q分解法。P-Q分解法在计算速度方面有显著的提高，迅速得到了推广。

牛顿法的特点是将非线性方程线性化。20世纪70年代后期，有人提出采用更精确的模型，即将泰勒级数的高阶项也包括进来，希望以此提高算法的性能，这便产生了保存非线性的潮流算法。另外，为了解决病态潮流计算，出现了将潮流计算表示为一个无约束非线性规划问题的模型，即非线性规划潮流算法。

近20多年来，潮流算法的研究仍然非常活泼，但是大多数研究都是围绕改良牛顿法和P-Q分解法进行的。此外，随着人工智能理论的开展，遗传算法、人工神经网络、模糊算法也逐渐被引入潮流计算。但是，到目前为止这些新的模型和算法还不能取代牛顿法和P-Q分解法的地位。由于电力系统规模的不断扩大，对

3

计算速度的要求不断提高，计算机的并行计算技术也将在潮流计算中得到广泛的应用，成为重要的研究领域。

1.4潮流计算的开展趋势

通过几十年的开展，潮流算法日趋成熟。近几年，对潮流算法的研究仍然是如何改善传统的潮流算法，即高斯-塞德尔法、牛顿法和快速解耦法。牛顿法，由于其在求解非线性潮流方程时采用的是逐次线性化的方法，为了进一步提高算法的收敛性和计算速度，人们考虑采用将泰勒级数的高阶项或非线性项也考虑进来，于是产生了二阶潮流算法。后来又提出了根据直角坐标形式的潮流方程是一个二次代数方程的特点，提出了采用直角坐标的保存非线性快速潮流算法。

4

2潮流计算的数学模型

2.1电力线路的数学模型及其应用

在电力系统稳态分析中的电力线路数学模型就是以电阻、电抗、电纳、电导表示的它们的等值电路。

r1

s式〔2.1〕

式中ρ——为导线材料的电阻率〔Ω•mm2/km〕；

s——为导线的额定截面积〔mm2〕。

x10.1445lgDm

0.0157r式〔2.2〕

式中r——为导线计算半径〔mm或cm〕；

Dm——为几何均距〔mm或cm〕，其单位应与r的相同。

b17.58106式〔2.3〕Dmlgr

Pgg12103式〔2.4〕U

式中b1——导线单位长度的电纳〔S/km〕；

g1——导线单位长度的电导〔S/km〕；

；Pg——三相线路泄漏和电晕损耗功率〔kW/km〕

U——线路线电压〔kV〕。

按上式求得单位长度导线的电阻、电抗、电纳、电导后，就可作最原始的电力线路等值电路图，如图2-1所示。这是单相等值电路。之所以可用单相等值电路代表三相，一方面由于

图2-1中等线路等值模型5

以单相等值电路代表三相虽已简化了不少计算，但由于电力线路的长度往往有数十乃至数百公里，如将每公里的电阻、电抗、电纳、电导都一一绘于图上，所得的等值电路仍十分复杂。何况，严格说来，电力线路的参数并不是均匀分布的，即使是极短的一段线段，都有相应大小的电阻、电抗、电纳、电导。换言之，即使是如此复杂的等值电路，也不能认为精确。但好在电力线路一般都不长，需分析的又往往只是它们的端点状况—两端电压、电流、功率，通常可不考虑线路的这种分布参数特性，只是在个别情况下才要用双曲函数研究具有均匀分布参数的线路。以下，先讨论一般线路的等值电路。

中等长度的线路通常指100km-300km之间的架空线路，这种线路的导纳一般不能略去，常用的是∏型等值电路。当线路长度为l(km)时：

Rrl1,Xx1l

Gg1l,Bb1l

2.2等值双绕组变压器模型及其应用

无论采用有名制或标幺制，凡涉及多电压级网络的计算都必须将网络中所有参数和变量归算至同一电压等级。这是因为以Γ型或T型等值电路做变压器模型时，这些等值电路模型并不能表达变压器实际具有的电压变换功能。以下将介绍另一种可等值的表达变压器电压变换功能的模型，它也是运用计算机进行电力系统分析时采用的变压器模型，虽然运用这种模型时并不排斥手算。既然这种模型可表达电压变换，在多电压级网络计算中采用这种变压器模型后，就可不必进行参数和变量的归算，这正是这种变压器模型的主要特点之一。以下就介绍这种变压器模型。如图2-2所示。

2

T

(b)

U1(a)21U1/k(c)

2

图2-2等值双绕组变压器

6

首先，从一个未作电压级归算的简单网络入手。设图中变压器的导纳或励磁支路和线路的导纳支路都可略去；设变压器两侧线路的阻抗都未经归算，即分别为高、低电压侧Ⅰ、Ⅱ侧线路的实际阻抗，变压器本身的阻抗那么归在低压侧；设变压器的变比k，其值为高、低压绕组电压之比。

显然，在这些假设条件下，如在变压器阻抗ZT的左侧串联以变比为K的理想

变压器如图(2-2c)所示，其效果就如同将变压器及其低压侧线路的阻抗都归算至高压侧，或将高压侧线路的阻抗归算至低压侧，从而实际上获得将所有参数和变量都归算在同一侧的等值网络，只要变压器的变比取的是实际变比，这一等值网络无疑是严格的。

由图(2-2c)可见流入理想变压器的功率为S1

率为S2*U1I1，流出理想变压器的功U1I2，流入流出变压器的功率应该相等，可得：

**

U1I1U1I2k式〔2.5〕从而有：

**

I1I2

另外由图2-2c可以直接得到：

U1kU2ZTI2

联立解方程组：式〔2.6〕式〔2.7〕

U1U2ZTI2I1I2k**

式〔2.8〕可得：112TT1UI1TT式〔2.9〕

即：

7

1Y11Zk2T1Y12ZTkY1

21ZTkY122ZTy12y10y12y12y20y12式〔2.10〕

y21y12的成立表达了无源电路的互易特性，然后令ZTYT，就可以作导纳支路表示的变压器模型如图(2-2e)所示以及以阻抗支路表示的变压器模型如图(2-2f)所示。

其中，y10(1k)(ZTk)，y20(k1)(ZTk)。

以下利用图2-2说明各种不同情况下等值变压器模型的应用，即多电压级网络中变压器和线路参数的计算，以及相应的理想变压器变比的取值。

(1)有名值、线路参数都归算到低压侧。据以图2-2的情况，由图可见，此时线路阻抗分别为上图中Z1,Z2，变压器阻抗那么由‘‘2

PkU22Uk%U22RT;XT；21000SN100SN

相应的理想变压器变比那么为k

是变压器实际变比。

(2)有名值、线路参数都归算到高压侧。这种情况下的线路阻抗分别为U12，这里所取得理想变压器的变比就Z1Z1’;Z2Z2’(

从而理想变压器变比为：

U1N2)U2NU2NU1kU1NU2

2.3电力网络的数学模型

有名制：所有参数和变量都以有名单位，如Ω、S、kV(V)、kA(A)、MVA(VA)等表示。

标幺制：所有参数和变量都以与他们同名基准值相对的标幺值表示，因此都没有单位。

对多电压级网络，变压器模型：采用等值变压器模型时，所有参数和变量可

8

不进行归算；采用有名制或标幺制取决于习惯。在我国，电力工程界使用标幺值已有多年；但在国外，有名制的使用也很普遍。至于变压器模型的使用范围，那么泾渭清楚。手算时，都是用Γ形或T型等值电路模型；计算机计算时，都是用等值变压器或Π型等值电路模型。

此外，在制定电力网络等值电路模型时，有时还同时作某些简化，常见的有：①线路的电导通常都被略去；②变压器的电导有时以具有定值的有功功率损耗的形式出现在电路中；③100km以下架空线路的电纳被略去；④100~300km架空线路或变压器的电纳有时以具有定值的容性或感性无功功率损耗的形式出现在电路中。有时，整个元件，甚至局部系统都可能不包括在等值电路中。例如，将某些发电厂的高压母线看作为可维持给定电压、输出给定功率的等值电源时，这些发电厂式〔2.11〕上式中，IB是节点注入电流的列向量，可理解为某个节点的电源电流与负荷电流之和，并规定电源流向网络的注入电流为正。因此，仅有负荷的负荷节点注入电流就具有负值。UB是节点电压的列向量。因通常以大地作参考节点，网络中有接地支路时，节点电压通常就指该节点的对地电压；网络中没有接地支路时，各节点电压可指各该节点与某一个被选定参考节点之间的电压差。YB是一个节点导纳矩阵，它的阶数n等于网络中除参考节点外的节点数。

它可展开为

I1Y11Y12Y13I2Y21Y22Y23

YY32Y3331I3

YYYn3n1n2In

2.4.1节点导纳矩阵的形成UY1n1Y2nU2Y3nU3式〔2.12〕YnnUn

根据定义直接求取节点导纳矩阵时，注意以下几点：

9

(1)节点导纳矩阵是方阵，其阶数就等于网络中除去参考节点外的节点数。参考节点一般取大地，编号为零。

(2)节点导纳矩阵是稀疏矩阵，其各行非零非对角元素就等于与该行相对应节点所连接的不接地支路数。

(3)节点导纳矩阵的对角元素就等于该节点所连接导纳的总和。因此，与没有接地支路的节点对应的行或列中，对角元素为非对角元素之和的负值。

(4)节点导纳矩阵的非对角元素等于连接节点i，j支路导纳的负值。因此，在一般情况下，节点导纳矩阵的对角元素往往大于非对角元素的负值。

(5)节点导纳矩阵一般是对称矩阵，这是网络的互易特性所决定的。从而，一般只要求取这个矩阵的上三角或下三角局部。

(6)网络中的变压器。

2.4.2节点导纳矩阵的修改

(1)从原有网络引出一支路，同时增加一节点。

设i为原有网络中的节点，j为新增加的节点，新增加支路导纳为

新增一节点，节点导纳矩阵将增加一阶。

新增的对角元Yij，由于在节点j上只有一个支路

(2)在原有网络的节点i、j之间增加一支路。

这时由于仅增加支路不增加节点，节点导纳矩阵阶数不变，但与节点i、j有关的元素应作一下修改，其增量为:

yii

切除一导纳为

yij。那么因yij，将为Yij=yij；新增的非对角元YiiYjiyij；原有矩阵中的对角元Yii将增加yii,yiiyij。yiiyij,yijyjiyij(3)在原有网络的节点i，j之间切除一支路。yij的支路，相当于增加一导纳为yij的支路，从而与节点i、j有关的元素应作如下修改：Yiiyij,Yjjyij,YijYijyij‘(4)原有网络的节点i、j之间的导纳由yij改变为yij。

‘yy这种情况相当于切除一导纳为ij的支路，并增加一导纳为ij的新支路。从

而与节点i、j有关的元素应作如下修改：

‘‘Yyy,Yyiiijijjjijyij,YijYijyijyij’(5)原有网络节点i、j之间变压器的变比由k改变为k。

‘‘这种情况相当于在i、j节点之间并联一个变比为k的变压器，再并联一个变比为k的变压器，即相当于修改变压器。修改前，i、j节点之间的自导纳和互导纳为：

10

Yii

yT(k1)y(1k)yTyT

yTyT,YjjTy,YyTijijkkkk2k2k

‘

修改后，引用“理想变压器〞的π型等值电路，变压器变比由k改变为k时，原网中与节点i、j有关的元素应作如下修改：

Yii0,Y

jj

1

’2k12y,YTk

ij

1Yjk1y’k

T

2.5潮流计算节点的类型

用一般的电路理论求解网络方程，目的是给出电压源〔或电流源〕研究网络内的电流(或电压)分布，作为根底的方程式，一般用线性代数方程式表示。然而在电力系统中，给出发电机或负荷连接母线上电压或电流(都是向量)的情况是很少的，一般是给出发电机母线上发电机的有功功率P和母线电压的幅值U，给出负荷母线上负荷消耗的有功功率P和无功功率Q。主要目的是由这些量去求电力系统内的各种电气量。所以，根据电力系统中各节点性质的不同，很自然地把节点分成三类：

(1)PQ节点

对这类节点，等值负荷功率PGi、QLi和等值电源功率PGi、QGi是给定的，从而注入功率Pi、Qi是给定的，待求的那么是节点电压的大小Ui和相位角i。属于这一类节点的有按给定有功无功功率发电的发电厂母线和没有其他电源的变电所母线。

(2)PV节点

对这类节点，等值负荷和等值电源的有功功率PLi、PGi是给定的，从而注入有功功率Pi是给定的。等值负荷的无功功率QLi和节点电压的大小Ui也是给定的。待求的那么是等值电源的无功功率QGi，从而注入无功功率Qi和节点电压的相位角i。有一定无功功率储藏的发电厂和一定无功功率电源的变电所母线都可选作为PV节点。

(3)平衡节点

潮流计算时，一般只设一个平衡节点。对这节点，等值负荷功率是给定的，节点电压的大小US和相位角S也是给定的，如给定US=1.0、S=0。待求的那么是等值电源功率PGs、QGs，从而注入功率Ps、Qs。担负调整系统频率任务的发电厂母线往往被选作为平衡节点。例如，为提高计算的收敛性。可以选择出线数多或者靠近电网中心的发电厂母线作平衡节点。

进行计算时，平衡节点是不可少的；PQ节点是大量的；PV节点较少，甚至

11

可能没有。

2.6节点功率方程

节点电压向量可以表示为极坐标的形式，也可以表示为直角坐标的形式，与此相对应，在潮流计算中节点功率方程也有两种形式。

节点功率可表示为：

n

YV〔i=1,2,…n〕PijQiV式〔2.13〕iijj

ji

如果上式中电压向量表示为极坐标的形式：V

iVjiie导纳矩阵中元素表示为：YijGijjBij因此：P

jijQi=Vie

i

(GijjBjj

ij)Vje

〔i=1,2,„n〕ji

又由e

j

cosjsin那么可以得到：PijQiViVj(GijjBij)(cosijjsinij)ji

式中：ijij为两个节点电压的相位差。

将上式按实部和虚部展开，得到：

n

PV

iVij(GijcosijBijsinij)ji

n

QiViVj(GijsinijBijcosij)

ji

这就是功率的极坐标方程式。

把上式中个节点的电压向量表示为直角坐标：

Vieijfi

eiVicosi，fiVisini代入式：PijQiV

i

n

Yi

V

j

ji

即可得到：

12

式

〔2.14〕式〔2.15〕

式〔2.16〕

式(2.17)

式〔2.18〕

式〔2.19〕式〔2.20〕

式(2.21)

式(2.22)

Pieiaifibi〔i=1,2,„n〕式(2.23)Qifiaieibi

(Geij

jijBijfi)ai

式(2.24)

式中(GjiijifBijej)bi

这就是功率的直角坐标方程式。

2·7潮流计算的约束条件

通过方程的求解所得到的计算结果代表了潮流方程在数学上的一组解答。但是，这组解答所反映的系统运行状态在工程上是否具有实际意义呢？这还要进行检验。因为电力系统运行时还必须满足一定技术上和经济上的要求。这些要求构成了潮流问题中某些变量的约束条件，常用约束条件有：

〔1〕所有节点电压必需满足：

ViminViVimax(i=1,2,3…n)

从保证电能质量和供电平安的要求来看，电力系统的所有电气设备必需运行在额定电压附近，PV节点的电压幅值必需按上述条件给定。因此，这一约束条件主要是对PV节点而言。

〔2〕所有电源节点的有功功率和无功功率必需满足：

QiminQiQimax

由于PQ节点的有功功率和无功功率以及PV节点的有功功率属于扰动变量不可控，对它们没有约束。对平衡节点的P和Q以及PV节点的Q应按上述条件进行检验。〔3〕某些节点之间电压的相位差应满足：PiminPiPimax

ijijmax

为了保证系统运行的稳定性，要求某些输电线路两端的电压相位差不超过一定的数值。因此，潮流计算可以归结为求解一组非线性方程组，并使其解满足一定的约束条件，如不满足，那么应修改某些变量的给定值，甚至修改系统运行方式，重新计算。

13

3牛顿－拉夫逊法潮流计算根本原理

3.1牛顿-拉夫逊法的根本原理

牛顿—拉夫逊迭代法是常用的解非线性方程组的方法，也是当前广泛采用的计算潮流的方法，其标准模式如下。

设有非线性方程组

f1x1,x2,,xny1

f2x1,x2,,xny2

式〔3.1〕

fnx1,x2,,xnyn

x1,x2,,xn

其近似解为。设近似解与精确解分别相差

x1,x2,,xn

，那么如下的

关式应该成立：

000

f1x1x,1x2x,,nxny12x

000f2x1x,1x2x,x,xy22nn

式〔3.2〕

000

fnx1x1,x2x,x,nxnyn2

上式中任何一式都可按泰勒级数展开，由此可得：

f1x1,x2,,xn

000

f2x1,x2,,xn

fnx1,x2,,xn

0

f1f1f1

x1x2xny1x10x20xn0

f2f2f2x1x2xny2x10x20xn0〔3.3〕式

fnfnfn

x1x2xnynx10x20xn0

f1f

x11x2x1x2

以第一式为例，

(0)(0)(0)(0)(0)(0)

f1(x1x1,x2x2,...,xnxn)f1(x1,x2,...,xn)

f1f1f1f1

，，…，,...,xn1y1，式子中：

x10x20xnxn

(0)(0)

分别表示以x1(0),x2,...,xn

带入这些偏导数表示式时的计算所得，1那么是一包含x1，x2，…，xn的高次方

14

与f1的高阶偏导数相乘的函数。如果xi(0)与精确解相差不大，那么xi的高次方可以略去，从而1也可以略去。

由此可得：

y1f1x1,x2,,xn

yfx0,x0,,x0

2n111

yfx0,x0,,x0

2n

111

000

f1f1f1

xn

x10x20

f1f1f1x10x20xnf1f1f1

xn

x10x20

0x1x2

式〔3.4〕0

xn0

或简写为：

的列向量。

将xi带入，可得f、J中的各个元素。然后运用任何一组解线性代数方程的方法，可求得xi

(0)

fJx式〔3.5〕

式中：J称函数fi的雅克比矩阵，x为由xi组成的列向量，f那么称不平衡向量

，从而球的经第一次迭代后xi的新值xi

(1)

(1)

xi(0)xi(0)。再将求

(2)

得的xi(1)代入，又可以求得f、J的新值，从而解得xi以及xi此循环而已，最后可获得足够精确的解。

xi(1)xi(1)。如

运用这种方法计算时，xi的初值要选择比拟接近他们的精确解，否那么迭代过程可能不收敛。将这种情况简单说明如下。设函数的图像如下图，运用这种方法解算f(x)y时的修正方程式为

yf(x(k))

df

x(k)dxk

按着修正方程式迭代求解过程就如图3-1中由x(0)求x(1)，x(2)……的过程。由图可见，如x的初值x(0)选择的接近其精确解，迭代过程将循序收敛；反之，将不收敛。正因为这样，某些运用牛顿拉夫逊计算潮流的程序中，第一、第二次迭代采用高斯赛德尔法，这是因为后者对xi的初值的选择没有严格要求。

15

图3-1牛顿－拉夫逊发的收敛过程

与运用高斯赛德尔法时不同，运用牛顿法拉夫逊法时，可以直接用以求解功率方程。

而为此需将Yij

jnj1

Ui

Y

ij

UjPijQi式〔3.6〕

GijjBij,Ueijfi代入

eijfiGijjBijejjfiPijQi式〔3.7〕

j1jn

并将实数局部和虚数局部分列

〔3.8a〕eGeiijjBijfjfiGijfjBijejPi式j1

jn

fGe

i

ij

j1

jn

j

〔3.8b〕BijfjeiGijfjBijejQi式

此外，由于系统中还有电压大小给定的PV节点，还应补充一组方程式

ei2fi2Ui2式〔3.8c〕

ei和fi分别表示迭代过程中求得的节点电压实部与虚部，Pi为PQ节点和PV

节点的注入有功功率，Qi为PQ节点的注入有功功率，Ui为PV节点的电压大小。

对照式(3.8)、式(3.1)可见，式(3-8)的右端项Pi、Qi、Ui分别是给定的注入功率和节点电压大小的平方值，他们就对应于式(3.1)右端项yi；式(3.2)的左端函数分别是由迭代过程求得的节点电压确定的注入功率和节点电压大小的平方值，它就对应于式(3.1)中的左端函数fi(x1,x2,...,xn)；于是，式(3.8)中的ei和fi...就对应

2

...,xn。至于修正方程式(3.4)中雅可比矩阵的各个元素，显然于式(3.1)中的x1、x2、

...的就是迭代过程中求得的注入功率各个节点电压大小的平方值相对应的ei、fi、

16

偏导数。

牛顿法的核心便是反复形成并求解修正方程。牛顿法当初始估计值和方程的精确解足够接近时，收敛速度非常快。

3.2牛顿-拉夫逊法潮流计算的修正方程

牛顿潮流计算的核心问题是修正方程式的建立和求解。为说明这一修正方程的建立过程，先对网络中个节点的编号作如下约定：

(1)网络中共有n个节点，编号为1,2,3,…,n，其中包含一个平衡节点，编号为s；

(2)网络中有m个PQ节点，编号为1,2,3,…,m，其中包含编号为s的平衡节点；(3)网络中有n-m个PV节点，编号为m+1,m+2,…,n。

据此，由式〔3-8a〕、(3-8b)、(3-8c)所组成的方程式组中共有2(n-1)个独立方程式。其中，式〔3-8a〕类型的有(n-1)个，包括除平衡节点为所有节点有功功率Pi的表示式，即i=1,2,3,…,n，is；式〔3-8b〕类型的有(m-1)个，包括所有PQ节点无功功率Qi的表示式，即i=1,2,3,…,n，is；式〔3-8c〕类型的有(n-1)-(m-1)=n-m个，包括所有PV节点电压Ui2的表达式，i=m+1,m+2,…,n。平衡节点s的功率和电压之所以不包括在这个方程式(3-10)NpnfpSpnepNnnfn

eSnnn

式中的Pi(0)、Qi(0)以及Ui(0)2分别为注入功率的节点电压平方的不平衡量。由式(3-8)可见，他们分别为：

17

〔3.11a〕PiPi[ei(GijejBijfj)fi(GijfjBijej)]式

j1

jn

QiQi[fi(GijejBijfj)ei(GijfjBijej)]

j1

jn

式〔3.11b〕

UU(efi)式(3.11c)

式子中的雅可比矩阵的各个元素分别为：

2

i2i2i

2

PPii

Hij;Nij

fjejQiQi

Jij;Lij

fjejRij

UiUi

;Sij

fjej

2

2

式(3.12)

为取这些偏导数，可将Pi、Qi以及Ui2分别展开如下：

〔3.13a)Piei(GiieiBiifi)fi(GiifiBiiei)[ei(GijejBijfj)fi(GijfjBijej)]式

j1

jijn

Qifi(GiieiBiifi)ei(GiifiBiiei)[fi(GijejBijfj)ei(GijfjBijej)]

jnj1ji

式(3.13b)

Ui2(ei2fi2)

式(3.13c)

当ji时，由于对特定的j，只有该特定节点的fi和ei是变量，由式(3.12)、式(3.13)可得：

PiPi

HijBijeiGijfi;NijGijeiBijfi

fjej

QiQi

BijfiGijeiNij;LijGijfiBijeiHij式(3.14a)Jijfjej

22

UiUi

Rij0;Sij0

fjej

当

ji时，为使这些偏导数的表示式更为简洁，先引入节点诸如电流的表示

式如下：

18

IYiiUiYijUj

j1ji

jn

[(GijeiBijfi)(GijejBijfj)]j[(GiifiBiiei)(GijfjBijej)]

j1ji

j1ji

jnjn

aiijbii

然后由式(3-12)、式(3-13)和上式可得

PiGijfiBijei

Hij

fjGiifiBiieibii

PiBijfiGijei

Nij

ejBijfiGijeiaii

ijij

ijij

ijQiBijfiGijei

Jij

BfGeafjijijiiiiji

ij

QiGijfiBijei

LijijGfBebejijiijiii

ijUi20Rij

ij2ffjiUi20

Sij

ej2ei

ij

式〔3.14b〕ij

由式(3.14a)可见，如果YijGijjBij0，即节点i、j之间无联系，这些元素

Hij

都等于零。从而，将雅可比矩阵分块，把每个2×2阶子阵(

Jij

NijLij

Hij

Rij

Nij

)Sij

作为分块矩阵的元素时，分块雅可比和节点导纳矩阵YB将有相同的结构，所以分块雅可比矩阵和节点导纳矩阵的结构相同是一个可以利用的特点。

仔细分析该修正方程可以发现如下特点。

〔1〕修正方程数目为2(n1)1个，在PV节点所占的比例不大时，修正方程的数目接近2(n1)个。

〔2〕雅可比矩阵的元素都是节点电压的函数，每次迭代，雅可比矩阵都需要重新形成。

〔3〕分析雅可比矩阵的非对角元素的表示式可见，某个非对角元素是否为零决定于相应节点导纳矩阵元素Yji是否为零。因此如将修正方程式按节点号的次序

19

Hij

排列，并将雅可比矩阵分块，把每个2×2阶子阵(如

MijNijLijHij

RijNij

)作Sij

为分块矩阵的元素，那么按节点号顺序而构成的分块雅可比矩阵将和节点导纳矩阵

具有相同的稀疏结构，是一个高度稀疏的矩阵。

〔4〕由于HijHji,NijNji，所以雅可比矩阵不是对称矩阵。

3.3潮流计算的根本特点

形成了雅可比矩阵并建立了修正方程式，运用牛顿-拉夫逊计算潮流的核心问题已经解决，已有可能列出根本计算步骤并编制流程图。

显然，虽修正方程有两种不同表示方法，但是牛顿拉夫逊潮流计算的根本步骤却总不外乎如下几步：

〔1〕形成节点导纳矩阵YB。〔2〕设个节点电压的初值ei(0)、fi(0)。

〔3〕将各节点电压的初值带入式中求修正方程式中的不平衡量Pi(0)、Qi(0)

以及Ui(0)2。

〔4〕将节点电压的初值代入式〔3.14b〕求修正方程式的系数矩阵，即雅可比矩阵的各个元素Hij、Nij、Jij、Lij以及Rij、Sij。

〔5〕解修正方程式，求各节点电压的变化量，即修正量ei(0)、fi(0)。〔6〕计算各节点电压的新值，即修正后值:

ei(1)ei(0)ei(0);fi(1)fi(0)fi(0)

〔7〕运用各节点电压侧初值自第三步开始进入下一次迭代。〔8〕计算平衡节点功率和线路功率。其中平衡节点功率为：

SsUsYsiUiPsjQs

j1

~

jn

式(3.15)

线路功率为：

SijUiIijUi[Uiyi0(UiUj)yij]PijjQijSjiUjIjiUj[Ujyj0(UjUi)yji]PjijQji

~

~

式(3.16)式(3.17)

从而，线路上的损耗功率为

SijSijSjiPijjQij

式(3.18)

20

图3-3线路上流通的电流和功率

3.4节点功率方程

本节主要讨论：复杂网络数学模型的建立，节点功率方程的计算。

在2.4节中我们得到了电力网络方程的系数矩阵即导纳矩阵。建立了节点导纳矩阵YB，就可以进行潮流分布计算。但由于工程实践中通常的是各节点的功率SB，实际计算时，几乎无一例外地要迭代解非线性的节点电压方程YBUBS。故应用联系节点电流和节点功率的关系式：

UB

PijQi

UiYijUjj1

n

这就是潮流问题最根本的方程式，是一个以节点电压Ui为变量的非线性方程组，并且通过迭代来求解

[e(Gei

iijjBijfj)fi(GijfjBijej)]PiBijfj)ei(GijfjBijej)]Qi

形式的潮流方程式。这两种形式的潮流方程统称为节点功率方程，是牛顿拉夫逊法的主要的数学模型。

节点功率方程可以通过牛顿拉夫逊法来有效的解算。我们对于不同的类型的节点，根据以上两式子得到牛顿拉夫逊算法的修正方程。在一点运用泰勒级数展开，略去二阶以上项可以得到雅可比矩阵个元素Hij、Nij、Jij、Lij、Rij、Sij。最后反复迭代形成并且求解修正方程式，从而得到节点功率方程较为精确的解。

[f(Geijj

21

4牛顿－拉夫逊法分解潮流程序

4·1牛顿－拉夫逊法分解潮流程序原理总框图

图4-1牛-拉法分解程序总框

对于图4-1，相关的计算公式如下；

(1)PiPiGijejBijfjeiGijfjBijejfi；

j=1j=n

22

(2)QiQiGijejBijfjfiGijfjBijejei；

j=1

j=n

(3)Ui2Ui2ei2fi2

(4)HijBijeiGijfi；NijGijeiBijfi；

JijNij；Lij=Hij；Rij0；Sij0

(5)HiibiiGiifiBiiei；NiiaiiGiieiBiifi；

JiiaiiGiieiBiifi；LiibiiGiifiBiiei；Rii2fi；Sii2ei

。

4.2形成节点导纳矩阵程序框图及代码

4.2。1形成节点导纳矩阵程序框图

图4-2形成节点导纳矩阵程序框图

对于图4-2相关的计算公式如下:

Y(p,p)ZLk2B2;Y(p,q)ZLk;Y(q,p)Y(p,q);Y(q,q)ZLB2

23

4.2.2形成节点导纳矩阵的程序代码

Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);

O=zeros(1,n);S1=zeros(nl);

fori=1:n

ifA(i,2)~=0;

p=A(i,1);

Y(p,p)=1./A(i,2);

end

end

fori=1:nl

ifB1(i,6)==0

p=B1(i,1);q=B1(i,2);

elsep=B1(i,2);q=B1(i,1);

end

Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5));

Y(q,p)=Y(p,q);

Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5))+B1(i,4)./2;

Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2;

end

%求导纳矩阵

disp(‘导纳矩阵Y=‘);

disp(Y);

G=real(Y);B=imag(Y);

其中，B1=由各支路参数形成的矩阵:;

%它以矩阵形式存贮支路的情况,每行存贮一条支路

%第一列存贮支路的一个端点

%第二列存贮支路的另一个端点

%第三列存贮支路的阻抗

%第四列存贮支路的对地导纳

%第五列存贮变压器的变比

%第六列存贮支路的序号

B2=由各节点参数形成的矩阵:;

%第一列为电源侧的功率

%第二列为负荷侧的功率

24

%第三列为该点的电压初始值

%第四列为该节点的幅值

%第五列为该节点对地导纳值

%第六列为该点的类型:1为平衡节点,2为PQ节点,3为PV节点

A=由节点号及其对地阻抗形成的矩阵：;

%第一列位节点号

%第二列位对地阻抗

4·3雅克比矩阵求取的程序框图及代码

4·3·1形成雅克比矩阵程序框图

4-3形成雅克比矩阵程序框图

4·3·2形成雅克比矩阵程序的代码

P=real(S);Q=imag(S);

k=0;IT=1;a=0;%k----迭代次数；;IT----没有到达精度要求个数

25

whileIT~=0

IT=0;a=a+1;num2=0;num1=0;

fori=1:n

ifi~=isb

C(i)=0;

D(i)=0;

forj=1:n

C(i)=C(i)+G(i,j)*e(j)-B(i,j)*f(j);

D(i)=D(i)+G(i,j)*f(j)+B(i,j)*e(j);

end

P1=C(i)*e(i)+f(i)*D(i);

Q1=f(i)*C(i)-D(i)*e(i);

V1=e(i)+f(i);

end

%针对PQ节点求雅克比矩阵各参数H,N,J,L

ifB2(i,6)~=1&B2(i,6)~=3%判断该节点是PQ节点：B2(i,6)=2num1=num1+1;%1---平衡节点,2---PQ节点,3---PV节点DP=P(i)-P1;

DQ=Q(i)-Q1;%有功功率偏移量开始

num1=num1+1;

DY(num1,1)=DQ;

disp(‘各节点功率不平衡量:DY=‘);

disp(DY);

forj=1:n

ifj==i&j~=ph

X1=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);%X1=H(i,i)X2=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);%X2=N(i,i)X3=-C(i)+G(i,i)*e(i)+B(i,i)*f(i);%X3=J(i,i)X4=D(i)+B(i,i)*e(i)-G(i,i)*f(i);%X4=L(i,i)p=2*i-1;q=2*j-1;

J(p,q)=X1;

m=q+1;

J(p,m)=X2;

26

DY(num1,1)=DP;%各节点功率偏移量形成的矩阵：10*1,先从PQ的

p=p+1;

J(p,q)=X3;

q=q+1;

J(p,q)=X4;

elseifj~=i&j~=ph

X1=B(i,j)*e(i)-G(i,j)*f(i);%X1=H(i,j)X2=-G(i,j)*e(i)-B(i,j)*f(i);%X2=N(i,j)X3=-X2;%X3=J(i,j)=-N(i,j)=-X2X4=X1;%X4=L(i,j)=H(i,j)=X1p=2*i-1;q=2*j-1;

J(p,q)=X1;

m=q+1;

J(p,m)=X2;

p=p+1;

J(p,q)=X3;

q=q+1;

J(p,q)=X4;

end

end

elseifB2(i,6)~=1&B2(i,6)~=2%判断该节点是PV节点：B2(i,6)=3%针对PV节点求雅克比矩阵各参数H,N,R,S

DP=P(i)-P1;

num1=2*i-1;

DY(num1,1)=DP;

DV=V(i)-V1;

num1=num1+1;

DY(num1,1)=DV;

forj=1:n

ifj~=ph&j~=i

X1=G(i,j)*e(i)-B(i,j)*f(i);%X1=H(i,j)X2=-B(i,j)*e(i)-G(i,j)*f(i);%X2=N(i,j)X5=0;%X5=R(i,j)X6=0;%X6=S(i,j)p=2*i-1;q=2*j-1;

27

J(p,q)=X5;

m=p+1;

J(m,q)=X1;

q=q+1;

J(p,q)=X6;J(m,q)=X2;

elseifj==i&j~=ph

X1=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);%X1=H(i,i)X2=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);%X2=N(i,i)X5=-2*e(i);%X5=R(i,i)X6=-2*f(i);%X6=S(i,i)p=2*i-1;q=2*j-1;

J(p,q)=X1;

m=p+1;

J(m,q)=X5;

q=q+1;

J(p,q)=X2;J(m,q)=X6;

end

end

end

end

disp(‘雅克比矩阵J=‘);

disp(J);

4·4求取Df、De的程序框图及代码

4·4·1求取Df、De的程序框图

如图4-4所示。

4·4·2求取Df、De的程序代码

%求解修正方程，即各节点电压修正量DX=[Df,De],DY=J*DX----->DX=inv〔J〕*DYDX=J\DY;

disp(‘各节点电压修正量:DX=‘)

disp(DX)

%求各节点电压新值

fori=1:(n-1)

28

p=2*i-1;

Df=DX(p);

q=2*i;

De=DX(q);

e(i)=e(i)-De;

f(i)=f(i)-Df;

end

图4-4求取Df、De的程序框图

5实例与分析

5.1一个6节点算例

29

(1)如下列图所示的6节点电力系统：

图5-1六节点电力系统图

该系统中，节点1是平衡节点，保持U1=1.05+j0为定值，节点6为PV节点，

其他四个节点都是PQ节点。给定的注入电压标幺值、线路阻抗标幺值、输出功率标幺值分别为表5-1、表5-2、表5-3中的数据。线路对地导纳标幺值之半Y0=j0.25

及线路阻抗标幺值、输出功率标幺值和变压器变比标幺值如图5-2所示的注释。

表5-1各节点电压标幺值参数

表5-2各节点电压标幺值参数

30

表5-3各节点电压标幺值参数

图5-2电力系统等值阻抗图

5·2根据算例输入相应节点的线路参数

1.n=6(节点数)；

2.n1=6(支路数)；

3.B1(由各支路参数形成的矩阵)，如表5-4所示；

4.输入修正值；ip=0.00001；

5.B2(由各节点参数形成的矩阵)，如表5-5所示；

6.A(节点对地导纳矩阵)，如表5-6所示；

31

表5-4由各支路参数形成的矩阵

表5-5由各支路参数形成的矩阵

注：1——平衡节点，2——PQ节点，3——PV节点.

表5-7由节点对地导纳矩阵形成的矩阵

但是这种将平衡节点编为1号节点，不利于程序的实现。为了方便编程，故对各节点重新编号。将原来的1号节点重新编为6号节点，其余五个节点号依次前移一位。即2号重新编为1号，3号重新编为2号，4号重新编为3号，5号重新编为4号，6号重新编为5号。

重新编号后，相应各支路参数、各节点参数都会发生变化。B1(由各支路参数形成的矩阵)将变为表5-8，B2(由各节点参数形成的矩阵)将变为表5-9，A(节点对地导纳矩阵)未发生变化。

表5-8由各支路参数形成的矩阵

32

表5-9由各支路参数形成的矩阵

5.3算例运行

5.3.1原始数据的输入

为了让程序更有通用性，以及方便初学用户使用，近几年，MATLAB编程大多都用到GUI功能。本次设计也不例外，用到了简单的人机对话界面作为原始数据输入界面。

在这次设计中，用到了一个Excel表格作为原始数据输入界面，通过该界面，用户不用依照矩阵的形式，将一连串的数据输入，而是，按照图表的提示，在图表中填入要求的电力系统节点、支路的参数，节点、支路个数以及要求精度即可。数据输入界面如图5-3所示。

在此，对数据的输入有以下几点说明：

①在节点信息里，节点电压为迭代计算时所设的初值。

②在节点信息里，节点类型一栏中，―1‖表示平衡节点，―2‖表示PQ节点，―3‖表示PV节点。

③KT一栏要求输入的是变压器的变比，非标准变比变压器，KT=k〔k1〕，标准变压器KT=1，假设该线路无变压器，KT=0。

④输入变压器电阻、电抗时，如无特殊说明，均采用归算到低压侧的数值，再进行计算。

33

⑤本条程序默认节点数，支路数均在在

100以下，可以解决绝大局部电力系统的潮流问题，假设遇到超大系统，可对程序做稍加调整，仍然适用。

图5-3Excel输入界面

5.3.2原始数据输入程序段

[x]=xlsread(‘sj.xls’,’A3:A3’)

[y]=xlsread(‘sj.xls’,’B3:B3’)

[z]=xlsread(‘sj.xls’,’C3:C3’)

pr=xlsread(‘sj.xls’,’B4:B4’)

[zhilu]=xlsread(‘sj.xls’,’E4:K100’)

[point]=xlsread(‘sj.xls’,’M4:V100’)

[duidi]=xlsread(‘sj.xls’,’X4:Y100’)

n=x;

n1=y;

ph=z;

B1=zhilu;

B2=point;

A=duidi;

T1=zeros(n,2);

T2=zeros(n1,3);

i=sqrt(-1);

34

forj=1:n

T1(j,1)=B1(j,3)+B1(j,4)*i;

T1(j,2)=B1(j,5)*i;

end

forj=1:n1

T2(j,1)=B2(j,2)+B2(j,3)*i;

T2(j,2)=B2(j,4)+B2(j,5)*i;

T2(j,3)=B2(j,6)+B2(j,7)*i;

end

B1=zeros(n,6);

B2=zeros(n1,6);

forj=1:n

B1(j,1)=zhilu(j,1);

B1(j,2)=zhilu(j,2);

B1(j,3)=T1(j,1);

B1(j,4)=T1(j,2);

B1(j,5)=zhilu(j,6);

B1(j,6)=zhilu(j,7);

end

forj=1:n1

B2(j,1)=T2(j,1);

B2(j,2)=T2(j,2);

B2(j,3)=T2(j,3);

B2(j,4)=point(j,8);

B2(j,5)=point(j,9);

B2(j,6)=point(j,10);

end

5.3.3导纳矩阵的形成

①求导纳矩阵Y中的非对角元元素Yij，假设无变压器,那么Yij直接为线路阻抗分之一取负值，假设有变压器，Yij为线路阻抗乘以KT后分之一再取负值。

②求导纳矩阵Y中对角元元素Yii，无变压器时Yii为Yij加上线路对地电导的一半乘j，有变压器时，对角元元素就与所输入的折算到哪一侧有关，如果支路起始端处于高压侧，支路起始节点的自导纳中要加上变压器等值导纳模型的对地支路的

35

(1-KT)/KT倍，支路终止节点的自导纳要加上变压器等值导纳模型的对地支路的(KT-1)/KT倍，如果支路起始端处于低压侧，情况正好相反支路起始节点的自导纳中要加上变压器等值导纳模型的对地支路的(KT-1)/KT倍，支路终止节点的自导纳要加上变压器等值导纳模型的对地支路的(1-KT)/KT倍。

5.3.4导纳矩阵Y为

5.3.5雅克比矩阵J〔k=0〕

5.3.6算例运行结果输出程序段

本程序在结束时，会通过fopen这个函数，在MATLAB运行程序所储存的文件夹内生成一个TXT文

fid=fopen(‘jisuanjieguo.txt’,’wt’);

fprintf(fid,’******************************电力系统08李进才121408131******************************\n’);

36

fprintf(fid,’*****************潮流计算输出结果*************\n’);fprintf(fid,’迭代次数k为:%d\n’,k);

fprintf(fid,’================================================\n’);fprintf(fid,’各节点的实际电压标么值U为(节点号从小到大排列):\n’);form=1:n

fprintf(fid,’第%d个实际电压标么值U:%f+i*(%f)\n’,m,real(U(m)),imag(U(m))’);

end

fprintf(fid,’================================================\n’);fprintf(fid,’各节点的电压大小V为(节点号从小到大排列):\n’);form=1:n

fprintf(fid,’第%d个节点的电压大小V为(节点号从小到大排列):%f\n’,m,V(m));

end

fprintf(fid,’================================================\n’);fprintf(fid,’各节点的电压相角shita为(节点号从小到大排列):\n’);form=1:n

fprintf(fid,’第%d个节点的电压相角时shita为(节点号从小到大排列):%f\n’,m,shita(m));

end

fprintf(fid,’================================================\n’);fprintf(fid,’各节点的功率S为(节点号从小到大排列):\n’);

form=1:n

fprintf(fid,’第%d个的复功率:%f+i*(%f)\n’,m,real(S(m)),imag(S(m))‘);

end

fprintf(fid,’================================================\n’);fprintf(fid,’节点有功功率P为:\n’);

form=1:n

fprintf(fid,’第