12.07.2015 Views

Dynamical Systems in Neuroscience:

Dynamical Systems in Neuroscience:

Dynamical Systems in Neuroscience:

SHOW MORE
SHOW LESS
  • No tags were found...

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

72 One-Dimensional <strong>Systems</strong>λ(V-Veq)VeqF(V)VFigure 3.21: Hartman-Grobman theorem:The non-l<strong>in</strong>ear system ˙V = F (V )is topologically equivalent to the l<strong>in</strong>earone ˙V = λ(V −V eq ) <strong>in</strong> the local (shaded)neighborhood of the hyperbolic equilibriumV eq .our analysis to a small neighborhood of the l<strong>in</strong>e R, e.g., a neighborhood of the rest<strong>in</strong>gstate or of the threshold, and study behavior locally <strong>in</strong> this neighborhood.An important tool <strong>in</strong> the local analysis of dynamical systems is the Hartman-Grobman theorem, which says that a non-l<strong>in</strong>ear one-dimensional system˙V = F (V )sufficiently near an equilibrium V = V eq is locally topologically equivalent to the l<strong>in</strong>earone˙V = λ(V − V eq ) (3.8)provided that the eigenvalueλ = F ′ (V eq )at the equilibrium is non-zero, i.e., the slope of F (V ) is non-zero. Such an equilibriumis called hyperbolic. Thus, nonl<strong>in</strong>ear systems near hyperbolic equilibria behave as ifthere were l<strong>in</strong>ear, as <strong>in</strong> Fig. 3.21.It is easy to f<strong>in</strong>d the exact solution of the l<strong>in</strong>earized system (3.8) with an <strong>in</strong>itialcondition V (0) = V 0 . It is V (t) = V eq + e λt (V 0 − V eq ) (check by differentiat<strong>in</strong>g).If the eigenvalue λ < 0, then e λt → 0 and V (t) → V eq as t → ∞, so that theequilibrium is stable. Conversely, if λ > 0, then e λt → ∞ mean<strong>in</strong>g that the <strong>in</strong>itialdisplacement, V 0 − V eq , grows with the time, and the equilibrium is unstable. Thus,the l<strong>in</strong>earization predicts qualitative dynamics at the equilibrium and quantitative rateof convergence/divergence to/from the equilibrium.If the eigenvalue λ = 0, then the equilibrium is non-hyperbolic, and analysis ofthe l<strong>in</strong>earized system ˙V = 0 cannot describe the behavior of the nonl<strong>in</strong>ear system.Typically, non-hyperbolic equilibria arise when the system undergoes a bifurcation,i.e., a qualitative change of behavior, which we consider next. To study stability, weneed to consider higher-order terms of the Taylor series of F (V ) at V eq .3.3.3 BifurcationsThe f<strong>in</strong>al and the most advanced step <strong>in</strong> the qualitative analysis of any dynamicalsystem is the bifurcation analysis. In general, a system is said to undergo a bifurcationwhen its phase portrait changes qualitatively. For example, the energy landscape <strong>in</strong>Fig. 3.22 changes so that the system is no longer bistable. The precise mathematicaldef<strong>in</strong>ition of a bifurcation will be given later.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!