一种基于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三角分解法快速求取电力系统节点阻抗矩 阵的方法。本专利技术不但克服了现有技术的不足，而且还利用单位矩阵E阵中Ek阵的结构特 点，省去了对A = LR三角分解法中间矩阵变量W阵的计算过程，直接利用R阵和Ek阵求取 Zk阵对角元及其以上的元素。 本专利技术是通过以下技术方案实现的，主要步骤为： 步骤1 :读取各支路数据，形成节点导纳矩阵Y ; 步骤2 :对Y阵进行A = LR三角分解，得到L、R两个子阵； 步骤2中具体实施过程如下： 本专利技术利用了比LDU三角分解法计算效率更高的LR三角分解法来求取Z阵元素。 因此先根据YZk = Ek，令Y = LR，得LRZk = Ek。再将LRZk = Ek进一步分解为LWk = Ek，RZk=Wk两个方程。Wk阵是仅用于中间计算的过渡变量矩阵。 步骤3 :用RZk = Ek直接求取Zk阵对角元Zkk及以上的非对角元素； 步骤3中具体实施过程如下： (1)对方程LWk = Ek求解Wk阵的过程 传统的LR三角分解法中对方程LWk = Ek要求解整个Wk阵。 本专利技术由于对Zk阵元素仅需求解其对角元Zkk及以上元素的计算方式，因此对方 程LWk = Ek也仅需求解Wk阵对角元及其以上的元素，对应的也仅需使用Ek阵第k行对角元 及以上的元素。而Ek阵元素的特点为：第k行的元素为1，其余元素全部为零。这种仅计算 Wk阵对角元及其以上元素的算法,使得Wk = Ek成立，即所求得的Wk阵对角元及其以上的元 素与Ek阵对角元及其以上的元素完全相同。因此本专利技术方法中对方程LWk = Ek的求解完 全可省去。 ⑵对方程RZk = Wk求解Zk阵的过程 如果仅考虑对角元及其以上的元素，则有Wk = Ek成立，因此对方程RZk = Wk的求 解可直接转化成对方程RZk = Ek中对角元及其以上的元素的求解。对方程RZk = Ek求解Zk阵的计算顺序为：Zn，-_-，Zk，---,Z1，且在计算各个Z k阵中，仅计算对角元Zkk及其以上的元 素，并根据Z阵元素的对称性得到对角元Zkk以左的元素。与传统的LDU三角分解法求取Zk阵元素的方式相比，利用对称性计算Z阵下三角的非对角元素可以减少50%非对角元素的 计算。当然，如果需要，完全可以不求Z阵下三角的非对角元素，以进一步加快计算速度。 步骤4 :将整个求出的Z阵写入数据文件以备后续程序使用。 考虑到程序的结构化，所形成的Z阵数据文件可由下一个程序调用执行。 本专利技术方法与LDU三角分解法相比主要具有下述三个优势： (1)选用了计算结构更简洁、效率更高、计算速度更快的LR三角分解法； (2)根据对称性直接得出Z阵下三角中的非对角元素，省去50%非对角元素的计 算； (3)利用单位矩阵Ek阵元素结构的特点，可省去对W阵元素的计算过程，只有直接 利用R阵和Ek阵求取Zk阵对角元及其以上元素的计算过程，而LDU三角分解法求取Zk阵 元素则有3个计算过程和W、H 2个中间变量矩阵。 【附图说明】 图1为本专利技术方法求取Z阵元素的计算流程图。 图2为传统的LDU三角分解法求取Z阵元素的计算流程图。 【具体实施方式】 本专利技术将通过以下实施例作进一步说明。 实施例1。 以nXn阶节点系统为例，分别比较传统的LDU三角分解法和本专利技术方法求解Z阵 元素过程的不同。比较结果如表1所示。 表1传统的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