一种计及传感器故障的分数阶网络系统状态估计方法技术方案

本发明专利技术公开了一种计及传感器故障的分数阶网络系统状态估计方法，用于分析分数阶网络系统因传感器故障而引起的数据丢包情况下的状态估计问题。该方法的具体步骤如下：首先对传感器故障情况下，发生数据随机丢包的系统模型进行了分析，建立了考虑数据随机丢包情况下的分数阶网络系统模型。接着，以传统的分数阶扩展卡尔曼滤波状态估计方法为基础，结合传感器故障情况下的分数阶网络系统模型，设计出了改进分数阶扩展卡尔曼状态估计方法。本方法适用于传感故障引起数据随机丢包情况下的分数阶网络系统状态估计问题，且流程简单易于实现。

【技术实现步骤摘要】
一种计及传感器故障的分数阶网络系统状态估计方法
本专利技术涉及一种计及传感器故障的分数阶网络系统状态估计方法，属于网络系统分析与控制

技术介绍
网络系统的分析与控制对于保证网络系统安全稳定的运行具有重要的意义。近年来，随着传感器技术的发展，网络系统的实时在线监测与控制成为了众多研究人员所关注的焦点。在现有的研究中，借助于传感器所获取的实时量测信息，通过设计动态的状态估计器，是实现网络系统实时分析与控制的主要途径。一般情况下，现场数据通过传感器进行量测，然后通过信息传输通道传到控制中心，但是需要注意的是，传感器所量测的信息并不总是真实的，因为其会受到外界的干扰，以及信号的衰减，甚至传感器故障的影响。所以，在进行网络系统的动态估计器时必须计及量测信号发生丢包的情况。分数阶网络系统由于可以更加精确的描述系统的结构，近年来被广泛应用于各个领域，如运用在电力系统网络中，可以更加精确的对电力系统中的节点电压和电流进行预测和估计。但是，在现有的分数阶网络系统研究中，计及传感器失败所引起的数据丢包现象主要集中于线性的分数阶网络，而对非线性分数阶网络计及传感器故障的分析和研究，国内外鲜有相关报道。为了进一步拓展分数阶网络的应用，本专利技术设计了计及传感器故障下的非线性分数阶网络系统的状态估计方法，并从理论上给予了证明。最后，实际的分数阶网络系统算例测试验证了本专利技术方法的有效性和实用性。
技术实现思路
专利技术目的：针对现有技术中存在的问题，本专利技术提供一种计及传感器失效的非线性分数阶网络状态估计方法。技术方案：一种计及传感器故障的分数阶网络系统状态估计方法，包括如下部分：1)计及传感器故障的分数阶网络系统建模对于传感器故障下，系统量测数据发生随机丢包的离散非线性分数阶网络系统，其状态方程xk+1和输出方程yk分别为：Δγxk+1＝f(xk)+wkyk＝Γkh(xk)+vk式中：xk+1表示k+1时刻的状态矢量，yk表示k时刻的输出矢量，f(·)和h(·)对应于两个可用泰勒级数展开的非线性函数，wk和vk分别为k时刻的系统噪声值和量测噪声值，二者相互独立无关，满足的协方差矩阵分别为Qk和Rk，式中γj和Γk计算公式如下式中n≥0是分数阶阶次，j≥0代表不同时刻，是符合伯努利分布的二进制标量，其取值为0或1；期望和方差分别为πi，πi(1-πi)，即满足(其中P(·)表示某件事发生的概率)在建立传感器网络失败所引起量测信号丢包的模型之后，则可以通过如下方法对量测信号数据丢包情况下的非线性分数阶网络系统进行状态估计。2)初始化k时刻的估计初始值和估计误差协方差Pk，估计时刻最大值N；式中E(·)表示对某变量进行求期望运算，(·)T表示求矩阵转置。3)计算k时刻的系统函数的雅克比矩阵，计算公式如下式中表示求函数在变量处的偏导。4)计算k时刻的反馈增益矩阵Kk，计算公式如下式中[·]-1为矩阵逆运算符，是阿达玛乘积算子，Υ式中计算公式如下5)计算k+1时刻的状态估计协方差矩阵Pk+1，计算公式如下6)计算k+1时刻的状态估计值计算公式如下7)若k+1＜N则进行下一时刻的迭代估计；反之，则结束迭代，输出估计结果。8)算法的证明过程如下证明：估计误差ek+1可以表示为基于泰勒级数展开可得到：所以估计误差ek+1可以近似等价为：对k+1时刻的状态估计协方差矩阵Pk+1推导化简可得到：式中简化上式中的一项为：并定义如下变量：则k+1时刻的状态估计协方差矩阵Pk+1可进一步简化为：通过完成观测增益的平方项，则：通过联立上述两式，可得到：进而可求得滤波反馈增益：当且仅当估计误差协方差矩阵Pk+1取得最小值，此时协方差矩阵Pk+1为所以k+1时刻的状态估计值可以通过如下公式求取附图说明图1为本专利技术实施例的方法流程图；图2为实施例的量测信号图；图3为实施例运用本专利技术和传统方法的状态估计结果图，其中(a)为本专利技术方法的状态估计结果图，(b)为传统方法的状态估计结果图；图4为实施例采用本专利技术和传统方法的估计误差图，其中(a)为本专利技术方法的估计误差图，(b)为传统方法的估计误差图。具体实施方式下面结合具体实施例，进一步阐明本专利技术，应理解这些实施例仅用于说明本专利技术而不用于限制本专利技术的范围，在阅读了本专利技术之后，本领域技术人员对本专利技术的各种等价形式的修改均落于本申请所附权利要求所限定的范围。如图1所示，计及传感器故障的分数阶网络系统状态估计方法，方法在计算机中是依次按照如下步骤实现的：1)初始化k时刻的估计初始值和估计误差协方差Pk，估计时刻最大值N；2)计算k时刻的系统函数的雅克比矩阵，计算公式如下3)计算k时刻的反馈增益矩阵Kk，计算公式如下4)计算k+1时刻的状态估计协方差矩阵Pk+1，计算公式如下5)计算k+1时刻的状态估计值计算公式如下6)若k+1＜N则进行下一时刻的迭代估计；反之，则结束迭代，输出估计结果。为了验证本专利技术方法的有效性，下面介绍本专利技术的一个实施例，考虑如下非线性分数阶网络系统Δγxk+1＝sin(xk)+wk式中分数阶阶次n＝0.9，因传感器失效所引起的量测信号丢包率为系统噪声wk和量测噪声vk所满足的协方差矩阵分别为在运用本专利技术方法对实施例非线性分数阶网络进行状态估计时，状态估计的初始值x0＝0.9；初始状态估计误差协方差矩阵为P0＝1，最大迭代估计时刻N＝100。分别运用本专利技术计及传感器故障的分数阶网络系统状态估计方法，以及传统分数阶卡尔曼滤波状态估计方法对实施例分数阶网络系统进行变量估计，不同算法的状态估计结果如图3所示，状态估计误差如图4所示。综合图3和图4所示的测试结果，可以得出如下结论：由于传感器失效会引起量测信号丢失，所以在对分数阶网络系统进行状态估计器设计时必须计及量测信号丢失的情形。本文档来自技高网...

【技术保护点】
一种计及传感器故障的分数阶网络系统状态估计方法，其特征在于，包括如下步骤：1)初始化k时刻的估计初始值

【技术特征摘要】
1.一种计及传感器故障的分数阶网络系统状态估计方法，其特征在于，包括如下步骤：1)初始化k时刻的估计初始值和估计误差协方差Pk，估计时刻最大值N；2)计算k时刻的系统函数的雅克比矩阵，计算公式如下

【专利技术属性】
技术研发人员：孙永辉，王义，张博文，卫志农，孙国强，翟苏巍，汪婧，
申请(专利权)人：河海大学，
类型：发明
国别省市：江苏,32