杂志的下一篇文章
变形过程的分数模型
下一篇特刊文章
随机延迟不确定CE151直升机系统基于采样数据的容错控制设计:Takagi-Sugeno模糊方法
期刊上的上一篇文章
多串分数阶记忆Hopfield神经网络的动态分析及其在物联网中的音频加密应用
 
 
订购文章重印
字体类型:
宋体 佐治亚州 宋体,Verdana
字体大小:
澳大利亚 澳大利亚 澳大利亚
行距:
列宽:
背景:
第条

用图论方法识别多分散随机多组模型的局部拓扑结构

1
西南交通大学数学学院,中国成都610031
2
日本名古屋丰田技术研究所计算智能实验室,邮编:468-8511
*
应向其发送信件的作者。
收到的提交文件:2022年5月23日/修订日期:2022年6月27日/接受日期:2022年6月28日/发布日期:2022年7月1日

摘要

:
本文研究了具有多重离散的随机多群模型的局部拓扑识别问题。基于自适应钉扎控制和图论方法,得到了具有多重离散的随机多群模型局部拓扑识别的一些充分判据。也就是说,可以成功地识别出未知的局部拓扑结构。最后,通过数值算例验证了理论结果的有效性。

1.简介

近年来,多群模型因其在生物学、流行病学等许多不同领域的广泛应用而受到越来越多的关注[1,2]. 多组模型的不同动力学行为已被广泛研究,参见[,4]为了全球稳定[5]用于同步和[6]用于平稳分布。
在许多多种群模型中,特别是在多径环境中,多重扩散具有很大的影响,并且物种之间总是存在不同的扩散。此外,自然界中的系统必须受到随机扰动的影响[7,8,9,10]. 具有多重分散性的随机多群模型是有效的数学模型,已引起越来越多的关注,参见[11,12,13]. 应该注意的是,许多研究中的拓扑结构[,4,5,6,8,9,10,11,12,13]已知。事实上,在许多实际应用中,拓扑结构通常是未知的或不确定的。因此,识别具有多重扩散的随机多群模型的未知拓扑结构非常重要。
在本文中,我们考虑以下具有多重分散的随机多群模型
d日 x个 k个 ( ) ( t吨 ) = ϕ k个 ( ) ( x个 k个 ( t吨 ) , t吨 ) + 小时 = 1 N个 k个 小时 ( ) H(H) k个 小时 ( ) ( x个 k个 ( ) ( t吨 ) , x个 小时 ( ) ( t吨 ) , t吨 ) d日 t吨 + ψ k个 ( ) ( x个 k个 ( ) ( t吨 ) , t吨 ) d日 W公司 ( t吨 ) , 1 , 1 k个 N个 ,
哪里 x个 k个 ( ) R(右) 是的状态向量-中的第个组件k个-第个组。我们定义 = = 1 对于 R(右) + . x个 k个 ( t吨 ) = ( ( x个 k个 ( 1 ) ( t吨 ) ) T型 , ( x个 k个 ( 2 ) ( t吨 ) ) T型 , , ( x个 k个 ( ) ( t吨 ) ) T型 ) T型 R(右) 表示的状态向量k个-第th组。 ϕ k个 ( ) ( x个 k个 ( t吨 ) , t吨 ) : R(右) × R(右) + R(右) 代表的性能-的第个组件k个-第个组。 ψ k个 ( ) ( x个 k个 ( ) ( t吨 ) , t吨 ) : R(右) × R(右) + R(右) 显示了上的扰动强度-的第个组件k个-第个组。 k个 小时 ( ) 是指-第个组件小时-第个组到k个-第个组。在这里, k个 小时 ( ) = 0 如果没有分散-第个组件小时-第个组到k个-第个组。 H(H) k个 小时 ( ) : R(右) × R(右) × R(右) + R(右) 是顶点的影响小时在顶点上k个, W公司 ( · ) 是一维布朗运动。为了更好地理解,我们画了一个由四个群和两个分散体组成的图,参见图1在这里, 11 ( 1 ) = 13 ( 1 ) = 14 ( 1 ) = 22 ( 1 ) = 24 ( 1 ) = 31 ( 1 ) = 33 ( 1 ) = 34 ( 1 ) = 44 ( 1 ) = 0 , 11 ( 2 ) = 14 ( 2 ) = 21 ( 2 ) = 22 ( 2 ) = 24 ( 2 ) = 33 ( 2 ) = 41 ( 2 ) = 42 ( 2 ) = 44 ( 2 ) = 0 .
在本文中,数学模型有若干组。如果我们为每个组添加一个控制器,将导致高控制成本。众所周知,固定控制策略是减少受控组数量的有效技术[14,15,16]. 另一方面,如果大量群具有未知的部分拓扑结构,或者我们只对部分拓扑结构感兴趣,那么钉扎控制是一种有效的技术。因此,我们将尝试使用固定控制机制来识别未知的部分拓扑结构(1).
值得一提的是,Lyapunov方法是研究局部拓扑识别的有效工具[17,18,19]. 然而,很难为(1)由于多组和分散。令人鼓舞的是,Li和他的合作者将图论和Lyapunov方法相结合,间接地为耦合系统构建了全局Lyapunow函数。同时,他们使用这种方法来研究耦合系统的全局稳定性[20]. 这种方法称为图形理论方法。据我们所知,利用图论方法对随机多群模型进行局部拓扑识别的工作很少。
基于上述讨论,本文旨在利用图理论方法识别具有多重离散的随机多群模型的未知局部拓扑结构。主要贡献如下。
  • 该数学模型是通用的,包括多重分散和随机扰动。
  • 钉扎控制器具有成本效益并且可以减少受控组。
  • 局部拓扑识别的图形理论方法是一种新颖的方法。
  • 成功地识别了具有多重扩散的随机多群模型的未知局部拓扑结构。
本文的其余部分组织如下。第2节,显示一些初步信息。主要结果介绍于第3节.英寸第4节,用两个数值算例验证了理论结果的有效性。结论见第5节.

2.前期工作

图论和随机微分方程的一些必要基础知识

有向图 M(M) = ( D类 , P(P) ) 包含一个集合 D类 = { 1 , 2 , , N个 } 顶点和集合P(P)弧的 ( , j个 ) ,其中 ( , j个 ) 是从初始顶点算起的弧到端子顶点j个.给定一个有向图 M(M) 具有N个顶点,我们定义加权矩阵 单位 = ( u个 j个 ) N个 × N个 谁的条目 u个 j个 等于弧的权重 ( j个 , ) 如果从顶点有弧到顶点j个,否则为0。带权矩阵的有向图单位描述为 ( M(M) , 单位 ) 拉普拉斯矩阵 属于 ( M(M) , 单位 ) 定义为 = ( j个 ) N个 × N个 ,其中 = j个 u个 j个 , j个 = u个 j个 ( j个 ) , 1 , j个 N个 .
假设系数 ϕ k个 ( ) ψ k个 ( ) 第页,共页(1)满足局部Lipschitz条件和线性增长条件[21,22]. 然后,对于初始值 x个 0 , (1)具有唯一的连续解,表示为 x个 ( t吨 ; x个 0 ) 此外,如果 ϕ k个 ( ) ( 0 , t吨 ) = 0 , ψ k个 ( ) ( 0 , t吨 ) = 0 H(H) k个 小时 ( ) ( 0 , 0 , t吨 ) = 0 ,然后(1)承认一个微不足道的解决方案 x个 ( t吨 ; 0 ) 0 .
我们定义微分算子 执行 V(V) k个 ( ) C类 2 , 1 ( R(右) × R(右) + ; R(右) + ) 以及(1)作为
V(V) k个 ( ) ( x个 k个 ( ) , t吨 ) = V(V) k个 ( ) ( x个 k个 ( ) , t吨 ) t吨 + V(V) k个 ( ) ( x个 k个 ( ) , t吨 ) x个 k个 ( ) ϕ k个 ( ) ( x个 k个 ( ) ( t吨 ) , t吨 ) + 小时 = 1 N个 k个 小时 ( ) H(H) k个 小时 ( ) ( x个 k个 ( ) ( t吨 ) , x个 小时 ( ) ( t吨 ) , t吨 ) + 1 2 跟踪 ψ k个 ( ) T型 ( x个 k个 ( ) ( t吨 ) , t吨 ) 2 V(V) k个 ( ) ( x个 k个 ( ) , t吨 ) x个 k个 ( ) 2 ψ k个 ( ) ( x个 k个 ( ) ( t吨 ) , t吨 ) ,
哪里
V(V) k个 ( ) ( x个 k个 ( ) , t吨 ) x个 k个 ( ) = V(V) k个 ( ) ( x个 k个 ( ) , t吨 ) x个 k个 ( 1 ) , V(V) k个 ( ) ( x个 k个 ( ) , t吨 ) x个 k个 ( 2 ) , , V(V) k个 ( ) ( x个 k个 ( ) , t吨 ) x个 k个 ( ) , 2 V(V) k个 ( ) ( x个 k个 ( ) , t吨 ) ( x个 k个 ( ) ) 2 = 2 V(V) k个 ( ) ( x个 k个 ( ) , t吨 ) x个 k个 ( 第页 ) x个 k个 ( q个 ) × .
柠檬 1
([21]).假设有一个函数 V(V) C类 2 , 1 ( R(右) N个 × R(右) + ; R(右) + ) ,一个函数 ϑ 1 ( R(右) + ; R(右) + ) 和一个连续函数 ξ : R(右) N个 R(右) + 这样的话
| x个 | inf公司 0 t吨 < V(V) ( x个 , t吨 ) = ,
和微分算子 作用于V的轨迹(1)满足
V(V) ( x个 , t吨 ) ϑ ( t吨 ) ξ ( x个 ) , ( x个 , t吨 ) R(右) N个 × R(右) + .
此外, ψ k个 被限定为 1 , 1 k个 N个 然后,对于每个初始值 x个 0 R(右) N个 , t吨 V(V) ( x个 ( t吨 ; x个 0 ) , t吨 ) 存在并且几乎肯定是有限的。此外,
t吨 ξ ( x个 ( t吨 ; x个 0 ) ) = 0 . .
柠檬 2
(中的定理2.2[20]).假设 N个 2 . ν k个 是第k个对角线元素的余因子 。然后,以下身份保持不变:
k个 , 小时 = 1 N个 ν k个 k个 小时 F类 k个 小时 ( x个 k个 , x个 小时 ) = ω ( ) ( , 第页 ) 电子 ( C类 ) F类 第页 ( x个 第页 , x个 ) .
在这里, F类 k个 小时 ( x个 k个 , x个 小时 ) , k个 , 小时 = 1 , 2 , , N个 是任意函数, 是的所有生成单圈图的集合 ( G公司 , A类 ) , ω ( ) 是的重量 、和 C类 表示的定向循环 .

3.主要成果

在本节中,我们将研究基于自适应钉扎同步和图形理论方法的具有多重分散的随机多组模型的局部拓扑识别问题。
为了识别具有多重分散的随机多群模型的局部拓扑结构,取(1)作为驱动系统,带有自适应固定控制器的响应系统描述为
d日 k个 ( ) ( t吨 ) = ϕ k个 ( ) ( k个 ( t吨 ) , t吨 ) + 小时 = 1 b条 k个 小时 ( ) H(H) k个 小时 ( ) ( k个 ( ) ( t吨 ) , 小时 ( ) ( t吨 ) , t吨 ) + 小时 = + 1 N个 b条 k个 小时 ( ) H(H) k个 小时 ( ) ( x个 k个 ( ) ( t吨 ) , x个 小时 ( ) ( t吨 ) , t吨 ) + u个 k个 ( ) ( t吨 ) d日 t吨 + ψ k个 ( ) ( k个 ( ) ( t吨 ) , t吨 ) d日 W公司 ( t吨 ) , 1 , 1 k个 .
ψ k个 ( ) ( k个 ( ) ( t吨 ) , t吨 ) ϕ k个 ( ) ( k个 ( t吨 ) , t吨 ) 扰动强度和函数控制-的第个组件k个-分别为响应网络中的第th组。 k个 ( ) = ( k个 1 ( ) , k个 2 ( ) , , k个 ( ) ) T型 R(右) 是的状态向量-中的第个组件k个-第个组。 k个 ( t吨 ) = ( ( k个 ( 1 ) ( t吨 ) ) T型 , ( k个 ( 2 ) ( t吨 ) ) T型 , , ( k个 ( ) ( t吨 ) ) T型 ) T型 R(右) 表示的状态向量k个-第个组。 B类 ( ) = ( b条 k个 小时 ( ) ) × N个 表示未知部分耦合矩阵的估计 A类 ( ) = k个 小时 ( ) × N个 ( 1 ) . u个 k个 ( ) 是通用的自适应固定控制器。在不损失一般性的情况下,识别由前沿组成的多分散随机多群模型的局部拓扑结构就足够了群体及其扩散,即, ( k个 小时 ( ) ) × N个 .让 e(电子) ( t吨 ) = ( t吨 ) x个 ( t吨 ) = ( e(电子) 1 T型 ( t吨 ) , e(电子) 2 T型 ( t吨 ) , , e(电子) T型 ( t吨 ) ) T型 是同步错误,其中 e(电子) k个 ( t吨 ) = k个 ( t吨 ) x个 k个 ( t吨 ) = ( ( e(电子) k个 ( 1 ) ( t吨 ) ) T型 , ( e(电子) k个 ( 2 ) ( t吨 ) ) T型 , , ( e(电子) k个 ( ) ( t吨 ) ) T型 ) T型 ( 1 k个 ) .我们表示 c(c) k个 小时 ( ) = b条 k个 小时 ( ) k个 小时 ( ) ( 1 , 1 k个 , 1 小时 N个 ) .系统间同步误差的动力系统(1)和()可以写为
d日 e(电子) k个 ( ) ( t吨 ) = [ ϕ k个 ( ) ( k个 ( t吨 ) , t吨 ) ϕ k个 ( ) ( x个 k个 ( t吨 ) , t吨 ) + 小时 = 1 b条 k个 小时 ( ) H(H) k个 小时 ( ) ( k个 ( ) ( t吨 ) , 小时 ( ) ( t吨 ) , t吨 ) + 小时 = + 1 N个 b条 k个 小时 ( ) H(H) k个 小时 ( ) ( x个 k个 ( ) ( t吨 ) , x个 小时 ( ) ( t吨 ) , t吨 ) 小时 = 1 N个 k个 小时 ( ) H(H) k个 小时 ( ) ( x个 k个 ( ) ( t吨 ) , x个 小时 ( ) ( t吨 ) , t吨 ) + u个 k个 ( ) ( t吨 ) ] d日 t吨 + [ ψ k个 ( ) ( k个 ( ) ( t吨 ) , t吨 ) ψ k个 ( ) ( x个 k个 ( ) ( t吨 ) , t吨 ) ] d日 W公司 ( t吨 ) = [ ϕ k个 ( ) ( k个 ( t吨 ) , t吨 ) ϕ k个 ( ) ( x个 k个 ( t吨 ) , t吨 ) + 小时 = 1 c(c) k个 小时 ( ) H(H) k个 小时 ( ) ( k个 ( ) ( t吨 ) , 小时 ( ) ( t吨 ) , t吨 ) + 小时 = 1 k个 小时 ( ) H(H) ¯ k个 小时 ( ) ( e(电子) k个 ( ) ( t吨 ) , e(电子) 小时 ( ) ( t吨 ) , t吨 ) + 小时 = + 1 N个 c(c) k个 小时 ( ) H(H) k个 小时 ( ) ( x个 k个 ( ) ( t吨 ) , x个 小时 ( ) ( t吨 ) , t吨 ) + u个 k个 ( ) ( t吨 ) ] d日 t吨 + [ ψ k个 ( ) ( k个 ( ) ( t吨 ) , t吨 ) ψ k个 ( ) ( x个 k个 ( ) ( t吨 ) , t吨 ] d日 W公司 ( t吨 ) , 1 , 1 k个 ,
哪里 H(H) ¯ k个 小时 ( ) ( e(电子) k个 ( ) ( t吨 ) , e(电子) 小时 ( ) ( t吨 ) , t吨 ) = H(H) k个 小时 ( ) ( k个 ( ) ( t吨 ) , 小时 ( ) ( t吨 ) , t吨 ) H(H) k个 小时 ( ) ( x个 k个 ( ) ( t吨 ) , x个 小时 ( ) ( t吨 ) , t吨 ) .
假设以下条件适用于每个 1 1 k个 .
  • (A1)有常数 α k个 ( ) j个 , 1 j个 如下不等式适用于任何 x个 k个 , k个 R(右) :
    e(电子) k个 ( ) T型 ϕ k个 ( ) ( k个 , t吨 ) ϕ k个 ( ) ( x个 k个 , t吨 ) j个 = 1 α k个 ( ) j个 e(电子) k个 ( j个 ) 2 .
  • (A2)假设 ψ k个 ( ) ( k个 ( ) , t吨 ) ψ k个 ( ) ( x个 k个 ( ) , t吨 ) 被限定为 x个 k个 ( ) , k个 ( ) R(右) 存在一个常数 β k个 ( ) 这样的话
    ψ k个 ( ) ( k个 ( ) , t吨 ) ψ k个 ( ) ( x个 k个 ( ) , t吨 ) β k个 ( ) k个 ( ) x个 k个 ( ) .
  • (A3)存在正常数 B类 k个 小时 ( ) D类 k个 小时 ( ) 这样的话
    | H(H) ¯ k个 小时 ( ) ( e(电子) k个 ( ) , e(电子) 小时 ( ) , t吨 ) | B类 k个 小时 ( ) | e(电子) k个 ( ) | + D类 k个 小时 ( ) | e(电子) 小时 ( ) | , e(电子) k个 ( ) , e(电子) 小时 ( ) R(右) .
  • (A4)假设每个 ( 1 ) k个 ( 1 k个 ) , { H(H) k个 小时 ( ) ( k个 ( ) , 小时 ( ) , t吨 ) } 小时 = 1 N个 在轨道上线性无关 { 小时 ( ) ( t吨 ) } 小时 = 1 N个 外部同步歧管的 { x个 小时 ( ) ( t吨 ) = 小时 ( ) ( t吨 ) } 小时 = 1 .
接下来,我们给出了一些自适应钉扎控制器和更新律 1 , 1 k个 .
u个 k个 ( ) ( t吨 ) = d日 k个 ( ) ( t吨 ) e(电子) k个 ( ) ( t吨 ) ,
d日 ˙ k个 ( ) ( t吨 ) = 第页 k个 ( ) e(电子) k个 ( ) ( t吨 ) T型 e(电子) k个 ( ) ( t吨 ) ,
b条 ˙ k个 小时 ( ) ( t吨 ) = δ k个 小时 ( ) e(电子) k个 ( ) ( t吨 ) T型 H(H) k个 小时 ( ) ( k个 ( ) ( t吨 ) , 小时 ( ) ( t吨 ) , t吨 ) , 1 小时 , η k个 小时 ( ) e(电子) k个 ( ) ( t吨 ) T型 H(H) k个 小时 ( ) ( x个 k个 ( ) ( t吨 ) , x个 小时 ( ) ( t吨 ) , t吨 ) , + 1 小时 N个 ,
哪里 第页 k个 ( ) , δ k个 小时 ( ) η k个 小时 ( ) 是任意正常数。
定理 1
如果(A1)–(A4)保持不变 ( M(M) , ( k个 小时 ( ) D类 k个 小时 ( ) ) × ) ( 1 ) 每个都紧密相连 ( 1 ) ,然后是未知的局部拓扑结构 ( M(M) , ( k个 小时 ( ) ) × N个 ) ( 1 ) 耦合网络的(1)可以通过以下方式识别 ( M(M) , ( b条 k个 小时 ( ) ) × N个 ) ( 1 ) 控制器下方(5)和更新法律(6)概率为1。也就是说,它适用于每个 ( 1 ) 那个
t吨 k个 = 1 小时 = 1 N个 | b条 k个 小时 ( ) ( t吨 ) k个 小时 ( ) | = 0 , . .
证明。 
我们定义
V(V) k个 ( ) ( e(电子) k个 ( ) , t吨 ) = 1 2 e(电子) k个 ( ) T型 e(电子) k个 ( ) + 1 2 小时 = 1 c(c) k个 小时 ( ) 2 δ k个 小时 ( ) + 1 2 小时 = + 1 N个 c(c) k个 小时 ( ) 2 η k个 小时 ( ) + 1 2 第页 k个 ( ) d日 k个 ( ) d日 * 2 ,
在哪儿 d日 * 是一个足够大的正数。然后,它从微分算子的定义出发 那个
V(V) k个 ( ) ( e(电子) k个 ( ) ( t吨 ) , t吨 ) = e(电子) k个 ( ) ( t吨 ) T型 [ ϕ k个 ( ) ( k个 ( t吨 ) , t吨 ) ϕ k个 ( ) ( x个 k个 ( t吨 ) , t吨 ) d日 k个 ( ) ( t吨 ) e(电子) k个 ( ) ( t吨 ) + 小时 = 1 c(c) k个 小时 ( ) H(H) k个 小时 ( ) ( k个 ( ) ( t吨 ) , 小时 ( ) ( t吨 ) , t吨 ) + 小时 = 1 k个 小时 ( ) H(H) ¯ k个 小时 ( ) ( e(电子) k个 ( ) ( t吨 ) , e(电子) 小时 ( ) ( t吨 ) , t吨 ) + 小时 = + 1 N个 c(c) k个 小时 ( ) H(H) k个 小时 ( ) ( x个 k个 ( ) ( t吨 ) , x个 小时 ( ) ( t吨 ) , t吨 ) ] + 1 2 跟踪 ψ k个 ( ) ( k个 ( ) ( t吨 ) , t吨 ) ψ k个 ( ) ( x个 k个 ( ) ( t吨 ) , t吨 ) T型 × ψ k个 ( ) ( k个 ( ) ( t吨 ) , t吨 ) ψ k个 ( ) ( x个 k个 ( ) ( t吨 ) , t吨 ) 小时 = 1 c(c) k个 小时 ( ) e(电子) k个 ( ) ( t吨 ) T型 H(H) k个 小时 ( ) ( k个 ( ) ( t吨 ) , 小时 ( ) ( t吨 ) , t吨 ) 小时 = + 1 N个 c(c) k个 小时 ( ) e(电子) k个 ( ) ( t吨 ) T型 H(H) k个 小时 ( ) ( x个 k个 ( ) ( t吨 ) , x个 小时 ( ) ( t吨 ) , t吨 ) + d日 k个 ( ) d日 * e(电子) k个 ( ) ( t吨 ) T型 e(电子) k个 ( ) ( t吨 ) j个 = 1 α k个 ( ) j个 e(电子) k个 ( j个 ) ( t吨 ) 2 + e(电子) k个 ( ) ( t吨 ) T型 小时 = 1 k个 小时 ( ) H(H) ¯ k个 小时 ( ) ( e(电子) k个 ( ) ( t吨 ) , e(电子) 小时 ( ) ( t吨 ) , t吨 ) d日 * e(电子) k个 ( ) ( t吨 ) T型 e(电子) k个 ( ) ( t吨 ) + 1 2 β k个 ( ) 2 e(电子) k个 ( ) ( t吨 ) 2 j个 = 1 α k个 ( ) j个 e(电子) k个 ( j个 ) ( t吨 ) 2 d日 * e(电子) k个 ( ) ( t吨 ) 2 + 1 2 β k个 ( ) 2 e(电子) k个 ( ) ( t吨 ) 2 + 小时 = 1 k个 小时 ( ) e(电子) k个 ( ) ( t吨 ) B类 k个 小时 ( ) e(电子) k个 ( ) ( t吨 ) + D类 k个 小时 ( ) e(电子) 小时 ( ) ( t吨 ) = j个 = 1 α k个 ( ) j个 e(电子) k个 ( j个 ) ( t吨 ) 2 + 1 2 小时 = 1 k个 小时 ( ) D类 k个 小时 ( ) e(电子) 小时 ( ) ( t吨 ) 2 e(电子) k个 ( ) ( t吨 ) 2 + 小时 = 1 k个 小时 ( ) B类 k个 小时 ( ) + D类 k个 小时 ( ) + 1 2 β k个 ( ) 2 d日 * e(电子) k个 ( ) ( t吨 ) 2 .
我们定义
V(V) ( e(电子) , t吨 ) = k个 = 1 = 1 λ k个 ( ) V(V) k个 ( ) ( e(电子) k个 ( ) , t吨 ) ,
哪里 λ k个 ( ) 是的辅因子k个-的拉普拉斯矩阵的第th个对角元素 M(M) , k个 小时 ( ) D类 k个 小时 ( ) × ( 1 ) 那么,不难推导出
V(V) ( e(电子) ( t吨 ) , t吨 ) = k个 = 1 = 1 λ k个 ( ) V(V) k个 ( ) ( e(电子) k个 ( ) ( t吨 ) , t吨 ) k个 = 1 = 1 λ k个 ( ) j个 = 1 α k个 ( ) j个 e(电子) k个 ( j个 ) ( t吨 ) 2 + 1 2 k个 = 1 = 1 λ k个 ( ) 小时 = 1 k个 小时 ( ) D类 k个 小时 ( ) e(电子) 小时 ( ) ( t吨 ) 2 e(电子) k个 ( ) ( t吨 ) 2 + k个 = 1 = 1 λ k个 ( ) 小时 = 1 k个 小时 ( ) B类 k个 小时 ( ) + D类 k个 小时 ( ) + 1 2 β k个 ( ) 2 d日 * e(电子) k个 ( ) ( t吨 ) 2 + + .
现在,我们计算两个部分, .
= k个 = 1 j个 = 1 λ k个 ( j个 ) = 1 α k个 ( j个 ) e(电子) k个 ( ) ( t吨 ) 2 = k个 = 1 = 1 λ k个 ( ) 1 λ k个 ( ) j个 = 1 λ k个 ( j个 ) α k个 ( j个 ) e(电子) k个 ( ) ( t吨 ) 2 .
根据引理2,它产生
k个 = 1 = 1 小时 = 1 λ k个 ( ) k个 小时 ( ) D类 k个 小时 ( ) F类 k个 小时 ( ) ( e(电子) k个 ( ) ( t吨 ) , e(电子) 小时 ( ) ( t吨 ) ) = = 1 ω ( ) ( v(v) , u个 ) 电子 ( C类 ) F类 u个 v(v) ( ) ( e(电子) u个 ( ) ( t吨 ) , e(电子) v(v) ( ) ( t吨 ) ) ,
哪里 是的所有生成单圈图的集合 M(M) , k个 小时 ( ) D类 k个 小时 ( ) × , ω ( ) 是的重量 C类 是指 .采取
F类 k个 小时 ( ) ( e(电子) k个 ( ) ( t吨 ) , e(电子) 小时 ( ) ( t吨 ) ) = 1 2 e(电子) 小时 ( ) ( t吨 ) 2 e(电子) k个 ( ) ( t吨 ) 2 ,
它认为
= 1 2 = 1 ω ( ) ( u个 , v(v) ) 电子 ( C类 ) e(电子) u个 ( ) ( t吨 ) 2 e(电子) v(v) ( ) ( t吨 ) 2 = 0 .
然后,可以获得
V(V) ( e(电子) ( t吨 ) , t吨 ) k个 = 1 = 1 λ k个 ( ) 1 λ k个 ( ) j个 = 1 λ k个 ( j个 ) α k个 ( j个 ) e(电子) k个 ( ) ( t吨 ) 2 + k个 = 1 = 1 λ k个 ( ) 小时 = 1 k个 小时 ( ) B类 k个 小时 ( ) + D类 k个 小时 ( ) + 1 2 β k个 ( ) 2 d日 * e(电子) k个 ( ) ( t吨 ) 2 = k个 = 1 = 1 λ k个 ( ) 1 λ k个 ( ) j个 = 1 λ k个 ( j个 ) α k个 ( j个 ) + 小时 = 1 k个 小时 ( ) B类 k个 小时 ( ) + D类 k个 小时 ( ) + 1 2 β k个 ( ) 2 d日 * e(电子) k个 ( ) ( t吨 ) 2 .
作为 d日 * 足够大,存在 σ k个 ( ) = j个 = 1 λ k个 ( j个 ) α k个 ( j个 ) / λ k个 ( ) + 小时 = 1 k个 小时 ( ) B类 k个 小时 ( ) + D类 k个 小时 ( ) + β k个 ( ) 2 / 2 d日 * < 0 因此,
V(V) ( e(电子) ( t吨 ) , t吨 ) k个 = 1 = 1 λ k个 ( ) σ k个 ( ) e(电子) k个 ( ) ( t吨 ) 2 ρ e(电子) T型 ( t吨 ) e(电子) ( t吨 ) ξ ( e(电子) ( t吨 ) ) ,
哪里 ρ > 0 是一个常数,由 λ k个 ( ) σ k个 ( ) , 1 , 1 k个 此外,上述分析表明
| e(电子) | inf公司 0 t吨 < V(V) ( e(电子) , t吨 ) = .
因此,从引理1, t吨 V(V) ( e(电子) ; t吨 ) 存在并且几乎肯定是有限的。它还认为 t吨 ξ ( e(电子) ( t吨 ) ) = 0 . . . 结合拉萨尔不变性原理、假设(A4)和误差系统(4),可以得到集合 M(M) = { e(电子) = 0 , b条 k个 小时 ( ) = k个 小时 ( ) , d日 k个 ( ) = d日 * , 1 , 1 k个 , 1 小时 N个 } 是最大的不变集 M(M) = { ξ ( e(电子) ) = 0 } = { e(电子) = 0 } 因此,对于误差系统的任何初始值(4),轨迹渐近收敛到 M(M) 概率为1[21]. 证明了具有多重扩散的随机多群模型(1)和()自适应控制器下渐近实现完全外同步(5)和更新法律(6). 此外,未知的多重拓扑结构 M(M) , k个 小时 ( ) × N个 ( 1 ) 已由成功识别 M(M) , b条 k个 小时 ( ) × N个 ( 1 ) 概率为1,这就完成了这个证明。
备注 1
众所周知,李亚普诺夫方法在局部拓扑识别的研究中起着重要的作用[17,18,19]. 然而,本文考虑了多重扩散和随机扰动,这使得数学模型更加复杂。因此,对于多群模型,很难直接构造全局Lyapunov函数。动机[20,23,24],我们通过顶点Lyapunov函数的加权求和来构造全局Lyapunow函数 V(V) k个 以…的形式 V(V) = k个 = 1 = 1 λ k个 V(V) k个 在这里, λ k个 ( ) 是拉普拉斯矩阵第k对角元素的余因子 M(M) , k个 小时 ( ) D类 k个 小时 ( ) × ( 1 ) 显然,这种方法与网络的拓扑结构密切相关。由于这种方法使用了图论的一些结果,所以通常称为图论方法。该方法可用于研究许多其他网络的动态行为。例如,在[20]利用图论方法研究了具有扩散的单种群生态模型的稳定性。他们获得的全球稳定结果比[25,26]. 此外,利用图理论方法还可以得到耦合振子和多匹配捕食-被捕食模型的全局渐近稳定性。
备注 2
虽然自适应控制方法可以实现多组模型的拓扑识别,但它需要为每个组添加一个控制器[27,28]. 然而,在本文中,该模型包含了大量组。将控制器添加到所有组有时很难实施,并且控制成本可能更高。因此,本文提出的钉扎控制是有用的。一方面,它可以降低控制成本,因为只需要控制一小部分群体。另一方面,更可行的方法是只向感兴趣的组添加控制器来识别相应的未知拓扑结构。
什么时候? = N个 ,我们可以获得整个拓扑识别。具体来说,驱动系统的相应响应系统(1)可以通过以下特征
d日 k个 ( ) ( t吨 ) = ϕ k个 ( ) ( k个 ( t吨 ) , t吨 ) + 小时 = 1 N个 b条 k个 小时 ( ) H(H) k个 小时 ( ) ( k个 ( ) ( t吨 ) , 小时 ( ) ( t吨 ) , t吨 ) + u个 k个 ( ) ( t吨 ) d日 t吨 + ψ k个 ( ) ( k个 ( ) ( t吨 ) , t吨 ) d日 W公司 ( t吨 ) , 1 , 1 k个 N个 .
然后,驱动系统之间的错误系统(1)和响应系统(10)可以表示为
d日 e(电子) k个 ( ) ( t吨 ) = [ ϕ k个 ( ) ( k个 ( t吨 ) , t吨 ) ϕ k个 ( ) ( x个 k个 ( t吨 ) , t吨 ) + 小时 = 1 N个 b条 k个 小时 ( ) H(H) k个 小时 ( ) ( k个 ( ) ( t吨 ) , 小时 ( ) ( t吨 ) , t吨 ) 小时 = 1 N个 k个 小时 ( ) H(H) k个 小时 ( ) ( x个 k个 ( ) ( t吨 ) , x个 小时 ( ) ( t吨 ) , t吨 ) + u个 k个 ( ) ( t吨 ) ] d日 t吨 + ψ k个 ( ) ( k个 ( ) ( t吨 ) , t吨 ) ψ k个 ( ) ( x个 k个 ( ) ( t吨 ) , t吨 ) d日 W公司 ( t吨 ) , 1 , 1 k个 N个 .
自适应控制器及其更新律 1 , 1 k个 N个 如下所示:
u个 k个 ( ) ( t吨 ) = d日 k个 ( ) ( t吨 ) e(电子) k个 ( ) ( t吨 ) , d日 ˙ k个 ( ) ( t吨 ) = 第页 k个 ( ) e(电子) k个 ( ) ( t吨 ) T型 e(电子) k个 ( ) ( t吨 ) ,
b条 ˙ k个 小时 ( ) ( t吨 ) = δ k个 小时 ( ) e(电子) k个 ( ) ( t吨 ) T型 H(H) k个 小时 ( ) ( k个 ( ) ( t吨 ) , 小时 ( ) ( t吨 ) , t吨 ) , 1 小时 N个 ,
哪里 第页 k个 ( ) n个 d日 δ k个 小时 ( ) 是任意正常数。然后,我们得到以下推论。
推论 1
如果(A1)–(A4)保持不变 ( M(M) , ( k个 小时 ( ) D类 k个 小时 ( ) ) N个 × N个 ) ( 1 ) 每个都紧密相连 ( 1 ) ,然后是未知的整体拓扑结构 ( M(M) , ( k个 小时 ( ) ) N个 × N个 ) ( 1 ) 耦合网络的(1)可以通过以下方式识别 ( M(M) , ( b条 k个 小时 ( ) ) N个 × N个 ) ( 1 ) 控制器下方(12)和更新法律(13). 也就是说,它适用于每个 ( 1 ) 那个
t吨 k个 = 1 N个 小时 = 1 N个 | b条 k个 小时 ( ) ( t吨 ) k个 小时 ( ) | = 0 , . .
这个证明类似于定理1,我们在这里省略了它。
如果 H(H) k个 小时 ( ) ( x个 k个 ( ) ( t吨 ) , x个 小时 ( ) ( t吨 ) , t吨 ) = Γ x个 小时 ( ) ( t吨 ) ,则可以在中找到相应的整个拓扑识别结果[27]. 因此,本文的数学模型更具一般性。此外,本文的理论结果是通用的,可以得到局部拓扑识别和整体拓扑识别。

4.仿真结果

在这一部分中,给出了两个仿真实例来验证理论结果的有效性。我们使用Lorenz系统来描述每个组的动态。显然,洛伦兹系统满足(A1)[29].
例子 1
我们考虑一个具有四个群和三种扩散的一般耦合系统。驱动系统可以采用以下形式:
d日 x个 k个 ( ) ( t吨 ) = ϕ k个 ( ) ( x个 k个 ( t吨 ) , t吨 ) + 小时 = 1 4 k个 小时 ( ) Γ x个 小时 ( ) ( t吨 ) d日 t吨 + ψ k个 ( ) ( x个 k个 ( ) ( t吨 ) , t吨 ) d日 W公司 ( t吨 ) , = 1 , 2 , , k个 = 1 , 2 , , 4 ,
哪里 x个 k个 ( t吨 ) = ( x个 k个 ( 1 ) ( t吨 ) , x个 k个 ( 2 ) ( t吨 ) , x个 k个 ( ) ( t吨 ) ) T型 ,   Γ 1 = Γ 2 = Γ = d日 0.1 , 0.1 , 0.1 .
  • ϕ k个 ( 1 ) ( x个 k个 ( t吨 ) , t吨 ) = 10 x个 k个 ( 1 ) ( t吨 ) x个 k个 ( 2 ) ( t吨 ) , ψ k个 ( 1 ) ( x个 k个 ( 1 ) ( t吨 ) , t吨 ) = 0.5 x个 k个 ( 1 ) ( t吨 ) ;
  • ϕ k个 ( 2 ) ( x个 k个 ( t吨 ) , t吨 ) = 28 x个 k个 ( 1 ) ( t吨 ) x个 k个 ( 2 ) ( t吨 ) x个 k个 ( 1 ) ( t吨 ) x个 k个 ( ) ( t吨 ) , ψ k个 ( 2 ) ( x个 k个 ( 2 ) ( t吨 ) , t吨 ) = 0.3 余弦 x个 k个 ( 2 ) ( t吨 ) ;
  • ϕ k个 ( ) ( x个 k个 ( t吨 ) , t吨 ) = x个 k个 ( 1 ) ( t吨 ) x个 k个 ( 2 ) ( t吨 ) 8 x个 k个 ( ) ( t吨 ) / , ψ k个 ( ) ( x个 k个 ( ) ( t吨 ) , t吨 ) = 0.1 余弦 x个 k个 ( ) ( t吨 ) .
一些砝码 k个 小时 ( ) 配置矩阵的 A类 ( ) = ( k个 小时 ( ) ) 4 × 4 ( = 1 , 2 , , k个 = 1 , 2 , , 4 ) 可以任意选择为
  • 12 ( 1 ) = 34 ( 1 ) = 43 ( 1 ) = 24 ( 2 ) = 31 ( 2 ) = 24 ( ) = 1 ; 13 ( 1 ) = 21 ( 1 ) = 34 ( 2 ) = 12 ( ) = 34 ( ) = 43 ( ) = 2 ;
  • 32 ( 1 ) = 41 ( 1 ) = 12 ( 2 ) = 43 ( 2 ) = 31 ( ) = ; 42 ( 1 ) = 21 ( 2 ) = 4 ;
    其他值设置为0。在不损失通用性的情况下,我们将控制器添加到前两个组中。然后,我们只需要识别部分权重配置矩阵 A类 ¯ ( ) = ( k个 小时 ( ) ) 2 × 4 ( = 1 , 2 , , k个 = 1 , 2 , 小时 = 1 , 2 , , 4 ) . 因此,具有自适应钉扎控制器的响应系统可以表示为
d日 k个 ( ) ( t吨 ) = ϕ k个 ( ) ( k个 ( t吨 ) , t吨 ) + 小时 = 1 2 b条 k个 小时 ( ) Γ 小时 ( ) ( t吨 ) + 小时 = 4 b条 k个 小时 ( ) Γ x个 小时 ( ) ( t吨 ) + u个 k个 ( ) ( t吨 ) d日 t吨 + ψ k个 ( ) ( k个 ( ) ( t吨 ) , t吨 ) d日 W公司 ( t吨 ) , = 1 , 2 , , k个 = 1 , 2 .
我们定义 e(电子) k个 ( ) ( t吨 ) = k个 ( ) ( t吨 ) x个 k个 ( ) ( t吨 ) ( = 1 , 2 , , k个 = 1 , 2 ) . 然后,误差系统(14)和(15)可以显示为
d日 e(电子) k个 ( ) ( t吨 ) = [ ϕ k个 ( ) ( k个 ( t吨 ) , t吨 ) ϕ k个 ( ) ( x个 k个 ( t吨 ) , t吨 ) + 小时 = 1 2 c(c) k个 小时 ( ) H(H) k个 小时 ( ) ( k个 ( ) ( t吨 ) , 小时 ( ) ( t吨 ) , t吨 ) + 小时 = 1 2 k个 小时 ( ) H(H) ¯ k个 小时 ( ) ( e(电子) k个 ( ) ( t吨 ) , e(电子) 小时 ( ) ( t吨 ) , t吨 ) + u个 k个 ( ) ( t吨 ) + 小时 = 4 c(c) k个 小时 ( ) H(H) k个 小时 ( ) ( x个 k个 ( ) ( t吨 ) , x个 小时 ( ) ( t吨 ) , t吨 ) ] d日 t吨 + [ ψ k个 ( ) ( k个 ( ) ( t吨 ) , t吨 ) ψ k个 ( ) ( x个 k个 ( ) ( t吨 ) , t吨 ] d日 W公司 ( t吨 ) , = 1 , 2 , , k个 = 1 , 2 .
可从定理1获得三个权重配置矩阵 A类 ¯ ( ) = ( k个 小时 ( ) ) 2 × 4 可以通过以下方式进行估算 B类 ( ) = ( b条 k个 小时 ( ) ) 2 × 4 根据以下控制器和更新规则 = 1 , 2 , , k个 = 1 , 2 .
u个 k个 ( ) ( t吨 ) = d日 k个 ( ) ( t吨 ) e(电子) k个 ( ) ( t吨 ) , d日 ˙ k个 ( ) ( t吨 ) = 第页 k个 ( ) e(电子) k个 ( ) ( t吨 ) T型 e(电子) k个 ( ) ( t吨 ) ,
b条 ˙ k个 小时 ( ) ( t吨 ) = δ k个 小时 ( ) e(电子) k个 ( ) ( t吨 ) T型 Γ 小时 ( ) ( t吨 ) , 小时 = 1 , 2 , η k个 小时 ( ) e(电子) k个 ( ) ( t吨 ) T型 Γ x个 小时 ( ) ( t吨 ) , 小时 = , 4 .
d日 1 ( ) = 1 , d日 2 ( ) = 2 , x个 k个 ( ) ( 0 ) = k个 ( ) ( 0 ) = 0 , 第页 k个 ( ) = δ k个 小时 ( ) = η k个 小时 ( ) = 1 , b条 k个 小时 ( ) ( 0 ) = 0.1 , ( = 1 , 2 , , k个 = 1 , 2 , 小时 = 1 , 2 , , 4 ) , 我们可以得到一些模拟结果。驱动系统同步错误(14)和响应系统(15)如所示图2.图2a显示第一组的同步错误图2b表示第二组的同步错误。显然,随着时间的推移,所有误差曲线都接近0,这意味着驱动响应系统可以达到同步。图3给出了时间演变 b条 k个 小时 ( ) ( = 1 , 2 , , k个 = 1 , 2 , 小时 = 1 , 2 , , 4 ) 在系统中(15).图3a显示了的识别结果 b条 k个 小时 ( 1 ) ( k个 = 1 , 2 , 小时 = 1 , 2 , , 4 ) 观察到曲线稳定在三个常数:2 b条 k个 小时 ( 1 ) b条 21 ( 1 ) ,1用于 b条 12 ( 1 ) 的和0 b条 11 ( 1 ) , b条 14 ( 1 ) , b条 22 ( 1 ) , b条 23 ( 1 ) , b条 24 ( 1 ) .图3b显示了 b条 k个 小时 ( 2 ) ( k个 = 1 , 2 , 小时 = 1 , 2 , , 4 ) 我们可以看到 b条 21 ( 2 ) 倾向于4, b条 12 ( 2 ) 倾向于3, b条 24 ( 2 ) 趋于1且 b条 11 ( 2 ) , b条 13 ( 2 ) , b条 14 ( 2 ) , b条 22 ( 2 ) , b条 23 ( 2 ) 趋于0。图3c显示了 b条 k个 小时 ( ) ( k个 = 1 , 2 , 小时 = 1 , 2 , , 4 ) 如我们所见,曲线 b条 12 ( ) 收敛到2,曲线 b条 24 ( ) 收敛到1和曲线 b条 11 ( ) , b条 13 ( ) , b条 14 ( ) , b条 21 ( ) , b条 22 ( ) , b条 23 ( ) 收敛到0。因此 b条 k个 小时 ( ) ( = 1 , 2 , , k个 = 1 , 2 , 小时 = 1 , 2 , , 4 ) 可以收敛到一个实值,这意味着局部拓扑识别是成功的。
例子 2
在这个例子中,我们假设网络结构遵循Newman和Watts提出的小世界算法[30]. 考虑系统(14)由八个组组成,固定组的数量与上面的示例相同。为了更好地查看,拓扑结构如所示图4,其中图4a遵守带参数的小世界算法 ( N个 = 8 , K = 2 , P(P) = 0.1 ) ,图4b遵守带参数的小世界算法 ( N个 = 8 , K = 2 , P(P) = 0.2 ) ,图4c遵守带参数的小世界算法 ( N个 = 8 , K = 2 , P(P) = 0.3 ) 其中,
12 ( 1 ) = 13 ( 1 ) = 17 ( 1 ) = 18 ( 1 ) = 21 ( 1 ) = 23 ( 1 ) = 24 ( 1 ) = 28 ( 1 ) = 31 ( 1 ) = 32 ( 1 ) = 34 ( 1 ) = 35 ( 1 ) = 36 ( 1 ) = 42 ( 1 ) = 43 ( 1 ) = 45 ( 1 ) = 46 ( 1 ) = 53 ( 1 ) = 54 ( 1 ) = 56 ( 1 ) = 57 ( 1 ) = 63 ( 1 ) = 64 ( 1 ) = 65 ( 1 ) = 67 ( 1 ) = 67 ( 1 ) = 68 ( 1 ) = 71 ( 1 ) = 75 ( 1 ) = 76 ( 1 ) = 78 ( 1 ) = 81 ( 1 ) = 82 ( 1 ) = 86 ( 1 ) = 87 ( 1 ) = 1 ; 12 ( 2 ) = 13 ( 2 ) = 14 ( 2 ) = 17 ( 2 ) = 18 ( 2 ) = 21 ( 2 ) = 24 ( 2 ) = 24 ( 2 ) = 28 ( 2 ) = 31 ( 2 ) = 32 ( 2 ) = 34 ( 2 ) = 35 ( 2 ) = 38 ( 2 ) = 41 ( 2 ) = 42 ( 2 ) = 43 ( 2 ) = 45 ( 2 ) = 46 ( 2 ) = 47 ( 2 ) = 53 ( 2 ) = 54 ( 2 ) = 56 ( 2 ) = 57 ( 2 ) = 64 ( 2 ) = 65 ( 2 ) = 67 ( 2 ) = 68 ( 2 ) = 71 ( 2 ) = 74 ( 2 ) = 75 ( 2 ) = 76 ( 2 ) = 78 ( 2 ) = 81 ( 2 ) = 82 ( 2 ) = 83 ( 2 ) = 86 ( 2 ) = 87 ( 2 ) = 1 ; 12 ( ) = 13 ( ) = 15 ( ) = 17 ( ) = 17 ( ) = 18 ( ) = 21 ( ) = 23 ( ) = 24 ( ) = 25 ( ) = 28 ( ) = 31 ( ) = 32 ( ) = 34 ( ) = 35 ( ) = 42 ( ) = 43 ( ) = 45 ( ) = 46 ( ) = 51 ( ) = 52 ( ) = 53 ( ) = 54 ( ) = 56 ( ) = 57 ( ) = 58 ( ) = 61 ( ) = 64 ( ) = 65 ( ) = 67 ( ) = 68 ( ) = 71 ( ) = 75 ( ) = 76 ( ) = 78 ( ) = 81 ( ) = 82 ( ) = 85 ( ) = 86 ( ) = 87 ( ) = 1 .
其他的 k个 小时 ( ) = 0 此外,一些参数可以任意设置如下: x个 1 ( ) ( 0 ) = 0.2 , x个 2 ( ) ( 0 ) = 0.3 , x个 ( ) ( 0 ) = 0.1 , x个 4 ( ) ( 0 ) = 0.6 , x个 5 ( ) ( 0 ) = 0.2 , x个 6 ( ) ( 0 ) = 0.2 , x个 7 ( ) ( 0 ) = 0.3 , x个 8 ( ) ( 0 ) = 0.4 , 1 ( ) ( 0 ) = 0.8 , 2 ( ) ( 0 ) = 0.7 , d日 1 ( ) = 1.5 , d日 2 ( ) = 1 , b条 k个 小时 ( ) ( 0 ) = 1 , ( = 1 , 2 , , k个 = 1 , 2 , 小时 = 1 , 2 , , 8 ) . 其他参数与上述示例相同。
在与上例相同的条件下,可以得出定理1的所有假设都满足。这意味着三个重量配置矩阵 A类 ¯ ( ) = ( k个 小时 ( ) ) 2 × 8 可以通过以下方式进行估算 B类 ( ) = ( b条 k个 小时 ( ) ) 2 × 8 在控制器和更新法则下(17)和(18)的 = 1 , 2 , , k个 = 1 , 2 , 小时 = 1 , 2 , , 8 .
图5显示了的仿真结果 B类 ( ) = ( b条 k个 小时 ( ) ) 2 × 8 ( = 1 , 2 , , k个 = 1 , 2 , 小时 = 1 , 2 , , 8 ) 很明显,所有曲线都稳定在实值,这意味着局部拓扑识别是成功的。

5.结论

总之,本文的内容有四个方面。
(1)
该模型考虑了多群模型、多重扩散和随机扰动。
(2)
利用图论方法,可以间接地构造具有多重离散的随机多群模型的全局李亚普诺夫函数。特别是,该方法可以用于研究大规模复杂网络的许多动力学行为,如耦合振子的稳定性和多匹配捕食者-食饵模型。
(3)
利用钉扎控制可以成功地识别随机多群模型的未知局部拓扑结构。
(4)
通过数值算例,可以看出本文得到的理论结果是有效的。
在本文中,噪声是白噪声。然而,在实际应用中也存在颜色噪声。因此,带彩色噪声的多组模型局部拓扑识别将是未来的工作之一。

作者贡献

概念化、C.Z.和D.X。;方法论、C.Z.和H.C。;验证,N.G。;形式分析,H.Y。;调查,R.L。;书面原稿编制,C.Z.和D.X。;写作审查和编辑,H.C.、H.Y.、R.L.和N.G.所有作者均已阅读并同意手稿的出版版本。

基金

这项工作得到了中央政府地方科技发展指导基金(No.2021ZYD0010)和中央高校基本科研业务费(No.2682020ZT109)的支持。

机构审查委员会声明

不适用。

知情同意书

不适用。

数据可用性声明

本文提供了本研究中的数据。

利益冲突

作者没有任何利益冲突。

术语

R(右) + 正实数。
R(右) n个 n个-维欧几里德空间。
R(右) × n个 一套×n个实矩阵。
电子 ( · ) 数学期望。
T型 矩阵或向量的转置。
第页 ( R(右) + ; R(右) + ) 的家庭 R(右) + -值随机变量y 电子 ( | | 第页 ) < .
C类 2 , 1 ( R(右) n个 × R(右) + ; R(右) + ) 所有非负函数族V(V)(x个,t吨)上的 R(右) n个 × R(右) + 在中连续两次可微的x个一次进去t吨.
( Ω , F类 , { F类 t吨 } t吨 0 , P(P) ) 带过滤的完全概率空间 { F类 t吨 } t吨 0 满足平常条件,即它是右连续的和F0包含所有P-null集。

工具书类

  1. 黄,G。;王J.L。;Zu,J.具有潜伏期分布的多组登革热疾病模型的全球动力学。数学。方法应用。科学。 2015,38, 2703–2718. [谷歌学者] [交叉参考]
  2. 刘振杰。;胡,J。;Wang,L.W.腮腺炎全球复发的建模与分析:一个具有无症状感染、一般接种和暴露分布的多组流行病模型。非线性分析。真实世界应用。 2017,37, 137–161. [谷歌学者] [交叉参考]
  3. 丁,D.Q。;Ding,X.H.具有时滞的多组疫苗接种流行病模型的全局稳定性。非线性分析。真实世界应用。 2011,12, 1991–1997. [谷歌学者] [交叉参考]
  4. 沈,M.W。;Xiao,Y.N.具有接种年龄和感染年龄的多组SVEIR流行病模型的全局稳定性。实际应用。数学。 2016,144, 137–157. [谷歌学者] [交叉参考]
  5. 刘,Y。;刘,A.R。;Li,W.X.具有扩散的随机多群模型的同步平稳分布。神经计算。申请。 2020,32, 5001–5013. [谷歌学者] [交叉参考]
  6. 刘,Y。;P.R.Yu。;朱棣文,D.H。;Su,H.具有散布和电报噪声的随机多群模型的平稳分布。非线性分析。混合系统。 2019,33, 93–103. [谷歌学者] [交叉参考]
  7. Chen,G.C.公司。;魏,F。;Wang,W.B.时滞随机区间系统的有限时间镇定及其在储能电路中的应用。数码产品 2019,8, 175. [谷歌学者] [交叉参考] [绿色版本]
  8. Wu,Y.B。;沈,B。;Ahn,C.K。;Li,W.X.随机复杂网络同步的间歇动态事件触发控制。IEEE传输。电路系统。我是Regul。巴普。 2021,68, 2639–2650. [谷歌学者] [交叉参考]
  9. Wu,Y.B。;陈,B.D。;Li,W.X.基于离散时间状态观测的反馈控制随机耦合系统同步。非线性分析。混合系统。 2017,26, 68–85. [谷歌学者] [交叉参考]
  10. Wu,Y.B。;Wang,Y。;Gunasekaran,N。;Vadivel,R.多智能体系统几乎可以肯定的共识:间歇性噪声。IEEE传输。电路系统。II快速简报 2022,69, 2897–2901. [谷歌学者] [交叉参考]
  11. 郭毅。;Li,Y.J。;Ding,X.H.基于图论方法的具有多重分散和随机扰动的多群模型的全局指数稳定性。菲洛马 2017,31, 5125–5138. [谷歌学者] [交叉参考] [绿色版本]
  12. Luo,T.J.基于离散时间状态观测的反馈控制对具有多重分散和随机扰动的多组模型的稳定性。申请。数学。计算。 2019,354, 396–410. [谷歌学者] [交叉参考]
  13. 郭毅。;丁,X.H。;Li,Y.J.具有分散和随机扰动的受电弓多群模型的随机稳定性。J.弗兰克尔。仪表工程应用。数学。 2016,353, 2980–2998. [谷歌学者] [交叉参考]
  14. 李小杰。;Yang,G.H.不确定复杂动态网络基于模糊逼近的全局钉扎同步控制。IEEE传输。赛博。 2017,47, 873–883. [谷歌学者] [交叉参考]
  15. 刘,H。;李,J。;李振英。;曾振国。;Lv,J.H.通过钉扎脉冲控制实现多路动态网络的层内同步。IEEE传输。赛博。 2022,52, 2110–2122. [谷歌学者] [交叉参考]
  16. 刘,L.F。;刘,K。;Xiang,H.Y。;Liu,Q.Pinning复杂动态网络的脉冲簇同步。物理A 2020,545, 123580. [谷歌学者] [交叉参考]
  17. 李,Z.M。;马,W.Y。;Ma,N.R.利用同步方法对回火分数阶复杂网络进行局部拓扑识别。数学。方法应用。科学。 2021. [谷歌学者] [交叉参考]
  18. 朱S.B。;周,J。;Lu,J.A.通过钉扎机制识别复杂动力网络的局部拓扑。混乱 2018,28, 043108. [谷歌学者] [交叉参考]
  19. 徐智强。;张J.X。;Tang,W.S.参数和通过钉扎控制的复杂延迟网络结构识别。事务处理。仪器测量。控制。 2013,35, 619–624. [谷歌学者] [交叉参考]
  20. 李,M.Y。;Shuai,Z.S.网络上耦合微分方程组的全局稳定性问题。J.差异。埃克。 2010,248, 1–20. [谷歌学者] [交叉参考] [绿色版本]
  21. Mao,X.R.拉萨尔定理的随机版本。J.差异。埃克。 1999,153, 175–195. [谷歌学者] [交叉参考] [绿色版本]
  22. X.R.毛。随机微分方程及其应用第2版。;霍伍德出版社:英国奇切斯特出版社,2007年。[谷歌学者]
  23. 张,C.M。;Shi,L.扩散网络上具有无限时滞的耦合捕食者-食饵系统周期性的图论方法。物理学。A类 2021,561, 125255. [谷歌学者] [交叉参考]
  24. Chen,H.L。;张,C.M。;徐,Q。;Feng,Y.L.基于自适应同步的时变延迟耦合随机多加权复杂网络拓扑识别的图论方法。神经过程。莱特。 2022,54, 181–205. [谷歌学者] [交叉参考]
  25. 贝瑞塔,E。;Takeuchi,Y.具有连续时滞的单种群扩散Volterra模型的全局稳定性。牛市。数学。生物。 1987,49, 431–448. [谷歌学者] [交叉参考]
  26. 卢,Z。;Takeuvhi,Y.单种群离散扩散系统的全局渐近行为。数学杂志。生物。 1993,32, 67–77. [谷歌学者] [交叉参考]
  27. 夏,D。;张,C.M。;姚明,X.P。;Chen,H.L.多重分散随机多群模型拓扑识别的图理论方法。J.上海标准。自然科学大学。 2021,20, 301–319. [谷歌学者]
  28. 姚,X.P。;夏,D。;Zhang,C.M.。基于自适应同步的多权重复杂网络拓扑识别:一种图形理论方法。数学。方法应用。科学。 2021,44, 1570–1584. [谷歌学者] [交叉参考]
  29. Lu,J.H。;陈,G.R.创造了一个新的混沌吸引子。国际法学分会。混乱 2002,12, 659–661. [谷歌学者] [交叉参考]
  30. 瓦茨,D。;Strogatz,S.“小世界”网络的集体动态。自然 1998,393, 440–442. [谷歌学者] [交叉参考]
图1。模型的四群有向图(1)带有两个分散器。
图1。模型的四群有向图(1)带有两个分散器。
分形06 00371 g001
图2。驱动系统同步误差的一些样本路径(14)和响应系统(15).
图2。驱动系统同步误差的一些样本路径(14)和响应系统(15).
分形06 00371 g002
图3。耦合构型矩阵的识别 A类 ¯ ( ) ( = 1 , 2 , ) 在(14).
图3。耦合构型矩阵的识别 A类 ¯ ( ) ( = 1 , 2 , ) 在(14).
分形06 00371 g003
图4。具有8个顶点的小世界网络的拓扑结构。
图4。具有8个顶点的小世界网络的拓扑结构。
分形06 00371 g004
图5。耦合构型矩阵的识别 A类 ¯ ( ) ( = 1 , 2 , ) .
图5。耦合构型矩阵的识别 A类 ¯ ( ) ( = 1 , 2 , ) .
分形06 00371 g005
出版商备注:MDPI对公布的地图和机构关联中的管辖权主张保持中立。

分享和引用

MDPI和ACS样式

张,C。;夏,D。;陈,H。;Yang,H.等人。;李,R。;北卡罗来纳州Gunasekaran。用图论方法识别多分散随机多群模型的局部拓扑结构。分形。 2022,6, 371.https://doi.org/10.3390/fractalfract6070371

AMA风格

Zhang C、Xia D、Chen H、Yang H、Li R、Gunasekaran N。用图论方法识别多分散随机多群模型的局部拓扑结构。分形和分数. 2022; 6(7):371.https://doi.org/10.3390/fractalfract6070371

芝加哥/图拉宾风格

张春梅、丹霞、陈慧玲、杨慧、李冉和纳拉潘·古纳塞卡兰。2022.“通过图论方法识别具有多重分散的随机多组模型的局部拓扑结构”分形和分数6,7号:371。https://doi.org/10.3390/fractalfract6070371

文章指标

返回页首顶部