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.

144 CHAPTER 8. BIFURCATIONSWith a = 0.7, b = 0.8, <strong>and</strong> c = 3, do <strong>the</strong> following:• Numerically obtain <strong>the</strong> equilibrium point <strong>of</strong> this model for several values <strong>of</strong> z,ranging between -2 <strong>and</strong> 0. There is only one real equilibrium point in thissystem.• Apply <strong>the</strong> result obtained above <strong>to</strong> <strong>the</strong> Jacobian matrix <strong>of</strong> <strong>the</strong> model, <strong>and</strong>numerically evaluate <strong>the</strong> stability <strong>of</strong> that equilibrium point for each value <strong>of</strong> z.• Estimate <strong>the</strong> critical thresholds <strong>of</strong> z at which a Hopf bifurcation occurs. Thereare two such critical thresholds.• Draw a series <strong>of</strong> its phase spaces with values <strong>of</strong> z varied from 0 <strong>to</strong> -2 <strong>to</strong>confirm your analytical prediction.8.4 Bifurcations in Discrete-Time ModelsThe bifurcations discussed above (saddle-node, transcritical, pitchfork, Hopf) are alsopossible in discrete-time dynamical systems with one variable:x t = F (x t−1 ) (8.35)The Jacobian matrix <strong>of</strong> this system is, again, a 1×1 matrix whose eigenvalue is its contentitself, which is given by dF/dx. Since this is a discrete-time model, <strong>the</strong> critical condition atwhich a bifurcation occurs is given bydF∣ dx ∣ = 1. (8.36)x=xeqLet’s work on <strong>the</strong> following example:x t = x t−1 + r − x 2 t−1 (8.37)This is a discrete-time analog <strong>of</strong> Eq. (8.3). Therefore, it has <strong>the</strong> same set <strong>of</strong> equilibriumpoints:x eq = ± √ r (8.38)Next, we calculate dF/dx as follows:dFdx = (r + x − x2 ) ′ = 1 − 2x (8.39)dF∣ dx ∣ = |1 ± 2 √ r| (8.40)√ x=± r

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

Saved successfully!

Ooh no, something went wrong!