一种基于对称稀疏矩阵技术的LR三角分解求取电力系统节点阻抗矩阵的方法技术方案

一种基于对称稀疏矩阵技术的LR三角分解求取电力系统节点阻抗矩阵Z的方法，属于电力系统分析计算领域。主要包括以下步骤：读入n节点系统各线路支路数据；形成节点导纳矩阵Y；根据Y阵元素的稀疏性和对称性对Y阵进行LR三角分解；按对称性回代求取Z阵元素；写Z阵数据到数据文件。本发明专利技术根据LR三角分解法元素结构的特点，提出LR合成矩阵的概念，方便对计算过程的理解而且节省存贮单元；采用对称稀疏矩阵技术不但省去了大量元素的计算，而且可简化所有lij元素的计算，继而大大提高了Z阵元素的求取速度。本发明专利技术方法原理简单、计算快捷。用本发明专利技术对IEEE-30、-57、-118节点系统进行验算，与不考虑稀疏性和对称性的LR三角分解法相比，计算速度提高约45％。

【技术实现步骤摘要】

本专利技术属于电力系统分析计算领域，涉及一种电力系统节点阻抗矩阵的方法。
技术介绍
在电力系统中用三角分解法求取节点阻抗矩阵Z时一般都使用A＝LDU三角分解法，但实际上，由于计算过程和计算变量的不同，LDU三角分解法的计算效率远低于LR三角分解法，因此用LR三角分解法求取Z阵元素应该是更好的选择。但在传统的LR三角分解法中，一般均未考虑利用元素的稀疏性和对称性，从而导致大量零元素和部分非零元素的计算，使计算效率大大降低。电力系统计算中稀疏矩阵技术运用很广，主要为省去大量零元素的存贮及计算，加快高斯消元法的计算速度。矩阵元素的存贮方案也很多，如按坐标存贮、按顺序存贮、按链表存贮等等。尽管这些存贮方式可以省去不少存贮单元，但计算速度并没有达到最优效果，而且这些存贮方式结构复杂，且对角元素与非对角元素分开存贮也使得存取过程繁琐，特别不利于对称矩阵中的数据处理。实际上，这些存贮方式主要为减少存贮单元，对存贮过程的简化或存贮速度的提高并没有特别优势。而且这些存贮方式主要用于高斯消元法中，极少用于三角分解法中。由于传统的稀疏矩阵技术一般不考虑矩阵元素结构的特点对非零元素进行存贮，这种存贮方式在进行LR三角分解计算时无法利用L、R因子阵元素的稀疏性、对称性及其相互间的关系等特点。因此如果将传统的稀疏矩阵技术用于三角分解法中，其计算过程较繁琐复杂、计算速度较慢、计算效果也并不理想。
技术实现思路
为了克服上述现有技术的不足，本专利技术提供了一种基于对称稀疏矩阵技术的LR三角分解快速求取电力系统节点阻抗矩阵的方法。本专利技术是通过以下技术方案实现的，主要包括以下步骤：步骤1：读入n节点系统各线路支路数据；步骤2：形成节点导纳矩阵Y；步骤3：根据Y阵元素及其在三角分解过程中的稀疏性和对称性结构进行A＝LR三角分解；步骤3中具体实施过程如下：(1)提出合成矩阵的概念。下图左侧为4阶Y阵，在此基础上建立L、R二个因子阵，其中L阵是单位下三角阵，lii＝1；R阵是上三角阵，rii≠1。根据其因子阵结构的特点，可建立下图右侧的4阶合成矩阵。Y11Y12Y13Y14Y21Y22Y23Y24Y31Y32Y33Y34Y41Y42Y43Y44,r11r12r13r14l21r22r23r24l31l32r33r34l41l42l43r44]]>(2)合成阵中，L、R阵各元素与Y阵元素的关系如下。r11=Y11r12=Y12r13=Y13r14=Y14l21=Y21/r11r22=Y22-l21r12r23=Y23-l21r13r24=Y24-l21r14l31=l32=r33=r34=Y31/r11(Y32-l31r12)/r22Y33-(l31r13+l32r23)Y34-(l31r14+l32r24)l41=l42=l43=r44=Y41/r11(Y42-l41r12)/r22[Y43-(l41r13+l42r23)]/r33Y44-(l41r14+l42r24+l43r34)]]>(3)从合成阵可以看出，L、R阵各元素之间在任何情况下均有规律性的对应关系：如lji＝rij/rii，这使得任一对角元以右的rij元素与对角元以下的lji元素尽管数值不等，但非零元素的位置对称。利用此特点，只需计算对角元以右非零的rij元素，就可根据lji＝rij/rii关系得到对角元以下非零的lji元素，从而大大简化对非零的lji元素的计算。(4)从合成阵可以看出，在同一行的lij、rii、rij元素之间，所有元素的计算都类似于对对角元左侧的lij元素逐个进行消元，只是消元过程结束后对角元左侧的各个lij元素均要除以其所在列的对角元素rjj，而对角元及其右侧的rij元素不用。(5)假设l31≠0，则对r11元素之下的l31元素要进行消元。先定义r11元素为对角元素，r11右侧的所有r1j元素均定义为交叉元素，l31元素定义为消元元素，则l31右侧的所有元素均定义为计算元素。消元过程中若不考虑元素的稀疏性，则要计算l31元素右侧的所有元素l32、r33、r34；若采用稀疏性技术，则仅需计算不为零的l31元素所在行与不为零的r1j元素所在列相交的计算元素，从而大大减少相应元素的计算。(6)当对角元左侧的l31、l32元素被消元后，将r33右侧非零的r34元素除以其所在行的对角元素r33，并直接赋值给r33以下的l43元素。此时(5)和(6)所采用的就是对称稀疏性技术。该技术不但利用元素的稀疏性省去了大量元素的计算，而且利用了元素的对称性，可简化所有lij元素的计算。依次循环，完成合成阵所有元素的计算，得到L、R因子阵的合成阵。步骤4：根据LRZk＝Ek回代求取Zk阵对角元Zkk及以上的非对角元素；步骤5：根据对称性求取对角元Zkk以左的非对角元素；步骤6：将Z阵写入数据文件以备后续程序使用。考虑到程序的结构化，所形成的Z阵数据文件可由下一个程序调用执行。本专利技术方法具有以下几点优点。(1)三角分解后的L、R阵放在一个合成阵中不但方便对计算过程的理解而且节省存贮单元。(2)三角分解过程中，合成阵中的L、R阵非零元素位置始终对称，可省去大量判断语句；且在求得对角元以右的rij元素后除以其所在行相应的对角元素rii，就可得对角元以下非零的lji元素，可简化所有lij元素的计算。(3)三角分解过程中，只需计算非零的lij元素所在行与非零的rij元素所在列相交的元素，从而大大减少计算元素的计算量，快速完成前代过程。附图说明图1为传统的不考虑元素稀疏性和对称性的LR三角分解法求取Z阵元素流程图图2为本专利技术方法求取Z阵元素流程图具体实施方式本专利技术将通过以下实施例作进一步说明。实施例1。分别用传统的不考虑元素稀疏性和对称性的LR三角分解法(图1)以及本专利技术方法(图2)对IEEE-30、-57、-118节点系统的Y阵求其Z阵元素，并比较其“分解”和“分解+回代”过程的平均计算时间。计算结果如表1所示。表1稀疏性及对称性对LR三角分解法“分解”和“分解+回代”过程时间的影响T1：传统的LR三角分解法在“分解”过程的平均计算时间T2：本专利技术方法在“分解”过程的平均计算时间T11：传统的LR三角分解法在“分解+回代”过程的平均计算时间T21：本专利技术方法在“分解+回代”过程的平均计算时间根据表1可以看出：无论是“分解”过程还是“分解+回代”过程，本专利技术方法的计算速度大大优于传统方法的计算速度。以IEEE-118节点系统为例，在“分解”过程中，本专利技术方法的平本文档来自技高网...

【技术保护点】
一种基于对称稀疏矩阵技术的LR三角分解求取电力系统节点阻抗矩阵的方法，其特征包括以下步骤：步骤1：读入n节点系统各线路支路数据；步骤2：形成节点导纳矩阵Y；步骤3：根据Y阵元素及其在三角分解过程中的稀疏性和对称性结构进行A＝LR三角分解；步骤4：根据LRZk＝Ek回代求取Zk阵对角元Zkk及以上的非对角元素；步骤5：根据对称性求取对角元Zkk以左的非对角元素；步骤6：将Z阵写入数据文件。

【技术特征摘要】
1.一种基于对称稀疏矩阵技术的LR三角分解求取电力系统节点阻抗矩阵的方法，
其特征包括以下步骤：
步骤1：读入n节点系统各线路支路数据；
步骤2：形成节点导纳矩阵Y；
步骤3：根据Y阵元素及其在三角分解过程中...

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