快速求取电力系统节点阻抗矩阵的高斯-约当消元新算法

摘要

本发明提供一种快速求取电力系统节点阻抗矩阵Z的高斯‑约当消元新算法，涉及电力系统分析计算领域，主要包括以下步骤：读取节点导纳矩阵Ｙ数据；将n阶Ｙ阵与n阶单位矩阵E构成n*2n阶特殊增广阵B＝YE；对B阵中Ｙ阵第1～n‑1列下三角元素消元，用四角规则计算Ｙ阵对角元素和上三角元素以及E阵对角元素和下三角元素，而通过每行规格化前的上三角元素赋值得到Ｙ阵的下三角元素，得B(n‑1)′＝Y(n‑1)′E(n‑1)′；用四角规则对B(n‑1)′＝Y(n‑1)′E(n‑1)′中Y(n‑1)′阵的上三角元素从第n～2列消元，仅计算E(n‑1)′阵的对角元素和下三角元素，得B2(n‑1)′＝Y2(n‑1)′E2(n‑1)′；根据E2(n‑1)′阵可直接得到Z阵的对角元素及下三角元素；根据对称性得Z阵的上三角元素；将Z阵写入数据文件。用本申请方法求取IEEE‑57～‑300节点系统Z阵，比传统高斯‑约当消元法计算速度高约65％。