12.07.2015 Views

Dynamical Systems in Neuroscience:

Dynamical Systems in Neuroscience:

Dynamical Systems in Neuroscience:

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

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

96 Two-Dimensional <strong>Systems</strong>def<strong>in</strong>es a constant vertical vector field depicted <strong>in</strong> Fig. 4.3b. The systemẋ = −x ,ẏ = −ydef<strong>in</strong>es a vector field that po<strong>in</strong>ts to the orig<strong>in</strong> (0, 0), as <strong>in</strong> Fig. 4.3c, and the systemẋ = −y , (4.3)ẏ = −x (4.4)def<strong>in</strong>es a saddle vector field, as <strong>in</strong> Fig. 4.3d. Vector fields provide geometrical <strong>in</strong>formationabout the jo<strong>in</strong>t evolution of state variables. For example, the vector field <strong>in</strong>Fig. 4.3d is directed rightward <strong>in</strong> the lower half-plane and leftward <strong>in</strong> the upper halfplane.Therefore the variable x(t) <strong>in</strong>creases when y < 0 and decreases otherwise, whichobviously follows from the equation (4.3). Quite often however geometrical analysis ofvector fields can provide <strong>in</strong>formation about the behavior of the system that may notbe obvious from the form of the functions f and g.4.1.1 Nullcl<strong>in</strong>esThe vector field <strong>in</strong> Fig. 4.3d is directed rightward (x <strong>in</strong>creases) or leftward (x decreases)<strong>in</strong> different regions of the phase plane. The set of po<strong>in</strong>ts where the vectorfield changes its horizontal direction is called x-nullcl<strong>in</strong>e, and it is def<strong>in</strong>ed by the equationf(x, y) = 0. Indeed, at any such po<strong>in</strong>t x neither <strong>in</strong>creases nor decreases becauseẋ = 0. The x-nullcl<strong>in</strong>e partitions the phase plane <strong>in</strong>to two regions where x moves <strong>in</strong>opposite directions. Similarly, the y-nullcl<strong>in</strong>e is def<strong>in</strong>ed by the equation g(x, y) = 0,and it denotes the set of po<strong>in</strong>ts where the vector field changes its vertical direction.This nullcl<strong>in</strong>e partitions the phase plane <strong>in</strong>to two regions where y either <strong>in</strong>creases ordecreases. The x- and y-nullcl<strong>in</strong>es partition the phase plane <strong>in</strong>to 4 different regions:(a) x and y <strong>in</strong>crease, (b) x decreases, y <strong>in</strong>creases, (c) x and y decrease, and (d) x<strong>in</strong>creases, y decreases, as we illustrate <strong>in</strong> Fig. 4.4.Each po<strong>in</strong>t of <strong>in</strong>tersection of the nullcl<strong>in</strong>es is an equilibrium po<strong>in</strong>t, s<strong>in</strong>ce f(x, y) =g(x, y) = 0 and hence ẋ = ẏ = 0. Conversely, every equilibrium of a two-dimensionalsystem is the po<strong>in</strong>t of <strong>in</strong>tersection of its nullcl<strong>in</strong>es. Because nullcl<strong>in</strong>es are so important,we consider two examples <strong>in</strong> detail below (the reader is urged to solve Ex. 1 at the endof this chapter).Let us determ<strong>in</strong>e nullcl<strong>in</strong>es of the system (4.3, 4.4) with the vector field shown <strong>in</strong>Fig. 4.3d. From (4.3) it follows that x-nullcl<strong>in</strong>e is the horizontal l<strong>in</strong>e y = 0, and from(4.4) it follows that y-nullcl<strong>in</strong>e is the vertical l<strong>in</strong>e x = 0. These nullcl<strong>in</strong>es (dashed l<strong>in</strong>es<strong>in</strong> Fig. 4.3d) partition the phase plane <strong>in</strong>to 4 quadrants, <strong>in</strong> each of which the vectorfield has a different direction. The <strong>in</strong>tersection of the nullcl<strong>in</strong>es is the equilibrium(0, 0). Later <strong>in</strong> this chapter we will study how to determ<strong>in</strong>e stability of equilibria <strong>in</strong>two-dimensional systems, though <strong>in</strong> this particular case one can easily guess that theequilibrium is not stable.

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

Saved successfully!

Ooh no, something went wrong!