1、第4 9卷 第3期2 0 2 3年9月延 边 大 学 学 报(自然科学版)J o u r n a l o f Y a n b i a n U n i v e r s i t y(N a t u r a l S c i e n c e E d i t i o n)V o l.4 9 N o.3S e p.2 0 2 3收稿日期:2 0 2 3 0 5 1 8基金项目:国家自然科学基金面上项目(1 1 6 7 2 0 7 4);福建省自然科学基金(2 0 2 2 J 0 1 6 5 7)第一作者:曾巧云(1 9 9 9),女,硕士研究生,研究方向为非线性神经动力学.通信作者:郑艳红(1 9 7 7
2、),女,博士,教授,研究方向为非线性神经动力学.文章编号:1 0 0 4-4 3 5 3(2 0 2 3)0 3-0 2 3 0-0 6一类具有时滞的帕金森病模型异常振荡的分岔分析曾巧云,郑艳红,易丹(福建师范大学 数学与统计学院;福建省分析数学及应用重点实验室:福州 3 5 0 1 1 7)摘要:研究了一个拓展的具有时滞的丘脑底核 苍白球网络模型存在振荡的理论条件,并利用R o u t h-H u r w i t z定理推导了该模型在平衡点的稳定性.数值模拟验证表明,所得的理论条件成立,且该模型产生分岔的临界点与理论的分岔点吻合度较高.当神经元集群间的传输时滞较小时,系统处于健康状态;当时滞
3、较大时,系统会发生过度的b e t a振荡,即系统会处于帕金森病状态;神经元集群间的连接权值也可对系统的振荡产生影响.关键词:帕金森病模型;振荡;H o p f分岔;时滞;连接权值;S TN-G P网络中图分类号:O 1 7 5;O 3 2 2 文献标志码:AB i f u r c a t i o n a n a l y s i s o f a b n o r m a l o s c i l l a t i o n s i n a c l a s s o f P a r k i n s o n s d i s e a s e m o d e l w i t h t i m e d e l a
4、yZ E NG Q i a o y u n,Z HE NG Y a n h o n g,Y I D a n(C o l l e g e o f M a t h e m a t i c s a n d S t a t i s t i c s,F u j i a n N o r m a l U n i v e r s i t y;F u j i a n K e y L a b o r a t o r y o f M a t h e m a t i c a l A n a l y s i s a n d A p p l i c a t i o n s:F u z h o u 3 5 0 1 1 7,
5、C h i n a)A b s t r a c t:I n t h i s p a p e r,t h e t h e o r e t i c a l c o n d i t i o n s o f o s c i l l a t i o n s i n a n e x t e n d e d s u b t h a l a m i c n u c l e u s-g l o b u s p a l l i d u s n e t w o r k(S T N-G P)m o d e l w i t h t i m e d e l a y w e r e s t u d i e d.A n d
6、t h e s t a b i l i t y o f t h e m o d e l a t e q u i l i b r i u m p o i n t w a s d e d u c e d b y R o u t h-H u r w i t z t h e o r e m.N u m e r i c a l s i m u l a t i o n r e s u l t s s h o w t h a t t h e t h e o r e t i c a l c o n d i-t i o n s a r e v a l i d,a n d t h e c r i t i c a
7、l p o i n t s o f b i f u r c a t i o n g e n e r a t e d b y t h e m o d e l a r e i n g o o d a g r e e m e n t w i t h t h e t h e o r e t i c a l b i f u r c a t i o n p o i n t s.T h e s y s t e m i s i n a h e a l t h y s t a t e w h e n t h e t r a n s m i s s i o n d e l a y b e t w e e n n
8、e u r o n c l u s t e r s i s s m a l l.T h e l a r g e r t i m e d e l a y m a k e s t h e s y s t e m h a v e e x c e s s i v e b e t a o s c i l l a t i o n,t h a t i s,t h e s y s t e m i s i n a P a r k i n s o n s s t a t e.T h e c o n n e c t i o n w e i g h t s b e t w e e n c l u s t e r s
9、o f n e u r o n s c a n a l s o a f f e c t t h e o s c i l l a t i o n o f t h e s y s t e m.K e y w o r d s:P a r k i n s o n s d i s e a s e m o d e l;o s c i l l a t i o n;H o p f b i f u r c a t i o n;t i m e d e l a y;c o n n e c t i o n w e i g h t;s u b t h a-l a m i c n u c l e u s-g l o b
10、u s p a l l i d u s n e t w o r k0 引言帕金森病(P D)是一种运动功能障碍疾病,其主要是由基底神经节中的丘脑底核(S T N)和苍白球(G P)之间的神经元因产生过度的b e t a振荡而引起的.为了研究该疾病,学者们提出了多种与P D密切 第3期曾巧云,等:一类具有时滞的帕金森病模型异常振荡的分岔分析相关的神经元振荡模型.例如:T e r m a n等1提出了一种S T N苍白球外侧(G P e)神经元网络模型,并对模型的不同活动模式进行了研究.H o l g a d o等2提出了一个平均放电率模型,并推导了该模型发生振荡的边界条件;其探究还发现,改变连接
11、权值可以使系统处于P D状态和出现H o p f分岔.H u等3-4对S T N-G P e网络的振荡频带进行了分析,研究显示连接权值和时滞都是系统产生振荡的关键因素,且选取不同的连接权值还会使系统产生双向H o p f分岔现象.赵静仪5将G P分为内侧(G P i)和G P e,并由此构建了一个新的P D模型,研究显示该模型中的连接权值和时滞与P D振荡密切相关.由于神经元集群具有异质性,因此一些学者对G P e进行了进一步分类,如G a s t等6将G P e分为典型细胞(G P e-p)和烷基苍白球细胞(G P e-a).此外,研究还显示大脑中的其他核团(如脚桥核)7也可对S T N-G
12、 P网络的振荡产生影响.基于上述研究,本文使用参数较少的平均放电率模型(将时滞作为分岔参数)研究了扩展的S T N-G P网络模型发生H o p f分岔的动力学机制.1 模型及其平衡点的稳定性图1 扩展的S T N-G P网络模型示意图扩展的S T N-G P网络模型如图1所示,其中实线箭头表示集群接收的兴奋性刺激,虚线箭头表示集群受到的抑制性刺激.S T N 1和S T N 2表示S T N在不同环境下的两个核团,其膜时间常数为S.基于模型本身固有的复杂性,本文将神经元集群间的传输时滞Ti j和膜时间常数i设为相同(即令Ti j=T,i=),其中i,j=S,G.图1中系统的动力学行为用平均放
13、电率方程可表示为:dS1dt=1(FS(-wG SG(t-T)+wC SC)-S1(t),dS2dt=1(FS(-wG SG(t-T)-S2(t),dGdt=1(FG(wS GS1(t-T)+wS GS2(t-T)-wG GG(t-T)-wX GX)-G(t).(1)其中:S1(t)、S2(t)和G(t)分别为S T N 1、S T N 2和G P的放电率;wi j为神经元集群i到j的连接权值;wC S和wX G分别表示皮层(C t x)对集群S T N 1和纹状体(S t r)对G P的输入强度;X和C分别表示纹状体和大脑皮层的输入常数;Fi(x)(Fi(x)0)为神经元集群i的激活函数,其
14、表达式为 Fi(x)=Mi1+Mi-BiBie x p(-4x/Mi),(2)其中Bi表示无外部输入时神经元集群i的放电率,Mi表示神经元集群i的最大放电率.假设系统的平衡点为(S01,S02,G0).下面讨论系统在该点的稳定性,即讨论系统存在H o p f分岔的条件.首先,对平衡点进行线性变换(S1=S1-S01,S2=S2-S02,G=G-G0),使系统的平衡点(S01,S02,G0)平移到原点(0,0,0)处,则系统(1)在平衡点处可线性化为如下矩阵形式:S 1(t)S 2(t)G(t)=B1S1(t)S2(t)G(t)+B2S1(t-T)S2(t-T)G(t-T).(3)其中:B1=-
15、1/000-1/000-1/,B2=00-a1 300-a2 3a3 1a3 2-a3 3 ,132延边大学学报(自然科学版)第4 9卷 a1 3=4wG S/(S01/MS-(S01/MS)2),a2 3=4wG S/(S02/MS-(S02/MS)2),a3 1=a3 2=4wS G/(G0/MG-(G0/MG)2),a3 3=4wG G/(G0/MG-(G0/MG)2).由拉普拉斯变换可知,系统(3)对应的特征方程为:I-B1-B2e-T=(+1)(+1)2+a3 3(+1)e-T+(a1 3a3 1+a2 3a3 2)e-2 T=0.(4)由式(4)和稳定性定理8可知,当系统的特征多项
16、式的特征根都有负实部时,该系统是稳定的.另外,由于式(4)有一个特征根(1=-1/)为负根,因此判断系统(1)是否稳定只需考虑如下方程(5)即可:(+1)2+a3 3(+1)e-T+(a1 3a3 1+a2 3a3 2)e-2 T=0.(5)令M=a1 3a3 1+a2 3a3 2,于是方程(5)可改写为:2+2+12+a3 3(+1)e-T+Me-2 T=0.(6)当T=0时,式(6)变为:2+(2+a3 3)+(12+a3 3+M)=0.(7)由R o u t h-H u r w i t z定理8可知,如果式(7)满足条件H 1(2/+a3 30,1/2+a3 3/+M0),则式(7)的根
17、都有负实部,故系统(1)在平衡点(S01,S02,G0)处是渐进稳定的.若条件H 1不成立,则式(7)至少有一个根没有负实部,此时系统(1)在平衡点(S01,S02,G0)处是不稳定的.当T0时,假定式(6)有一对复共轭纯虚根=i(0).将这对复共轭纯虚根代入式(6)中可得:(M-2+12)c o s(T)-2s i n(T)=-a3 3,2c o s(T)+(-M-2+12)s i n(T)=-a3 3.(8)情形1 当方程组(8)的系数行列式为0,即 M-2+12-22-M-2+12=0(9)时,如果条件H 2(M-1/20)成立,则由式(9)可得1=M-1/2.由于方程组(8)中的2个式
18、子是线性相关的,因此只需将1代入方程组(8)中的任一个式子即可.本文选取第1个式子.将第1个式子与c o s2(1T)+s i n2(1T)=1进行联立和整理可得:(221+1)c o s2(1T)+a3 3c o s(1T)+(a23 324-221)=0.(1 0)由于1=M-1/2为实数,且c o s(1T)为实值,因此式(1 0)须满足条件H 3(4(221+1)-a23 320).此时,对式(1 0)进行求解可得:c o s(1T)=-a3 3 14(221+1)-a23 322(221+1).(1 1)再由式(1 1)可知Tk1,2=11a r c c o s(-a3 3 14(2
19、21+1)-a23 322(221+1)+2k1,k=0,1,n.由上述可知,当方程组(8)的系数行列式为0时,系统(1)产生H o p f分岔的临界时滞为T0=m i nTkl|l=1,2;232 第3期曾巧云,等:一类具有时滞的帕金森病模型异常振荡的分岔分析k=0,1,n,且此时式(6)对应的纯虚根为i0.情形2 当方程组(8)的系数行列式不为0时,对其求解可得:c o s(T)=-a3 31+M 2+22,s i n(T)=a3 321+M 2+22.(1 2)将式(1 2)代入c o s2(T)+s i n2(T)=1中可得:4+(2M-a23 3+22)2+(M2+2M-a23 32
20、+14)=0.(1 3)令p=2,则式(1 3)可变为:p2+(2M-a23 3+22)p+(M2+2M-a23 32+14)=0.(1 4)假设条件H 4(M+1/)2a23 34M且M-1/20)可得:2,3=(a23 3-2M)2a3 32a23 3-4M-222.(1 5)再由式(1 2)可知Tkj+1=1ja r c c o s(-a3 31+M 2+22)+2kj,k=0,1,n;j=2,3.由上述可知,当方程组(8)的系数行列式不为0时,系统(1)产生H o p f分岔的临界时滞为T0=m i nTkl|l=3,4;k=0,1,n,且方程(6)对应的纯虚根为i0.下面对参数T的横
21、截条件进行验证.将隐函数(T)代入式(6)中后对T进行求导可得:ddT -1=(2 2+2)e T+a3 32 M 2e-T-(32+22+)e T-T.由上式及横截条件的定义可知,如果条件H 5(R eddT -1=i0 0)成立,则参数T的横截条件成立.定理1当系统(1)满足情形1时,若条件H 1H 3、H 5成立(或当系统(1)满足情形2时,若条件H 1、H 4、H 5成立),则系统(1)在平衡点(S01,S02,G0)处有如下结论:1)当T0,T0)时,系统(1)在平衡点(S01,S02,G0)处是渐近稳定的;2)当TT0时,系统(1)在平衡点(S01,S02,G0)处是不稳定的;3)
22、当T=T0时,系统(1)在平衡点(S01,S02,G0)处产生H o p f分岔.2 数值模拟取系统(1)的参数为:MS=3 0 0s p k/s,MG=4 0 0s p k/s,BS=1 7s p k/s,BG=7 5s p k/s,C=2 7s p k/s,X=2 s p k/s,=0.0 0 6 s,wS G=1 9,wG S=1.1 2,wG G=6.6,wC S=2.4 2,wX G=1 5.1.对系统(1)进行计算可得:当T=0时,条件H 1成立;当T0时,方程组(8)的系数行列式为0,且情形1中的条件H 2、H 3、H 5都成立.由此可得系统(1)在理论上产生H o p f分岔的
23、临界时滞为T0=0.0 0 7 6 s.2.1 时滞对系统振荡的影响图2为G P放电率与时滞T的关系.由图2可以看出:当T逐渐增大到0.0 0 7 4 s时,G P的放电率开始发生分岔,且该临界点(0.0 0 7 4 s)与理论上发生分岔的临界点高度吻合.当神经元集群间的传输时滞未达到临界点时,系统处于稳定状态(S S);当时滞大于产生分岔的临界点时,系统处于振荡状态(O S),且时滞越大系统的振幅也越大.332延边大学学报(自然科学版)第4 9卷 图3为神经元集群G P的主峰频率(D F)与T的关系.由图3可以看出,当T0,T0)时,系统对应的主峰频率为0;当TT0时,系统的振荡频率主要发生
24、在1 33 0H z的范围内.由此可知,该现象为典型的b e t a振荡.图2 G P放电率与时滞T的关系 图3 G P的主峰频率与T的关系 图4 不同传输时滞T所对应的时间序列图4为不同传输时滞T所对应的时间序列.由图4(a)可以看出,当T=0.0 0 5 s时,系统(1)在初始阶段产生了轻微振荡,随后逐渐进入了稳定状态,这表明系统的正平衡点是一个渐近稳定点.由图4(b)可以看出,当T=0.0 0 8 s时,系统出现了高频的周期振荡,这表明此时系统处于P D状态.由图4(c)可以看出,当T=0.0 1 5s时,系统出现了低频的周期振荡,且振幅较大,这表明增大神经核团间的传输时滞会增加P D的
25、振荡程度.2.2 连接权值对系统振荡的影响图5为不同的wG G对系统发生振荡的影响.由图5(a)可以看出,当wG G取值在1.3,6.9 范围内时,wG G值越大系统发生振荡所需的时滞越大,即此时wG G值越大系统的稳定区间越大.由图5(b)可以看出,当取T=0.0 0 8 s和wG G1.3时,G P的振幅随wG G值的增大而增大;但当wG G值到达1.3时,G P的振幅出现减小趋势.这说明,较小或较大的wG G值有助于降低系统的P D振荡程度.(a)振荡边界 (b)振荡振幅图5 不同的wG G对系统振荡的影响图6为G P的振幅和频率在wS G和wG S共同作用下的变化情况.由图6(a)可以
26、看出:当wG S较大、wS G较小时,G P的振幅很小,且其放电率几乎保持不变,即此时系统基本处于稳定状态;当wG S越小、wS G越大时,G P的振幅越大,即此时系统处于振荡不断加强的状态.由图6(b)可以看出,系统在wS G和wG S的共同作用下,其发生的振荡主要为2 53 0H z的b e t a振荡.432 第3期曾巧云,等:一类具有时滞的帕金森病模型异常振荡的分岔分析(a)G P的振幅变化情况 (b)G P的频率变化情况图6 G P的振幅和频率在wS G和wG S共同作用下的变化情况3 结论本文研究了一个扩展的具有传输时滞的S T N-G P网络模型,给出了该模型正平衡点稳定性的判别
27、方法和模型发生P D振荡的临界条件.研究表明,时滞会影响系统的稳定性.当时滞T在0,T0)范围内变化时,系统处于稳态;当T=T0时,系统产生H o p f分岔,即此时模型开始从健康状态转变为帕金森状态.通过调优神经元集群间的连接权值可以减小系统的振荡振幅,即可以改善系统的P D振荡程度.本文研究结果有助于人们更好地了解P D振荡的发病机制,并可为临床诊断提供参考.在今后的研究中,我们将考虑具有不同时滞的神经元集群对系统P D振荡的影响.参考文献:1 T E RMAN D,R U B I N J E,Y EW A C,e t a l.A c t i v i t y p a t t e r n s
28、 i n a m o d e l f o r t h e s u b t h a l a m o p a l l i d a l n e t w o r k o f t h e b a s a l g a n g l i aJ.T h e J o u r n a l o f N e u r o s c i e n c e,2 0 0 2,2 2(7):2 9 6 3-2 9 7 6.2 HO L GA D O A J N,T E R R Y J R,B OG A C Z R.C o n d i t i o n s f o r t h e g e n e r a t i o n o f b e
29、t a o s c i l l a t i o n s i n t h e s u b t h a l a m i c n u c l e u s-g l o b u s p a l l i d u s n e t w o r kJ.T h e J o u r n a l o f N e u r o s c i e n c e,2 0 1 0,3 0(3 7):1 2 3 4 0-1 2 3 5 2.3 HU B,XU M B,WANG Z Z,e t a l.T h e t h e o r e t i c a l m e c h a n i s m o f P a r k i n s o n
30、s o s c i l l a t i o n f r e q u e n c y b a n d s:a c o m p u-t a t i o n a l m o d e l s t u d yJ.C o g n i t i v e N e u r o d y n a m i c s,2 0 2 1,1 5(4):7 2 1-7 3 1.4 HU B,XU M B,Z HU L Y,e t a l.A b i d i r e c t i o n a l H o p f b i f u r c a t i o n a n a l y s i s o f P a r k i n s o ns
31、o s c i l l a t i o n i n a s i m p l i f i e d b a s a l g a n g l i a m o d e lJ.J o u r n a l o f T h e o r e t i c a l B i o l o g y,2 0 2 2,5 3 6:1 1 0 9 7 9.5 赵静仪.帕金森病的基底神经节网络动力学研究D.呼和浩特:内蒙古大学,2 0 2 0.6 G A S T R,G ON G R,S CHM I D T H,e t a l.O n t h e r o l e o f a r k y p a l l i d a l a n
32、d p r o t o t y p i c a l n e u r o n s f o r p h a s e t r a n s i t i o n s i n t h e e x t e r n a l p a l l i d u mJ.T h e J o u r n a l o f N e u r o s c i e n c e,2 0 2 1,4 1(3 1):6 6 7 3-6 6 8 3.7 陈国泰,郑艳红.丘脑底核 苍白球外侧 脚桥核网络的振荡动力学分析J.福建师范大学学报(自然科学版),2 0 2 2,3 8(3):1 1-1 6.8 马知恩,周义仓,李承治.常微分方程定性与稳定性方法M.第2版.北京:科学出版社,2 0 1 5:7 4-8 1.532