一种随机通讯协议下复杂网络的状态估计方法技术

一种随机通讯协议下复杂网络的状态估计方法，它用于控制系统的网络状态估计技术领域。本发明专利技术解决了现有的状态估计方法不能够同时处理具有随机内耦合以及测量丢失现象的复杂网络的状态估计的问题。本发明专利技术同时考虑了具有未知概率的测量丢失现象以及随机内耦合对状态估计性能的影响，利用扩展卡尔曼滤波方法全面考虑了估计误差协方差矩阵的有效信息，达到了抗扰动的目的；与现有的通讯协议下复杂网络状态估计方法相比，本发明专利技术的方法可以将估计误差控制在极小的范围内，在易于求解的同时，可以将估计的精确度提高10％以上。本发明专利技术可以应用于网络状态估计技术领域用。

【技术实现步骤摘要】
一种随机通讯协议下复杂网络的状态估计方法
本专利技术属于控制系统的网络状态估计
，具体涉及一种随机通讯协议下复杂网络的状态估计方法。
技术介绍
复杂网络的状态估计问题一直是控制系统中一种重要的研究问题，在工程、电网、社交网络等领域的信号估计任务中获得广泛应用。对于随机通讯协议下具有随机内耦合以及测量丢失现象的非线性复杂网络，由于现有的状态估计方法不能同时处理具有此类现象的复杂网络的状态估计问题，因此，这些现象一直影响着现有状态估计方法的状态估计性能。
技术实现思路
本专利技术的目的是为解决现有的状态估计方法不能够同时处理具有随机内耦合以及测量丢失现象的复杂网络的状态估计的问题。本专利技术为解决上述技术问题采取的技术方案是：一种随机通讯协议下复杂网络的状态估计方法，该方法的具体步骤为：步骤一、在随机通讯协议下，建立具有测量丢失以及随机内耦合现象的复杂网络动态模型；步骤二、对于步骤一建立的复杂网络动态模型进行增广，并构造状态估计器；利用状态估计器对增广后复杂网络动态模型进行状态估计；步骤三、计算增广后复杂网络动态模型状态估计的一步预测误差协方差矩阵上界∑k+1|k；步骤四、根据步骤三中计算出的增广后复杂网络动态模型状态估计的一步预测误差协方差矩阵上界∑k+1|k，计算出增广后复杂网络动态模型在k+1时刻的估计增益矩阵Kk+1；步骤五、将步骤四中获得的增广后复杂网络动态模型在k+1时刻的估计增益矩阵Kk+1代入步骤二中的构造状态估计器，得到增广后复杂网络动态模型在k+1时刻的状态估计判断k+1是否达到复杂网络总时长M，若k+1<M，则执行步骤六，若k+1＝M，则结束对复杂网络的状态估计；步骤六、根据步骤四中计算出的增广后复杂网络动态模型在k+1时刻的估计增益矩阵Kk+1，计算出增广后复杂网络动态模型在k+1时刻的估计误差协方差矩阵上界∑k+1|k+1；令k＝k+1，执行步骤三，直至满足k+1＝M。本专利技术的有益效果是：本专利技术提供了一种随机通讯协议下复杂网络的状态估计方法，本专利技术的随机通讯协议下复杂网络的状态估计方法同时考虑了具有未知概率的测量丢失现象以及随机内耦合对状态估计性能的影响，利用扩展卡尔曼滤波方法全面考虑了估计误差协方差矩阵的有效信息，达到了抗扰动的目的；与现有的通讯协议下复杂网络状态估计方法相比，本专利技术的方法可以将估计误差控制在极小的范围内，在易于求解的同时，可以将估计的精确度提高10％以上。附图说明图1为本专利技术所述的一种随机通讯协议下复杂网络的状态估计方法的流程图；图2为本专利技术第一个节点的第一个变量在k时刻的状态轨迹及其估计的对比图；图3为本专利技术第一个节点的第二个变量在k时刻的状态轨迹及其估计的对比图；图4为本专利技术第二个节点的第一个变量在k时刻的状态轨迹及其估计的对比图；图5为本专利技术第二个节点的第二个变量在k时刻的状态轨迹及其估计的对比图；图6为本专利技术第三个节点的第一个变量在k时刻的状态轨迹及其估计的对比图；图7为本专利技术第三个节点的第二个变量在k时刻的状态轨迹及其估计的对比图；图8为本专利技术的γk值随时刻的变化曲线图；图9为本专利技术的测量丢失概率分别为0.25、0.55、0.85和1下的，估计误差协方差矩阵上界的迹tr(∑k+1|k+1)随时刻变化的对比图；具体实施方式下面结合附图对本专利技术的技术方案作进一步的说明，但并不局限于此，凡是对本专利技术技术方案进行修改或者等同替换，而不脱离本专利技术技术方案的精神和范围，均应涵盖在本专利技术的保护范围中。具体实施方式一：结合图1说明本实施方式。本实施方式所述的一种随机通讯协议下复杂网络的状态估计方法，该方法的具体步骤为：步骤一、在随机通讯协议下，建立具有测量丢失以及随机内耦合现象的复杂网络动态模型；步骤二、对于步骤一建立的复杂网络动态模型进行增广，并构造状态估计器；利用状态估计器对增广后复杂网络动态模型进行状态估计；步骤三、计算增广后复杂网络动态模型状态估计的一步预测误差协方差矩阵上界∑k+1|k；步骤四、根据步骤三中计算出的增广后复杂网络动态模型状态估计的一步预测误差协方差矩阵上界∑k+1|k，计算出增广后复杂网络动态模型在k+1时刻的估计增益矩阵Kk+1；步骤五、将步骤四中获得的增广后复杂网络动态模型在k+1时刻的估计增益矩阵Kk+1代入步骤二中的构造状态估计器，得到增广后复杂网络动态模型在k+1时刻的状态估计判断k+1是否达到复杂网络总时长M，若k+1<M，则执行步骤六，若k+1＝M，则结束对复杂网络的状态估计；步骤六、根据步骤四中计算出的增广后复杂网络动态模型在k+1时刻的估计增益矩阵Kk+1，计算出增广后复杂网络动态模型在k+1时刻的估计误差协方差矩阵上界∑k+1|k+1；令k＝k+1，执行步骤三，直至满足k+1＝M。所述复杂网络可以为卫星构成的网络、机器人构成的网络、航天器构成的网络或雷达构成的网络。具体实施方式二：本实施方式对实施方式一所述的一种随机通讯协议下复杂网络的状态估计方法进行进一步的限定，所述步骤一的具体过程为：在随机通讯协议下，建立具有测量丢失以及随机内耦合现象的复杂网络动态模型，所述复杂网络动态模型的状态空间形式为：yi，k＝λi，kCi，kxi，k+vi，k(2)其中：xi，k代表第i个节点在第k时刻的状态变量，i＝1，2，…，N，N代表节点的个数，；xi，k+1代表第i个节点在第k+1时刻的状态变量，fi(xi，k)是第i个节点在第k时刻的状态变量对应的非线性函数，且非线性函数fi(xi，k)对xi，k是可微的，j表示与i耦合的节点，j＝1，2，…，N；wij代表第i个节点与第j个节点的耦合强度；代表n维欧氏空间，Γ代表内耦合矩阵，代表内耦合不确定矩阵；xj，k代表第j个节点在第k时刻的状态变量；yi，k为第i个节点在第k时刻的测量输出，且代表P维欧氏空间；是第i个节点在第k时刻的均值为零、方差为Qi，k的过程噪声，且代表q维欧氏空间，是第i个节点在第k时刻的均值为零、方差为Ri，k的测量噪声；Bi，k为第i个节点在第k时刻的噪声分布矩阵，Ci，k为第i个节点在第k时刻的测量矩阵；ξi，k为服从均值为零的单位方差的高斯白噪声，用来刻画随机内耦合现象，λi，k为服从伯努利分布的随机变量，且随机变量λi，k＝1时的概率Prob{λi，k＝1}和随机变量λi，k＝0时的概率Prob{λi，k＝0}分别表示如下：其中，为第i个节点在k时刻未发生测量丢失现象的期望概率，是已知的常数，Δλi，k代表未知概率，满足|Δλi，k|≤Qi，k刻画概率的未知性，是未知概率的上界；γk代表k时刻选中的通过复杂网络传输的节点，则经过随机通讯协议后，第i个节点在第k时刻的实际测量输出为：定义中间变量则将(4)式表示为其中，δ(γk-i)是δ函数，δ(γk-i)定义为：γk＝i的概率Prob{γk＝i}为pi，pi表示第i个节点被选中传输到复杂网络的概率值。为了避免数据传输过程中由于网络拥堵造成的数据冲突现象，引入随机通讯协议(RAP)。状态变量的形式根据网络的输入确定。具体实施方式三：本实施方式对实施方式二所述的一种随机通讯协议下复杂网络的状态估计方法进行进一步的限定，本实施方式中的所述步骤二的具体过程为：对于复杂网络动态模型的每本文档来自技高网...

【技术保护点】
1.一种随机通讯协议下复杂网络的状态估计方法，其特征在于，该方法的具体步骤为：步骤一、在随机通讯协议下，建立具有测量丢失以及随机内耦合现象的复杂网络动态模型；步骤二、对于步骤一建立的复杂网络动态模型进行增广，并构造状态估计器；利用状态估计器对增广后复杂网络动态模型进行状态估计；步骤三、计算增广后复杂网络动态模型状态估计的一步预测误差协方差矩阵上界∑k+1|k；步骤四、根据步骤三中计算出的增广后复杂网络动态模型状态估计的一步预测误差协方差矩阵上界∑k+1|k，计算出增广后复杂网络动态模型在k+1时刻的估计增益矩阵Kk+1；步骤五、将步骤四中获得的增广后复杂网络动态模型在k+1时刻的估计增益矩阵Kk+1代入步骤二中的构造状态估计器，得到增广后复杂网络动态模型在k+1时刻的状态估计

【技术特征摘要】
1.一种随机通讯协议下复杂网络的状态估计方法，其特征在于，该方法的具体步骤为：步骤一、在随机通讯协议下，建立具有测量丢失以及随机内耦合现象的复杂网络动态模型；步骤二、对于步骤一建立的复杂网络动态模型进行增广，并构造状态估计器；利用状态估计器对增广后复杂网络动态模型进行状态估计；步骤三、计算增广后复杂网络动态模型状态估计的一步预测误差协方差矩阵上界∑k+1|k；步骤四、根据步骤三中计算出的增广后复杂网络动态模型状态估计的一步预测误差协方差矩阵上界∑k+1|k，计算出增广后复杂网络动态模型在k+1时刻的估计增益矩阵Kk+1；步骤五、将步骤四中获得的增广后复杂网络动态模型在k+1时刻的估计增益矩阵Kk+1代入步骤二中的构造状态估计器，得到增广后复杂网络动态模型在k+1时刻的状态估计判断k+1是否达到复杂网络总时长M，若k+1<M，则执行步骤六，若k+1＝M，则结束对复杂网络的状态估计；步骤六、根据步骤四中计算出的增广后复杂网络动态模型在k+1时刻的估计增益矩阵Kk+1，计算出增广后复杂网络动态模型在k+1时刻的估计误差协方差矩阵上界∑k+1|k+1；令k＝k+1，执行步骤三，直至满足k+1＝M。2.根据权利要求1所述的一种随机通讯协议下复杂网络的状态估计方法，其特征在于，所述步骤一的具体过程为：在随机通讯协议下，建立具有测量丢失以及随机内耦合现象的复杂网络动态模型，所述复杂网络动态模型的状态空间形式为：yi,k＝λi,kCi,kxi,k+νi,k(2)其中：xi,k代表第i个节点在第k时刻的状态变量，i＝1,2,…,N，N代表节点的个数，；xi,k+1代表第i个节点在第k+1时刻的状态变量，fi(xi,k)是第i个节点在第k时刻的状态变量对应的非线性函数，j表示与i耦合的节点，j＝1,2,…,N；wij代表第i个节点与第j个节点的耦合强度；代表n维欧氏空间，Γ代表内耦合矩阵，代表内耦合不确定矩阵；xj,k代表第j个节点在第k时刻的状态变量；yi,k为第i个节点在第k时刻的测量输出；是第i个节点在第k时刻的均值为零的过程噪声，且代表q维欧氏空间，νi,k是第i个节点在第k时刻的均值为零的测量噪声；Bi,k为第i个节点在第k时刻的噪声分布矩阵，Ci,k为第i个节点在第k时刻的测量矩阵；ξi,k为服从均值为零的单位方差的高斯白噪声，λi,k为服从伯努利分布的随机变量，且随机变量λi,k＝1时的概率Prob{λi,k＝1}和随机变量λi,k＝0时的概率Prob{λi,k＝0}分别表示如下：其中，为已知的常数，Δλi,k代表未知概率；γk代表k时刻选中的通过复杂网络传输的节点，则经过随机通讯协议后，第i个节点在第k时刻的实际测量输出为：定义中间变量则将(4)式表示为其中，δ(γk-i)是δ函数，δ(γk-i)定义为：γk＝i的概率Prob{γk＝i}为pi，pi表示第i个节点被选中传输到复杂网络的概率值。3.根据权利要求2所述的一种随机通讯协议下复杂网络的状态估计方法，其特征在于，所述步骤二的具体过程为：对于复杂网络动态模型的每个节点i进行增广得到公式(5)：其中：xk为增广后复杂网络动态模型在第k时刻的状态变量；x1,k为第1个节点在第k时刻的状态变量，为x1,k的转置，xk+1为增广后复杂网络动态模型在第k+1时刻的状态变量；f(xk)为增广...

【专利技术属性】
技术研发人员：胡军，张红旭，赵文杰，杨昱，宋诗宇，高岩，
申请(专利权)人：哈尔滨理工大学，
类型：发明
国别省市：黑龙江,23