牛顿拉夫逊法.pptVIP

（一）牛顿-拉夫逊法在数学上是求解非线性代数方程式的有效方法。其要点是把非线性方程式的求解过程变成反复地对相应的线性方程式进行求解的过程，即通常所称的逐次线性化过程。 将非线性代数方程组 (1-22) 在待求量 的某一个初始估计值 附近，展开成泰勒级数并略去二阶及以上的高阶项，得到线性化方程组 (1-24) 称为牛顿法的修正方程式。 由上式根据初值 可求得第一次迭代的修正量 (1-25) 将 和 相加，得到变量的第一次改进值 。 因此，应用牛顿法求解的迭代格式为 (1-26) (1-27) 上两式中： 是函数 对于 的一阶偏导数矩阵，即雅可比矩阵 ，为迭代次数。 牛顿法当初值 和方程的精确解足够接近时，具有平方收敛特性。 （二）牛顿潮流算法的修正方程式 将牛顿法用于求解电力系统潮流计算问题时，由于所采用的数学表达式以及复电压变量采用的坐标形式的不同，可以形成牛顿潮流算法的不同形式。 以下讨论用得最广泛的 采用功率方程式模型，而电压变量分别采用极坐标和直角坐标的两种形式。 1 极坐标形式 令 ，对每个 节点及 节点，根据式(1-13)，有 (1-28) 对每个 节点，根据式(1-14)，有 (1-29) 将上述方程式在某个近似解附近用泰勒级数展开，略去二阶及以上的高阶项后，得到以矩阵形式表示的修正方程式 (1-30) 式中： 为节点个数， 为 节点数，雅可比矩阵是 阶非奇异方阵。 雅可比矩阵各元素的表示式如下： 2 直角坐标形式 令 ，此时每个节点，都有两个方程式。因此共有 个方程式。 对每个PQ 节点，根据式(1-11)和式(1-12)有： (1-39) (1-40) 对每个 节点，除了有与式(1-39)相同的有功功率方程式之外，还有 (1-41) 采用直角坐标形式的修正方程式为 (1-42) 雅可比矩阵各元素的表示式如下： 分析以上两种类型的修正方程式，可以看出两者具有以下的共同特点。 (1) 修正方程式的数目分别为 及 个，在 节点比例不大时，两者的方程式数目基本接近 个。 (2) 雅可比矩阵的元素都是节点电压的函数，每次迭代，雅可比矩阵都需要重新形成。 (3) 从雅可比阵非对角元素的表示式可见，某个非对角元素是否为零决定于相应的节点导纳矩阵元素 是否为零。如将修正方程式按节点号的次序排列，并将雅可比矩阵分块，把每个 阶子阵 作为分块矩阵的 元素，则按节点号顺序而构成的分块雅可比矩阵将和节点导纳矩阵具有同样的稀疏结构，是一个高度稀疏的矩阵。 (4) 和节点导纳矩阵具有相同稀疏结构的分块雅可比矩阵在位置上对称，但由于 ，所以雅可比矩阵不是对称阵。 修正方程式的这些特点决定了牛顿法潮流程序特点，在设计算法时应重点考虑。 （三）修正方程式的处理和求解 有效地处理修正方程式是提高牛顿法潮流程序计算速度并降低内存需求量的关键。 结合修正方程式的求解，目前实用的牛顿法潮流程序的程序特点主要有以下三个方面，这些程序特点对牛顿法潮流程序性能的提高起着决定性的作用。 1 对于稀疏矩阵，在计算机中只储存其非零元素，且只有非零元素才参加运算。 * 牛顿拉夫逊法 牛顿法解非线性方程 原理：将非线性方程线性化 —— Taylor 展开 取 x0 ? x*，将 f (x)在 x0 做一阶Taylor展开: ，? 在 x0 和 x* 之间