时滞作用下的化学突触耦合的Hindmarsh-Rose神经元同步研究

本文通过数值模拟的方法研究了化学突触耦合的Hindmarsh-Rose神经元的同步问题,以及时滞对耦合的Hindmarsh-Rose神经元同步的影响.研究发现:耦合强度能够影响化学突触的耦合神经元系统的放电活动以及耦合神经元之间的同步状态.通过调整化学突触耦合的两个Hindmarsh-Rose神经元之间的耦合强度,可以使得神经元系统从非同步状态变为同步状态.同时也发现,适当的时滞会使得两个Hindmarsh-Rose化学突触耦合神经元系统从同步状态变为非同步状态,从而破坏系统的同步活动.     

具有时滞的帕金森模型的振荡动力学分析

研究大脑基底神经节中产生异常β振荡的起源有助于分析帕金森病的致病机理.本文系统地研究了改进的皮质-基底神经节(E-I-STN-GPe-GPi)共振模型的振荡动力学.首先,通过Routh-Hurwitz准则和稳定性理论获得了该模型局部平衡点处的稳定性与Hopf分岔发生的条件,并且推导出该共振模型存在Hopf分岔的时滞参数范围.研究发现,增加突触传输时滞能够使模型产生Hopf分岔,并且诱导β振荡的产生,使系统在健康和帕金森病这两个状态之间相互转换.其次,揭示了β振荡的产生与丘脑底核相关的突触连接强度有关.数值模拟发现,当丘脑底核同时受到兴奋性神经元集群和苍白球外侧较强的促进作用时,丘脑底核产生振荡.最后,分析了与苍白球内侧有关的参数对其产生振荡的影响,研究结果发现,当较小的苍白球外侧突触连接强度和较大的突触传输时滞共同作用时,苍白球内侧更容易发生振荡,且振幅越来越大.希望本文对E-I-STN-GPe-GPi共振模型的动力学特征的研究有助于人们理解帕金森病的致病机理和揭示帕金森病异常β振荡的来源.     

电磁场下神经元模型中混合簇的分岔机制

Pre-B?tzinger复合体是新生哺乳动物呼吸节律起源的关键部位, 是呼吸节律产生的中枢. 忆阻器的功能类似于神经元突触的可塑性, 可用其模拟磁通量.本文在Butera动力学模型的基础上引入刺激电流和磁通控制忆阻器, 分别研究这两个因素对单个pre-B?tzinger复合体神经元中混合簇放电模式的影响.通过无量纲化的方法对变量进行时间尺度分析, 结果表明, 模型包含3个不同的时间尺度.通过快慢分解和分岔分析研究了神经元混合簇放电产生和转迁的动力学机制.电流和磁通量都可以影响混合簇中胞体簇的个数, 减小电流和磁通量的值, 混合簇中胞体簇的个数也会相应减少, 并使簇的类型由"fold/homoclinic"型簇放电转迁为经由"fold/homoclinic"滞后环的"Hopf/Hopf"型簇放电.双参数分岔分析表明, 随着钙离子浓度的逐渐增加, 全系统轨线在鞍结分岔曲线和同宿轨分岔曲线之间来回跃迁, 是混合簇的产生分岔机制.全系统轨线在鞍结分岔曲线和同宿轨分岔曲线之间跃迁的次数, 与混合簇中胞体簇的个数相对应.      

轮对系统的Hopf分岔研究

针对轮对系统中的非线性动力学问题, 本文基于Hopf分岔代数判据得到考虑陀螺效应的轮对系统Hopf分岔点解析表达式, 即轮对系统蛇形失稳的线性临界速度解析表达式. 基于分岔理论得到轮对系统的第一、第二Lyapunov系数表达式, 并结合打靶法分别得到不同纵向刚度下, 考虑陀螺效应与不考虑陀螺效应的轮对系统分岔图. 通过对比有无陀螺效应的轮对系统分岔图发现, 在同一纵向刚度下, 考虑陀螺效应的轮对系统线性临界速度和非线性临界速度均大于不考虑陀螺效应的轮对系统, 即陀螺效应可以提高轮对系统的运动稳定性. 基于Bautin分岔理论, 以纵向刚度和纵向速度作为参数, 分别得到考虑陀螺效应和不考虑陀螺效应的轮对系统, 从亚临界Hopf分岔到超临界Hopf分岔, 再从超临界Hopf分岔到亚临界Hopf分岔的迁移机理拓扑图. 通过对比有、无陀螺效应的轮对系统Bautin分岔拓扑图发现, 陀螺效应将改变轮对系统的退化Hopf分岔点, 但对于轮对系统Bautin分岔拓扑图的影响不大.      

一类忆阻神经元的电活动多模振荡及Hamilton 能量反馈控制

根据法拉第电磁感应定律,在离子穿越细胞膜或者在外界电磁辐射下,细胞内外的电生理环境会产生电磁感应效应,继而会影响神经元的电活动行为. 基于此,本文考虑电磁感应影响下的 Hindmarsh-Rose (HR) 神经元模型,研究了其混合模式振荡放电特征,并设计一个 Hamilton 能量反馈控制器,将其控制到不同的周期簇放电状态. 首先,通过理论分析发现磁通 HR 神经元系统的 Hopf 分岔使其平衡点的稳定性发生了改变,并产生极限环,进而研究了 Hopf 分岔点附近膜电压的放电特征. 基于双参数数值仿真发现该系统具有丰富的分岔结构,在不同的参数平面上存在倍周期分岔、伴有混沌的加周期分岔、无混沌的加周期分岔以及共存的混合模式振荡. 最后,为了有效控制膜电压的混合模式振荡,利用亥姆霍兹理论计算出磁通 HR 神经元系统的 Hamilton 能量函数并设计 Hamilton 能量反馈控制器,通过数值仿真分析了膜电压在不同反馈增益下的簇放电状态,发现该控制器能够有效地控制膜电压到不同的周期簇放电模式. 本文的研究结果为探究电磁感应下神经元的分岔结构及其能量控制领域提供了有用的理论支撑.      

MULTI-MODE OSCILLATIONS AND HAMILTON ENERGY FEEDBACK CONTROL OF A CLASS OF MEMRISTOR NEURON$^{\bf 1)}$

根据法拉第电磁感应定律,在离子穿越细胞膜或者在外界电磁辐射下,细胞内外的电生理环境会产生电磁感应效应,继而会影响神经元的电活动行为. 基于此,本文考虑电磁感应影响下的 Hindmarsh-Rose (HR) 神经元模型,研究了其混合模式振荡放电特征,并设计一个 Hamilton 能量反馈控制器,将其控制到不同的周期簇放电状态. 首先,通过理论分析发现磁通 HR 神经元系统的 Hopf 分岔使其平衡点的稳定性发生了改变,并产生极限环,进而研究了 Hopf 分岔点附近膜电压的放电特征. 基于双参数数值仿真发现该系统具有丰富的分岔结构,在不同的参数平面上存在倍周期分岔、伴有混沌的加周期分岔、无混沌的加周期分岔以及共存的混合模式振荡. 最后,为了有效控制膜电压的混合模式振荡,利用亥姆霍兹理论计算出磁通 HR 神经元系统的 Hamilton 能量函数并设计 Hamilton 能量反馈控制器,通过数值仿真分析了膜电压在不同反馈增益下的簇放电状态,发现该控制器能够有效地控制膜电压到不同的周期簇放电模式. 本文的研究结果为探究电磁感应下神经元的分岔结构及其能量控制领域提供了有用的理论支撑.     

通过分岔控制改变神经元的兴奋性类型

提出一种通过分岔控制改变神经元兴奋性类型的方法.采用一个基于washout滤波器的动态反馈控制实现对一个二维的Hindmarsh-Rose类的模型神经元的分岔动力学控制.这一模型神经元从静息态到峰放电态跨越一个不变圆上鞍结分岔(saddle-node on invariant circle,SNIC),呈现出第一类兴奋性.在该SNIC分岔前所期望的参数值处产生一个Hopf分岔,然后通过选择适当的控制器参数调节Hopf分岔的临界性.这样,模型神经元就呈现为第二类兴奋性,因此神经元兴奋性就从第一类改变成第二类.在这个控制器中,线性控制增益决定着Hopf分岔的位置,而非线性增益决定着Hopf分岔的临界性.     

时滞Rucklidge系统Hopf分岔分析与电路实现

以一类新的单时滞Rucklidge系统为分析对象,通过计算时滞系统的平衡点,分析该系统在各平衡点的稳定性和Hopf分岔的存在性,得到其发生Hopf分岔的条件。Matlab多组数值仿真验证了理论分析的正确性。基于此设计了一种可切换时滞与非时滞的混沌电路,并运用Multisim14.0进行仿真,实验结果表明,该电路可行且有效。     

实验性神经起步点自发放电的整数倍节律和随机自共振

本文介绍了实验中发现的无外界周期刺激的神经起步点放电节律随[Ca＋＋]。变化产生的整数倍节律，并用描写神经放电的理论模型（Chay模型）进行数值模拟。结果发现：在相应的参数区间，确定性模型为－Hopf分岔，无整数倍节律；在随机模型中，在Hopf分岔点附近，整数倍节律产生，该整数倍节律是通过随机自共振产生的。实验中与模型的整数倍节律处于桢的参数区间，位于周期1和阈下振荡之间：并且有相同的特征；其峰峰间期处于一个基本峰峰新时期的整数倍，峰峰间期出现频率随峰峰新时期增加呈现出指数降低。这提示，实验中的整数倍节律是通过随机自共振产生的。     

时滞Rucklidge系统Hopf分岔分析与电路实现

以一类新的单时滞Rucklidge系统为分析对象,通过计算时滞系统的平衡点,分析该系统在各平衡点的稳定性和Hopf分岔的存在性,得到其发生Hopf分岔的条件。Matlab多组数值仿真验证了理论分析的正确性。基于此设计了一种可切换时滞与非时滞的混沌电路,并运用Multisim14.0进行仿真,实验结果表明,该电路可行且有效。     

Stability and bifurcation analysis in tri-neuron model with time delay

A simple delayed neural network model with three neurons is considered. By constructing suitable Lyapunov functions, we obtain sufficient delay-dependent criteria to ensure global asymptotical stability of the equilibrium of a tri-neuron network with single time delay. Local stability of the model is investigated by analyzing the associated characteristic equation. It is found that Hopf bifurcation occurs when the time delay varies and passes a sequence of critical values. The stability and direction of bifurcating periodic solution are determined by applying the normal form theory and the center manifold theorem. If the associated characteristic equation of linearized system evaluated at a critical point involves a repeated pair of pure imaginary eigenvalues, then the double Hopf bifurcation is also found to occur in this model. Our main attention will be paid to the double Hopf bifurcation associated with resonance. Some Numerical examples are finally given for justifying the theoretical results.     

The Bifurcation Study of 1:2 Resonance in a Delayed System of Two Coupled Neurons

In this paper, we consider a delayed system of differential equations modeling two neurons: one is excitatory, the other is inhibitory. We study the stability and bifurcations of the trivial equilibrium. Using center manifold theory for delay differential equations, we develop the universal unfolding of the system when the trivial equilibrium point has a double zero eigenvalue. In particular, we show a universal unfolding may be obtained by perturbing any two of the parameters in the system. Our study shows that the dynamics on the center manifold are characterized by a planar system whose vector field has the property of 1:2 resonance, also frequently referred as the Bogdanov–Takens bifurcation with $Z_2$ symmetry. We show that the unfolding of the singularity exhibits Hopf bifurcation, pitchfork bifurcation, homoclinic bifurcation, and fold bifurcation of limit cycles. The symmetry gives rise to a “figure-eight” homoclinic orbit.     

Complex dynamics of a four neuron network model having a pair of short-cut connections with multiple delays

This paper deals with dynamic behaviors on Hopfield type of ring neural network of four neurons having a pair of short-cut connections with multiple time delays. By suitable transformation and under certain assumptions on multiple time delays, the model is reduced to four dimensional nonlinear delay differential equations with three delays. Regarding these time delays as parameters, delay independent sufficient conditions for no stability switches of the trivial equilibrium of the linearized system are derived. Conditions for stability switching with respect to one delay parameter which is not associated with short-cut connection are obtained. Hopf bifurcations with respect to two other delays which are associated with short-cut connection are also obtained. Using the normal form method and center manifold theory, the direction of the Hopf bifurcation, stability and the properties of Hopf-bifurcating periodic solutions are determined. Using numerical simulations of the nonlinear model, different rich dynamical behaviors such as quasiperiodicity, torus attractor and chaotic-bands are also observed for suitable range of three delay parameters. Lyapunov exponents are also calculated using the AnT 4.669 tool for verification of chaotic dynamics.     

Multiple scales analysis for double Hopf bifurcation with?1:3?resonance

The dynamical behavior of a general n-dimensional delay differential equation (DDE) around a 1:3 resonant double Hopf bifurcation point is analyzed. The method of multiple scales is used to obtain complex bifurcation equations. By expressing complex amplitudes in a mixed polar-Cartesian representation, the complex bifurcation equations are again obtained in real form. As an illustration, a system of two coupled van der Pol oscillators is considered and a set of parameter values for which a 1:3 resonant double Hopf bifurcation occurs is established. The dynamical behavior around the resonant double Hopf bifurcation point is analyzed in terms of three control parameters. The validity of analytical results is shown by their consistency with numerical simulations.     

Using delay to quench undesirable vibrations

A van der Pol type system with delayed feedback is explored by employing the two variable expansion perturbation method. The perturbation scheme is based on choosing a critical value for the delay corresponding to a Hopf bifurcation in the unperturbed ε=0 system. The resulting amplitude–delay relation predicts two Hopf bifurcation curves, such that in the region between these two curves oscillations will be quenched. The perturbation results are verified by comparison with numerical integration.     

Study of mixed-mode oscillations in a parametrically excited van der Pol system

This paper investigates the effects of slowly varying parametric excitation on the dynamics of van der Pol system. Periodic bifurcation delay behaviors are exhibited when the parametric excitation slowly passes through Hopf bifurcation value of the controlled van der Pol system. The first bifurcation delay behavior relies on initial conditions, while the bifurcation delay behaviors that follow the first one are immune to initial conditions. These bifurcation delay behaviors result in a hysteresis loop between the spiking attractor and the rest state, which is responsible for the generation of mixed-mode oscillations. Then an approximate calculation for the number of spikes in each cluster of repetitive spiking of mixed-mode oscillations is explored based on bifurcation delay behaviors. Theoretical results agree well with numerical simulations.     

Double Hopf bifurcation of time-delayed feedback control for maglev system

This paper undertakes an analysis of a double Hopf bifurcation of a maglev system with time-delayed feedback. At the intersection point of the Hopf bifurcation curves in velocity feedback control gain and time delay space, the maglev system has a codimension 2 double Hopf bifurcation. To gain insight into the periodic solution which arises from the double Hopf bifurcation and the unfolding, we calculate the normal form of double Hopf bifurcation using the method of multiple scales. Numerical simulations are carried out with two pairs of feedback control parameters, which show different unfoldings of the maglev system and we verify the theoretical analysis.     

Resonant double Hopf bifurcation in a diffusive Ginzburg–Landau model with delayed feedback

We investigate the resonant double Hopf bifurcation in a diffusive complex Ginzburg–Landau model with delayed feedback and phase shift. The conditions for the existence of resonant double Hopf bifurcation are obtained by analyzing the roots’ distribution of the characteristic equation, and a general formula to determine the bifurcation point is given. For the cases of 1:2 and 1:3 resonance, we choose time delay, feedback strength and phase shift as bifurcation parameters and derive the normal forms which are proved to be the same as those in non-resonant cases. The impact of cubic terms on the unfolding types is discussed after obtaining the normal form till 3rd order. By fixing phase shift, we find that varying time delay and feedback strength simultaneously can induce the coexistence of two different periodic solutions, the existence of quasi-periodic solutions and strange attractors. Also, the effects on the existence of transient quasi-periodic solution exerted by the phase shift are illustrated.

     

Multiple scales and normal forms in a ring of delay coupled oscillators with application to chaotic Hindmarsh–Rose neurons

We study the appearance and stability of spatiotemporal periodic patterns like phase-locked oscillations, mirror-reflecting waves, standing waves, in-phase or antiphase oscillations, and coexistence of multiple patterns, in a ring of bidirectionally delay coupled oscillators. Hopf bifurcation, Hopf–Hopf bifurcation, and the equivariant Hopf bifurcation are studied in the viewpoint of normal forms obtained by using the method of multiple scales which is a kind of perturbation technique, thus a clear bifurcation scenario is depicted. We find time delay significantly affects the dynamics and induces rich spatiotemporal patterns. With the help of the unfolding system near Hopf–Hopf bifurcation, it is confirmed in some regions two kinds of stable oscillations may coexist. These phenomena are shown for the delay coupled limit cycle oscillators as well as for the delay coupled chaotic Hindmarsh–Rose neurons.     

Bifurcation dynamics of the modified physiological model of artificial pancreas with insulin secretion delay

In this paper, we modify the original physiological model of artificial pancreas by introducing the insulin secretion time delay. The non-resonant double Hopf bifurcation is analyzed by the Center Manifold Theorem and Normal Form Method. Numerical results supporting the theoretical analysis are presented in some typical parameter regions. It is shown that the critical value of technological delay and the area of death island of the non-resonant double Hopf bifurcation in the modified model are far less than those in the original model. This implies that when the secretion delay appears, the smaller technological delay can induce the double Hopf bifurcation. In addition, the region IV with complex coexisting bi-stability also decreases sharply. Furthermore, the rich dynamics such as various period, quasi-period and chaotic behaviors are found when some key parameters are changed. The obtained results can have important theoretical guidance for the diagnosis and treatment of diabetes patients.