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.

152 CHAPTER 8. BIFURCATIONSExercise 8.5 Logistic map Period-doubling bifurcations <strong>and</strong> chaos are not justfor abstract, contrived ma<strong>the</strong>matical equations, but <strong>the</strong>y can occur in various models<strong>of</strong> real-world biological, ecological, social, <strong>and</strong> engineering phenomena. Thesimplest possible example would be <strong>the</strong> logistic map we introduced in Section 5.5:x t = rx t−1 (1 − x t−1 ) (8.42)This is a ma<strong>the</strong>matical model <strong>of</strong> population dynamics, where x t represents <strong>the</strong>population <strong>of</strong> a species that reproduce in discrete (non-overlapping) generations.This model was used by British/Australian ma<strong>the</strong>matical ecologist Robert May inhis influential 1976 Nature paper [32] <strong>to</strong> illustrate how a very simple ma<strong>the</strong>maticalmodel could produce as<strong>to</strong>nishingly complex behaviors.• Conduct a bifurcation analysis <strong>of</strong> this model <strong>to</strong> find <strong>the</strong> critical thresholds <strong>of</strong> rat which bifurcations occur.• Study <strong>the</strong> stability <strong>of</strong> each equilibrium point in each parameter range <strong>and</strong> summarize<strong>the</strong> results in a table.• Simulate <strong>the</strong> model with several selected values <strong>of</strong> r <strong>to</strong> confirm <strong>the</strong> results <strong>of</strong>analysis.• Draw a bifurcation diagram <strong>of</strong> this model for 0 < r < 4.Exercise 8.6 Stability analysis <strong>of</strong> periodic trajec<strong>to</strong>ries The stability <strong>of</strong> aperiod-2 trajec<strong>to</strong>ry <strong>of</strong> a discrete-time model x t = F (x t−1 ) can be studied by <strong>the</strong>stability analysis <strong>of</strong> ano<strong>the</strong>r model made <strong>of</strong> a composite function <strong>of</strong> F :y τ = G(y τ−1 ) = F (F (y τ−1 )) (8.43)This is because <strong>the</strong> period-2 trajec<strong>to</strong>ry in F corresponds <strong>to</strong> one <strong>of</strong> <strong>the</strong> equilibriumpoints <strong>of</strong> G(◦) = F (F (◦)). If such an equilibrium point <strong>of</strong> G is being destabilizedso that dG/dx ≈ −1, it means that <strong>the</strong> period-2 trajec<strong>to</strong>ry in question is losing <strong>the</strong>stability, <strong>and</strong> thus ano<strong>the</strong>r period-doubling bifurcation in<strong>to</strong> a period-4 trajec<strong>to</strong>ry isabout <strong>to</strong> occur. Using this technique, analytically obtain <strong>the</strong> critical threshold <strong>of</strong> rin Eq. (8.37) at which <strong>the</strong> second period-doubling bifurcation occurs (from period 2<strong>to</strong> period 4). Then, draw cobweb plots <strong>of</strong> y τ = G(y τ−1 ) for several values <strong>of</strong> r near<strong>the</strong> critical threshold <strong>to</strong> see what is happening <strong>the</strong>re.

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

Saved successfully!

Ooh no, something went wrong!