一、潮流计算概述、基本方法.ppt

一、潮流计算概述、基本方法.ppt

《一、潮流计算概述、基本方法.ppt》由会员分享，可在线阅读，更多相关《一、潮流计算概述、基本方法.ppt（109页珍藏版）》请在冰豆网上搜索。

电气工程学院电气工程学院电力系统稳态分析电气工程学院电气工程学院0、课程简介电气工程学院电气工程学院n电力系统潮流计算n基础知识n概述、潮流问题的数学模型nGeuss-Seidal法，N-R法n线性稀疏方程的解法nFDLF法n保留非线性潮流算法n最小化潮流算法n最优潮流问题n几个特殊性质的潮流计算简介主要内容

（一）主要内容

（一）主要内容

（一）主要内容

（一）电气工程学院电气工程学院n电力系统状态估计n概述n电力系统运行状态的表征与可观察性n最小二乘估计n不良数据的检测、不良数据的辩识n非二次准则的电力系统状态估计方法简介主要内容

（二）主要内容

（二）主要内容

（二）主要内容

（二）电气工程学院电气工程学院n电力系统静态安全分析n概述n电力系统静态等值n支络开断模拟n发电机开断模拟n预想事故的自动选择主要内容（三）主要内容（三）主要内容（三）主要内容（三）电气工程学院电气工程学院n电力系统复杂故障分析n简单故障的分析n用于故障分析的两口网络方程n复杂故障分析主要内容（四）主要内容（四）主要内容（四）主要内容（四）电气工程学院电气工程学院参考书目参考书目参考书目参考书目n电力系统分析诸骏伟水利电力n现代电力系统分析王锡凡方万良杜正春科学n电力系统稳态分析陈珩水利电力n电力系统静态安全分析吴际舜上海交大n电子数字计算机的应用电力系统计算西安交大等六院校合编水利电力n高等电力网络分析张伯明、陈寿荪清华大学n电力系统状态估计于尔铿水利电力n稀疏矩阵：

算法及程序实现杨绍祺等高等教育n线性优化及其扩展：

理论与方法方述诚等科学电气工程学院电气工程学院n中国电机工程学报n电力系统自动化n电网技术n电力自动化设备n继电器参考文献来源参考文献来源参考文献来源参考文献来源（国内国内国内国内）电气工程学院电气工程学院nIEEE（InstituteofElectricalandElectronicsEngineers）nPES（PowerEngineeringSociety）nIEEEtransactionsonPowerApparatusandSystems（PAS）1986年止，分为：

nIEEEtransactionsonPowerSystemsnIEEEtransactionsonPowerDeliverynIEEEtransactionsonEnergyConversionnIEE（TheinstitutionofElectricalEngineers）nPICA（PowerIndustryComputerApplication）nCIGRE（InternationalCouncilonLargeElectricSystems）参考文献来源（国际）参考文献来源（国际）参考文献来源（国际）参考文献来源（国际）电气工程学院电气工程学院电气工程学院电气工程学院电气工程学院电气工程学院一潮流计算概述、基本方法电气工程学院电气工程学院内容提要内容提要nn基础知识基础知识nn节点方程节点方程nn变压器等值电路变压器等值电路nn移相器数学模型移相器数学模型nn节点导纳矩阵节点导纳矩阵nn潮流计算概述潮流计算概述nn潮流计算问题的数学模型潮流计算问题的数学模型nn高斯赛德尔法高斯赛德尔法nn牛顿拉夫逊法牛顿拉夫逊法nn快速分解法快速分解法电气工程学院电气工程学院一、基一、基一、基一、基础础础础知知知知识识识识

（一）节点方程n分析交流电路有两种方法：

节点电压法和回路电流法n节点电压法比较普遍n以图示的两个电源，一个等值负荷系统为例说明节点方程n系统是5节点6支路n以地为参考，根据基尔霍夫第一定律，得到电气工程学院电气工程学院n以基尔霍夫第一定律可以列出节点电流方程：

n按节点电压整理后得到：

n左式中，左端是由各节点流出的电流，右端是向各节点注入的电流。

n左式可以表示为规范的形式电气工程学院电气工程学院n前述式子表示为规范形式如下：

n可以看出，其中的元素如下；n左式中，即为相应节点间的自导纳及互导纳。

其余节点间互导纳为零。

电气工程学院电气工程学院n重写规范形式如下：

n上式为电力网络的节点方程节点方程。

n在求出节点电压后，就可以求出各支路电流，从而使网络变量得以求解。

n节点方程反映了各节点电压与注入电流间的关系。

在此例中，除节点4、5外，其余节点注入电流均为0。

电气工程学院电气工程学院n一般情况下，如果电力网络有n个节点，则有节点方程：

n式中：

nY是导纳矩阵，对角元是节点i的自导纳，非对角元是节点间的互导纳。

分别是节点注入电流列向量及节点电压列向量节点方程节点方程节点方程节点方程电气工程学院电气工程学院n反映了电力网络的参数及接线情况n由导纳矩阵所构成的节点方程式是电力网络广泛应用的一种数学模型。

（二）节点导纳矩阵

（二）节点导纳矩阵

（二）节点导纳矩阵

（二）节点导纳矩阵电气工程学院电气工程学院n节点导纳物理意义：

n如果在节点i加一单位电压，而把其余节点全部接地节点导纳矩阵物理意义节点导纳矩阵物理意义节点导纳矩阵物理意义节点导纳矩阵物理意义电气工程学院电气工程学院n则上述节点方程式成为n节点自导纳Yii节点i加单位电压，其它节点接地时，节点i向电网注入的电流。

n节点互导纳Yji节点i加单位电压，其它节点接地时，节点j向电网注入的电流。

电气工程学院电气工程学院n例，有以下三节点网络n导纳矩阵有如下形式，现考虑如何求其中各元素举例举例电气工程学院电气工程学院n从图中可以看出：

n形成导纳阵第一列元素Y11，Y21，Y31。

应在节点1加单位电压，节点2、3接地。

举例电气工程学院电气工程学院n从图中可以看出：

n形成导纳阵第二列元素Y12，Y22，Y32。

应在节点2加单位电压，节点1、3接地。

举例电气工程学院电气工程学院n从图中可以看出：

n形成导纳阵第三列元素Y13，Y23，Y33。

应在节点3加单位电压，节点1、2接地。

举例电气工程学院电气工程学院n最后，得到该网络的导纳矩阵举例电气工程学院电气工程学院导纳矩阵物理意义（续）导纳矩阵物理意义（续）n令nk=i时，上式说明，当网络中除节点i以外所有节点都接地时，从节点i注入网络的电流同施加于节点i的电压之比，即节点自导纳Yii。

节点i加单位电压，其它节点接地时，节点i向电网注入的电流。

n自导纳Yii是节点i以外的所有节点都接地时节点i对地的总导纳。

显然，应等于与节点i相接的各支路导纳之和。

得电气工程学院电气工程学院导纳矩阵物理意义（续）导纳矩阵物理意义（续）nki时，上式说明，当网络中除节点k以外所有节点都接地时，从节点i注入网络的电流同施加于节点k的电压之比，即节点互导纳Yik。

节点k加单位电压，其它节点接地时，节点i向电网注入的电流。

n此时节点i的电流实际上是自网络流出并进入地中的电流，所以互导纳Yik应等于节点i，k间的支路导纳的负值。

电气工程学院电气工程学院n特点：

n当不含移相器时，导纳阵为对称矩阵n导纳矩阵为稀疏矩阵n出线数24条，每行非对角元中仅有24个非零元n例如，节点数分别10，1000的两个网络，平均出线为3n前者非零元40个，占总数40。

n后者非零元4000个，占总数0.4。

n计算时充分利用对称及稀疏性节点导纳矩阵特点与构成节点导纳矩阵特点与构成节点导纳矩阵特点与构成节点导纳矩阵特点与构成电气工程学院电气工程学院n阶数等于网络节点数n各行非对角元中非零元个数等于对应节点所连的不接地支路数n各对角元，即各节点的自导纳，等于相应节点所连支路的导纳之和n导纳矩阵非对角元素Yij等于节点i与j之间的支路导纳负值。

导纳矩阵导纳矩阵导纳矩阵导纳矩阵特点与特点与特点与特点与构成构成构成构成电气工程学院电气工程学院（三）变压器等值电路（三）变压器等值电路（三）变压器等值电路（三）变压器等值电路n忽略变压器励磁回路或作为负荷或阻抗单独处理时，变压器可以用漏抗串联一个无损耗理想变压器来模拟。

电气工程学院电气工程学院n由上式解得：

n写成：

电气工程学院电气工程学院n得变压器等值电路：

n或用相应导纳表示：

n其中，yT=1/zT电气工程学院电气工程学院n前述漏抗zT是放在变比为1的一侧n思考：

n如漏抗zT放在变比为K的一侧时，如何建立相应模型？

电气工程学院电气工程学院n漏抗zT放在变比为1的一侧时：

n如漏抗zT放在变比为K的一侧时，可以用：

n变化成：

电气工程学院电气工程学院（四）移相器（四）移相器（四）移相器（四）移相器n与变压器不同的是，移相器改变电压相位，因此，变比K是复数。

电气工程学院电气工程学院移相器数学模型移相器数学模型移相器数学模型移相器数学模型n要知道和的关系，要利用功率守恒原理。

n式中，、分别是和的共轭，从上式得到：

电气工程学院电气工程学院n最终有：

n其中：

n由于变比K为复数，造成Yij与Yji不等，因此，移相器没有相应等值电路。

而且，含有移相器的电力网络的导纳矩阵不对称。

电气工程学院电气工程学院二、潮流计算概述二、潮流计算概述二、潮流计算概述二、潮流计算概述n电力系统常规潮流计算：

根据给定的网络结构及运行条件，求出整个网络的运行状态。

n运行状态包括：

母线的电压、网络中的功率分布及功率损耗等。

电气工程学院电气工程学院n潮流计算分离线计算及在线计算n离线计算：

安排运行方式、规划供电方案、故障分析、优化计算n在线计算：

安全分析等潮流计算的应用场合潮流计算的应用场合潮流计算的应用场合潮流计算的应用场合电气工程学院电气工程学院n本质上是求解非线性代数方程本质上是求解非线性代数方程n在在数数学学上上一一般般都都是是采采用用非非线线性性代代数数方方程程求求解解，须须采用迭代法。

采用迭代法。

潮流计算的性质潮流计算的性质潮流计算的性质潮流计算的性质电气工程学院电气工程学院n潮流计算的基本要求：

潮流计算的基本要求：

（1）计算速度；）计算速度；（22）计算机内存使用量；）计算机内存使用量；（33）算法的收敛可靠性；）算法的收敛可靠性；（44）程序设计的方便性及算法扩充移植等的通用灵活性。

）程序设计的方便性及算法扩充移植等的通用灵活性。

潮流计算的基本要求潮流计算的基本要求潮流计算的基本要求潮流计算的基本要求电气工程学院电气工程学院n基本潮流算法：

基本潮流算法：

n高斯高斯-塞德尔法塞德尔法n牛顿法牛顿法n快速解耦法。

快速解耦法。

n几种改进算法：

几种改进算法：

n引入泰勒级数的高阶项，提高精度引入泰勒级数的高阶项，提高精度-保留非线性的保留非线性的潮流计算；潮流计算；n解决病态潮流解决病态潮流-最小化潮流计算法。

最小化潮流计算法。

n最优潮流问题：

兼顾电力系统的经济性、安全性和电能质最优潮流问题：

兼顾电力系统的经济性、安全性和电能质量。

量。

几种主要的潮流计算方法几种主要的潮流计算方法几种主要的潮流计算方法几种主要的潮流计算方法电气工程学院电气工程学院三、三、三、三、潮流计算问题的数学模型潮流计算问题的数学模型潮流计算问题的数学模型潮流计算问题的数学模型潮流计算所用的电力系统由变压器、输电线路、电容潮流计算所用的电力系统由变压器、输电线路、电容器、电抗器等静止线性元件所构成，并用集中参数表器、电抗器等静止线性元件所构成，并用集中参数表示的串联或并联等值支路来模拟。

示的串联或并联等值支路来模拟。

潮流计算普遍采用潮流计算普遍采用节点法节点法电气工程学院电气工程学院n上四式中：

上四式中：

Y，Z，Yij,Zij分别是节点导纳矩阵、节点阻分别是节点导纳矩阵、节点阻抗矩阵及其相应的元素；抗矩阵及其相应的元素；n为电力系统节点数。

为电力系统节点数。

或用用节点法节点法，节点电压与节点电流之间的关系，节点电压与节点电流之间的关系：

其展开式分别是：

电气工程学院电气工程学院n实际中只知道节点注入功率，节点电流和节点功率的关系式：

n将上式代入式（1-3）、（1-4）得到：

n或电气工程学院电气工程学院n重写（1-6）（1-7）式。

这就是潮流计算问题最基本的方程式，是一个以节点电压为变量的非线性代数方程组。

n由此可见，采用节点功率作为节点注入量是造成方程组呈非线性的根本原因。

电气工程学院电气工程学院n对于电力系统中的每个节点，要确定其运行状态，需要有四个变量；n有功注入P、无功注入Q、电压模值U及电压相角。

nn个节点总共有4n个运行变量要确定。

n式（1-6）或（1-7），总共包括2n个实数方程式，由此仅可以解得2n个未知运行变量。

n为此在潮流计算前必须将另外的2n个变量作为已知量而预先指定。

潮流计算中的方程个数与变量个数潮流计算中的方程个数与变量个数潮流计算中的方程个数与变量个数潮流计算中的方程个数与变量个数电气工程学院电气工程学院nPQ节点：

给出运行参数（P,Q），待求（V,）。

通常有变电所母线，某些出力P、Q给定的发电厂。

nPV节点：

给出（P,V），待求（Q,）。

必须有可调节无功电源，用于维持电压值。

通常选有一定无功功率储备的发电厂母线。

或有无功补偿设备的变电所。

nV节点或平衡节点：

系统中一般只设一个。

待求P,Q。

选调频发电厂母线，也可以为提高收敛性而选择出线最多的发电厂母线为平衡节点。

潮流计算中节点分类潮流计算中节点分类潮流计算中节点分类潮流计算中节点分类电气工程学院电气工程学院实际电力系统中的节点类型实际电力系统中的节点类型实际电力系统中的节点类型实际电力系统中的节点类型123454.过渡节点：

过渡节点：

PQ为为0的给定的给定PQ节点，如图中的节点，如图中的51.负荷节点：

给定功率负荷节点：

给定功率P、Q如图中的如图中的3、4节点节点2.发电机节点：

发电机节点：

如图中的节点如图中的节点1，可能有两种，可能有两种情况：

情况：

给定给定P、Q运行，给定运行，给定P、V运行运行3.负荷发电机混合节点负荷发电机混合节点：

PQ节点，如图中的节点，如图中的2发电机节点发电机节点负荷节点负荷节点负荷节点负荷节点混合节点混合节点过渡节点过渡节点电气工程学院电气工程学院潮流计算中节点类型的划分潮流计算中节点类型的划分潮流计算中节点类型的划分潮流计算中节点类型的划分3.平衡节点：

已知V、也称为松弛节点，摇摆节点12345平衡节点平衡节点PQ节点节点PQ节点节点PV节点节点PQ节点节点1.PQ节点：

已知P、Q负荷、过渡节点，PQ给定的发电机节点，大部分节点2.PV节点：

已知P、V给定PV的发电机节点，具有可调电源的变电所，少量节点电气工程学院电气工程学院PQV节点节点P节点节点4.P节点：

已知P5.PQV节点：

已知P、Q、VASVG6.V节点：

节点：

已知V8.PQV：

已知P、Q、V、7.Q节点：

节点：

已知Q电气工程学院电气工程学院例题：

例题：

例题：

例题：

IEEE22IEEE22IEEE22IEEE22节点类型划分节点类型划分节点类型划分节点类型划分平衡节点：

平衡节点：

PV节点：

节点：

PQ节点：

1）平衡节点从发电机节点中选择）平衡节点从发电机节点中选择2）除平衡机以外的发电机节点一般选作）除平衡机以外的发电机节点一般选作PV节点，节点，装有无功补偿装置的中间节点也可选作装有无功补偿装置的中间节点也可选作PV节点节点3）负荷节点和其它中间节点一般选作）负荷节点和其它中间节点一般选作PQ节点节点电气工程学院电气工程学院n交流电力系统中的复数电压变量可以用两种坐标形式来表示：

n或n而复数导纳为n由上述几式带入（1-6）可得到以下两种潮流方程。

电气工程学院电气工程学院n潮流方程的直角坐标形式:

n潮流方程的极坐标形式:

（i=1,2,3,n）（i=1,2,3,n）其中，表示标号为j的节点与i直接相连，并且包括j=i电气工程学院电气工程学院n而由某个电源发出的有功,无功功率则是由运行人员控制,是自变量或称为控制变量.n各个节点的电压模值或相角,则属于随着控制变量的改变而变化的因变量或状态变量.n若以p,u,x分别表示扰动变量、控制变量、状态变量,则潮流方程可以用下式表示f（x,u,p）=0（115）n根据上式,潮流计算的含义就是针对某个扰动变量p,根据给定的控制变量u,求出相应的状态变量x。

n每个节点的注入功率是该节点的电源输入功率和负荷需求功率的代数和.负荷需求的功率是取决于用户,称之为不可控变量或扰动变量.电气工程学院电气工程学院四、四、四、四、高斯高斯高斯高斯-塞德尔法塞德尔法塞德尔法塞德尔法n由式（1-6）可以得:

n式中:

Pis、Qis为节点给定的注入有功、无功功率。

n假定节点1为平衡节点,其给定电压为。

平衡节点不参加迭代。

于是对应这种情况的高斯-塞德尔迭代格式为:

计算Ui（k+1）时，用到了（2,i-1）的Uj（k+1），以及（i+1,n）的Uj（k）。

（i=2,3,.,n）电气工程学院电气工程学院n从一组假定的初值出发,依次进行迭代计算,迭代收敛的判据是（118）n当系统存在PV节点时，对应于这类节点的电压不修正。

并根据对应PV节点电压修正注入功率。

电气工程学院电气工程学院n高斯-塞德尔算法的优点:

n原理简单,程序设计十分容易。

线性非线性方程组均适用。

n导纳矩阵是一个对称且高度稀疏的矩阵。

n因此占用内存非常节省。

n每次迭代的计算量也小。

是各种潮流算法中最小的。

n高斯-塞德尔算法的缺点:

n收敛速度很慢。

（松散耦合）n迭代次数将随所计算网络节点数的增加而直线上升n病态条件的系统，计算往往会发生收敛困难。

电气工程学院电气工程学院为提高算法收敛速度,常用的方法是在迭代过程中加入加速因子a,即取式中：

是通过式1-16求得节点i电压的第（k+1）次迭代值；则是实际采用的节点i电压的第（k+1）次迭代值；a为加速因子，一般取1a2。

（119）电气工程学院电气工程学院n病态条件的系统。

n节点间相位角差很大的重负荷系统；n包含有负阻抗支路的系统；n具有较长的辐射型线路的系统；n长线路与短线路接在同一节点上，而且长短线路的长度比值又很大的系统。

病态条件的系统病态条件的系统病态条件的系统病态条件的系统n此外，平衡节点所在位置的不同选择也会影响收敛性能。

此外，平衡节点所在位置的不同选择也会影响收敛性能。

电气工程学院电气工程学院n为克服这些缺点，提出了基于节点阻抗矩阵的高斯赛德尔迭代法。

其迭代公式为:

基于节点阻抗矩阵的高斯赛德尔法基于节点阻抗矩阵的高斯赛德尔法基于节点阻抗矩阵的高斯赛德尔法基于节点阻抗矩阵的高斯赛德尔法（120）（121）电气工程学院电气工程学院n节点阻抗矩阵的高斯赛德尔法优点：

n算法的收敛速度比较快。

（紧密耦合）n达到收敛所需迭代次数与网络规模关系不大n主要缺点n阻抗矩阵所占用的内存量大n每次迭代的计算量也很大。

电气工程学院电气工程学院（一一）牛顿拉夫逊法的一般概念牛顿拉夫逊法的一般概念n要点：

把非线性方程的求解过程变成反复对相应的线性方程进行求解的过程，即通常所称的逐次线性化过程。

n对于非线性代数方程组f（x）=0即fi（x1,x2,xn）=0线性化方法：

在待求量x的某一个初始估计值x（0）附近，将上式展开成泰勒级数并略去二阶及以上的高阶项五、牛顿拉夫逊法五、牛顿拉夫逊法五、牛顿拉夫逊法五、牛顿拉夫逊法（122）（123）电气工程学院电气工程学院n对f（x）=0在x（0）附近泰勒展开级数并略去二阶及以上的高阶项，得到：

f（x（0）+f（x（0）x（0）=0（124）n上式称为牛顿法的修正方程式。

上式称为牛顿法的修正方程式。

n由此得到第一次迭代的修正量：

由此得到第一次迭代的修正量：

x（0）=f（x（0）-1f（x（0）（125）电气工程学院电气工程学院n从一定的初值x（0）出发，应用牛顿法求解的迭代格式为:

f（x（k）x（k）=f（x（k）x（k+1）=x（k）+x（k）nf（x）是函数f（x）对于变量x的一阶偏导数矩阵，即雅可比矩阵J；k为迭代次数。

电气工程学院电气工程学院几何认识几何认识几何认识几何认识又称切线法。

平方收敛性。

又称切线法。

平方收敛性。

下一步下一步迭代迭代第第k+1k+1步步迭代迭代电气工程学院电气工程学院例题例题电气工程学院电气工程学院以下讨论f（x）采用功率方程式模型，电压变量则采用极坐标和直角坐标的两种形式。

（1）修正方程的极坐标形式令，则采用极坐标形式的潮流方程是：

对每个PQ节点及PV节点，根据式1-13，有对每个PQ节点，根据式1-14，有（二二二二）牛顿潮流算法牛顿潮流算法牛顿潮流算法牛顿潮流算法的修正方程式的修正方程式的修正方程式的修正方程式电气工程学院电气工程学院n将上述方程式在某个近似解附近用泰勒级数展开，并略去二阶及以上的高阶项后，得到以矩阵形式表示的修正方程式为（注意下式子右侧有负号）n式中：

n为节点总数；m为PV节点数，雅可比矩阵是（2n-m-2）阶非奇异方阵。

（130）除以U只是为了雅可比矩阵元素的表达一致性好些电气工程学院电气工程学院计算计算时雅可比矩阵各元素时雅可比矩阵各元素极坐标雅可比矩阵各元素的表示式极坐标雅可比矩阵各元素的表示式极坐标雅可比矩阵各元素的表示式极坐标雅可比矩阵各元素的表示式（131）（132）（133）（134）电气工程学院电气工程学院计算计算i=j时雅可比矩阵各元素时雅可比矩阵各元素雅可比矩阵各元素的表示式雅可比矩阵各元素的表示式雅可比矩阵各元素的表示式雅可比矩阵各元素的表示式电气工程学院电气工程学院雅可比矩阵各元素的表示式雅可比矩阵各元素的表示式雅可比矩阵各元素的表示式雅可比矩阵各元素的表示式电气工程学院电气工程学院

（2）

（2）修正方程的修正方程的修正方程的修正方程的直角坐标形式直角坐标形式直角坐标形式直角坐标形式n对每个PQ节点，根据式1-11和1-12有n对每个PV节点，还有n令令，对每个节点都有二个方程式，所以在不计平衡节点，对每个节点都有二个方程式，所以在不计平衡节点方程式的情况下，总共有方程式的情况下，总共有2（n-1）个方程式。

个方程式。

电气工程学院电气工程学院（142）采用直角坐标形式的修正方程式为采用直角坐标形式的修正方程式为采用直角坐标形式的修正方程式为采用直角坐标形式的修正方程式为电气工程学院电气工程学院n雅可比矩阵各元素的表示式如下雅可比矩阵各元素的表示式雅可比矩阵各元素的表示式雅可比矩阵各元素的表示式雅可比矩阵各元素的表示式（143.48）电气工程学院电气工程学院n雅可比矩阵各元素的表示式如下雅可比矩阵各元素的表示式雅可比矩阵各元素的表示式雅可比矩阵各元素的表示式雅可比矩阵各元素的表示式电气工程学院电气工程学院雅可比矩阵各元素的表示式雅可比矩阵各元素的表示式雅可比矩阵各元素的表示式雅可比矩阵各元素的表示式n雅可比矩阵各元素的表示式如下（149.54）电气工程学院电气工程学院雅可比矩阵各元素的表示式雅可比矩阵各元素的表示式雅可比矩阵各元素的表示式雅可比矩阵各元素的表示式电气工程学院电气工程学院电气工程学院电气工程学院n以上极坐标及直角类型的修正方程式，有以下特点：

n修正方程式的数目分别为2（n-1）-m个及2（n-1）个，在PV节点所占的比例不大时，两者的方程式数目基本接近2（n-1）个。

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

n雅可比矩阵的非对角元是否为零决定于相应的节点导纳阵元素Yij是否为零。

n和节点导纳矩阵具有相同稀疏结构的分块雅可比矩阵在位置上对称，但雅可比矩阵不对称。

牛顿潮流算法修正方程式特点牛顿潮流算法修正方程式特点牛顿潮流算法修正方程式特点牛顿潮流算法修正方程式特点电气工程学院电气工程学院牛顿潮流算法修正方程式示例牛顿潮流算法修正方程式示例牛顿潮流算法修正方程式示例牛顿潮流算法修正方程式示例n示例系统：

6节点系统，3为PV节点，6为平衡节点。

n导纳矩阵结构：

电气工程学院电气工程学院n按节点号顺序而构成的分块雅可比矩阵和节点导纳矩阵有同样的稀疏结构。

电气工程学院电气工程学院n实用的牛顿法潮流程序中的程序特点主要有以下三方面n稀疏矩阵n“压缩”方式只储存其非零元素n只有非零元素才参加运算，n修正方程式的求解过程，采用对包括常数项的增广矩阵以按行消去法进行消元运算。

对增广矩阵边形成、边消元、边存储。

所需存储量是消元运算