一类分数阶中立型时滞神经网络的稳定性
Stability of a Class of Fractional-Order Neutral Time-Delay Neural Networks
摘要: 本文考虑了一类具有时变时滞的分数阶中立型神经网络模型。根据同胚映射原理证明了该系统平衡点的存在唯一性。然后通过构造Lyapunov泛函以及利用Lyapunov稳定性理论,证明了系统渐近稳定的充分条件。
Abstract: In this paper, a class of fractional-order neutral neural network models with time-varying delays was considered. The uniqueness of the existence of the equilibrium point of the system is proved by the principle of homeomorphic mapping. Then, by constructing Lyapunov functional and using Lyapunov stability theory, the sufficient conditions for asymptotic stability of the system are proved.
文章引用:曹行. 一类分数阶中立型时滞神经网络的稳定性[J]. 理论数学, 2025, 15(1): 63-68. https://doi.org/10.12677/pm.2025.151007

1. 引言

近年来,随着科技的飞速发展,神经网络模型已经广泛应用在很多智能计算的学科领域上[1] [2],许多学者指出,分数阶微积分非常适合描述各种材料和过程的记忆性和遗传性[3]。在此基础上,越来越多的学者关注分数阶神经网络模型,研究分析其动力学行为[4] [5]。稳定性问题是动力系统的首要问题。因此,分数阶神经网络系统的稳定性是一个重要的研究领域,近年来取得了一系列有意义的结果[6]-[12]

目前,关于分数阶神经网络的稳定性分析主要涉及有限时间稳定、全局渐近稳定、Mittag-Leffler稳定等。其中,文献[9]利用分数阶Gronwall不等式,研究了分数阶神经网络的有限时间稳定性问题。文献[10]利用比较原理、不动点定理,研究了在Caputo导数意义下,一类分数阶神经网络的渐近稳定性问题。文献[11]运用Lyapunov泛函方法,研究了Riemann-Liouville分数阶神经网络的渐近稳定性。最近,文献[12]利用同胚映射原理、Lyapunov泛函方法,研究了一类整数阶中立型时滞神经网络模型的渐近稳定性。文献[13]利用压缩映像原理、Lyapunov泛函方法,研究了一类具有分布时滞的分数阶神经网络模型的渐近稳定性。

受上述讨论的启发,我们将一类整数阶神经网络推广到分数阶情形,提出了一类Riemann-Liouville分数阶中立型时滞神经网络,利用同胚映射定理证明了平衡点的存在唯一性,得到了该分数阶神经网络渐近稳定性的一个充分条件。同胚映射定理是常用的、巧妙的证明平衡点的存在唯一性的方法之一。

本文的主要结构如下:在第2节中介绍了Riemann-Liouville导数的定义、性质和一些重要结论。在第3节中,研究证明了平衡点的存在唯一性,同时得到渐近稳定性的一个充分条件。在第4节中,通过一个数值例子来验证主要结果的有效性。最后,给出全文的主要结论。

2. 模型描述与预备知识

2.1. 模型描述

本文考虑了如下分数阶中立型神经网络系统:

D t 0 RL t α x i ( t )= a i x i ( t )+ j=1 n b ij f j ( x j ( t ) )+ j=1 n c ij g j ( x j ( tτ( t ) ) )+ j=1 n d ij ( D t 0 RL t α x i ( tδ ) ) + U i (1)

其中, D t 0 RL t α x i ( t ) 表示Riemann-Liouville分数阶导数, i=1,2,n,0<α<1 x i ( t ) 表示第 i 个神经元在 t 时刻的状态, n 为神经元的个数, A=diag( a 1 , a 2 ,, a n ) ,常数 a i >0 是神经元的自调节参数, τ( t ) 是连续的有界时变时滞函数, 0<τ( t ) τ 0 <1 ,连接权重矩阵 B= ( b ij ) m×n ,C= ( c ij ) m×n ,D=diag( d 1 , d 2 ,, d n ) δ 为常数时滞, U i 表示恒定的外部输入, f,g 表示激活函数,系统(1)的初始条件为: D t 0 RL t α1 x i ( t )= ω i ( t ),t[ τ 0 ,0 ]

2.2. 预备知识

定义1 函数 y(t) p 阶Riemann-Liouville分数阶积分为 D t 0 RL t p y( t )= 1 Γ( p ) t 0 t ( ts ) p1 y( s ) ds,p>0 ,其中 Γ( p ) 表示 Gamma 函数,且 Γ( p )= t 0 t p1 e t dt

定义2 函数 y(t) q 阶Riemann-Liouville分数阶导数为:

D t 0 RL t q y( t )= 1 Γ( nq ) d n d t n t 0 t ( ts ) nq1 y( s ) ds,n1q<n Z +

本文给出如下假设与引理:

A1 :假设激活函数 f( x ) g( x ) 满足 Lipschitz 条件,即存在 l i 1 >0,s.t. f i ( u ) f i ( v ) l i 1 | uv | l j 2 >0,s.t. g j ( u ) g j ( v ) l j 2 | uv | 对任意的 u,vR,i=1,2,n

引理1 [12] 任意向量 x= ( x 1 , x 2 ,, x n ) T ,y= ( y 1 , y 2 , y n ) T 以及任意的常数 p( p>0 ) ,则下列不等式成立: ±2 x T yp x T x+ p 1 y T y

引理2 [12] 对任意实数矩阵 A= ( a ij ) n×n 和向量 x= ( x 1 , x 2 ,, x n ) T ,有: x T A T Ax i=1 n a i ~ x i 2 ,其中 a i ~ = j=1 n | k=1 n a ki a kj | ,i=1,2,,n

引理3 [12] f(x) 为连续函数,且 f R n 上的单射,若 x + 时, f( x ) + ,则 f( x ) R n 上的同胚映射。

引理4 [13] α>β>0 ,有以下等式成立: D t 0 RL t α ( D t 0 RL t β y( t ) )= D t 0 RL t αβ y( t )

引理5 [13] y( t ) R n 是一个可微的向量函数,则有以下不等式:

D t 0 RL t α ( y T ( t )Py( t ) )2 y T ( t )P D t 0 RL t α y( t )

3. 主要结果

3.1. 平衡点的存在唯一性

定理1 若激活函数满足条件假设 A1 ,且存在正数 p 1 p 2 ,使得 φ=ψ L 2 γ>0 ,则系统(1)存在唯一的平衡点,其中:

ψ=diag( ζ 1 , ζ 2 ,, ζ n ), ζ i =( 2 p 1 p 2 ) d i 2 ,γ=diag( γ 1 , γ 2 ,, γ n ), γ i = p 1 1 b i + p 2 1 c i

b i = j=1 n | k=1 n b ki b kj | , c i = j=1 n | k=1 n c ki c kj |, i=1,2,,n

证明:对于系统(1),写成向量形式为

D t 0 RL t α x( t )=Ax( t )+Bf( x( t ) )+Cg( x( tτ( t ) ) )+D( D t 0 RL t α x( tδ ) )+U (2)

假设 x 0 = ( x 1 , x 2 ,, x n ) T 是系统(1)的一个平衡点,则有

A x 0 ( t )+Bf( x 0 ( t ) )+Cg( x 0 ( tτ( t ) ) )+U=0 (3)

定义映射: ϕ( x )=Ax+Bf( x )+Cf( x )+U

下面证明 ϕ(x) 是一个单射。对任意的 x,y R n ,且 xy

那么: ϕ( x )ϕ( y )=A( xy )+B(f( ( x )f( y ) )+C( f( x )f( y ) )

在上式两边左乘 2 ( xy ) T A 得到:

2 ( xy ) T A( ϕ( x )ϕ( y ) )=2 ( xy ) T AA( xy ) +2 ( xy ) T AB( f( x )f( y ) )+2 ( xy ) T AC( f( x )f( y ) ) (4)

由引理1可知:

2 ( xy ) T AA( xy )=2 ( xy ) T A 2 ( xy ) (5)

2 ( xy ) T AB( f( x )f( y ) ) p 1 ( xy ) T A 2 ( xy )+ p 1 1 ( f( x )f( y ) ) T B T B( f( x )f( y ) ) (6)

2 ( xy ) T AC( f( x )f( y ) ) p 2 ( xy ) T A 2 ( xy )+ p 2 1 ( f( x )f( y ) ) T C T C( f( x )f( y ) ) (7)

将(5)~(7)式代入(4)式得到:

2 ( xy ) T A( ϕ( x )ϕ( y ) )2 ( xy ) T A 2 ( xy )+( p 1 + p 2 ) ( xy ) T A 2 ( xy )+

p 1 1 ( f( x )f( y ) ) T B T B( f( x )f( y ) )+ p 2 1 ( f( x )f( y ) ) T C T C( f( x )f( y ) )

由引理2与假设 A1 可得:

2 ( xy ) T A( ϕ( x )ϕ( y ) )2 ( xy ) T A 2 ( xy )+( p 1 + p 2 ) ( xy ) T A 2 ( xy )+

p 1 1 i=1 n b i ~ l i 2 ( x i y i ) 2 + p 2 1 i=1 n c i ~ l i 2 ( x i y i ) 2

2 ( xy ) T A( ϕ( x )ϕ( y ) ) i=1 n [ ( 2 p 1 p 2 ) a i 2 p 1 1 l i 2 b i ~ p 2 1 l i 2 c i ~ ] ( x i y i ) 2

i=1 n [ ( 2 p 1 p 2 ) a i 2 p 1 1 l i 2 b i ~ p 2 1 l i 2 c i ~ ] ( x i y i ) 2 = i=1 n [ ( 2 p 1 p 2 ) a i 2 l i 2 ( p 1 1 b i ~ + p 2 1 c i ~ ) ] ( x i y i ) 2

= i=1 n ( ζ i ~ l i 2 γ i ) ( x i y i ) 2 = ( xy ) T ( ψ L 2 γ )( xy )= ( xy ) T φ( xy )

所以, 2 ( xy ) T A( ϕ( x )ϕ( y ) ) ( xy ) T φ( xy ) ,由于 φ=ψ L 2 γ>0 ,那么对于任意的 xy ,都有 ϕ( x )ϕ( y ) ,即 ϕ( x ) 是一个单射。

另外,在上述证明过程中令 y=0 ,则有: 2 x T A( ϕ( x )ϕ( 0 ) ) x T φx λ min ( φ ) x T x λ min ( φ ) x 2 ,由 Schwarz 不等式可知:

2 x T A ϕ( x )ϕ(0) λ min ( φ ) x 2

x 0 时, ϕ( x )ϕ( 0 ) λ min ( φ ) x 2 A

那么 x + 时, ϕ( x ) + ,所以 ϕ( x ) 是一个同胚映射。

由上述证明与引理3可得, ϕ(x) 是同胚映射,由同胚映射定理知,系统(1)存在唯一的平衡点。

3.2. 渐近稳定性

定理2 假设条件 A1 与定理1成立,若存在正定矩阵 P Q R ,使得以下不等式成立,

Ω=[ 2PA+ L 2 T Q L 2 ( PB+ L 1 T P L 1 B ) PC PD * 2P L 1 A( L 1 1 ) P L 1 C P L 1 D * * Q 0 * * * R ]<0 (8)

则系统(1)的解是渐近稳定的。其中, L 1 =diag[ l 1 1 , l 2 1 ,, l n 1 ], L 2 =diag[ l 1 2 , l 2 2 ,, l n 2 ]

证明:构造 Lyapunov 泛函 V( t )= i=1 4 V i ( t ) ,其中:

V 1 ( t )= D t 0 RL t α1 x T ( t )Px( t ) V 2 ( t )= D t 0 RL t α1 f T ( x( t ) )Pf( x( t ) )

V 3 ( t )= tτ(t) t g T ( x( s ) )Qg( x( s ) )ds V 4 ( t )= δ t ( D t 0 RL t α x( tδ ) ) T R D t 0 RL t α x( tδ ) dt

V( t )= i=1 4 V i ( t ) 求导,并且由引理4、5可得:

V 1 ( t )= D t 0 RL t α x T ( t )Px( t )2 x T ( t )P D t 0 RL t α x( t )

=2 x T ( t )PAx( t )+2 x T ( t )PBf( x( t ) )+2 x T ( t )PCg( x( tτ( t ) ) )+2 x T ( t )PD D t 0 RL t α x( tδ ) (9)

V 2 ( t )= D t 0 RL t α f T ( x( t ) )Pf( x( t ) )2 f T ( x( t ) )P D t 0 RL t α f( x( t ) )

2 f T ( x( t ) )P L 1 Ax( t )+2 f T ( x( t ) )P L 1 Bf( x( t ) )+ 2 f T ( x( t ) )P L 1 Cg( x( tτ( t ) ) )+2 f T ( x( t ) )P L 1 D D t 0 RL t α x( tδ ) (10)

V 3 ( t )= g T ( x( t ) )Qg( x( t ) ) g T ( x( tτ( t ) ) )Qg( x( tτ( t ) ) ) (11)

V 4 ( t )= ( D t 0 RL t α x( tδ ) ) T R D t 0 RL t α x( tδ ) (12)

根据上述(9)~(12)式有:

V ( t )= V 1 ( t )+ V 2 ( t )+ V 3 ( t )+ V 4 ( t ) 2 x T ( t )PAx( t )+2 x T ( t )PBf( x( t ) )+2 x T ( t )PCg( x( tτ( t ) ) )+2 x T ( t )PD D t 0 RL t α x( tδ ) 2 f T ( x( t ) )P L 1 Ax( t )+2 f T ( x( t ) )P L 1 Bf( x( t ) )+2 f T ( x( t ) )P L 1 Cg( x( tτ( t ) ) )+2 f T ( x( t ) )P L 1 D D t 0 RL t α x( tδ )+ g T ( x( t ) )Qg( x( t ) ) g T ( x( tτ( t ) ) )Qg( x( tτ( t ) ) ) ( D t 0 RL t α x( tδ ) ) T R D t 0 RL t α x( tδ )

由假设 A1 激活函数 f( x ) g( x ) 满足 Lipschitz 条件,所以有:

V ( t )= V 1 ( t )+ V 2 ( t )+ V 3 ( t )+ V 4 ( t ) x T ( t )( 2PA L 2 T Q L 2 )x( t )+ x T ( t )( 2PB+2 L 1 T P L 1 B )f( x( t ) )+ x T ( t )( 2PC )g( x( tτ( t ) ) )+ x T ( t )( 2PD ) D t 0 RL t α x( tδ )2 f T ( x( t ) )P L 1 A( L 1 1 )f( x( t ) )+ f T ( x( t ) )( 2P L 1 C )g( x( tτ( t ) ) )+ f T ( x( t ) )( 2P L 1 D ) D t 0 RL t α x( tδ ) g T ( x( tτ( t ) ) )Qg( x( tτ( t ) ) ) ( D t 0 RL t α x( tδ ) ) T R D t 0 RL t α x( tδ )

= [ x( t ) f( x( t ) ) g( x( tτ( t ) ) ) D t 0 RL t α x( tδ ) ] T [ 2PA+ L 2 T Q L 2 ( PB+ L 1 T P L 1 B ) PC PD * 2P L 1 A( L 1 1 ) P L 1 C P L 1 D * * Q 0 * * * R ][ x( t ) f( x( t ) ) g( x( tτ( t ) ) ) D t 0 RL t α x( tδ ) ]

根据公式(8),可以得到 V ( t ) <0 ,因此系统(1)的解是渐近稳定的。

4. 数值例子

例1考虑如下的二维中立型神经网络模型:

D t 0 RL t α x( t )=Ax( t )+Bf( x( t ) )+Cg( x( tτ( t ) ) )+D( D t 0 RL t α x( tδ ) )+U (13)

其中, α( 0,1 ),x( t )= ( cos( x 1 ( t ) ),cos( x 2 ( t ) ) ) T U= ( 0,0 ) T

f( x( t ) )=(2cos( x 1 ( t ) ),2cos ( x 2 ( t ) ) T g( x( tτ( t ) ) )= ( 1 2 cos( x 1 ( tτ( t ) ) ), 1 2 cos( x 2 ( tτ( t ) ) ) ) T

A=[ 3 1 0 3 ],B=[ 1 1 1 1 ],C=[ 1 1 1 1 ],D=[ 2 0 0 2 ]

l i 1 = l j 2 =1 ( i=1,2,j=1,2 ) ,那么 L 1 =[ 1 0 0 1 ], L 2 =[ 1 0 0 1 ]

LMI 得可行解:

P=[ 0.2027 0.0148 0.0148 0.2321 ],Q=[ 0.5466 0.0695 0.0695 0.5723 ],R=[ 1.0392 0.0305 0.0305 1.0459 ]

根据定理2可知,神经网络模型(13)是渐近稳定的。

5. 结语

本文研究一类分数阶时滞中立型神经网络模型。利用同胚映射原理证明平衡点的存在唯一性问题。结合Lyapunov稳定性方法,得到神经网络模型渐近稳定的条件,例子表明结果是可行有效的。相对于其他神经网络模型的研究,本文的研究突出于对分数阶中立型神经网络模型的分析,对于具有时变时滞的分数阶中立型神经网络而言,其他的动力学性质仍值得我们分析。

参考文献

[1] 白琮, 黄玲, 陈佳楠, 等. 面向大规模图像分类的深度卷积神经网络优化[J]. 软件学报, 2018, 29(4): 1029-1038.
[2] 刘涛. 基于粒计算-神经网络的汽车电子控制系统故障分析[J]. 佳木斯大学学报(自然科学版), 2021, 39(4): 59-67.
[3] Podlubny, I. (1999) Fractional Differential Equations. Academic Press, 62-75.
[4] Pahnehkolaei, S.M.A., Alfi, A. and Tenreiro Machado, J.A. (2017) Dynamic Stability Analysis of Fractional Order Leaky Integrator Echo State Neural Networks. Communications in Nonlinear Science and Numerical Simulation, 47, 328-337.
https://doi.org/10.1016/j.cnsns.2016.11.013
[5] 张欣欣. Caputo型分数阶神经网络的稳定性分析[D]: [硕士学位论文]. 秦皇岛: 燕山大学, 2015.
[6] Du, F. and Lu, J. (2020) Finite-time Stability of Neutral Fractional Order Time Delay Systems with Lipschitz Nonlinearities. Applied Mathematics and Computation, 375, Article ID: 125079.
https://doi.org/10.1016/j.amc.2020.125079
[7] Wang, Q., Lu, D. and Fang, Y. (2015) Stability Analysis of Impulsive Fractional Differential Systems with Delay. Applied Mathematics Letters, 40, 1-6.
https://doi.org/10.1016/j.aml.2014.08.017
[8] 韩小花. 分数阶微分方程的若干稳定性问题[D]: [硕士学位论文]. 南宁: 广西民族大学, 2016.
[9] Du, F. and Lu, J. (2021) New Criteria for Finite-Time Stability of Fractional Order Memristor-Based Neural Networks with Time Delays. Neurocomputing, 421, 349-359.
https://doi.org/10.1016/j.neucom.2020.09.039
[10] 王利敏, 宋乾坤, 赵振江. 基于忆阻的分数阶时滞复值神经网络的全局渐近稳定性[J]. 应用数学和力学, 2017, 38(3): 333-346.
[11] Zhang, H., Ye, R., Cao, J., Ahmed, A., Li, X. and Wan, Y. (2017) Lyapunov Functional Approach to Stability Analysis of Riemann‐Liouville Fractional Neural Networks with Time‐Varying Delays. Asian Journal of Control, 20, 1938-1951.
https://doi.org/10.1002/asjc.1675
[12] 田宝单, 钱紫舒, 徐鑫. 一类中立型时滞神经网络的稳定性分析[J]. 海南大学学报自然科学版, 2024, 42(3): 244-250.
[13] 刘雪铃, 黄静. 分数阶时滞神经网络渐近稳定与同步分析[J]. 淮北师范大学学报自然科学版, 2023, 44(1): 16-20.