局部耦合的双层忆阻神经元网络的嵌合态和同步

第２１卷第９期２０２３年９月
动力学与控制学报
J O U R N A L O FD Y N AM I C SA N DC O N T R O L
V o l ．２１N o ．９
S e p
．２０２３文章编号:１６７２Ｇ６５５３Ｇ２０２３Ｇ２１(９)Ｇ０６７Ｇ００７
D O I :１０．６０５２/１６７２Ｇ６５５３Ｇ２０２２Ｇ０３３
㊀２０２２Ｇ０５Ｇ１５收到第１稿,２０２２Ｇ０６Ｇ２３收到修改稿．∗国家自然科学基金资助项目(１２２７２０９２,１１９７２１１５),N a t i o n a lN a t u r a l S c i e n c eF o u n d a t i o no fC h i n a (１２２７２０９２,１１９７２１１５)．
†通信作者E Ｇm a i l :y a d i a h a n ＠d h u ．e d u ．c n
局部耦合的双层忆阻神经元网络的嵌合态和同步∗
许一凡１㊀袁园园１㊀韩芳２†
㊀王直杰２
(东华大学信息科学与技术学院,上海㊀２０１６２０
)摘要㊀嵌合态被发现存在于神经系统并且可能在神经元节律㊁大脑的睡眠和记忆等诸多神经过程中发挥重要作用．本文考虑神经元交互中的电磁感应现象,建立了以H i n d m a r s h ＧR o s e 神经元为节点的局部耦合双层忆阻神经元网络,研究其嵌合态时空动力学模式及产生机理．结果表明,改变层内㊁层间突触耦合强度会使网络产生移动和不完美移动嵌合态等多种类型的嵌合模式,其中不完美移动嵌合态中不相干的区域会扩展到网络的相干域．特别地,在特定耦合强度下,存在一种新的嵌合态活动模式,即一部分神经元处于嵌合态,另一部分神经元处于移动嵌合态．考虑神经元突触的忆阻特性,发现忆阻参数的增加能够使处于嵌合态的神经元网络转变为同步态,且耦合强度越大,达到同步态所需要的忆阻参数值越小．进一步探究双层网络的同步性,发现层间耦合强度和忆阻参数的增大有助于网络达到更好的同步．研究结果表明神经元之间的相互作用可以激发双层神经元网络产生多种嵌合态模式,电磁感应可以促进网络由嵌合态向同步态转迁,这些结果有助于理解人脑中复杂的神经放电过程和信息处理机制,并为可能的类脑装置应用提供参考．关键词㊀嵌合态,㊀双层神经元网络,㊀忆阻突触,㊀局部耦合中图分类号:Q ４２
文献标志码:A
C h e m i r a S t a t e s a n dS t n c h r o n i z a t i o n i naT w o ＧL a y e rL o c a l l y C o u p
l e d M e m r i s t i v eH i n d m a r s h ＧR o s eN e u r o n a lN e t w o r k
∗
X uY i f a n １㊀Y u a nY u a n y u a n １㊀H a nF a n g ２†㊀W a n g Z h i j
i e ２(D e p a r t m e n t o f I n f o r m a t i o nS c i e n c e a n dT e c h n o l o g y ,D o n g h u aU n i v e r s i t y ,S h a n g
h a i ㊀２０１６２０,C h i n a )A b s t r a c t ㊀C h i m e r a s t a t e s h a v eb e e n f o u n d i nn e r v o u s s y s t e m s a n dm a yp l a y a n i m p o r t a n t r o l e i n m a n y
n e u r a l p r o c e s s e s i n c l u d i n g n e u r o n a l r h y t h m s ,s l e e p a n dm e m o r y i n t h e b r a i n ．C o n s i d e r i n g t h e p h e n o m e Ｇn o no f e l e c t r o m a g n e t i c i n d u c t i o n i n t h e n e u r o n i n t e r a c t i o n s ,a t w o Ｇl a y
e rm e m r i s t i v e n e u r o n a l n e t w o r ko
f l o c a l l y c o u p l e dH i n d m a r s h ＧR o s e n e u r o n sw a s e s t a b l i s h e d t o s t u d y t h e s p a t i o t e m p o r a l d y n a m i c a l p a t t e r n s o f c h i m e r a s t a t e s a n d t h em e c h a n i s m s o f t h e i r
g e n e r a t i o n ．I t i s f o u n d t
h a t c h a n g
i n g t
h e i n t r a Ｇa n d i n t e r Ｇl a y e r s y n a p t i c c o u p l i n g s t r e n g t h s c a n i n d u c e s e v e r a l t y p e s o f c h i m e r a s t a t e s i nn e t w o r k s ,i n c l u d i n g t r a v e Ｇl i n g a n d i m p e r f e c t t r a v e l i n g c h i m e r a s t a t e sw h e r e t h e i n c o h e r e n t t r a v e l i n g d o m a i n s p r e a d s i n t o t h e c o h e r Ｇe n t d o m a i no f t h en e t w o r k ．I n p a r t i c u l a r ,a p h e n o m e n o n t h a t s e v e r a l n e u r o n s a r e i nc h i m e r as t a t e s a n d o t h e r s a r e i n t r a v e l i n g c h i m e r a s t a t e s i sw i t n e s s e d a t s p e c i f i c c o u p l i n g s t r e n g
t h s ．M o r e o v e r ,i t i s a l s o o b Ｇs e r v e d t h a t i n c r e a s i n g t h em e m r i s t i v e p a r a m e t e r c a n i n d u c e t h e n e u r o n a l n e t w o r k t o c h a n g
e f r o mc h i m e r a s t a t e s t o s y n c h r o n o u s s t a t e s ,a n d t h e s t r o n g e r t h e c o u p l i n g s t r e n g t h s ,t h e s m a l l e r t h e v a l u e o f t h em e m Ｇr i s t i v e p a r a m e t e r r e q u i r e d t o r e a c hs y n c h r o n o u s s t a t e s ．C o n s i d e r i n g t h e i n t e r Ｇl a y e r s y
n c h r o n i z a t i o n ,t h e i n c r e a s eo f t h e c o u p l i n g s t r e n g t h s a n d t h em e m r i s t i v e p a r a m e t e r h e l p s t h e n e t w o r k t o a c h i e v e b e t t e r s y
n Ｇ
动㊀力㊀学㊀与㊀控㊀制㊀学㊀报２０２３年第２１卷c h r o n i s a t i o n．T h e r e s u l t s d e m o n s t r a t e t h a t t h e i n t e r a c t i o n sb e t w e e nn e u r o n so n t w oＧl a y e rn e u r o n a l n e tＧw o r k s c a n s t i m u l a t e a v a r i e t y o f p a t t e r n s o f c h i m e r a s t a t e s,a n d t h a t t h e e f f e c t o f e l e c t r o m a g n e t i c i n d u cＧt i o n c a n f a c i l i t a t e t h e t r a n s i t i o n f r o mc h i m e r a s t a t e s t oc o h e r e n t s t a t e s．T h e s e r e s u l t s a r eu s e f u l f o ru nＧd e r s t a n d i n g t h ec o m p l e xn e u r a l f i r i n gp r o c e s s e sa n d i n f o r m a t i o n p r o c e s s i n g m e c h a n i s m s i nt h eh u m a n b r a i na n dm a yp r o v i d e a r e f e r e n c e f o r t h e a p p l i c a t i o no f t h eb r a i nＧl i k e d e v i c e s．
K e y w o r d s㊀c h i m e r a s t a t e s,㊀t w oＧl a y e r n e u r o n a l n e t w o r k,㊀m e m r i s t i v e s y n a p s e s,㊀l o c a l c o u p l i n g
引言
人脑中有大量神经元,这些神经元组成复杂的网络结构参与神经信息处理以及大脑的其他活动,通过建模研究神经元网络的动态活动机制具有十分重要的意义．２００２年,K u r a m o t o在由相同振荡
器通过非局部耦合形成的系统中发现了一种同步与非同步态共存的复杂动态模式[１],后来A b r a m s 等人[２]引用古希腊神话中的神兽将这种奇特的现象命名为 C h i m e r a ．二十年来,人们对嵌合态的研究逐渐深入,在生物㊁物理和化学等多学科领域进行了广泛探索,证实它不仅可以存在于相位振荡器,还存在于混沌系统[３,４]和H o p f范式[５,６]等诸多系统中．在神经系统中,自然界中的多种神经过程都证明了嵌合态的存在,如一些鸟类和水生哺乳动物的单半球慢波睡眠现象和人类的首晚现象[７],这种同步与不同步共存的现象还与帕金森症㊁癫痫发作等多种神经疾病密切相关[８,９]．因此,神经元网络中的嵌合态获得了广泛的关注,在H i n d m a r s hＧR o s e[９,１０]㊁F i t z H u g hＧN a g u m o[１１]和M o r r i sＧL eＧc a r[１２]等多种神经元模型构成的多种网络拓扑结构中都观察到了嵌合态,H u o等人[１３]研究了随机网络㊁全连接网络和无标度网络中的嵌合态特征,证实了网络结构会影响嵌合态出现的可能性．H i z a n iＧd i s等人[１４]建立了基于秀丽隐杆线虫土壤蠕虫连接组的神经元网络模型,证实了模块化网络也会出现类嵌合态现象．为了区分嵌合态的多种特征,人们根据时空动力学将其划分为移动嵌合态㊁暂态嵌合态㊁不完美嵌合态[１５,１６]等类型,同时对可能影响嵌合态存在的耦合强度㊁网络拓扑和时延[１７]等多种因素进行了深入研究．
近些年来,对复杂网络集体行为的研究在生物科学㊁计算机技术和统计物理等[１８,１９]多个领域引起了广泛重视．局部耦合网络是常见的规则网络连接方式,它在感知分组和模式分割中有很大潜力[２０],且局部耦合的非线性特征能在嵌合态的产
生中发挥重要作用[２１]．大脑皮层网络十分复杂,单层网络不足以说明系统功能,需要用多层网络结构模拟更加真实的系统．多层结构势必会对神经元网络的动力学行为产生影响．特别地,在多层网络中也发现了嵌合态现象,R u z z e n e等人[２２]通过构建由F i t z H u g hＧN a g u m o神经元组成的多层网络,研究了利用复用结构远程控制嵌合态的可能性,为控制嵌合态提供了新思路．M a j h i等人[２３]研究了多层神经元网络中突触和信息传输延迟对嵌合态和同步簇模式的影响．
神经元之间的信息传递主要通过电突触和化学突触两种功能性连接进行,已经在电和化学突触耦合的多种神经元网络中观察到嵌合态现象．生物神经网络中信号传输和离子交换可能会产生复杂的电磁场,电磁场与膜电位之间存在忆阻联系,并且神经元突触连接本质上是有记忆特性的．因此选择忆阻突触构建神经元网络能够有效模拟生物突触动态行为,也可以帮助解释神经系统的信息处理和调控机制．王松等[２４]研究含有忆阻突触的两耦合神经元系统,证实了信号传输中的时滞可以影响系统稳定性和放电模式．马杰等[２５]研究了含有时滞的忆阻突触耦合的四维H o p f i e l d神经网络,发现时滞和忆阻耦合强度会使系统产生多种时空动力学行为．L i u等[２６]构建由不同神经元组成的层间忆阻突触连接的多层神经元网络,探究其动态鲁棒性及放电模式的转换,发现弱电磁感应会破坏动态鲁棒性．但是在生物神经元网络中突触的记忆性广泛存在,探究神经元之间的忆阻特性与复杂网络的时空动力学特征之间的关系具有十分重要的意义．本文着力于探究忆阻突触耦合的神经元网络中的嵌合态时空动力学行为,构建了一个基于H i n dＧm a r s hＧR o s e神经元的局部耦合双层神经元网络,层
８６
第９期
许一凡等:局部耦合的双层忆阻神经元网络的嵌合态和同步
内和层间节点都用忆阻突触连接．首先,研究在不同耦合强度下神经元网络出现的多种嵌合态动力学行为．然后,探讨忆阻参数对神经元网络动力学行为的影响和不同参数与网络层间同步性的关系．
１㊀模型描述与分析方法
１．１㊀神经元网络模型
本文中,神经元网络为双层结构,单个振子与其在层内两个方向上最近邻和另一层网络对应位置上的神经元连接(称为局部耦合),用H i n d Ｇ
m a r s h ＧR o s e 神经元作为网络的节点,用忆阻突触模拟生物突触连接层内和层间神经元．网络结构如图１所示(红色和紫色的圆分别代表网络第一层和第二层中的神经元,蓝色的长方形代表忆阻突触),网络动力学用如下方程描述
:
图１㊀局部耦合的双层忆阻神经元网络结构图
F i g ．１㊀S c h e m a t i c d i a g r a mo f t h e t w o Ｇl a y
e rm e m r i s t i v e n e u r a l n e t w o r ko
f l o c a l l y c o u p
l e d x
i ,１
＝a x ２i ,１－x ３
i ,１－y i ,１－z i ,１＋
㊀ε１[M (φi －１,１)(x i －１,１－x i ,１)＋M (φi ,１)(x i ＋１,１－㊀x i ,１)]－λ３M (φ
i ,１２)(x i ,１－x i ,２)y i ,１
＝(a ＋α)x ２
i ,１－y i ,１
z
i ,１＝w (b x i ,１－z i ,１＋
c )φ
i ,１
＝x i ,１－x i ＋１,１－δ１φ
i ,１x
i ,２＝a x ２i ,２－x ３
i ,２－y i ,２－z i ,２＋
㊀ε２[M (φi －１,２)(x i －１,２－x i ,２)＋M (φi ,２)(x i ＋１,２－㊀x i ,２)]＋λ３M (φ
i ,１２)(x i ,１－x i ,２)y
i ,２
＝(
a ＋α)x ２
i ,２
－y i ,２
z i ,２＝w (b x i ,２－z i ,２＋c )φ
i ,２＝x i ,２－x i ＋１,２－δ２φ
i ,２φ
i ,１２＝x i ,１－x i ,２－δ３φ
i ,１２(１)式中,x ㊁y 和z 分别表示单个神经元的膜电位㊁
恢复变量和与钠㊁钾离子有关的慢适应电流．N 表示每层网络中的神经元数量,x j ０＝x j N ,x j N ＋１＝x j １,j ＝
１,２表示两层结构,式中的δ表示忆阻器的遗忘参数,参数ε代表每层神经元之间忆阻突触的耦合强度,λ３代表层间忆阻突触的耦合强度,变量φi ,j 和
M (φi ,j )分别表示第j 行第i 个忆阻突触的磁通量和忆感,M (φi ,j )是q (φ
i ,j )电荷相对于磁通量φi ,j 的导数,使用三阶磁通控制忆阻突触,其电荷方程表示为:
q (φi ,j )＝σφi ,j ＋θφ
３
i ,j (２)对磁通量求导得到忆阻的忆感表示为:
M (φi ,j )＝σ＋３θφ
２
i ,j (３
)式中,σ和θ是忆阻器的参数,
它们的取值与环境和忆阻器自身条件有关,将神经元的参数设置为a ＝
１．４５,α＝１．６,u ＝０．００１,b ＝９,c ＝５,δj ＝０．５,δ３＝０．５,θ＝０．０２．
对式(１
)中每个神经元的变量设置初始值,参考文献[２７],设置前N /２个神经元为有规律的初始值,后N /２个神经元的初始条件加入小的随机波动,表示为:当i ＝１,２, ,N /２时,x i ０,j ＝０
．０１(i －N /２),y i ０,j ＝０．０２(i －N /２),z i ０,j ＝０．０３(i －N /２),对其他神经元设置为x i ０,j ＝０．１(N /２－i ),y i ０,j ＝
０．１２(N /２－i ),z i ０,j ＝０．２１(N /２－i )．１．２㊀分析方法
为了更好地区分嵌合态与其他状态,引入不相
干强度(S I ,t h e s t r e n g
t h o f i n c o h e r e n c e )[２８
]这一统计变量,首先对相邻神经元膜电位x i ,j 求取差值,
得到w i ,j ＝x i ,j －x i ＋１,j ,w i ,j ＝(w １,i ,j ,w ２,i ,j , ,w d ,i ,j )
ɪR d
,为了更加清楚地量化同步状态,将每层神经元分为M 组大小为n ＝N /M 的部分,因此局部标准变量(l o c a l s t a n d a r dd e v i a t i o n )ηl (m )可以描述为:
㊀ηi ,j (m )＝‹１
n ðm n
k ＝n (m －１)＋
１[w l ,k ,j －‹w i ,j ›
]›t (４
)m ＝１,２, ,M ,l ＝１,２, ,d ,i ＝１,２, ,N ,d 为
神经元模型的维度,‹w i ,j ›＝１
n ðN
i ＝１w l ,i ,j (t ),‹ ›t 代表一段时间内的平均值,η
i ,j (m )可以计算每n 个９
６
动㊀力㊀学㊀与㊀控㊀制㊀学㊀报２０２３年第２１卷
连续神经元的局部标准偏差,不相干强度描述为:
㊀S I j ＝
１－ðM
m ＝１
s m ,j M
,s m ,j ＝Θ[δ０－η
i ,j (m )](５)其中Θ( )是阶跃函数,δ０是提前设定的阈值并通常很小,S I ＝０,S I ＝１和０＜S I ＜１分别表示同步,不同步和嵌合态．另外,使用局部序参数(l o c a lo r d e r p a r a m e Ｇ
t e r
)的统计量[２９]
更加直观地辨别嵌合态㊁不同步态
和同步态,局部序参数代表神经元网络的局部有序性,表示不相干的程度．其公式描述如下:
L i ,j ＝１２p ð|i －k |ɤp
e η
Φk (６
)式中η＝
－１,
i ＝１,２, ,N 表示神经元个数,k ＝１,２, ,N ,j ＝１,２．参数p 指每个神经元两边最近邻的数量,Φk ＝a r c t a n (y k /
x k )是第k 个神经元的几何相位[３０
]．L i ,j ʈ０表明第
j 层第i 个神经元属于不相干态的节点,L i ,j ʈ１表明第j 层第i 个神经元属于相干部分,L i ,j ＝
１代表完全同步．为了进一步探讨两层网络动态行为之间的关
系,引入层间同步误差[３１
],其公式描述为:
E ＝l i m
ΔT ➝ɕ１
ΔT
ʏ
ΔT
１
(M －１)N ðM －１j ＝１ðN
i ＝１
x i ,j －x i ,j ＋１d t (７
)M 和N 分别代表网络中的层数和每一层神经元的数量．E ＝０代表网络达到最佳同步状态．
２㊀结果与分析
２．１㊀耦合强度诱发多种嵌合态
首先,将层内和层间突触耦合强度参数设置为ε１＝５．５㊁λ３＝０．１㊁ε２＝
０,第一层和第二层层内耦合强度较大的差值使网络的两层呈现出不同的状态,此时忆阻参数固定在σ＝０．１２[３２]
．
图１给出了双层网络时空响应图㊁神经元时间序列㊁神经元群在某
一时刻的膜电位分布和局部序参量时空图．从图２
(b
)可以看出第一层网络中一些神经元同步,而其余的神经元不同步,这是嵌合态的特征．图２(a )显示了第一层所有神经元的膜电位在一段时间内随时间演化过程,神经元的同步状态随时间发生变化,呈现出移动嵌合态,图２(c )中的局部序参量的时空响应也与之相符．第二层网络的神经元在移动
过程中一些相干态的神经元被转化为不相干态,这种非相干域扩展到相干域的特殊移动嵌合态是不完美移动嵌合态．图２(d )清楚地显示了在时间和空间移动的过程中宽度变化的时空图,从时刻所有神经元膜电位分布的快照可以观察到相干与不相干的神经元共存,有三组不相干簇,放电方式为高原簇放电,如图２(e )所示．计算出每个神经元的局部序参量,如图２(f )所示,黄色的区域表示相干节点,两个连续的相干域之间蓝色的不相干域在时空中传播时,其宽度会发生变化,出现不相干域向相干域扩展的情况,这是不完美移动嵌合态的明确标志．
图２㊀在ε１＝５．５㊁λ３＝０．１㊁ε２＝０下双层忆阻神经网络的
时空响应图㊁膜电位分布情况及时间序列㊁局部序参量
F i g ．２㊀S p a t i o t e m p o r a l p
l o t s ,t i m e s e r i e s a n d t h e d i s t r i b u t i o n ,t h e l o c a l o r d e r p a r a m e t e r i n t w o Ｇl a y
e rm e m r i s t i v e n e u r o n a l n e t w o r kw i t h ε１＝５．５,λ３＝０
．１a n d ε２＝０特别地,当层间耦合强度较大但层内耦合强度都很小,即ε１＝０．０２㊁λ３＝５㊁ε２＝０．０２时,两层网络都出现了嵌合态与移动嵌合态共存的现象．以第一层网络为例,神经元膜电位时空分布图呈现出规则的条纹状,而右半部分为倾斜的条纹,表示神经元的状态随着时间推移发生变化,如图３(a
)所示．图３(b )表示t ＝３９００m s 时所有神经元的膜电位分布情况,以i ＝４０为界分为两部分,第Ⅰ部分神经元处于无序嵌合态,第Ⅱ部分的神经元被划分为几个相干簇,一部分神经元放电,另一部分属于静息,少０
７
第９期许一凡等:局部耦合的双层忆阻神经元网络的嵌合态和同步
数神经元处于不相干态,且每个神经元的同步状态会随着时间的推移发生改变,处于移动嵌合态．局部序参量的时空响应图更加清楚地表现这一现象,如图３(c)所示,第Ⅰ部分中相干域与不相干域交替出现且保持静态,属于嵌合态,第Ⅱ部分中神经元相干与不相干的状态随时间发生变化,属于移动嵌合态．由此得出,耦合强度对多种类型嵌合态的产生有重要作用．
图３㊀在ε１＝０．０２㊁λ３＝５㊁ε２＝０．０２下忆阻神经网络第一层的时空响应图㊁膜电位分布情况及时间序列㊁局部序参量
F i g．３㊀S p a t i o t e m p o r a l p l o t s,t i m e s e r i e s a n d t h e d i s t r i b u t i o n, t h e l o c a l o r d e r p a r a m e t e r i n t h e f i r s t l a y e r o f t h em e m r i s t i v e
n e u r o n a l n e t w o r kw i t hε１＝０．０２,λ３＝５a n dε２＝０．０２２．２㊀忆阻参数促进嵌合态向同步转迁
调整层内和层间耦合强度,神经元网络中没有出现同步态现象．当ε１＝１．５㊁λ３＝０．８㊁ε２＝０时,神经元网络表现为移动嵌合态,此时忆阻参数σ＝０．１２．在耦合强度不变的条件下增大忆阻参数,到σ＝４．２时所有神经元的膜电位显示出周期放电的模式,且神经元间膜电位彼此关联,如图４(a)㊁４(b)所示,神经元处于完全同步态．图４(c)计算出一段时间内的局部序参量,所有神经元在时空内都保持同步状态．由此看出,忆阻参数σ的增大能打破神经元网络的嵌合态,最终使所有神经元均处于同步态．采用不相干强度S I来定量分析双层网络的动力学行为,图５(a)㊁５(b)分别是两层网络中设置λ３＝０．８㊁ε２＝０而ε１取１．５㊁２㊁３不同值时,所有神经元的状态随忆阻参数变化的情况．可以看出不同层间耦合强度值ε１＝１．５㊁２㊁３,网络先处于嵌合态(０＜S I＜１),忆阻参数增大到σ＝４．５㊁３．５㊁３之后网络呈现同步态(S I＝０),且层内耦合越大,想要达到同步态所需要的忆阻参数值越小．同样地,图５(c)和５(d)分别表示两层网络在不同强度层间作用时,神经元状态随忆阻参数变化的情况,此时ε１＝３㊁ε２＝０．随着层间作用的增强,神经元网络可以在较小的忆阻参数时就从嵌合态转变为同步态．由此可以得出,层内㊁层间相互作用和忆阻参数的增加都可促进网络向同步态的转迁．
图４㊀在ε１＝１．５㊁λ３＝０．８㊁ε２＝０㊁σ＝４．２下忆阻神经网络第一层的时空响应图㊁膜电位分布情况及时间序列㊁局部序参量F i g．４㊀S p a t i o t e m p o r a l p l o t s,t i m e s e r i e s a n d t h e d i s t r i b u t i o n,t h e l o c a l o r d e r p a r a m e t e r i n t h e f i r s t l a y e r o f t h em e m r i s t i v e n e u r o n a l n e t w o r kw i t hε１＝１．５,λ３＝０．８,ε２＝０,σ＝４．
２
图５㊀双层神经元网络中不同耦合强度下
不相干强度随忆阻参数变化曲线
F i g．５㊀V a r i a t i o no f t h e s t r e n g t ho f i n c o h e r e n c ew i t h t h e
m e m r i s t i v e p a r a m e t e r f o r d i f f e r e n t v a l u e s o f c o u p l i n g
s t r e n g t h s i n t h e t w oＧl a y e r n e u r o n a l n e t w o r k ２．３㊀不同参数对层间同步性的影响
为了探究层内㊁层间相互作用和忆阻突触对层间同步动态行为的影响,设置㊁,作出神经元网络在不同忆阻参数下层间同步误差随层间耦合强度变化情况的曲线图．从图６(a)可以看出,随着忆阻参
１７
动㊀力㊀学㊀与㊀控㊀制㊀学㊀报２０２３年第２１卷
数的增大层间同步误差有微小的下降,没有层间作用时层间误差很大,一旦层间作用增大,误差就迅速减小,并且层间耦合强度持续增大层间误差一直保持在较小的值．图６(b )是ε１＝１㊁σ＝０．１２且层间耦合强度取不同值时,层间误差随层内耦合变化的情况,ε２ɤ１时层间误差随层内耦合的增大而减小,在ε２＝１时两层网络层内耦合强度相同,每个神经元受到来自其他神经元的输入电流相同,误差为零,网络达到完全同步状态,而ε２＞１时误差随层内耦合的增大而增大,但是在层间作用较大时这种变化并不明显．综上,层间作用和忆阻参数的增大均会促进双层网络的层间同步
．
图６㊀层间同步误差变化曲线
F i g ．６㊀V a r i a t i o n c u r v e o f i n t e r Ｇl a y e r s y
n c h r o n i z a t i o ne r r o r ３㊀结论
本文构建了一个局部耦合的忆阻突触连接的H i n d m a r s h ＧR o s e 双层神经元网络,研究不同耦合强度和忆阻参数下神经元网络的嵌合态动力学．首先调节层内和层间耦合强度,出现两层网络分别表现不同状态的情况,一层网络中的神经元状态为随时间发生变化的移动嵌合态,另一层出现移动过程中非相干域扩展到相干域的现象,为不完美移动嵌合态．继续改变耦合强度值,发现两层都出现了一部分神经元处于嵌合态,另一部分神经元处于移动嵌合态的现象,说明耦合强度可以激发神经元网络出现多种嵌合态行为．然后,增大忆阻突触参数,观察到了调整耦合强度时没有出现的同步态,说明忆阻参数的增加能促使网络由嵌合态转迁为同步态．进一步发现层内和层间耦合强度增大,忆阻参数在更小值就能使网络转变为同步态．最后探究了不同参数对层间同步性的影响,结果表明增大层间耦合强度和忆阻参数会促进双层网络层间同步,且在两层网络层内耦合强度相同时达到最佳同步状态．上述研究对探索神经系统中的信息处理和放电节律
有关机制有重要帮助,可能为类脑装置的实现提供了新思路．
参考文献
[１
]K U R AMO T O Y ,B A T T O G T O K H D．C o e x i s t e n c e
o f c o h e r e n c ea n di n c o h e r e n c ei nn o n l o c a l l y c o u p l e d p
h a s e o s c i l l a t o r s [J ]．N o n l i n e a rP h e n o m e n a i nC o m Ｇp l e xS y
s t e m s ,２００２,５(４):３８０－３８５．[２
]A B R AM SD M ,S T R O G A T ZS H．C h i m e r as t a t e s
f o r c o u p l e d o s c i l l a t o r s [J ]．P h y s i c a l R e v i e wL e t t e r s ,２００４,９３(１７):１７４１０２．
[３]
OM E L C H E N K O I ,MA I S T R E N K O Y ,HÖV E L
P ,e t a l ．L o s so f c o h e r e n c e i nd y
n a m i c a l n e t w o r k s :s p a t i a l c h a o sa n dc h i m e r as t a t e s [J ]．P h y
s i c a lR e Ｇv i e wL e t t e r s ,２０１１,１０６(２３):２３４１０２．
[４]
P A V L O V A O N ,P A V L O V A N．P r e d i c t i o n o f
c o m p l e xo s c i l l a t i o n s i n t h e
d y n a m i c s o f c o u p l
e d c h a Ｇo t i c s y s t e m s u s i n g t r a n s i e n t s [J ]．P h y
s i c aA :S t a t i s Ｇt i c a l M e c h a n i c sa n dI t s A p p
l i c a t i o n s ,２０２０,５４５:１２３８１８．
[５
]V ÜL L I N G SA ,H I Z A N I D I S J ,OM E L C H E N K OI
,e t a l ．C l u s t e r e dc h i m e r as t a t e s i ns y s t e m so f t y p e ＧI e x c i t a b i l i t y [J ]．N e wJ o u r n a l o fP h y s i c s ,２０１４,１６(１２):１２３０３９．
[６
]S C HM I D T L ,S C HÖN L E B E R K ,K R I S C H E R K ,
e t a l ．C o e x i s t e n c eo fs y n c h r o n y a n di n c o h e r e n c ei n o s c i l l a t o r y m e d i a u n d e r n o n l i n e a r g l o b a lc o u p l i n g
[J ]．C h a o s :A nI n t e r d i s c i p l i n a r y J
o u r n a l o fN o n l i n Ｇe a r S c i e n c e ,２０１４,２４(１):０１３１０２．
[７]T AMA K IM ,B A N GJ W ,WA T A N A B E T ,e t a l ．N i g h tw a t c hi no n eb r a i nh e m i s p h e r ed u r i n g s l e e p
a s s o c i a t e dw i t h t h e f i r s t Ｇn i g h t e f f e c t i nh u m a n s [J ]．C u r r e n tB i o l o g y ,２０１６,２６(９):１１９０－１１９４．[８]UH L HA A SPJ ,S I N G E R W．N e u r a l s y n c h r o n y i
n b r a i nd i s o r d e r s :r e l e v a n c e f o r c o g n i t i v ed y
s f u n c t i o n s a n d p a t h o p h y s i o l o g y [J ]．N e u r o n ,２００６,５２(１):１５５－１６８．
[９]F A ND G ,W A N G Q Y ．S y n c h r o n i z a t i o na n db u r s t i n g
t r a n s i t i o n o f t h e c o u p l e dH i n d m a r s h ＧR o s e s y s t e m sw i t h a s y m m e t r i c a l t i m e Ｇd e l a y s [J ]．S c i e n c e C h i n aT e c h n o l o g
i Ｇc a lS c i e n c e ,２０１７,６０(７):１０１９－１０３１．
[１０]Y A N G Y ,X I A N GC ,D A IX ,e t a l ．C h i m e r a s t a t e s
a n d c l u s t e r s o l u t i o n s i nH i n d m a r s h ＧR o s e n e u r a l n e t Ｇw o r k s w i t hs t a t er e s e t t i n g p r o c e s s [J ]．C o g
n i t i v e ２
７
第９期许一凡等:局部耦合的双层忆阻神经元网络的嵌合态和同步
N e u r o d y n a m i c s,２０２２,１６(１):２１５－２２８．[１１]Z HA N G Y,WA N GP,Z HA N G M,e t a l．C h i m e r a
d y n a m i c s i nn o n l o c a l l y c o u p l
e db i c o m p o n e n t o s c i l l aＧ
t o r s[J]．E u r o p h y s i c sL e t t e r s,２０２１,１３５:４０００４．[１２]C A L I M A,T O R R E SJ J,O Z E R M,e t a l．C h i m e r a s t a t e si n h y b r i d c o u p l e d n e u r o n p o p u l a t i o n s[J]．
N e u r a lN e t w o r k s,２０２０,１２６:１０８－１１７．
[１３]HU OS,T I A N C,K A N G L,e t a l．C h i m e r a s t a t e s o f n e u r o n n e t w o r k s w i t h a d a p t i v e c o u p l i n g[J]．
N o n l i n e a rD y n a m i c s,２０１９,９６(１):７５－８６．[１４]H I Z A N I D I SJ,K O U V A R I S N E,Z A M O R AＧLóP E Z G,e ta l．C h i m e r aＧl i k es t a t e si n m o d u l a rn e u r a ln e tＧ
w o r k s[J]．S c i e n t i f i cR e p o r t s,２０１６,６(１):１－１１．[１５]B E R AB K,G HO S H D,B A N E R J E E T．I m p e r f e c t t r a v e l i n g c h i m e r as t a t e si n d u c e d b y l o c a ls y n a p t i c
g r a d i e n t c o u p l i n g[J]．P h y s i c a lR e v i e w E,２０１６,９４
(１):０１２２１５．
[１６]S H E P E L E VIA,B U K H A V,S T R E L K O V A GI,
e t a l．C h i m e r as t a t e si ne n s e m b l e so fb i s t a b l ee l eＧ
m e n t sw i t hr e g u l a r a n dc h a o t i cd y n a m i c s[J]．N o nＧ
l i n e a rD y n a m i c s,２０１７,９０(４):２３１７－２３３０．[１７]K A S A T K I NDV,Y A N C HU KS,S C H L LE,e t a l．S e l fＧo r g a n i z e d e m e r g e n c e o fm u l t i l a y e r s t r u c t u r e a n d
c h i m e r a s t a t e s i n
d y n a m i c a l n
e t w o r k sw i t ha d a p t i v e
c o u p l i n g s[J]．P h y s i c a lR e v i e w E,２０１７,９６(６):
０６２２１１．
[１８]A L B E R TR,B A R A BÁS IAL．S t a t i s t i c a lm e c h a n i c s o f c o m p l e xn e t w o r k s[J]．R e v i e w s o fM o d e r nP h y sＧ
i c s,２００２,７４(１):４７．
[１９]N E WMA N M EJ．T h es t r u c t u r ea n df u n c t i o no f
c o m p l e xn e t w o r k s[J]．S I AM R e v i e w,２００３,４５
(２):１６７－２５６．
[２０]WA N GDL．E m e r g e n t s y n c h r o n y i n l o c a l l y c o u p l e d n e u r a l o s c i l l a t o r s[J]．I E E E T r a n s a c t i o n so nN e u r a l
N e t w o r k s,１９９５,６(４):９４１－９４８．
[２１]B E R A B K,G HO S H D．C h i m e r as t a t e s i n p u r e l y l o c a l d e l a yＧc o u p l e do s c i l l a t o r s[J]．P h y s i c a lR e v i e w
E,２０１６,９３(５):０５２２２３．
[２２]R U Z Z E N EG,OM E L C H E N K OI,S AW I C K I J,e t
a l．R e m o t e p a c e m a k e rc o n t r o lo fc h i m e r as t a t e si n
m u l t i l a y e r n e t w o r k s o f n e u r o n s[J]．P h y s i c a l R e v i e w
E,２０２０,１０２(５):０５２２１６．
[２３]MA J H I S,P E R C M,G HO S H D．C h i m e r a s t a t e s i n u n c o u p l e dn e u r o n s i n d u c e db y am u l t i l a y e r s t r u c t u r e
[J]．S c i e n t i f i cR e p o r t s,２０１６,６:３９０３３．[２４]王松,茅晓晨．含时滞的忆阻耦合H R神经元的复杂放电行为[J]．动力学与控制学报,２０２０,１８(１):３３－
３９．
WA N GS,MA OXC．C o m p l e x f i r i n g b e h a v i o r s o f a
m e m r i s t i v eH Rn e u r o n sw i t h t i m ed e l a y s[J]．J o u rＧ
n a l o fD y n a m i c sa n dC o n t r o l,２０２０,１８(１):３３－３９．
(i nC h i n e s e)
[２５]马杰,高洁,独盟盟等．时滞下忆阻突触耦合
H o p f i e l d神经网络的动力学行为分析[J]．动力学与
控制学报,２０２１,１９(６):５９－６６．
MAJ,G A OJ,D U M M,e t a l．D y n a m i cb e h a v i o r
a n a l y s i so f m e m r i s t i v es y n a p t i cc o u p l i n g H p o f i l e d
n e u r a l n e t w o r ku n d e r t i m e d e l a y[J]．J o u r n a l o fD yＧ
n a m i c s a n dC o n t r o l,２０２１,１９(６):５９－６６．(i n C h iＧ
n e s e)
[２６]L I U Y,S U NZ,Y A N GX,e t a l．D y n a m i c a l r o b u s tＧn e s s a n d f i r i n g m o d e s i nm u l t i l a y e rm e m r i s t i v en e uＧ
r a ln e t w o r k so fn o n i d e n t i c a ln e u r o n s[J]．A p p l i e d
M a t h e m a t i c s a n dC o m p u t a t i o n,２０２１,４０９:１２６３８４．[２７]X UF,Z HA N GJ,J I N M,e t a l．C h i m e r a s t a t e s a n d s y n c h r o n i z a t i o n b e h a v i o ri n m u l t i l a y e r m e m r i s t i v e
n e u r a l n e t w o r k s[J]．N o n l i n e a rD y n a m i c s,２０１８,９４
(２):１－９．
[２８]G O P A LR,C HA N D R A S E K A R V K,V E N K A T EＧS A N A,e ta l．O b s e r v a t i o na n dc h a r a c t e r i z a t i o no f
c h i m e r as t a t e si nc o u p l e
d d y n a m i c a ls y s t
e m s w i t h
n o n l o c a l c o u p l i n g[J]．P h y s i c a lR e v i e wE,２０１４,８９
(５):０５２９１４．
[２９]O M E L C H E N K OI,O M E L C H E N K O E,HÖV E LP,
e ta l．W h e nn o n l o c a l c o u p l i n g b e t w e e no s c i l l a t o r sb eＧ
c o m e s s t r o n g e r:p a t c h e
d s y n c h r o n y o r m u l t i c h i m
e r a
s t a t e s[J]．P h y s i c a l R e v i e wL e t t e r,２０１３,１１０:２２４１０１．[３０]OM E L C H E N K O I,MA I S T R E N K O Y,HÖV E L P,e t a l．L o s so f c o h e r e n c e i nd y n a m i c a l n e t w o r k s:
s p a t i a l c h a o sa n dc h i m e r as t a t e s[J]．P h y s i c a lR eＧ
v i e wL e t t e r s,２０１１,１０６(２３):２３４１０２．
[３１]R A K S H I TS,MA J H IS,B E R A B K,e t a l．T i m eＧv a r y i n g m u l t i p l e xn e t w o r k:I n t r a l a y e r a n d i n t e r l a y e r
s y n c h r o n i z a t i o n[J]．P h y s i c a lR e v i e w E,２０１７,９６
(６):０６２３０８．
[３２]Y U A N Y Y,Y A N G H,HA N F,e t a l．T r a v e l i n g
c h i m e r as t a t e s i nl o c a l l y c o u p l e
d m
e m r i s t i v e H i n dＧ
m a r s hＧR o s en e u r o n a ln e t w o r k sa n dc i r c u i ts i m u l aＧ
t i o n[J]．S c i e n c e C h i n a T e c h n o l o g i c a l S c i e n c e s,
２０２２,６５(７):１４４５－１４５５．
３７。