一种用于电力系统线性化模型降阶的RNN矢量化方法

一种用于电力系统线性化模型降阶的rnn矢量化方法
技术领域
1.本发明涉及电力系统降阶领域，具体涉及一种用于电力系统线性化模型降阶的rnn矢量化方法。

背景技术：

2.电力系统规模庞大，对其直接进行研究与分析，既是不方便也是不必要的，因此，电力系统模型降阶一直以来都颇受关注。目前，线性系统的模型降阶方法已较为成熟，其中，平衡截断法能在很大程度上保证降阶系统模型的输入输出特性与原系统模型的输入输出特性一致，且降阶系统模型能保持原系统稳定性，在频域上也具有全局的降阶误差，这些优势使得平衡截断的模型降阶方法受到广泛的关注(朱泽翔.电力系统动态参数辨识与模型降阶研究[d].浙江大学,2018.)。然而，平衡截断法的计算速度受制于lyapunov方程的求解效率，这也一定程度上限制了其应用范围。
[0003]
关于lyapunov方程的高效求解方法，学界已有丰富的研究成果，譬如经典的bartels-stewart方法(bartels r.,stewart g.(1972).solution of the equation a x+x b＝cm acm,15(9),820
–
826.doi:10.1145/361573.361582)、最小残差方法(lin y.,simoncini v.(2013)minimal residual methods for large scale lyapunov equations.appl.numer.math.,72,52
–
71.doi:10.1016/j.apnum.2013.04.004)、低秩迭代法(stykel t.(2008).low-rank iterative methods for projected generalized lyapunov equations.electron.trans.numer.anal.,30(1),187
–
202.doi:10.1080/14689360802423530)等。这些研究都是针对软件算法，而近年来随着循环神经网络(rnn)的发展，基于rnn的lyapunov方程在线求解器也引发关注(yi c.,chen y.,lu z.(2011).improved gradient-based neural networks for online solution of lyapunov matrix equation.inform.process.lett.,111(16),780
–
786.doi:10.1016/j.ipl.2011.05.010，yi c.,chen y.,lan x.(2013)parison on neural solvers for the lyapunov matrix equation with stationary&nonstationary coeffificients.appl.math.model.,37(4),2495
–
2502.doi:10.1016/j.apm.2012.06.022，xiao l.,liao b.(2016).a convergence-accelerated zhang neural network and its solution application to lyapunov equation,neurocomputing,193,213
–
218.doi:10.1016/j.neucom.2016.02.021)。与传统的软件算法不同，基于rnn的求解方法，充分利用rnn的并行处理能力和易于硬件实现的特点，进一步提高了lyapunov方程的求解效率。
[0004]
rnn被用来求解lyapunov方程，其最终目的是构建出lyapunov方程的在线求解器，因此，对在线计算而言，提高计算效率是最关键的任务。不少学者通过提高rnn模型的收敛性来提高计算效率(xiao l.,zhang y.,hu z.,dai j.(2019).performance benefits of robust nonlinear zeroing neural network for finding accurate solution of lyapunov equation in presence of various noises.ieee trans.ind.inf.,15(9),
5161-5171.doi:10.1109/tii.2019.2900659)，而事实上，rnn模型的矢量化运算效率也还有提升的空间。另外，rnn模型的硬件电路接线复杂、元件数量多，不利于减小硬件体积及降低故障率。

技术实现要素：

[0005]
针对求解lyapunov方程的实时性需求及rnn模型的硬件电路接线复杂、元件数量多等问题，本发明以rnn的矢量化运算为研究对象，发现矢量化后的系数矩阵具有固定的元素构成，因此提出一种用于电力系统线性化模型降阶的rnn矢量化方法，比传统的基于克罗内克积的矢量化方法更简单且计算效率更高，在软件耗时上有明显优势，还可降低硬件电路接线复杂度及故障率，有利于减小硬件体积。
[0006]
本发明的目的至少通过如下技术方案之一是实现。
[0007]
一种用于电力系统线性化模型降阶的rnn矢量化方法，包括以下步骤：
[0008]
s1、基于电力系统线性化模型得到lyapunov方程并传送参数矩阵到矢量化功能模块；
[0009]
s2、根据矢量化矩阵的结构特点，完成矢量化功能模块的搭建，并将输出信号传送到循环神经网络(rnn)功能模块；
[0010]
s3、搭建循环神经网络(rnn)功能模块；
[0011]
s4、确定计算参数并驱动循环神经网络(rnn)功能模块，从而得到lyapunov方程的解信号；
[0012]
s5、基于步骤s4的解信号，得到电力系统的线性化降阶模型。
[0013]
进一步地，步骤s1具体如下：
[0014]
电力系统的小信号模型为：
[0015][0016]
其中，a、b、c、d均为电力系统线性化模型的常数矩阵，δxo、分别表示电力系统线性化模型内部的状态变量及其微分量，δu表示电力系统线性化模型的输入变量，δyo表示电力系统线性化模型的输出变量，n表示状态变量及其微分量的维数，m、p分别表示输入变量和输出变量的维数，表示实数域；δxo一般包括线性化的发电机转子角、发电机角频率及动态元件的状态量，δu和δyo需要根据具体研究对象来确定。
[0017]
基于平衡截断法原理得到电力系统的对偶lyapunov方程：
[0018][0019]
其中，g、x分别为可控性矩阵和可观性矩阵；
[0020]
对于可观性矩阵，lyapunov方程如下：
[0021]at
x(t)+x(t)a＝-q
ꢀꢀꢀ
(3)
[0022]
其中，q是常系数对称正定矩阵且q＝c
t
c，x(t)是可观性lyapunov方
程的解，t为时间变量；
[0023]
需要说明的是，公式(3)与公式(2)的可观性矩阵在形式上有不同，由于本发明采用基于循环神经网络的迭代方法求解可观性矩阵，所以公式(3)的解矩阵x(t)含时间变量。
[0024]
对于可控性矩阵，lyapunov方程如下：
[0025]
ag(t)+g(t)a
t
＝-w
ꢀꢀꢀ
(4)
[0026]
其中，w是常系数对称正定矩阵且w＝bb
t
，g(t)是可控性lyapunov方程的解，
[0027]
将参数矩阵a、q和w传送到矢量化功能模块。
[0028]
进一步地，步骤s2中，矢量化功能模块针对的是参数矩阵矢量化，其核心是计算为方便描述，令为方便描述，令表示克罗内克积运算，i表示n阶单位矩阵；经研究，的固有元素构成为：
[0029][0030]
其中，a
ij
表示矩阵a第i行第j列的元素，1≤i≤n，1≤j≤n。
[0031]
因此，直接根据的结构特点求取
[0032]
需要指出，矢量化功能模块的具体形式可以是软件算法，也可以是硬件电路，视需求而定。
[0033]
若采用软件算法搭建首先根据分别根据和的元素构成将矩阵a的元素a
ij
赋值给矩阵m；因此，本发明的矢量化方法用少量加法运算及赋值运算取代了传统克罗内克积法的大量加法与乘法运算，具有更高的计算效率。
[0034]
若采用硬件电路搭建根据式(5)，把a
11
,
…
,a
1n
,
…
,a
n1
,
…
,a
nn
传送给对应的其中，为m的对角线元素，需要配置加法器，而m的非对角线元素则直接用接线相连。因此，本发明的矢量化方法通过少量加法器即可得到的所有元素，而传统的克罗内克积法需要大量乘法器和加法器，增加了硬件电路接线复杂度及故障率，且不利于减小硬件体积。
[0035]
进一步地，步骤s2中，矢量化功能模块还需要将n阶矩阵q和w展开为n2阶的一维向量vecq和vecw；
[0036]
若采用软件算法搭建vecq和vecw，则根据vec(q)＝[q
11
,
…
,q
1n
,q
21
,
…
,q
n1
,
…
,qnn
]
t
和vec(w)＝[w
11
,
…
,w
1n
,w
21
,
…
,w
n1
,
…
,w
nn
]
t
分别对q和w进行矩阵维数变换即可，式中的q
ij
和w
ij
分别表示矩阵q和w中第i行第j列的元素，且1≤i≤n，1≤j≤n；
[0037]
若采用硬件电路搭建vecq和vecw，则令若采用硬件电路搭建vecq和vecw，则令分别把q
11
,
…
,q
nn
、w
11
,
…
,w
nn
通过接线传送给对应的即可；
[0038]
根据接收的参数矩阵a、q和w，将矢量化结果传送到神经网络功能模块。
[0039]
进一步地，步骤s3中，由于rnn主要包括梯度神经网络(gnn)和张氏神经网络(znn)，根据两者的定义构建如下模型：
[0040]
对于可观性矩阵，矢量化后的梯度神经网络(gnn)模型为：
[0041][0042]
其中，其中，其中，表示克罗内克积运算，i表示n阶单位矩阵，vec(x(t))＝[x
11
(t),
…
,x
1n
(t),x
21
(t),
…
,x
n1
(t),
…
,x
nn
(t)]
t
，其中，x
ij
(t)表示矩阵x(t)的元素；γ是大于0的常数，是x(t)关于时间的导数，表示gnn的激活函数；
[0043]
矢量化后的张氏神经网络(znn)模型为：
[0044][0045]
其中，是znn的激活函数；
[0046]
同理，对于可控性矩阵，矢量化后的梯度神经网络(gnn)模型为：
[0047][0048]
矢量化后的张氏神经网络(znn)模型为：
[0049][0050]
矢量化的rnn模型构成神经网络功能模块；
[0051]
需要说明的是，rnn的矢量化包括参数矩阵矢量化和网络结构矢量化，前者是步骤s2矢量化功能模块的任务，后者则是步骤s3神经网络功能模块的任务。
[0052]
另外，与矢量化功能模块类似，神经网络功能模块的具体形式可以是软件算法，也
可以是硬件电路。采用软件算法搭建神经网络功能模块，就是根据式(6-9)构建微分方程组模型。采用硬件电路搭建循环神经网络(rnn)功能模块，则是根据式(6-9)，通过积分器、加法器、乘法器等搭建电路模型。
[0053]
接收矢量化功能模块的输出，从而完成rnn的矢量化。
[0054]
进一步地，步骤s4具体如下：
[0055]
s4.1、分别设置解向量vecx(t)和vecg(t)的初值为vecx(0)和vecg(0)，选择合适的激活函数和常数γ；
[0056]
s4.2、驱动循环神经网络(rnn)功能模块得到lyapunov方程的解。
[0057]
进一步地，步骤s5具体如下：
[0058]
s5.1、根据步骤s2-s4，得到可控性矩阵g和可观性矩阵x，对可控性矩阵g和可观性矩阵x进行cholesky分解，得到：
[0059][0060]
其中，zg和z
x
都是下三角矩阵；
[0061]
然后，再进行svd分解，得到：
[0062][0063]
其中，u和v是列正交矩阵，∑＝diag(σ1,σ2…
σn)，σ1≥σ2…
≥σi…
≥σn，σi表示hankel奇异值矩阵第i行的对角线元素；
[0064]
这样，左右投影矩阵通过下式获得：
[0065][0066]
其中，ur和vr分别表示矩阵u和v前r列，
[0067]
根据左右投影矩阵，得到降阶系统的系数矩阵ar、br、cr、dr：
[0068][0069]
s5.2、进一步地，得到降阶系统表达式：
[0070][0071]
其中，且和分别表示降阶系统的状态变量及其微分量。
[0072]
本发明相对于现有技术，具有如下的优点及效果：
[0073]
1、本发明首次以rnn的矢量化运算为切入点，探索出了更高效更简单的rnn矢量化方法。
[0074]
2、本发明所提出的高效矢量化方法，无论是软件仿真还是硬件实现，都比传统的克罗内克积方法更简单、更快速。具体来说，在软件仿真耗时上有明显优势，还可降低硬件电路接线复杂度及故障率，有利于减小硬件体积。
附图说明
[0075]
图1是本发明实施例中一种用于电力系统线性化模型降阶的rnn矢量化方法的流程图。
[0076]
图2是本发明实施例中的gnn硬件电路图。
[0077]
图3是本发明实施例中的97阶算例降阶前后系统时域曲线对比图。
[0078]
图4是本发明实施例中的97阶算例降阶前后系统频域曲线对比图。
具体实施方式
[0079]
下面结合具体实施例对本发明作进一步具体详细描述。
[0080]
实施例：
[0081]
一种用于电力系统线性化模型降阶的rnn矢量化方法，如图1所示，包括以下步骤：
[0082]
s1、基于电力系统线性化模型得到lyapunov方程并传送参数矩阵到矢量化功能模块，具体如下：
[0083]
电力系统的小信号模型为：
[0084][0085]
其中，a、b、c、d均为电力系统线性化模型的常数矩阵，其中，a、b、c、d均为电力系统线性化模型的常数矩阵，分别表示电力系统线性化模型内部的状态变量及其微分量，表示电力系统线性化模型的输入变量，表示电力系统线性化模型的输出变量，n表示状态变量及其微分量的维数，m、p分别表示输入变量和输出变量的维数，表示实数域；δxo一般包括线性化的发电机转子角、发电机角频率及动态元件的状态量，δu和δyo需要根据具体研究对象来确定。
[0086]
基于平衡截断法原理得到电力系统的对偶lyapunov方程：
[0087][0088]
其中，g、x分别为可控性矩阵和可观性矩阵；
[0089]
对于可观性矩阵，lyapunov方程如下：
[0090]at
x(t)+x(t)a＝-q
ꢀꢀ
(3)
[0091]
其中，是常系数对称正定矩阵且q＝c
t
c，是可观性lyapunov方程的解，t为时间变量；
[0092]
需要说明的是，公式(3)与公式(2)的可观性矩阵在形式上有不同，由于本发明采用基于循环神经网络的迭代方法求解可观性矩阵，所以公式(3)的解矩阵x(t)含时间变量。
[0093]
对于可控性矩阵，lyapunov方程如下：
[0094]
ag(t)+g(t)a
t
＝-w
ꢀꢀ
(4)
[0095]
其中，是常系数对称正定矩阵且w＝bb
t
，是可控性lyapunov方程的解；
[0096]
将参数矩阵a、q和w传送到矢量化功能模块。
[0097]
本实施例中，采用pst软件工具箱的16机系统和48机系统算例，分别为35阶算例和97阶算例。其中，δxo包括线性化的发电机转子角、发电机角频率和发电机的励磁电压，δu包括发电机角速度差，δyo包括励磁电压辅助信号。
[0098]
s2、根据矢量化矩阵的结构特点，完成矢量化功能模块的搭建，并将输出信号传送到循环神经网络(rnn)功能模块；
[0099]
矢量化功能模块针对的是参数矩阵矢量化，其核心是计算矢量化功能模块针对的是参数矩阵矢量化，其核心是计算为方便描述，令为方便描述，令表示克罗内克积运算，i表示n阶单位矩阵；经研究，的固有元素构成为：
[0100][0101]
其中，a
ij
表示矩阵a的元素。
[0102]
因此，直接根据的结构特点求取
[0103]
需要指出，矢量化功能模块的具体形式可以是软件算法，也可以是硬件电路，视需求而定。
[0104]
本实施例中，采用软件算法搭建首先根据首先根据分别根据和的元素构成将矩阵a的元素a
ij
赋值给矩阵m。因此，本发明的矢量化方法用少量加法运算及赋值运算取代了传统克罗内克积法的大量加法与乘法运算，具有更高的计算效率。
[0105]
为了比较本发明提出的参数矩阵矢量化方法与传统的克罗内克积方法的性能差异，对前述的35阶、97阶算例进行仿真，将两种方法计算的耗时记录在表1中。
[0106]
表1两种矢量化方法的耗时对比
[0107]
[0108]
表1中，方法a表示本文所提矢量化方法，方法b表示克罗内克积法，耗时比值是方法a的耗时除以方法b的耗时。
[0109]
从表1的数据可知，两个不同算例下，本发明所提矢量化方法的计算效率均高于传统的克罗内克积方法。
[0110]
进一步地，矢量化功能模块还需要将n阶矩阵q和w展开为n2阶的一维向量vecq和vecw；
[0111]
本实施例中，采用软件算法搭建vecq和vecw，根据vec(q)＝[q
11
,
…
,q
1n
,q
21
,
…
,q
n1
,
…
,q
nn
]
t
和vec(w)＝[w
11
,
…
,w
1n
,w
21
,
…
,w
n1
,
…
,w
nn
]
t
分别对q和w进行矩阵维数变换即可，式中的q
ij
和w
ij
分别表示矩阵q和w的元素，且1≤i≤n，1≤j≤n；
[0112]
根据接收的参数矩阵a、q和w，将矢量化结果传送到神经网络功能模块。
[0113]
s3、搭建循环神经网络(rnn)功能模块；
[0114]
由于rnn主要包括梯度神经网络(gnn)和张氏神经网络(znn)，根据两者的定义可构建如下模型：
[0115]
对于可观性矩阵，矢量化后的梯度神经网络(gnn)模型为：
[0116][0117]
其中，其中，其中，表示克罗内克积运算，i表示n阶单位矩阵，vec(x(t))＝[x
11
(t),
…
,x
1n
(t),x
21
(t),
…
,x
n1
(t),
…
,x
nn
(t)]
t
，其中，x
ij
(t)表示矩阵x(t)的元素；γ是大于0的常数，是x(t)关于时间的导数，表示gnn的激活函数；
[0118]
矢量化后的张氏神经网络(znn)模型为：
[0119][0120]
其中，φ是znn的激活函数；
[0121]
同理，对于可控性矩阵，矢量化后的梯度神经网络(gnn)模型为：
[0122][0123]
矢量化后的张氏神经网络(znn)模型为：
[0124]
[0125]
矢量化的rnn模型构成神经网络功能模块。
[0126]
需要说明的是，rnn的矢量化包括参数矩阵矢量化和网络结构矢量化，前者是步骤s2矢量化功能模块的任务，后者则是步骤s3神经网络功能模块的任务。
[0127]
另外，与矢量化功能模块类似，神经网络功能模块的具体形式可以是软件算法，也可以是硬件电路。采用软件算法搭建神经网络功能模块，就是根据式(6-9)构建微分方程组模型。采用硬件电路搭建循环神经网络(rnn)功能模块，则是根据式(6-9)，通过积分器、加法器、乘法器等搭建电路模型。
[0128]
接收矢量化功能模块的输出，从而完成rnn的矢量化。
[0129]
s4、确定计算参数并驱动循环神经网络(rnn)功能模块，从而得到lyapunov方程的解信号，具体如下：
[0130]
(1)分别设置解向量vecx(t)和vecg(t)的初值为vecx(0)和vecg(0)，选择合适的激活函数和常数γ。
[0131]
(2)驱动循环神经网络(rnn)功能模块得到lyapunov方程的解。
[0132]
在本实施例中，以97阶算例的可观性矩阵为例，vecx(0)为全零矩阵，激活函数均选取线性激活函数，γ为100，若采用软件求解，则调用matlab软件的ode45函数。若采用硬件电路，以gnn为例，可利用图2所示的电路图进行求解[5]，其中件电路，以gnn为例，可利用图2所示的电路图进行求解[5]，其中
[0133]
采用软件求解znn的结果如表2所示。
[0134]
表2 97阶算例下的znn求解精度
[0135] 线性激活函数znn范数0.0055
[0136]
表格中的范数是指znn收敛后的||a
t
x+xa+c||f，表征的是znn求解精度。从表2可知，范数非常接近0，说明所得lyapunov解的精度很高。
[0137]
s5、基于步骤s4的解信号，得到电力系统的线性化降阶模型，具体如下：
[0138]
(1)根据步骤s2-s4，得到可控性矩阵g和可观性矩阵x，对可控性矩阵g和可观性矩阵x进行cholesky分解，得到：
[0139][0140]
其中，zg和z
x
都是下三角矩阵；
[0141]
然后，再进行svd分解，得到：
[0142][0143]
其中，u和v是列正交矩阵，∑＝diag(σ1,σ2…
σn)，σ1≥σ2…
≥σi…
≥σn，σi表示hankel奇异值矩阵第i行的对角线元素；
[0144]
这样，左右投影矩阵可以通过下式获得：
[0145][0146]
其中，ur和vr分别表示矩阵u和v前r列，
[0147]
根据左右投影矩阵，就可得到降阶系统的系数矩阵ar、br、cr、dr：
[0148][0149]
(2)更进一步地，可以得到降阶系统表达式：
[0150][0151]
其中，且和分别表示降阶系统的状态变量及其微分量。
[0152]
本实施例中，降阶系统阶数r取为70.
[0153]
为对比降阶前后的系统动态特性，采用0.01p.u的阶跃扰动，比较了原系统模型和降阶系统模型的时域响应和频域响应，分别如图3和图4所示。
[0154]
从图3和图4可知，降阶系统模型可以很好地近似于原系统。
[0155]
实施例2：
[0156]
与实施例1的不同之处在于，本实施例中，采用硬件电路搭建vecq和vecw，令分别把q
11
,
…
,q
nn
、w
11
,
…
,w
nn
通过接线传送给对应的即可。
[0157]
实施例3：
[0158]
与实施例1的不同之处在于，本实施例中，采用硬件电路搭建根据式(5)，把a
11
,
…
,a
1n
,
…
,a
n1
,
…
,a
nn
传送给对应的其中，为m的对角线元素，需要配置加法器，而m的非对角线元素则直接用接线相连。因此，本发明的矢量化方法通过少量加法器即可得到的所有元素，而传统的克罗内克积法需要大量乘法器和加法器，增加了硬件电路接线复杂度及故障率，且不利于减小硬件体积。
[0159]
本领域的技术人员容易理解，上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受上述实施例的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化，均应为等效的置换方式，都包含在本发明的保护范围之内。

技术特征：

1.一种用于电力系统线性化模型降阶的rnn矢量化方法，其特征在于，包括以下步骤：s1、基于电力系统线性化模型得到lyapunov方程并传送参数矩阵到矢量化功能模块；s2、根据矢量化矩阵的结构特点，完成矢量化功能模块的搭建，并将输出信号传送到循环神经网络功能模块；s3、搭建循环神经网络功能模块；s4、确定计算参数并驱动循环神经网络功能模块，从而得到lyapunov方程的解信号；s5、基于步骤s4的解信号，得到电力系统的线性化降阶模型。2.根据权利要求1所述的一种用于电力系统线性化模型降阶的rnn矢量化方法，其特征在于，步骤s1具体如下：电力系统的小信号模型为：其中，a、b、c、d均为电力系统线性化模型的常数矩阵，δx
o
、分别表示电力系统线性化模型内部的状态变量及其微分量，δu表示电力系统线性化模型的输入变量，δy
o
表示电力系统线性化模型的输出变量，n表示状态变量及其微分量的维数，m、p分别表示输入变量和输出变量的维数，表示实数域；δx
o
一般包括线性化的发电机转子角、发电机角频率及动态元件的状态量，δu和δy
o
需要根据具体研究对象来确定；基于平衡截断法原理得到电力系统的对偶lyapunov方程：其中，g、x分别为可控性矩阵和可观性矩阵；对于可观性矩阵，lyapunov方程如下：a
t
x(t)+x(t)a＝-q
ꢀꢀꢀꢀꢀꢀ
(3)其中，q是常系数对称正定矩阵且q＝c
t
c，x(t)是可观性lyapunov方程的解，t为时间变量；对于可控性矩阵，lyapunov方程如下：ag(t)+g(t)a
t
＝-w
ꢀꢀꢀꢀꢀ
(4)其中，w是常系数对称正定矩阵且w＝bb
t
，g(t)是可控性lyapunov方程的解，将参数矩阵a、q和w传送到矢量化功能模块。3.根据权利要求1所述的一种用于电力系统线性化模型降阶的rnn矢量化方法，其特征在于，步骤s2中，矢量化功能模块针对的是参数矩阵矢量化，其核心是计算为方便描述，令为方便描述，令为方便描述，令表示克罗内克积运算，i表示n阶单位矩阵；经研究，的固有元素构成为：
其中，a
ij
表示矩阵a第i行第j列的元素，1≤i≤n，1≤j≤n；因此，直接根据的结构特点求取4.根据权利要求3所述的一种用于电力系统线性化模型降阶的rnn矢量化方法，其特征在于，若采用软件算法搭建首先根据首先根据分别根据和的元素构成将矩阵a的元素a
ij
赋值给矩阵m；若采用硬件电路搭建根据式(5)，把a
11
,
…
,a
1n
,
…
,a
n1
,
…
,a
nn
传送给对应的其中，为m的对角线元素，需要配置加法器，而m的非对角线元素则直接用接线相连。5.根据权利要求2所述的一种用于电力系统线性化模型降阶的rnn矢量化方法，其特征在于，步骤s2中，矢量化功能模块还需要将n阶矩阵q和w展开为n2阶的一维向量vecq和vecw；若采用软件算法搭建vecq和vecw，则根据vec(q)＝[q
11
,
…
,q
1n
,q
21
,
…
,q
n1
,
…
,q
nn
]
t
和vec(w)＝[w
11
,
…
,w
1n
,w
21
,
…
,w
n1
,
…
,w
nn
]
t
分别对q和w进行矩阵维数变换即可，式中的q
ij
和w
ij
分别表示矩阵q和w中第i行第j列的元素，且1≤i≤n，1≤j≤n；若采用硬件电路搭建vecq和vecw，则令若采用硬件电路搭建vecq和vecw，则令分别把q
11
,
…
,q
nn
、w
11
,
…
,w
nn
通过接线传送给对应的即可；根据接收的参数矩阵a、q和w，将矢量化结果传送到神经网络功能模块。6.根据权利要求1所述的一种用于电力系统线性化模型降阶的rnn矢量化方法，其特征在于，步骤s3中，由于循环神经网络主要包括梯度神经网络和张氏神经网络，根据两者的定义构建如下模型：对于可观性矩阵，矢量化后的梯度神经网络模型为：
其中，其中，其中，表示克罗内克积运算，i表示n阶单位矩阵，vec(x(t))＝[x
11
(t),
…
,x
1n
(t),x
21
(t),
…
,x
n1
(t),
…
,x
nn
(t)]
t
，其中，x
ij
(t)表示矩阵x(t)的元素；γ是大于0的常数，是x(t)关于时间的导数，表示梯度神经网络的激活函数；矢量化后的张氏神经网络模型为：其中，是张氏神经网络的激活函数。7.根据权利要求6所述的一种用于电力系统线性化模型降阶的rnn矢量化方法，其特征在于，对于可控性矩阵，矢量化后的梯度神经网络模型为：矢量化后的张氏神经网络模型为：矢量化的循环神经网络模型构成神经网络功能模块。8.根据权利要求7所述的一种用于电力系统线性化模型降阶的rnn矢量化方法，其特征在于，若采用软件算法搭建神经网络功能模块，根据式(6)～式(9)构建微分方程组模型；若采用硬件电路搭建循环神经网络功能模块，则是根据式(6)～式(9)搭建电路模型，电路模型包括积分器、加法器和乘法器；接收矢量化功能模块的输出，从而完成循环神经网络的矢量化。9.根据权利要求1～8任一项所述的一种用于电力系统线性化模型降阶的rnn矢量化方法，其特征在于，步骤s4具体如下：s4.1、分别设置解向量vecx(t)和vecg(t)的初值为vecx(0)和vecg(0)，选择合适的激活函数和常数γ；s4.2、驱动循环神经网络功能模块得到lyapunov方程的解。10.根据权利要求9所述的一种用于电力系统线性化模型降阶的rnn矢量化方法，其特征在于，步骤s5具体如下：s5.1、根据步骤s2-s4，得到可控性矩阵g和可观性矩阵x，对可控性矩阵g和可观性矩阵x进行cholesky分解，得到：其中，z
g
和z
x
都是下三角矩阵；然后，再进行svd分解，得到：
其中，u和v是列正交矩阵，∑＝diag(σ1,σ2…
σ
n
)，σ1≥σ2…
≥σ
i
…
≥σ
n
，σ
i
表示hankel奇异值矩阵第i行的对角线元素；这样，左右投影矩阵通过下式获得：其中，u
r
和v
r
分别表示矩阵u和v前r列，根据左右投影矩阵，得到降阶系统的系数矩阵a
r
、b
r
、c
r
、d
r
：s5.2、进一步地，得到降阶系统表达式：其中，且和分别表示降阶系统的状态变量及其微分量。

技术总结

本发明公开了一种用于电力系统线性化模型降阶的RNN矢量化方法。所述方法包括以下步骤：基于电力系统线性化模型得到Lyapunov方程并传送参数矩阵到矢量化功能模块；根据矢量化矩阵的结构特点，完成矢量化功能模块的搭建，并将输出信号传送到循环神经网络功能模块；搭建循环神经网络功能模块；确定计算参数并驱动循环神经网络功能模块，从而得到Lyapunov方程的解信号；基于解信号，得到电力系统的线性化降阶模型。本发明相比传统的基于克罗内克积的矢量化方法，在软件耗时上有明显优势，还可降低硬件电路接线复杂度及故障率，有利于减小硬件体积。件体积。件体积。