Pure Mathematics
Vol. 11  No. 01 ( 2021 ), Article ID: 39932 , 11 pages
10.12677/PM.2021.111012

改进的一重积分不等式在时滞神经网络 中的应用

柳莹莹,熊良林*,邹梅,蔡丽,王珍

云南民族大学数学与计算机科学学院,云南 昆明

收稿日期:2020年12月14日;录用日期:2021年1月15日;发布日期:2021年1月22日

摘要

本文主要研究一类一重积分不等式的改进并将其用于分析时滞神经网络稳定性的问题。首先,借助交互凸组合不等式和多个零等式,改进一类一重积分不等式。其次,根据研究的时变时滞系统的属性,构造一类新的含有更多时滞信息的增广型Lyapunov泛函。然后,利用改进的积分不等式及其它分析技巧,获得保守性较低的稳定性判据。最后,通过一个数值例子验证所得结果的有效性和优越性。

关键词

交互凸组合不等式,增广的Lyapunov泛函,零等式,改进的一重积分不等式

The Application of Improved Single Integral Inequality in Neural Networks with Time-Delayed

Yingying Liu, Lianglin Xiong*, Mei Zou, Li Cai, Zhen Wang

School of Mathematics and Computer Science, Yunnan Minzu University, Kunming Yunnan

Received: Dec. 14th, 2020; accepted: Jan. 15th, 2021; published: Jan. 22nd, 2021

ABSTRACT

This paper is concerned with an improved single integral inequality with application in the time-delayed neural network. First of all, an improved single integral inequality is proposed by virtue of the reciprocally convex combination inequality and zero-qualities. Secondly, relay on the characteristics of the system more information about time-delay is considered in constructing the augmented Lyapunov functional. Then, a less conservative stability condition is obtained through utilizing the improved integral inequality and other analysis techniques. In the end, to verify the effectiveness and superiority of the derived results, a numerical example is provided.

Keywords:Reciprocally Convex Combination Inequality, Augmented Lyapunov Functional, Zero-Qualities, An Improved Single Integral Inequality

Copyright © 2021 by author(s) and Hans Publishers Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

1. 引言

时滞神经网络的稳定性问题是神经网络处理联想记忆、信号处理和一些优化问题等多个实际问题的基础 [1] [2] [3]。一方面,在应用神经网络之前,应确保稳定性;另一方面,在神经网络的实现中,由于神经元的有限切换速率就不可避免的出现时滞,很多相关研究表明时滞会导致不良的动态行为,这些行为会使得系统的性能下降甚至不稳定 [4],因此研究时滞神经网络的稳定性是非常有意义的。

到目前为止,学者们针对时滞神经网络系统的稳定性研究主要采用Lyapunov泛函方法。通过利用Lyapunov稳定性理论 [5],构造恰当的Lyapunov泛函以及尽可能精细地估计泛函的导数来得到低保守性的稳定性准则。对于构造Lyapunov泛函方法,目前比较常见的有构造带有增广的Lyapunov泛函 [6],利用时滞分割方法构造泛函 [7] 以及构造带有多重积分项的泛函 [8]。对于估计泛函的导数,主要是估计求导后产生的积分项,力求获得求导后更严格的上界。为了处理这些积分,学者们主要通过不等式放缩的方法,例如Jensen不等式 [9],基于Wirtinger的积分不等式 [10],基于自由矩阵的积分不等式 [11],改进的自由矩阵积分不等式 [12] 以及互凸不等式 [13] 等等。这些不等式有效地降低了时滞神经网络的时滞依赖稳定性条件,不过仍然有改进的空间。

文献 [14] 中通过引入增广的Lyapunov泛函并应用于Wirtinger不等式中,得到了保守性较低的时滞相关稳定性标准,文献 [15] 利用零等式改进具有增广的一重积分不等式,与直接使用Jensen不等式放缩相比,运用零等式可以巧妙加入自由矩阵,从而提供更大的自由度进而得到更紧的时滞上界。然而,文献 [15] 中不等式的放缩只带有一重积分信息,具有一定的保守性。众所周知,高阶积分项在降低时滞系统的保守性中扮演着十分重要的作用,例如在估计一重积分不等式中引入二重积分项,这时二重积分状态信息会与一重状态信息以及系统的状态信息交叉,从而得到更大的时滞上界,因此建立一个新的带有自由矩阵和二重积分项的具有增广的一重积分不等式对降低时滞神经网络系统的保守性具有深远的意义。此外,值得一提的是在增广向量中引入s相关项能有效克服条件的保守性,因为它将一重积分放入增广向量中并把这个增广向量放入带有一重积分的Lyapunov泛函中,在对这个泛函求导后能够获得更多状态信息。因此,本文拟通过新的积分不等式处理技巧结合更符合系统的Lyapunov泛函,获得保守性更低的稳定性条件。

为了简化表述,我们有必要做如下的符号说明: R n 表示n维欧式空间; R n × m 代表 n × m 的矩阵; sym { X } = X + X T diag { } 表示对角矩阵; col { } 列向量 X T ,和 X 1 分别矩阵的转置和逆矩阵; * 表示矩阵的对称项。

2. 预备知识

2.1. 系统描述

考虑一般的带有时变时滞的神经网络系统

{ x ˙ ( t ) = A x ( t ) + W 1 f ( W 0 x ( t ) ) + W 2 f ( W 0 x ( t h ( t ) ) ) , x ( θ ) = ϕ ( θ ) , θ [ h , 0 ] . (1)

其中 x ( t ) R n 是系统的状态向量 A R n × n W 0 R n × n W 1 R n × n ,和 W 2 R n × n 是已知的实数矩阵,初始条件 φ ( θ ) 是连续的微分函数,h和u是已知的常数,且满足如下条件:

0 < h ( t ) < h , h ˙ ( t ) < u . (2)

f ( W 0 x ( t ) ) = col { f ( W 01 x ( t ) ) , f ( W 02 x ( t ) ) , , f ( W 0 n x ( t ) ) } 是神经元的激活函数 W 0 i W 0 的第i行,假设 f ( 0 ) = 0 且满足Lipschitz条件

l i f i ( a ) f i ( b ) a b l i + , (3)

其中, l i l i + 是一些常数,且 L = diag { l 1 , l 2 , , l n } L + = diag { l 1 + , l 2 + , , l n + }

2.2. 引理

在分析系统(1)的时滞依赖稳定性判据之前,本节首先引入如下一些重要的引理:

引理1. [16] 给定一个二次函数 f ( y ) = a 2 y 2 + a 1 y + a 0 ,当 β [ 0 , 1 ] 时,如果有

T i = f ( h i ) < 0 ( i = 1 , 2 ) , T 3 = β 2 h 12 2 a 2 + f ( h 1 ) < 0 , T 4 = ( 1 β ) 2 h 12 2 a 2 + f ( h 2 ) < 0 ,

那么 f ( y ) < 0 ( h 1 y h 2 ) 成立,其中 h 12 = h 2 h 1

引理2. [12] 设 w ( t ) 是一个可微函数: [ a , b ] R n ,那么存在对称正定矩阵 H R n × n ,对于任意具有恰当维度的矩阵L,M,有下列的不等式成立:

a b w ˙ T ( s ) H w ˙ ( s ) d s sym { σ 0 T L σ 1 + σ 0 T M σ 2 } ( b a ) σ 0 T ( 3 L H 1 L T + M H 1 M T 3 ) σ 0

其中 σ 1 = x ( b ) x ( a ) σ 2 = x ( b ) + x ( a ) 2 b a a b x ( s ) d s σ 0 是任意向量。

引理3. [17] 若对称正定矩阵 R > 0 x ( t ) 是一个可微函数: [ a , b ] R n ,则下列不等式成立

a b x T ( s ) R x ( s ) d s 1 b a ( a b x ( s ) d s ) T R ( a b x ( s ) d s ) + 3 b a Ω 1 T R Ω 1 ,

a b x ˙ T ( s ) R x ˙ ( s ) d s 1 b a Ω 2 T R Ω 2 + 3 b a Ω 3 T R Ω 3 ,

其中

Ω 1 = a b x ( s ) d s 2 b a a b β b x ( s ) d s d β , Ω 2 = x ( b ) x ( a ) , Ω 3 = x ( b ) + x ( a ) 2 b a a b x ( s ) d s .

为了能得到保守性更低的稳定性条件,本文首先引入一个改进的一重积分不等式如下:

引理4. 若半正定对称矩阵 Y i R n × n ( i = 1 , 3 ) X i R n × n ( i = 1 , 2 ) ,任意矩阵 Y 2 R n × n X 3 R n × n X 4 R 2 n × 2 n x ( t ) 是一个可微函数: [ a , b ] R n ,满足:

Ω 0 = [ Y 1 Y 2 * Y 3 ] 0 , Ω 1 = [ Y 1 Y 2 + X 1 * Y 3 ] 0 , Ω 2 = [ Y 1 Y 2 + X 2 * Y 3 ] 0 ,

Ω 3 = [ Y 3 X 3 * Y 3 ] 0 , Ω 4 = [ Ω 1 X 4 * Ω 2 ] 0 ,

那么,对任意向量 c [ a , b ] ,有下面的不等式成立

Φ = a b W T ( s ) Ω 0 W ( s ) d s ξ T Ω ξ . (4)

其中

W ( s ) = col { x ( s ) , x ˙ ( s ) } , ξ = col { x ( b ) , x ( c ) , x ( a ) , w 1 , w 2 , v 1 , v 2 } , w 1 = 1 b c c b x ( s ) d s ,

w 2 = 1 c a a c x ( s ) d s , v 1 = 1 b c c b θ b x ( s ) d s d θ , v 2 = 1 c a a c θ c x ( s ) d s d θ ,

e i = [ 0 n , ( i 1 ) n I n 0 n , ( 7 i ) n ] ( i = 1 , 2 , , 7 ) , ζ 1 = col { e 1 e 2 , e 2 e 3 } ,

ζ 2 = col { ( b c ) e 4 2 e 6 , 2 e 4 e 1 e 2 , ( c a ) e 5 2 e 7 , 2 e 5 e 2 e 3 } ,

Ω = sym { e 4 T ( Y 2 + X 1 ) ( e 1 e 2 ) e 5 T ( Y 2 + X 2 ) ( e 2 e 3 ) } ( b c ) e 4 T Y 1 e 4 ( c a ) e 5 T Y 1 e 5 e 2 T ( X 1 X 2 ) e 2 + e 1 T X 1 e 1 e 3 T X 2 e 3 1 b a ζ 1 T Ω 3 ζ 1 1 b a ζ 2 T Ω 4 ζ 2 .

证明:对于半正定对称矩阵 X i R n × n ( i = 1 , 2 ) ,有下面的零等式成立

0 = 2 c b x T ( s ) X 1 x ˙ ( s ) d s + x T ( b ) X 1 x ( b ) x T ( c ) X 1 x ( c ) ,

0 = 2 a c x T ( s ) X 2 x ˙ ( s ) d s + x T ( c ) X 2 x ( c ) x T ( a ) X 2 x ( a ) .

将上面两个零等式加入二重积分,合并同类项后可以得到

Φ = a b W T ( s ) Ω 0 W ( s ) d s 2 c b x T ( s ) X 1 x ˙ ( s ) d s 2 a c x T ( s ) X 2 x ˙ ( s ) d s x T ( c ) ( X 1 X 2 ) x ( c ) + x T ( b ) X 1 x ( b ) x T ( a ) X 2 x ( a ) = c b W T ( s ) Ω 1 W ( s ) d s a c W T ( s ) Ω 2 W ( s ) d s x T ( c ) ( X 1 X 2 ) x ( c ) + x T ( b ) X 1 x ( b ) x T ( a ) X 2 x ( a )

利用引理3中的不等式放缩,展开计算如下

Φ 1 b c ( c b W ( s ) d s ) T Ω 1 ( c b W ( s ) d s ) 1 c a ( a c W ( s ) d s ) T Ω 2 ( a c W ( s ) d s ) 3 b c ( c b W ( s ) d s 2 b c c b β b W ( s ) d s d β ) T Ω 1 ( c b W ( s ) d s 2 b c c b β b W ( s ) d s d β ) 3 c a ( a c W ( s ) d s 2 c a a c β c W ( s ) d s d β ) T Ω 2 ( a c W ( s ) d s 2 c a a c β c W ( s ) d s d β ) x T ( c ) ( X 1 X 2 ) x ( c ) + x T ( b ) X 1 x ( b ) x T ( a ) X 2 x ( a )

= ( b c ) w 1 T Y 1 w 1 2 w 1 T ( Y 2 + X 1 ) ( x ( b ) x ( c ) ) 1 b c ( x ( b ) x ( c ) ) T Y 3 ( x ( b ) x ( c ) ) ( c a ) w 2 T Y 1 w 2 2 w 2 T ( Y 2 + X 2 ) ( x ( c ) x ( a ) ) 1 c a ( x ( c ) x ( a ) ) T Y 3 ( x ( c ) x ( a ) ) 3 b c [ ( b c ) w 1 2 v 1 2 w 1 x ( b ) x ( c ) ] T [ Y 1 Y 2 + X 1 * Y 3 ] [ ( b c ) w 1 2 v 1 2 w 1 x ( b ) x ( c ) ] 3 c a [ ( c a ) w 2 2 v 2 2 w 2 x ( c ) x ( a ) ] T [ Y 1 Y 2 + X 2 * Y 3 ] [ ( c a ) w 2 2 v 2 2 w 2 x ( c ) x ( a ) ] x T ( c ) ( X 1 X 2 ) x ( c ) + x T ( b ) X 1 x ( b ) x T ( a ) X 2 x ( a )

利用交互凸组合不等式可得

Φ ( b c ) w 1 T Y 1 w 1 2 w 1 T ( Y 2 + X 1 ) ( x ( b ) x ( c ) ) ( c a ) w 2 T Y 1 w 2 2 w 2 T ( Y 2 + X 2 ) ( x ( c ) x ( a ) ) x T ( c ) ( X 1 X 2 ) x ( c ) + x T ( b ) X 1 x ( b ) x T ( a ) X 2 x ( a ) 1 b a [ x ( b ) x ( c ) x ( c ) x ( a ) ] T [ Y 3 X 3 * Y 3 ] [ x ( b ) x ( c ) x ( c ) x ( a ) ] 3 b a [ ( b c ) w 1 2 v 1 2 w 1 x ( b ) x ( c ) ( c a ) w 2 2 v 2 2 w 2 x ( c ) x ( a ) ] T [ Ω 1 X 4 * Ω 2 ] [ ( b c ) w 1 2 v 1 2 w 1 x ( b ) x ( c ) ( c a ) w 2 2 v 2 2 w 2 x ( c ) x ( a ) ] = ξ T Ω ξ

结论成立,证毕。

3. 理论结果及证明

在本节当中,基于上面提到的一重积分不等式引理以及广义的自由矩阵积分不等式引理,我们通过恰当地构造Lyapunov泛函分析得到关于系统(1)的新的时滞依赖条件。

为了简化向量和矩阵表示给出以下术语:

ϖ 1 ( t ) = 1 h ( t ) t h ( t ) t x ( s ) d s , ϖ 2 ( t ) = 1 t h ( t ) t h t h ( t ) x ( s ) d s , ϑ 1 ( t ) = 1 h ( t ) t h ( t ) t θ t x ( s ) d s d θ ,

ϑ 2 ( t ) = 1 t h ( t ) t h t h ( t ) θ t h ( t ) x ( s ) d s d θ , η 2 ( t ) = col { x ( t ) , x ˙ ( t ) , f ( W 0 x ( t ) ) } ,

η 1 ( t , s ) = col { x ( s ) , x ˙ ( s ) , f ( W 0 x ( s ) ) , x ( t ) , x ( t h ) , s t x ( θ ) d θ } , η 3 ( t ) = col { x ( t ) , x ˙ ( t ) } ,

ξ 1 ( t ) = col { x ( t ) , x ( t h ) , ϖ 1 ( t ) , ϑ 1 ( t ) , t h t x ( s ) d s , t h ( t ) t f ( W 0 x ( s ) ) d s , t h t f ( W 0 x ( s ) ) d s } ,

ξ ( t ) = col { x ( t ) , x ( t h ( t ) ) , x ( t h ) , ϖ 1 ( t ) , ϖ 2 ( t ) , ϑ 1 ( t ) , ϑ 2 ( t ) , x ˙ ( t h ( t ) ) , x ˙ ( t h ) , f ( W 0 x ( t ) ) , f ( W 0 x ( t h ( t ) ) ) , t h ( t ) t f ( W 0 x ( s ) ) d s , t h t h ( t ) f ( W 0 x ( s ) ) d s } ,

g i 1 ( s ) = f i ( s ) l i s , g i 2 ( s ) = l i + s f i ( s ) .

定理1 对于给定的实数h和u,若存在对称正定矩阵 P R 7 n × 7 n Q 1 R 5 n × 5 n Q 2 R 5 n × 5 n R R 3 n × 3 n Z = [ Z 1 Z 2 * Z 3 ] R 2 n × 2 n N i R n × n ( i = 1 , 2 ) ,半正定对角矩阵 T 12 R n × n T 13 R n × n T 23 R n × n V k R n × n ( k = 1 , 2 , 3 , 4 ) D i R n × n 任意矩阵 M k R n × n S 1 R 3 n × 3 n S 2 R n × n S 3 R 2 n × 2 n 得下列线性矩阵不等式成立

Γ m = [ Z 1 Z 2 + N i * Z 3 ] > 0 ( m = 1 , 2 ) , Γ 3 = [ Z 3 S 2 * Z 3 ] > 0 Γ 4 = [ Γ 1 S 3 * Γ 2 ] > 0 ,(5)

Γ 5 = [ R S 1 * R ] > 0 , (6)

Γ 6 ( h ( t ) = 0 ) = ( δ ) T [ Ξ h M 3 h M 4 * H 0 * * 3 H ] ( δ ) < 0 , (7)

Γ 7 ( h ( t ) = h ) = ( δ ) T [ Ξ h M 1 h M 2 * H 0 * * 3 H ] ( δ ) < 0 , (8)

Γ 8 = ( δ ) T ( γ 2 h 2 J 2 + Γ 6 ) ( δ ) < 0 , Γ 9 = ( δ ) T ( ( 1 γ ) 2 h 2 J 2 + Γ 7 ) ( δ ) < 0 ,(9)

则时滞神经网络系统(1)是渐近稳定的。

其中

Ξ ( h ( t ) ) = n = 1 8 Π n , Π 1 = sym { G 1 T P G 2 } , Π 3 = h 2 G 9 T R G 9 G 10 T Γ 5 G 10 ,

Π 2 = G 3 T Q 1 G 3 ( 1 u ) G 4 T Q 1 G 4 + G 3 T Q 2 G 3 G 5 T Q 2 G 5 + sym { G 7 T Q 1 G 6 + G 8 T Q 2 G 6 } ,

Π 4 = h G 11 T Z G 11 1 h G 12 T Γ 3 G 12 1 h G 13 T Γ 4 G 13 + sym { e 4 T ( Y 2 + X 1 ) ( e 1 e 2 ) e 5 T ( Y 2 + X 2 ) ( e 2 e 3 ) } h ( t ) e 4 T Y 1 e 4 ( h h ( t ) ) e 5 T Y 1 e 5 e 2 T ( X 1 X 2 ) e 2 + e 1 T X 1 e 1 e 3 T X 2 e 3 ,

Π 5 = h f 0 H f 0 + sym { G 14 T M 1 G 15 + G 14 T M 2 G 16 + G 14 T M 3 G 17 + G 14 T M 4 G 18 } ,

Π 6 = sym { [ ( e 10 L 2 W 0 e 1 ) T V 1 + ( L 1 W 0 e 1 e 10 ) T V 2 ] W 0 f 0 } + sym { [ ( e 12 L 2 W 0 e 3 ) T V 3 + ( L 1 W 0 e 3 e 12 ) T V 4 ] W 0 e 9 } ,

Π 7 = sym { [ e 10 e 11 L 2 W 0 ( e 1 e 2 ) ] T T 12 [ L 1 W 0 ( e 1 e 2 ) ( e 10 e 11 ) ] } + sym { [ e 10 e 12 L 2 W 0 ( e 1 e 3 ) ] T T 13 [ L 1 W 0 ( e 1 e 3 ) ( e 10 e 12 ) ] } + sym { [ e 11 e 12 L 2 W 0 ( e 2 e 3 ) ] T T 23 [ L 1 W 0 ( e 2 e 3 ) ( e 11 e 12 ) ] } + sym { ( e 10 L 2 W 0 e 1 ) T T 1 ( L 1 W 0 e 1 e 10 ) } + sym { ( e 11 L 2 W 0 e 2 ) T T 2 ( L 1 W 0 e 2 e 11 ) } + sym { ( e 12 L 2 W 0 e 3 ) T T 3 ( L 1 W 0 e 3 e 12 ) } ,

Π 8 = sym { ( e 13 L 2 e 4 ) T D 1 ( L 1 e 4 e 13 ) } + sym { ( e 14 L 2 e 5 ) T D 2 ( L 1 e 5 e 14 ) } ,

ϒ = h ( t ) G 14 T ( 3 M 1 H 1 M 1 T + M 2 H 1 M 2 T 3 ) G 14 + ( h h ( t ) ) G 14 T ( 3 M 3 H 1 M 3 T + M 4 H 1 M 4 T 3 ) G 14 ,

e i = [ 0 n , ( i 1 ) n I n 0 n , ( 9 i ) n ] ( i = 1 , 2 , , 9 ) ,

G 1 = col { e 1 , e 3 , h ( t ) e 4 , h ( t ) e 6 , h ( t ) e 4 + ( h h ( t ) ) e 5 , e 13 , e 13 + e 14 } , G 4 = col { e 2 , e 8 , e 11 , e 1 , e 3 , h ( t ) e 4 } ,

G 2 = col { f 0 , e 9 , e 1 ( 1 u ) e 2 , h ( t ) ( e 1 ( 1 u ) e 4 ) , e 1 e 3 , e 10 ( 1 u ) e 11 , e 10 e 12 } ,

G 3 = col { e 1 , f 0 , e 10 , e 1 , e 3 , e 0 } , G 5 = col { e 3 , e 9 , e 12 , e 1 , e 3 , h ( t ) e 4 + ( h h ( t ) ) e 5 } ,

G 6 = col { e 0 , e 0 , e 0 , f 0 , e 9 , e 1 } , G 7 = col { h ( t ) e 4 , e 1 e 2 , e 13 , h ( t ) e 1 , h ( t ) e 3 , h ( t ) e 6 } ,

G 8 = col { h ( t ) e 4 + ( h h ( t ) ) e 5 , e 1 e 3 , e 13 + e 14 , h e 1 , h e 3 , h ( t ) e 6 + ( h h ( t ) ) e 7 + h ( t ) ( h h ( t ) ) e 4 } ,

G 9 = col { e 1 , f 0 , e 10 } , G 10 = col { h ( t ) e 4 , e 1 e 2 , e 13 , ( h h ( t ) ) e 5 , e 2 e 3 , e 14 } , G 11 = col { e 1 , f 0 } ,

G 12 = col { e 1 e 2 , e 2 e 3 } , G 13 = col { h ( t ) e 4 2 e 6 , 2 e 4 e 1 e 2 , ( h h ( t ) ) e 5 2 e 7 , 2 e 5 e 2 e 3 } ,

G 14 = col { e 1 , e 2 , e 3 , e 4 , e 5 , e 6 , e 7 , e 8 , e 9 , e 10 , e 11 , e 12 , e 13 , e 14 } , G 15 = e 1 e 2 , G 16 = e 1 + e 2 2 e 4 ,

G 17 = e 2 e 3 , G 18 = e 2 + e 3 2 e 5 , G 19 = col { e 0 , e 0 , e 4 , e 6 , e 4 e 5 , e 0 , e 0 } ,

G 20 = col { e 0 , e 0 , e 0 , e 1 ( 1 u ) e 4 , e 0 , e 0 , e 0 } , G 21 = col { e 0 , e 0 , e 0 , e 0 , e 0 , e 4 } ,

G 22 = col { e 0 , e 0 , e 0 , e 0 , e 0 , e 4 e 5 } , G 23 = col { e 4 , e 0 , e 0 , e 5 , e 0 , e 0 } , G 24 = col { e 4 , e 0 , e 5 , e 0 } .

证明:构造如下Lyapunov泛函

V ( t ) = d = 1 6 V d ,

V 1 ( t ) = ξ 1 T ( t ) P ξ 1 ( t ) , V 2 ( t ) = t h ( t ) t η 1 T ( t , s ) Q 1 η 1 ( t , s ) d s + t h t η 1 T ( t , s ) Q 2 η 1 ( t , s ) d s ,

V 3 ( t ) = h t h t θ t η 2 T ( s ) R η 2 ( s ) d s d θ , V 4 ( t ) = t h t θ t η 3 T ( s ) Z η 3 ( s ) d s d θ ,

V 5 ( t ) = t h t θ t x ˙ T ( s ) H x ˙ ( s ) d s d θ ,

V 6 ( t ) = 2 i = 1 n 0 W 0 i x ( t ) [ V 1 i g i 1 ( s ) + V 2 i g i 2 ( s ) ] d s + 2 i = 1 n 0 W 0 i x ( t h ) [ V 3 i g i 1 ( s ) + V 4 i g i 2 ( s ) ] d s .

V ( t ) 求导,可得

V ˙ 1 ( t ) = 2 ξ 1 T ( t ) P ξ ˙ 1 ( t ) , V ˙ 3 ( t ) = h 2 η 2 T ( t ) R η 2 ( t ) t h t η 2 T ( s ) R η 2 ( s ) d s ,

V ˙ 2 ( t ) = η 1 T ( t , t ) Q 1 η 1 ( t , t ) ( 1 u ) η 1 T ( t , t h ( t ) ) Q 1 η 1 ( t , t h ( t ) ) + η 1 T ( t , t ) Q 2 η 1 ( t , t ) η 1 T ( t , t h ) Q 2 η 1 ( t , t h ) + ξ T ( t ) sym { G 7 T Q 1 G 6 + G 8 T Q 2 G 6 } ξ ( t ) ,

V ˙ 4 ( t ) = h η 3 T ( t ) Z η 3 ( t ) t h t η 3 T ( s ) Z η 3 ( s ) d s , V ˙ 5 ( t ) = h x ˙ T ( t ) H x ˙ ( t ) t h t x ˙ T ( t ) H x ˙ ( t ) d s ,

V ˙ 6 ( t ) = 2 i = 1 n [ V 1 i g i 1 ( W 0 i x ( t ) ) + V 2 i g i 2 ( W 0 i x ( t ) ) ] W 0 i x ˙ ( t ) + 2 i = 1 n [ V 3 i g i 1 ( W 0 i x ( t h ) ) + V 4 i g i 2 ( W 0 i x ( t h ) ) ] W 0 i x ˙ ( t h ) .

为了考虑更多的积分信息,利用时滞分割方法将一重积分区间 [ 0 , h ] 划分成 [ 0 , h ( t ) ] [ h ( t ) , h ] 两个部分,运用(6)中的条件和Jensen不等式,对 V ˙ 3 ( t ) 中的积分放缩得

t h t η 2 T ( s ) R η 2 ( s ) d s = t h ( t ) t η 2 T ( s ) R η 2 ( s ) d s t h t h ( t ) η 2 T ( s ) R η 2 ( s ) d s ξ T ( t ) G 10 T Γ 5 G 10 ξ ( t ) .

运用(5)中的条件和引理4提出的不等式对 V ˙ 4 ( t ) 中的积分项进行放缩,可以得到

t h t η 3 T ( s ) Z η 3 ( s ) d s ξ T ( t ) { sym { e 4 T ( Y 2 + X 1 ) ( e 1 e 2 ) e 5 T ( Y 2 + X 2 ) ( e 2 e 3 ) } h ( t ) e 4 T Y 1 e 4 ( h h ( t ) ) e 5 T Y 1 e 5 e 2 T ( X 1 X 2 ) e 2 + e 1 T X 1 e 1 e 3 T X 2 e 3 1 h G 12 T Γ 3 G 12 1 h G 13 T Γ 4 G 13 } ξ ( t ) .

运用引理2中的自由矩阵不等式将 V ˙ 5 ( t ) 中的积分项进行放缩,可以得到

t h t x ˙ T ( t ) H x ˙ ( t ) d s = t h ( t ) t x ˙ T ( t ) H x ˙ ( t ) d s t h t h ( t ) x ˙ T ( t ) H x ˙ ( t ) d s ,

t h ( t ) t x ˙ T ( t ) H x ˙ ( t ) d s sym { τ 0 T M 1 δ 1 + τ 0 T M 2 τ 2 } + h ( t ) τ 0 T ( 3 M 1 H 1 M 1 T + M 2 H 1 M 2 T 3 ) τ 0 ,

t h t h ( t ) x ˙ T ( t ) H x ˙ ( t ) d s sym { τ 0 T M 3 τ 3 + τ 0 T M 4 τ 4 } + ( h h ( t ) ) τ 0 T ( 3 M 3 H 1 M 3 T + M 4 H 1 M 4 T 3 ) τ 0 ,

其中

τ 0 = G 14 ξ ( t ) , τ 1 = x ( t ) x ( t h ( t ) ) , τ 2 = x ( t ) + x ( t h ( t ) ) 2 ϖ 1 , τ 3 = x ( t h ( t ) ) x ( t h ) ,

τ 4 = x ( t h ( t ) ) + x ( t h ) 2 ϖ 2 .

另外,由(3)中的Lipchitz条件可知,对任意的正定对角矩阵 D i 0 T i j 0 ( i = 1 , 2 ; j = 2 , 3 ) T k 0 ( k = 1 , 2 , 3 , 4 ) ,有下面的条件成立

2 [ f ( s k ) L s k ] T T k [ f ( s k ) L + s k ] 0 ,

2 [ f ( s i ) f ( s j ) L ( s i s j ) ] T T i j [ f ( s i ) f ( s j ) L + ( s i s j ) ] 0 ,

2 [ t h ( t ) t ( f ( W 0 x ( s ) ) L 1 x ( s ) ) d s ] T D 1 [ t h ( t ) t ( f ( W 0 x ( s ) ) L 2 x ( s ) ) d s ] 0 ,

2 [ t h t h ( t ) ( f ( W 0 x ( s ) ) L 1 x ( s ) ) d s ] T D 2 [ t h t h ( t ) ( f ( W 0 x ( s ) ) L 2 x ( s ) ) d s ] 0 .

Ψ ( h ( t ) ) = Ξ ( h ( t ) ) + ϒ ,由于 h ( t ) 在是二次的,因此可以改写成

Ψ ( h ( t ) ) = J 2 h 2 ( t ) + J 1 h ( t ) + J 0 ,

其中 J 0 J 1 J 2 和都是常数矩阵,且

J 2 = sym { e 19 T P e 20 ( L e 4 ) T D 1 ( L + e 4 ) ( L e 5 ) T D 2 ( L + e 5 ) } ( 1 u ) G 21 T Q 1 G 21 G 22 T Q 2 G 22 G 23 T Γ 5 G 23 1 h G 24 T Γ 4 G 24 .

因此,通过引理1,如果(7)~(9)的条件满足,那么可以得到 Ψ ( h ( t ) ) < 0 成立。

利用Finsler’s引理 [18],让 δ ξ ( t ) = 0 ,如果(5)~(9)的条件满足,我们可以得到 ( δ ) T Ψ ( h ( t ) ) ( δ ) < 0 成立,这就意味着 ξ T ( t ) Ψ ( h ( t ) ) ξ ( t ) < 0 成立。从而定理1的结论成立,证毕。

4. 数值实例

本节结合前面的定理,利用数值实例证明本文所提方法的有效性和优越性。

例:考虑时滞系统(1),其中

A = [ 1.5 0 0 0.7 ] , W 0 = [ 0.0503 10.0454 0.0987 0.2075 ] , W 1 = [ 0.2381 0.9320 0.0388 0.5062 ] ,

W 2 = diag { 1 , 1 } , L 1 = diag { 0.3 , 0.8 } , L 2 = 0 .

运用定理1,结合Matlab计算,当u取不同值时,通过调节参数 γ 可以得到不同的时滞上界。如表1所示,当 u = 0.55 γ = 0.5 时,我们得到的时滞上界 h = 12.4199 ,它比文献 [1], [13], [2], [8], [6], [12], [3] 以及 [19] 分别大265.02%,264.11%,215.47%,260.68%,100.01%,82.22%,41.01%,40.86%。另外,令

ϖ ( t ) = 12.4199 2 + 12.4199 2 ( sin ( 1.1 t 12.4199 ) ) , f ( t ) = [ 0.3 tanh ( x 1 ( t ) ) , 0.8 tanh ( x 2 ( t ) ) ] T

这时,取初始条件 x ( θ ) = [ 0.7 , 0.3 ] T ,从图1的状态运行状况,我们可以清楚地看到,时滞神经网络系统是渐进稳定的。

Table 1. The upper bound of the maximum allowable time delay of h when u takes different values

表1. 当u取不同值时h最大允许的时滞上界

5. 结论

本文探讨了时滞神经网络的稳定性问题,首先通过两个零等式和交互凸组合不等式引入自由矩阵改进具有增广的一重积分不等式,然后在充分考虑各个状态向量的基础上构造带有增广的Lyapunov泛函,并利用提出的一重积分不等式以及广义的自由矩阵积分不等式处理泛函的导数。最后,通过一个数值例子来说明所得结论的有效性和优越性。本文的方法还可以扩展到具有马尔科夫跳跃的时滞神经网络系统的稳定性以及控制问题的研究中。

Figure 1. State responses of the neural network

图1. 神经网络的状态响应

基金项目

国家自然科学资金(12061088);云南省教育厅科学研究资金项目(2020Y0236)。

文章引用

柳莹莹,熊良林,邹 梅,蔡 丽,王 珍. 改进的一重积分不等式在时滞神经网络中的应用
The Application of Improved Single Integral Inequality in Neural Networks with Time-Delayed[J]. 理论数学, 2021, 11(01): 79-89. https://doi.org/10.12677/PM.2021.111012

参考文献

  1. 1. Zhang, C.K., He, Y., Jiang, L., Wu, Q.H. and Wu, M. (2013) Delay-Dependent Stability Criteria for Generalized Neural Networks with Two Delay Components. IEEE Transactions on Neural Networks and Learning Systems, 25, 1263-1276. https://doi.org/10.1109/TNNLS.2013.2284968

  2. 2. Kwon, O.M., Park, M.J., Lee, S.M., Park, J.H. and Cha, E.-J. (2012) Stability for Neural Networks with Time-Varying Delays via Some New Approaches. IEEE Transactions on Neural Networks and Learning Systems, 24, 181-193. https://doi.org/10.1109/TNNLS.2012.2224883

  3. 3. Yang, B., Wang, J. and Liu, X.D. (2017) Improved De-lay-Dependent Stability Criteria for Generalized Neural Networks with Time-Varying Delays. Information Sciences, 420, 299-312. https://doi.org/10.1016/j.ins.2017.08.072

  4. 4. Zhang, X.M. and Han, Q.L. (2018) State Estimation for Static Neural Networks With Time-Varying Delays Based on an Improved Reciprocally Convex Inequality. IEEE Transactions on Neural Networks & Learning Systems, 29, 1376-1381. https://doi.org/10.1109/TNNLS.2017.2661862

  5. 5. 廖晓昕. 稳定性理论、方法和应用[M]. 武汉: 华中科技大学出版社, 2010.

  6. 6. Zhang, C.K., He, Y., Jiang, L. and Wu, M. (2015) Stability Analysis for Delayed Neural Net-works Considering both Conservativeness and Complexity. IEEE Transactions on Neural Networks and Learning Sys-tems, 27, 1486-1501. https://doi.org/10.1109/TNNLS.2015.2449898

  7. 7. Ge, C., Hua, C.C. and Guan, X.P. (2017) New De-lay-Dependent Stability Criteria for Neural Networks With Time-Varying Delay Using Delay-Decomposition Approach. IEEE Transactions on Neural Networks & Learning Systems, 25, 1378-1383. https://doi.org/10.1109/TNNLS.2013.2285564

  8. 8. Zhang, H.G., Yang, F.S., Liu, X.D. and Zhang, Q.L. (2013) Stability Analysis for Neural Networks with Time-Varying Delay Based on Quadratic Convex Combination. IEEE Transactions on Neural Networks and Learning Systems, 24, 513-521. https://doi.org/10.1109/TNNLS.2012.2236571

  9. 9. Gu, K., Kharitonov, V. L. and Chen, J. (2013) Stability of Time-Delay Systems. Birkhäuser, Boston.

  10. 10. Seuret, A. and Gouaisbaut, F. (2013) Wirtinger-Based Integral Ine-quality: Application to Time-Delay Systems. Automatica, 49, 2860-2866. https://doi.org/10.1016/j.automatica.2013.05.030

  11. 11. Zeng, H.B., He, Y., Wu, M. and She, J.H. (2015) Free-Matrix-Based Integral Inequality for Stability Analysis of Systems with Time-Varying Delay. IEEE Transactions on Automatic Control, 60, 2768-2772. https://doi.org/10.1109/TAC.2015.2404271

  12. 12. Zhang, C.K., He, Y., Jiang, L., Lin, W.-J. and Wu, M. (2017) Delay-Dependent Stability Analysis of Neural Networks with Time-Varying Delay: A Generalized Free-Weighting-Matrix Approach. Applied Mathematics and Computation, 294, 102-120. https://doi.org/10.1016/j.amc.2016.08.043

  13. 13. Li, T., Wang, T., Song, A.G. and Fei, S.M. (2013) Combined Convex Technique on Delay-Dependent Stability for Delayed Neural Networks. IEEE Transactions on Neural Networks and Learning Systems, 24, 1459-1466. https://doi.org/10.1109/TNNLS.2013.2256796

  14. 14. Kwon, O.M., Park, M.J., Park, J.H., Lee, S.M. and Cha, E.J. (2014) Improved Results on Stability of Linear Systems with Time-Varying Delays via Wirtinger-Based Integral Ine-quality. Journal of the Franklin Institute, 351, 5386-5398. https://doi.org/10.1016/j.jfranklin.2014.09.021

  15. 15. Zhang, R., Zeng, D., Park, J.H., Zhong, S.M., Liu, Y.J. and Zhou, X. (2019) New Approaches to Stability Analysis for Time-Varying Delay Systems. Journal of the Franklin In-stitute, 356, 4174-4189. https://doi.org/10.1016/j.jfranklin.2019.02.029

  16. 16. Zhang, C.K., Long, F., He, Y., Yao, W., Jiang, L. and Wu, M. (2020) A Relaxed Quadratic Function Negative-Determination Lemma and Its Application to Time-Delay Systems. Automatica, 113, Article ID: 108764. https://doi.org/10.1016/j.automatica.2019.108764

  17. 17. Park, P.G., Lee, W.I. and Lee, S.Y. (2015) Auxiliary Function-Based Integral Inequalities for Quadratic Functions and Their Applications to Time-Delay Systems. Journal of the Franklin Institute, 352, 1378-1396. https://doi.org/10.1016/j.jfranklin.2015.01.004

  18. 18. Skelton, R.E., Iwasaki, T. and Grigoriadi, K. (1997) A Unifisfied Algebraic Approach to Control Design. CRC Press, Boca Raton.

  19. 19. Zhang, R.M., Zeng, D.Q., Liu, X.Z., Zhong, S.M. and Cheng, J. (2019) New Results on Stability Analysis for Delayed Markovian Generalized Neural Net-works with Partly Unknown Transition Rates. IEEE Transactions on Neural Networks and Learning Systems, 30, 3384-3395. https://doi.org/10.1109/TNNLS.2019.2891552

  20. NOTES

    *通讯作者。

期刊菜单