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.

120 CHAPTER 7. CONTINUOUS-TIME MODELS II: ANALYSISExercise 7.9 Simplify <strong>the</strong> following two-dimensional differential equation modelby variable rescaling:dx= ax(1 − x) − bxydt(7.38)dy= cy(1 − y) − dxydt(7.39)7.4 Asymp<strong>to</strong>tic Behavior <strong>of</strong> Continuous-Time Linear Dynamical<strong>Systems</strong>A general formula for continuous-time linear dynamical systems is given bydx= Ax, (7.40)dtwhere x is <strong>the</strong> state vec<strong>to</strong>r <strong>of</strong> <strong>the</strong> system <strong>and</strong> A is <strong>the</strong> coefficient matrix. As discussedbefore, you could add a constant vec<strong>to</strong>r a <strong>to</strong> <strong>the</strong> right h<strong>and</strong> side, but it can always be convertedin<strong>to</strong> a constant-free form by increasing <strong>the</strong> dimensions <strong>of</strong> <strong>the</strong> system, as follows:⎛ ⎞y = ⎝ x 1⎠ (7.41)⎛dydt = ⎝ A a0 0⎞ ⎛⎠ ⎝ x 1⎞⎠ = By (7.42)Note that <strong>the</strong> last-row, last-column element <strong>of</strong> <strong>the</strong> exp<strong>and</strong>ed coefficient matrix is now 0,not 1, because <strong>of</strong> Eq. (6.25). This result guarantees that <strong>the</strong> constant-free form given inEq. (7.40) is general enough <strong>to</strong> represent various dynamics <strong>of</strong> linear dynamical systems.Now, what is <strong>the</strong> asymp<strong>to</strong>tic behavior <strong>of</strong> Eq. (7.40)? This may not look so intuitive, butit turns out that <strong>the</strong>re is a closed-form solution available for this case as well. Here is <strong>the</strong>solution, which is generally applicable for any square matrix A:x(t) = e At x(0) (7.43)Here, e X is a matrix exponential for a square matrix X, which is defined as∞∑e X X k=k! , (7.44)k=0

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

Saved successfully!

Ooh no, something went wrong!