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

本发明属于控制系统的网络状态估计技术领域，具体涉及一种随机通讯协议下复杂网络的状态估计方法。

背景技术

复杂网络的状态估计问题一直是控制系统中一种重要的研究问题，在工程、电网、社交网络等领域的信号估计任务中获得广泛应用。

对于随机通讯协议下具有随机内耦合以及测量丢失现象的非线性复杂网络，由于现有的状态估计方法不能同时处理具有此类现象的复杂网络的状态估计问题，因此，这些现象一直影响着现有状态估计方法的状态估计性能。

技术实现要素：

本发明的目的是为解决现有的状态估计方法不能够同时处理具有随机内耦合以及测量丢失现象的复杂网络的状态估计的问题。

本发明为解决上述技术问题采取的技术方案是：

一种随机通讯协议下复杂网络的状态估计方法，该方法的具体步骤为：

步骤一、在随机通讯协议下，建立具有测量丢失以及随机内耦合现象的复杂网络动态模型；

步骤二、对于步骤一建立的复杂网络动态模型进行增广，并构造状态估计器；利用状态估计器对增广后复杂网络动态模型进行状态估计；

步骤三、计算增广后复杂网络动态模型状态估计的一步预测误差协方差矩阵上界∑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)。状态变量的形式根据网络的输入确定。

具体实施方式三：本实施方式对实施方式二所述的一种随机通讯协议下复杂网络的状态估计方法进行进一步的限定，本实施方式中的所述步骤二的具体过程为：

对于复杂网络动态模型的每个节点i进行增广得到公式(5)：

其中：xk为增广后复杂网络动态模型在第k时刻的状态变量；x1，k为第1个节点在第k时刻的状态变量，为x1，k的转置，xk+1为增广后复杂网络动态模型在第k+1时刻的状态变量；f(xk)为增广后复杂网络动态模型在第k时刻的状态变量对应的非线性函数，f1(x1，k)是第1个节点在第k时刻的状态变量对应的非线性函数，为f1(x1，k)的转置；w＝[wij]n×n是耦合矩阵，代表w和γ的克罗内克(kronecker)积；代表w和的克罗内克(kronecker)积；

yk为增广后的复杂网络在第k时刻的测量输出，

ξk代表随机变量ξi，k增广后的矩阵，diag{ξ1，k，ξ2，k，…，ξn，k}代表由{ξ1，k，ξ2，k，…，ξn，k}构成的对角矩阵；

bk代表由bi，k组成的对角矩阵，bk＝diag{b1，k，b2，k，…，bn，k}；

代表增广后复杂网络动态模型在第k时刻的过程噪声；

ck代表由ci，k组成的对角矩阵，ck＝diag{c1，k，c2，k，…，cn，k}；

vk代表增广后复杂网络动态模型在第k时刻的测量噪声；

λk为由服从伯努利分布的随机变量λi，k增广后的矩阵，；ip为p维的单位矩阵，diag{λ1，k，λ2，k，…，λn，k}代表由{λ1，k，λ2，k，…，λn，k}构成的对角矩阵，代表克罗内克(kronecker)积；{*}^t代表元素{*}的转置；

利用公式(5)式构造状态估计器：

式中：是增广后复杂网络动态模型在k时刻的一步预测，为k时刻对增广后复杂网络动态模型状态的估计，为节点1在k时刻对增广后复杂网络动态模型状态的估计，是的转置矩阵；为对应的非线性函数，为的转置；为k+1时刻对增广后复杂网络动态模型状态的估计，为增广后复杂网络动态模型在k+1时刻的实际测量输出，代表协议判别矩阵，且γk+1代表k+1时刻选中的通过复杂网络传输的节点；为期望概率矩阵，且ck+1代表由ci，k+1组成的对角矩阵，ci，k+1为第i个节点在第k+1时刻的测量矩阵，kk+1为增广后复杂网络动态模型在k+1时刻的估计增益矩阵，为增广后的复杂网络在第k+1时刻的实际测量输出，

具体实施方式四：本实施方式对实施方式三所述的一种随机通讯协议下复杂网络的状态估计方法进行进一步的限定，本实施方式中的步骤三的具体过程为：

按照下式计算出增广后复杂网络动态模型状态估计的一步预测误差协方差矩阵上界∑k+1|k：

式中，ak为f(xk)在点的泰勒展开系数矩阵，且∑k|k为增广后复杂网络动态模型在k时刻的估计误差协方差矩阵的上界，为∑k|k取逆；hk为f(xk)在点泰勒展开后的误差矩阵，ηk为大于零的时变实数且满足i是单位矩阵，为ηk的逆；lk为f(xk)在点泰勒展高阶后的无穷小误差矩阵，qk是增广后的过程噪声协方差矩阵(也是由过程噪声的方差为qi，k组成的对角矩阵，qk＝diag{q1，k，q2，k，…，qn，k})，是bk的转置矩阵；为过渡矩阵，且1n×n为n×n维元素全都为1的矩阵，in为n阶单位矩阵；为中间变量，且ε1为大于零的已知权重，是哈达玛积(hadamardproduct)。

具体实施方式五：本实施方式对实施方式四所述的一种随机通讯协议下复杂网络的状态估计方法进行进一步的限定，本实施方式中的步骤四的具体过程为：

按如下公式计算增广后复杂网络动态模型在k+1时刻的估计增益矩阵kk+1：

式中：ε2为大于零的常数，为ck+1的转置矩阵，为的转置矩阵，代表概率判别矩阵，φi＝diag{0，0，…，ip，…，0}(对角线第i块为p维单位矩阵，其余块全为p维的全0矩阵)代表第i个节点的判别矩阵，diag{·}为由{·}构成的对角矩阵，пk+1为中间变量矩阵，且其中：

为n×p维的单位矩阵，是随机变量λi，k的协方差矩阵的增广形式，是中间变量，且ε3为大于零的常数，υk+1代表增广后复杂网络动态模型在k+1时刻的未知概率的上界矩阵，rk+1是增广后复杂网络动态模型的测量噪声在k+1时刻的方差；tr{*}表示矩阵{*}的迹，为中间变量矩阵。

qi，k代表未知概率的上界，min为取小函数。

具体实施方式六：本实施方式对实施方式五所述的一种随机通讯协议下复杂网络的状态估计方法进行进一步的限定，本实施方式中步骤六根据步骤四中计算出的增广后复杂网络的k+1时刻的估计增益矩阵kk+1，计算出增广后复杂网络在k+1时刻的估计误差协方差矩阵上界∑k+1|k+1的具体过程为：

采用如下公式计算增广后复杂网络动态模型在k+1时刻的估计误差协方差矩阵上界∑k+1|k+1：

其中：为中间变量，为的转置矩阵；为φi的转置矩阵，为kk+1的转置矩阵；∑k+1k+1为增广后复杂网络在k+1时刻的估计误差协方差矩阵上界。

此外，可以证明在每个节点的测量丢失概率都确定且都一致的情况下，即随着测量丢失概率的增加，估计误差协方差矩阵上界∑k+1|k+1的迹是不增的。

实施例

采用本发明所述方法进行仿真：

系统参数：

b1，k＝[0.110.14]^t，b2，k＝[0.18+0.1sin(2k)0.12]^t，b3，k＝[0.140.12]^t，

c1，k＝[0.650.75]，c2，k＝[0.80.75]，c3，k＝[0.60.65+0.1cos(7k)]，γ＝diag{0.9，0.9}，

其它仿真初始值选取如下：

(i＝1，2，3)，q1，k＝0.15，q2，k＝0.25，q3，k＝0.1，r1，k＝0.15，r2，k＝0.15，r3，k＝0.2，ε1＝1，ε2＝0.1，ε3＝1，∑0|0＝10i6，p1＝0.3，p2＝0.3，p3＝0.4，l＝0.01i6，h＝0.01i6。

状态估计器效果：

如图2为第一个节点的第一个变量在k时刻的状态轨迹及其估计的对比图；如图3为第一个节点的第二个变量在k时刻的状态轨迹及其估计的对比图；如图4为第二个节点的第一个变量在k时刻的状态轨迹及其估计的对比图；如图5为第二个节点的第二个变量在k时刻的状态轨迹及其估计的对比图；如图6为第三个节点的第一个变量在k时刻的状态轨迹及其估计的对比图；如图7为第三个节点的第二个变量在k时刻的状态轨迹及其估计的对比图；如图8为本发明的γk值随时刻的变化曲线图；如图9为测量丢失概率分别为0.25、0.55、0.85和1下的，估计误差协方差矩阵上界σk+1k+1的迹随时刻变化的对比图；

由图2-图7可见，针对随机通讯协议下具有测量丢失现象的非线性复杂网络，所发明的状态估计器设计方法可有效地估计出目标状态。

由图9可以看出，随着概率的增加，估计误差协方差矩阵上界∑k+1|k+1的迹是不增的。