一种考虑量测信号丢失的动态振荡信号参数辨识方法
【专利摘要】本发明公开了一种考虑量测信号丢失的动态振荡信号参数辨识方法。该方法首先采用二项分布建立了量测信号丢失的数学模型；然后，在此基础上设计出了适用于量测信号丢失情形下动态振荡信号参数辨识的方法，即改进的扩展卡尔曼滤波算法。该算法在设计时不仅考虑了量测信号丢失的情况，而且引入了系统噪声和量测噪声协方差矩阵的动态计算方法，可以有效的避免因噪声协方差矩阵设置不当而引起的算法失效，提高了算法的效率。该算法因考虑了实际工程背景，且简单方便，具有较高的工程应用价值。
【专利说明】
-种考虑量测信号丢失的动态振荡信号参数辨识方法
技术领域
[0001] 本发明设及一种考虑量测信号丢失的动态振荡信号参数辨识方法，属于信号分析 与参数辨识技术领域。
【背景技术】
[0002] 近年来，随着现代电网规模不断扩大，电网互联程度的日益提高，系统受到大、小 扰动后产生的动态振荡已经成为制约电网安全稳定运行的最主要的因素之一。由于运些动 态振荡信号可W提供关于电力系统运行模式的重要信息，所W发现并准确掌握运些振荡信 号特征对于电力系统安全稳定运行具有重要意义。
[0003] 鉴于动态振荡信号辨识的重要性，研究人员提出许多辨识方法，如最大似然法、普 龙尼法，卡尔曼滤波方法等。但是，运些方法在设计时未考虑量测信号发生数据丢包情况。 因此，在量测信号发生丢包时，运些方法无法有效的对动态振荡信号的参数进行有效的辨 识。所W，研究量测信号丢包的动态振荡信号参数辨识方法具有重要现实意义和工程应用 价值。
【发明内容】
[0004] 发明目的：针对现有技术中存在的问题，为了有效地解决量测信号丢包时动态振 荡信号的参数辨识，克服已有辨识算法的缺点，本发明设计了一种考虑量测信号丢失的动 态振荡信号参数辨识方法，有效的实现了量测信号丢包情况下的动态振荡信号参数辨识。
[0005] 技术方案:一种考虑量测信号丢失的动态振荡信号参数辨识方法，该方法在计算 机中是依次按照如下步骤实现的：
[0006] (1)、获取量测信号丢失的离散状态空间表达式，模型公式为：
[0007]
[000引巧中，X化)表不k时刻的状态问量，y化)表示k时刻的输出量测向量，f( ?)和K ?) 是对应于具体问题中的非线性函数，W化)和V化)分别是系统噪声和量测噪声，曰化)是符合 二项分布的白噪声随机序列，即，S(k)=diag{ 丫 i(k), 丫 2(k),..., 丫 m(k)}，丫 i(k)(i = 1.. .m)是不相关的随机变量，且与W化)和V化），X(O)不相关。
[0009] (2)、丫 1化)的取值为0或者1，且满足如下计算公式：
[0010] prob( Ti(k)=0) = I-Ui(k)
[0011] prob(丫 i(k) = l)=]ii化）
[0012] 式中，prob(丫 1化）=0) = 1-化化)表示的是第i个量测量的丢失率。
[0013] (3)、初始化，包括:设定参数辨识的初值4，初始参数辨识误差协方差与，W及过程 噪声和量测噪声所满足的初始协方差矩阵Qo和Ro,整体算法迭代次数最大值S;
[0014] (4)、获取量测数据丢包情况下的量测数据y化）；
[0015] (5)、由已知条件，计算k时刻的系统噪声和量测噪声所满足的动态协方差矩阵Q 化)和R化），计算公式为：
[0016]
[0017]
[001 引
[0019]
[0020] 式中k代表第k时刻迭代，N是动态估计窗口值，其是一个常数，G化)是k时刻的卡尔 曼滤波增益，是k时刻的状态估计协方差，
良示非线性函数M ?)在馬 处的雅克比矩阵，其中M ?)对应具体问题输出方程中的非线性函数，上标T表示矩阵的装 置。系数頁巧)的计算公式如下：
[0021]
[0022] (6)、由已经得到的k-1时刻的状态估计值和状态估计误差协方差，利用本发明算 法的预测步，得到k时刻的状态预测值和状态预测误差协方差，计算公式为：
[0023]
[0024]
[0025] 式中，-*(*)表示k时刻的状态预测值，f( ?)表示非线性函数，；作-1)表示k-1时刻 的状态估计向量。表示k时刻的状态预测误差协方1
t示非线性函 数f (X)在免批-1)处的雅克比矩阵，户表示k-1时刻的状态估计误差协方差，上标T表示 转置，Q化-1)是系统噪声k-1时刻所满足的动态协方差矩阵。
[0026] (7)、由已知条件计算k时刻的中间迭代变量M化），其计算公式为：
[0027]
[002引式中O符号表示化damard乘积，其定义为[AoB]ij=[AijXBij]，式中系数ミ(A)的取 值计算规则为：
[0029；
[0030；
[0031] (8)、利用k时刻的中间迭代值M化)和k时刻的预测误差协方差计算k时刻的最优滤 波增益，计算步骤为：
[0032]
[0033] 式中，G化)表示k时刻的最优滤波增益，#(!；>表示k时刻的状态预测误差协方差，上 标T表示转置：
g示非线性函数Mx)在处的雅克比矩阵。
[0034] (9)、利用k时刻的最优滤波增益和k时刻的预测误差协方差计算k时刻的估计误差 协方差，计算步骤为：
[0035]
[0036] (10)、结合已得到的k时刻的最优滤波增益和k时刻的状态预测值计算k时刻的状 态估计值.1+當击驢九，
[0037]
[0038] (11)、按照上述步骤，进行多次迭代辨识，若k《S，则迭代继续，若k>S，则迭代结 束，输出辨识结果。
【附图说明】
[0039] 图1为本发明实施例的方法流程图；
[0040] 图2实施例的动态振荡信号；
[0041 ]图3为实施例采用本发明方法^及邸。算法的信号频率辨识结果对比；
[0042] 图4为实施例采用本发明方法W及邸。算法的信号阻尼因子辨识结果对比；
[0043] 图5为实施例采用本发明方法信号频率和阻尼因子的辨识误差。
【具体实施方式】
[0044] 下面结合具体实施例，进一步阐明本发明，应理解运些实施例仅用于说明本发明 而不用于限制本发明的范围，在阅读了本发明之后，本领域技术人员对本发明的各种等价 形式的修改均落于本申请所附权利要求所限定的范围。
[0045] 如图1所示，动态振荡信号参数辨识方法，其包含如下步骤：
[0046] (1)、获取量测信号丢失的离散状态空间表达式；
[0047] (2)、初始化，包括:设定参数辨识的初值^，初始参数辨识误差协方差皆，W及过程 噪声和量测噪声所满足的初始协方差矩阵Qo和Ro,整体算法迭代次数最大值S;
[0048] (3)、获取量测数据丢包情况下的量测数据y化）；
[0049] (4)、由已知条件，计算k时刻的系统噪声和量测噪声所满足的协方差矩阵Q化)和R (k)；
[0050] (5)、由已经得到的k-1时刻的状态估计值和状态估计误差协方差，利用本发明算 法的预测步，得到k时刻的状态预测值和状态预测误差协方差；
[0051 ] (6)、由已知条件计算k时刻的中间迭代变量M化）；
[0052] (7)、利用k时刻的中间迭代值M化)和k时刻的预测误差协方差计算k时刻的最优滤 波增益；
[0053] (8)、利用k时刻的最优滤波增益和k时刻的预测误差协方差计算k时刻的估计误差 协方差；
[0054] (9)、结合已得到的k时刻的最优滤波增益和k时刻的状态预测值计算k时刻的状态 估计值；
[0055] (10)、按照上述步骤，进行多次迭代辨识，若k《S，则迭代继续，若k>S，则迭代结 束，输出辨识结果。
[0056] -般情况下动态振荡信号可W表示为多个指数衰减的正弦信号的之和，可W描述 为如下形式：
[0化7]
[005引式中，Al,Si,Wi, (61是实数的未知参数，n(t)是一个零均值白噪声。其中，Si称为动 态信号的阻尼因子，Wi是动态信号的频率，其中Wi, Si为待估参数。经过推理可W得到动态信 号的状态变量分量中包含待估参数的离散状态空间模型。考虑由N个指数衰减的正弦信号 总和组成的动杰倍号.Jt4N个献杰巧畳形式可W表达如下：
[0化9]
[0060]
[0061]
[0062]
[0063]式中i代表运些变量和参数是属于动态信号的第i个衰减正弦信号。k代表时刻，fs 代表采样频率。根据推理可得到k+1时刻的状态分量：
[0066] X4i-1, k+1 = X4i-1, k+W4i-l, k[0067] X4i,k+l = X4i,k+W4i,k[0068] 则考虑量测数据丢包的输出方程为：
[0064
[00化
[0069 V (-A
/
[0070] 式中，k2i-i = cos( 4 i)，k2i = -sin( 4 i)，nk为均值为零的白噪声，刮k)为对应维度 的二项分布随机序列。
[0071] 所W，考虑量测数据丢包的动态振荡信号的状态空间模型一般可W表示为：
[0072]
[0073] 式中，f( ?)和K ?)代表可W根据泰勒级数展开进行线性化的非线性函数，w(k) 和V化)是均值为零的高斯白噪声序列，分别满足动态协方差矩阵Q化)和R化)。具体而言：
[0074]
[0075] 而函数h(xk)可W表示为如下形式：
[0076] H=化iksOO... ,k2i-ik2i00... ,k2N-ik2N〇0)
[0077] hUk)=Hxk
[0078] 至此，考虑量测数据丢包的动态振荡信号状态空间表达式已建立。在此基础上，贝U 可W运用本发明所介绍的方法，进行量测信号丢包情况下的动态振荡信号参数辨识，得到 辨识结果。
[0079] 下面介绍本发明的一个实施例：
[0080] 考虑动态振荡信号为：
[0081 ] y(k)=巧k)e_〇'〇ik sin(0.4k)+nk，0《k《300
[0082] 式中k是信号采样时刻，nk是高斯白噪声。如图2所示，该动态信号是由一个指数衰 减的正弦信号组成，该动态振荡信号的频率为W = O.4,阻尼因子为5 = 0.01。在运用本发明 所提出的方法进行动态信号参数辨识时，所采用的相关初始参数值为：
[0083]
[0084]
[0085] 在对本算例进行辨识时，假设量测信号的丢包率11 = 0.1，在计算噪声所满足的动 态误差协方差矩阵时，动态估计窗口值N=2。
[0086] 图1为实施例所用的算法流程图，图2为实施例的动态振荡信号，图3为实施例采用 本发明方法W及邸尸算法的信号频率W辨识结果对比，图4为实施例采用本发明方法W及邸F 算法的信号阻尼因子S辨识结果对比，图5为实施例采用本发明方法信号频率和阻尼因子的 辨识误差。对比结果表明，本发明所提的方法有效的实现了量测信号丢失情形下的动态振 荡信号参数辨识，且较邸F算法具有更好的收敛性和精度。
【主权项】
1. 一种考虑量测信号丢失的动态振荡信号参数辨识方法，其特征在于，包含如下步骤： (1 )、获取量测信号丢失的离散状态空间表达式； (2) 、初始化，包括:设定参数辨识的初值，初始参数辨识误差协方差為.，以及过程噪声 和量测噪声所满足的初始协方差矩阵qq和r q，整体算法迭代次数最大值S; (3) 、获取量测数据丢包情况下的量测数据y(k); (4) 、由已知条件，计算k时刻的系统噪声和量测噪声所满足的协方差矩阵Q(k)和R(k); (5) 、由已经得到的k-1时刻的状态估计值和状态估计误差协方差，利用本发明算法的 预测步，得到k时刻的状态预测值和状态预测误差协方差； (6) 、由已知条件计算k时刻的中间迭代变量M(k); (7) 、利用k时刻的中间迭代值M(k)和k时刻的预测误差协方差计算k时刻的最优滤波增 益； (8) 、利用k时刻的最优滤波增益和k时刻的预测误差协方差计算k时刻的估计误差协方 差； (9) 、结合已得到的k时刻的最优滤波增益和k时刻的状态预测值计算k时刻的状态估计 值； (10) 、按照上述步骤，进行多次迭代辨识，若S，则迭代继续，若k> S，则迭代结束，输 出辨识结果。2. 如权利要求1所述的考虑量测信号丢失的动态振荡信号参数辨识方法，其特征在于， 获取量测信号丢失的离散状态空间表达式，模型公式为：式中，x(k)表示k时刻的状态向量，y(k)表示k时刻的输出量测向量，f( ·)和h( ·)是对 应于具体问题中的非线性函数，w(k)和v(k)分别是系统噪声和量测噪声，Ξ(1〇是符合二项 分布的白噪声随机序列，BP，S(k)=diag{ γι(1〇, Y2(k)，···，Ym(k)}，Yi(k)(i = l. · .m)是 不相关的随机变量，且与w(k)和v(k)，x(0)不相关。 (2)、丫1(1〇的取值为0或者1，且满足如下计算公式： prob( Yi(k)=0) = l-yi(k) prob( γ i(k) = 1) =yi(k) 式中，prob( γ i(k)=0) = l-yi(k)表示的是第i个量测量的丢失率。3. 如权利要求2所述的考虑量测信号丢失的动态振荡信号参数辨识方法，其特征在于， 由已知条件，计算k时刻的系统噪声和量测噪声所满足的动态协方差矩阵Q(k)和R(k)，计算 公式为：式中k代表第k时刻迭代，N是动态估计窗口值，其是一个常数，G (k)是k时刻的卡尔曼滤 波增益，#0；)是k时刻的状态估计协方差表示非线性函数h( ·)在毛处 的雅克比矩阵，其中h( ·)对应具体问题输出方程中的非线性函数，上标T表示矩阵的装置。 系数g(fc)的计算公式如下：4. 如权利要求3所述的考虑量测信号丢失的动态振荡信号参数辨识方法，其特征在于， 由已经得到的k-Ι时刻的状态估计值和状态估计误差协方差，利用本发明算法的预测步，得 到k时刻的状态预测值和状态预测误差协方差，计算公式为：式中，50)表示k时刻的状态预测值，f( ·)表示非线性函数，表示k-Ι时刻的状 态估计向量;表示k时刻的状态预测误差协方差，示非线性函数f (X)在;^-1)处的雅克比矩阵，#0-1)表示k-1时刻的状态估计误差协方差，上标T表示转 置，Q(k-l)是系统噪声k-Ι时刻所满足的动态协方差矩阵； 由已知条件计算k时刻的中间迭代变量M(k)，其计算公式为：式中〇符号表示Hadamard乘积，其定义为[A〇B]ij = [Aij XBij]，式中系数的取值计 算规则为：5. 如权利要求4所述的考虑量测信号丢失的动态振荡信号参数辨识方法，其特征在于， 利用k时刻的中间迭代值M(k)和k时刻的预测误差协方差计算k时刻的最优滤波增益，计算 步骤为：式中，G(k)表示k时刻的最优滤波增益，户("表示k时刻的状态预测误差协方差，上标T 表示转置表示非线性函数h(x)在处的雅克比矩阵； 利用k时刻的最优滤波增益和k时刻的预测误差协方差计算k时刻的估计误差协方差， 计算步骤为：结合已得到的k时刻的最优滤波增益和k时刻的状态预测值计算k时刻的状态估计值， 计算步骤为：
【文档编号】G06K9/00GK105956565SQ201610302188
【公开日】2016年9月21日
【申请日】2016年5月9日
【发明人】王 义, 孙永辉, 卫志农, 孙国强, 张世达, 郭敏
【申请人】河海大学