潮流计算的基本算法及使用方法 - 数据结构与算法.docx

多平衡机等运行参数。平衡发电机是电气岛内的电压相角参考点，当素可做如下讨论：ii元素表示为当ji时，雅可比矩阵中各对角元 多平衡机等运行参数。平衡发电机是电气岛内的电压相角参考点，当 素可做如下讨论：ii元素表示为当ji时，雅可比矩阵中各对角元 的修正方程节点电压以直角坐标表示时，令UeYGjB，则式（1 不需要参加迭代求节点电压。因此，对于n个节点的系统只能列出2 i 1 2 n 一、 潮流计算的基本算法 1. 牛顿－拉夫逊法 1．1 概述 牛顿－拉夫逊法是目前求解非线性方程最好的一种方法。这种方法的特点就是把对非线 性方程的求解过程变成反复对相应的线性方程求解的过程，通常称为逐次线性化过程，就 是牛顿－拉夫逊法的核心。 牛顿-拉夫逊法的基本原理是在解的某一邻域内的某一初始点出发，沿着该点的一阶偏 导数——雅可比矩阵，朝减小方程的残差的方向前进一步，在新的点上再计算残差和雅可 矩阵继续前进，重复这一过程直到残差达到收敛标准，即得到了非线性方程组的解。因为 越靠近解，偏导数的方向越准，收敛速度也越快，所以牛顿法具有二阶收敛特性。而所谓 “某一邻域”是指雅可比方向均指向解的范围，否则可能走向非线性函数的其它极值点， 一般来说潮流由平电压即各母线电压( 相角为0，幅值为1) 启动即在此邻域内。 1．2 一般概念 对于非线性代数方程组 即 f x ,x , ,x 0 i 1,2, n (1 － 1) 在待求量 x 的某一个初始计算值x 0 附件，将上式展开泰勒级数并略去二阶及以上的高 阶项，得到如下的线性化的方程组 f x 0 f x 0 x 0 0 (1 －2) 上式称之为牛顿法的修正方程式。由此可以求得第一次迭代的修正量 x 0 f x 0 1 f x 0 (1 －3) 将 x 0 和x 0 相加，得到变量的第一次改进值x 1 。接着再从x 1 出发，重复上述计算过 程。因此从一定的初值x 0 出发，应用牛顿法求解的迭代格式为 f x k x k f x k (1 －4) 解。因为越靠近解，偏导数的方向越准，收敛速度也越快，所以牛顿VB2－5）（2－6）（2－7）（2－8。2.快速分解法2．1概述快速分解法的基本思想是：把节点功率和K略去不计，即认为它们的元素都等于零。这样，n1m阶的方程n 解。因为越靠近解，偏导数的方向越准，收敛速度也越快，所以牛顿 VB2－5）（2－6）（2－7）（2－8 。2.快速分解法2．1概述快速分解法的基本思想是：把节点功率 和K略去不计，即认为它们的元素都等于零。这样，n1m阶的方程 n i i n i i ij ij ij i i i i ij jBij i j j j e jf 0 （ 1－7） 上两式中： f x 是函数 f x 对于变量 x 的一阶偏导数矩阵，即雅可比矩阵 J ；k 为迭代次 数。 由式（ 1－4 ）和式子（ 1－5 ）可见，牛顿法的核心便是反复形成求解修正方程式。牛 顿法当初始估计值 x 0 和方程的精确解足够接近时，收敛速度非常快，具有平方收敛特性。 1．3 潮流计算的修正方程 运用牛顿－拉夫逊法计算潮流分布时， 首先要找出描述电力系统的非线性方程。 这里仍 从节点电压方程入手，设电力系统导纳矩阵已知，则系统中某节点（ i 节点）电压方程为 从而得 S U Y U i i ij j j 1 进而有 P jQ U Y U i ij j j 1 0 （ 1－6） 式（ 1－6）中， 左边第一项为给定的节点注入功率，第二项为由节点电压求得的节点注 入功率。他们二者之差就是节点功率的不平衡量。 现在有待解决的问题就是各节点功率的 不平衡量都趋近于零时， 各节点电压应具有的价值。 由此可见， 如将式（ 1－6）作为牛顿－拉夫逊中的非线性函数 F X 0 ，其中节点电压 就相当于变量 X 。建立了这种对应关系， 就可列出修正方程式， 并迭代求解。 但由于节点 电压可有两种表示方式——以直角做表或者极坐标表示，因而列出的迭代方程相应地也有 两种，下面分