一类带时滞的四元数神经网络的全局指数反同步
Global Exponential Anti-Synchronization for a Class of Quaternion-Valued Neural Networks with Time Delays

作者: 陈星宇 , 朱培勇 :电子科技大学数学科学学院,四川 成都;

关键词: 四元数神经网络指数反同步时滞Quaternion-Valued Neural Networks Exponential Anti-Synchronization Time Delays

摘要:
本文就一类带时滞的四元数值神经网络(QVNNs)的指数反同步问题进行研究。首先依据Hamilton法则,将QVNNs分解为4个等价的实值神经网络(RVNNs)。然后通过构造合适的Lyapunov-Krasovskii泛函,利用驱动–响应同步方法和矩阵不等式理论,得到了四元数值时滞神经网络关于指数反同步的一个结果。最后,通过数值仿真算例验证了所提出方法的有效性。

Abstract: In this paper, the problem of global exponential anti-synchronization for a class of quaterni-on-valued neural networks (QVNNs) with time delays is studied. Firstly, in light of Hamilton’s law, the QVNNs were decomposed into four real-valued parts to compose the equivalent real-valued sys-tems. Next, by constructing appropriate Lyapunov-Krasovskii functionals, utilizing drive-response synchronization method and matrix inequality theory, a result about global exponential an-ti-synchronization of QVNNs with time delays is realized. Finally, a numerical example was given to verify the effectiveness of the proposed method.

1. 引言

自1990年Pecora和Caroll [1] 引入驱动–响应同步概念提出混沌同步原理以来,因其在保密通信领域的应用价值,诸多学者对神经网络的同步性问题进行了研究,提出了多种同步方法 [2] [3] [4] [5] [6]。同时同步的概念也得到了拓宽,如完全同步 [7]、延迟同步 [8]、反同步 [9] [10] 等。

四元数是复数在某种意义上的扩展,在处理风速预报、彩色图像、阵列信号等实值神经网络和复值神经网络无法处理的问题上具有显著优势。此外由于四元数乘法不可交换的特性,四元数神经网络动力学特性比实值神经网络和复值神经网络更复杂。然而,文献 [2] - [10] 主要集中于实值神经网络之间的同步问题,因此,一个自然的问题是:四元数神经网络反同步是否有相应的结果?

本文就上述四元数神经网络问题进行研究,主要考虑一类具有时滞的四元数神经网络的指数反同步问题,利用李亚普洛夫理论和驱动–响应同步方法探讨,给出了四元数值时滞神经网络实现指数反同步的一个判定准则和反同步控制器设计方法,通过数值仿真验证了所提出方法的有效性。

2. 预备知识

2.1. 四元数代数

一个四元数 q 可被表示如下形式:

q = q 1 + q 2 i + q 3 j + q 4 k ,

其中 q 1 , q 2 , q 3 , q 4 R i , j , k 是虚数单位,且满足

i 2 = j 2 = k 2 = i j k = 1 , i j = j i = k , j k = k j = i , k i = i k = j ,

q 1 为四元数 q 的实部,记为 R e ( q ) = q 1 q 2 i + q 3 j + q 4 k 称为四元数 q 的虚部,记为 I m ( q ) = q 2 i + q 3 j + q 4 k 。由此我们可知四元数乘法不可交换。定义另外一个四元数

p = p 1 + p 2 i + p 3 j + p 4 k ,

则四元数 p , q 的和与乘积分别定义为

p + q = ( p 1 + q 1 ) + ( p 2 + q 2 ) i + ( p 3 + q 3 ) j + ( p 4 + q 4 ) k ,

p q = ( p 1 q 1 p 2 q 2 p 3 q 3 p 4 q 4 ) + ( p 1 q 2 + p 2 q 1 + p 3 q 4 p 4 q 3 ) i + ( p 1 q 3 + p 3 q 1 + p 4 q 2 p 2 q 4 ) j + ( p 1 q 4 + p 4 q 1 + p 2 q 3 p 3 q 2 ) k .

p = p ¯ = p 1 p 2 i p 3 j p 4 k 表示 p 的共轭转置,并且 p 的模表示为

| p | = p p * = p 1 2 + p 2 2 + p 3 2 + p 4 2 .

u = ( u 1 , u 2 , , u n ) T n ,则 u 的模可表示为 | u | = ( | u 1 | , | u 2 | , , | u n | ) T u 的范数表示为

u = i = 1 n | u i | 2 .

注:本文用 R m × n m × n m × n 分别表示 m × n 维实值、复值、四元数值矩阵, A B ( > ) 0 A B 是半正定(正定)矩阵。

2.2. 系统表示及基本引理

本文考虑如下一类四元数值时滞神经网络模型:

s ˙ ( t ) = A s ( t ) + B f ( s ( t ) ) + C g ( s ( t τ ) ) + J , (1)

其中 s ( t ) = ( s 1 ( t ) , , s n ( t ) ) T n 是神经元的状态变量,对角矩阵 A n × n 表示神经元的自反馈系数矩阵且 A > 0 ,矩阵 B , C n × n 分别表示神经元之间的与时滞无关、与时滞有关的权值矩阵, f ( s ( ) ) = ( f 1 ( s 1 ( ) ) , , f n ( s n ( ) ) ) T n g ( s ( ) ) = ( g 1 ( s 1 ( ) ) , , g n ( s n ( ) ) ) T n n 分别表示没有时滞和有时滞的激活函数, J = ( J 1 , , J n ) T n 表示系统外部输入, τ 是信号传输时滞且 τ 0

系统(1)的初始值是 s ( t ) = ϕ ( v ) ϕ ( v ) v [ τ , 0 ] 上有界且连续。

以系统(1)为驱动系统,设计相应的响应系统如下所示:

s ˜ ˙ ( t ) = A s ˜ ( t ) + B f ( s ˜ ( t ) ) + C g ( s ˜ ( t τ ) ) + J + U ( t ) , (2)

其中 s ˜ ( t ) = ( s ˜ 1 ( t ) , , s ˜ n ( t ) ) T n 表示相应系统神经元状态变量, U ( t ) = ( U 1 ( t ) , U 2 ( t ) , , U n ( t ) ) n 是实现驱动系统(1)和响应系统(2)指数反同步的控制器。

系统(2)的初始值是 s ˜ ( t ) = ϕ ˜ ( v ) ϕ ˜ ( v ) v [ τ , 0 ] 上有界且连续。

本文给出下列假设:

假设1 为方便计算,设激活函数 f ( ) g ( ) x Q n 可以分解成如下形式:

f ( x ) = f R ( x R ) + f I ( x I ) i + f J ( x J ) j + f K ( x K ) k ,

g ( x ) = g R ( x R ) + g I ( x I ) i + g J ( x J ) j + g K ( x K ) k ,

其中 f μ ( x μ ) : n n ( μ = R , I , J , K )

假设2 对 s 1 , s 2 n ,激活函数 f ( ) , g ( ) 满足Lipschitz条件,即存在常数 l p ( μ ) > 0 m p ( μ ) > 0 使得:

| f p ( μ ) ( s 1 ( μ ) ) + f p ( μ ) ( s 2 ( μ ) ) | l p ( μ ) | s 1 ( μ ) + s 2 ( μ ) | , | g p ( μ ) ( s 1 ( μ ) ) + g p ( μ ) ( s 2 ( μ ) ) | m p ( μ ) | s 1 ( μ ) + s 2 ( μ ) | , (3)

成立,其中 p = 1 , 2 , , n μ = r , i , j , k

引理1 [6] 假设存在正数 k 1 k 2 满足 k 1 > k 2 > 0 p ( t ) 是非负函数且满足

D + p ( t ) k 1 p ( t ) + k 2 p ˜ ( t ) ,

其中 p ˜ ( t ) = sup t τ s t { p ( s ) } ,则

p ( t ) p ˜ ( t 0 ) e λ ( t t 0 ) , t t 0 ,

成立,其中 λ 是方程 λ = k 1 k 2 e λ τ 的唯一解。

引理2 [11] 假设 Ω 1 Ω 2 Ω 3 均是任意维数的实矩阵,且 Ω 3 > 0 ,则对 x , y R n

2 x T Ω 1 T Ω 2 y x T Ω 1 T Ω 3 Ω 1 x + y T Ω 2 T Ω 3 1 Ω 2 y

成立。

定义系统反同步误差信号为 e ( t ) = s ( t ) + s ˜ ( t ) ,因此有

e ˙ ( t ) = A e ( t ) + B f ˜ ( e ( t ) ) + C g ˜ ( e ( t τ ( t ) ) ) + U ( t ) , (4)

其中 f ˜ ( e ( t ) ) = f ( s ˜ ( t ) ) + f ( s ( t ) ) , g ˜ ( e ( t τ ( t ) ) ) = g ( s ˜ ( t τ ( t ) ) ) + g ( s ( t τ ( t ) ) ) , e ( v ) = ϕ ˜ ( v ) + ϕ ( v ) 是误差系统(4)的初始值。

类比文献 [12] 中指数反同步的定义,引入下面定义:

定义1 如果存在常数 M > 0 κ > 0 使得误差满足下列不等式:

e ( t ) M sup τ s 0 e ( s ) exp ( κ t ) ,

则称驱动系统(1)和响应系统(2)是全局指数反同步的。

利用四元数运算法则,式(4)所示误差系统可以分解为如下等价的实值系统:

e ˙ R ( t ) = A e R + B R f ˜ R ( e R ( t ) ) B I f ˜ I ( e I ( t ) ) B J f ˜ J ( e J ( t ) ) B K f ˜ K ( e K ( t ) ) + C R g ˜ R ( e R ( t τ ) ) C I g ˜ I ( e I ( t τ ) ) C J g ˜ J ( e J ( t τ ) ) C K g ˜ K ( e K ( t τ ) ) + U R ( t ) , e ˙ I ( t ) = A e I + B R f ˜ I ( e I ( t ) ) + B I f ˜ R ( e R ( t ) ) + B J f ˜ K ( e K ( t ) ) B K f ˜ J ( e J ( t ) ) + C R g ˜ I ( e I ( t τ ) ) + C I g ˜ R ( e R ( t τ ) ) + C J g ˜ K ( e K ( t τ ) ) C K g ˜ J ( e J ( t τ ) ) + U I ( t ) , e ˙ J ( t ) = A e J + B R f ˜ J ( e J ( t ) ) + B J f ˜ R ( e R ( t ) ) B I f ˜ K ( e K ( t ) ) + B K f ˜ I ( e I ( t ) ) + C R g ˜ J ( e J ( t τ ) ) + C J g ˜ R ( e R ( t τ ) ) C I g ˜ K ( e K ( t τ ) ) + C K g ˜ I ( e I ( t τ ) ) + U J ( t ) , e ˙ K ( t ) = A e K + B R f ˜ K ( e K ( t ) ) + B K f ˜ R ( e R ( t ) ) + B I f ˜ J ( e J ( t ) ) B J f ˜ I ( e I ( t ) ) + C R g ˜ K ( e K ( t τ ) ) + C K g ˜ R ( e R ( t τ ) ) + C I g ˜ J ( e J ( t τ ) ) C J g ˜ I ( e I ( t τ ) ) + U K ( t ) . (5)

3. 主要结果

定理1 在假设1和假设2被满足的基础上,设计控制器的形式为:

U u ( t ) = ζ p u e p u ( t ) , (6)

其中 ζ p u > 0 ζ p u R p = 1 , 2 , , n u = R , I , J , K

如果存在矩阵 L ¯ u M ¯ u > 0 ( u = R , I , J , K ) 和控制增益 ζ p u 使得下式成立:

ζ u > λ max ( Ψ + 2 L ¯ u ) , min { ζ u λ max ( Ψ + 2 L ¯ u ) } > max { λ max ( M ¯ R ) , λ max ( M ¯ I ) , λ max ( M ¯ J ) , λ max ( M ¯ K ) } > 0 , (7)

Ψ = A + 1 2 { B R ( B R ) T + B I ( B I ) T + B J ( B J ) T + B K ( B K ) T + C R ( C R ) T + C I ( C I ) T + C J ( C J ) T + C K ( C K ) T } ,

其中 ζ R = min 1 p n { ζ p R } , ζ I = min 1 p n { ζ p I } , ζ J = min 1 p n { ζ p J } , ζ K = min 1 p n { ζ p K }

则在控制器(6)下驱动系统(1)和响应系统(2)是全局指数反同步的。

证明:构造Lyapunov函数形式如下:

V ( t ) = 1 2 ( ( e R ( t ) ) T e R ( t ) + ( e I ( t ) ) T e I ( t ) + ( e J ( t ) ) T e J ( t ) + ( e K ( t ) ) T e K ( t ) ) (8)

对(8)进行求导,其时间导数如下:

V ˙ ( t ) = [ e ˙ R ( t ) ] T e R ( t ) + [ e ˙ I ( t ) ] T e I ( t ) + [ e ˙ J ( t ) ] T e J ( t ) + [ e ˙ K ( t ) ] T e K ( t ) (9)

由假设2可知,一定存在对角矩阵 L ¯ , M ¯ 使得下式成立:

[ f ˜ ( μ ) ( e μ ( t ) ) ] T f ˜ ( μ ) ( e μ ( t ) ) [ e μ ( t ) ] T L ¯ μ e μ ( t ) , [ g ˜ ( μ ) ( e μ ( t τ ) ) ] T g ˜ ( μ ) ( e μ ( t τ ) ) [ e μ ( t τ ) ] T M ¯ μ e μ ( t τ ) , (10)

其中 L ¯ μ = d i a g ( ( l 1 μ ) 2 , , ( l n μ ) 2 ) M ¯ μ = d i a g ( ( m 1 μ ) 2 , , ( m n μ ) 2 ) μ = R , I , J , K

由引理2可得下列不等式

[ e α ( t ) ] T B β f ˜ γ ( e γ ( t ) ) 1 / 2 [ e α ( t ) ] T B β ( B β ) T e α ( t ) + 1 / 2 ( f ˜ γ ( e γ ( t ) ) ) T f ˜ γ ( e γ ( t ) ) , [ e α ( t ) ] T C β g ˜ γ ( e γ ( t τ ) ) 1 / 2 [ e α ( t ) ] T C β ( C β ) T e α ( t ) + 1 / 2 ( g ˜ γ ( e γ ( t τ ) ) ) T g ˜ γ ( e γ ( t τ ) ) (11)

成立,其中 α , β , γ = R , I , J , K

由式(9)~(11)有:

V . ( t ) [ e R ( t ) ] T { A ζ R + 2 L ¯ R + 1 2 { B R ( B R ) T + B I ( B I ) T + B J ( B J ) T + B K ( B K ) T + C R ( C R ) T + C I ( C I ) T + C J ( C J ) T + C K ( C K ) T } } e R ( t ) + [ e I ( t ) ] T { A ζ I + 2 L ¯ I + 1 2 { B R ( B R ) T + B I ( B I ) T + B J ( B J ) T + B K ( B K ) T + C R ( C R ) T + C I ( C I ) T + C J ( C J ) T + C K ( C K ) T } } e I ( t )

+ [ e J ( t ) ] T { A ζ J + 2 L ¯ J + 1 2 { B R ( B R ) T + B I ( B I ) T + B J ( B J ) T + B K ( B K ) T + C R ( C R ) T + C I ( C I ) T + C J ( C J ) T + C K ( C K ) T } } × e J ( t ) + [ e K ( t ) ] T { A ζ K + 2 L ¯ K + 1 2 { B R ( B R ) T + B I ( B I ) T + B J ( B J ) T + B K ( B K ) T + C R ( C R ) T + C I ( C I ) T + C J ( C J ) T + C K ( C K ) T } } e K ( t )

+ 2 { [ e R ( t τ ) ] T × M ¯ R e R ( t τ ) + [ e I ( t τ ) ] T M ¯ I e I ( t τ ) + [ e J ( t τ ) ] T M ¯ J e J ( t τ ) + [ e K ( t τ ) ] T × M ¯ K e K ( t τ ) } α V ( t ) + β V ( t τ ) . (12)

由引理1很容易得到:

V ( t ) max τ θ 0 V ( θ ) e κ t ,

| e R ( t ) | 2 + | e I ( t ) | 2 + | e J ( t ) | 2 + | e K ( t ) | 2 max τ θ 0 ( | e R ( θ ) | 2 + | e I ( θ ) | 2 + | e J ( θ ) | 2 + | e K ( θ ) | 2 ) e κ t ,

e ( t ) 2 max τ θ 0 ( e ( θ ) 2 ) e κ t ,

其中 κ 是方程 α β e κ τ κ = 0 的解。因此驱动系统(1)和响应系统(2)在控制器(6)下实现了全局指数反同步。

4. 仿真示例

例1 考虑如下时滞具有常时滞的四元数值神经网络:

s ˙ i ( t ) = A i j s j ( t ) + B i j f j ( s j ( t ) ) + C i j g j ( s j ( t τ j ) ) + J i , (13)

为了实现反同步,设计如下响应系统:

s ˙ i ( t ) = A i j s j ( t ) + B i j f j ( s j ( t ) ) + C i j g j ( s j ( t τ j ) ) + J i + u i ( t ) , (14)

其中

A = ( 0.035 0 0 0.09 ) ,

B = ( 1.5 0.3575 i 0.525 j 3.15 k 1.25 + 0.2145 i + 0.98 j + 3.6 k 1.25 + 0.3575 i + 0.035 j + 5.4 k 1.5 0.2145 i + 0.56 j 7.2 k ) ,

C = ( 0.75 + 0.2295 i + 0.27 j + 0.675 k 1.5 0.3825 i 0.27 j 0.6 k 0.315 + 0.612 i + 0.135 j 0.45 k 0.75 0.459 i 0.675 j 0.375 k ) ,

J = ( 0.2 0.15 i 0.14 j 0.12 k 0.5 + 0.19 i 0.3 j + 0.16 k ) ,

f 1 ( s 1 ) = 1.53 tan h ( s 1 R ) + 1.53 tan h ( s 1 I ) i + 1.53 tan h ( s 1 J ) j + 1.53 tan h ( s 1 K ) k , f 2 ( s 2 ) = 0.85 tan h ( s 2 R ) + 0.85 tan h ( s 2 I ) i + 0.85 tan h ( s 2 J ) j + 0.85 tan h ( s 2 K ) k ,

g 1 ( s 1 ) = 0.22 tan h ( s 1 R τ 1 ) + 0.22 tan h ( s 1 I τ 1 ) i + 0.22 tan h ( s 1 J τ 1 ) j + 0.22 tan h ( s 1 K τ 1 ) k , g 2 ( s 2 ) = 0.42 tan h ( s 2 R τ 2 ) + 0.42 tan h ( s 2 I τ 2 ) i + 0.42 tan h ( s 2 J τ 2 ) j + 0.42 tan h ( s 2 K τ 2 ) k .

对于 τ 1 = 0.1 τ 2 = 0.2 ,式(6)所示控制器参数取为 ζ i R = 52.10 , ζ i I = 50.16 , ζ i J = 50.16 , ζ i K = 57.36 ( i = 1 , 2 ) ,很容易验证其满足定理1条件,因此可以得出系统(13)和系统(14)在控制器下是指数反同步的结论。

随机选取10组数据。图1为系统运行时的误差变化曲线,显然驱动响应系统实现了反同步。

Figure 1. The anti-synchronization error curves between s i s ˜ i ( i = 1 , 2 )

图1. s i s ˜ i ( i = 1 , 2 ) 之间的反同步误差变化曲线

令驱动系统(13)和响应系统(14)的初始条件分别为 s 1 = 0.4 + 0.3 i + 0.5 j + 0.8 k s 2 = 0.2 + 0.2 i + 0.3 j + 1.1 k s ˜ 1 = 0.65 + 0.5 i 0.7 j + 0.1 k s ˜ 2 = 0.2 0.1 i + 0.2 j + 0.9 k ,系统(13)-(14)的状态变量的实部和虚部随时间的变化如图2所示。

Figure 2. Variation curves of the real and imaginary parts of a variable

图2. 变量的实部和虚部变化曲线

5. 结束语

利用Hamilton法则将具有常时滞的四元数值神经网络系统转化为等价的实值神经网络系统,设计一个简单易于实现的控制器,给出了驱动系统和响应系统的实现全局指数反同步的判别准则。最后通过数值仿真,验证了提出的设计方案的有效性,丰富了四元数值时滞神经网络动力学行为分析的相关研究。

文章引用: 陈星宇 , 朱培勇 (2021) 一类带时滞的四元数神经网络的全局指数反同步。 理论数学, 11, 1263-1270. doi: 10.12677/PM.2021.117140

参考文献

[1] Pecora, L.M. and Carroll, T.L. (1990) Synchronization in Chaotic Systems. Physical Review Letters, 64, 821-824.
https://doi.org/10.1103/PhysRevLett.64.821

[2] 张玮玮. 时滞混沌神经网络的指数同步[J]. 安庆师范学院学报:自然科学版, 2013, 19(3): 30-33.

[3] 任玲, 刘丽霞, 贝晓萌. 一类混沌系统的同步、反同步及其相关问题研究[J]. 齐鲁工业大学学报:自然科学版, 2015, 29(2): 56-60.

[4] 陈学菲, 刘辉昭. 参数不确定非自治混沌系统的自适应指数同步[J]. 应用数学和力学, 2021, 42(3): 316-322.

[5] 闫欢, 赵振江, 宋乾坤. 具有泄漏时滞的复值神经网络的全局同步性[J]. 应用数学和力学, 2016, 37(8): 832-841.

[6] Gopalsamy, K. (1992) Stability and Oscillations in Delay Differential Equations of Population Dynamics. Kluwer, Dordrecht, The Netherlands.
https://doi.org/10.1007/978-94-015-7920-9

[7] 范云, 张荣, 过榴晓. 时变耦合网络的完全同步[J]. 江南大学学报(自然科学版), 2011, 10(3): 258-262.

[8] 贾飞蕾, 徐伟. 一类参数不确定混沌系统的延迟同步[J]. 物理学报, 2007, 56(6): 3101-3106.

[9] 吕翎, 夏晓岚. 非线性耦合时空混沌系统的反同步研究[J]. 物理学报, 2009, 58(2): 814-818.

[10] Wu, Z., Wang, Z., Zhou, T. and Tan, S. (2020) Global Synchronization and Anti-Synchronization of Fractional-Order Complex-Valued Gene Regulatory Networks with Time-Varying Delays. IEEE Access, 8, 150555-150572.
https://doi.org/10.1109/ACCESS.2020.3016706

[11] 吴桂华, 廖晓峰. 脉冲时滞Hopfield神经网络的全局指数稳定性[J]. 计算机技术与发展, 2009, 19(9): 25-27.

[12] 王兰梓, 房辉. 混沌系统异结构之间的指数反同步[J]. 西南民族大学学报(自然科学版), 2010, 36(5): 721-726.

分享
Top