发明名称 一种基于对称稀疏矩阵技术的改进LR三角分解求取电力系统节点阻抗矩阵的方法
摘要 一种基于对称稀疏矩阵技术的改进LR三角分解求取电力系统节点阻抗矩阵Z的方法,属于电力系统分析计算领域。包括以下步骤:读取数据文件;形成节点导纳矩阵Y;根据稀疏性和对称性对Y阵进行LR三角分解求R阵元素;根据R阵元素稀疏性求Z<sub>k</sub>阵对角元Z<sub>kk</sub>及以上元素;按对称性求Z<sub>kk</sub>以左元素;写Z阵数据到数据文件。本发明方法根据稀疏性和对称性按过程法对Y阵进行LR三角分解,仅求上三角的R阵元素,大幅提高三角分解速度;利用单位矩阵E阵结构特点省略W阵元素的计算,利用R阵元素的稀疏性直接解方程RZ<sub>k</sub>=E<sub>k</sub>,大幅提高回代求解速度。用本发明方法对IEEE-30、-57、-118节点系统进行验算,与传统的LR三角分解法相比,计算速度可提高约83~98%。
申请公布号 CN104572584A 申请公布日期 2015.04.29
申请号 CN201410790069.9 申请日期 2014.12.17
申请人 南昌大学 发明人 陈恳;万新儒;刘单;刘阳涵
分类号 G06F17/16(2006.01)I 主分类号 G06F17/16(2006.01)I
代理机构 南昌新天下专利商标代理有限公司 36115 代理人 施秀瑾
主权项 一种基于对称稀疏矩阵技术的改进LR三角分解求取电力系统节点阻抗矩阵的方法,其特征包括以下步骤:步骤1:读入n节点系统各线路支路数据文件;步骤2:形成节点导纳矩阵Y;步骤3:根据稀疏性和对称性对Y阵进行LR三角分解,但仅求R阵元素,并记录其非零元素位置;步骤4:根据RZ<sub>k</sub>=E<sub>k</sub>应用R阵元素的稀疏性回代求Z<sub>k</sub>阵对角元Z<sub>kk</sub>及以上的非对角元素;步骤5:根据对称性求对角元Z<sub>kk</sub>以左的非对角元素;步骤6:将Z阵写入数据文件。
地址 330031 江西省南昌市红谷滩新区学府大道999号