现代电力系统分析(2011-3).pdfVIP

研究生学位课： 现代电力系统分析 任课教师：王继东 四、最优潮流计算的牛顿算法 最优潮流作为一个非线性规划问题，可以 利用非线性规划的各种方法来求解，更由 于结合了电力系统的固有物理特性，在变 量的划分、等式及不等式约束条件的处 理、有功与无功的分解、变量修正方向的 决定、甚至基本潮流计算方法的选择等等 方面，都可以有各种不同的方案。为此即 使是采用非线性规划方法，也曾出现过为 数甚多的最优潮流算法。其中，在1984年 由参考文献(Sun D I)提出的最优潮流牛顿 算法，得到了国内外学者的高度评价。 (一)牛顿法的基本原理 如同上面提到的梯度法或最速下降法， 牛顿法是另一种求无约束极值的方法。 设无约束最优化问题 f minx ( )x R ∈ n 其极值存在的必要条件 ▽f(x)=0，在一般 况下为一个非线性代数方程组。 现在用牛顿法对非线性代数方程组求解， 于是得到优化的迭代格式为 式中▽f(x(k))为目标函数f(x)的梯度向 量 ；k为迭代次数；H(x)= ▽2f(x)为目标函 数f(x)的海森矩阵，是目标函数对于x的二 阶导数，故牛顿法又称为海森矩阵法。 算法的收敛判据是，||▽f(x(k))|| ε 。 牛顿法在按上述的基本格式进行迭代时，其搜索 方向为 (( ))k 1 ( )k − k s H [x ( )]f x ( − ) ∇ 可见这种方法与最速下降法比较，除了利用了目 标函数的一阶导数之外，还利用了目标函数的二 阶导数，考虑了梯度变化的趋势，因此所得到的 搜索方向比最速下降法好，能较快地找到最优 点。 牛顿法在有一个较好的初值，并且H(x(k))为正定 的情况下，收敛速度极快，具有二阶收敛速度， 这是该法的突出优点。 但是牛顿法的使用也受到一些限制 ： (1)要求f(x)二阶连续可微； (2)每一步都要计算海森矩阵及其逆阵，内存量 和计算工作量都很大。为此,对于变量维数很高的 优化计算，实用上往往被迫转而采用不必直接求H 及其逆阵的拟牛顿法(变尺度法) 。 但是在有些情况下，海森矩阵是一个稀疏阵，于 是可以采用结合了稀疏矩阵技术的高斯消去法等 一整套极其有效的方法，直接求解修正方程以求 得△x，其计算效率极高。 而在电力系统最优潮流计算问题中，通过模型的 适当建立，相应的海森矩阵可以是一个高度稀疏 的矩阵，从而使海森矩阵法这种收敛速度极快的 方法完全可以在最优潮流计算这样的大规模非线 性规划问题中得到应用。而这正是下面要介绍的 牛顿最优潮流算法的最基本特色， (二)最优潮流牛顿算法 在最优潮流牛顿算法中，对变量不再区 分为控制变量及状态变量，而统一写为x， 这样便于构造稀疏的海森矩阵，优化是在 全空间中进行的。 于是最优潮流计算归结为如下非线性规 划问题 min f (x ) ⎫ ⎪⎪ . . ( ) 0 s t g x ⎬ ⎪ ( ) 0 h x ≤ ⎭⎪