基于对称稀疏矩阵技术的分段对称反向高斯-约当消元法制造技术

本发明专利技术提出基于对称稀疏矩阵技术的分段对称反向高斯‑约当消元法，属于电力系统分析计算领域。主要步骤包括：读取Y阵及系统数据；构成n*2n阶特殊增广阵[YE]；判断并记录Y阵第i行非零元素信息于数组S

【技术实现步骤摘要】
基于对称稀疏矩阵技术的分段对称反向高斯-约当消元法
本专利技术属于电力系统分析计算领域，具体涉及一种基于对称稀疏矩阵技术的分段对称反向高斯-约当消元法，应用于求解电力系统节点阻抗矩阵，可用于电力系统等工程领域的相关计算。
技术介绍
传统高斯-约当消元法(约当法)派生于传统高斯消元法(高斯法)，均用于求解变系数而不是常系数方程，两者计算原理类似，但计算过程不同。理论计算中的约当法同时对上下三角元素消元，没有回代直接得到方程解，高斯法一般对下三角元素消元，用上三角元素回代得到方程解。表面上看，约当法的计算过程似乎更简单，但实际上由于约当法对上三角元素的消元计算比高斯法用上三角元素的回代计算更为复杂，因此约当法的计算速度总是低于高斯法。也由于约当法是同时对上下三角元素消元，因此对对称矩阵下三角元素消元时难以利用元素的对称性；对上三角元素消元时有大量的无效计算；对上下三角元素消元时均须应用计算公式，不利于计算过程的理解和编程；未将取倒后的对角元素作为规格化因子，造成大量的除法计算；对上下三角元素消元时均无法利用矩阵稀疏性等等。此外，其求解过程与高斯法一本文档来自技高网...

【技术保护点】
1.基于对称稀疏矩阵技术的分段对称反向高斯-约当消元法，其特征在于，包括以下步骤：/n步骤1：打开数据文件，读取Y阵及系统数据到相关数组；/n步骤2：将n阶Y阵与n阶E阵构成n*2n阶特殊增广阵[YE]；/n步骤3：判断Y阵第i行规格化前对角元以右的非零元素，在非零元素信息数组S

【技术特征摘要】
1.基于对称稀疏矩阵技术的分段对称反向高斯-约当消元法，其特征在于，包括以下步骤：
步骤1：打开数据文件，读取Y阵及系统数据到相关数组；
步骤2：将n阶Y阵与n阶E阵构成n*2n阶特殊增广阵[YE]；
步骤3：判断Y阵第i行规格化前对角元以右的非零元素，在非零元素信息数组Sij中记录非零元素列号和个数sij，并按对称性将第i行规格化前的非零元素赋值给第i列对角元以下的相应元素，再将对角元素取倒，规格化Y阵第i行对角元以右的非零元素和E阵第i行、第1～i列的元素；
(1)判断Y阵第i行规格化前对角元以右的非零元素，在Sij数组中记录其列号和个数sij，按对称性将第i行规格化前的非零元素赋值第i列对角元以下的相应元素；
(2)将对角元素取倒，规格化Y阵第i行对角元以右的非零元素和E阵第i行、第1～i列的元素；
步骤4：对Y阵下三角元素按第1～n-1列、从上往下正向按列消元，根据数组Sij所记录的第i行对角元以右非零元素的列号和个数sij，用四角规则计算对Y阵第i列...

【专利技术属性】
技术研发人员：陈恳，熊哲浩，郭甲宝，魏艺君，廖嘉文，熊守江，
申请(专利权)人：南昌大学，
类型：发明
国别省市：江西;36