一种考虑量测信号丢失的改进集合卡尔曼滤波估计方法技术

本发明专利技术公开了一种考虑量测信号丢失的改进集合卡尔曼滤波估计方法。该方法首先采用二项分布建立了量测信号丢失的数学模型；然后，通过对现有的集合卡尔曼滤波状态估计方法进行改进，设计出了改进集合卡尔曼滤波状态估计算法。该算法因考虑了实际工程背景中广泛存在的量测信号丢失的情况，且简单易于实现，有效的解决了量测信号丢失下的非线性系统的状态估计和参数辨识问题，具有一定的工程应用价值。

【技术实现步骤摘要】

本专利技术涉及一种考虑量测信号丢失的改进集合卡尔曼滤波估计方法，属于信号分析与状态估计

技术介绍
近年来，随着现代电网规模不断扩大，电网互联程度的日益提高，为了能够更好的研究电力系统的运行情况，如何快速有效的实现对电力系统的动态状态估计问题显得尤为重要，因为，通过分析这些状态信号，可以有效的分析电网运行的状态，对于电力系统安全稳定运行具有重要意义。然而，在运用最小二乘法，卡尔曼滤波等方法对电网的运行状态进行估计时，大都忽略了对采样信号数据发生丢包情形的考虑。在实际应用中，把传感器采集的现场信号传送到中央控制室，不可避免的要发生数据在信道的丢失情形。为了更加准确有效的实现对电网状态信号的状态估计，有必要对考虑量测信号丢失情况下的状态估计算法进行研究。
技术实现思路
专利技术目的：基于以上分析，为了有效的解决量测信号丢失情形下的状态估计问题，本专利技术设计了考虑量测信号丢失情形下的状态估计方法。本专利技术首先采用二项分布序列，实现了对量测信号丢失情形的状态空间模型建模，然后在此基础上以传统集合卡尔曼滤波为基本工具，推导出了用于量测信号丢失情形下解决状态估计问题的改进集合卡尔曼滤波。技术方案：一种考虑量测信号丢失的改进集合卡尔曼滤波状态估计方法，该方法在计算机中是依次按照如下步骤实现的：(1)、获取量测信号丢失的离散状态空间表达式，模型公式为： x k = f ( x k - 1 ) + w k y k = Ξ ( k ) h ( x k ) + v k ]]>式中，xk表示k时刻的状态向量，yk表示k时刻的输出量测向量，f(·)和h(·)是对应于具体问题中的非线性函数，wk和vk分别是系统噪声和量测噪声，Ξ(k)是符合二项分布的白噪声随机序列，即，
γi(k)(i＝1…m)是不相关的随机变量，且与wk和vk，x(0)不相关。(2)、γi(k)的取值为0或者1，且满足如下计算公式：prob(γi(k)＝0)＝1-μi(k)prob(γi(k)＝1)＝μi(k)式中，1-μi(k)是第i个量测量的丢失率。(3)、初始化，包括：设定丢包率的取值，样本采样数q，状态估计的初始值以及过程噪声和量测噪声所满足的协方差矩阵Q和R，整体算法的最大迭代次数L。(4)、由已知的k-1时刻的状态估计值，利用改进集合卡尔曼滤波的预测步，得到k时刻的q个采样集合的状态预测值计算公式为： x k f i = f ( x k a i , u k ) + w k i ]]> X k f = ( x k f 1 , ... , x k f q ) ]]>式中，表示k时刻的第i个样本的状态预测值。(5)、在上一步的基础上，计算样本采样预测均值和样本采样预测输出值，计算公式如下： x ‾ k f = 1 / q Σ i = 1 q x k f i ]]> y k f i = Ξ ‾ ( k ) h ( x k f i ) ]]> y ‾ k f = 1 / q Σ i = 1 q y k f i ]]>式中，是q个样本k时刻预测值的均值，是第i个样本在k时刻丢包情形下的量测输出值，是q个样本k时刻量测输出值的均值。上式中系数的取值遵循如下规律： Ξ ‾ ( k ) 本文档来自技高网...

【技术保护点】
一种考虑量测信号丢失的改进集合卡尔曼滤波估计方法，其特征在于，包含如下步骤：(1)、获取量测信号丢失的离散状态空间表达式，模型公式为：xk=f(xk-1)+wkyk=Ξ(k)h(xk)+vk]]>式中，xk表示k时刻的状态向量，yk表示k时刻的输出量测向量，f(·)和h(·)是对应于具体问题中的非线性函数，wk和vk分别是系统噪声和量测噪声，Ξ(k)是符合二项分布的白噪声随机序列，即，Ξ(k)＝diag{γ1(k)，γ2(k)，·，·γm(k)}，γi(k)(i＝1...m)是不相关的随机变量，且与wk和vk，x(0)不相关。(2)、γi(k)的取值为0或者1，且满足如下计算公式：prob(γi(k)＝0)＝1‑μi(k)prob(γi(k)＝1)＝μi(k)式中，1‑μi(k)是第i个量测量的丢失率。(3)、初始化，包括：设定丢包率的取值，样本采样数q，状态估计的初始值以及过程噪声和量测噪声所满足的协方差矩阵Q和R，整体算法的最大迭代次数L。(4)、由已知的k‑1时刻的状态估计值，利用改进集合卡尔曼滤波的预测步，得到k时刻的q个采样集合的状态预测值计算公式为：xkfi=f(xkai,uk)+wki]]>Xkf=(xkf1,...,xkfq)]]>式中，表示k时刻的第i个样本的状态预测值。(5)、在上一步的基础上，计算样本采样预测均值和样本采样预测输出值，计算公式如下：x‾kf=1/qΣi=1qxkfi]]>ykfi=Ξ‾(k)h(xkfi)]]>y‾kf=1/qΣj=1qykfi]]>式中，是q个样本k时刻预测值的均值，是第i个样本在k时刻丢包情形下的量测输出值，是q个样本k时刻量测输出值的均值。上式中系数的取值遵循如下规律：Ξ‾(k)=E{Ξ(k)}=diag{μ1(k)μ2(k)...μm(k)}]]>(6)、在上一步基础上，利用改进集合卡尔曼滤波的预测步，计算k时刻的集合状态误差矩阵以及集合输出误差矩阵计算公式为：Ekf=[xkf1-x‾kf...xkfq-x‾kf]]]>Eykf=[ykf1-y‾kf...ykfq-y‾kf]]]>(7)、计算k时刻的交叉协方差矩阵和计算公式如下：pxykf=1q-1Ekf(Eykf)T]]>pyykf=1q-1Eykf(Eykf)T]]>上式中括号外上标T表示矩阵的转置。(8)、依据改进集合卡尔曼滤波的估计步，计算k时刻的状态估计均值其计算步骤为：Kk=pxykf(pxykf)-1]]>xkai=xkfi+Kk(yk-Ξ‾(k)h(xkfi))]]>x‾ka=1/qΣi=1qxkai]]>式中，Kk是k时刻的改进卡尔曼滤波增益，是第i个样本k时刻的状态估计值，是k时刻q个样本状态估计均值，其值作为k时刻对系统的状态估计值。(9)、按照上述步骤，进行多次迭代估计，至k＞L，迭代停止，输出状态估计结果。...

【技术特征摘要】
1.一种考虑量测信号丢失的改进集合卡尔曼滤波估计方法，其特征在于，包含如下步骤：(1)、获取量测信号丢失的离散状态空间表达式，模型公式为： x k = f ( x k - 1 ) + w k ...

【专利技术属性】
技术研发人员：沈谋全，周玉庭，孙世荣，严沈，陈爱华，
申请(专利权)人：南京工业大学，
类型：发明
国别省市：江苏;32