[发明专利]一种基于A＝LDU三角分解法求解电力系统节点阻抗矩阵的方法

摘要

一种基于A＝LDU三角分解法求解电力系统节点阻抗矩阵的方法，属于电力系统分析计算领域。主要包括以下步骤：形成节点导纳矩阵Y，对Y阵进行LDU三角分解；对DH_k＝W_k(E_k)方程求取h_kk元素；对UZ_k＝H_k求取Z_k阵对角元Z_kk及以上的非对角元素；根据对称性求对角元Z_kk以左的非对角元素；写Z阵数据到数据文件。本发明方法主要利用了单位矩阵E阵中E_k阵的结构特点、Z_k阵元素的计算顺序以及Z阵元素的对称性，省去了W阵的计算，将对H阵的计算简化为仅对其对角元素h_kk的计算，对Z阵可省去50％非对角元素的计算，大幅度提高了Z阵元素的计算速度。用本发明方法对IEEE-57、-118、-300节点系统进行验算，与传统的LDU三角分解法相比，本发明方法计算速度可提高约35～45％。

主权项

一种基于A＝LDU三角分解法求解电力系统节点阻抗矩阵的方法，其特征是包括以下步骤：步骤1：读取各支路数据，形成节点导纳矩阵Y；步骤2：对节点导纳矩阵Y进行LDU三角分解，得到L、D、U三个因子阵：根据YZ_k＝E_k，令Y＝LDU，得LDUZ_k＝E_k，再将LDUZ_k＝E_k进一步分解为LW_k＝E_k，DH_k＝W_k，UZ_k＝H_k三个方程；步骤3：对DH_k＝W_k(E_k)方程，仅求取H_k阵中的对角元素h_kk，通过UZ_k＝H_k求取Z_k阵对角元Z_kk及以上的非对角元素，利用对称性得到对角元Z_kk以左的非对角元素；步骤4：将Z阵写入数据文件。