一种基于稀疏矩阵技术快速形成及读写电力系统节点导纳矩阵数据的方法技术方案

一种基于稀疏矩阵技术快速形成及读写电力系统节点导纳矩阵数据的方法，包括以下步骤：构建节点导纳矩阵Y阵的数组Y(n,22)；形成并读取系统各线路支路数据；计算各节点的自导纳Yii、Yjj和上三角中i-j节点间非零元素的互导纳Yij，形成仅含对角元素和上三角非零元素的Y阵；将Y阵写入数据文件。本发明专利技术按规定方式形成、写入或读取、计算和存贮Y阵上三角中对角元素的行下标、非零的非对角元素的列下标以及相应的自导纳和互导纳值，不但大大减少了大量零元素和下三角元素的存贮空间，因而大大提高了Y阵的形成及Y阵数据的读写速度，而且存贮方式简单明了，非常易于后续程序对Y阵存贮数据的计算及应用。

【技术实现步骤摘要】
一种基于稀疏矩阵技术快速形成及读写电力系统节点导纳矩阵数据的方法
本专利技术属于电力系统分析计算领域，涉及一种快速形成及读写电力系统节点导纳矩阵数据的方法。
技术介绍
在大型电力系统节点导纳矩阵的形成、存贮及读入过程中，如不考虑节点导纳矩阵Y阵元素的稀疏性和对称性，会导致大量零元素和对称元素的存贮以及对称元素的计算，从而造成形成Y阵所需时间较长、所需存贮空间极大、读写相应Y阵的数据文件耗费时间较长。传统的不考虑元素稀疏性的Y阵数组形式为Y(n,2n)，可存贮Y阵全部元素的数值，元素的行、列下标直接由元素在Y阵中的位置确定。这种数组形式简单直观，方便对Y阵数据的处理，但由于大量零元素的存贮而占据极大的存贮空间。传统的考虑元素稀疏性Y阵的数组形式多样，可以省略大量零元素的存贮单元，但Y阵元素的存贮方式及存贮结构复杂，有时也未利用Y阵的对称性，尽管元素存贮效力较高，但对Y阵数据的处理并不方便。因此考虑元素稀疏性及对称性的存贮方式不但可大幅节省存贮单元，而且可减少形成Y阵所需时间、以及对Y阵数据文件的读写时间。目前，电力系统计算中稀疏矩阵技术运用很广，矩阵的存贮方案也很多，如按坐标存贮、按顺序存贮、按链表存贮等等。尽管这些存贮方式可以省去不少存贮单元，但还没有达到最优效果，而且这些存贮方式结构复杂，且对角元素与非对角元素分开存贮也使得存取过程繁琐，不利于对Y阵的数据处理。
技术实现思路
本专利技术的目的是克服现有技术的不足，提供一种基于稀疏矩阵技术快速形成及读写电力系统节点导纳矩阵数据的方法。本专利技术是通过以下技术方案实现的。本专利技术包括以下步骤：步骤1：构建节点导纳矩阵Y阵的数组Y(n,22)。本专利技术根据电力系统每个节点一般最多和6条支路相连的原则，同时考虑Y阵数组元素的对称性，定义Y阵数组为Y(n,22)，并按规定方式计算和存贮Y阵上三角中对角元素的行下标、非零的非对角元素的列下标、以及相应的自导纳和互导纳值、各节点连接的支路数三组数据，不但省去了大量非零元素和下三角元素的存贮，并且存贮方式简单直观，所建立的数据文件的读写速度极快，也可直接写入传统的不考虑元素稀疏性的Y阵数组进行后续消元等计算。定义Y阵数组为Y(n,22)，主要是考虑任一节点可能最多与6条支路相连，将其列数分为3组。第1组：“支路累加计数列S”，提供与相应节点连接的支路数，位于第1列，其数值Si由程序自动累加，以保证准确读取与相应节点连接的支路数，进一步节省存贮空间；第2组：“父节点组”，提供与相应节点对应的自导纳的行号i和参数，共3列，位于第2～4列；第3组：“子节点组”，提供与相应父节点连接的i<j的互导纳的列号j和参数，最多提供6条与“父节点”相连的“子节点”的列号和参数，共3*6＝18列，位于第5～22列，其中第5、8、11、14、17、20列存放与“父节点组”节点相连接的子节点的列号，其余存放相应的互导纳参数。“父节点”有几条支路与“子节点”相连则，则“子节点组”就有几组数据。如“子节点组”有3条支路，则“子节点组”就只有3*3＝9列数据。所以“子节点组”的数据一般未存满18列，而“支路累加计数列S”的作用就是帮助提高对数据的读取效率。步骤2：按规定方式形成并读取系统各线路支路数据；形成和读取支路数据I、J、R、X、K时要求按i<j以及j1<j2<j3<j4<j5<j6的方式，以便确保仅形成Y阵的对角元素和上三角元素。步骤3：计算各节点的自导纳Yii、Yjj和上三角中i-j对应节点间的非零元素的互导纳Yij，形成仅含对角元素和上三角非零元素的Y阵；步骤3中具体实施过程如下：(1)Y阵中第1列、“支路累加计数列S”的S值帮助确定第i行对应的节点有几条支路相连，以便准确读写与第i行对应的节点有关的互导纳Yij，但不包括Yji。(2)根据1条i-j对应节点间的支路数据求取i、j对应节点的自导纳Yii、Yjj，对角元的行下标i、j作为“父节点”始终存放在Y阵第i、j行的第2列，相应的电导gii、gjj及电纳bii、bjj存放在第3～4列；如果在第i、j行的第2列上已经存有i、j值，则相应的电导gii、gjj及电纳bii、bjj要与原来存放在第3～4列的电导g′ii、g′jj及电纳b′ii、b′jj进行累加。(3)根据该支路数据求取i-j对应节点间的互导纳Yij，第1次出现的j1对应节点下标存放在Y阵第i行的第5列，相应的电导gi,j1、电纳bi,j1存放在第6～7列；第2次出现的j2对应节点下标存放在第8列，相应的电导gi,j2、电纳bi,j2存放在第9～10列；依此类推。(4)计算完成Y阵中所有对角元素和上三角的非零元素，得到仅含对角元素和上三角非零元素的Y阵。表1为本专利技术考虑元素稀疏性时Y阵元素的形成及存贮形式。表1本专利技术考虑元素稀疏性时Y阵元素的形成及存贮形式步骤4：将Y阵写入数据文件。考虑到程序的结构化，形成Y阵程序到此结束，而所形成的Y阵数据文件的调用则由下一个程序执行。在新的计算程序中打开按上述方式形成并存储的Y阵数据文件，可将Y阵数据文件数据根据需要直接读入Y(n,2n)或Y(n,22)数组，以便进行后续计算。在新的计算程序中读入Y(n,22)数据文件数据比读入Y(n,2n)数据文件数据所需时间要少得多，且可直接将Y(n,22)数据文件的数据写入Y(n,22)数组或Y(n,2n)数组以便后续程序计算。本专利技术按规定方式形成、写入或读取、计算和存贮Y阵上三角中对角元的行下标、非零的非对角元的列下标以及相应的自导纳和互导纳值以及1列“支路累加计数列S”的数据。这种存贮方式不但大大减少了Y阵的存贮空间，并且简单明了，易于数据的检索、修改和应用，且大大降低了Y阵的形成及数据文件的读取时间。本专利技术与传统的不考虑元素稀疏性的Y阵元素存贮形式相比，不但大大减少了大量零元素和下三角元素的存贮空间，因而大大提高了Y阵的形成及Y阵数据的读写速度。本专利技术与各种传统的、考虑元素稀疏性的Y阵元素存贮方案，如按坐标存贮、按顺序存贮、按链表存贮等等方案相比，还可进一步减少存贮空间、提高存贮Y阵数据的读写速度，而且存贮方式简单明了，非常易于后续程序对Y阵存贮数据的计算及应用。附图说明图1为本专利技术的流程图。图2为不考虑元素稀疏性及对称性时形成Y阵的流程图。图3为本专利技术考虑元素稀疏性及对称性时形成Y阵的流程图。具体实施方式本专利技术将通过以下实施举例作进一步说明。实施例1。以IEEE-300节点系统为例，分别比较传统的不考虑元素稀疏性时Y阵元素的存贮形式和本专利技术Y阵元素的存贮形式。(1)传统的不考虑元素稀疏性时Y阵元素的存贮形式传统的不考虑元素稀疏性时Y阵元素的存贮包括了所有的非零元素和下三角元素，其Y阵元素的存贮形式如表2。表2不考虑元素稀疏性时IEEE-300节点系统Y阵元素的存贮形式根据表2可以看出，不考虑稀疏性时IEEE-300节点系统Y阵实部、虚部元素的存贮需300*600个单元。(2)本专利技术对IEEE-300节点系统Y阵元素的存贮形式根据每个节点一般最多和6条支路相连的原则，且仅存贮Y阵的对角元素和上三角的非零元素，此时Y阵元素的存贮形式如表3。表3本专利技术对IEEE-300节点系统Y阵元素的存贮形式根据表3可以看本文档来自技高网...

【技术保护点】
一种基于稀疏矩阵技术快速形成及读写电力系统节点导纳矩阵数据的方法，其特征是包括以下步骤：步骤1：构建节点导纳矩阵Y阵的数组Y(n,22)；定义Y阵数组为Y(n,22)，将其列数分为3组，计算和存贮Y阵上三角中对角元素的行下标、非零的非对角元素的列下标、以及相应的自导纳和互导纳值、各节点连接的支路数三组数据；第1组：“支路累加计数列S”，提供与相应节点连接的支路数，位于第1列，其数值Si由程序自动累加；第2组：“父节点组”，提供自导纳的行号i和参数，共3列，位于第2～4列；第3组：“子节点组”，提供i<j的互导纳的列号j和参数，最多提供6条与“父节点”相连的“子节点”的列号和参数，共3*6＝18列，位于第5～22列，其中第5、8、11、14、17、20列存放与“父节点组”节点相连接的子节点的列号，其余存放相应的互导纳参数；步骤2：形成并读取系统各线路支路数据；形成和读取支路数据(I、J、R、X、K)时按节点号i<j以及j1<j2<j3<j4<j5<j6的方式，以确保仅形成Y阵的对角元素和上三角元素；步骤3：计算各节点的自导纳Yii、Yjj和上三角中i‑j节点间非零元素的互导纳Yij，形成仅含对角元素和上三角非零元素的Y阵；具体过程如下：(1)Y阵中第1列、“支路累加计数列S”的Si值帮助确定第i行(或第i节点)有几条支路相连，以便准确读写与i节点有关的互导纳Yij；(2)根据1条i‑j节点间的支路数据求取i、j节点的自导纳Yii、Yjj，对角元的行下标i、j作为“父节点”始终存放在Y阵第i、j行的第2列，相应的电导gii、gjj及电纳bii、bjj存放在第3～4列；如果在第i、j行的第2列上已经存有i、j值，则相应的电导gii、gjj及电纳bii、bjj要与原来存放在第3～4列的电导g′ii、g′jj及电纳b′ii、b′jj进行累加；(3)根据该支路数据求取i‑j节点间的互导纳Yij，第1次出现的j1节点下标存放在Y阵第i行的第5列，相应的电导gi,j1、电纳bi,j1存放在第6～7列；第2次出现的j2节点下标存放在第8列，相应的电导gi,j2、电纳bi,j2存放在第9～10列；依此类推；(4)计算完成Y阵中所有对角元素和上三角的非零元素，得到仅含对角元素和上三角非零元素的Y阵；步骤4：将Y阵写入数据文件。...

【技术特征摘要】
1.一种基于稀疏矩阵技术快速形成及读写电力系统节点导纳矩阵数据的方法，其特征是包括以下步骤：步骤1：构建节点导纳矩阵Y阵的数组Y(n,22)；定义Y阵数组为Y(n,22)，将其列数分为3组，计算和存贮Y阵上三角中对角元素的行下标、非零的非对角元素的列下标、以及相应的自导纳和互导纳值、各节点连接的支路数三组数据；第1组：“支路累加计数列S”，提供与相应节点连接的支路数，位于第1列，其数值Si由程序自动累加；第2组：“父节点组”，提供与相应节点对应的自导纳的行号i和参数，共3列，位于第2～4列；第3组：“子节点组”，提供与相应父节点连接的i<j的互导纳的列号j和参数，最多提供6条与“父节点”相连的“子节点”的列号和参数，共3*6＝18列，位于第5～22列，其中第5、8、11、14、17、20列存放与“父节点组”节点相连接的子节点的列号，其余存放相应的互导纳参数；步骤2：形成并读取系统各线路支路数据；形成和读取支路数据I、J、R、X、K时按i<j以及j1<j2<j3<j4<j5<j6的方式，以确保仅形成Y阵的对角元素和上三角元素；步骤3：计算各节点的自导纳Yii、Yj...

【专利技术属性】
技术研发人员：陈恳，万新儒，汪亚茜，彭丽君，
申请(专利权)人：南昌大学，
类型：发明
国别省市：江西;36