Dynamical Systems in Neuroscience:

Dynamical Systems in Neuroscience: Dynamical Systems in Neuroscience:

ib.cnea.gov.ar
from ib.cnea.gov.ar More from this publisher
12.07.2015 Views

106 Two-Dimensional SystemsFor example, is the equilibrium in Fig. 4.4 stable? What about the equilibrium inFig. 4.10? The vector fields in the neighborhoods of the two equilibria exhibit subtledifferences that are difficult to spot without the help of analytical tools, which wediscuss next.4.2.2 Local linear analysisBelow we remind the reader some basic concepts of linear algebra, assuming that thereader has some familiarity with matrices, eigenvectors and eigenvalues. Consider atwo-dimensional dynamical systemẋ = f(x, y) (4.5)ẏ = g(x, y) (4.6)having an equilibrium point (x 0 , y 0 ). The nonlinear functions f and g can be linearizednear the equilibrium; i.e., written in the formf(x, y) = a(x − x 0 ) + b(y − y 0 ) + higher-order terms,g(x, y) = c(x − x 0 ) + d(y − y 0 ) + higher-order terms,where higher-order terms include (x − x 0 ) 2 , (x − x 0 )(y − y 0 ), (x − x 0 ) 3 , etc., anda = ∂f∂x (x 0, y 0 ),c = ∂g∂x (x 0, y 0 ),b = ∂f∂y (x 0, y 0 ),d = ∂g∂y (x 0, y 0 ),are the partial derivatives of f and g with respect of the state variables x and yevaluated at the equilibrium (x 0 , y 0 ) (first, evaluate the derivatives, then substitutex = x 0 and y = y 0 ). Many questions regarding the stability of the equilibrium can beanswered by considering the corresponding linear system˙u = au + bw , (4.7)ẇ = cu + dw , (4.8)where u = x − x 0 and w = y − y 0 are the deviations from the equilibrium, and thehigher-order terms, u 2 , uw, w 3 , etc., are neglected. We can write this system in thevector form( ) ( ) ( )˙u a b u=.ẇ c d wThe linearization matrix( ) a bL =c d

Two-Dimensional Systems 107is called the Jacobian matrix of the system (4.5, 4.6) at the equilibrium (x 0 , y 0 ). Forexample, the Jacobian matrix of the system (4.3, 4.4) at the origin is( ) 0 −1. (4.9)−1 0It is important to remember that Jacobian matrices are defined for equilibria, andthat a nonlinear system can have many equilibria and hence many different Jacobianmatrices.4.2.3 Eigenvalues and eigenvectorsA non-zero vector v ∈ R 2 is said to be an eigenvector of the matrix L correspondingto the eigenvalue λ ifLv = λv (matrix notation) .For example, the matrix (4.9) has two eigenvectors( ) 1v 1 =and v12 =( 1−1corresponding to the eigenvalues λ 1 = −1 and λ 2 = 1, respectively. Any textbook onlinear algebra explains how to find eigenvectors and eigenvalues of an arbitrary matrix.It is important for the reader to get comfortable with these notions, since they are usedextensively in the rest of the book.Eigenvalues play important role in analysis of stability of equilibria. To find theeigenvalues of a 2×2-matrix L, one solves the characteristic equation( )a − λ bdet= 0 .c d − λThis equation can be written in the polynomial form (a − λ)(d − λ) − bc = 0 orwhereλ 2 − τλ + ∆ = 0 ,τ = tr L = a + d and ∆ = det L = ad − bcare the trace and the determinant of the matrix L, respectively.polynomial has two solutions of the form)Such a quadraticλ 1 = τ + √ τ 2 − 4∆2andλ 2 = τ − √ τ 2 − 4∆2(4.10)and they are either real (when τ 2 −4∆ ≥ 0) or complex-conjugate (when τ 2 −4∆ < 0).What can you say about the case τ 2 = 4∆?

Two-Dimensional <strong>Systems</strong> 107is called the Jacobian matrix of the system (4.5, 4.6) at the equilibrium (x 0 , y 0 ). Forexample, the Jacobian matrix of the system (4.3, 4.4) at the orig<strong>in</strong> is( ) 0 −1. (4.9)−1 0It is important to remember that Jacobian matrices are def<strong>in</strong>ed for equilibria, andthat a nonl<strong>in</strong>ear system can have many equilibria and hence many different Jacobianmatrices.4.2.3 Eigenvalues and eigenvectorsA non-zero vector v ∈ R 2 is said to be an eigenvector of the matrix L correspond<strong>in</strong>gto the eigenvalue λ ifLv = λv (matrix notation) .For example, the matrix (4.9) has two eigenvectors( ) 1v 1 =and v12 =( 1−1correspond<strong>in</strong>g to the eigenvalues λ 1 = −1 and λ 2 = 1, respectively. Any textbook onl<strong>in</strong>ear algebra expla<strong>in</strong>s how to f<strong>in</strong>d eigenvectors and eigenvalues of an arbitrary matrix.It is important for the reader to get comfortable with these notions, s<strong>in</strong>ce they are usedextensively <strong>in</strong> the rest of the book.Eigenvalues play important role <strong>in</strong> analysis of stability of equilibria. To f<strong>in</strong>d theeigenvalues of a 2×2-matrix L, one solves the characteristic equation( )a − λ bdet= 0 .c d − λThis equation can be written <strong>in</strong> the polynomial form (a − λ)(d − λ) − bc = 0 orwhereλ 2 − τλ + ∆ = 0 ,τ = tr L = a + d and ∆ = det L = ad − bcare the trace and the determ<strong>in</strong>ant of the matrix L, respectively.polynomial has two solutions of the form)Such a quadraticλ 1 = τ + √ τ 2 − 4∆2andλ 2 = τ − √ τ 2 − 4∆2(4.10)and they are either real (when τ 2 −4∆ ≥ 0) or complex-conjugate (when τ 2 −4∆ < 0).What can you say about the case τ 2 = 4∆?

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

Saved successfully!

Ooh no, something went wrong!