一种基于A＝LR三角分解法快速求取电力系统节点阻抗矩阵的方法技术方案

一种基于LR三角分解法快速求取电力系统节点阻抗矩阵的方法，属于电力系统分析计算领域，主要步骤如下：形成节点导纳矩阵Y；对Y阵进行LR三角分解；用RZk＝Ek直接求取Zk阵对角元Zkk及以上的非对角元素；根据对称性求取对角元Zkk以左的非对角元素；写Z阵数据到数据文件。本发明专利技术利用了比LDU三角分解法计算效率更高的LR三角分解法来求取Z阵元素；利用单位矩阵E阵中Ek阵的结构特点，省去对W阵的计算过程，直接利用R阵求取Zk阵对角元及其以上的元素；并根据Z阵元素的对称性直接得到Z阵的下三角元素，进而大幅提高Z阵元素的求取速度。用本发明专利技术方法求解IEEE-30、-57、-118节点系统的Z阵元素，与传统的LDU三角分解法相比，计算速度大幅度提高。

【技术实现步骤摘要】
-种基于A = LR三角分解法快速求取电力系统节点阻抗矩 阵的方法
本专利技术属于电力系统分析计算领域，涉及一种求取电力系统节点阻抗矩阵的方 法。
技术介绍
节点阻抗矩阵Z在电力系统中应用十分广泛且具有重要的作用。传统的求取Z阵 的方法有A = LDU三角分解法、支路追加法、导纳矩阵Y消元求逆法等。在现有方法中，LDU 三角分解法在求取Z阵元素时计算速度最快，因而使用最多。该方法将nXn阶的Z阵分成 n个Zk阵，依次对每一列Zk阵的全部元素进行求解，在计算非对角元素时未利用Z阵元素 的对称性，导致大量计算时间的浪费。 此外，在各种三角分解法中，A = LR或A =⑶三角分解法形成因子阵所需计算元 素的个数最少，且只有1个中间矩阵变量W。而LDU三角分解法形成因子阵所需计算元素的 个数较多，且有2个中间矩阵变量W、H。因此LDU三角分解法的计算效率远比LR或⑶三角 分解法的计算效率要低。相应地，其计算速度比LR或CU三角分解法的计算速度约慢30% 以上。因此用LDU三角分解法求取Z阵元素并非是对三角分解法应用的最佳选择。 【专利技术内容】 本专利技本文档来自技高网...

【技术保护点】
一种基于A＝LR三角分解法快速求取电力系统节点阻抗矩阵的方法，其特征包括以下步骤：步骤1：读取各支路数据，形成节点导纳矩阵Y；步骤2：对Y阵进行A＝LR三角分解，得到L、R两个子阵：先根据YZk＝Ek，令Y＝LR，得LRZk＝Ek。再将LRZk＝Ek进一步分解为LWk＝Ek，RZk＝Wk二个方程；步骤3：用RZk＝Ek直接求取Zk阵对角元Zkk及以上的非对角元素；步骤4：将整个求出的Z阵写入数据文件。

【技术特征摘要】
1. 一种基于A = LR三角分解法快速求取电力系统节点阻抗矩阵的方法，其特征包括以 下步骤： 步骤1 :读取各支路数据，形成节点导纳矩阵Y ; 步骤2 :对Y阵进行A = LR三角分解，得到L、R两个子阵： 先根据 YZk = Ek，...

【专利技术属性】
技术研发人员：陈恳，席小青，万新儒，罗仁露，
申请(专利权)人：南昌大学，
类型：发明
国别省市：江西;36