CN107831823B - 一种用于分析和优化电网拓扑结构的高斯消元方法 - Google Patents
一种用于分析和优化电网拓扑结构的高斯消元方法 Download PDFInfo
- Publication number
- CN107831823B CN107831823B CN201711345248.1A CN201711345248A CN107831823B CN 107831823 B CN107831823 B CN 107831823B CN 201711345248 A CN201711345248 A CN 201711345248A CN 107831823 B CN107831823 B CN 107831823B
- Authority
- CN
- China
- Prior art keywords
- matrix
- data
- generalized
- bit
- result
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 52
- 230000008030 elimination Effects 0.000 title claims abstract description 44
- 238000003379 elimination reaction Methods 0.000 title claims abstract description 44
- 230000003287 optical effect Effects 0.000 claims abstract description 85
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 43
- 238000012545 processing Methods 0.000 claims abstract description 23
- 238000004458 analytical method Methods 0.000 claims abstract description 20
- 238000005516 engineering process Methods 0.000 claims abstract description 11
- 239000011159 matrix material Substances 0.000 claims description 158
- 238000006243 chemical reaction Methods 0.000 claims description 41
- 238000006467 substitution reaction Methods 0.000 claims description 22
- 238000004364 calculation method Methods 0.000 claims description 11
- 230000009466 transformation Effects 0.000 claims description 9
- 238000013461 design Methods 0.000 claims description 8
- 230000005540 biological transmission Effects 0.000 claims description 7
- 238000002360 preparation method Methods 0.000 claims description 6
- 230000008901 benefit Effects 0.000 abstract description 5
- 230000006872 improvement Effects 0.000 abstract description 3
- 238000005457 optimization Methods 0.000 abstract description 3
- 238000010937 topological data analysis Methods 0.000 abstract 1
- 239000004973 liquid crystal related substance Substances 0.000 description 8
- 230000008569 process Effects 0.000 description 6
- 238000010276 construction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06E—OPTICAL COMPUTING DEVICES; COMPUTING DEVICES USING OTHER RADIATIONS WITH SIMILAR PROPERTIES
- G06E3/00—Devices not provided for in group G06E1/00, e.g. for processing analogue or hybrid data
- G06E3/001—Analogue devices in which mathematical operations are carried out with the aid of optical or electro-optical elements
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06E—OPTICAL COMPUTING DEVICES; COMPUTING DEVICES USING OTHER RADIATIONS WITH SIMILAR PROPERTIES
- G06E3/00—Devices not provided for in group G06E1/00, e.g. for processing analogue or hybrid data
- G06E3/008—Matrix or vector computation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Mathematical Physics (AREA)
- Health & Medical Sciences (AREA)
- Nonlinear Science (AREA)
- Development Economics (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Water Supply & Treatment (AREA)
- Public Health (AREA)
- Computing Systems (AREA)
- Optics & Photonics (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Signal Processing For Digital Recording And Reproducing (AREA)
- Compression Or Coding Systems Of Tv Signals (AREA)
Abstract
本发明公开了一种用于分析和优化电网拓扑结构的高斯消元方法,该方法以提高电网拓扑的分析效率与优化电网结构为出发点,利用可重构光学处理器的资源优势,构造出大位数的广义加法器和广义乘法器,借助数据拼接和剪辑技术,将高斯消元法一次迭代内部多行、多列的运算作全并行处理,将传统算法中分解、前代、回代的三重循环改进为一重循环完成,以降低高斯消元法分析网络拓扑的延时,提高对复杂网络数据的分析和处理能力。本发明致力于提升高斯消元法动态实时处理复杂电网拓扑结构的能力,本发明能有效应对电网拓扑分析面临的海量数据问题,提高电网拓扑结构分析的效率。
Description
技术领域
本发明涉及一种用于分析和优化电网拓扑结构的高斯消元方法。
背景技术
现实中电力拓扑网络由场站设备和输电线路等组成,可以用无向图G、关联矩阵A等表示节点之间的逻辑关系,抽象出数学模型,以布尔代数的相关运算法则为基础,分析拓扑结构中任意两个结点之间的连通状态。2007年,山东大学的韩学山团队将电网拓扑分析转化为求A的传递闭包阵Ac,并且推导出如下结论:Ac=An-1,Ac=AAc+I,可见Ac的计算转化为类似线性方程组的求解,通过高斯消元法完成。传统的高斯消元法分析n结点的拓扑结构时,自乘运算量为n2(n-2)(2n-1),算法复杂度为O(2n4),韩学山团队将算法复杂度改进至O(n3)。尽管学者们致力于拓扑分析算法的改进研究,但是现有算法的局限性较大,加之电网拓扑结构复杂多变,拓扑分析的实时性、效率、计算耗时面临着严峻的挑战,设计一种更高效、计算延时更小的高斯消元法实现方案是迫在眉睫需要解决的问题。
2000年,上海大学的金翊教授及其科研团队提出三值光学计算机的概念和结构,采用无光态W和偏振方向相互正交的两个偏振光态H、V表达信息,以三值逻辑、MSD数字***为基础完成运算,其核心优势在于光学处理器是可重构的,即可以根据计算需求,随时构造出用户或者常用工程算法需要的二值或三值逻辑运算器,并且借助光学运算器数据位众多的特点,合理划分出多个运算分区,并行完成计算任务,达到提高算法效率、降低延时的目的。
发明内容
本发明的目的在于提出一种用于分析和优化电网拓扑结构的高斯消元方法,该方法以三值光学计算机核心计算技术为基础,以提高电网拓扑结构分析和优化的效率为目的。
本发明为了实现上述目的,采用如下技术方案:
一种用于分析和优化电网拓扑结构的高斯消元方法,以三值光学计算机的核心计算技术为基础,具体包括如下步骤:
s1.在三值光学计算机的用户界面上输入关联矩阵A、矩阵维数n以及运算法则,点击“确定”按钮,将运算请求送入到三值光学计算机内部;
s2.三值光学计算机生成SZG文件,该SZG文件中包含关联矩阵A的矩阵元素aij、矩阵维数n、广义乘法以及广义加法的运算法则;
s3.任务调度模块解析SZG文件获取n,用来申请和分配可重构光学处理器的数据位数;
s4.数据位管理模块计算需要的数据位数总量VT=2n(n-1),查询当前可重构光学处理器上空闲的数据位区段h~h+VT-1,并将数据位区段h~h+VT-1分配使用,任务调度模块整理关联矩阵A的矩阵元素aij,生成对应的操作数编码,并进行格式标准化处理;
s5.任务调度模块将分配的数据位区段h~h+VT-1、关联矩阵A的矩阵元素aij以及重构命令码发送至三值光学计算机的底层控制软件;
s6.重构模块根据广义乘法和广义加法的光路设计生成重构命令码,并执行这些重构命令码,构造出一个VM=n(n-1)位的广义乘法器M和一个VN=n(n-1)位的广义加法器N;
s7.底层控制软件实施高斯消元法的分解、前代以及回代的迭代计算步骤;
在每次运算时,将一次迭代内部的运算数据拼接并一次性送入可重构光学处理器进行变换,结束时通过剪辑技术生成每一个变换的结果;
回代运算输出的结果,记入矩阵Ac,就是电网拓扑的传递闭包阵;
s8.任务调度模块将运算输出矩阵Ac写入SZG结果文件中,经过文件解析步骤转化为十进制表示形式并返回该结果;
其中,任务调度模块、数据位管理模块和重构模块为三值光学计算机的软件模块。
优选地,所述步骤s7中,高斯消元法的分解、前代以及回代的迭代计算步骤具体为:
s71高斯消元法分解运算的实施步骤为:
s711.获取关联矩阵A的矩阵元素aij及矩阵维数n,完成迭代变量i、C11、O11、C12、O12、S1、R1的定义和赋值运算,令i的初始值为1;
s712.准备分解步骤的广义乘法运算的输入数据
将矩阵A第i行的矩阵元素ai,i+1...ain先拼接,复制n-i份,再拼接成一个(n-i)2的数据,末位补零成n(n-1)位的数据,送入到编码器生成广义乘法器的控制光路编码信息,记入变量C11;当i=1时,拼接复制补零后,得到数据a12...a1na12...a1n...a12...a1n0...0,将其送入编码器;
将矩阵A第i列的矩阵元素ai+1,i...ani先逐位复制n-i份,再拼接成一个(n-i)2的数据,末位补零成n(n-1)位的数据,送入到编码器生成广义乘法器的主光路编码信息,记入变量O11;当i=1时,将数据a21...a21a31...a31...an1...an10...0送入编码器;
s713.实施分解步骤的广义乘法运算
将C11和O11送入n(n-1)位广义乘法器中完成变换,解码器获取并保存结果,记入变量S1;
s714.准备分解步骤的广义加法运算的输入数据
将矩阵A元素ai+1,i+1~ann以行为单位进行拼接,得到一个(n-i)2位的数据,末位补零成n(n-1)位的数据,送入到编码器生成广义加法器的控制光路编码信息,记入变量C12;当i=1时,拼接补零后得到的数据为a22...a2na32...a3n...an2...ann;
将S1生成广义加法器的主光路编码信息,记入变量O12;
s715.实施分解步骤的广义加法运算
将C12和O12送入n(n-1)位广义加法器中进行变换,解码器获取变换结果,记入变量R1;
s716.解码器对结果R1进行按位剪辑,前(n-i)2位作为矩阵新元素ai+1,i+1~ann更新矩阵A;
s717.i增加1,解码器将上一轮迭代得到的矩阵进行回馈处理,重复分解步骤s712~s716,直到i=n-1,最终结果仍保存至矩阵A,A=LU,分解运算结束;
s72高斯消元法前代运算的实施步骤为:
s721.获取分解运算的结果矩阵A、L,完成迭代变量j、C21、O21、C22、O22、S2、R2的定义和赋值运算,矩阵B赋初值为单位矩阵,令j的初始值为1;
s722.准备前代步骤的广义乘法运算的输入数据
将矩阵B对角线以下的元素bj1...bjjbj1...bj,j+1...bj1...bj,n-1拼接成一个的数据,末位补零成n(n-1)位的数据,送入到编码器生成广义乘法器的控制光路编码信息,记入变量C21;当j=1时,拼接复制补零后,得到数据b11b11b12b11b12b13...b11b12...b1,n-10...0,将其送入编码器;
对应地,复制矩阵L的第j列元素,拼接成的数据lj+1,j...lj+1,jlj+2, j...lj+2,j...lnj...lnj,末位补零成n(n-1)位的数据,送入到编码器生成广义乘法器的主光路编码信息,记入变量O21;当j=1时,将数据l21l31l31l41l41l41...ln1...ln10...0送入编码器;
s723.实施前代步骤的广义乘法运算
将C21和O21送入n(n-1)位广义乘法器中完成变换,解码器获取并保存结果,记入变量S2;
s724.准备前代步骤的广义加法运算的输入数据
将矩阵B对角线以下的元素拼接成位数据bj+1,1...bj+1,jbj+2, 1...bj+2,j+1...bn1...bn,n-1,末位补零成n(n-1)位数据,送入到编码器生成广义加法器的控制光路编码信息,记入变量C22;当j=1时,拼接补零后的数据是b21b31b32b41b42b43...bn1bn2...bn,n-10...0;
将S2生成广义加法器的主光路编码信息,记入变量O22;
s725.实施前代步骤的广义加法运算
将C22和O22送入n(n-1)位广义加法器中进行变换,解码器获取变换结果,记入变量R2;
s727.j增加1,解码器对矩阵B元素进行回馈处理,重复前代步骤s722~s726,直到j=n-1,变换后的结果仍保存至矩阵B,前代运算结束;
s73高斯消元法回代运算的实施步骤为:
s731.获取分解、前代运算的结果矩阵U和B,完成迭代变量k、C31、O31、C32、O32、S3、R3的定义和赋值运算,矩阵Ac赋初值为0,令k的初始值为1;
s732.准备回代步骤的广义乘法运算的输入数据
将矩阵B第n-k+1行元素执行拼接操作,复制(n-k)份,再拼接成一个n(n-k)位数据,末位补零成n(n-1)位数据,送入到编码器生成广义乘法器的控制光路编码信息,记入变量C31;当k=1时,拼接补零后,得到数据bn1...bnnbn1...bnn...bn1...bnn0...0,将其送入编码器;
对应地,将矩阵U的第n-k+1列前n-k个元素复制n份拼接成n(n-k)位数据,末位补零成n(n-1)位的数据,送入到编码器生成广义乘法器的主光路编码信息,记入变量O31;当k=1时,将数据u1n...u1nu2n...u2n...un-1,n...un-1,n0...0送入编码器;
s733.实施回代步骤的广义乘法运算
将C31和O31送入n(n-1)位广义乘法器中完成变换,解码器获取并保存结果,记入变量S3;
s734.准备回代步骤的广义加法运算的输入数据
将矩阵B前n-k行元素顺序拼接成n(n-k)位数据b11...b1nb21...b2n...bn-k, 1...bn-k,n,末位补零成n(n-1)位数据,送入到编码器生成广义加法器的控制光路编码信息,记入变量C32,当k=1时,拼接补零后的数据是b11...b1nb21...b2n...bn-1,1...bn-1,n0...0;
将S3生成广义加法器的主光路编码信息,记入变量O32;
s735.实施回代步骤的广义加法运算
将C32和O32送入n(n-1)位广义加法器中进行变换,解码器获取变换结果,记入变量R3;
s736.解码器对结果R3进行按位剪辑,前n(n-k)位作为新元素更新矩阵B;
s737.k增加1,解码器对矩阵B元素进行回馈处理,重复回代步骤s732~s736,直到k=n-1,矩阵B变换后的结果保存至Ac,回代运算结束;
本发明具有如下优点:
本发明方法以提高电网拓扑的分析效率与结构优化为出发点,利用可重构光学处理器的资源优势,构造出大位数的广义加法器和广义乘法器,借助数据拼接和剪辑技术,将高斯消元法内部的一次迭代内部多行、多列的运算作全并行处理,将三重循环改进为一重循环实现,以降低网络拓扑分析延时、提高对复杂网络的数据分析和处理能力。
附图说明
图1为本发明中广义乘法器的光路设计结构图;
图2为本发明中广义加法器的光路设计结构图;
图3为本发明中高斯消元法分解运算的并行实现方案流程图。
具体实施方式
本发明的基本思想为:以提高电网拓扑的分析效率与结构优化为出发点,基于可重构光学处理器、数据拼接与剪辑等光学并行计算技术,设计高斯消元法的高效实现方案。
下面结合附图以及具体实施方式对本发明作进一步详细说明:
求解传递闭包阵Ac的过程分为分解、前代以及回代三个计算步骤。
为了实现本发明的目的,按照如下步骤优化技术方案:
1、分析迭代步骤的运算量和需要的运算器数量及规模
高斯消元法的内部步骤均为循环计算过程。
首次迭代需要的运算器数量最多,随着迭代控制变量的增大,矩阵中的数据处理规模降低,用到的运算器数量减少,故运算器的分析以首次迭代中用到的数量为准。
(1)分解运算
(2)前代运算
(3)回代运算
2、为高斯消元法设计可重构光学处理器的构造方案
将三个步骤所需要的运算器数量及规模列举表1中,经过分析可见,一次迭代过程中最多需要的运算器数量为n(n-1)个,每次运算处理对象都是1位,均为二值逻辑变换,并且变换过程中没有数据位数膨胀的问题。于是,为提高可重构光学处理器计算资源效率和改良算法耗时,每次运算时,尽可能将多个1位的数据拼接成一个大数据,一次性送入到可重构光学处理器发生变换,结束后再根据剪辑技术生成每一个变换的结果。
表1高斯消元法各步骤需要的处理器类型及规模
为高斯消元法设计的可重构光学处理器的构造策略如下:
将矩阵规模n作为申请可重构光学处理器数据位数的参数,分配一段连续的、长度为VT=2n(n-1)的空闲数据位数。
构造出一个VM=n(n-1)位的广义乘法器M和一个VN=n(n-1)位的广义加法器N。
3、设计运算器的光路结构,生成可重构光学处理器的重构命令码
根据广义加法和广义乘法的运算法则,写出对应的真值表,如表2和表3:
在三值光学计算机内部,用垂直偏振光V表示拓扑结构中节点的连通状态,用无光态W表示节点的断开状态,将真值表转化为用V和W表达的物理状态迁移表。
确定构造广义加法器和广义乘法器的光学基本处理单元并完成光路设计,如图1和图2所示,生成对应运算器的重构指令编码,构造出相应的运算器。
在图1和图2中,b1、b2、b3为控制光路输入,a1、a2、a3为主光路输入,v1、v2、v3、v4、v5、v6、v7为垂直偏振片,只能允许垂直偏振光穿过,h1为水平偏振片,只能允许水平偏振光穿过,d1、d2、d3为液晶的控制信号端,Lc1为常旋光液晶,Lc2、Lc3为常不旋光液晶,c1、c2为输出光信号。
广义乘法器的工作原理如下:主光路输入信号a1穿过垂直偏振片v2到达常旋光液晶Lc1,同时控制光路输入b1穿过垂直偏振片v1,当b1有信号时,该信号通过控制信号端d1作用到Lc1上,使Lc1不旋转主光路信号,当b1无信号时,d1不产生作用,Lc1将主光路信号的方向旋转90°,Lc1输出的光信号穿过垂直偏振片v3后,产生变换结果,输出光信号c1。
广义加法器的工作原理如下:主光路输入a2穿过垂直偏振片v4到达常不旋光液晶Lc2,同时控制光路输入b2穿过水平偏振片h1生成控制信号,该控制信号通过d2作用于常不旋光液晶Lc2,进而对穿过Lc2的信号产生旋转作用,之后该光信号穿过垂直偏振片v5,产生第一条光路的变换结果;主光路输入a3穿过垂直偏振片v6到达常不旋光液晶Lc3,同时控制光路输入b3生成常不旋光液晶Lc3的控制信号,该控制信号通过d3对穿过Lc3的输入信号产生旋转作用,之后该光信号穿过垂直偏振片v7产生第二条光路的变换结果;
第一条光路和第二条光路的变换结果叠加到一起,输出光信号c2。
本发明将用于电网拓扑结构中分析节点连通状态的高斯消元法,首次移植到三值光学计算平台上,设计出一种高斯消元法的高效方案,该技术方法的显著优势在于:
(1)首次根据电网拓扑节点状态分析的特点,分析高斯消元法的计算需求,设计出一种大位数的广义加法器和广义乘法器,借助光路结构实现广义加法和广义乘法运算。
(2)对传统高斯消元法的分解、前代以及回代计算步骤做出改进,设计了一种在可重构光学处理器上实现的创新方法。该方法将分解、前代和回代步骤的三重循环改进为一重循环实现,一次迭代中将矩阵多行、多列的布尔逻辑运算一并送入可重构光学处理器,执行一次逻辑变换完成,将耗时缩短至一个时钟周期。
(3)结合广义乘法、广义加法运算数量多、位数小、运算零散的特点,引入数据拼接与剪辑技术,使运算效率得到显著的提升。
(4)将分解、前代和回代运算总量从2n2(n-1)降低至3n(n-1),算法复杂度降低至O(n2)。
(5)本发明中的高斯消元法以关联矩阵A的维数n来申请可重构光学处理器的数据位数,当电网拓扑节点增多、关联矩阵规模增大时,可以分配更多的数据位数,构造规模更大的广义加法器和广义乘法器,可以有效提高高斯消元法对复杂电网拓扑结构的处理能力。
为了对本发明进行详细说明,首先对三值光学计算机的内部软硬件组成进行介绍:
三值光学计算机是一种新型计算机,与传统计算机相比,最显著的硬件优势是可重构光学处理器,它拥有数量众多的数据位数,可以根据运算的需求对数据位区段进行划分,构造出不同种类、不同规模的多个运算器,并行处理运算请求。另外,三值光学计算机内还包含着用户界面、任务调度模块、重构模块以及底层控制软件等一系列软件资源,不同软件模块执行不同的流程,软、硬件资源协同合作完成用户请求。其中,用户界面负责接收运算请求、数据和返回运算结果,SZG文件负责传输命令和数据,任务调度模块负责协调和处理运算,数据位管理模块负责管理和分配可重构光学处理器的数据位数,重构模块负责具体运算器的构造,底层控制软件负责完成迭代运算,编码器生成可重构光学处理器的输入数据编码,可重构光学处理器完成具体的数据变换,解码器负责处理可重构光学处理器输出的结果。
高斯消元法在内部实施时被划分为多个运算模块,由多个软、硬件模块相互协调配合完成。一种用于分析和优化电网拓扑结构的高斯消元方法包括如下步骤:
s1.在三值光学计算机的用户界面上输入关联矩阵A、矩阵维数n以及运算法则,点击“确定”按钮,将运算请求送入到三值光学计算机内部;
s2.三值光学计算机生成SZG文件,该SZG文件中包含关联矩阵A的矩阵元素aij、矩阵维数n、广义乘法以及广义加法的运算法则;
s3.任务调度模块解析SZG文件获取n,用来申请和分配可重构光学处理器的数据位数;
s4.数据位管理模块计算需要的数据位数总量VT=2n(n-1),查询当前可重构光学处理器上空闲的数据位区段h~h+VT-1,并将数据位区段h~h+VT-1分配使用,任务调度模块整理关联矩阵A的矩阵元素aij,生成对应的操作数编码,并进行格式标准化处理;
s5.任务调度模块将分配的数据位区段h~h+VT-1、关联矩阵A的矩阵元素aij以及重构命令码发送至三值光学计算机的底层控制软件;
s6.重构模块根据广义乘法和广义加法的光路设计生成重构命令码,并执行这些重构命令码,构造出一个VM=n(n-1)位的广义乘法器M和一个VN=n(n-1)位的广义加法器N;
s7.底层控制软件实施高斯消元法的分解、前代以及回代的迭代计算步骤;
在每次运算时,将一次迭代内部的运算数据拼接并一次性送入可重构光学处理器进行变换,结束时通过剪辑技术生成每一个变换的结果;
回代运算输出的结果,记入矩阵Ac,就是电网拓扑的传递闭包阵。
下面对高斯消元法的分解、前代以及回代步骤分别进行如下说明:
s71如图3所示,高斯消元法分解运算的实施步骤为:
s711.获取关联矩阵A的矩阵元素aij及矩阵维数n,完成迭代变量i、C11、O11、C12、O12、S1、R1的定义和赋值运算,令i的初始值为1。
s712.准备分解步骤的广义乘法运算的输入数据
将矩阵A第i行的矩阵元素ai,i+1...ain先拼接,复制n-i份,再拼接成一个(n-i)2的数据,末位补零成n(n-1)位的数据,送入到编码器生成广义乘法器的控制光路编码信息,记入变量C11;当i=1时,拼接复制补零后,得到数据a12...a1na12...a1n...a12...a1n0...0,将其送入编码器。
将矩阵A第i列的矩阵元素ai+1,i...ani先逐位复制n-i份,再拼接成一个(n-i)2的数据,末位补零成n(n-1)位的数据,送入到编码器生成广义乘法器的主光路编码信息,记入变量O11;当i=1时,将数据a21...a21a31...a31...an1...an10...0送入编码器。
s713.实施分解步骤的广义乘法运算
将C11和O11送入n(n-1)位广义乘法器中完成变换,解码器获取并保存结果,记入变量S1:
s714.准备分解步骤的广义加法运算的输入数据
将矩阵A元素ai+1,i+1~ann以行为单位进行拼接,得到一个(n-i)2位的数据,末位补零成n(n-1)位的数据,送入到编码器生成广义加法器的控制光路编码信息,记入变量C12;当i=1时,拼接补零后得到的数据为a22...a2na32...a3n...an2...ann。
将S1生成广义加法器的主光路编码信息,记入变量O12。
s715.实施分解步骤的广义加法运算
将C12和O12送入n(n-1)位广义加法器中进行变换,解码器获取变换结果,记入变量R1:
s716.解码器对结果R1进行按位剪辑,前(n-i)2位作为矩阵新元素ai+1,i+1~ann更新矩阵A:
s717.i增加1,解码器将上一轮迭代得到的矩阵进行回馈处理,重复分解步骤s712~s716,直到i=n-1,最终结果仍保存至矩阵A,A=LU,分解运算结束。
s72高斯消元法前代运算的实施步骤为:
s721.获取分解运算的结果矩阵A、L,完成迭代变量j、C21、O21、C22、O22、S2、R2的定义和赋值运算,矩阵B赋初值为单位矩阵,令j的初始值为1。
s722.准备前代步骤的广义乘法运算的输入数据
将矩阵B对角线以下的元素bj1...bjjbj1...bj,j+1...bj1...bj,n-1拼接成一个的数据,末位补零成n(n-1)位的数据,送入到编码器生成广义乘法器的控制光路编码信息,记入变量C21;当j=1时,拼接复制补零后,得到数据b11b11b12b11b12b13...b11b12...b1,n-10...0,将其送入编码器。
对应地,复制矩阵L的第j列元素,拼接成的数据lj+1,j...lj+1,jlj+2, j...lj+2,j...lnj...lnj,末位补零成n(n-1)位的数据,送入到编码器生成广义乘法器的主光路编码信息,记入变量O21;当j=1时,将数据l21l31l31l41l41l41...ln1...ln10...0送入编码器。
s723.实施前代步骤的广义乘法运算
将C21和O21送入n(n-1)位广义乘法器中完成变换,解码器获取并保存结果,记入变量S2:
s724.准备前代步骤的广义加法运算的输入数据
将矩阵B对角线以下的元素拼接成位数据bj+1,1...bj+1,jbj+2,1...bj+2,j+1...bn1...bn,n-1,末位补零成n(n-1)位数据,送入到编码器生成广义加法器的控制光路编码信息,记入变量C22;当j=1时,拼接补零后的数据是b21b31b32b41b42b43...bn1bn2...bn,n-10...0。
将S2生成广义加法器的主光路编码信息,记入变量O22。
s725.实施前代步骤的广义加法运算
将C22和O22送入n(n-1)位广义加法器中进行变换,解码器获取变换结果,记入变量R2:
s727.j增加1,解码器对矩阵B元素进行回馈处理,重复前代步骤s722~s726,直到j=n-1,变换后的结果仍保存至矩阵B,前代运算结束。
s73高斯消元法回代运算的实施步骤为:
s731.获取分解、前代运算的结果矩阵U和B,完成迭代变量k、C31、O31、C32、O32、S3、R3的定义和赋值运算,矩阵Ac赋初值为0,令k的初始值为1。
s732.准备回代步骤的广义乘法运算的输入数据
将矩阵B第n-k+1行元素执行拼接操作,复制(n-k)份,再拼接成一个n(n-k)位数据,末位补零成n(n-1)位数据,送入到编码器生成广义乘法器的控制光路编码信息,记入变量C31;当k=1时,拼接补零后,得到数据bn1...bnnbn1...bnn...bn1...bnn0...0,将其送入编码器。
对应地,将矩阵U的第n-k+1列前n-k个元素复制n份拼接成n(n-k)位数据,末位补零成n(n-1)位的数据,送入到编码器生成广义乘法器的主光路编码信息,记入变量O31;当k=1时,将数据u1n...u1nu2n...u2n...un-1,n...un-1,n0...0送入编码器。
s733.实施回代步骤的广义乘法运算
将C31和O31送入n(n-1)位广义乘法器中完成变换,解码器获取并保存结果,记入变量S3:
s734.准备回代步骤的广义加法运算的输入数据
将矩阵B前n-k行元素顺序拼接成n(n-k)位数据b11...b1nb21...b2n...bn-k, 1...bn-k,n,末位补零成n(n-1)位数据,送入到编码器生成广义加法器的控制光路编码信息,记入变量C32,当k=1时,拼接补零后的数据是b11...b1nb21...b2n...bn-1,1...bn-1,n0...0。
将S3生成广义加法器的主光路编码信息,记入变量O32。
s735.实施回代步骤的广义加法运算
将C32和O32送入n(n-1)位广义加法器中进行变换,解码器获取变换结果,记入变量R3:
s736.解码器对结果R3进行按位剪辑,前n(n-k)位作为新元素更新矩阵B:
s737.k增加1,解码器对矩阵B元素进行回馈处理,重复回代步骤s732~s736,直到k=n-1,矩阵B变换后的结果保存至Ac,回代运算结束。
s8.任务调度模块将运算输出矩阵Ac写入SZG结果文件中,经过文件解析步骤转化为十进制表示形式并返回该结果。
本发明能有效应对电网拓扑分析面临的海量数据问题,提高电网拓扑结构分析的效率。
当然,以上说明仅仅为本发明的较佳实施例,本发明并不限于列举上述实施例,应当说明的是,任何熟悉本领域的技术人员在本说明书的教导下,所做出的所有等同替代、明显变形形式,均落在本说明书的实质范围之内,理应受到本发明的保护。
Claims (2)
1.一种用于分析和优化电网拓扑结构的高斯消元法,其特征在于,该高斯消元法以三值光学计算机的核心计算技术为基础,具体包括如下步骤:
s1.在三值光学计算机的用户界面上输入关联矩阵A、矩阵维数n以及运算法则,点击“确定”按钮,将运算请求送入到三值光学计算机内部;
s2.三值光学计算机生成SZG文件,该SZG文件中包含关联矩阵A的矩阵元素aij、矩阵维数n、广义乘法以及广义加法的运算法则;
s3.任务调度模块解析SZG文件获取n,用来申请和分配可重构光学处理器的数据位数;
s4.数据位管理模块计算需要的数据位数总量VT=2n(n-1),查询当前可重构光学处理器上空闲的数据位区段h~h+VT-1,并将数据位区段h~h+VT-1分配使用,任务调度模块整理关联矩阵A的矩阵元素aij,生成对应的操作数编码,并进行格式标准化处理;
s5.任务调度模块将分配的数据位区段h~h+VT-1、关联矩阵A的矩阵元素aij以及重构命令码发送至三值光学计算机的底层控制软件;
s6.重构模块根据广义乘法和广义加法的光路设计生成重构命令码,并执行这些重构命令码,构造出一个VM=n(n-1)位的广义乘法器M和一个VN=n(n-1)位的广义加法器N;
s7.底层控制软件实施高斯消元法的分解、前代以及回代的迭代计算步骤;
在每次运算时,将一次迭代内部的运算数据拼接并一次性送入可重构光学处理器进行变换,结束时通过剪辑技术生成每一个变换的结果;
回代运算输出的结果,记入矩阵Ac,就是电网拓扑的传递闭包阵;
s8.任务调度模块将运算输出矩阵Ac写入SZG结果文件中,经过文件解析步骤转化为十进制表示形式并返回该结果;
其中,任务调度模块、数据位管理模块和重构模块为三值光学计算机的软件模块。
2.根据权利要求1所述的一种用于分析和优化电网拓扑结构的高斯消元法,其特征在于,所述步骤s7中,高斯消元法的分解、前代以及回代的迭代计算步骤具体为:
s71高斯消元法分解运算的实施步骤为:
s711.获取关联矩阵A的矩阵元素aij及矩阵维数n,完成迭代变量i、C11、O11、C12、O12、S1、R1的定义和赋值运算,令i的初始值为1;
s712.准备分解步骤的广义乘法运算的输入数据
将矩阵A第i行的矩阵元素ai,i+1...ain先拼接,复制n-i份,再拼接成一个(n-i)2的数据,末位补零成n(n-1)位的数据,送入到编码器生成广义乘法器的控制光路编码信息,记入变量C11;当i=1时,拼接复制补零后,得到数据a12...a1na12...a1n...a12...a1n0...0,将其送入编码器;
将矩阵A第i列的矩阵元素ai+1,i...ani先逐位复制n-i份,再拼接成一个(n-i)2的数据,末位补零成n(n-1)位的数据,送入到编码器生成广义乘法器的主光路编码信息,记入变量O11;当i=1时,将数据a21...a21a31...a31...an1...an10...0送入编码器;
s713.实施分解步骤的广义乘法运算
将C11和O11送入n(n-1)位广义乘法器中完成变换,解码器获取并保存结果,记入变量S1;
s714.准备分解步骤的广义加法运算的输入数据
将矩阵A元素ai+1,i+1~ann以行为单位进行拼接,得到一个(n-i)2位的数据,末位补零成n(n-1)位的数据,送入到编码器生成广义加法器的控制光路编码信息,记入变量C12;当i=1时,拼接补零后得到的数据为a22...a2na32...a3n...an2...ann;
将S1生成广义加法器的主光路编码信息,记入变量O12;
s715.实施分解步骤的广义加法运算
将C12和O12送入n(n-1)位广义加法器中进行变换,解码器获取变换结果,记入变量R1;
s716.解码器对结果R1进行按位剪辑,前(n-i)2位作为矩阵新元素ai+1,i+1~ann更新矩阵A;
s717.i增加1,解码器将上一轮迭代得到的矩阵进行回馈处理,重复分解步骤s712~s716,直到i=n-1,最终结果仍保存至矩阵A,A=LU,分解运算结束;
s72高斯消元法前代运算的实施步骤为:
s721.获取分解运算的结果矩阵A、L,完成迭代变量j、C21、O21、C22、O22、S2、R2的定义和赋值运算,矩阵B赋初值为单位矩阵,令j的初始值为1;
s722.准备前代步骤的广义乘法运算的输入数据
将矩阵B对角线以下的元素bj1...bjjbj1...bj,j+1...bj1...bj,n-1拼接成一个的数据,末位补零成n(n-1)位的数据,送入到编码器生成广义乘法器的控制光路编码信息,记入变量C21;当j=1时,拼接复制补零后,得到数据b11b11b12b11b12b13...b11b12...b1,n- 10...0,将其送入编码器;
对应地,复制矩阵L的第j列元素,拼接成的数据lj+1,j...lj+1,jlj+2, j...lj+2,j...lnj...lnj,末位补零成n(n-1)位的数据,送入到编码器生成广义乘法器的主光路编码信息,记入变量O21;当j=1时,将数据l21l31l31l41l41l41...ln1...ln10...0送入编码器;
s723.实施前代步骤的广义乘法运算
将C21和O21送入n(n-1)位广义乘法器中完成变换,解码器获取并保存结果,记入变量S2;
s724.准备前代步骤的广义加法运算的输入数据
将矩阵B对角线以下的元素拼接成位数据bj+1,1...bj+1,jbj+2,1...bj+2,j+ 1...bn1...bn,n-1,末位补零成n(n-1)位数据,送入到编码器生成广义加法器的控制光路编码信息,记入变量C22;当j=1时,拼接补零后的数据是b21b31b32b41b42b43...bn1bn2...bn,n- 10...0;
将S2生成广义加法器的主光路编码信息,记入变量O22;
s725.实施前代步骤的广义加法运算
将C22和O22送入n(n-1)位广义加法器中进行变换,解码器获取变换结果,记入变量R2;
s727.j增加1,解码器对矩阵B元素进行回馈处理,重复前代步骤s722~s726,直到j=n-1,变换后的结果仍保存至矩阵B,前代运算结束;
s73高斯消元法回代运算的实施步骤为:
s731.获取分解、前代运算的结果矩阵U和B,完成迭代变量k、C31、O31、C32、O32、S3、R3的定义和赋值运算,矩阵Ac赋初值为0,令k的初始值为1;
s732.准备回代步骤的广义乘法运算的输入数据
将矩阵B第n-k+1行元素执行拼接操作,复制n-k份,再拼接成一个n(n-k)位数据,末位补零成n(n-1)位数据,送入到编码器生成广义乘法器的控制光路编码信息,记入变量C31;当k=1时,拼接补零后,得到数据bn1...bnnbn1...bnn...bn1...bnn0...0,将其送入编码器;
对应地,将矩阵U的第n-k+1列前n-k个元素复制n份拼接成n(n-k)位数据,末位补零成n(n-1)位的数据,送入到编码器生成广义乘法器的主光路编码信息,记入变量O31;当k=1时,将数据u1n...u1nu2n...u2n...un-1,n...un-1,n0...0送入编码器;
s733.实施回代步骤的广义乘法运算
将C31和O31送入n(n-1)位广义乘法器中完成变换,解码器获取并保存结果,记入变量S3;
s734.准备回代步骤的广义加法运算的输入数据
将矩阵B前n-k行元素顺序拼接成n(n-k)位数据b11...b1nb21...b2n...bn-k,1...bn-k,n,末位补零成n(n-1)位数据,送入到编码器生成广义加法器的控制光路编码信息,记入变量C32,当k=1时,拼接补零后的数据是b11...b1nb21...b2n...bn-1,1...bn-1,n0...0;
将S3生成广义加法器的主光路编码信息,记入变量O32;
s735.实施回代步骤的广义加法运算
将C32和O32送入n(n-1)位广义加法器中进行变换,解码器获取变换结果,记入变量R3;
s736.解码器对结果R3进行按位剪辑,前n(n-k)位作为新元素更新矩阵B;
s737.k增加1,解码器对矩阵B元素进行回馈处理,重复回代步骤s732~s736,直到k=n-1,矩阵B变换后的结果保存至Ac,回代运算结束;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711345248.1A CN107831823B (zh) | 2017-12-15 | 2017-12-15 | 一种用于分析和优化电网拓扑结构的高斯消元方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711345248.1A CN107831823B (zh) | 2017-12-15 | 2017-12-15 | 一种用于分析和优化电网拓扑结构的高斯消元方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107831823A CN107831823A (zh) | 2018-03-23 |
CN107831823B true CN107831823B (zh) | 2021-03-12 |
Family
ID=61644348
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711345248.1A Active CN107831823B (zh) | 2017-12-15 | 2017-12-15 | 一种用于分析和优化电网拓扑结构的高斯消元方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107831823B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109063291B (zh) * | 2018-07-20 | 2021-07-13 | 西安交通大学 | 机电装备冷却通道结构智能化拓扑优化设计方法 |
CN110244817B (zh) * | 2019-05-18 | 2023-06-13 | 南京惟心光电***有限公司 | 一种基于光电计算阵列的偏微分方程求解器及其方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102012892B (zh) * | 2010-12-13 | 2012-05-02 | 上海大学 | 可重构的三值光学处理器 |
-
2017
- 2017-12-15 CN CN201711345248.1A patent/CN107831823B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN107831823A (zh) | 2018-03-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108009106A (zh) | 神经网络计算模组 | |
CN112162723A (zh) | 一种量子加法运算方法、装置、电子装置及存储介质 | |
CN110163354A (zh) | 一种计算装置及方法 | |
CN112114776A (zh) | 一种量子乘法运算方法、装置、电子装置及存储介质 | |
CN107341133B (zh) | 基于任意维数矩阵lu分解的可重构计算结构的调度方法 | |
CN105373517A (zh) | 基于Spark的分布式稠密矩阵求逆并行化运算方法 | |
CN109934336A (zh) | 基于最优结构搜索的神经网络动态加速平台设计方法及神经网络动态加速平台 | |
CN110276447A (zh) | 一种计算装置及方法 | |
CN107831823B (zh) | 一种用于分析和优化电网拓扑结构的高斯消元方法 | |
CN108491924B (zh) | 一种面向人工智能计算的神经网络数据串行流水处理装置 | |
CN110163350A (zh) | 一种计算装置及方法 | |
CN105227259A (zh) | 一种m序列并行产生方法和装置 | |
CN111461335A (zh) | 基于mpi多进程的含噪声单量子逻辑门实现方法及装置 | |
CN105322973A (zh) | 一种rs码编码器及编码方法 | |
CN112836813A (zh) | 一种用于混合精度神经网络计算的可重构脉动阵列*** | |
CN114511094B (zh) | 一种量子算法的优化方法、装置、存储介质与电子装置 | |
CN111461336B (zh) | 基于mpi多进程的含噪声双量子逻辑门实现方法及装置 | |
CN107368459B (zh) | 基于任意维数矩阵乘法的可重构计算结构的调度方法 | |
CN106873942B (zh) | 结构量计算机的msd乘法计算的方法 | |
Kai et al. | Design of MSD multiplier for ternary optical computer processor based on minimum module | |
Peng et al. | Realization of a tri-valued programmable cellular automata with ternary optical computer | |
CN209708122U (zh) | 一种计算单元、阵列、模块、硬件*** | |
CN112214200A (zh) | 一种量子减法运算方法、装置、电子装置及存储介质 | |
CN103902762A (zh) | 一种针对正定对称矩阵进行最小二乘方程求解的电路结构 | |
JP7227322B2 (ja) | 位相演算を行うための方法および装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |