CN104714928B - 一种基于对称稀疏矩阵技术的高斯消元法求取电力系统节点阻抗矩阵的方法 - Google Patents
一种基于对称稀疏矩阵技术的高斯消元法求取电力系统节点阻抗矩阵的方法 Download PDFInfo
- Publication number
- CN104714928B CN104714928B CN201510026804.3A CN201510026804A CN104714928B CN 104714928 B CN104714928 B CN 104714928B CN 201510026804 A CN201510026804 A CN 201510026804A CN 104714928 B CN104714928 B CN 104714928B
- Authority
- CN
- China
- Prior art keywords
- battle array
- sparsity
- diagonal
- symmetrical
- diagonal element
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Complex Calculations (AREA)
Abstract
一种基于对称稀疏矩阵技术的高斯消元法求取电力系统节点阻抗矩阵的方法,属于电力系统分析计算领域。主要包括以下步骤:形成节点导纳矩阵Y;将Y阵与En阵形成增广阵Bn=[YEn];根据对称稀疏性对Bn阵消元得Bn (n‑1)′=[Y(n‑1)′En (n‑1)′];根据Y(n‑1)′Zn=En (n‑1)′、稀疏性、对称性求Zn阵对角元Znn以上和以左元素;根据Y(n‑1)′阵得Y(k‑1)′阵;根据Y(k‑1)′Zk=Ek (k‑1)′、稀疏性、对称性求Zk阵对角元Zkk以上和以左元素。本发明方法利用对称稀疏性,避免了前代过程的所有无效计算,减少了约50%非零元素的计算;利用E阵元素结构的特点以及上三角元素的稀疏性按对称方式回代求取Zk阵元素,大大加速了回代计算。用本发明方法对IEEE‑30、‑57、‑118节点等系统进行验算,与传统的高斯消元法和LDU三角分解法相比,对IEEE‑118节点系统计算速度可提高96~97%。
Description
技术领域
本发明属于电力系统分析计算领域,涉及一种求取电力系统节点阻抗矩阵的方法。
背景技术
电力系统中一般都用传统的LDU三角分解法、个别文献也介绍无规格化的高斯消元法求取节点阻抗矩阵Z。这2种方法求解Z阵时一般是将对n*n阶Z阵的求解转换成对n个Zk阵(Z1~Zn)整列元素的求解,未用到Z阵元素的对称性进行计算,即未仅计算Zk阵中对角元Zkk及以上的元素Zk-1,k~Z1k,再根据Zk-1,k~Z1k直接获得对角元以左的元素Zk,k-1~Zk1。因此计算整列Zk阵元素的方式比用对称性计算要多计算约50%的Zk阵元素。
LDU三角分解法适合于求解常系数方程,但计算原理及计算过程均比高斯消元法复杂。由于LDU三角分解法的计算过程中含有规格化的因子矩阵,因此求解Z阵元素时其计算效率高于不含规格化的高斯消元法。实际上,含规格化的高斯消元法比不含规格化的高斯消元法的计算速度可高约30%。因此,如果用含规格化的高斯消元法将对一个n*n阶Z阵的求解转换成对n个Zk阵的求解,则其计算速度应优于LDU三角分解法。
此外,由于节点导纳矩阵Y阵是对称稀疏矩阵,传统方法在消元过程中未有效地利用Y阵的对称性和稀疏性,在回代过程中未考虑E阵元素结构的特点和消元后矩阵上三角元素的稀疏性,从而使计算效率更是大大降低。
电力系统计算中稀疏矩阵技术运用很广,主要为省去大量零元素的存贮及计算,加快高斯消元法的计算速度。矩阵元素的存贮方案也很多,如按坐标存贮、按顺序存贮、按链表存贮等等。尽管这些存贮方式可以省去不少存贮单元,但计算速度并没有达到最优效果,而且这些存贮方式结构复杂,且对角元素与非对角元素分开存贮也使得存取过程繁琐,不利于对称矩阵中数据的处理。实际上,这些存贮方式主要为减少存贮单元,对计算过程的简化或计算速度的提高并无优势。而且这些存贮方式主要用于高斯消元法中,在三角分解法中的应用较为复杂。且由于传统的稀疏矩阵技术一般不考虑矩阵元素结构的特点对非零元素进行存贮,需形成另外的存储矩阵,因此原理复杂,计算速度慢。
发明内容
为了克服上述现有技术的不足,本发明提供了一种基于对称稀疏矩阵技术的高斯消元法求取电力系统节点阻抗矩阵的方法。
本发明是通过以下技术方案实现的,主要包括以下步骤:
步骤1:读入n节点系统各线路支路数据;
步骤2:形成节点导纳矩阵Y;
步骤3:Y阵与E阵最后一列形成增广阵Bn=[Y En];
步骤4:根据对称稀疏性对Bn阵进行n-1次含规格化的高斯消元得Bn (n-1)′=[Y(n-1)′En (n-1)′],并记录Y(n-1)′阵上三角中非零元素的位置;
步骤4中具体实施过程如下:
(1)用对称稀疏矩阵技术快速消元,同时记录其上三角非零元素的位置。
消元前的Y阵如左下式,进行含规格化的消元后,Y阵变成Y(n-1)′阵如右下式。
假设Y31≠0,则对Y31元素要进行消元。定义:Y11元素为对角元素;Y11右侧的所有元素Y1j均为交叉元素;Y31元素为消元元素;Y31右侧的所有元素Y3j为计算元素。
1)根据对角元以右非零的交叉元素的值和位置用对称性获得对角元以下非零的消元元素的值和位置,并仅计算其行列交互点上的计算元素。
Y(n-1)′阵中,每行元素在消元前后、规格化之前,其对角元以右非零的交叉元素与对角元以下非零的消元元素数值相等、位置对称。而在规格化之后,对角元以右非零的交叉元素均除以了对角元素,此时与对角元以下非零的消元元素数值不等(只相差对角元的数值)、但位置仍然对称。
因此可根据对角元以右非零的交叉元素的值和位置用对称性获得对角元以下非零的消元元素的值和位置,并利用非零的交叉元素所在行以及非零的消元元素所在列的交互点来确定需计算的计算元素。这种用非零的交叉元素根据对称性来确定非零的消元元素、再用它们行与列的交互点来确定所需计算的计算元素的方式可大大减少对元素的无效减少,加快高斯消元的计算速度。
2)按对称稀疏性,仅计算非零的消元元素所在行其对角元及其以右和非零的交叉元素所在列交互点上的计算元素。
假设Y31≠0。
不考虑元素的稀疏性:要计算Y31元素右侧所有的交叉元素Y32、Y33、…;
考虑元素的稀疏性,仅计算Y31元素所在行之右、与非零的交叉元素Y1j所在列相交的所有计算元素Y3j(1<j<n),这可大大减少零元素的计算;
考虑元素的对称稀疏性:仅仅计算Y31元素所在行的对角元素Y33及其以右、与非零的交叉元素Y1j所在列相交的所有计算元素Y3j(3≤j<n),这可再减少约50%非零元素的计算。因为对角元Y33以下的非零元素均可根据对称性从Y33以右的非零元素中获得。
3)将规格化前对角元以右的非零元素赋值给对角元以下的相应元素
在Y33的前一个元素Y32消元完成后,先将对角元Y33以右的非零元素Y3j的值按对称性赋给Y33以下的非零元素Yj3,然后再对Y33以右的非零元素Y3j规格化。
4)记录上三角矩阵中非零元素的位置以便回代过程中应用上三角元素的稀疏性。
步骤5:根据Y(n-1)′Zn=En (n-1)′以及上三角元素的稀疏性求Zn阵元素;
步骤6:根据对称性得其对角元Znn以左的所有元素;
步骤7:根据Y(n-1)′阵得Y(k-1)′阵,并根据Y(k-1)′Zk=Ek (k-1)′以及上三角元素的稀疏性等求第n-1~1列Zk阵的对角元Zkk及以上的非对角元素;
由于Ek阵元素结构的特点,在对第k行元素规格化后,Ek阵的对角元从ekk=1变化
为
步骤8:根据对称性得Zkk以左的所有元素;
本发明方法对Y(k-1)′Zk=Ek (k-1)′求解Zk阵元素的计算顺序为:Zn,…,Zk,…,Z1,且在计算各个Zk阵过程中,仅计算Zk阵对角元Zkk及其以上的非对角元素,即计算Zkk,Zk-1,k,Zk-2,k,…,Z1k,再根据对称性得Zkk以左的所有元素Zk,k-1,Zk,k-2,…,Zk1。这种计算方式具有较高的计算效率。
步骤9:将Z阵写入数据文件以备后续程序使用。
考虑到程序的结构化,所形成的Z阵数据文件可由下一个程序调用执行。
本发明方法具有以下几点优点。
(1)可根据元素本身的对称性和稀疏性进行计算,不但大大减少消元过程中对元素无效计算,而且可减少约50%非零元素的计算,从而大幅提高前代计算速度。
(2)消元过程中下三角的非零元素可通过与之对应的、规格化前的上三角的非零元素获得。
(3)回代过程利用单位矩阵E元素结构特点,可将Ek (k-1)′阵的求解转换成对其对角元素的计算,使得用对角元素和上三角元素就可完成相应计算,大大减少回代过程元素的计算量。
(4)回代过程利用E阵元素结构的特点以及上三角元素的稀疏性按对称方式回代求取Zk阵对角元Zkk及以上元素、根据对称性求Zkk以左元素,也可大幅提高回代计算速度。
本发明对传统的含规格化的高斯消元法在前代过程中利用了Y(n-1)′阵元素的对称稀疏性、E阵元素结构的特点,回代过程中改变了Zk阵的求解顺序和求解方式,利用了Y(n-1)′阵上三角元素的稀疏性等等,可大大提高计算速度。
附图说明
图1无规格化的高斯消元法对Y阵求Z阵的计算流程图。
图2LDU三角分解法对Y阵求Z阵的计算流程图。
图3本发明方法对Y阵求Z阵的计算流程图。
具体实施方式
本发明将通过以下实施例作进一步说明。
实施例1。
不考虑元素的稀疏性、考虑元素的稀疏性以及考虑元素的对称稀疏性计算过程的比较。进行含规格化的消元后,Y阵变成如下Y(n-1)′阵。
假设Y31≠0,则对Y31元素要进行消元。定义:Y11元素为对角元素;Y11右侧的所有元素Y1j均为交叉元素;Y31元素为消元元素;Y31右侧的所有元素Y3j为计算元素。
(1)不考虑元素的稀疏性,要计算Y31元素右侧所有的计算元素Y32、Y33、…。
(2)考虑元素的稀疏性,仅计算Y31元素所在行之右、与非零的交叉Y1j元素所在列交互的所有计算元素Y3j(1<j<n)。与(1)相比,大大减少了零元素计算。
(3)考虑元素的对称稀疏性,仅仅计算Y31元素所在行的对角元素Y33及其以右、与非零的交叉Y1j元素所在列交互的所有计算元素Y3j(3<j<n)。与(2)相比,再减少50%非零元素的计算。对Y33以左的元素无需进行消元计算,只是在对Y33左侧元素Y32消元完成后,需将Y33以右的非零元素按对应关系赋值给Y33以下的非零元素。本发明方法按(3)进行。
实施例2。
分别用传统的无规格化的高斯消元法(图1)、LDU三角分解法(图2)以及本发明方法(图3)对IEEE-30、-57、-118节点系统的Y阵求其Z阵元素,并比较其平均计算时间。计算结果如表1所示。
表1 无规格化的高斯消元法、LDU三角分解法与本发明求解Z阵计算时间的比较
T1:无规格化高斯消元法平均计算时间;
T2:LDU三角分解法平均计算时间;
T3:本发明方法平均计算时间;
T2/T1:LDU三角分解法与无规格化高斯消元法平均计算时间百分比;
T3/T1:本发明方法与无规格化高斯消元法平均计算时间百分比;
T3/T2:本发明方法与LDU三角分解法平均计算时间百分比。
根据表1可以看出:
LDU三角分解法的计算时间比无规格化的高斯消元法快30%~40%,但随着节点数越多,两者计算时间的差距在减小。
而本发明方法比无规格化的高斯消元法快80%~97%,比LDU三角分解法快65%~96%,且随着节点数越多,本发明方法计算时间的优势越明显。如对IEEE-118节点系统,与传统的高斯消元法和LDU三角分解法相比,本发明方法计算速度可提高约96~97%。
上述计算结果表明,本发明提出的基于对称稀疏矩阵技术的高斯消元法求取电力系统节点阻抗矩阵的方法与传统的无规格化的高斯消元法和LDU三角分解法相比,计算速度大大加快。
本方法可以采用任何一种编程语言和编程环境实现,这里采用C++编程语言,开发环境是Visual C++。
Claims (1)
1.一种基于对称稀疏矩阵技术的高斯消元法求取电力系统节点阻抗矩阵的方法,其特征包括以下步骤:
步骤1:读入n节点系统各线路支路数据;
步骤2:形成节点导纳矩阵Y;
步骤3:Y阵与E阵最后一列形成增广阵Bn=[Y En];
步骤4:根据对称稀疏性对Bn阵进行n-1次含规格化的高斯消元得Bn (n-1)′=[Y(n-1)′ En (n -1)′],并记录Y(n-1)′阵上三角中非零元素的位置;
步骤4具体实施过程如下:
(1)判断对角元以右非零的交叉元素用对称性确定对角元以下非零的消元元素;
(2)以按列消元方式,计算非零的交叉元素和非零的消元元素交互点上、并位于各行对角元以右的计算元素;
(3)将规格化前的对角元以右的非零元素赋值给对角元以下的相应元素;
(4)记录第一次消元过程中上三角非零元素的位置以便在后续的反复回代和前代过程中应用;
步骤5:根据Y(n-1)′Zn=En (n-1)′以及上三角元素的稀疏性求Zn阵元素;
步骤5具体实施过程如下:
(1)规定Zk阵的求取顺序按从第n~1列、Zk阵元素的求取顺序从Zkk~Z1k;
(2)对方程Y(n-1)′Zn=En (n-1)′中仅求取En阵中的对角元素
(3)根据步骤4记录的Y(n-1)′阵上三角元素的稀疏性,对方程Y(n-1)′Zn=En (n-1)′的求解,仅计算Zn阵对角元Znn及其以上的非对角元素;
步骤6:根据对称性得其对角元Znn以左的所有元素;
步骤7:根据Y(n-1)′阵得Y(k-1)′阵,并根据Y(k-1)′Zk=Ek (k-1)′以及上三角元素的稀疏性等求第n-1~1列Zk阵的对角元Zkk及以上的非对角元素;
(1)根据Y(n-1)′阵第k行及以上的元素得Y(k-1)′阵;
(2)对方程YZk=Ek中仅求取Ek阵中的对角元素
(3)根据步骤4记录的Y(n-1)′阵第k行及以上上三角元素的稀疏性,对方程Y(k-1)′Zk=Ek (k-1)′的求解,仅计算Zk阵对角元Zkk及其以上的非对角元素;
步骤8:根据对称性得Zkk以左的所有元素;
步骤9:将Z阵写入数据文件。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510026804.3A CN104714928B (zh) | 2015-01-20 | 2015-01-20 | 一种基于对称稀疏矩阵技术的高斯消元法求取电力系统节点阻抗矩阵的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510026804.3A CN104714928B (zh) | 2015-01-20 | 2015-01-20 | 一种基于对称稀疏矩阵技术的高斯消元法求取电力系统节点阻抗矩阵的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104714928A CN104714928A (zh) | 2015-06-17 |
CN104714928B true CN104714928B (zh) | 2018-10-19 |
Family
ID=53414275
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510026804.3A Expired - Fee Related CN104714928B (zh) | 2015-01-20 | 2015-01-20 | 一种基于对称稀疏矩阵技术的高斯消元法求取电力系统节点阻抗矩阵的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104714928B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105354422B (zh) * | 2015-11-12 | 2018-07-20 | 南昌大学 | 一种基于对称稀疏矩阵技术快速求取极坐标牛顿-拉夫逊法潮流的方法 |
CN105786769B (zh) * | 2016-02-15 | 2021-03-26 | 南昌大学 | 一种基于快速数据读取及对称稀疏因子表法在极坐标pq分解法潮流中的应用 |
CN105786984B (zh) * | 2016-02-15 | 2019-02-01 | 南昌大学 | 一种基于稀疏技术可快速读写的电力系统pq分解法潮流数据存贮方法 |
CN109284476B (zh) * | 2018-08-27 | 2023-05-02 | 南昌大学 | 非零元素随机存放和随机对称消元求取电力系统节点阻抗的方法 |
CN110717145A (zh) * | 2019-09-19 | 2020-01-21 | 南昌大学 | 基于对称稀疏矩阵技术的分段对称反向高斯-约当消元法 |
-
2015
- 2015-01-20 CN CN201510026804.3A patent/CN104714928B/zh not_active Expired - Fee Related
Non-Patent Citations (4)
Title |
---|
电力系统潮流分析软件开发;邱莉;《中国优秀硕士学位论文全文数据库-信息科学辑》;20091115;第2009年卷(第11期);第3章 * |
矩量法中阻抗矩阵的优化填充技术;韩国栋;《微波学报》;20100831;15-19 * |
配电网潮流计算的实用算法;李伟等;《南昌大学学报(工科版)》;20110930;第33卷(第3期);303-306 * |
配电网的回路阻抗法潮流计算与仿真可视化;王文举等;《同济大学学报(自然科学版)》;20120331;第40卷(第3期);459-467 * |
Also Published As
Publication number | Publication date |
---|---|
CN104714928A (zh) | 2015-06-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104714928B (zh) | 一种基于对称稀疏矩阵技术的高斯消元法求取电力系统节点阻抗矩阵的方法 | |
CN104317553B (zh) | 一种基于稀疏矩阵技术快速形成及读写电力系统节点导纳矩阵数据的方法 | |
CN104933528B (zh) | 一种基于稀疏矩阵技术快速形成电力系统潮流计算中雅可比矩阵的方法 | |
US11010130B2 (en) | Floating point processor prototype of multi-channel data | |
CN104317776A (zh) | 一种基于稀疏矩阵技术求解电力系统节点阻抗矩阵的方法 | |
WO2020103883A1 (zh) | 执行矩阵乘法运算的方法、电路及soc | |
CN105354422B (zh) | 一种基于对称稀疏矩阵技术快速求取极坐标牛顿-拉夫逊法潮流的方法 | |
CN104715422A (zh) | 一种基于对称稀疏矩阵技术的因子表法求取电力系统节点阻抗矩阵的方法 | |
CN110826186B (zh) | 基于对称稀疏矩阵技术和非零元素随机存放的lr三角分解法 | |
CN115146226B (zh) | 基于张量压缩方法的流数据处理方法、装置及设备 | |
CN105786769B (zh) | 一种基于快速数据读取及对称稀疏因子表法在极坐标pq分解法潮流中的应用 | |
CN104657337B (zh) | 一种基于cu三角分解求取电力系统节点阻抗矩阵的方法 | |
CN104391824A (zh) | 一种基于a=lr三角分解法快速求取电力系统节点阻抗矩阵的方法 | |
CN104391823A (zh) | 一种基于a=ldu三角分解法求解电力系统节点阻抗矩阵的方法 | |
CN104408026A (zh) | 一种基于对称稀疏矩阵技术的lr三角分解求取电力系统节点阻抗矩阵的方法 | |
CN114237548A (zh) | 基于非易失性存储器阵列的复数点乘运算的方法及系统 | |
CN109241492A (zh) | 快速求取电力系统节点阻抗矩阵的高斯-约当消元新算法 | |
CN105703359B (zh) | 一种对称稀疏因子表法在直角坐标pq分解法潮流计算中的应用 | |
CN110717145A (zh) | 基于对称稀疏矩阵技术的分段对称反向高斯-约当消元法 | |
CN109191016A (zh) | 快速求取电力系统节点阻抗矩阵的高斯-约当因子表法 | |
CN104572584B (zh) | 基于稀疏技术的lr分解求取电力系统节点阻抗矩阵方法 | |
Huerta et al. | Determining work partitioning on closely coupled heterogeneous computing systems using statistical design of experiments | |
CN109831235B (zh) | 一种信道矩阵的svd分解方法及装置 | |
CN113379046B (zh) | 卷积神经网络的加速计算方法、存储介质和计算机设备 | |
CN114268091B (zh) | 基于综合阻抗矩阵的故障处理方法、装置、设备及介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20181019 Termination date: 20220120 |
|
CF01 | Termination of patent right due to non-payment of annual fee |