[发明专利]基于递推随机子空间的电力系统低频振荡在线辨识方法

权利要求书

1.一种基于递推随机子空间的电力系统低频振荡在线辨识方法，利用随机子空间方法辨识电力系统低频振荡模态信息，并引入遗忘因子对Hankel协方差矩阵更新，基于投影近似子空间跟踪方法对子空间进行递推，避免SVD运算，降低了算法复杂度、提升辨识实时性，其特征在于：包括如下步骤，

步骤1：由有限长度的振荡电气量实测数据构造Hankel协方差矩阵：

式中，为振荡电气量实测时间序列数据，，，为数据协方差，为数据窗的长度；

步骤2：对矩阵H做一次SVD分解计算：

将奇异值从大到小排列：

令，，和分别为左、右奇异向量，、和、分别为相应信号子空间和噪声子空间的左右奇异向量，是包含全部奇异值的对角矩阵，包含个主奇异值，为张成的信号子空间，为张成的噪声子空间；

步骤3：进行以下计算，其结果作为步骤4随机子空间迭代递归运算的初始值：

其中，符号表示广义逆；为权值矩阵，为协方差矩阵估计；

步骤4：利用时间序列数据、，广义可观测矩阵的计算更新可由以下基于投影近似子空间跟踪方法的递归计算实现：

其中，遗忘因子可取值0.95-0.99；

步骤5:依据系统的最小实现理论，系统的广义可观测矩阵可表示为：        

其中，C为线性状态空间输出矩阵，为系统矩阵，矩阵由下式求出：

上式中，和分别表示矩阵去除最后一行和第一行后的矩阵；

步骤6:求取矩阵的特征值，设采样时间间隔为，可辨识出各振荡模态的频率衰减因子和阻尼比为：

步骤7:随机子空间模态辨识方法由下式可以直接确定振荡模态振型：

上式中，和分别表示算法模型振型和系统振型，也是矩阵特征值对应的右特征向量；

步骤8:通过最小二乘法求取各模态的幅值和相角；对于个采样数据构造向量，模态幅值和相角辨识通过求解以下线性系统方程：

其中，，最小二乘解为：

由此，各振荡模态幅值和相位为：

步骤9：通过时间序列数据的更新，进行步骤4-8递推计算实现振荡模态信息的在线辨识。