A Model of Evolutionary Dynamics with Quasiperiodic Forcing
|
|
- Marilynn Griffin
- 5 years ago
- Views:
Transcription
1 paper presented at Society for Experimental Mechanics (SEM) IMAC XXXIII Conference on Structural Dynamics February 2-5, 205, Orlando FL A Model of Evolutionary Dynamics with Quasiperiodic Forcing Elizabeth Wesson Center for Applied Mathematics Cornell University Ithaca, NY 4853, USA Richard Rand Professor of Mathematics and Mech. & Aero. Eng. Cornell University Ithaca, NY 4853, USA October 23, 204 Abstract Evolutionary dynamics combines game theory and nonlinear dynamics to model competition in biological and social situations. The replicator equation is a standard paradigm in evolutionary dynamics. The growth rate of each strategy is its excess fitness: the deviation of its fitness from the average. The game-theoretic aspect of the model lies in the choice of fitness function, which is determined by a payoff matrix. Previous work by Ruelas and Rand investigated the Rock-Paper-Scissors replicator dynamics problem with periodic forcing of the payoff coefficients. This work extends the previous to consider the case of quasiperiodic forcing. This model may find applications in biological or social systems where competition is affected by cyclical processes on different scales, such as days/years or weeks/years. We study the quasiperiodically forced Rock-Paper-Scissors problem using numerical simulation, and Floquet theory and harmonic balance. We investigate the linear stability of the interior equilibrium point; we find that the region of stability in frequency space has fractal boundary. Keywords Replicator equation, quasiperiodic forcing, Floquet theory, harmonic balance, numerical integration Introduction The field of evolutionary dynamics combines game theory and nonlinear dynamics to model population shifts due to competition in biological and social situations. One standard paradigm [9, 3 uses the replicator equation, ẋ i = x i (f i (x) φ), i =,..., n () where x i is the frequency, or relative abundance, of strategy i; the unit vector x is the vector of frequencies; f i (x) is the fitness of strategy i; and φ is the average fitness, defined by φ = i x i f i (x). (2) The replicator equation can be derived [0 from an exponential model of population growth, ξ i = ξ i f i, i =,..., n (3) where ξ i is the population of strategy i, assuming that f i depends only on the frequencies: f i = f i (x). The derivation consists of a simple change of variables: x i ξ i /p where p = i ξ i is the total population. The game-theoretic component of the replicator model lies in the choice of fitness functions. Define the payoff matrix A = (a ij ) where a ij is the expected reward for a strategy i individual vs. a strategy j individual. We assume
2 the population is well-mixed, so that any individual competes against each strategy at a rate proportional to that strategy s frequency in the population. Then the fitness f i is the total expected payoff for strategy i vs. all strategies: f i (x) = (Ax) i = j a ij x j. (4) In this work, we generalize the replicator model to systems in which the payoff coefficients are quasiperiodic functions of time. Previous work by Ruelas and Rand [7, 8 investigated the Rock-Paper-Scissors replicator dynamics problem with periodic forcing of the payoff coefficients. We also consider a forced Rock-Paper-Scissors system. The quasiperiodically forced replicator model may find applications in biological or social systems where competition is affected by cyclical processes on different scales, such as days/years or weeks/years. 2 The model 2. Rock-Paper-Scissors games with quasiperiodic forcing Rock-Paper-Scissors (RPS) games are a class of three-strategy evolutionary games in which each strategy is neutral vs. itself, and has a positive expected payoff vs. one of the other strategies and a negative expected payoff vs. the remaining strategy. The payoff matrix is thus A = 0 b 2 a a 2 0 b 3 b a 3 0. (5) We perturb off of the canonical case, a = = b 3 =, by taking 0 F (t) + F (t) A = 0 0 (6) where the forcing function F is given by F (t) = ɛ(( δ) cos ω t + δ cos ω 2 t). (7) For ease of notation, write (x, x 2, x 3 ) = (x, y, z). The dynamics occur in the simplex S {(x, y, z) R x, y, z [0, } (8) but since x, y, z are the frequencies of the three strategies, and hence x + y + z =, we can eliminate z using z = x y. Therefore, the region of interest is T, the projection of S into the x y plane: See Figure. Thus the replicator equation () becomes Note that ẋ = 0 when x = 0, ẏ = 0 when y = 0, and T {(x, y) R x, y, x + y [0, }. (9) ẋ = x(x + 2y )( + (x )F (t)) (0) ẏ = y(2x + y x(x + 2y )F (t)) () ẋ + ẏ = (x + y )(xf (t)(x + 2y ) x + y) (2) so that ẋ+ẏ = 0 when x+y =, which means that x+y = is an invariant manifold. This shows that the boundary of T is invariant, so trajectories cannot escape the region of interest. It is known [4 that in the unperturbed case (ɛ = 0) there is an equilibrium point at (x, y) = ( 3, 3 ), and the interior of T is filled with periodic orbits. We see from equations (0)-() that this interior equilibrium point persists when ɛ 0. Numerical integration suggests that the Lyapunov stability of motions around the equilibrium point depends sensitively on the values of ω and ω 2. See Figure 2. We investigate the stability of the interior equilibrium using Floquet theory and harmonic balance, as well as by numerical methods.
3 z 0.0 x 0.0 y 0.0 Figure : A curve in S and its projection in T 2.2 Linearization To study the linear stability of the equilibrium point, we set x = u + 3, y = v + 3, substitute these into (0)-() and linearize, to obtain u = (u + 2v)(3 + 2F (t)) (3) 9 v = (F (t)(u + 2v) + 3(2u + v)). (4) 9 The linearized system (3)-(4) can also be written [6 as a single second-order equation on u, by differentiating (3) and substituting in expressions for v from (4) and v from (3). This gives us where g(t)ü ġ(t) u 9 g2 (t)u = 0 (5) g(t) = 3 2F (t) = 3 2ɛ(( δ) cos ω t + δ cos ω 2 t). (6) Now that we have a linear system with coefficients that are functions of time, we use Floquet theory to determine the stability of the origin. 3 Floquet theory Floquet theory is concerned with systems of differential equations of the form dx dt = M(t)x, M(t + T ) = M(t). (7) We have the system (3)-(4), which can be written as [ u = [ g(t) 2g(t) v 9 (9 g(t)) g(t) 2 [ u v [ u B(t) v (8) where g(t) is as in (6). In general, B(t) is not periodic, since ω and ω 2 are rationally independent. However, the set of points for which ω and ω 2 are rationally dependent is dense in the ω ω 2 plane, and solutions of (8) must vary continuously with ω and ω 2, so it is reasonable to consider only the case that F (t), and hence g(t) and B(t), are in fact periodic.
4 x x t t (a) ω = ω 2 =.2 (b) ω = ω 2 = 0.9 Figure 2: Numerical solutions for x(t) with identical initial conditions x(0) = y(0) = 0.33 and parameters ɛ = 0.9, δ = 0.6, but with different ω, ω 2 Assume that ω 2 = a b ω in lowest terms, where a and b are relatively prime integers. Then we can make the change of variables τ = ω t, so ω 2 t = a b τ. Since a and b are relatively prime, we see that F, and hence g and B, have period T = 2πb in τ. Thus (8) becomes [ u v = [ u B(τ), B(τ + 2πb) = B(τ) (9) ω v where u indicates du/dτ. This has the same form as (7), so we can apply the results of Floquet theory. Suppose that there is a fundamental solution matrix of (9), [ u (τ) u X(τ) = 2 (τ) v (τ) v 2 (τ) where [ u (0) v (0) = [ 0, [ u2 (0) v 2 (0) = [ 0 (20). (2) Then the Floquet matrix is C = X(T ) = X(2πb), and stability is determined by the eigenvalues of C: We can show [5 that det C =, as follows. Define the Wronskian λ 2 (trc)λ + det C = 0. (22) W (τ) = det X(τ) = u (τ)v 2 (τ) u 2 (τ)v (τ). (23) Notice that W (0) = det X(0) =. Then taking the time derivative of W and using (9) gives dw dτ = u (τ)v 2 (τ) + u (τ)v 2(τ) u 2(τ)v (τ) u 2 (τ)v (τ) = ( g(τ)(u + 2v )v 2 + 9ω 2 u (9u 2 (u 2 + 2v 2 )g(τ)) g(τ)(u 2 + 2v 2 )v ) 2 u 2(9u (u + 2v )g(τ)) = 0. (24) This shows that W (τ) = for all τ, and in particular W (T ) = det C =. Therefore, λ = trc ± trc 2 4 (25) 2 which means [5 that the transition between stable and unstable solutions occurs when trc = 2, and this corresponds to periodic solutions of period T = 2πb or 2T = 4πb. Given the period of the solutions on the transition curves in the ω ω 2 plane, we use harmonic balance to approximate those transition curves.
5 4 Harmonic balance We seek solutions to (5) of period 4πb in τ: u = α k cos k=0 ( ) kτ + β k sin 2b ( ) kτ. (26) 2b Since ω 2 = a b ω where a and b are relatively prime, any integer k can be written as na + mb for some integers n and m [, 2. That is, there is a one-to-one correspondence between integers k and ordered pairs (m, n). We can therefore write the solution as u = = m=0 n= m=0 n= ( ma + nb α mn cos 2b ( mω2 + nω α mn cos t 2 ) ( ma + nb τ + β mn sin 2b ) τ (27) ) ( ) mω2 + nω + β mn sin t. (28) 2 We substitute a truncated version of (28) into (5), expand the trigonometric functions and collect like terms. This results in cosine terms whose coefficients are functions of the α mn, and sine terms whose coefficients are functions of the β mn. Let the coefficient matrices of these two sets of terms be Q and R, respectively. In order for a nontrivial solution to exist, the determinants of both coefficient matrices must vanish [5. We solve the equations det Q = 0 and det R = 0 for relations between ω and ω 2. This gives the approximate transition curves seen in Figure 3. It has been shown [7, 6 that in a periodically forced RPS system (i.e. δ = 0 in our model) there are tongues of instability emerging from ω = 2/n 3 in the ω ɛ plane. Our harmonic balance analysis is consistent with this: we observe bands of instability around ω = 2/ 3 and ω 2 = 2/ 3, which get broader as ɛ increases. We also see narrower regions of instability along the lines nω + mω 2 = 2/ 3, for each n, m used in the truncated solution (28). Thus the boundary of the region of instability exhibits self-similarity when we consider ω, ω 2 [0, 2 k for k = 0,, Numerical integration In order to check the results of the harmonic balance method, we generate an approximate stability diagram by numerical integration of the linearized system (5). For randomly chosen parameters (ω, ω 2 ) [0, 2, we choose random initial conditions (u(0), u(0)) on the unit circle - since the system is linear, the amplitude of the initial condition needs only to be consistent between trials. We then integrate the system for 000 time steps using ode45 in Matlab. This is an explicit Runge-Kutta (4,5) method that is recommended in the Matlab documentation for most non-stiff problems. We considered a motion to be unstable if max u(t) > 0. The set of points (ω, ω 2 ) corresponding to unstable motions were plotted using matplotlib.pyplot in Python. See Figure 4. Each plot in Figure 4 contains approximately points. We note that the unstable regions given by numerical integration appear to be consistent with the transition curves predicted by harmonic balance (Figure 3). The regions of instability around ω = 2/ 3 and ω 2 = 2/ 3 are visible for all tested values of ɛ and δ, and as ɛ increases, more tongues of the form nω + mω 2 = 2/ 3 become visible. 6 Lyapunov exponents A second, and more informative, numerical approach for determining stability is the computation of approximate Lyapunov exponents. This is a measure of a solution s rate of divergence from the equilibrium point [, and is defined as λ = lim sup t ln u(t). (29) t If the limit is finite, then u(t) e λt or smaller as t. A positive Lyapunov exponent indicates that the solution is unstable.
6 2.0 Ω (a) ɛ =, δ = 0.6 Ω Ω (b) ɛ = 0.9, δ = 0.6 Ω Ω (c) ɛ =.3, δ = 0.6 Figure 3: Transition curves predicted by harmonic balance with 5 m 5, 0 n 5 for various values of ɛ
7 (a) ɛ =, δ = 0.6 (b) ɛ = 0.9, δ = 0.6 (c) ɛ =.3, δ = 0.6 (d) Detail view: ɛ =.3, δ = 0.6 Figure 4: Plots of unstable points in the (ω, ω 2 ) plane for various values of ɛ
8 We do not find any negative Lyapunov exponents, but note [5 that the system (5) can be converted to a Hill s equation ( 4g(t) ġ(t) 2 ) 8g(t) g(t) z z 36g(t) 2 = 0 (30) by making the change of variables u = g(t)z. Since g(t) is bounded, u is bounded if and only if z is bounded. And since there is no dissipation in (30), stable solutions correspond to λ = 0. We approximate the Lyapunov exponents numerically by integrating as above, and taking λ sup ln u(t). (3) 900<t<000 t See Figure 5. The shape of the unstable region is the same as in Figure 4, but this method allows us to see a sharp increase in unstable solutions rate of growth along the line ω = ω 2. 7 Conclusion The replicator equation with quasiperiodic perturbation may be used to model biological or social systems where competition is affected by cyclical processes on different scales. We have investigated the linear stability of the interior equilibrium point for RPS systems with quasiperiodic perturbation, using Floquet theory and harmonic balance, as well as numerical integration and numerical computation of Lyapunov exponents. We find that stability depends sensitively on the frequencies ω and ω 2, and that the region of instability in the ω ω 2 plane exhibits self-similarity. References [ Guckenheimer, J. and Holmes, P. Nonlinear Oscillations, Dynamical Systems, and Bifurcations of Vector Fields. Springer-Verlag, New York, 983. [2 Herstein, I. Topics in Algebra, 2nd ed., John Wiley & Sons, New York, 975. [3 Hofbauer, J. and Sigmund, K. Evolutionary Games and Population Dynamics. Cambridge University Press, Cambridge, UK, 998. [4 Nowak, M., Evolutionary Dynamics. Belknap Press of Harvard Univ. Press, Cambridge, MA. [5 Rand, R. Lecture Notes on Nonlinear Vibrations. Published online by The Internet-First University Press, [6 Rand, R., Yazhbin, M., and Rand, D. Evolutionary dynamics of a system with periodic coefficients. Comm Nonlinear Sci Numer Simulat, 6: , 20. [7 Ruelas, R., Rand, D., and Rand, R. Nonlinear parametric excitation of an evolutionary dynamical system. J. Mechanical Engineering science (Proceedings Inst. Mech. Eng., Part C), 226:92-920, 202. [8 Ruelas, R., Rand, D., and Rand, R. Parametric Excitation and Evolutionary Dynamics. J. Applied Mechanics 80, 0503, 203. [9 Sigmund, K. Introduction to Evolutionary Game Theory. In Evolutionary Game Dynamics, K. Sigmund, ed., Vol. 69 of Proceedings of Symposia in Applied Mathematics. American Mathematical Society, Chap., pp. -26, 20. [0 Taylor, P. and Jonker, L. Evolutionarily Stable Strategies and Game Dynamics. Mathematical Biosciences 40(-2), pp , 978. [ Zounes, R. and Rand, R. Transition Curves for the Quasi-periodic Mathieu Equation. Siam J. Appl. Math, 40(4), pp , 998.
9 (a) ɛ =, δ = 0.6. Contours between λ = 0 and λ = (b) ɛ = 0.9, δ = 0.6. Contours between λ = 0 and λ = (c) ɛ =.3, δ = 0.6. Contours between λ = 0 and λ = 0.2. (d) Detail view: ɛ =.3, δ = 0.6. Contours as in (c). Figure 5: Contour plot of Lyapunov exponents in the (ω, ω 2 ) plane for various values of ɛ
Available online at ScienceDirect. Procedia IUTAM 19 (2016 ) IUTAM Symposium Analytical Methods in Nonlinear Dynamics
Available online at www.sciencedirect.com ScienceDirect Procedia IUTAM 19 (2016 ) 11 18 IUTAM Symposium Analytical Methods in Nonlinear Dynamics A model of evolutionary dynamics with quasiperiodic forcing
More informationJournal of Applied Nonlinear Dynamics
Journal of Applied Nonlinear Dynamics 4(2) (2015) 131 140 Journal of Applied Nonlinear Dynamics https://lhscientificpublishing.com/journals/jand-default.aspx A Model of Evolutionary Dynamics with Quasiperiodic
More informationCommun Nonlinear Sci Numer Simulat
Commun Nonlinear Sci Numer Simulat 16 (011) 887 895 Contents lists available at ScienceDirect Commun Nonlinear Sci Numer Simulat journal homepage: www.elsevier.com/locate/cnsns Short communication Evolutionary
More information2:2:1 Resonance in the Quasiperiodic Mathieu Equation
Nonlinear Dynamics 31: 367 374, 003. 003 Kluwer Academic Publishers. Printed in the Netherlands. ::1 Resonance in the Quasiperiodic Mathieu Equation RICHARD RAND Department of Theoretical and Applied Mechanics,
More informationCoupled Parametrically Driven Modes in Synchrotron Dynamics
paper presented at Society for Experimental Mechanics SEM) IMAC XXXIII Conference on Structural Dynamics February -5, 015, Orlando FL Coupled Parametrically Driven Modes in Synchrotron Dynamics Alexander
More information2:1:1 Resonance in the Quasi-Periodic Mathieu Equation
Nonlinear Dynamics 005) 40: 95 03 c pringer 005 :: Resonance in the Quasi-Periodic Mathieu Equation RICHARD RAND and TINA MORRION Department of Theoretical & Applied Mechanics, Cornell University, Ithaca,
More informationHigher Order Averaging : periodic solutions, linear systems and an application
Higher Order Averaging : periodic solutions, linear systems and an application Hartono and A.H.P. van der Burgh Faculty of Information Technology and Systems, Department of Applied Mathematical Analysis,
More informationEvolutionary Games and Periodic Fitness
Evolutionary Games and Periodic Fitness Philippe Uyttendaele Frank Thuijsman Pieter Collins Ralf Peeters Gijs Schoenmakers Ronald Westra March 16, 2012 Abstract One thing that nearly all stability concepts
More informationCoexistence phenomenon in autoparametric excitation of two
International Journal of Non-Linear Mechanics 40 (005) 60 70 wwwelseviercom/locate/nlm Coexistence phenomenon in autoparametric excitation of two degree of freedom systems Geoffrey Recktenwald, Richard
More informationEntrance Exam, Differential Equations April, (Solve exactly 6 out of the 8 problems) y + 2y + y cos(x 2 y) = 0, y(0) = 2, y (0) = 4.
Entrance Exam, Differential Equations April, 7 (Solve exactly 6 out of the 8 problems). Consider the following initial value problem: { y + y + y cos(x y) =, y() = y. Find all the values y such that the
More informationTwo models for the parametric forcing of a nonlinear oscillator
Nonlinear Dyn (007) 50:147 160 DOI 10.1007/s11071-006-9148-3 ORIGINAL ARTICLE Two models for the parametric forcing of a nonlinear oscillator Nazha Abouhazim Mohamed Belhaq Richard H. Rand Received: 3
More information8 Example 1: The van der Pol oscillator (Strogatz Chapter 7)
8 Example 1: The van der Pol oscillator (Strogatz Chapter 7) So far we have seen some different possibilities of what can happen in two-dimensional systems (local and global attractors and bifurcations)
More informationNonlinear parametric excitation of an evolutionary dynamical system
Special Issue Article Nonlinear parametric excitation of an evolutionary dynamical system Proc IMechE Part C: J Mechanical Engineering Science 226(8) 1912 1920! IMechE 2011 Reprints and permissions: sagepub.co.uk/journalspermissions.nav
More informationLinear and Nonlinear Oscillators (Lecture 2)
Linear and Nonlinear Oscillators (Lecture 2) January 25, 2016 7/441 Lecture outline A simple model of a linear oscillator lies in the foundation of many physical phenomena in accelerator dynamics. A typical
More information2.10 Saddles, Nodes, Foci and Centers
2.10 Saddles, Nodes, Foci and Centers In Section 1.5, a linear system (1 where x R 2 was said to have a saddle, node, focus or center at the origin if its phase portrait was linearly equivalent to one
More informationCALCULATION OF NONLINEAR VIBRATIONS OF PIECEWISE-LINEAR SYSTEMS USING THE SHOOTING METHOD
Vietnam Journal of Mechanics, VAST, Vol. 34, No. 3 (2012), pp. 157 167 CALCULATION OF NONLINEAR VIBRATIONS OF PIECEWISE-LINEAR SYSTEMS USING THE SHOOTING METHOD Nguyen Van Khang, Hoang Manh Cuong, Nguyen
More information6 Linear Equation. 6.1 Equation with constant coefficients
6 Linear Equation 6.1 Equation with constant coefficients Consider the equation ẋ = Ax, x R n. This equating has n independent solutions. If the eigenvalues are distinct then the solutions are c k e λ
More informationCOMPLEX DYNAMICS AND CHAOS CONTROL IN DUFFING-VAN DER POL EQUATION WITH TWO EXTERNAL PERIODIC FORCING TERMS
International J. of Math. Sci. & Engg. Appls. (IJMSEA) ISSN 0973-9424, Vol. 9 No. III (September, 2015), pp. 197-210 COMPLEX DYNAMICS AND CHAOS CONTROL IN DUFFING-VAN DER POL EQUATION WITH TWO EXTERNAL
More information7 Planar systems of linear ODE
7 Planar systems of linear ODE Here I restrict my attention to a very special class of autonomous ODE: linear ODE with constant coefficients This is arguably the only class of ODE for which explicit solution
More informationChapter #4 EEE8086-EEE8115. Robust and Adaptive Control Systems
Chapter #4 Robust and Adaptive Control Systems Nonlinear Dynamics.... Linear Combination.... Equilibrium points... 3 3. Linearisation... 5 4. Limit cycles... 3 5. Bifurcations... 4 6. Stability... 6 7.
More informationMathematical Foundations of Neuroscience - Lecture 7. Bifurcations II.
Mathematical Foundations of Neuroscience - Lecture 7. Bifurcations II. Filip Piękniewski Faculty of Mathematics and Computer Science, Nicolaus Copernicus University, Toruń, Poland Winter 2009/2010 Filip
More informationLagrangian Analysis of 2D and 3D Ocean Flows from Eulerian Velocity Data
Flows from Second-year Ph.D. student, Applied Math and Scientific Computing Project Advisor: Kayo Ide Department of Atmospheric and Oceanic Science Center for Scientific Computation and Mathematical Modeling
More informationChapter III. Stability of Linear Systems
1 Chapter III Stability of Linear Systems 1. Stability and state transition matrix 2. Time-varying (non-autonomous) systems 3. Time-invariant systems 1 STABILITY AND STATE TRANSITION MATRIX 2 In this chapter,
More informationApplied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.
Printed Name: Signature: Applied Math Qualifying Exam 11 October 2014 Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. 2 Part 1 (1) Let Ω be an open subset of R
More informationDISCRETIZED BEST-RESPONSE DYNAMICS FOR THE ROCK-PAPER-SCISSORS GAME. Peter Bednarik. Josef Hofbauer. (Communicated by William H.
Journal of Dynamics and Games doi:0.94/jdg.207005 American Institute of Mathematical Sciences Volume 4, Number, January 207 pp. 75 86 DISCRETIZED BEST-RESPONSE DYNAMICS FOR THE ROCK-PAPER-SCISSORS GAME
More informationHalf of Final Exam Name: Practice Problems October 28, 2014
Math 54. Treibergs Half of Final Exam Name: Practice Problems October 28, 24 Half of the final will be over material since the last midterm exam, such as the practice problems given here. The other half
More informationDynamics of Two Coupled van der Pol Oscillators with Delay Coupling Revisited
Dynamics of Two Coupled van der Pol Oscillators with Delay Coupling Revisited arxiv:1705.03100v1 [math.ds] 8 May 017 Mark Gluzman Center for Applied Mathematics Cornell University and Richard Rand Dept.
More informationIdentification Methods for Structural Systems. Prof. Dr. Eleni Chatzi Lecture March, 2016
Prof. Dr. Eleni Chatzi Lecture 4-09. March, 2016 Fundamentals Overview Multiple DOF Systems State-space Formulation Eigenvalue Analysis The Mode Superposition Method The effect of Damping on Structural
More informationMath Ordinary Differential Equations
Math 411 - Ordinary Differential Equations Review Notes - 1 1 - Basic Theory A first order ordinary differential equation has the form x = f(t, x) (11) Here x = dx/dt Given an initial data x(t 0 ) = x
More informationLogistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations
Logistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations S. Y. Ha and J. Park Department of Mathematical Sciences Seoul National University Sep 23, 2013 Contents 1 Logistic Map 2 Euler and
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More information2:1 Resonance in the delayed nonlinear Mathieu equation
Nonlinear Dyn 007) 50:31 35 DOI 10.1007/s11071-006-916-5 ORIGINAL ARTICLE :1 Resonance in the delayed nonlinear Mathieu equation Tina M. Morrison Richard H. Rand Received: 9 March 006 / Accepted: 19 September
More informationDynamics of a mass-spring-pendulum system with vastly different frequencies
Dynamics of a mass-spring-pendulum system with vastly different frequencies Hiba Sheheitli, hs497@cornell.edu Richard H. Rand, rhr2@cornell.edu Cornell University, Ithaca, NY, USA Abstract. We investigate
More informationIrrational behavior in the Brown von Neumann Nash dynamics
Irrational behavior in the Brown von Neumann Nash dynamics Ulrich Berger a and Josef Hofbauer b a Vienna University of Economics and Business Administration, Department VW 5, Augasse 2-6, A-1090 Wien,
More information3 Stability and Lyapunov Functions
CDS140a Nonlinear Systems: Local Theory 02/01/2011 3 Stability and Lyapunov Functions 3.1 Lyapunov Stability Denition: An equilibrium point x 0 of (1) is stable if for all ɛ > 0, there exists a δ > 0 such
More informationModeling and Analysis of Dynamic Systems
Modeling and Analysis of Dynamic Systems Dr. Guillaume Ducard Fall 2017 Institute for Dynamic Systems and Control ETH Zurich, Switzerland G. Ducard c 1 / 57 Outline 1 Lecture 13: Linear System - Stability
More informationChapter 14 Three Ways of Treating a Linear Delay Differential Equation
Chapter 14 Three Ways of Treating a Linear Delay Differential Equation Si Mohamed Sah and Richard H. Rand Abstract This work concerns the occurrence of Hopf bifurcations in delay differential equations
More information2D-Volterra-Lotka Modeling For 2 Species
Majalat Al-Ulum Al-Insaniya wat - Tatbiqiya 2D-Volterra-Lotka Modeling For 2 Species Alhashmi Darah 1 University of Almergeb Department of Mathematics Faculty of Science Zliten Libya. Abstract The purpose
More informationL = 1 2 a(q) q2 V (q).
Physics 3550, Fall 2011 Motion near equilibrium - Small Oscillations Relevant Sections in Text: 5.1 5.6 Motion near equilibrium 1 degree of freedom One of the most important situations in physics is motion
More informationGeometry of Resonance Tongues
Geometry of Resonance Tongues Henk Broer with Martin Golubitsky, Sijbo Holtman, Mark Levi, Carles Simó & Gert Vegter Instituut voor Wiskunde en Informatica Rijksuniversiteit Groningen Resonance p.1/36
More informationNONSTANDARD NUMERICAL METHODS FOR A CLASS OF PREDATOR-PREY MODELS WITH PREDATOR INTERFERENCE
Sixth Mississippi State Conference on Differential Equations and Computational Simulations, Electronic Journal of Differential Equations, Conference 15 (2007), pp. 67 75. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu
More informationDETC DELAY DIFFERENTIAL EQUATIONS IN THE DYNAMICS OF GENE COPYING
Proceedings of the ASME 2007 International Design Engineering Technical Conferences & Computers and Information in Engineering Conference IDETC/CIE 2007 September 4-7, 2007, Las Vegas, Nevada, USA DETC2007-3424
More informationStability Analysis of a Hydrodynamic Journal Bearing With Rotating Herringbone Grooves
G. H. Jang e-mail: ghjang@hanyang.ac.kr J. W. Yoon PREM, Department of Mechanical Engineering, Hanyang University, Seoul, 33-79, Korea Stability Analysis of a Hydrodynamic Journal Bearing With Rotating
More informationClassification of Phase Portraits at Equilibria for u (t) = f( u(t))
Classification of Phase Portraits at Equilibria for u t = f ut Transfer of Local Linearized Phase Portrait Transfer of Local Linearized Stability How to Classify Linear Equilibria Justification of the
More informationHardening nonlinearity effects on forced vibration of viscoelastic dynamical systems to external step perturbation field
Int. J. Adv. Appl. Math. and Mech. 3(2) (215) 16 32 (ISSN: 2347-2529) IJAAMM Journal homepage: www.ijaamm.com International Journal of Advances in Applied Mathematics and Mechanics Hardening nonlinearity
More informationThe Big Picture. Discuss Examples of unpredictability. Odds, Stanisław Lem, The New Yorker (1974) Chaos, Scientific American (1986)
The Big Picture Discuss Examples of unpredictability Odds, Stanisław Lem, The New Yorker (1974) Chaos, Scientific American (1986) Lecture 2: Natural Computation & Self-Organization, Physics 256A (Winter
More informationSTABILITY. Phase portraits and local stability
MAS271 Methods for differential equations Dr. R. Jain STABILITY Phase portraits and local stability We are interested in system of ordinary differential equations of the form ẋ = f(x, y), ẏ = g(x, y),
More informationCDS 101 Precourse Phase Plane Analysis and Stability
CDS 101 Precourse Phase Plane Analysis and Stability Melvin Leok Control and Dynamical Systems California Institute of Technology Pasadena, CA, 26 September, 2002. mleok@cds.caltech.edu http://www.cds.caltech.edu/
More informationEvolutionary Dynamics of Lewis Signaling Games: Signaling Systems vs. Partial Pooling
This version: October 8, 2006. Introduction Evolutionary Dynamics of Lewis Signaling Games: Signaling Systems vs. Partial Pooling Simon Huttegger Brian Skyrms Rory Smead Kevin Zollman In Lewis signaling
More informationDynamicsofTwoCoupledVanderPolOscillatorswithDelayCouplingRevisited
Global Journal of Science Frontier Research: F Mathematics and Decision Sciences Volume 7 Issue 5 Version.0 Year 07 Type : Double Blind Peer Reviewed International Research Journal Publisher: Global Journals
More informationProf. Krstic Nonlinear Systems MAE281A Homework set 1 Linearization & phase portrait
Prof. Krstic Nonlinear Systems MAE28A Homework set Linearization & phase portrait. For each of the following systems, find all equilibrium points and determine the type of each isolated equilibrium. Use
More informationB5.6 Nonlinear Systems
B5.6 Nonlinear Systems 1. Linear systems Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Linear systems 1.1 Differential Equations 1.2 Linear flows 1.3 Linear maps
More informationPoincaré Map, Floquet Theory, and Stability of Periodic Orbits
Poincaré Map, Floquet Theory, and Stability of Periodic Orbits CDS140A Lecturer: W.S. Koon Fall, 2006 1 Poincaré Maps Definition (Poincaré Map): Consider ẋ = f(x) with periodic solution x(t). Construct
More informationEvolutionary Dynamics and Extensive Form Games by Ross Cressman. Reviewed by William H. Sandholm *
Evolutionary Dynamics and Extensive Form Games by Ross Cressman Reviewed by William H. Sandholm * Noncooperative game theory is one of a handful of fundamental frameworks used for economic modeling. It
More informationWHAT IS A CHAOTIC ATTRACTOR?
WHAT IS A CHAOTIC ATTRACTOR? CLARK ROBINSON Abstract. Devaney gave a mathematical definition of the term chaos, which had earlier been introduced by Yorke. We discuss issues involved in choosing the properties
More informationQuasiperiodic phenomena in the Van der Pol - Mathieu equation
Quasiperiodic penomena in te Van der Pol - Matieu equation F. Veerman and F. Verulst Department of Matematics, Utrect University P.O. Box 80.010, 3508 TA Utrect Te Neterlands April 8, 009 Abstract Te Van
More informationDRIVEN and COUPLED OSCILLATORS. I Parametric forcing The pendulum in 1:2 resonance Santiago de Compostela
DRIVEN and COUPLED OSCILLATORS I Parametric forcing The pendulum in 1:2 resonance Santiago de Compostela II Coupled oscillators Resonance tongues Huygens s synchronisation III Coupled cell system with
More informationarxiv: v1 [quant-ph] 9 Sep 2011
Stability analysis of surface ion traps Arkadas Ozakin and Fayaz Shaikh Quantum Information Systems Group Georgia Tech Research Institute March 24, 2018 arxiv:1109.2160v1 [quant-ph] 9 Sep 2011 Abstract
More informationFeedback control for a chemostat with two organisms
Feedback control for a chemostat with two organisms Patrick De Leenheer and Hal Smith Arizona State University Department of Mathematics and Statistics Tempe, AZ 85287 email: leenheer@math.la.asu.edu,
More informationSymmetry Properties of Confined Convective States
Symmetry Properties of Confined Convective States John Guckenheimer Cornell University 1 Introduction This paper is a commentary on the experimental observation observations of Bensimon et al. [1] of convection
More informationLMI Methods in Optimal and Robust Control
LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 15: Nonlinear Systems and Lyapunov Functions Overview Our next goal is to extend LMI s and optimization to nonlinear
More informationUnderstand the existence and uniqueness theorems and what they tell you about solutions to initial value problems.
Review Outline To review for the final, look over the following outline and look at problems from the book and on the old exam s and exam reviews to find problems about each of the following topics.. Basics
More informationComplicated behavior of dynamical systems. Mathematical methods and computer experiments.
Complicated behavior of dynamical systems. Mathematical methods and computer experiments. Kuznetsov N.V. 1, Leonov G.A. 1, and Seledzhi S.M. 1 St.Petersburg State University Universitetsky pr. 28 198504
More informationContinuous Threshold Policy Harvesting in Predator-Prey Models
Continuous Threshold Policy Harvesting in Predator-Prey Models Jon Bohn and Kaitlin Speer Department of Mathematics, University of Wisconsin - Madison Department of Mathematics, Baylor University July
More informationApplied Dynamical Systems
Applied Dynamical Systems Recommended Reading: (1) Morris W. Hirsch, Stephen Smale, and Robert L. Devaney. Differential equations, dynamical systems, and an introduction to chaos. Elsevier/Academic Press,
More informationODE Homework 1. Due Wed. 19 August 2009; At the beginning of the class
ODE Homework Due Wed. 9 August 2009; At the beginning of the class. (a) Solve Lẏ + Ry = E sin(ωt) with y(0) = k () L, R, E, ω are positive constants. (b) What is the limit of the solution as ω 0? (c) Is
More informationCDS 101/110a: Lecture 2.1 Dynamic Behavior
CDS 11/11a: Lecture 2.1 Dynamic Behavior Richard M. Murray 6 October 28 Goals: Learn to use phase portraits to visualize behavior of dynamical systems Understand different types of stability for an equilibrium
More informationCalculus and Differential Equations II
MATH 250 B Second order autonomous linear systems We are mostly interested with 2 2 first order autonomous systems of the form { x = a x + b y y = c x + d y where x and y are functions of t and a, b, c,
More informationKey words. quasi-periodic, Floquet theory, Hill s equation, Mathieu equation, perturbations, stability
SIAM J. APPL. MATH. c 1998 Society for Industrial and Applied Mathematics Vol. 58, No. 4, pp. 1094 1115, August 1998 004 TRANSITION CURVES FOR THE QUASI-PERIODIC MATHIEU EQUATION RANDOLPH S. ZOUNES AND
More informationAn Introduction to Evolutionary Game Theory: Lecture 2
An Introduction to Evolutionary Game Theory: Lecture 2 Mauro Mobilia Lectures delivered at the Graduate School on Nonlinear and Stochastic Systems in Biology held in the Department of Applied Mathematics,
More informationMATH 23 Exam 2 Review Solutions
MATH 23 Exam 2 Review Solutions Problem 1. Use the method of reduction of order to find a second solution of the given differential equation x 2 y (x 0.1875)y = 0, x > 0, y 1 (x) = x 1/4 e 2 x Solution
More informationThree ways of treating a linear delay differential equation
Proceedings of the 5th International Conference on Nonlinear Dynamics ND-KhPI2016 September 27-30, 2016, Kharkov, Ukraine Three ways of treating a linear delay differential equation Si Mohamed Sah 1 *,
More informationDYNAMICS OF THREE COUPLED VAN DER POL OSCILLATORS WITH APPLICATION TO CIRCADIAN RHYTHMS
Proceedings of IDETC/CIE 2005 ASME 2005 International Design Engineering Technical Conferences & Computers and Information in Engineering Conference September 24-28, 2005, Long Beach, California USA DETC2005-84017
More informationMCE693/793: Analysis and Control of Nonlinear Systems
MCE693/793: Analysis and Control of Nonlinear Systems Systems of Differential Equations Phase Plane Analysis Hanz Richter Mechanical Engineering Department Cleveland State University Systems of Nonlinear
More information1 Introduction. Minho Lee 1 Jihoon Lee 1 Gunhee Jang 1
DOI 10.1007/s005-015-5-5 TECHNICAL PAPER Stability analysis of a whirling rigid rotor supported by stationary grooved FDBs considering the five degrees of freedom of a general rotor bearing system Minho
More informationMechanisms of Chaos: Stable Instability
Mechanisms of Chaos: Stable Instability Reading for this lecture: NDAC, Sec. 2.-2.3, 9.3, and.5. Unpredictability: Orbit complicated: difficult to follow Repeatedly convergent and divergent Net amplification
More informationMath 331 Homework Assignment Chapter 7 Page 1 of 9
Math Homework Assignment Chapter 7 Page of 9 Instructions: Please make sure to demonstrate every step in your calculations. Return your answers including this homework sheet back to the instructor as a
More informationTHREE DIMENSIONAL SYSTEMS. Lecture 6: The Lorenz Equations
THREE DIMENSIONAL SYSTEMS Lecture 6: The Lorenz Equations 6. The Lorenz (1963) Equations The Lorenz equations were originally derived by Saltzman (1962) as a minimalist model of thermal convection in a
More informationThe goal of this chapter is to study linear systems of ordinary differential equations: dt,..., dx ) T
1 1 Linear Systems The goal of this chapter is to study linear systems of ordinary differential equations: ẋ = Ax, x(0) = x 0, (1) where x R n, A is an n n matrix and ẋ = dx ( dt = dx1 dt,..., dx ) T n.
More informationHILL S EQUATION WITH RANDOM FORCING TERMS. Fred C. Adams and Anthony M. Bloch Also: Suzie Butler, Jeff Druce and Jake Ketchum
1 HILL S EQUATION WITH RANDOM FORCING TERMS Fred C. Adams and Anthony M. Bloch Also: Suzie Butler, Jeff Druce and Jake Ketchum 2 Orbit problems in astrophysics Solutions to Hill s equation with forcing
More informationnonlinear oscillators. method of averaging
Physics 4 Spring 7 nonlinear oscillators. method of averaging lecture notes, spring semester 7 http://www.phys.uconn.edu/ rozman/courses/p4_7s/ Last modified: April, 7 Oscillator with nonlinear friction
More information2012 NCTS Workshop on Dynamical Systems
Barbara Gentz gentz@math.uni-bielefeld.de http://www.math.uni-bielefeld.de/ gentz 2012 NCTS Workshop on Dynamical Systems National Center for Theoretical Sciences, National Tsing-Hua University Hsinchu,
More informationMassachusetts Institute of Technology Physics Department
Massachusetts Institute of Technology Physics Department Physics 8.32 Fall 2006 Quantum Theory I October 9, 2006 Assignment 6 Due October 20, 2006 Announcements There will be a makeup lecture on Friday,
More informationIrrational behavior in the Brown von Neumann Nash dynamics
Games and Economic Behavior 56 (2006) 1 6 www.elsevier.com/locate/geb Irrational behavior in the Brown von Neumann Nash dynamics Ulrich Berger a,, Josef Hofbauer b a Vienna University of Economics and
More informationApplied Differential Equation. November 30, 2012
Applied Differential Equation November 3, Contents 5 System of First Order Linear Equations 5 Introduction and Review of matrices 5 Systems of Linear Algebraic Equations, Linear Independence, Eigenvalues,
More informationStudies on Rayleigh Equation
Johan Matheus Tuwankotta Studies on Rayleigh Equation Intisari Dalam makalah ini dipelajari persamaan oscilator tak linear yaitu persamaan Rayleigh : x x = µ ( ( x ) ) x. Masalah kestabilan lokal di sekitar
More informationControl Systems Design
ELEC4410 Control Systems Design Lecture 13: Stability Julio H. Braslavsky julio@ee.newcastle.edu.au School of Electrical Engineering and Computer Science Lecture 13: Stability p.1/20 Outline Input-Output
More informationTranscendental cases in stability problem. Hamiltonian systems
of Hamiltonian systems Boris S. Bardin Moscow Aviation Institute (Technical University) Faculty of Applied Mathematics and Physics Department of Theoretical Mechanics Hamiltonian Dynamics and Celestial
More informationLecture 9. Systems of Two First Order Linear ODEs
Math 245 - Mathematics of Physics and Engineering I Lecture 9. Systems of Two First Order Linear ODEs January 30, 2012 Konstantin Zuev (USC) Math 245, Lecture 9 January 30, 2012 1 / 15 Agenda General Form
More informationLecture 10: Singular Perturbations and Averaging 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.243j (Fall 2003): DYNAMICS OF NONLINEAR SYSTEMS by A. Megretski Lecture 10: Singular Perturbations and
More informationThe Big Picture. Python environment? Discuss Examples of unpredictability. Chaos, Scientific American (1986)
The Big Picture Python environment? Discuss Examples of unpredictability Email homework to me: chaos@cse.ucdavis.edu Chaos, Scientific American (1986) Odds, Stanislaw Lem, The New Yorker (1974) 1 Nonlinear
More informationB5.6 Nonlinear Systems
B5.6 Nonlinear Systems 5. Global Bifurcations, Homoclinic chaos, Melnikov s method Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Motivation 1.1 The problem 1.2 A
More informationCDS 101/110a: Lecture 2.1 Dynamic Behavior
CDS 11/11a: Lecture.1 Dynamic Behavior Richard M. Murray 6 October 8 Goals: Learn to use phase portraits to visualize behavior of dynamical systems Understand different types of stability for an equilibrium
More informationControl Systems. Time response. L. Lanari
Control Systems Time response L. Lanari outline zero-state solution matrix exponential total response (sum of zero-state and zero-input responses) Dirac impulse impulse response change of coordinates (state)
More informationSystem Control Engineering 0
System Control Engineering 0 Koichi Hashimoto Graduate School of Information Sciences Text: Nonlinear Control Systems Analysis and Design, Wiley Author: Horacio J. Marquez Web: http://www.ic.is.tohoku.ac.jp/~koichi/system_control/
More informationTwo-Body Problem. Central Potential. 1D Motion
Two-Body Problem. Central Potential. D Motion The simplest non-trivial dynamical problem is the problem of two particles. The equations of motion read. m r = F 2, () We already know that the center of
More informationEN Nonlinear Control and Planning in Robotics Lecture 3: Stability February 4, 2015
EN530.678 Nonlinear Control and Planning in Robotics Lecture 3: Stability February 4, 2015 Prof: Marin Kobilarov 0.1 Model prerequisites Consider ẋ = f(t, x). We will make the following basic assumptions
More informationStability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games
Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,
More informationLecture 4. Chapter 4: Lyapunov Stability. Eugenio Schuster. Mechanical Engineering and Mechanics Lehigh University.
Lecture 4 Chapter 4: Lyapunov Stability Eugenio Schuster schuster@lehigh.edu Mechanical Engineering and Mechanics Lehigh University Lecture 4 p. 1/86 Autonomous Systems Consider the autonomous system ẋ
More information6 Parametric oscillator
6 Parametric oscillator 6. Mathieu equation We now study a different kind of forced pendulum. Specifically, imagine subjecting the pivot of a simple frictionless pendulum to an alternating vertical motion:
More information