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.

5.6. ASYMPTOTIC BEHAVIOR OF DISCRETE-TIME LINEAR DYNAMICAL... 85The first array shows <strong>the</strong> list <strong>of</strong> eigenvalues (it surely includes <strong>the</strong> golden ratio), while <strong>the</strong>second one shows <strong>the</strong> eigenvec<strong>to</strong>r matrix (i.e., a square matrix whose column vec<strong>to</strong>rsare eigenvec<strong>to</strong>rs <strong>of</strong> <strong>the</strong> original matrix). The eigenvec<strong>to</strong>rs are listed in <strong>the</strong> same order aseigenvalues. Now we need <strong>to</strong> interpret this result. The eigenvalues are 1.61803399 <strong>and</strong>-0.61803399. Which one is dominant?The answer is <strong>the</strong> first one, because its absolute value is greater than <strong>the</strong> secondone’s. This means that, asymp<strong>to</strong>tically, <strong>the</strong> system’s behavior looks like this:x t ≈ 1.61803399 x t−1 (5.52)Namely, <strong>the</strong> dominant eigenvalue tells us <strong>the</strong> asymp<strong>to</strong>tic ratio <strong>of</strong> magnitudes <strong>of</strong> <strong>the</strong>state vec<strong>to</strong>rs between two consecutive time points (in this case, it approaches <strong>the</strong> goldenratio). If <strong>the</strong> absolute value <strong>of</strong> <strong>the</strong> dominant eigenvalue is greater than 1, <strong>the</strong>n <strong>the</strong> systemwill diverge <strong>to</strong> infinity, i.e., <strong>the</strong> system is unstable. If less than 1, <strong>the</strong> system will eventuallyshrink <strong>to</strong> zero, i.e., <strong>the</strong> system is stable. If it is precisely 1, <strong>the</strong>n <strong>the</strong> dominant eigenvec<strong>to</strong>rcomponent <strong>of</strong> <strong>the</strong> system’s state will be 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 sameinterpretation can 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 (given byits corresponding eigenvec<strong>to</strong>r) grows or shrinks over time. For discrete-time models:• |λ| > 1 means that <strong>the</strong> component is growing.• |λ| < 1 means that <strong>the</strong> component is shrinking.• |λ| = 1 means that <strong>the</strong> component is conserved.For discrete-time models, <strong>the</strong> absolute value <strong>of</strong> <strong>the</strong> dominant eigenvalue λ d determines<strong>the</strong> stability <strong>of</strong> <strong>the</strong> whole system as follows:• |λ d | > 1: The system is unstable, diverging <strong>to</strong> infinity.• |λ d | < 1: The system is stable, converging <strong>to</strong> <strong>the</strong> origin.• |λ d | = 1: 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.

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

Saved successfully!

Ooh no, something went wrong!