基于H∞扩展卡尔曼滤波的低频振荡信号参数辨识方法技术

本发明专利技术提供了一种基于H∞扩展卡尔曼滤波的低频振荡信号参数辨识方法，利用H∞滤波理论，在低频振荡信号参数辨识时，有效的计及模型不确定性的影响，避免了由于模型参数不确定性所引起的参数辨识误差，且由于采用噪声协方差矩阵自适应技术，动态调整协方差矩阵，从而使得所提方法具有更强的鲁棒性，利于获取更准确的低频振荡信号参数辨识结果。

Parameter identification method of low frequency oscillating signal based on H infinity extended Calman filter

The present invention provides an extended Calman filter based on H for low frequency oscillation signal parameter identification method, using H filtering theory, the low frequency oscillation signal parameter identification, considering the effect of model uncertainty, to avoid the error of parameter identification of model parameter uncertainty caused by, and due to the adaptive noise the dynamic adjustment of the covariance matrix, the covariance matrix, so that the robustness of the proposed method is more conducive to obtain the low frequency oscillation signal, the parameter identification results more accurate.

【技术实现步骤摘要】
基于H∞扩展卡尔曼滤波的低频振荡信号参数辨识方法
本专利技术涉及一种电力系统，具体涉及一种电力系统低频振荡信号提取方法。
技术介绍
近年来，在全国电网互联与西电东输的过程中，电力交换更加频繁，电力系统的安全稳定问题大多表现为低频振荡。因此，如何有效的提取电力系统低频振荡信号所表征的信息，对于电力系统的安全稳定分析具有重要意义。在目前的研究中，利用现场实测数据分析处理信号，得到振荡特征参数是研究电力系统低频振荡的一种有效途径。其常用的方法主要包括实时快速傅里叶变换(fastFouriertransform，FFT)、小波算法、Prony算法和扩展卡尔曼滤波算法等。实时FFT的精度受数据窗限制，不能反映振荡的阻尼特性；小波算法可以反映信号的时变特性，但存在小波基难以选取的问题；Prony算法能直接提取幅值、相位、频率和衰减因子，算法简便，因此被广泛用于电力系统低频振荡模式的识别。但是，Prony算法对噪声较敏感，识别含噪低频振荡信号时的误差较大；当振荡模式为多阶且采样率增大时，识别振荡幅值和初相的计算量呈指数增加，矩阵求逆运算困难。基于扩展卡尔曼滤波(extendedKalmanfilter，EKF)的低频振荡信号参数辨识方法，不仅具有在线辨识功能，而且计算占用内存低，因此应用较广。但是，需要注意的是，EKF方法无法考虑振荡信号建模过程中所引入的不确定性，并且其辨识结果易受噪声初始方差矩阵的影响。
技术实现思路
专利技术目的：本专利技术的目的在于针对现有技术存在的问题，提出了一种基于H∞扩展卡尔曼滤波的低频振荡信号参数辨识方法，其有效的计及模型不确定性的影响，且动态调整协方差矩阵，实现电力系统低频振荡信号参数的准确辨识。技术方案：本专利技术提供了一种基于H∞扩展卡尔曼滤波的低频振荡信号参数辨识方法，包括以下步骤：(1)设定滤波相关的初始值，包括k＝0时刻的状态估计初始值状态估计误差协方差系统噪声和量测噪声协方差矩阵的初始值Q0和R0、移动窗口值L以及最大估计时刻N；(2)获取电力系统低频振荡信号量测序列输入值yk，其量测函数定义如下：yk＝h(xk)+vk式中h(·)表示已知的量测函数，xk为k时刻的参数真值，vk表示k时刻的量测噪声，其满足的协防差矩阵为Rk；(3)计算k时刻的参数预测值计算公式如下：式中f(·)表示已知的系统函数，为k-1时刻的参数估计值；(4)计算k时刻的参数预测误差协方差计算公式如下：式中表示非线性函数f(·)在处的雅克比矩阵，Qk-1表示k-1时刻的系统噪声协方差矩阵；(5)计算k时刻H∞扩展卡尔曼滤波增益Gk，计算公式如下：式中(·)-1为求矩阵的逆运算，表示的非线性函数h(·)在处的雅克比矩阵；(6)计算k时刻的估计误差协方差计算公式如下：式中I为对应维度的单位矩阵，Re，k计算公式如下：其中参数γ设置的计算公式为：式中λ为待设置的大于1的正参数，电力系统参数辨识时取值区间为[1，30]，eig(·)表示取相应矩阵的特征值，max(·)表示取最大值；(7)计算k时刻的参数估计值计算公式如下：式中yk为k时刻的量测值；(8)计算新息序列，计算公式如下：(9)取移动窗口大小为L时，计算窗口内新息序列sk的平均值，即新息矩阵Cvk，其计算公式如下：(10)在上一步的基础上，动态计算k+1时刻系统噪声协方差矩阵Qk和量测噪声协防差矩阵Rk，计算公式如下：(11)按照(3)-(10)步骤依据时间序列进行电力系统低频振荡信号参数辨识，直至k+1＞N时迭代停止，输出参数辨识结果。有益效果：本专利技术利用H∞滤波理论，在低频振荡信号参数辨识时，有效的计及模型不确定性的影响，避免了由于模型参数不确定性所引起的参数辨识误差，且由于采用噪声协方差矩阵自适应技术，动态调整协方差矩阵，从而使得所提方法具有更强的鲁棒性，利于获取更准确的低频振荡信号参数辨识结果。附图说明图1为本专利技术方法的流程图；图2为实施例的电力系统低频振荡信号量测值；图3为实施例对低频振荡信号频率参数w辨识结果；图4为实施例对低频振荡信号阻尼因子参数δ辨识结果；图5为实施例对低频振荡信号参数辨识的均方根偏差。具体实施方式下面对本专利技术技术方案进行详细说明，但是本专利技术的保护范围不局限于所述实施例。一般情况下电力系统低频振荡信号可以表示为多个指数衰减的正弦信号的之和，可以描述为如下形式：式中，Ai，δi，wi，φi是实数的未知参数，N是构成一个振荡信号衰减的正弦信号总数，下标i表示相关参数属于构成低频振荡信号的第i个衰减的余弦信号，n(t)是一个零均值白噪声。其中，δi称为低频振荡信号的阻尼因子，wi表示低频振荡信号的频率，wi，δi为待估参数，经过推理可以得到低频振荡信号的状态变量分量中包含待估参数的离散状态空间模型。考虑由N个指数衰减的正弦信号总和组成的电力系统低频振荡信号，其4N个状态变量形式可以表达如下：x4i-1，k＝wix4i，k＝δi式中i代表这些变量和参数是属于电力系统低频振荡信号的第i个衰减正弦信号(i＝1…N)，k代表时刻，fs代表采样频率。根据推理可得到k+1时刻的状态分量：x4i-1，k+1＝x4i-1，k+w4i-1，kx4i，k+1＝x4i，k+w4i，k则其输出方程为：式中，k2i-1＝cos(φi)，k2i＝-sin(φi)，nk为均值为零的白噪声，所以，电力系统低频振荡信号的状态空间模型一般可以表示为：式中，f(·)和h(·)代表可以根据泰勒级数展开进行线性化的非线性函数，xk+1表示k+1时刻的状态量和待辨识参数分量，yk为k时刻的低频振荡信号量测序列，wk和vk是均值为零的高斯白噪声序列，分别满足协方差矩阵Qk和Rk。具体而言，电力系统低频振荡信号中：而函数h(xk)可以表示为如下形式：H＝(k1k200…，k2i-1k2i00…，k2N1k2N00)h(xk)＝βxk式中β表示已知常系数矩阵。至此，状态变量分量中包含电力系统低频振荡信号模型待估参数的状态空间模型已经建立，在此基础之上，则可以运用本专利技术所介绍的方法，对电力系统低频振荡信号参数辨识。实施例：为了验证本专利技术方法的有效性和实用性，本实施例选取如下电力系统低频振荡信号进行参数辨识分析y(t)＝e-δtCos(wt+φ)+n(t)该低频振荡信号是由一个指数衰减的正弦信号组成。该低频振荡信号待辨识的参数：阻尼因子δ＝0.01，频率w＝0.5rad/s，φ＝0，n(t)是高斯白噪声，其所满足的协方差矩阵为r＝10-5，取采样时间为T＝1s，本实施例在进行仿真实验时取前300采样时刻量测值进行算法验证，即N为300。在运用本专利技术所提出的方法对实施例低频振荡信号进行参数辨识时，所取滤波初始估计误差协方差和系统噪声协方差矩阵初始值为：参数辨识的初始值选取为量测噪声协方差矩阵初始值设置为真实值的103倍，即R0＝10-2；过程噪声动态估计窗口值L取为10，λ取值为20。如图1所示，运用本专利技术方法对电力系统低频振荡信号参数辨识方法，其包含如下步骤：(1)设定滤波相关的初始值，包括k＝0时刻的状态估计初始值状态估计误差协方差系统噪声和量测噪声协方差矩阵的初始值Q0和R0、移动窗口值L以及最大估计时刻N；(2)获取电力系统低频振荡信号量测序列输入值yk，其量测函数定义如下y本文档来自技高网...

【技术保护点】
一种基于H∞扩展卡尔曼滤波的低频振荡信号参数辨识方法，其特征在于：包括以下步骤：(1)设定滤波相关的初始值，包括k＝0时刻的状态估计初始值

【技术特征摘要】
1.一种基于H∞扩展卡尔曼滤波的低频振荡信号参数辨识方法，其特征在于：包括以下步骤：(1)设定滤波相关的初始值，包括k＝0时刻的状态估计初始值状态估计误差协方差系统噪声和量测噪声协方差矩阵的初始值Q0和R0、移动窗口值L以及最大估计时刻N；(2)获取电力系统低频振荡信号量测序列输入值yk，其量测函数定义如下：yk＝h(xk)+vk式中h(·)表示已知的量测函数，xk为k时刻的参数真值，vk表示k时刻的量测噪声，其满足的协防差矩阵为Rk；(3)计算k时刻的参数预测值计算公式如下：式中f(·)表示已知的系统函数，为k-1时刻的参数估计值；(4)计算k时刻的参数预测误差协方差计算公式如下：式中表示非线性函数f(·)在处的雅克比矩阵，Qk-1表示k-1时刻的系统噪声协方差矩阵；(5)计算k时刻H∞扩展卡尔曼滤波增益Gk，计算公式如下：式中(·)-1为求矩阵的逆运算，表示的非线性函数h(·)在处的雅克比矩阵；(6)计算k时刻的估计误差协方差计算公式如下：式中I为对应维度的单位矩阵，Re，k...

【专利技术属性】
技术研发人员：王义，钟永洁，孙永辉，武小鹏，吕欣欣，翟苏巍，
申请(专利权)人：河海大学，
类型：发明
国别省市：江苏,32