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.

11.5. EXAMPLES OF BIOLOGICAL CELLULAR AUTOMATA MODELS 203R aR iFigure 11.10: CA model <strong>of</strong> Turing patterns. Left: Schematic illustration <strong>of</strong> <strong>the</strong> model’sstate-transition function. Right: A sample simulation result.This kind <strong>of</strong> spatial dynamics, driven by propagation <strong>of</strong> states between adjacent cellsthat are physically <strong>to</strong>uching each o<strong>the</strong>r, are called contact processes. This model <strong>and</strong> <strong>the</strong>following two are all examples <strong>of</strong> contact processes.A stylized CA model <strong>of</strong> this excitable media can be developed as follows. Assume atwo-dimensional space made <strong>of</strong> cells where each cell takes ei<strong>the</strong>r a normal (0; quiescent),excited (1), or refrac<strong>to</strong>ry (2, 3, . . . k) state. A normal cell becomes excited s<strong>to</strong>chasticallywith a probability determined by a function <strong>of</strong> <strong>the</strong> number <strong>of</strong> excited cells in its neighborhood.An excited cell becomes refrac<strong>to</strong>ry (2) immediately, while a refrac<strong>to</strong>ry cell remainsrefrac<strong>to</strong>ry for a while (but its state keeps counting up) <strong>and</strong> <strong>the</strong>n it comes back <strong>to</strong> normalafter it reaches k. Figure 11.11 shows <strong>the</strong> schematic illustration <strong>of</strong> this state-transitionfunction, as well as a sample simulation result you can get if you implement this modelsuccessfully. These kinds <strong>of</strong> uncoordinated traveling waves <strong>of</strong> excitation (including spirals)are actually experimentally observable on <strong>the</strong> surface <strong>of</strong> a human heart under cardiac arrhythmia.Exercise 11.7following:Implement a CA model <strong>of</strong> excitable media in Python. Then try <strong>the</strong>

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

Saved successfully!

Ooh no, something went wrong!