[发明专利]一种快速求解电力系统节点阻抗矩阵的方法

说明书

技术领域

本发明涉及电力系统分析计算领域，主要涉及一种快速求取节点阻抗矩阵的计算方法。

背景技术

在电力系统分析计算中，经常会用到节点阻抗矩阵，传统方法中常用LDU三角分解法和高斯消元法求解节点阻抗矩阵。

各种三角分解法均类似于因子表法，比较适合于系数矩阵不变方程的反复求解。许多文献均介绍用LDU三角分解法求解Y矩阵的逆矩阵Z，主要为了把对Z矩阵的求解转换成对Z_k矩阵的求解。由于三角分解法的计算公式均与高斯消元计算公式紧密关联，如果不考虑用于多个逆矩阵Z的反复求解(大多数情况下只求1个逆矩阵)，三角分解法在原理上不如高斯消元法简单直接、计算速度快，更不会比含规格化的高斯-约当消元法理想。只是由于传统的高斯消元法用于求解逆矩阵时一般未涉及规格化计算(计算流程图见图1)，而LDU三角分解法又包括了规格化计算(计算流程图见图2)，因此LDU三角分解法比不含规格化的高斯消元法计算速度更快。所以更多的文献在介绍求取逆矩阵时用LDU三角分解法，而不是高斯消元法。但传统的LDU三角分解法和高斯消元法都存在原理复杂，计算时间长等问题。

(1)高斯消元法或LDU三角分解法求取逆矩阵的方式

YZ＝ZY＝E  (1)

根据上式可得

YZ_k＝E_k  (k＝1,2,……,n)(2)

传统的高斯消元法或LDU三角分解法求取逆矩阵的方式上是用式(2)，即通过求取第1～n列的Z_k矩阵来获取Z矩阵，而不是同时完整的求出整个Z矩阵。由于LDU三角分解法含规格化计算，因此LDU三角分解法较高斯消元法具有明显计算速度上的优势。计算流程图分别见图1和图2。

发明内容

本发明的目的是提供了一种快速求解电力系统节点阻抗矩阵新方法，能提高电力系统分析计算中的计算速度。

本发明是通过以下技术方案实现的。

首先将节点导纳矩阵和单位矩阵形成特殊的增广矩阵，之后再对此增广矩阵进行含规格化的高斯-约当消元，即可求出方程的最终解，其基本步骤如下：

步骤1：输入节点导纳矩阵数据；

步骤2：节点导纳矩阵Y阵、单位方阵E阵构建增广矩阵B＝[Y E]；

步骤3：对B阵规格化及n次高斯-约当消元得B⁽ⁿ⁾″＝[E E⁽ⁿ⁾″]；

步骤4：得逆矩阵Z＝E⁽ⁿ⁾″；

步骤5：输出结果。

对增广矩阵B进行n次含规格化的高斯-约当后得到的新的矩阵E⁽ⁿ⁾″即为节点阻抗矩阵Z。

本发明所述的步骤2中，与传统的LDU三角分解法和高斯消元法求逆矩阵的方式不同。传统方法消元过程中形成的增广矩阵为B＝[Y E_k]，然后一列列的求取Z_k矩阵。本发明是利用高斯-约当消元法的计算原理对所需求逆阵的矩阵和单位矩阵形成特殊的增广矩阵B＝[Y E]，消元完成后可直接获得完整的逆矩阵Z，且无回代过程。

本发明所述的步骤3和步骤4中，对形成的特殊增广矩阵B进行含规格化的高斯-约当消元后可得B⁽ⁿ⁾″＝[Y⁽ⁿ⁾″ E⁽ⁿ⁾″]。此时Y阵变成了n阶单位方阵Y⁽ⁿ⁾″，而n阶单位方阵E变成了n阶方阵E⁽ⁿ⁾″阵，E⁽ⁿ⁾″阵就是所要求的节点阻抗矩阵Z。

由于本发明的整个Z矩阵是同时求得，因此无需用主要针对常系数方程多次求解的LDU三角分解法以及计算效率较低的不含规格化的高斯消元法。因此本发明可以快速求解节点阻抗矩阵，且原理简单，程序编写方便。用本发明对IEEE-57、-118、-300节点等系统进行计算，与传统的LDU三角分解法和高斯消元法求逆矩阵的方法相比，计算速度可分别提高约25～50％(见实施例1)。

本发明提出的一种快速求解电力系统节点阻抗矩阵新方法，利用了高斯-约当消元法在计算速度上的优势，对所求逆的矩阵和单位矩阵形成特殊的增广矩阵并进行消元。与传统的LDU三角分解法和高斯消元法求逆矩阵的方式不同，不是一列列的求取Z_k矩阵，而是同时解出整个Z矩阵，即本发明是直接用式(1)而不是用式(2)来完成求取逆矩阵，且无回代过程。这些特点决定了本方法的速度优势。计算流程图见见图3。