具有耦合时滞的模块神经元网络簇同步的抑制方法技术

本发明专利技术公开了一种具有耦合时滞的模块神经元网络簇同步的抑制方法，包括如下步骤：建立模块神经元网络模型；利用二维离散映射描述模块神经元网络中节点上的局部神经元模型；利用单个神经元的快变量和整个网络平均场的时间历程，生动地描述神经元的放电动力学；计算序参数R，确定耦合强度对诱导簇同步的影响；提出利用微分反馈控制技术去抑制簇同步，计算抑制参数S，确定微分反馈控制技术对减小或抑制簇同步的有效性。本发明专利技术发现较大的耦合强度能够诱导模块神经元网络的簇同步，并且当微分反馈控制的控制参数位于控制域内时，微分反馈控制能够有效地抑制簇同步。当引入耦合时滞时，较小的耦合时滞能够提高微分反馈控制抑制簇同步的效果。

【技术实现步骤摘要】
具有耦合时滞的模块神经元网络簇同步的抑制方法
本专利技术涉及生物领域，具体涉及一种具有耦合时滞的模块神经元网络簇同步的抑制方法。
技术介绍
神经元的同步放电活动对生物神经信息处理具有关键性作用，且与神经系统许多正常的生理机制相关。然而，最近的医学实验表明并不是神经元所有的同步模式都能促进脑功能的正常运转。神经元异常的同步化放电活动会严重破坏大脑正常的功能和结构，进而导致某些神经性生理疾病，如帕金森症、癫痫症、手的颤抖等。因此，研究如何抑制神经元异常的同步放电活动具有重要的理论意义和潜在的应用价值。临床上，消除神经元病态同步有效的一种方法是深度脑刺激(DBS)。该方法在大脑的特定区域中植入电极，通过施加高频脉冲刺激，从而控制被刺激区域神经元的同步放电活动。然而，此方法治疗神经系统疾病的生物物理机制仍不十分明确，并且该技术在应用中存在一定的局限性。由于DBS技术存在以上不足，根据非线性动力学和统计物理的控制理论和方法，近年来学者们提出了一些消除神经元间病态同步的其它方法。例如，非线性时滞反馈技术，已经实现了由非线性微分方程组成的神经元网络的同步抑制，其效果不依赖于刺激参数的改变；滤波器辅助的平均场反馈方法，可以有效地控制神经元网络的同步活动；局部周期信号刺激方法，也可用于抑制神经元网络中的病态同步放电活动。在上述抑制神经元同步的研究过程中，如何在保持神经元本身放电特性不变的条件下，有效地消除神经系统的病态同步放电，这是一直令人困扰的一个问题。德国学者Rosenblum和Pikowsky根据Pyragas的微分反馈控制思想，提出了利用线性时滞反馈信号来消除大脑中的病态同步节律。该方法的优点是不需要知道各个神经元的确切参数，只需测得神经元集群的节律，将延迟信号反馈给被刺激神经元即可，并且这种微分反馈控制方法是一种非侵害性控制方法，它并不改变单个神经元的动力学性质。目前，微分反馈控制技术初步用于全局耦合神经元网络和小世界神经元网络的同步抑制中。最新关于猕猴和猫大脑皮层的解剖实验证实，具有相似连接特征或功能特征的神经元往往聚集在同一个模块中，这表明脑皮层的神经元网络具有模块结构。同时，由于突触间隙的存在及动作电位沿轴突的有限传播速度，时滞在神经元网络中是普遍存在的。因而，具有耦合时滞的模块神经元网络是对真实大脑神经系统更贴切的描述。然而，如何抑制耦合时滞作用下模块神经元网络的病态同步，目前还没有见到报道。发展抑制含有耦合时滞的模块神经元网络簇同步的方法与技术迫在眉睫。
技术实现思路
为解决上述问题，本专利技术鉴于微分反馈控制技术的优点，提供了一种对具有耦合时滞的模块神经元网络簇同步的抑制方法。为实现上述目的，本专利技术采取的技术方案为：具有耦合时滞的模块神经元网络簇同步的抑制方法，其特征在于，包括如下步骤：S1、建立模块神经元网络模型；S11、假设模块网络存在M个子网络，每个子网络含有N个节点，每个子网络是小世界网络，每个子网络中的每个神经元与它最邻近的和次邻近的节点相连；S12、利用一个给定的概率pintra增加子网络内节点间的连接；S13、选择来自第I(I＝1，...，M)和第J(J≠I，J＝1，...，M)个子网络中的一些节点对，并且用子网络间连接概率pinter连接选定的节点对，生成一个模块神经元网络；S2、利用二维离散映射描述模块神经元网络中节点上的局部神经元模型，神经元网络的表达式如下：yI，i(n+1)＝yI，i(n)+c[xI，i(n)-E]；此时u(n)是微分反馈控制项，形式为：u(n)＝K(X(n-τ)-X(n))，其中K是反馈强度，τ是反馈时滞，X是网络的平均场，具体为是耦合项，形式为：式中，(I，i)为第I个子网络中的第i个神经元，i＝1，...，N；I＝1，...，M；x为神经元膜电位的动力学；y为所有的离子电流通过神经元膜渗入神经元内的集体行为，具有修复作用；ε为不同神经元间的全局耦合强度，并且τ1是耦合时滞；AI，J(i，j)为网络的连接矩阵；当第I个子网络中的第i个神经元与第J个子网络中的第j个神经元相连接时，AI，J(i，j)＝AJ，I(j，i)＝1，否则，AI，J(i，j)＝AJ，I(j，i)＝0，并且AI，I(i，i)＝0；S3、利用单个神经元的快变量和整个网络平均场的时间历程，生动地描述神经元网络的动力学；S4、通过以下公式计算序参数R，确定耦合强度对簇同步的影响：式中，表示第J个子网络中第j个神经元在时刻n处的簇相位，并且簇相位的描述如下：式中，nJ，j，k是第J个子网络中第j个神经元的第k个簇开始放电的时刻；一般地，R的值越大，表明神经元网络簇同步的程度越大，特别地，当所有神经元簇放电达到簇同步时，簇相位几乎一致，从而R接近于1；相反地，当所有簇放电神经元处于完全不簇同步状态时，簇相位几乎不相关，从而R≈0；S5、提出利用微分反馈控制技术去抑制模块神经元网络的簇同步，通过以下公式计算抑制参数S，确定微分反馈控制机制对减小或抑制簇同步的有效性；式中，X和Xf分别是不存在和存在微分反馈控制时平均场的值。一般地，S值越大，微分反馈控制抑制簇同步越有效。本具体实施中，假设S＞3的参数区域为有效抑制簇同步的控制域。本专利技术发现较大的耦合强度能够诱导模块神经元网络的簇同步，并且当微分反馈控制的控制参数在控制域内时，本专利技术提出的微分反馈控制技术能够有效地抑制簇同步。当引入耦合中的时滞时，我们观察到较小值的耦合时滞能够促进微分反馈控制技术抑制簇同步，并且存在最优的耦合时滞使得这种抑制效果能够达到最佳。这些结果有助于理解真实神经系统的簇放电同步，并为治疗簇同步引发的一些神经性生理疾病提供一定的理论指导意义。附图说明图1为本专利技术实施例中由M＝3个小世界子网络构成的模块网络的图例，并且每个子网络中有N＝12个节点。图2为本专利技术实施例中当ε＝0[(a)和(c)]时和当ε＝0.05[(b)和(d)]时，两个随机选择的神经元膜电位x的时间历程示意图；网络模型的参数值为M＝3和N＝30。图3为本专利技术实施例中当(a)ε＝0和(b)ε＝0.05时，模块神经元网络的时空图，网络模型的参数值为M＝3和N＝30。图4(a)为本专利技术实施例中未耦合时(ε＝0)的平均场的时间历程，网络模型的参数值为M＝3和N＝30。图4(b)为强耦合(ε＝0.05)时平均场的时间历程，网络模型的参数值为M＝3和N＝30。图5为本专利技术实施例中取子网络个数不同时，序参数R随着耦合强度ε变化的曲线。图6(a)为本专利技术实施例中当存在微分反馈控制u(n)时，网络平均场随着时间的变化。图6(b)为本专利技术实施例中微分反馈控制信号u(n)的时间历程。图7为本专利技术实施例中当存在微分反馈控制时，两个随机选择的神经元膜电位x的时间历程，网络模型的参数为M＝3，N＝30，K＝0.15，τ＝120。图8为本专利技术实施例中在不同个数的子网络构成的模块网路中，抑制参数S随着控制参数组合(K和τ)的变化；图中：(a)M＝3；(b)M＝5；(c)M＝7。图9为本专利技术实施例中在不同个数的子网络构成的模块网络中，当控制时滞τ＝120时，对于不同的耦合强度，抑制参数S随着耦合时滞τ1的变化曲线；图中：(a)M＝3；(b)M＝5；(c)M＝7。图10为本专利技术实施例中当τ＝120本文档来自技高网...

【技术保护点】
具有耦合时滞的模块神经元网络簇同步的抑制方法，其特征在于，包括如下步骤：S1、建立模块神经元网络模型；S11、假设模块网络存在M个子网络，每个子网络含有N个节点，每个子网络是小世界网络，每个子网络中的每个神经元与它最邻近的和次邻近的节点相连；S12、利用一个给定的概率Pintra增加子网络内节点间的连接；S13、选择来自第I(I＝1，…，M)和第J(J≠I，J＝1，...，M)个子网络中的一些节点对，并且用子网络间连接概率Pinter连接选定的节点对，生成一个模块神经元网络；S2、利用二维离散映射描述模块神经元网络中节点上的局部神经元模型，神经元网络的表达式如下：xI,i(n+1)=xI,i(n)+F[xI,i(n)]-yI,i(n)-βH[xI,i(n)-d]+IIjsyn(n)+u(n);]]>yI，i(n+1)＝yI，i(n)+c[xI，i(n)‑E]；H(n)=1x≥00x<0;]]>F(x)=-m0xx≤Eminm1(x-a)Emin<-m0(x-1)x≥Emaxx<Emax,]]>此时Emin=am1m0+m1,Emax=m0+am1m0+m1.]]>u(n)是微分反馈控制项，形式为：u(n)＝K(X(n‑τ)‑X(n))，其中K是反馈强度，τ是反馈时滞，X是网络的平均场，具体为是耦合项，形式为：II,isyn(n)=ϵΣJΣjAI,J(i,j)(xJ,j(n-τ1)-xI,i(n));]]>式中，(I，i)为第I个子网络中的第i个神经元，i＝1，...，N；I＝1，...，M；x和y分别代表神经元的膜电位和恢复变量；ε为不同神经元间的全局耦合强度，τ1是耦合时滞；AI，J(i，j)为网络的连接矩阵；当第I个子网络中的第i个神经元与第J个子网络中的第j个神经元相连接时，AI，J(i，j)＝AJ，I(j，i)＝1，否则，AI，J(i，j)＝AJ，I(j，i)＝0，并且AI，I(i，i)＝0；S3、利用单个神经元的快变量和整个网络平均场的时间历程，生动地描述神经元网络的动力学；S4、通过以下公式计算序参数R，确定耦合强度对簇同步的影响：式中，表示第J个子网络中第j个神经元在时刻n处的簇相位，并且簇相位的描述如下：式中，nJ，j，k是第J个子网络中第j个神经元的第k个簇开始放电的时刻；S5、提出利用微分反馈控制技术去抑制模块神经元网络的簇同步，通过以下公式计算抑制参数S，确定微分反馈控制机制对减小或抑制簇同步的有效性；S=Var(X)Var(Xf),]]>式中，X和Xf分别是不存在和存在微分反馈控制时平均场的值。...

【技术特征摘要】
1.具有耦合时滞的模块神经元网络簇同步的抑制方法，其特征在于，包括如下步骤：S1、建立大脑神经系统的模块神经元网络模型；S11、假设模块网络存在M个子网络,每个子网络含有N个节点，每个子网络是小世界网络,每个子网络中的每个神经元与它最邻近的和次邻近的节点相连；S12、利用一个给定的概率pintra增加子网络内节点间的连接；S13、选择来自第I(I＝1,...,M)和第J(J≠I,J＝1,...,M)个子网络中的一些节点对,并且用子网络间连接概率pinter连接选定的节点对，生成一个模块神经元网络；S2、利用二维离散映射描述模块神经元网络中节点上的局部神经元模型，神经元网络的表达式如下：yI,i(n+1)＝yI,i(n)+c[xI,i(n)-E]；此时u(n)是微分反馈控制项，形式为：u(n)＝K(X(n-τ)-X(n))，其中K是反馈强度，τ是反馈时滞，X是网络的平均场，是所述模块神经元网络的输出,具体为是耦合项,形式为：...

【专利技术属性】
技术研发人员：杨晓丽，胡丽萍，
申请(专利权)人：陕西师范大学，
类型：发明
国别省市：陕西;61