一种快速存贮及读取电力系统稀疏矩阵数据的方法，属于电力系统分析计算领域。本发明专利技术对导纳矩阵Y仅存贮主节点的行号、与主节点连接的子节点数S、主节点和子节点非零元素的列号和参数，大大减少了存贮单元和数据文件的读取时间，对Y阵数据的检索、修改、应用及对形成雅克比矩阵、计算节点电流、功率等较为简单直观、方便快捷。本发明专利技术对IEEE-57～-300节点系统进行验算完全可行，对IEEE-300节点系统与传统方法中不考虑稀疏性时对Y阵元素的存贮方式相比，存贮单元仅为传统方法的4.33％，写入或读取数据文件的时间分别为传统方法的6.55％和6.39％，且节点数越多，优势越明显。本发明专利技术可用于各工程领域中对称或不对称稀疏矩阵数据的快速存贮及读取。
【技术实现步骤摘要】
本专利技术属于电力系统分析计算领域。
技术介绍
电力系统分析计算中涉及到许多对称或不对称的极度稀疏矩阵。这些稀疏矩阵数 据在存贮、读取及计算过程中，若不考虑稀疏性，会导致大量零元素的存贮、读取及计算，从 而造成存贮空间的极大浪费、读写数据文件及计算过程时间较长。因此，考虑对称或不对称 稀疏矩阵元素的存贮方式不但可大幅节省存贮单元，而且可大大减少数据文件的存贮、读 取及计算过程的时间。 传统的不考虑元素稀疏性的Y阵数组形式为Y (η，2η)，存贮Y阵全部元素的数值， 元素的行、列下标直接由元素在Y阵中的位置确定。这种数组形式简单直观，方便对Y阵数 据的处理，但大量零元素的存贮导致存贮空间的极大浪费。且在对如节点电流等连续求和 计算过程中，如不引入非零判断语句，其计算效率较低；如引入非零判断语句，则计算速度 将受到影响。 电力系统计算中稀疏矩阵技术运用很广，矩阵的存贮方案也很多，如按坐标存贮、 按顺序存贮、按链表存贮等等。尽管这些存贮方式虽可节省不少存贮单元，但由于存贮方式 和结构复杂，对角元与非对角元分开存贮使得数据的读取过程繁琐，不利于数据的计算处 理，计算过程达不到最优效果。 虽然对称的Y阵可以按一定方式仅存贮其上三角的非零元素，根据对称性获得下 三角的非零元素，但这种计算方式在程序计算中，角标的变换、元素之间的赋值其实也占用 很多时间，而且有时也并不方便。
技术实现思路
为了克服上述现有技术的不足，本专利技术提出一种快速存贮及读取电力系统稀疏矩 阵数据的方法。 本专利技术将Y阵中的节点分为主节点和子节点。仅存贮主节点的行号、与主节点连 接的子节点数S(包括主节点）、主节点以及具有非零元素子节点的列号和参数（自导纳和 互导纳）。这种存贮方式不但大大减少了 Y阵的存贮空间，且简单明了，易于数据的检索、修 改和计算应用，还大大减少了 Y阵数据的存贮单元，以及Y阵数据文件的存贮和读取时间。 本专利技术是通过以下技术方案实现的。 本专利技术所述的包括以下步 骤： 步骤1 :定义Y阵数组Y (η，26)。 本专利技术先确定相应系统中连接子节点数最多的主节点，假设为7个，加主节点共8 个，即S mx= 8。因此，一行中主节点和子节点的列号及参数（实部、虚部）需24个单元，加 行号和节点数二个存贮单元，共需26个存贮单元。因此，定义Y阵数组为Y (η，26)，将其每 行分为3组，分别存贮该行的行号、与行号对应的节点数、主节点和所有非零元素子节点的 列号及参数，具体方式如下。 第1组："主节点行号i"，位于第1列，为主节点对应的行号，为检验数据之用； 第2组："节点数S1 "，位于第2列，为主节点及与该主节点连接的子节点数之和，S 值由程序自动累加以保证快速读取对应主节点和子节点的参数； 第3组："主节点、子节点列号j及参数"，位于第3~26列；按递增顺序存贮主节 点及与该主节点连接的子节点的列号j，分别位于第3*k列；相应主节点或子节点的自导纳 和互导纳（均分实部、虚部）分别位于3*k+l、3*k+2列，k的取值为：1~S 1; 本专利技术方法Y阵元素的存贮方式如下表所示。 步骤2 :从数据文件中读取全部支路数据； 步骤3 :计算各节点的自导纳和互导纳； 步骤4 :形成Y (n，26)形式的Y阵； 步骤5 :将Y (n，26)写入数据文件以备后续程序使用。 考虑到程序的模块化，形成Y阵程序到此结束，而Y阵数据文件的调用则由下一个 程序执行。 在新的计算程序中打开按上述方式存贮的Y(n，26)数据文件，将数据文件数据根 据需要直接读入Υ(η，2η)数组（如用于消元计算等）或Υ(η，26)数组（如形成J阵或计算 功率等），以便后续应用。读入Y (η，26)数据文件的数据比读入Y (η，2η)数据文件的数据所 需时间要少得多。 本专利技术方法可用于对称或不对称稀疏矩阵数据的存贮。与传统的不考虑元素稀疏 性的Y阵元素存贮形式相比，不但大大减少了大量零元素存贮空间，且大大提高了 Y阵数据 的读写速度。与按坐标存贮、按顺序存贮、按链表存贮等等方案相比，仍较多地减少了存贮 空间，提高了数据的读写速度，且后续数据的计算处理等更为方便。与仅存贮上三角非零元 素的方式相比，几乎不需增加任何存贮单元，且在形成J阵、计算IpPpQJt，同时使用上下 三角的非零元素，比根据上三角的非零元素按对称性获得下三角的非零元素更方便、也更 快捷。而且本专利技术方法存贮方式简单直观，数据文件的存贮及读取速度快，可直接用于相应 的后续计算。【附图说明】 图1为本专利技术方法的流程图。【具体实施方式】 本专利技术将通过以下实施例作进一步说明。 本专利技术方法对于对称及不对称稀疏矩阵都具有同样的计算效率。为举例方便，以 对称的节点导纳矩阵Y为例。 实施例1。对IEEE-300节点系统分别比较传统的不考虑元素稀疏性时Y阵元素的 存贮方式和本专利技术方法对Y阵元素的存贮方式。 (1)传统的不考虑元素稀疏性时Y阵元素的存贮方式如表1。 表1不考虑元素稀疏性时IEEE-300节点系统Y阵元素的存贮方式 根据表1可以看出，不考虑稀疏性时IEEE-300节点系统Y阵元素的存贮单元数为 300*600 个。 (2)本专利技术方法对IEEE-300节点系统Y阵元素的存贮形式 对IEEE-300节点系统，主节点连接子节点数最多为7个，加主节点共8个，因此定 义Y阵数组为Y (n，26)，本专利技术方法对Y阵元素的存贮方式如表2。 表2本专利技术方法对IEEE-300节点系统Y阵元素的存贮方式 根据表2可以看出，本专利技术方法对Y阵的存贮单元数为300*26个，仅为不考虑元 素稀疏性时Y阵元素存贮单元的4. 3 %。且节点数越大，节省的单元数越多。 由于本专利技术只存贮Y阵上下三角的非零元素，大大减少了 Y阵零元素的存贮单元， 并大大加快了 Y阵数据文件的存贮和读取速度，且用Y (n，26)数组比Y (η，2η)在形成J阵、 计算L、Pp QJt计算效率更高。本专利技术方法与仅存贮Y阵上三角非零元素的方法相比，存 IC单元也并未增加，且在形成J阵、计算Ii、Pi、Qi时，由于无需考虑对称性进行赋值，计算速 度反而略有增加。 实施例2。分别用传统的不考虑稀疏性时Y阵的存贮方法以及本专利技术方法对 IEEE-57、-118、-300节点系统的Y阵进行存贮，比较其数据文件读写过程的平均迭代时间 及所需的存贮单元数。 表3给出了用传统方法与本专利技术方法对上述系统Y阵读写数据过程的平均时间及 存贮单元数的比较。 表3传统方法与本方法对IEEE系统Y阵数据文件读写时间及存贮单元数的比较 V Tn、U1:传统方法对Y阵数据文件的平均读写时间、所需的存贮单元数； T2、T21、U2:本专利技术方法对Y阵数据文件的平均读写时间、所需的存贮单元数； W T21/Tn、U2AJ1:本专利技术方法与传统方法的平均读写时间、所需存贮单元数的 百分比。 可以看出： 1.无论读或写数据文件，本专利技术方法的存贮单元数和速度大大优于传统方法。 以IEEE-300节点系统为例，本专利技术方法写数据文件的时间仅为传统方法的 6. 55% ;本专利技术方法读数据文件的时间仅为传统方法的6. 39% ;本专利技术方法所需存贮单元 仅为传统方法的4. 33%。 2.电力系统节点数越多，本本文档来自技高网...
【技术保护点】
一种快速存贮及读取电力系统稀疏矩阵数据的方法，其特征包括以下步骤：步骤1：定义Y阵数组Y(n,26)；定义Y阵数组为Y(n,26)，将其每行分为3组，分别存贮该行的行号、与行号对应的节点数、主节点和所有非零元素子节点的列号及参数，具体方式如下：第1组：“主节点行号i”，位于第1列，为主节点对应的行号，为检验数据之用；第2组：“节点数Si”，位于第2列，为主节点及与该主节点连接的子节点数之和，S值由程序自动累加以保证快速读取对应主节点和子节点的参数；第3组：“主节点、子节点列号j及参数”，位于第3～26列；按递增顺序存贮主节点及与该主节点连接的子节点的列号j，分别位于第3*k列；相应主节点或子节点的自导纳和互导纳分别位于3*k+1、3*k+2列，k的取值为：1～Si；步骤2：从数据文件中读取全部支路数据；步骤3：计算各节点的自导纳和互导纳；步骤4：形成Y(n,26)形式的Y阵；步骤5：将Y(n,26)写入数据文件以备后续程序使用。
【技术特征摘要】
【专利技术属性】
技术研发人员：陈恳，万新儒，邵尉哲，席小青，
申请(专利权)人：南昌大学，
类型：发明
国别省市：江西;36
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。