电力网络求解方法.docVIP

三角分解在电力网络方程求解中的应用 张庆华 （控制科学与工程学院，模式识别与智能系统专业，2011010236） 摘要：从电气角度来看，无论电力网络如何复杂，原则上都可以首先做出它的等值电路，然后用主流电路理论进行分析计算。电力网络方程主要用高斯消去法求解。高斯消去法求解线性方程组又由消去运算和回代运算两部分组成。应用矩阵论中的三角分解理论简化消去过程，提高了整体的计算速度。 电力网络方程 电力网络的数学模型是现代电力系统分析的基础。正常情况下的电力潮流和优化潮流分析、故障情况下短路电流计算以及电力系统静态安全分析和动态稳走性的评估，都离不开电力网络的数学模型。这里所谓电力网络，是指由输电线路、电力变压器、并(串)联电容器等静止元件所构成的整体【1】。从电气角度来看，无论电力网络如何复杂，原则上都可以首先做出它的等值电路，然后用主流电路理论进行分析计算。 电力网络通常是由相应的节点导纳矩阵或节点阻抗拒阵来描述的。在现代电力系统分析中，我们需要面对成千上万个节点及电力网络所连接的电力系统。对电力网络的描述和处理往往成为解决有关问题的关键【2】。应用矩阵论中的三角分解理论简化节点矩阵，能高速有效地处理电力网络方程，是电路分析计算的关键。 三角分解的应用 电力网络方程的形成 一般情况下，如果电力网络有n个节点则可用矩阵的形式表示为I=YV。 其中：I：节点注入电流列向量； Y：节点导纳矩阵； V：节点电压列向量； 式I=YV反映了各节点电压与注入电流的关系。左端为各节点流出的电流，右端为向各节点注入的电流。Y为nn的电力网络节点导纳矩阵。 高斯消去法 目前，对于已知节点导纳矩阵与节点电压列向量求解节点注入电流的电力网络方程主要用高斯消去法求解【3】。高斯消去法求解线性方程组由消去运算和回代运算两部分组成。消去运算又叫前代后算，可以按行进行，也可按列进行。同样，回代运算可以按行进行，也可按列进行。通常采用“消去运算按列进行，回代运算按行进行”的方式较多。 设有n阶线性方程组AX=B。其中距阵A和向量B的分别代表节点导纳矩阵与节点电压列向量，其可以是实数或复数。由于消去运算只对A和B进行，因此，为了算法叙述方便，把B作为第n+1列附在A之后，形成n×(n+1)阶增广矩阵。 经过对增广矩阵的n次消去运算，使矩阵A对角线以下的元素全部化为零，从而得到增广距阵： （1） 与之对应的方程是，它与原方程AX=B同解。从而求得各注入电流。 因子表与三角分解 但在实际计算中，经常遇到这种情况:对于方程组需要多次求解，每次仅改变其常数B，即测得的节点电压不同。而系数矩阵A（电力网络节点导纳矩阵）是不变的。这时，为了提高计算速度，可以利用因子表求解。 因子表可以理解为高斯消去法解线性方程组的过程中时常数项B全部运算的一种记录表格。如前所述，高斯消去法分为消去过程和回代过程。回代过程的运算由对系数矩阵进行消去运算后得到的上三角矩阵元素确定，见式(1)。为了对常数项进行消去运算(又叫前代运算)，还必须记录消去过程运算所需要的运算因子。消去过程中的运算又分规格化运算和消去运算，消去过程中对常数项B中的第i个元素的运算包括: （i=1,2,……,n） (2) (k=1,2,……,i-1) (3) 将上式中的运算因子,,……，及逐行放在下三角部分和式（1）的上三角矩阵元素合在一起，就得到了因子表： 厖  EMBED Equation.KSEE3 \* MERGEFORMAT 厖  EMBED Equation.KSEE3 \* MERGEFORMAT 厖  EMBED Equation.KSEE3 \* MERGEFORMAT 厖  EMBED Equation.KSEE3 \* MERGEFORMAT 厖 厖 厖 厖 厖 厖  EMBED Equation.KSEE3 \* MERGEFORMAT 厖 其中下三角元素用来对常数项B进行消去(前代）运算，上三角矩阵元素用来进行消去回代运算。因子表也可以表示为如下形式：