15.08.2015 Views

Introduction to the Modeling and Analysis of Complex Systems

introduction-to-the-modeling-and-analysis-of-complex-systems-sayama-pdf

introduction-to-the-modeling-and-analysis-of-complex-systems-sayama-pdf

SHOW MORE
SHOW LESS

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

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

122 CHAPTER 7. CONTINUOUS-TIME MODELS II: ANALYSISSimilar <strong>to</strong> discrete-time models, <strong>the</strong> dominant eigenvalues <strong>and</strong> eigenvec<strong>to</strong>rs tell us <strong>the</strong>asymp<strong>to</strong>tic behavior <strong>of</strong> continuous-time models, but with a little different stability criterion.Namely, if <strong>the</strong> real part <strong>of</strong> <strong>the</strong> dominant eigenvalue is greater than 0, <strong>the</strong>n <strong>the</strong> systemdiverges <strong>to</strong> infinity, i.e., <strong>the</strong> system is unstable. If it is less than 0, <strong>the</strong> system eventuallyshrinks <strong>to</strong> zero, i.e., <strong>the</strong> system is stable. If it is precisely 0, <strong>the</strong>n <strong>the</strong> dominant eigenvec<strong>to</strong>rcomponent <strong>of</strong> <strong>the</strong> system’s state is conserved with nei<strong>the</strong>r divergence nor convergence,<strong>and</strong> thus <strong>the</strong> system may converge <strong>to</strong> a non-zero equilibrium point. The same interpretationcan be applied <strong>to</strong> non-dominant eigenvalues as well.An eigenvalue tells us whe<strong>the</strong>r a particular component <strong>of</strong> a system’s state (givenby its corresponding eigenvec<strong>to</strong>r) grows or shrinks over time. For continuous-timemodels:• Re(λ) > 0 means that <strong>the</strong> component is growing.• Re(λ) < 0 means that <strong>the</strong> component is shrinking.• Re(λ) = 0 means that <strong>the</strong> component is conserved.For continuous-time models, <strong>the</strong> real part <strong>of</strong> <strong>the</strong> dominant eigenvalue λ d determines<strong>the</strong> stability <strong>of</strong> <strong>the</strong> whole system as follows:• Re(λ d ) > 0: The system is unstable, diverging <strong>to</strong> infinity.• Re(λ d ) < 0: The system is stable, converging <strong>to</strong> <strong>the</strong> origin.• Re(λ d ) = 0: The system is stable, but <strong>the</strong> dominant eigenvec<strong>to</strong>r component isconserved, <strong>and</strong> <strong>the</strong>refore <strong>the</strong> system may converge <strong>to</strong> a non-zero equilibriumpoint.Here is an example <strong>of</strong> a general two-dimensional linear dynamical system in continuoustime (a.k.a. <strong>the</strong> “love affairs” model proposed by Strogatz [29]):( )dx a bdt = x = Ax (7.52)c dThe eigenvalues <strong>of</strong> <strong>the</strong> coefficient matrix can be obtained by solving <strong>the</strong> following equationfor λ:( )a − λ bdet= 0 (7.53)c d − λ

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

Saved successfully!

Ooh no, something went wrong!