下载文档原格式
电力系统中的潮流计算与优化方法潮流计算是电力系统运行和规划中的重要环节，它用于计算电力系统中各节点的电压、相角、有功、无功功率以及线路、变压器等的潮流分布情况。
对电力系统进行潮流计算可以帮助电力系统运行人员了解系统的稳定性、可靠性以及容载能力，也可以为电力系统规划提供数据支持。
本文将介绍电力系统潮流计算的基本方法与优化技术。
一、潮流计算的基本方法1.1 普通潮流计算方法潮流计算的基本方法是牛顿-拉夫逊迭代法（Newton-Raphson Iteration Method）和高尔顿法（Gauss-Seidel Method）。
牛顿-拉夫逊迭代法主要是通过不断迭代求解雅可比矩阵的逆，直到迭代误差小于给定阀值时停止迭代；高尔顿法则是逐一更新所有节点的电压与相角，直至所有节点的迭代误差都小于给定阀值。
1.2 快速潮流计算方法在大型电力系统中，普通的潮流计算方法计算速度较慢。
因此，研究人员提出了一些针对快速潮流计算的方法，如快速牛顿-拉夫逊法（Fast Newton-Raphson Method）和DC潮流计算方法。
快速牛顿-拉夫逊法通过简化牛顿-拉夫逊法的迭代公式，减少计算量，提高计算速度；DC潮流计算方法则是将潮流计算问题转化为一个线性方程组的求解问题，进一步提升计算效率。
二、潮流计算的优化技术2.1 改进的潮流计算算法为了提高潮流计算的准确性和收敛速度，研究人员提出了一些改进的潮流计算算法。
其中，改进的牛顿-拉夫逊法（Improved Newton-Raphson Method）是一种结合牛顿-拉夫逊法和割线法的算法，通过混合使用这两种方法，实现在减小迭代误差的同时加快计算速度。
此外，基于粒子群优化算法（Particle Swarm Optimization）和遗传算法（Genetic Algorithm）的潮流计算算法也得到了广泛研究和应用。
2.2 潮流优化潮流计算不仅可以用于分析电力系统的工作状态，还可以作为优化问题的约束条件。
潮流计算的基本算法及使用方法一、二、潮流计算的基本算法1.牛顿－拉夫逊法1．1 概述牛顿－拉夫逊法是目前求解非线性方程最好的一种方法。
这种方法的特点就是把对非线性方程的求解过程变成反复对相应的线性方程求解的过程，通常称为逐次线性化过程，就是牛顿－拉夫逊法的核心。
牛顿-拉夫逊法的基本原理是在解的某一邻域内的某一初始点出发，沿着该点的一阶偏导数——雅可比矩阵，朝减小方程的残差的方向前进一步，在新的点上再计算残差和雅可矩阵继续前进，重复这一过程直到残差达到收敛标准，即得到了非线性方程组的解。
因为越靠近解，偏导数的方向越准，收敛速度也越快，所以牛顿法具有二阶收敛特性。
而所谓“某一邻域”是指雅可比方向均指向解的范围，否则可能走向非线性函数的其它极值点，一般来说潮流由平电压即各母线电压(相角为0，幅值为1)启动即在此邻域内。
1．2一般概念对于非线性代数方程组 即()0,,,21=n i x x x f ()n i ,2,1=(1－1)在待求量x 的某一个初始计算值()0x 附件，将上式展开泰勒级数并略去二阶及以上的高阶项，得到如下的线性化的方程组()()()()()0000=∆'+x x f x f(1－2)上式称之为牛顿法的修正方程式。
由此可以求得第一次迭代的修正量()()()[]()()0100x f x f x -'-=∆ (1－3)将()0x ∆和()0x 相加，得到变量的第一次改进值()1x 。
接着再从()1x 出发，重复上述计算过程。
因此从一定的初值()0x 出发，应用牛顿法求解的迭代格式为()()()()()k k k x f x x f -=∆' (1－4)()()()k k k x x x ∆+=+1(1－5)上两式中：()x f '是函数()x f 对于变量x 的一阶偏导数矩阵，即雅可比矩阵J ；k 为迭代次数。
由式（1－4）和式子（1－5）可见，牛顿法的核心便是反复形成求解修正方程式。
电力系统潮流计算电力系统潮流计算是电力系统运行分析中的重要环节。
它通过对电力系统中各节点的电压、相角以及功率等参数进行计算和分析，从而得出电力系统的稳态运行状态。
本文将从潮流计算的基本原理、计算方法、应用及其发展等方面进行阐述。
一、潮流计算的基本原理电力系统潮流计算的基本原理是基于潮流方程建立的。
潮流方程是一组非线性的方程，描述了电力系统中各节点的电压、相角以及功率之间的关系。
潮流计算的目的就是求解这组非线性方程，以确定电力系统的电压幅值、相角及有功、无功功率的分布情况。
二、潮流计算的基本方法潮流计算的基本方法主要有直接法、迭代法以及牛顿-拉夫逊法。
直接法是通过直接求解潮流方程得到电力系统的潮流状况，但对于大规模复杂的电力系统来说，直接法计算复杂度高。
迭代法是通过对电力系统的节点逐个进行迭代计算，直到满足预设的收敛条件。
牛顿-拉夫逊法是一种较为高效的迭代法，它通过近似潮流方程的雅可比矩阵，实现了计算的高效和稳定。
三、潮流计算的应用潮流计算在电力系统运行与规划中起着重要作用。
首先，潮流计算可以用于电力系统的稳态分析，确定电力系统在各种工况下的电压、相角等参数，以判断电力系统是否存在潮流拥挤、电压失调等问题。
其次，潮流计算还可以用于电力系统的优化调度，通过调整电力系统的发电机出力、负荷组织等参数，以改善电力系统的经济性和可靠性。
此外，潮流计算还可以用于电力系统规划，通过对电力系统进行潮流计算，可以为新建电源、输电线路以及变电站等设备的规划和选择提供科学依据。
四、潮流计算的发展随着电力系统的规模不断扩大和复杂度的提高，潮流计算技术也得到了迅速的发展。
传统的潮流计算方法在计算效率和计算精度上存在一定的局限性。
因此，近年来研究者提出了基于改进的迭代方法、高精度的求解算法以及并行计算等技术，以提高潮流计算的速度和准确性。
此外，随着可再生能源的不断融入电力系统，潮流计算还需要考虑多种能源的互联互通问题，这对潮流计算提出了新的挑战，需要进一步的研究和改进。
潮流计算的基本算法及使用方法Company number：【0089WT-8898YT-W8CCB-BUUT-202108】潮流计算的基本算法及使用方法一、 潮流计算的基本算法1.牛顿－拉夫逊法1．1 概述牛顿－拉夫逊法是目前求解非线性方程最好的一种方法。
这种方法的特点就是把对非线性方程的求解过程变成反复对相应的线性方程求解的过程，通常称为逐次线性化过程，就是牛顿－拉夫逊法的核心。
牛顿-拉夫逊法的基本原理是在解的某一邻域内的某一初始点出发，沿着该点的一阶偏导数——雅可比矩阵，朝减小方程的残差的方向前进一步，在新的点上再计算残差和雅可矩阵继续前进，重复这一过程直到残差达到收敛标准，即得到了非线性方程组的解。
因为越靠近解，偏导数的方向越准，收敛速度也越快，所以牛顿法具有二阶收敛特性。
而所谓“某一邻域”是指雅可比方向均指向解的范围，否则可能走向非线性函数的其它极值点，一般来说潮流由平电压即各母线电压(相角为0，幅值为1)启动即在此邻域内。
1．2 一般概念对于非线性代数方程组即 ()0,,,21=n i x x x f ()n i ,2,1= (1－1)在待求量x 的某一个初始计算值()0x 附件，将上式展开泰勒级数并略去二阶及以上的高阶项，得到如下的线性化的方程组()()()()()0000=∆'+x x f x f (1－2)上式称之为牛顿法的修正方程式。
由此可以求得第一次迭代的修正量()()()[]()()0100x f x f x -'-=∆ (1－3)将()0x ∆和()0x 相加，得到变量的第一次改进值()1x 。
接着再从()1x 出发，重复上述计算过程。
因此从一定的初值()0x 出发，应用牛顿法求解的迭代格式为()()()()()k k k x f x x f -=∆' (1－4)()()()k k k x x x ∆+=+1 (1－5)上两式中：()x f '是函数()x f 对于变量x 的一阶偏导数矩阵，即雅可比矩阵J ；k 为迭代次数。
潮流计算的基本算法及使用方法一、二、潮流计算的基本算法1.牛顿－拉夫逊法1．1 概述牛顿－拉夫逊法是目前求解非线性方程最好的一种方法。
这种方法的特点就是把对非线性方程的求解过程变成反复对相应的线性方程求解的过程，通常称为逐次线性化过程，就是牛顿－拉夫逊法的核心。
牛顿-拉夫逊法的基本原理是在解的某一邻域内的某一初始点出发，沿着该点的一阶偏导数——雅可比矩阵，朝减小方程的残差的方向前进一步，在新的点上再计算残差和雅可矩阵继续前进，重复这一过程直到残差达到收敛标准，即得到了非线性方程组的解。
因为越靠近解，偏导数的方向越准，收敛速度也越快，所以牛顿法具有二阶收敛特性。
而所谓“某一邻域”是指雅可比方向均指向解的范围，否则可能走向非线性函数的其它极值点，一般来说潮流由平电压即各母线电压(相角为0，幅值为1)启动即在此邻域内。
1．2一般概念对于非线性代数方程组即 ()0,,,21=n i x x x f ()n i ,2,1= (1－1)在待求量x 的某一个初始计算值()0x 附件，将上式展开泰勒级数并略去二阶及以上的高阶项，得到如下的线性化的方程组()()()()()0000=∆'+x x f x f (1－2)上式称之为牛顿法的修正方程式。
由此可以求得第一次迭代的修正量()()()[]()()0100x f x f x -'-=∆ (1－3)将()0x ∆和()0x 相加，得到变量的第一次改进值()1x 。
接着再从()1x 出发，重复上述计算过程。
因此从一定的初值()0x 出发，应用牛顿法求解的迭代格式为()()()()()k k k x f x x f -=∆' (1－4)()()()k k k x x x ∆+=+1 (1－5)上两式中：()x f '是函数()x f 对于变量x 的一阶偏导数矩阵，即雅可比矩阵J ；k 为迭代次数。
由式（1－4）和式子（1－5）可见，牛顿法的核心便是反复形成求解修正方程式。
潮流计算的基本算法及使用方法一、二、潮流计算的基本算法1.牛顿－拉夫逊法1．1 概述牛顿－拉夫逊法是目前求解非线性方程最好的一种方法。
这种方法的特点就是把对非线性方程的求解过程变成反复对相应的线性方程求解的过程，通常称为逐次线性化过程，就是牛顿－拉夫逊法的核心。
牛顿-拉夫逊法的基本原理是在解的某一邻域内的某一初始点出发，沿着该点的一阶偏导数——雅可比矩阵，朝减小方程的残差的方向前进一步，在新的点上再计算残差和雅可矩阵继续前进，重复这一过程直到残差达到收敛标准，即得到了非线性方程组的解。
因为越靠近解，偏导数的方向越准，收敛速度也越快，所以牛顿法具有二阶收敛特性。
而所谓“某一邻域”是指雅可比方向均指向解的范围，否则可能走向非线性函数的其它极值点，一般来说潮流由平电压即各母线电压(相角为0，幅值为1)启动即在此邻域内。
1．2一般概念对于非线性代数方程组即 ()0,,,21=n i x x x f ()n i ,2,1= (1－1)在待求量x 的某一个初始计算值()0x 附件，将上式展开泰勒级数并略去二阶及以上的高阶项，得到如下的线性化的方程组()()()()()0000=∆'+x x f x f (1－2)上式称之为牛顿法的修正方程式。
由此可以求得第一次迭代的修正量()()()[]()()0100x f x f x -'-=∆ (1－3)将()0x ∆和()0x 相加，得到变量的第一次改进值()1x 。
接着再从()1x 出发，重复上述计算过程。
因此从一定的初值()0x 出发，应用牛顿法求解的迭代格式为()()()()()k k k x f x x f -=∆' (1－4)()()()k k k x x x ∆+=+1 (1－5)上两式中：()x f '是函数()x f 对于变量x 的一阶偏导数矩阵，即雅可比矩阵J ；k 为迭代次数。
由式（1－4）和式子（1－5）可见，牛顿法的核心便是反复形成求解修正方程式。
潮流计算的基本算法及使用方法一、欧阳歌谷（2021.02.01）二、潮流计算的基本算法1.牛顿－拉夫逊法1．1 概述牛顿－拉夫逊法是目前求解非线性方程最好的一种方法。
这种方法的特点就是把对非线性方程的求解过程变成反复对相应的线性方程求解的过程，通常称为逐次线性化过程，就是牛顿－拉夫逊法的核心。
牛顿-拉夫逊法的基本原理是在解的某一邻域内的某一初始点出发，沿着该点的一阶偏导数——雅可比矩阵，朝减小方程的残差的方向前进一步，在新的点上再计算残差和雅可矩阵继续前进，重复这一过程直到残差达到收敛标准，即得到了非线性方程组的解。
因为越靠近解，偏导数的方向越准，收敛速度也越快，所以牛顿法具有二阶收敛特性。
而所谓“某一邻域”是指雅可比方向均指向解的范围，否则可能走向非线性函数的其它极值点，一般来说潮流由平电压即各母线电压(相角为0，幅值为1)启动即在此邻域内。
1．2一般概念对于非线性代数方程组即()0,,,21=nixxxf ()ni,2,1= (1－1)在待求量x 的某一个初始计算值()0x 附件，将上式展开泰勒级数并略去二阶及以上的高阶项，得到如下的线性化的方程组()()()()()0000=∆'+x x f x f (1－2)上式称之为牛顿法的修正方程式。
由此可以求得第一次迭代的修正量()()()[]()()0100x f x f x -'-=∆ (1－3)将()0x ∆和()0x 相加，得到变量的第一次改进值()1x 。
接着再从()1x 出发，重复上述计算过程。
因此从一定的初值()0x 出发，应用牛顿法求解的迭代格式为()()()()()k k k x f x x f -=∆' (1－4)()()()k k k x x x ∆+=+1 (1－5)上两式中：()x f '是函数()x f 对于变量x 的一阶偏导数矩阵，即雅可比矩阵J ；k 为迭代次数。
由式（1－4）和式子（1－5）可见，牛顿法的核心便是反复形成求解修正方程式。
电力系统的潮流计算电力系统的潮流计算是电力系统分析中的基础工作，主要用于计算电力系统中各节点的电压和功率流动情况。
通过潮流计算可以得到电力系统的电压、功率、功率因数等关键参数，为电力系统的运行和规划提供有效的参考依据。
本文将介绍电力系统潮流计算的基本原理、计算方法和应用。
一、电力系统潮流计算的基本原理电力系统潮流计算基于电力系统的能量守恒原理和基尔霍夫电流定律，通过建立电力系统的节点电压和功率平衡方程组来描述系统中各节点间的电压和功率流动关系。
潮流计算的基本原理可简述为以下三个步骤：1.建立节点电压方程：根据基尔霍夫电流定律，将电力系统中各节点的电流状况表达为节点电压和导纳矩阵之间的乘积关系。
2.建立功率平衡方程：根据能量守恒原理，将电力系统中各支路的功率流动表达为节点电压和导纳矩阵之间的乘积关系。
3.解算节点电压：通过求解节点电压方程组，得到系统中各节点的电压值。
二、电力系统潮流计算的常用方法电力系统潮流计算常用的方法有高斯-赛德尔迭代法、牛顿-拉夫逊迭代法和快速潮流法等。
其中，高斯-赛德尔迭代法是一种基于节点电压的迭代算法，通过在每一次迭代中更新节点电压值来逐步逼近系统潮流平衡状态。
牛顿-拉夫逊迭代法是一种基于节点电压和节点功率的迭代算法，通过在每一次迭代中同时更新节点电压和节点功率值来逼近系统潮流平衡状态。
快速潮流法则是一种通过行列式运算直接求解节点电压的方法，对于大规模复杂的电力系统具有较高的计算效率和精度。
三、电力系统潮流计算的应用电力系统潮流计算在电力系统的规划和运行中有广泛应用。
具体应用包括：1.电力系统规划：通过潮流计算可以预测系统中各节点的电压和功率流动情况，为电力系统的设计和扩建提供参考依据。
2.电力系统稳定性分析：潮流计算可以帮助分析系统中节点电压偏差、功率瓶颈等问题，为系统的稳态和暂态稳定性分析提供基础数据。
3.运行状态分析：潮流计算可以实时监测系统中各节点的电压和功率流动情况，为电力系统的运行调度提供参考。
潮流计算的基本算法及使用方法一、二、潮流计算的基本算法1.牛顿－拉夫逊法1．1 概述牛顿－拉夫逊法是目前求解非线性方程最好的一种方法。
这种方法的特点就是把对非线性方程的求解过程变成反复对相应的线性方程求解的过程，通常称为逐次线性化过程，就是牛顿－拉夫逊法的核心。
牛顿-拉夫逊法的基本原理是在解的某一邻域内的某一初始点出发，沿着该点的一阶偏导数——雅可比矩阵，朝减小方程的残差的方向前进一步，在新的点上再计算残差和雅可矩阵继续前进，重复这一过程直到残差达到收敛标准，即得到了非线性方程组的解。
因为越靠近解，偏导数的方向越准，收敛速度也越快，所以牛顿法具有二阶收敛特性。
而所谓“某一邻域”是指雅可比方向均指向解的范围，否则可能走向非线性函数的其它极值点，一般来说潮流由平电压即各母线电压(相角为0，幅值为1)启动即在此邻域内。
1．2一般概念对于非线性代数方程组即 ()0,,,21=n i x x x f ()n i ,2,1= (1－1) 在待求量x 的某一个初始计算值()0x 附件，将上式展开泰勒级数并略去二阶及以上的高阶项，得到如下的线性化的方程组()()()()()0000=∆'+x x f x f (1－2)上式称之为牛顿法的修正方程式。
由此可以求得第一次迭代的修正量()()()[]()()0100x f x f x -'-=∆ (1－3)将()0x ∆和()0x 相加，得到变量的第一次改进值()1x 。
接着再从()1x 出发，重复上述计算过程。
因此从一定的初值()0x 出发，应用牛顿法求解的迭代格式为()()()()()k k k x f x x f -=∆' (1－4)()()()k k k x x x ∆+=+1 (1－5)上两式中：()x f '是函数()x f 对于变量x 的一阶偏导数矩阵，即雅可比矩阵J ；k 为迭代次数。
由式（1－4）和式子（1－5）可见，牛顿法的核心便是反复形成求解修正方程式。
潮流计算的基本算法及使用方法一、 潮流计算的基本算法1. 牛顿－拉夫逊法1．1 概述牛顿－拉夫逊法是目前求解非线性方程最好的一种方法。
这种方法的特点就是把对非线性方程的求解过程变成反复对相应的线性方程求解的过程，通常称为逐次线性化过程，就是牛顿－拉夫逊法的核心。
牛顿-拉夫逊法的基本原理是在解的某一邻域内的某一初始点出发，沿着该点的一阶偏导数——雅可比矩阵，朝减小方程的残差的方向前进一步，在新的点上再计算残差和雅可矩阵继续前进，重复这一过程直到残差达到收敛标准，即得到了非线性方程组的解。
因为越靠近解，偏导数的方向越准，收敛速度也越快，所以牛顿法具有二阶收敛特性。
而所谓“某一邻域”是指雅可比方向均指向解的范围，否则可能走向非线性函数的其它极值点，一般来说潮流由平电压即各母线电压(相角为0，幅值为1)启动即在此邻域内。
1．2 一般概念对于非线性代数方程组()0=x f即 ()0,,,21=n i x x x f ()n i ,2,1= (1－1)在待求量x 的某一个初始计算值()0x附件，将上式展开泰勒级数并略去二阶及以上的高阶项，得到如下的线性化的方程组()()()()()0000=∆'+x x f x f (1－2)上式称之为牛顿法的修正方程式。
由此可以求得第一次迭代的修正量()()()[]()()0100x f x f x -'-=∆ (1－3)将()0x ∆和()0x相加，得到变量的第一次改进值()1x 。
接着再从()1x 出发，重复上述计算过程。
因此从一定的初值()0x出发，应用牛顿法求解的迭代格式为()()()()()k k k x f x x f -=∆' (1－4)()()()k k k x x x ∆+=+1 (1－5)上两式中：()x f '是函数()x f 对于变量x 的一阶偏导数矩阵，即雅可比矩阵J ；k 为迭代次数。
由式（1－4）和式子（1－5）可见，牛顿法的核心便是反复形成求解修正方程式。
潮流计算的基本算法及使用方法一、 潮流计算的基本算法1. 牛顿－拉夫逊法1．1 概述牛顿－拉夫逊法是目前求解非线性方程最好的一种方法。
这种方法的特点就是把对非线性方程的求解过程变成反复对相应的线性方程求解的过程，通常称为逐次线性化过程，就是牛顿－拉夫逊法的核心。
牛顿-拉夫逊法的基本原理是在解的某一邻域内的某一初始点出发，沿着该点的一阶偏导数——雅可比矩阵，朝减小方程的残差的方向前进一步，在新的点上再计算残差和雅可矩阵继续前进，重复这一过程直到残差达到收敛标准，即得到了非线性方程组的解。
因为越靠近解，偏导数的方向越准，收敛速度也越快，所以牛顿法具有二阶收敛特性。
而所谓“某一邻域”是指雅可比方向均指向解的范围，否则可能走向非线性函数的其它极值点，一般来说潮流由平电压即各母线电压(相角为0，幅值为1)启动即在此邻域内。
1．2 一般概念 对于非线性代数方程组()0=x f即 ()0,,,21=n i x x x f Λ ()n i Λ,2,1= (1－1)在待求量x 的某一个初始计算值()0x附件，将上式展开泰勒级数并略去二阶及以上的高阶项，得到如下的线性化的方程组()()()()()0000=∆'+x x f x f (1－2)上式称之为牛顿法的修正方程式。
由此可以求得第一次迭代的修正量()()()[]()()0100x f x f x -'-=∆ (1－3)将()0x ∆和()0x相加，得到变量的第一次改进值()1x 。
接着再从()1x 出发，重复上述计算过程。
因此从一定的初值()0x出发，应用牛顿法求解的迭代格式为()()()()()k k k x f x x f -=∆' (1－4)()()()k k k x x x ∆+=+1 (1－5)上两式中：()x f '是函数()x f 对于变量x 的一阶偏导数矩阵，即雅可比矩阵J ；k 为迭代次数。
由式（1－4）和式子（1－5）可见，牛顿法的核心便是反复形成求解修正方程式。
牛顿法当初始估计值()0x和方程的精确解足够接近时，收敛速度非常快，具有平方收敛特性。
1．3 潮流计算的修正方程运用牛顿－拉夫逊法计算潮流分布时，首先要找出描述电力系统的非线性方程。
这里仍从节点电压方程入手，设电力系统导纳矩阵已知，则系统中某节点（i 节点）电压方程为∑=**•⎪⎪⎪⎭⎫ ⎝⎛=nj i i j ij U S U Y 1从而得∑=**••=nj j ij i i U Y U S 1进而有()01=-+*=*•∑jnj ij ii i UY U jQ P（1－6）式（1－6）中，左边第一项为给定的节点注入功率，第二项为由节点电压求得的节点注入功率。
他们二者之差就是节点功率的不平衡量。
现在有待解决的问题就是各节点功率的不平衡量都趋近于零时，各节点电压应具有的价值。
由此可见，如将式（1－6）作为牛顿－拉夫逊中的非线性函数()0=X F ，其中节点电压就相当于变量X 。
建立了这种对应关系，就可列出修正方程式，并迭代求解。
但由于节点电压可有两种表示方式——以直角做表或者极坐标表示，因而列出的迭代方程相应地也有两种，下面分别讨论。
1．3．1 直角坐标表示的修正方程节点电压以直角坐标表示时，令i i i jf e U +=•、j j j jf e U +=•，且将导纳矩阵中元素表示为ij ij ij jB G Y +=，则式（1－7）改变为 ()()()()01=--+-+∑=nj j j ij iji i i i jf e jB Gjf e jQ P（1－7）再将实部和虚部分开，可得()()[]()()[]⎪⎪⎭⎪⎪⎬⎫=+---=++--∑∑==0011nj j ij j ij i j ij j ij i i nj j ij j ij i j ij j ij i i e B f G e f B e G f Q e B f G f f B e G e P （1－8）这就是直角坐标下的功率方程。
可见，一个节点列出了有功和无功两个方程。
对于PQ 节点（1,,21-=m i Λ，），给定量为节点注入功率，记为i P '、i Q '，则由式（2－8）可得功率的不平衡量，作为非线性方程()()[]()()[]⎪⎪⎭⎪⎪⎬⎫+---'=∆++--'=∆∑∑==nj j ij j ij i j ij j ij i i i nj j ij j ij i j ij j ij i i i e B f G e f B e G f Q Q e B f G f f B e G e P P 11 （1－9）式中i P ∆、i Q ∆——分别表示第i 节点的有功功率的不平衡量和无功功率的不平衡量。
对于PV 节点（n m m i ,,2,1Λ++=），给定量为节点注入有功功率及电压数值，记为i P '、i U '，因此，可以利用有功功率的不平衡量和电压的不平衡量表示出非线性方程，即有()()[]()⎪⎭⎪⎬⎫+-'=∆++--'=∆∑=22221i i i i nj j ij j ij i j ij j ij i i i f e U U e B f G f f B e G e P P（1－10）式中i U ∆为电压的不平衡量。
对于平衡节点（m i =），因为电压数值及相位角给定，所以S s S jf e U +=•也确定，不需要参加迭代求节点电压。
因此，对于n 个节点的系统只能列出()12-n 个方程，其中有功功率方程()1-n 个，无功功率方程()1-m 个，电压方程()m n -个。
将式（1－9）、式（1－10） 非线性方程联立，称为n 个节点系统的非线性方程组，且按泰勒级数在()0i f 、()0i e （m i n i ≠=,,,2,1Λ）展开，并略去高次项，得到以矩阵形式表示的修正方程如下⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡∆∆∆∆∆∆∆∆⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡∆∆∆∆∆∆∆∆n n p p nn nnnp npn n n n nn nn np np n n n n pn pn pp pp p p p p pn pn pp pp p p p p n n p p n n p p n n p p n n p p n n pp e f e f e f e f S R S R S R S R N H N H N H N H S R S R S R S R N H N H N H N H L J L J L J L J N H N H N H N H L J L J L J L J N H N H N H N H U P U P Q P Q P M M M M 22112211221122112211222222222121222222222121111112121111111112121111222211 （1－11） 上式中雅可比矩阵的各个元素则分别为jiij f P H ∂∆∂=jiij e P N ∂∆∂=jiij f Q J ∂∆∂=jiij e Q L ∂∆∂=j i ij f U R ∂∆∂=2ji ij e U S ∂∆∂=2将（1－11）写成缩写形式[]⎥⎦⎤⎢⎣⎡∆∆=⎥⎦⎤⎢⎣⎡∆∆⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡∆∆∆e f e f S RL J N HU Q P J 2 （1－12）对雅可比矩阵各元素可做如下讨论：当i j ≠时，对于特定的j ，只有该特定点的i f 和i e 是变量，于是雅可比矩阵中各非对角元素表示为i ij i ij jiij f G e B f P H -=∂∆∂=i ij i ij jiij f B e G e P N --=∂∆∂=i ij i ij jiij e G f B f Q J +=∂∆∂= i ij i ij jiij e B f G e Q L +-=∂∆∂=02=∂∆∂=j i ij f U R02=∂∆∂=ji ij e U S当i j =时，雅可比矩阵中各对角元素的表示式为()i ii i ii nj j ij j ij j iij e B f G e B f G f P H +-+-=∂∆∂=∑=1()i ii i ii n j j ij j ij j iij f B e G f B e G e P N ----=∂∆∂=∑=1()f B e G f B e G f Q J ii i ii n j j ij j ij j iij ++--=∂∆∂=∑=1()i ii i ii n j j ij j ij j iij e B f G e B f G e Q L +-++=∂∆∂=∑=1i j i ij f f U R 22-=∂∆∂=i ji ij e e U S 22-=∂∆∂=由上述表达式可知，直角坐标的雅可比矩阵有以下特点：1)雅可比矩阵是()12-n 阶方阵，由于ji ij H H ≠、ji ij N N ≠等等，所以它是一个不对称的方阵。
2) 雅可比矩阵中诸元素是节点电压的函数，在迭代过程中随电压的变化而不断地改变。
3)雅可比矩阵的非对角元素与节点导纳矩阵B Y 中对应的非对角元素有关，当B Y 中的ij Y 为零时，雅可比矩阵中相应的ij H 、ij N 、ij J 、ij L 也都为零，因此，雅可比矩阵也是一个稀疏矩阵。
1．3．2 极坐标表示的修正方程在牛顿－拉夫逊计算中，选择功率方程∑=**•=-+nj jij ii i UY U jQ P 10作为非线性函数方程，把式中电压向量表示为极坐标形式()i i i j i i j U e U U i δδδsin cos +==•()j j j j j j j U eU U jδδδsin cos +==•则节点功率方程变为()()()0sin cos sin cos 1=--+-+∑=j j j nj ij ij i i i i i j U jB G j U jQ P δδδδ将上式分解成实部和虚部()0sin cos 1=+-∑=nj ij ij ij ij j i i B G U U P δδ()0cos sin 1=--∑=nj ij ij ij ij j i i B G U U Q δδ这就是功率方程的极坐标形式，由此可得到描述电力系统的非线性方程。
对于PQ 节点，给定了()()⎪⎪⎭⎪⎪⎬⎫--'=+-'=∆∑∑==nj ij ij ij ij j i i i nj ij ij ij ij j i i i B G U U Q Q B G U U P P 11cos sin sin cos δδδδ ()121-=m i Λ、、 （1－13）对于PV 节点，给定了i P '、i U '，而i Q '未知，式（1－13）中i Q ∆将失去作用，于是PV节点仅保留i P ∆方程，以求得电压的相位角。
