具有切换拓扑结构的复杂网络的状态估计方法技术

本发明专利技术公开了具有切换拓扑结构的复杂网络的状态估计方法。本发明专利技术方法首先基于二进制编码方案，建立具有时滞的离散非线性复杂网络模型，对其进行状态估计，得到状态估计误差，进而得到状态估计增广系统；根据李亚普诺夫稳定性定理，得到状态估计器增益矩阵；将状态估计器增益矩阵带入状态估计器模型，完成对考虑切换拓扑结构的时滞非线性复杂网络系统基于二进制编码的状态估计。本发明专利技术方法考虑了具有耦合性、时滞以及非线性的复杂网络，能够用来建模许多工程中的复杂网络。本发明专利技术方法能够减少传输过程能量的消耗与网络带宽资源的占用，对保障网络安全平稳运行有着重要应用。保障网络安全平稳运行有着重要应用。保障网络安全平稳运行有着重要应用。

【技术实现步骤摘要】
具有切换拓扑结构的复杂网络的状态估计方法


[0001]本专利技术属于计算机
，特别是状态估计、控制方法
，具体涉及一种具有切换拓扑结构的时滞非线性复杂网络基于二进制编码的状态估计方法。

技术介绍

[0002]复杂网络由于能够模拟社交网络、电网、通信网络等各种大型复杂系统，而受到了广大学者的关注，而现实生活中的电网和通信网络的动态变化使分析动态问题带来了巨大挑战。众所周知，状态信息在复杂系统的动态分析中起着至关重要的作用。然而，在大多数情况下，各节点的状态信息是未知的，所以状态估计已成为一种热门的研究方法，通过已知可测量的信息来进行状态估计。由于通信网络的飞速发展，二进制编码方案成为数字通信中最流行的方案之一由于其传输的鲁棒性和易实现性，在二进制编码方案下，信号首先应该被编码成一组二进制位，通过二进制对称信道进行传输。由于信道噪声的影响，此传输过程不可避免地会产生随机比特误差，使得解码得到的测量信号明显偏离原始信号。近年来，人们对于二进制编码方案进行了初步的研究，然而在二进制编码下，还未有针对于具有切换拓扑结构的复杂网络状态估计方法。

技术实现思路

[0003]本专利技术的目的在于提供一种具有切换拓扑结构的复杂网络的状态估计方法，以解决目前现有的状态估计方法不能同时考虑马尔科夫参数、传感器非线性、时滞、随机干扰，进而影响状态估计性能的问题。
[0004]本专利技术方法具体是：
[0005]步骤(1)基于二进制编码方案，建立具有时滞的离散非线性复杂网络模型；
[0006]步骤(2)基于二进制编码方案，对具有时滞的离散非线性复杂网络模型进行状态估计，得到状态估计误差；
[0007]步骤(3)根据状态估计误差，得到状态估计增广系统；
[0008]步骤(4)利用状态估计增广系统，根据李亚普诺夫稳定性定理，得到状态估计器增益矩阵；
[0009]步骤(5)将状态估计器增益矩阵带入步骤(2)中的状态估计器模型，完成对考虑切换拓扑结构的时滞非线性复杂网络系统基于二进制编码的状态估计。
[0010]进一步，步骤(1)具体是：
[0011]所述的离散非线性复杂网络模型状态空间形式为：
[0012][0013][0014]其中，表示复杂网络动态模型中第i个节点在时刻的状态向量，表示第i个节点在时刻的状态向量,表示第i个节点在时刻从系统中接收到的测量
输出；f(
·
)表示满足扇形有界条件的非线性激励函数；和表示期望为零、且满足方差有界的噪声；表示复杂网络动态系统的时滞满足上下界为正标量τ
m
和τ
M
的范围，表示第j个节点在时刻的状态向量；在时刻的模态为随机变量，是一个齐次马尔可夫链在有限的状态空间范围{1,2,
…
,S}，模态序号s＝1,2,
…
,S，S为模态数量，Γ表示内耦合矩阵是一个对角阵，中间函数δ(
·
)定义为)定义为表示复杂网络动态模型中第s种模态下第i个节点和第j个节点的加权邻接参数，i,j＝1,2,
…
,N，且i≠j，N为节点数量，加权邻接[
·
]N
×
N
表示N
×
N维矩阵；A
i
、C
i
、D
i
为已知的参数矩阵；
[0015]测量输出首先被编码到一个有限长度二进制字符串中，然后通过一个无记忆的二进制对称信道传输到估计器进行处理，具体过程如下：
[0016]首先测量输出通过编码器被编码成一个长度为L的二进制字符串，将的振幅范围[
‑
h,h]范围划分为2
L
‑
1段，且均匀间隔长度
[0017]使用随机截断函数对测量输出进行预处理，得到截断信号具体是：当J
l
表示设定的第l个二进制字符的数值，l＝1,2,
…
,L，得到：的概率的概率的概率概率值截断误差是一个随机变量，的期望的期望的方差
[0018]将截断信号编码成一个长度为L的二进制位字符串且其中的字符其中的字符通过得到。
[0019]二进制位字符串通过无记忆的二进制对称信道进行传输，以交叉概率进行翻转，翻转后字符串且其中的字符且其中的字符且其中的字符为伯努利变量，满足为伯努利变量，满足交叉概率α∈[0,1]。
[0020]根据接收到的翻转后字符串解码后恢复的信号
[0021]的期望方差中间值
[0022]为了补偿失真，进行恢复测量，恢复后测量输出
[0023]设定比特反转误差为一个随机变量，其期望方差
[0024]最后得到的离散非线性复杂网络模型的测量输出
[0025]再进一步，步骤(2)中状态估计器模型为：
[0026][0027]其中，表示状态估计器模型中第i个节点在时刻对的估计值，表示第i个节点在时刻对的估计值，表示与模态相关的待求状态估计器增益矩阵。
[0028]状态估计误差方程为：
[0029]其中，K
i,s
表示与第s个模态相关的待求状态估计器增益矩阵。
[0030]第i个节点在时刻的状态估计误差第i个节点在时刻的状态估计误差
[0031]更进一步，步骤(3)中状态估计增广系统表示为：
[0032][0033]其中，时刻增广后状态估计误差col
N
{
·
}表示N维列向量，增广后待求状态估计器增益矩阵K
s
＝diag
N
{K
i,s
}，diag
N
{
·
}表示N维对角阵；增广后的参数矩阵A＝col
N
{A
i
},D＝col
N
{D
i
},C＝diag
N
{C
i
}，增广后的截断误差增广后的比特反转误差增广后的加权邻接矩阵中间函数函数表示克罗内克积运算。
[0034]又进一步，步骤(4)具体是：
[0035]由得到矩阵P
S
,X
S
，通过公式求解出状态估计器增益矩阵K
S
；其中，P
S
,X
S
为待求矩阵，I为单位矩阵，γ，π均为未知的正常数，T表示转置操作，*表示矩阵中的对称项，非线性系数矩阵非线性参数矩阵Q表示未知对称正定矩阵，U1，U2为已知的常数矩阵，待求的变异矩阵转移概率值满足0＜λ
s
′
s
＜1且
[0036]所述的李亚普诺夫稳定性定理为：
[0037][0038]其中：
[0039]函数第一分量函数第一分量表示与相关的未知对称正定矩阵；
[0040]函数第二分量t表示t时刻；
[0041]函数第三分量t
′...

【技术保护点】

【技术特征摘要】
1.具有切换拓扑结构的复杂网络的状态估计方法，其特征在于：步骤(1)基于二进制编码方案，建立具有时滞的离散非线性复杂网络模型；步骤(2)基于二进制编码方案，对具有时滞的离散非线性复杂网络模型进行状态估计，得到状态估计误差；步骤(3)根据状态估计误差，得到状态估计增广系统；步骤(4)利用状态估计增广系统，根据李亚普诺夫稳定性定理，得到状态估计器增益矩阵；步骤(5)将状态估计器增益矩阵带入步骤(2)中的状态估计器模型，完成对考虑切换拓扑结构的时滞非线性复杂网络系统基于二进制编码的状态估计。2.如权利要求1所述的具有切换拓扑结构的复杂网络的状态估计方法，其特征在于：步骤(1)中所述的离散非线性复杂网络模型的状态空间形式为：骤(1)中所述的离散非线性复杂网络模型的状态空间形式为：其中，表示复杂网络动态模型中第i个节点在时刻的状态向量，表示第i个节点在时刻的状态向量,表示第i个节点在时刻从系统中接收到的测量输出；f(
·
)表示满足扇形有界条件的非线性激励函数；和表示期望为零、且满足方差有界的噪声；表示复杂网络动态系统的时滞满足上下界为正标量τ
m
和τ
M
的范围，表示第j个节点在时刻的状态向量；在时刻的模态为随机变量，是一个齐次马尔可夫链在有限的状态空间范围{1,2,
…
,S}，模态序号s＝1,2,
…
,S，S为模态数量，Γ表示内耦合矩阵是一个对角阵，中间函数δ(
·
)定义为)定义为表示复杂网络动态模型中第s种模态下第i个节点和第j个节点的加权邻接参数，i,j＝1,2,
…
,N，且i≠j，N为节点数量，加权邻接[
·
]
N
×
N
表示N
×
N维矩阵；A
i
、C
i
、D
i
为已知的参数矩阵；测量输出首先被编码到一个有限长度二进制字符串中，然后通过一个无记忆的二进制对称信道传输到估计器进行处理，具体过程如下：首先测量输出通过编码器被编码成一个长度为L的二进制字符串，将的振幅范围[
‑
h,h]范围划分为2
L
‑
1段，且均匀间隔长度使用随机截断函数对测量输出进行预处理，得到截断信号具体是：当J
l
表示设定的第l个二进制字符的数值，l＝1,2,
…
,L，得到：的概率的概率概率值截断误差是一个随机变量，的期望
的方差将截断信号编码成一个长度为L的二进制位字符串且其中的字符通过得到；二进制位...

【专利技术属性】
技术研发人员：李琦，陈凯歌，谭海龙，盛伟国，
申请(专利权)人：杭州师范大学，
类型：发明
国别省市：