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.

14.3. LINEAR STABILITY ANALYSIS OF CONTINUOUS FIELD MODELS 281dynamical system. What we have done is <strong>to</strong> constrain <strong>the</strong> shape <strong>of</strong> <strong>the</strong> small perturbations(i.e., deviations from <strong>the</strong> homogeneous equilibrium state) <strong>to</strong> a certain eigenfunction<strong>to</strong> eliminate spatial effects, <strong>and</strong> <strong>the</strong>n ignore any higher-order terms <strong>to</strong> linearize <strong>the</strong> dynamics.Each point in this new (∆a, ∆c) phase space still represents a certain spatialconfiguration <strong>of</strong> <strong>the</strong> original model, as illustrated in Fig. 14.3. Studying <strong>the</strong> stability <strong>of</strong> <strong>the</strong>origin (∆a, ∆c) = (0, 0) tells us if <strong>the</strong> Keller-Segel model can remain homogeneous or if itundergoes a spontaneous pattern formation.aDcaa eqcxa eqcxc eqxc eqxDaaaa eqcxa eqcxc eqxc eqxhomogeneousequilibrium stateFigure 14.3: Visual illustration <strong>of</strong> how <strong>to</strong> interpret each point in <strong>the</strong> (∆a, ∆c) phasespace. Each call-out shows <strong>the</strong> actual state <strong>of</strong> <strong>the</strong> system that corresponds <strong>to</strong> eachpoint in <strong>the</strong> phase space.The only thing we need <strong>to</strong> do is <strong>to</strong> calculate <strong>the</strong> eigenvalues <strong>of</strong> <strong>the</strong> matrix in Eq. (14.65)<strong>and</strong> check <strong>the</strong> signs <strong>of</strong> <strong>the</strong>ir real parts. This may be cumbersome, but let’s get it done.Here is <strong>the</strong> calculation process, where λ is <strong>the</strong> eigenvalue <strong>of</strong> <strong>the</strong> matrix:

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

Saved successfully!

Ooh no, something went wrong!