A note on Hopf bifurcations in a delayed diffusive Lotka-Volterra predator-prey system

The diffusive Lotka-Volterra predator-prey system with two delays is reconsidered here. The stability of the coexistence equilibrium and associated Hopf bifurcation are investigated by analyzing the characteristic equations, and our results complement earlier ones. We also obtain that in a special case, a Hopf bifurcation of spatial inhomogeneous periodic solutions occurs in the system.     

Categorization of neural excitability using threshold models

A classification of spiking neurons according to the transition from quiescence to periodic firing of action potentials is commonly used. Nonbursting neurons are classified into two types, type I and type II excitability. We use simple phenomenological spiking neuron models to derive a criterion for the determination of the neural excitability based on the afterpotential following a spike. The crucial characteristic is the existence for type II model of a positive overshoot, that is, a delayed after depolarization, during the recovery process of the membrane potential. Our prediction is numerically tested using well-known type I and type II models including the Connor, Walter, & McKown (1977) model and the Hodgkin-Huxley (1952) model.     

Interspike interval statistics in the stochastic Hodgkin-Huxley model: coexistence of gamma frequency bursts and highly irregular firing

When the classical Hodgkin-Huxley equations are simulated with Na- and K-channel noise and constant applied current, the distribution of interspike intervals is bimodal: one part is an exponential tail, as often assumed, while the other is a narrow gaussian peak centered at a short interspike interval value. The gaussian arises from bursts of spikes in the gamma-frequency range, the tail from the interburst intervals, giving overall an extraordinarily high coefficient of variation--up to 2.5 for 180,000 Na channels when I approximately 7 microA/cm(2). Since neurons with a bimodal ISI distribution are common, it may be a useful model for any neuron with class 2 firing. The underlying mechanism is due to a subcritical Hopf bifurcation, together with a switching region in phase-space where a fixed point is very close to a system limit cycle. This mechanism may be present in many different classes of neurons and may contribute to widely observed highly irregular neural spiking.     

转子-密封系统中气流激振力的非线性动力学特性分析

在高参数汽轮机组和航空发动机等旋转机械中，转子-密封中的气流激振力对转子非线性动力学特性的影响不容忽视．本研究中建立了转子-密封系统三维流场模型，应用计算流体动力学（CFD）软件对可压缩气流流场进行模拟计算，获得了密封流场特性．由流场计算结果进一步获得了Muszynska气流激振力模型中的相关经验系数，使得此模型更加适用于气流激振力的计算．在对转子一密封系统进行非线性动力学分析过程中应用幂级数展开形式建立了系统幂级数模型．利用平均法得到气流激振力的1：2亚谐共振分岔方程，进一步应用奇异性理论和Hopf分岔理论研究了系统1：2亚谐共振的转迁集和系统超临界Hopf分岔与亚临界Hopf分岔的存在条件．通过参数控制方法抑制了转子-密封系统出现亚临界分岔的出现，使得系统稳定性提高．本文的分析结果对工程设计和操作具有一定的指导作用和意义．     

On the phase reduction and response dynamics of neural oscillator populations

We undertake a probabilistic analysis of the response of repetitively firing neural populations to simple pulselike stimuli. Recalling and extending results from the literature, we compute phase response curves (PRCs) valid near bifurcations to periodic firing for Hindmarsh-Rose, Hodgkin-Huxley, FitzHugh-Nagumo, and Morris-Lecar models, encompassing the four generic (codimension one) bifurcations. Phase density equations are then used to analyze the role of the bifurcation, and the resulting PRC, in responses to stimuli. In particular, we explore the interplay among stimulus duration, baseline firing frequency, and population-level response patterns. We interpret the results in terms of the signal processing measure of gain and discuss further applications and experimentally testable predictions.     

分数阶计算机病毒模型的稳定性与分岔分析

随着有线和无线通信网络的普及,计算机病毒已经成为当代信息社会的一大威胁,单纯依靠杀毒软件已经无法彻底清除病毒,而通过对其在互联网上的传播机制的分析,以及对其模型的研究,可以找到有效的防范计算机病毒的对策。因此,基于非线性动力学与分数阶系统理论,建立了一类具有饱和发生率的分数阶时滞SIQR计算机病毒模型。计算出模型的平衡点,并通过分析相应的特征方程研究了时滞对平衡点稳定性的影响。选择时滞作为分岔参数,得到了发生Hopf分岔的时滞临界值。研究发现,系统的动力学行为依赖于分岔的临界值,同时给出了系统局部稳定和产生Hopf分岔的条件。在此基础上,研究了分数阶阶次的变化对分岔阈值的影响。最后,通过数值模拟验证了理论分析的正确性。     

A Hopf bifurcation theorem for singular differential–algebraic equations

We prove a Hopf bifurcation result for singular differential–algebraic equations (DAE) under the assumption that a trivial locus of equilibria is situated on the singularity as the bifurcation occurs. The structure that we need to obtain this result is that the linearisation of the DAE has a particular index-2 Kronecker normal form, which is said to be simple index-2. This is so-named because the nilpotent mapping used to define the Kronecker index of the pencil has the smallest possible non-trivial rank, namely one. This allows us to recast the equation in terms of a singular normal form to which a local centre-manifold reduction and, subsequently, the Hopf bifurcation theorem applies.     

Computing bifurcation points via characteristics gain loci

The authors show how to check the crossing on the imaginary axis by the eigenvalues of the linearized system of differential equations depending on a real parameter μ via feedback system theory. E. Hopf's theorem (1942) refers to a system of ordinary differential equations depending on the real parameter μ in which, when a single pair of complex conjugate eigenvalues of the linearized equations cross the imaginary axis under the parameter vibration, near this critical condition periodic orbits appear. The authors present simple formulas for both static (one eigenvalue zero) and dynamic or Hopf (a single pure imaginary pair) bifurcations, and show some singular conditions (degeneracies) by continuing the bifurcation curves in the steady-state manifold. The bifurcation curves and singular sets of an interesting chemical reactor which possesses multiplicity in the equilibrium solutions and in the Hopf bifurcation points are described     

Numerical analysis of bifurcations by continuation methods

A numerical calculation technique for computing bifurcation values of interconnected dynamical systems is presented. The technique is based on continuation methods in which the bifurcation value of interconnected dynamical systems can be calculated from the bifurcation value of a subsystem using a set of coupled differential equations. As an example, the value of the Hopf bifurcation of an interconnected dynamical system is calculated.     

机床无刷直流电机系统的分岔分析与控制

文章主要研究了机床无刷直流电机系统的Hopf分岔控制问题.首先,对系统进行分岔分析,通过计算极限环曲率系数判定系统的Hopf分岔类型;然后设计Washout滤波器对系统进行分岔控制,根据Hopf分岔理论给出使原系统Hopf分岔位置发生改变的参数条件,利用Normal Form方法计算出受控系统的Hopf分岔正规型,根据正规型的实部大小判定Hopf分岔类型,给出使原系统Hopf分岔类型发生改变的参数条件;并借助MATLAB软件对理论结果进行数值仿真,理论结果和数值仿真表明:控制器中的线性增益能使系统在所期望的参数值处发生Hopf分岔,甚至消除Hopf分岔,控制器中的非线性增益能改变原系统的Hopf分岔类型及极限环幅值的大小.研究结果对无刷直流电动机系统的工程实际具有一定的指导意义.     

规范形Qi系统的Hopf分岔分析及控制

文章研究了一个Qi系统的Hopf分岔控制问题.根据计算的极限环曲率系数,判定原系统的Hopf分岔类型,并采用washout滤波器控制该系统的分岔行为.首先讨论了控制器的线性增益对Hopf分岔点位置的影响,然后引入规范形计算方法,求出受控系统的Hopf分岔规范形.分析了规范形中系数对控制参数的选择原则所产生之影响,以及对Hopf分岔类型及极限环幅值的影响.理论和仿真结果表明,控制器的线性增益能使原系统的Hopf分岔点延迟或消失,而非线性增益能则改变极限环的稳定性和极限环幅值的大小.最后把washout滤波器和线性控制器的控制效果作了对比,发现washout滤波器比之线性控制器具有一定的优势.     

传输时延环境下信息物理融合系统中恶意病毒传播的稳定性与分岔分析

本文考虑到恶意病毒在信息物理融合系统中的传播具有时延性,基于非线性动力学理论建立了一类更具一般性的含有时滞的恶意病毒传播模型.通过选取时滞作为分岔参数,并讨论相关的特征方程,研究了时滞对系统局部稳定性和Hopf分岔的影响.研究发现,系统的动力学行为依赖于分岔的临界值.此外,给出了保证系统稳定性和产生Hopf分岔的条件....     

Toric ideals and graph theory to analyze Hopf bifurcations in mass action systems

A family of polynomial differential systems describing the behavior of a chemical reaction network with generalized mass action kinetics is investigated. The coefficients and monomials are given by graphs. The aim of this investigation is to clarify the algebraic-discrete aspects of a Hopf bifurcation in these special differential equations. We apply concepts from toric geometry and convex geometry. As usual in stoichiometric network analysis we consider the solution set as a convex polyhedral cone and we intersect it with the deformed toric variety of the monomials. Using Gröbner bases the polynomial entries of the Jacobian are expressed in different coordinate systems. Then the Hurwitz criterion is applied in order to determine parameter regions where a Hopf bifurcation occurs. Examples from chemistry illustrate the theoretical results.     

Control bifurcations

A parametrized nonlinear differential equation can have multiple equilibria as the parameter is varied. A local bifurcation of a parametrized differential equation occurs at an equilibrium where there is a change in the topological character of the nearby solution curves. This typically happens because some eigenvalues of the parametrized linear approximating differential equation cross the imaginary axis and there is a change in stability of the equilibrium. The topological nature of the solutions is unchanged by smooth changes of state coordinates so these may be used to bring the differential equation into Poincare/spl acute/ normal form. From this normal form, the type of the bifurcation can be determined. For differential equations depending on a single parameter, the typical ways that the system can bifurcate are fully understood, e.g., the fold (or saddle node), the transcritical and the Hopf bifurcation. A nonlinear control system has multiple equilibria typically parametrized by the set value of the control. A control bifurcation of a nonlinear system typically occurs when its linear approximation loses stabilizability. The ways in which this can happen are understood through the appropriate normal forms. We present the quadratic and cubic normal forms of a scalar input nonlinear control system around an equilibrium point. These are the normal forms under quadratic and cubic change of state coordinates and invertible state feedback. The system need not be linearly controllable. We study some important control bifurcations, the analogues of the classical fold, transcritical and Hopf bifurcations.     

The Hopf bifurcation theorem for quasilinear differential-algebraic equations

We prove a generalization of the Hopf bifurcation theorem for quasilinear differential equations (DAEs), i.e. equations of the form A(μ, χ)χ = G(μ, χ) where the matrix A(χ, μ) has constant but not full rank and hence the system cannot be made into an explicit ODE. The paper includes an appendix by J. Ernsthausen addressing the numerical calculation of the Hopf points in the DAE setting.     

Bifurcation analysis in a silicon neuron

In this paper, we describe an analysis of the nonlinear dynamical phenomenon associated with a silicon neuron. Our silicon neuron in Very Large Scale Integration (VLSI) integrates Hodgkin?CHuxley (HH) model formalism, including the membrane voltage dependency of temporal dynamics. Analysis of the bifurcation conditions allow us to identify different regimes in the parameter space that are desirable for biasing our silicon neuron. This approach of studying bifurcations is useful because it is believed that computational properties of neurons are based on the bifurcations exhibited by these dynamical systems in response to some changing stimulus. We describe numerical simulations of the Hopf bifurcation which is characteristic of class 2 excitability in the HH model. We also show experimental measurements of an observed phenomenon in biological neurons and termed excitation block, firing rate and effect of current impulses. Hence, by showing that this silicon neuron has similar bifurcations to a certain class of biological neurons, we can claim that the silicon neuron can also perform similar computations.     

具有时滞和饱和接触率的SIRS模型的Hopf分支

应用微分方程分支理论,研究了具有时滞和饱和接触率的SIRS模型,以时滞[τ]为分支参数,运用Hopf分支理论,得到当时滞[τ]充分小时正平衡点是局部渐近稳定的,当[τ]经过一系列临界值时模型出现Hopf分支。用Matlab软件进行数值仿真验证了结论的正确性。     

A modeling study predicts the presence of voltage gated Ca2+ channels on myelinated central axons

The objective of this current study was to investigate whether voltage gated Ca(2+) channels are present on axons of the adult rat optic nerve (RON). Simulations of axonal excitability using a Hodgkin-Huxley based one-compartment model incorporating I(Na), I(K) and leak currents were used to predict conditions under which the potential contribution of a Ca(2+) current to an evoked action potential could be measured. Under control conditions the inclusion of a high threshold Ca(2+) current (I(Ca)) in the model had a negligible effect on the action potential. Reducing I(K), by decreasing the value of g(K), elongated the repolarizing phase of the action potential, increasing its duration. Subsequent incorporation of I(Ca) in the model revealed a significant I(Ca) contribution to the repolarizing phase of the action potential. The simulation thus suggests that Ca(2+) channels may be present on RON axons, but that pharmacological intervention is required to unmask their presence. Experiments based on the simulations revealed that there was no significant contribution of I(Ca) to the control action potential. However, as predicted by the simulation, reducing the repolarizing effect of I(K) by adding the K(+) channel blocker 4-AP revealed a Ca(2+) component on the repolarizing phase of the action potential that was blocked by the Ca(2+) channel inhibitor nifedipine.     

A numerical algorithm coupling a bifurcating indicator and a direct method for the computation of Hopf bifurcation points in fluid mechanics

This paper deals with the computation of Hopf bifurcation points in fluid mechanics. This computation is done by coupling a bifurcation indicator proposed recently (Cadou et al., 2006) [1] and a direct method (Jackson, 1987; Jepson, 1981) [2] and [3] which consists in solving an augmented system whose solutions are Hopf bifurcation points. The bifurcation indicator gives initial critical values (Reynolds number, Strouhal frequency) for the direct method iterations. Some classical numerical examples from fluid mechanics, in two dimensions, are studied to demonstrate the efficiency and the reliability of such an algorithm.