具有时变时滞的Cohen-Grossberg神经网络的广义同步
Cohen-Grossberg Neural Networks with Time-Varying Delays Generalized Synchronization
摘要: 本文利用辅助系统方法研究了一类具有时变时滞的Cohen-Grossberg型神经网络的广义同步问题。首先,基于李雅普诺夫稳定性理论,得到了保证响应系统和辅助系统之间实现指数同步的充分条件。其次,利用线性矩阵不等式,得到了响应系统和辅助系统有限时间混合外同步的充分条件。随后,根据辅助系统方法得出驱动系统与响应系统广义同步。上述结果同样适用于延迟细胞神经网络,结果具有一般性。最后,给出相应的数值模拟来验证所得结论的有效性。
Abstract: This paper investigates the generalized synchronization problem of a class of Cohen-Grossberg type neural networks with time-varying delays using the auxiliary system method. Firstly, based on Lyapunov stability theory, sufficient conditions are obtained to ensure exponential synchronization between the response system and the auxiliary system. Secondly, using linear matrix inequality, sufficient conditions for finite time mixed external synchronization between the response system and the auxiliary system are obtained, According to the auxiliary system method, the generalized synchronization between the driving system and the response system is obtained. The above re-sults are also applicable to delayed cellular neural networks, and the results are general. Finally, corresponding numerical simulations are provided to verify the effectiveness of the obtained con-clusions.
文章引用:李疏桐, 蔡旭颖, 胡湘涵, 李婷, 张川. 具有时变时滞的Cohen-Grossberg神经网络的广义同步[J]. 应用数学进展, 2024, 13(2): 832-847. https://doi.org/10.12677/AAM.2024.132080

1. 引言

脑网络是由大量神经元相互连接而成的复杂结构。这些连接形成了广泛的通信网络,允许神经元之间传递电信号和化学信号,从而实现各种生理和认知功能。因此,研究脑网络对我们认识不同网络功能障碍在神经异常和神经障碍中所起的作用提供了新的见解。Cohen-Grossberg型神经网络(CGNN)是Cohen和Grossberg首次在 [1] 中提出的一种神经网络,在保密通信,模型识别,优化问题等领域有广泛应用 [2] [3] 。它同时也是用于描述神经元之间相互作用和信息处理的数学模型,对脑网络功能的实现不可或缺。由于信息处理的切换速度有限,在研究神经网络稳定性,收敛性等特性的过程中,时间延迟尤其是时变延迟是不可避免的,因此在神经网络的建模中应该考虑时变时滞 [4] [5] [6] 。

如果复杂网络中各个节点的动态都相同,那么称这样的网络为同构网络。但相比同构网络,异构网络即网络中的节点具有不同的动态,更能代表一般情况。以生态群落中捕食者–猎物的相互作用为例,捕食者和猎物相互影响彼此的进化,每个个体的行为和另一个个体总是不同。此外,捕食者群落的拓扑结构与被捕食者群落的拓扑结构也不同 [7] 。

同步是许多复杂网络中一种典型的集体行为。它关注的是网络中每个节点达到相同稳态的一致行为。直到最近几年,人们的兴趣才转向多层网络的同步 [8] 。根据同步方式的不同可分为完全同步,指数同步,聚类同步,反同步,混合同步等 [9] [10] [11] [12] 。在许多实际问题中,完全同步以难实现。例如,在人类疾病网络中,疾病层和人类层无法实现完全同步。然而,疾病可以和宿主实现和谐共存,这种现象被称为“广义同步” [13] [14] [15] 。广义同步在工程网络,生物系统,社会活动和许多其他领域发挥了重要作用。现实世界中的大量现象暗示着广义同步的出现。

广义同步是完全同步的扩展,它意味着响应系统的状态通过非线性光滑的函数映射与驱动系统的状态同步,但是在现实操作中,该映射是难以寻找的。因此,我们可以用辅助系统方法来巧妙地避开这一难题,当响应系统与辅助系统达到完全同步时,对于任意初始状态不同于辅助系统的响应系统,响应系统将和驱动系统到达广义同步。目前,辅助系统方法已经在实际应用中取得了一定成果,例如,Abarbanel H D I等人用辅助系统方法检测电路之间的混沌同步,展示了电路同步和非同步运动的现象 [4] 。

通常,很多注意力都集中在网络内部的同步上,这被称为“内部同步”。除了内部同步之外,我们还可以观察到复杂网络的其他类型的同步行为,比如两个耦合网络之间的“外部同步”,即无论内部网络是否同步,耦合网络的对应节点都将实现同步。这种现象跨越了多个领域,从果蝇时钟神经元群体之间的集体行为,到生态学中捕食者–猎物群落之间的平衡,再到研究人员网络和教育工作者网络之间的拥塞。最近,人们报道了各种外部同步现象,包括完全外部同步(COS)、自适应外部同步(AOS)、逆外部同步(IOS)和广义外部同步(GOS)。在当前的工作中,我们提出了一种新的外部同步行为,即混合外部同步(MOS)。Wang等 [16] 研究了具有相同拓扑结构和耦合延迟的两个复杂网络之间的混合外部同步。Zheng等 [17] 通过具有脉冲控制效应的自适应反馈控制,分析了两个拓扑结构和输出耦合不相同的复杂网络之间的混合外同步。

此外,在实际工程中,我们感兴趣的是在固定有限时间内的有界同步。因此,引入了有限时间稳定性或同步的概念。存在两个有限时间稳定性和同步的概念。第一个概念是指在初始条件下,对于给定的边界,系统的状态在固定的时间间隔内保持在规定的边界内,这也称为有限时间有界性。第二个概念指系统在有限时间内趋于0。复杂网络的有限时间有界性在实际生活中应用广泛,例如,电网无法实现完全同步,因此希望在给定范围内获得条件,使发动机之间的转子相位差保持不变 [18] 。近年来,复杂动态网络的有界控制发展非常迅速,是控制理论和应用中的重要问题之一 [6] 。在现有的研究结果中,大部分的同步都是定义在无限时间上的,然而,研究CGNN型神经网络的有限时间同步更有意义,也更合理。据我们了解,关于CGNN网络的有限时间稳定性或同步的研究很少。在本文中我们研究有关具有时变时滞的CGNN网络的有限时间混合外同步,填补了这一方面的空白。

本文研究了一类具有时变时滞的Cohen-Grossberg型神经网络的指数同步和有限时间混合外同步问题。本文的创新点在于:

(1) 提出的模型更具有一般性。相比于文献 [8] ,本文研究了两层异构神经网络,更具有普遍性。同时利用了辅助系统方法,将证明驱动–响应系统广义同步这一复杂抽象的问题转化为证明响应–辅助系统的完全同步,使问题的解决变得可行,提供了一种新的研究思路。

(2) 基于李雅普诺夫稳定性理论,研究了具有时变时滞的Cohen-Grossberg型神经网络的指数同步,推导了该网络的指数同步的充分条件。

(3) 利用线性矩阵不等式研究有限时间混合外同步,推导了上述网络有限时间混合外同步的新的稳定性判据。

本文结构如下:第1节介绍了Cohen-Grossberg型神经网络模型,指数同步与有限时间混合外同步的定义,假设和引理。第2节基于李雅普诺夫理论和同步定义,在理论上证明了指数同步和有限时间混合外同步。作为特例,给出了延迟细胞神经网络同步的几个充分条件。第3节通过几个数值例子证明了结论的有效性和可行性。最后,我们在第4节对全文进行总结。

符号说明: R n 表示n维欧氏空间。 R n × n 是所有 n × n 阶实矩阵。 A T 表示矩阵A的转置。 w ˙ ( t ) 代表 w ( t ) 相对于时间t的导数。给定向量 x = [ x 1 T , x 2 T , , x n T ] T R n x 表示x的范数。 I = { 1 , 2 , , n } e ε 表示以e为底的指数。

2. 网络模型的建立

考虑一类具有时变时滞的Cohen-Grossberg神经网络,其动力学方程为

x ˙ i ( t ) = p i ( x i ( t ) ) [ a i ( x i ( t ) ) j = 1 n b i j f j ( x j ( t ) ) j = 1 n c i j f j ( x j ( t τ ( t ) ) ) I i ] . (1)

其中 i Ι = { 1 , 2 , , n } n 2 表示神经网络中节点的数目, x i ( t ) 表示第i个神经元在时刻t的状态变量。 p i ( ) 表示放大函数, a i ( ) 是一个适当的函数, b i j c i j 表示第j个神经元与第i个神经元的连接强度。 f j ( x j ( t ) ) f j ( x j ( t τ ( t ) ) ) 表示第j个神经元在时刻t和 t τ ( t ) 的输出, τ ( t ) 表示时间延迟,并且满足 0 τ ( t ) τ ¯ 0 τ ˙ ( t ) μ < 1 I i 是来自网络外部的输入。

本文将(1)作为驱动系统,对应的响应系统描述如下

y ˙ i ( t ) = d i ( y i ( t ) ) [ m i ( y i ( t ) ) j = 1 n h i j g j ( y j ( t ) ) j = 1 n s i j g j ( y j ( t τ ( t ) ) ) I i ] k i 1 ( y i ( t ) ϕ ( x i ) ) k i 2 ( y i ( t τ ( t ) ) ) ϕ ( x i ( t τ ( t ) ) ) (2)

其中 k i 1 k i 2 为反馈增益, ϕ ( ) 是一个变换。

则辅助系统为

z ˙ i ( t ) = d i ( z i ( t ) ) [ m i ( z i ( t ) ) j = 1 n h i j g j ( z j ( t ) ) j = 1 n s i j g j ( z j ( t τ ( t ) ) ) I i ] k i 1 ( z i ( t ) ϕ ( x i ) ) k i 2 ( z i ( t τ ( t ) ) ) ϕ ( x i ( t τ ( t ) ) ) (3)

定义1 [5] 对于系统(2)和系统(3),如果存在 M > 0 σ > 0 使得

i = 1 n | y i ( t ) z i ( t ) | 2 M e σ t φ y φ z , t > 0 .

其中 φ y φ z = sup τ t 0 i = 1 n | φ y i ( t ) φ z i ( t ) | 2 ,则称系统(2)和(3)实现指数同步。

定义2 [19] 响应系统(2)和辅助系统(3)被认为是相对于 c 1 , c 2 , T 实现了有限时间混合外同步,如果对于

sup τ ¯ θ 0 e ( θ ) 2 c 1 ,

e ( t ) 2 c 2 , t [ 0 , T ] .

其中 c 1 c 2 ,且 e ( t ) = [ e 1 T ( t ) , e 2 T ( t ) , , e N T ( t ) ] T

定义3 [20] [21] 对于驱动系统和响应系统的解 x R n y R n ,如果存在映射 ϕ ( x ) : R n R m ,使得

lim t y ( t ) ϕ ( x ( t ) ) = 0 ,

则说明驱动系统和响应系统达到了广义同步。设广义同步误差系统 E ( t ) = y ( t ) ϕ ( x ( t ) ) ,如果存在一个时间常量 T > 0 ,使得 lim t T E ( t ) = 0 ,且当 t T 时有 E ( t ) 0 ,则称驱动系统和响应系统关于向量函数 ϕ ( x ) 在T时刻达到有限时间广义同步,其中T称为同步时间。

假设1 [22] d i ( ) 是一个可微函数,存在正数 d ¯ i d _ i ,使得 0 < d _ i d i ( ) d ¯ i | d ˙ i ( ) | d ˜ i

假设2 [23] 非线性函数 g i ( ) 满足全局Lipschitz条件,即存在一个常数 L i 1 > 0 ,使得

| g i ( α ) g i ( β ) | L i 1 | α β | , α , β R .

假设3 [24] 存在常数 M i ,使得

| g i ( z ) | M i , z R , i = 1 , 2 , , n .

假设4 [25] 对于 i = 1 , 2 , , n ,函数 d i ( ) 满足全局Lipschitz条件,即存在常数 L i 2 > 0 ,有

| d i ( α ) d i ( β ) | L i 2 | α β | , α , β R .

假设5 [26] 对于函数 w i ( ) = d i ( ) m i ( ) ,存在 γ i > 0 ,使得

w i ( y ) w i ( x ) y x γ i .

其中 x , y R ,且 x y

假设6 I i ( t ) 是一个有界函数,即存在正数 I i + > 0 ,使得 | I i ( t ) | I i +

引理1 [19] (柯西不等式)对于任意对称正定矩阵 M R n × n x , y R n ,有

± 2 x T y x T M x + y T M 1 y .

3. 主要结果

对于响应系统与辅助系统,本文定义误差向量 e i ( t ) = y i ( t ) z i ( t ) ,根据(2)~(3)式,同步误差动力学为

e ˙ i ( t ) = [ d i ( y i ( t ) ) m i ( y i ( t ) ) d i ( z i ( t ) ) m i ( z i ( t ) ) ] + d i ( y i ( t ) ) [ j = 1 n h i j g i ( y j ( t ) ) j = 1 n h i j g i ( z j ( t ) ) + j = 1 n s i j g j ( y j ( t τ ( t ) ) ) j = 1 n s i j g j ( z j ( t τ ( t ) ) ) ] + [ d i ( y i ( t ) ) d i ( z i ( t ) ) ] [ j = 1 n h i j g i ( z j ( t ) ) + j = 1 n s i j g j ( z j ( t τ ( t ) ) ) + I i ] k i 1 e i ( t ) k i 2 e i ( t τ ( t ) ) (4)

3.1. Cohen-Grossberg神经网络的指数同步

定理1:在假设1~3和假设5的条件下,响应系统(2)和辅助系统(3)实现指数同步,由辅助系统方法可知驱动系统(1)和响应系统(2)实现广义同步,如果存在正数 k i 1 , k i 2 > 0 ,且满足

j = 1 n d ¯ j | s j i | L i 1 + k i 2 j = 1 n | q i j | ( 1 μ ) e ε τ < 0 ,

ε 2 r i + P i 2 k i 1 + k i 2 + j = 1 n | q i j | < 0 .

其中 P i = d ¯ i j = 1 n | h i j | L j + j = 1 n d ¯ j | h j i | L i + d ¯ i j = 1 n | s i j | L j + 2 d ¯ i [ j = 1 n ( | h i j | + | s i j | ) M j + I i + ]

证明:设Lyapunov函数为

V ( e ( t ) ) = i = 1 n ( e ε t e i T ( t ) e i ( t ) + j = 1 n | q i j | t τ ( t ) t e ε s e i T ( s ) e i ( s ) d s ) .

V ( e ( t ) ) 沿误差系统关于时间t求导得

V ˙ ( e ( t ) ) = i = 1 n [ ε e ε t e i T ( t ) e i ( t ) + 2 e ε t e i T ( t ) e ˙ i ( t ) + j = 1 n | q i j | ( e ε t e i T ( t ) e i ( t ) ( 1 τ ˙ ( t ) ) e ε ( t τ ( t ) ) e i T ( t τ ( t ) ) e i ( t τ ( t ) ) ]

= i = 1 n [ ε e ε t e i T ( t ) e i ( t ) + 2 e ε t e i T ( t ) [ [ d i ( y i ( t ) ) m i ( y i ( t ) ) d i ( z i ( t ) ) m i ( z i ( t ) ) ] + d i ( y i ( t ) ) [ j = 1 n h i j g j ( y j ( t ) ) j = 1 n h i j g j ( z j ( t ) ) + j = 1 n s i j g j ( y j ( t τ ( t ) ) ) j = 1 n s i j g j ( z j ( t τ ( t ) ) ) ] k i 1 e i ( t ) k i 2 e i ( t τ ( t ) ) + [ d i ( y i ( t ) ) d i ( z i ( t ) ) ] [ j = 1 n h i j g j ( z j ( t ) ) + j = 1 n s i j g j ( z j ( t τ ( t ) ) ) + I i ] ] + j = 1 n | q i j | ( e ε t e i T ( t ) e i ( t ) ( 1 τ ˙ ( t ) ) e ε ( t τ ( t ) ) e i T ( t τ ( t ) ) e i ( t τ ( t ) ) ]

由假设5知

[ d i ( y i ( t ) ) m i ( y i ( t ) ) d i ( z i ( t ) ) m i ( z i ( t ) ) ] γ i e i ( t ) .

由假设1和假设2知

d i ( y i ( t ) ) [ j = 1 n h i j ( g j ( y j ( t ) ) g j ( z j ( t ) ) ) + j = 1 n s i j ( g j ( y j ( t τ ( t ) ) ) g j ( z j ( t τ ( t ) ) ) ) ] d ¯ i [ j = 1 n | h i j | L j 1 | e j ( t ) | + j = 1 n | s i j | L j 1 | e j ( t τ ( t ) ) | ]

根据假设3与假设6可得

V ˙ ( e ( t ) ) i = 1 n [ ε e ε t e i T ( t ) e i ( t ) 2 γ i e ε t e i T ( t ) e i ( t ) + 2 d ¯ i j = 1 n | h i j | L j 1 e ε t e i T ( t ) | e j ( t ) | + 2 d ¯ i j = 1 n | s i j | L j 1 e ε t e i T ( t ) | e j ( t τ ( t ) ) | + 2 d ˜ i [ j = 1 n ( | h i j | + | s i j | ) M j + I i + ] e ε t e i T ( t ) | e i ( t ) | 2 k i 1 e ε t e i T ( t ) e i ( t ) 2 k i 2 e ε t e i T ( t ) e i ( t τ ( t ) ) + j = 1 n [ | q i j | e ε t e i T ( t ) e i ( t ) | q i j | ( 1 τ ˙ ( t ) ) e ε ( t τ ( t ) ) e i T ( t τ ( t ) ) e i ( t τ ( t ) ) ] ]

由引理1知

V ˙ ( e ( t ) ) i = 1 n [ ε e ε t e i T ( t ) e i ( t ) 2 γ i e ε t e i T ( t ) e i ( t ) + d ¯ i j = 1 n | h i j | L j 1 e ε t e i T ( t ) e i ( t ) + d ¯ i j = 1 n | h i j | L j 1 e ε t e j T ( t ) e j ( t ) + d ¯ i j = 1 n | s i j | L j 1 e ε t e i T ( t ) e i ( t ) + d ¯ i j = 1 n | s i j | L j 1 e ε t e j T ( t τ ( t ) ) e j ( t τ ( t ) ) + [ 2 d ˜ i [ j = 1 n ( | h i j | + | s i j | ) M j + I i + ] e ε t e i T ( t ) e i ( t ) 2 k i 1 e ε t e i T ( t ) e i ( t ) + k i 2 e ε t e i T ( t ) e i ( t ) + k i 2 e ε t e i T ( t τ ( t ) ) e i ( t τ ( t ) )

+ j = 1 n [ | q i j | e ε t e i T ( t ) e i ( t ) | q i j | ( 1 μ ) e ε ( t τ ) e i T ( t τ ( t ) ) e i ( t τ ( t ) ) ] ] i = 1 n [ ε 2 r i + P i 2 k i 1 + k i 2 + j = 1 n | q i j | ] e ε t e i T ( t ) e i ( t ) + i = 1 n [ j = 1 n d ¯ j | s j i | L i 1 + k i 2 j = 1 n | q i j | ( 1 μ ) e ε τ ] e ε t e i T ( t τ ( t ) ) e i ( t τ ( t ) ) < 0

其中 P i = d ¯ i j = 1 n | h i j | L j 1 + j = 1 n d ¯ j | h j i | L i 1 + d ¯ i j = 1 n | s i j | L j 1 + 2 d ˜ i [ j = 1 n ( | h i j | + | s i j | ) M j + I i + ] 。所以对任意的 t 0 ,有

V ( 0 ) = i = 1 n ( e i T ( 0 ) e i ( 0 ) + j = 1 n | q i j | τ ( 0 ) 0 e ε t e i T ( s ) e i ( s ) d s ) m ϕ φ

其中 ϕ φ = sup τ t 0 j = 1 n | ϕ j ( t ) φ j ( t ) | 2 m = 1 + τ ¯ max 1 i n [ j = 1 n | q i j | ] ,故 i = 1 n | e i ( t ) | 2 m e ε t ϕ φ 。因此系统(2)和(3)实现指数同步,由辅助系统方法和定义3可知系统(1)和(2)实现广义同步。

3.2. Cohen-Grossberg神经网络的有限时间混合外同步

定理2:对于给定的标量 τ ¯ > 0 ,在假设1~5的条件下,响应系统(2)和辅助系统(3)可以实现相对于 c 1 , c 2 , T ( c 1 < c 2 ) 的有限时间混合外同步。由辅助系统方法可知驱动系统(1)和响应系统(2)实现广义同步如果存在控制增益常数 k i 1 k i 2 α i > 0 α = max { α i } ,使得以下成立

2 γ i + j = 1 n d ¯ i 2 ( h i j 2 + s i j 2 ) + L i 1 2 + 2 L i 2 [ j = 1 n ( | h i j | + | s i j | ) M j + | I i | ] 2 k i 1 + k i 2 + 1 α i < 0 ,

e α T c 1 ( 1 + τ ¯ ) c 2 ,

L i 1 2 + k i 2 ( 1 μ ) 0 .

其中 i = 1 , 2 , , n

证明:设Lyapunov函数为

V ( e ( t ) ) = i = 1 n e i T ( t ) e i ( t ) + i = 1 n t τ ( t ) t e i T ( s ) e i ( s ) d s .

V ( e ( t ) ) 沿着局部误差系统关于时间t求导得

V ˙ ( e ( t ) ) = 2 i = 1 n e i T ( t ) e ˙ i ( t ) + i = 1 n e i T ( t ) e i ( t ) + i = 1 n ( 1 τ ˙ ( t ) ) e i T ( t τ ( t ) ) e i ( t τ ( t ) ) .

i = 1 n e i T ( t ) e ˙ i ( t ) = i = 1 n e i T ( t ) [ d i ( y i ( t ) ) m i ( y i ( t ) ) d i ( z i ( t ) ) m i ( z i ( t ) ) ] i = 1 n e i T ( t ) d i ( y i ( t ) ) [ j = 1 n h i j g j ( y j ( t ) ) + j = 1 n h i j g j ( z j ( t ) ) j = 1 n s i j g j ( y j ( t τ ( t ) ) ) + j = 1 n s i j g j ( z j ( t τ ( t ) ) ) ] + i = 1 n e i T ( t ) [ d i ( y i ( t ) ) d i ( z i ( t ) ) ] [ j = 1 n h i j g j ( z j ( t ) ) + j = 1 n s i j g j ( z j ( t τ ( t ) ) ) + I i ] i = 1 n e i T ( t ) k i 1 e i ( t ) i = 1 n e i T ( t ) k i 2 e i ( t τ ( t ) ) .

由假设5可得

i = 1 n e i T ( t ) [ d i ( y i ( t ) ) m i ( y i ( t ) ) d i ( z i ( t ) ) m i ( z i ( t ) ) ] i = 1 n e i T ( t ) γ i e i ( t ) .

由假设1、2和引理1得到

i = 1 n e i T ( t ) d i ( y i ( t ) ) [ j = 1 n h i j ( g j ( z j ( t ) ) g j ( y j ( t ) ) ) + j = 1 n s i j ( g j ( z j ( t τ ( t ) ) ) g j ( y j ( t τ ( t ) ) ) ) ] i = 1 n | e i T ( t ) | d ¯ i ( j = 1 n | h i j | L j 1 | e j ( t ) | + j = 1 n | s i j | L j 1 | e j ( t τ ( t ) ) | ) 1 2 i = 1 n j = 1 n [ | e i T ( t ) | d ¯ i 2 h i j 2 | e i ( t ) | + | e j T ( t ) | L j 1 2 | e j ( t ) | + | e i T ( t ) | d ¯ i 2 s i j 2 | e i ( t ) | + | e j T ( t τ ( t ) ) | L j 1 2 | e j ( t τ ( t ) ) | ] .

由假设3、4可知

i = 1 n e i T ( t ) [ d i ( y i ( t ) ) d i ( z i ( t ) ) ] [ j = 1 n h i j g j ( z j ( t ) ) + j = 1 n s i j g j ( z j ( t τ ( t ) ) ) + I i ] i = 1 n e i T ( t ) L i 2 e i ( t ) [ j = 1 n ( | h i j | + | s i j | ) M j + | I i | ] .

由引理1得到

i = 1 n e i T ( t ) k i 1 e i ( t ) i = 1 n e i T ( t ) k i 2 e i ( t τ ( t ) ) i = 1 n e i T ( t ) k i 1 e i ( t ) + 1 2 i = 1 n k i 2 [ e i T ( t ) e i ( t ) + e i T ( t τ ( t ) ) e i ( t τ ( t ) ) ] .

综上所述,

V ˙ ( e ( t ) ) 2 i = 1 n e i T ( t ) γ i e i ( t ) + i = 1 n j = 1 n [ | e i T ( t ) | d ¯ i 2 h i j 2 | e i ( t ) | + | e j T ( t ) | L j 1 2 | e j ( t ) | + | e i T ( t ) | d ¯ i 2 s i j 2 | e i ( t ) | ] + i = 1 n j = 1 n | e j T ( t τ ( t ) ) | L j 1 2 | e j ( t τ ( t ) ) | + 2 i = 1 n e i T ( t ) L i 2 e i ( t ) [ j = 1 n ( | h i j | + | s i j | ) M j + | I i | ] 2 i = 1 n e i T ( t ) k i 1 e i ( t ) + i = 1 n k i 2 e i T ( t ) e i ( t ) + i = 1 n k i 2 e i T ( t τ ( t ) ) e i ( t τ ( t ) ) + i = 1 n e i T ( t ) e i ( t ) + i = 1 n ( 1 τ ˙ ( t ) ) e i T ( t τ ( t ) ) e i ( t τ ( t ) )

= i = 1 n e i T ( t ) { 2 γ i + j = 1 n d ¯ i 2 ( h i j 2 + s i j 2 ) + L i 1 2 + 2 L i 2 [ j = 1 n ( | h i j | + | s i j | ) M j + | I i | ] 2 k i 1 + k i 2 + 1 } e i ( t ) + i = 1 n ( L i 1 2 + k i 2 ( 1 τ ˙ ( t ) ) ) e i T ( t τ ( t ) ) e i ( t τ ( t ) )

则当 L i 1 2 + k i 2 ( 1 μ ) 0 时,

V ˙ ( e ( t ) ) i = 1 n e i T ( t ) Ω i e i ( t )

其中 Ω i = 2 γ i + j = 1 n d ¯ i 2 ( h i j 2 + s i j 2 ) + L i 1 2 + 2 L i 2 [ j = 1 n ( | h i j | + | s i j | ) M j + | I i | ] 2 k i 1 + k i 2 + 1

注意到

V ( e ( t ) ) i = 1 n e i T ( t ) e i ( t ) = e ( t ) 2 .

根据

Ω i α i < 0 ,

Ω i < α i .

α = max { α i } ,

其中 i = 1 , 2 , , n

则有

V ˙ ( e ( t ) ) < i = 1 n e i T ( t ) Ω i e i ( t ) i = 1 n e i T ( t ) α i e i ( t ) α V ( e ( t ) ) .

上式左右分别乘 e α t ,并从0到t积分, t [ 0 , T ] ,则有

e α t V ( e ( t ) ) < V ( e ( 0 ) ) .

由于

V ( e ( 0 ) ) = i = 1 n e i T ( 0 ) e i ( 0 ) + i = 1 n τ ( t ) 0 e i T ( s ) e i ( s ) d s e ( 0 ) 2 + τ ¯ sup τ ¯ θ 0 e ( θ ) 2 c 1 ( 1 + τ ¯ ) .

因此,

e ( t ) 2 e α T c 1 ( 1 + τ ¯ ) .

由定理2知,对于所有的 t [ 0 , T ] ,存在

e ( t ) 2 < c 2 .

于是,响应系统(2)和辅助系统(3)相对于 c 1 , c 2 , T ( c 1 < c 2 )实现了有限时间混合外同步,由辅助系统方法以及定义3可知系统(1)和(2)实现有限时间广义同步。

如果放大函数 p i ( ) = 1 ,则驱动系统(1)变为

x ˙ i ( t ) = a i ( x i ( t ) ) + j = 1 n b i j f j ( x j ( t ) ) + j = 1 n c i j f j x j ( t τ ( t ) ) + I i . (5)

响应系统(2)与辅助系统(3)分别变为

y ˙ i ( t ) = m i ( y i ( t ) ) + j = 1 n h i j g j ( y j ( t ) ) + j = 1 n s i j g j y j ( t τ ( t ) ) + I i k i 1 ( y i ( t ) ϕ ( x i ) ) k i 2 ( y i ( t τ ( t ) ) ) ϕ ( x i ( t τ ( t ) ) ) (6)

z ˙ i ( t ) = m i ( z i ( t ) ) + j = 1 n h i j g j ( z j ( t ) ) + j = 1 n s i j g j z j ( t τ ( t ) ) + I i k i 1 ( z i ( t ) ϕ ( x i ) ) k i 2 ( z i ( t τ ( t ) ) ) ϕ ( x i ( t τ ( t ) ) ) (7)

根据定理1可得以下推论:

推论1:在假设1~3和假设5的条件下,系统(6)与系统(7)实现指数同步,如果

j = 1 n d ¯ j | s j i | L i 1 + k i 2 j = 1 n | q i j | ( 1 μ ) e ε τ < 0 ,

ε 2 γ i + P i 2 k i 1 + k i 2 + j = 1 n | q i j | < 0 .

其中 P i = d ¯ i j = 1 n | h i j | L j + j = 1 n d ¯ j | h j i | L i + d ¯ i j = 1 n | s i j | L j + 2 d ¯ i [ j = 1 n ( | h i j | + | s i j | ) M j + I i + ]

推论2:对于给定的标量 τ ¯ > 0 ,在假设5~10的条件下,系统(6)和系统(7)可以实现相对于 c 1 , c 2 , T ( c 1 < c 2 )的有限时间混合外同步。如果存在控制增益常数 k i 1 k i 2 α i > 0 α = max { α i } ,使得以下成立

2 γ i + j = 1 n ( h i j 2 + s i j 2 ) + L i 1 2 2 k i 1 + k i 2 + 1 α i < 0 ,

e α T c 1 ( 1 + τ ¯ ) c 2 ,

L i 2 2 + k i 2 1 + μ 0 .

其中 i = 1 , 2 , , n

4. 数值仿真

考虑下面的具有时变时滞Cohen-Grossberg型神经网络

y ˙ i ( t ) = d i ( y i ( t ) ) [ m i ( y i ( t ) ) j = 1 5 h i j g j ( y j ( t ) ) j = 1 5 s i j g j y j ( t τ ( t ) ) I i ] k i 1 ( y i ( t ) ϕ ( x i ) ) k i 2 ( y i ( t τ ( t ) ) ) ϕ ( x i ( t τ ( t ) ) ) (8)

其中 i = 1 , 2 , 3 , 4 , 5

辅助系统描述为

z ˙ i ( t ) = d i ( z i ( t ) ) [ m i ( z i ( t ) ) j = 1 5 h i j g j ( z j ( t ) ) j = 1 5 s i j g j z j ( t τ ( t ) ) I i ] k i 1 ( z i ( t ) ϕ ( x i ) ) k i 2 ( z i ( t τ ( t ) ) ) ϕ ( x i ( t τ ( t ) ) ) (9)

4.1. Cohen-Grossberg神经网络的指数同步数值仿真

为了验证定理一结论的正确性,令 g i ( ) = 0.4 tanh ( ) m 1 ( ) = 1.4 ( ) m 2 ( ) = 0.4 ( ) m 3 ( ) = 1.4 ( ) m 4 ( ) = 0.4 ( ) m 5 ( ) = 1.4 ( ) h 11 = 1.7 h 12 = 0.6 h 13 = 0.5 h 14 = 2.5 h 15 = 0.2 h 21 = 0.2 h 22 = 0.14 h 23 = 0.14 h 24 = 0.13 h 25 = 0.1 h 31 = 1.8 h 32 = 0.2 h 33 = 2 h 34 = 0.4 h 35 = 1.7 h 41 = 0.5 h 42 = 0.2 h 43 = 0.17 h 44 = 0.16 h 45 = 0.09 h 51 = 0.1 h 52 = 0.1 h 53 = 0.2 h 54 = 0.14 h 55 = 0.16 s 11 = 0.2 s 12 = 0.14 s 13 = 0.14 s 14 = 0.5 s 15 = 0.2 s 21 = 0.13 s 22 = 0.25 s 23 = 0.1 s 24 = 0.16 s 25 = 0.17 s 31 = 0.1 s 32 = 0.1 s 33 = 0.19 s 34 = 0.1 s 35 = 0.14 s 41 = 1.8 s 42 = 0.1 s 43 = 2 s 44 = 0.4 s 45 = 1.7 s 51 = 0.5 s 52 = 0.2 s 53 = 0.09 s 54 = 0.16 s 55 = 0.2 I i ( t ) = 0

τ ( t ) = e t 1 + e t d i ( ) = 0.7 + 0.1 1 + ( ) 2

所以,可以得到,对于 i = 1 , 2 , 3 , 4 , 5 ,有

0.7 d i ( ) 0.8 , | d ˙ i ( ) | 0.1 ,

另一方面,对于任意的 x , y R x y ,有

d 1 ( y ) m 1 ( y ) d 1 ( x ) m 1 ( x ) y x = 0.98 + 0.14 1 x y ( 1 + x 2 ) ( 1 + y 2 ) 0.98 0.07 ( x 2 + y 2 ) ( 1 + x 2 ) ( 1 + y 2 ) 0.91 .

从而 γ 1 = 0.91 。同理可得 γ 2 = 0.26 γ 3 = 0.91 γ 4 = 0.26 γ 5 = 0.91

经过计算得到

| g i ( ) | 0.4

并且

| g i ( x ) g i ( y ) | 0.4 | x y | .

进一步可得

L i 1 = 0.4 M i = 0.4 d ¯ i = 0.8 d _ i = 0.7 d ˜ i = 0.1 τ ¯ = 1 μ = 0.25 i = 1 , 2 , 3 , 4 , 5 .

选择 k 11 = 4.88 k 12 = 2.08 k 21 = 4.88 k 22 = 2.08 k 31 = 4.88 k 32 = 2.08 k 41 = 5.08 k 42 = 1.08 k 51 = 4.88 k 52 = 2.08 q i j = 1 i = 1 , 2 , 3 , 4 , 5 j = 1 , 2 , 3 , 4 , 5 ε = 0.1 。通过计算得到

ε 2 γ 1 + d ¯ 1 j = 1 5 | h 1 j | L j 1 + j = 1 5 d ¯ j | h j 1 | L 11 + d ¯ 1 j = 1 5 | s 1 j | L 11 + 2 d ˜ 1 [ j = 1 5 ( | h 1 j | + | s 1 j | ) M j + I 1 + ] 2 k 11 + k 12 + j = 1 5 | q 1 j | = 0.352 .

ε 2 γ 2 + d ¯ 2 j = 1 5 | h 2 j | L j 1 + j = 1 5 d ¯ j | h j 2 | L 21 + d ¯ 2 j = 1 5 | s 2 j | L j 1 + 2 d ˜ 2 [ j = 1 5 ( | h 2 j | + | s 2 j | ) M j + I 2 + ] 2 k 21 + k 22 + j = 1 5 | q 2 j | = 2.0952 .

ε 2 γ 3 + d ¯ 3 j = 1 5 | h 1 j | L j 1 + j = 1 5 d ¯ j | h j 3 | L 31 + d ¯ 3 j = 1 5 | s 3 j | L j 1 + 2 d ˜ 3 [ j = 1 5 ( | h 3 j | + | s 3 j | ) M j + I 3 + ] 2 k 31 + k 32 + j = 1 5 | q 3 j | = 0.7448 .

ε 2 γ 4 + d ¯ 4 j = 1 5 | h 4 j | L j 1 + j = 1 5 d ¯ j | h j 4 | L 41 + d ¯ 4 j = 1 5 | s 4 j | L j 1 + 2 d ˜ 4 [ j = 1 5 ( | h 4 j | + | s 4 j | ) M j + I 4 + ] 2 k 41 + k 42 + j = 1 5 | q 4 j | = 0.5864 .

ε 2 γ 5 + d ¯ 5 j = 1 5 | h 5 j | L j 1 + j = 1 5 d ¯ j | h j 5 | L 51 + d ¯ 5 j = 1 5 | s 5 j | L j 1 + 2 d ˜ 5 [ j = 1 5 ( | h 5 j | + | s 5 j | ) M j + I 5 + ] 2 k 51 + k 52 + j = 1 5 | q 5 j | = 4.34 .

d ¯ j j = 1 5 | s 1 j | L j 1 + k 12 j = 1 5 | q 1 j | ( 1 μ ) e ε τ = 0.4395 .

d ¯ j j = 1 5 | s 2 j | L j 1 + k 22 j = 1 5 | q 2 j | ( 1 μ ) e ε τ = 1.0603 .

d ¯ j j = 1 5 | s 3 j | L j 1 + k 32 j = 1 5 | q 3 j | ( 1 μ ) e ε τ = 0.5067 .

d ¯ j j = 1 5 | s 4 j | L j 1 + k 42 j = 1 5 | q 4 j | ( 1 μ ) e ε τ = 1.8907 .

d ¯ j j = 1 5 | s 5 j | L j 1 + k 52 j = 1 5 | q 5 j | ( 1 μ ) e ε τ = 0.5419 .

在定理1的条件下,(8)和(9)是指数同步的,同步误差 e i ( t ) 的时间演化如图1图2所示。因此,由辅助系统方法知驱动系统(1)和响应系统(2)实现广义同步。

Figure 1. Synchronization error e i 1

图1. 同步误差 e i 1

Figure 2. Synchronization error e i 2

图2. 同步误差 e i 2

4.2. Cohen-Grossberg神经网络的有限时间混合外同步数值仿真

为了验证定理2结论的正确性,本文依然选择选择5维的CGNN,其中 h 11 = 1.7 h 12 = 0.6 h 13 = 2 h 14 = 0.2 h 15 = 1 h 21 = 0.5 h 22 = 2.5 h 23 = 0.5 h 24 = 0.5 h 25 = 3 h 31 = 1 h 32 = 0.1 h 33 = 1.5 h 34 = 0.8 h 35 = 1 h 41 = 0.2 h 42 = 0.8 h 43 = 2.5 h 44 = 1.8 h 45 = 0.1 h 51 = 2 h 52 = 0.4 h 53 = 1.7 h 54 = 0.6 h 55 = 0.5 s 11 = 1.8 s 12 = 0.1 s 13 = 2.5 s 14 = 1.7 s 15 = 0.5 s 21 = 2 s 22 = 0.4 s 23 = 0.6 s 24 = 2.5 s 25 = 1.25 s 31 = 3 s 32 = 0.7 s 33 = 1 s 34 = 1.2 s 35 = 1.6 s 41 = 1.3 s 42 = 2 s 43 = 0.5 s 44 = 0.3 s 45 = 1.1 s 51 = 0.3 s 52 = 1 s 53 = 0.6 s 54 = 1.1 s 55 = 0.4 m 1 ( ) = 1.4 ( ) m 2 ( ) = 0.4 ( ) m 3 ( ) = 1.6 ( ) m 4 ( ) = 1.9 ( ) m 5 ( ) = 2 ( ) d i = 0.7 + 0.1 1 + ( ) 2 g i = 0.4 tanh ( ) I i = 0 τ ( t ) = e t 1 + e t i = 1 , 2 , 3 , 4 , 5

所以,我们可以得到:

0.7 < d i 0.8

并且

| g j ( μ ) g j ( ν ) | 0.4 | μ ν | .

进一步可得 L i 1 = 0.4 L i 2 = 0.1 M i = 0.4 g ¯ = 0.4 d ¯ i = 0.8 d _ i = 0.7 τ ¯ = 1 μ = 1 3

φ 1 ( s ) = 0.1 φ 2 ( s ) = 0.2 φ 3 ( s ) = 0.6 φ 4 ( s ) = 0.3 φ 5 ( s ) = 0.5 ϕ 1 ( s ) = 0.1 ϕ 2 ( s ) = 0.4 ϕ 3 ( s ) = 0.5 ϕ 4 ( s ) = 0.1 ϕ 5 ( s ) = 0.2 作为系统的初值,其中, s [ 1 , 0 ]

c 1 = 2.2 c 2 = 5 T = 17 ,利用MATLAB中的LMI控制工具箱对定理2中的LMI进行求解,得到一组可行解:

α 1 = 0.0057 α 2 = 0.0063 α 3 = 0.0065 α 4 = 0.0054 α 5 = 0.0039 ,令 α = max { α i } i = 1 , 2 , 3 , 4 , 5 k 11 = 288.0601 k 12 = 0.2533 k 12 = 350.5336 k 22 = 0 .2533 k 31 = 293 .3588 k 32 = 0.2533 k 41 = 418 .3834 k 42 = 0.2533 k 51 = 37 .6107 k 52 = 0.2536

因此,响应网络(2)与辅助网络(3)可以实现相对于 c 1 , c 2 , T 的有限时间混合外同步,数值模拟如图所示

Figure 3. Synchronization error e i 3

图3. 同步误差 e i 3

Figure 4. Synchronization error e i 4

图4. 同步误差 e i 4

采用欧拉法进行数值模拟,从图3图4中不难看出,误差系统可以渐近收敛于零,这意味着响应系统(2)和辅助系统(3)已经实现了有限时间混合外同步。在图3图4中,我们注意到,响应系统(2)和辅助系统(3)具有快速完美地实现了关于 c 1 , c 2 , T 的有限时间混合外同步,有限时间混合外同步的时间T约为17 s。

5. 结论

本文针对Cohen-Grossberg型神经网络,考虑时变时滞因素,研究其实现广义同步问题。为了得到响应系统和辅助系统实现指数同步的充分条件,基于李雅普诺夫稳定性理论,构造李雅普诺夫函数,从理论上证明了两层网络同步。再结合线性矩阵不等式,得到了响应系统和辅助系统有限时间混合外同步的充分条件,最后通过数值仿真验证了理论的有效性和可行性。此外,本文研究的主要是线性系统问题,在未来的工作中,我们将研究非线性系统的广义同步问题。

基金项目

省级大学生创新训练项目(S202310446045)。

NOTES

*第一作者。

#通讯作者。

参考文献

[1] Cohen, M.A. and Grossberg, S. (1983) Absolute Stability of Global Pattern Formation and Parallel Memory Storage by Competitive Neural Networks. IEEE Transactions on Systems, Man, and Cybernetics, SMC-13, 815-826.
https://doi.org/10.1109/TSMC.1983.6313075
[2] 张雅美, 郝涛, 尹四倍, 等. 具有时延和随机扰动的未知C-G神经网络的有限时间函数投影同步及其在保密通信中的应用[J]. 应用数学和力学, 2020, 41(12): 1405-1416.
[3] Nagamani, G. and Radhika, T. (2016) A Quadratic Convex Combination Approach on Robust Dissi-pativity and Passivity Analysis for Takagi-Sugeno Fuzzy Cohen-Grossberg Neural Networks with Time-Varying De-lays. Mathematical Methods in the Applied Sciences, 39, 3880-3896.
https://doi.org/10.1002/mma.3835
[4] Abarbanel, H.D.I., Rulkov, N.F. and Sushchik, M.M. (1996) Generalized Synchronization of Chaos: The Auxiliary System Approach. Physical Review E, 53, 4528-4535.
https://doi.org/10.1103/PhysRevE.53.4528
[5] Zhang, C., Deng, F., Zhao, X., et al. (2016) P-Th Exponential Synchronization of Cohen-Grossberg Neural Network with Mixed Time-Varying Delays and Unknown Parameters Us-ing Impulsive Control Method. Neurocomputing, 218, 432-438.
https://doi.org/10.1016/j.neucom.2016.09.002
[6] Xu, Y., Zhang, J., Zhou, W., et al. (2017) Finite-Time Bounded Synchronization of the Growing Complex Network with Nondelayed and Delayed Coupling. Discrete Dynamics in Na-ture and Society, 2017, Article ID: 6501583.
https://doi.org/10.1155/2017/6501583
[7] Wu, X., Zheng, W.X. and Zhou, J. (2009) Generalized Outer Synchro-nization between Complex Dynamical Networks. Chaos: An Interdisciplinary Journal of Nonlinear Science, 19, Article ID: 013109.
https://doi.org/10.1063/1.3072787
[8] Huang, J., Wu, X. and Tu, Z. (2023) Exponential Synchroni-zation of Cohen-Grossberg Neural Networks with Delays. IEEE Access, 11, 49992-49999.
https://doi.org/10.1109/ACCESS.2023.3277986
[9] Liang, K. and Wanli, L. (2019) Exponential Synchronization in Inertial Cohen-Grossberg Neural Networks with Time Delays. Journal of the Franklin Institute, 356, 11285-11304.
https://doi.org/10.1016/j.jfranklin.2019.07.027
[10] Jalan, S. and Singh, A. (2016) Cluster Synchronization in Mul-tiplex Networks. Europhysics Letters, 113, 30002.
https://doi.org/10.1209/0295-5075/113/30002
[11] Sundarapandian, V. (2013) Analysis and Anti-Synchronization of a Novel Chaotic System via Active and Adaptive Controllers. Journal of Engineering Science and Technology Review, 6, 45-52.
https://doi.org/10.25103/jestr.064.06
[12] Ma, T. and Zhang, J. (2015) Hybrid Synchronization of Cou-pled Fractional-Order Complex Networks. Neurocomputing, 157, 166-172.
https://doi.org/10.1016/j.neucom.2015.01.022
[13] 韦相, 赵军产, 胡春华. 两个异构复杂网络的广义同步与参数识别[J]. 自动化学报, 2017, 43(4): 595-603.
[14] Zhou, J., Chen, J., Lu, J., et al. (2016) On Applicability of Auxil-iary System Approach to Detect Generalized Synchronization in Complex Network. IEEE Transactions on Automatic Control, 62, 3468-3473.
https://doi.org/10.1109/TAC.2016.2615679
[15] Bao, H. and Cao, J. (2016) Finite-Time Generalized Synchroniza-tion of Nonidentical Delayed Chaotic Systems. Nonlinear Analysis: Modelling and Control, 21, 306-324.
https://doi.org/10.15388/NA.2016.3.2
[16] Wang, J.W., Ma, Q., Zeng, L., et al. (2011) Mixed Outer Synchroniza-tion of Coupled Complex Networks with Time-Varying Coupling Delay. Chaos: An Interdisciplinary Journal of Non-linear Science, 21, Article ID: 013121.
https://doi.org/10.1063/1.3555836
[17] Zheng, S. and Shao, W. (2013) Mixed Outer Synchronization of Dynamical Networks with Nonidentical Nodes and Output Coupling. Nonlinear Dynamics, 73, 2343-2352.
https://doi.org/10.1007/s11071-013-0945-1
[18] Dorfler, F. and Bullo, F. (2012) Synchronization and Transient Stability in Power Networks and Nonuniform Kuramoto Oscillators. SIAM Journal on Control and Optimization, 50, 1616-1642.
https://doi.org/10.1137/110851584
[19] He, P., Ma, S.H. and Fan, T. (2012) Finite-Time Mixed Outer Synchronization of Complex Networks with Coupling Time-Varying Delay. Chaos: An Interdisciplinary Journal of Nonlinear Science, 22, Article ID: 043151.
https://doi.org/10.1063/1.4773005
[20] 王智琦. 多层网络的广义同步控制方法[D]: [硕士学位论文]. 长春: 吉林大学, 2022.
[21] 朱泽飞, 涂俐兰, 吴泽虎. 异维混沌动力系统的有限时间广义同步[J]. 数学杂志, 2017, 37(2): 365-375.
[22] Hu, C., Yu, J. and Jiang, H. (2014) Finite-Time Synchronization of Delayed Neural Networks with Cohen-Grossberg Type Based on Delayed Feedback Control. Neurocomputing, 143, 90-96.
https://doi.org/10.1016/j.neucom.2014.06.016
[23] Zhang, Z. and Cao, J. (2018) Novel Finite-Time Synchroniza-tion Criteria for Inertial Neural Networks with Time Delays via Integral Inequality Method. IEEE Transactions on Neural Networks and Learning Systems, 30, 1476-1485.
https://doi.org/10.1109/TNNLS.2018.2868800
[24] Yang, X., Cao, J. and Liang, J. (2016) Exponential Synchro-nization of Memristive Neural Networks with Delays: Interval Matrix Method. IEEE Transactions on Neural Networks and Learning Systems, 28, 1878-1888.
https://doi.org/10.1109/TNNLS.2016.2561298
[25] Yu, J., Hu, C., Jiang, H., et al. (2011) Exponential Synchroni-zation of Cohen-Grossberg Neural Networks via Periodically Intermittent Control. Neurocomputing, 74, 1776-1782.
https://doi.org/10.1016/j.neucom.2011.02.015
[26] Chinnathambi, R., Rihan, F.A. and Shanmugam, L. (2018) Sta-bilization of Delayed Cohen-Grossberg BAM Neural Networks. Mathematical Methods in the Applied Sciences, 41, 593-605.
https://doi.org/10.1002/mma.4630