Nonlinear dynamics & chaos BECS


 Ira Willis
 2 years ago
 Views:
Transcription
1 Nonlinear dynamics & chaos BECS
2 Phase portraits Focus: nonlinear systems in two dimensions General form of a vector field on the phase plane: Vector notation:
3 Phase portraits Solution x(t) describes a trajectory on the phase plane The whole plane is filled with (nonintersecting) trajectories For nonlinear systems, there is no hope to find trajectories analytically Our approach: determining the qualitative behavior of the solutions!
4 Phase portraits Huge variety of possible phase portraits Salient features: 1) Fixed points (A, B, C): equilibria of the system 2) Closed orbits (D): periodic solutions 3) Arrangement of trajectories near fixed points and closed orbits 4) Stability or instability of the fixed points and closed orbits
5 Numerical computation of RungeKutta method phase portraits A stepsize Δt=0.1 usually provides sufficient accuracy
6 Example I Fixed points Stability Exponentially growing solutions: the fixed point is unstable
7 Example I Phase portrait: nullclines The nullclines are the curves where On the nullclines the flow is either purely horizontal or purely vertical
8 Example I The fixed point is a nonlinear version of a saddle point
9 Existence, uniqueness and topological consequences Are we sure that there is a solution? Existence and Uniqueness Theorem: let f be continuous and all its partial derivatives be continuous for x in some open connected set. Then for the initial value problem has a solution x(t) on some time interval (τ, τ) about t=0, and the solution is unique
10 Existence, uniqueness and topological consequences Corollary: different trajectories never intersect! If two trajectories did intersect there would be two solutions starting from the same point (the crossing point)
11 Existence, uniqueness and topological consequences Consequence: any trajectory starting from inside a closed orbit will be trapped inside it forever! What will happen in the limit of t >? The trajectory will either converge to a fixed point or to a periodic orbit!
12 Fixed points and linearization Fixed point (x *, y * ) Does the disturbance from the fixed point grow or decay?
13 Fixed points and linearization Likewise:
14 Fixed points and linearization Jacobian matrix at the fixed point (x *, y * ), multivariate analog of the derivative f (x*) for 1dimensional systems If we neglect secondorder terms and higher, we obtain the linearized system
15 Fixed points and linearization Dynamics can be analyzed with the methods used for linear systems Effect of small nonlinear terms If the fixed point is not one of the borderline cases (centers, degenerate nodes, stars, nonisolated fixed points) the predicted type of the linearized system is the correct one
16 Example I Fixed points: (0,0), (1,0), (1,0) (0,0) > (±1,0) > Stable node Saddle points
17 Example I 1) x and y equations are uncoupled 2) ydirection: trajectories decay exponentially to y=0 3) xdirection: trajectories are attracted to x=0 and repelled by x= ± 1 4) Vertical lines x=0 and x = ± 1 are invariant: if the particle is there it will always stay there 5) The axis y=0 is invariant too x
18 Example I Phase portrait is symmetric with respect to the x and the y axes
19 Example II (0,0) is a fixed point > linearization τ = 0, Δ=1 > 0 > the fixed point is a center (borderline case) Is that so? Let us switch to polar coordinates
20 Example II
21 Example II
22 Example II Radial and angular motions are independent The fixed point is a spiral (stable for a < 0, unstable for a > 0) Centers are delicate: one needs that the orbit closes perfectly after one cycle, the slightest perturbation turns it into a spiral
23 Fixed points and linearization Stars and degenerate nodes can be altered by small nonlinearities, but their stability does not change! Stars and degenerate nodes stay in fact well within regions of stability or instability: small perturbations will leave them in those areas.
24 Fixed points and linearization Robust cases 1) Repellers (or sources): both eigenvalues have positive real part 2) Attractors (or sinks): both eigenvalues have negative real part 3) Saddles: one eigenvalue is positive, the other is negative Marginal cases 1) Centers: both eigenvalues are purely imaginary 2) Higherorder and nonisolated fixed points: at least one eigenvalue is zero Marginal cases are those where at least one eigenvalue satisfies Re(λ) = 0
25 Fixed points and linearization If Re(λ) 0 for both eigenvalues, the fixed point is called hyperbolic: in this case its type is predicted by the linearization Same thing in higherorder systems HartmanGrobman Theorem: the local phase portrait near a hyperbolic fixed point is topologically equivalent to the phase portrait of the linearization
26 Fixed points and linearization Intuition: two phase portraits are equivalent if one is a distorted version of the other What type of distortion? Bending, warping, but not ripping Consequences: closed orbits stay closed, trajectories connecting saddle points must not be broken, etc. A phase portrait is structurally stable if its topology cannot be changed by an arbitrarily small perturbation to the vector field Examples: the phase portrait of a saddle point is structurally stable, that of a center is not, as a small perturbation converts the center into a spiral
27 Rabbits versus Sheep
28 Rabbits versus Sheep LotkaVolterra model of competition between two species Rabbits and sheep competing for the same limited resource (e.g. grass): no predators, seasonal effects, etc. 1) Each species would grow to its carrying capacity in the absence of the other > logistic growth 2) When rabbits and sheep encounter each other, trouble starts: sheep kick rabbits away > conflicts occur at a rate proportional to the size of each population, reducing the growth rate for each species 3) Rabbits reproduce faster but they are more penalized by conflicts
29 Rabbits versus Sheep x(t) 0 > population of rabbits y(t) 0 > population of sheep Fixed points (0,0), (0,2), (3,0), (1,1)
30 Rabbits versus Sheep (0,0) Eigenvalues are λ= 2, 3 > the origin is an unstable node Trajectories near a node are tangential to the slower eigendirection (here the yaxis, for which λ = 2 < 3)
31 Rabbits versus Sheep (0,2) Eigenvalues are λ= 1, 2 > (0,2) is a stable node Trajectories near a node are tangential to the slower eigendirection [here v = (1, 2), for which λ = 1 > 1 <  2 ]
32 Rabbits versus Sheep (3,0) Eigenvalues are λ= 3, 1 > (3,0) is a stable node Trajectories near a node are tangential to the slower eigendirection [here v = (3, 1), for which λ = 1 > 1 <  3 )]
33 Rabbits versus Sheep (1,1) Eigenvalues are λ= > (1,1) is a saddle point
34 Rabbits versus Sheep The xy axes contain horizontal/vertical trajectories
35 Rabbits versus Sheep Biological interpretation: one species drives the other to extinction!
36 Rabbits versus Sheep Principle of competitive exclusion: two species competing for the same limited resource typically cannot coexist Basin of attraction of an attracting fixed point: set of initial conditions leading to that fixed point as t >
37 Conservative systems Equation of motion of a mass subject to a nonlinear force F(x): F(x) has no dependence on the velocity or time > no damping or friction, no timedependent driving force Energy is conserved V(x) is the potential energy
38 Conservative systems Systems with a conserved quantity are called conservative General definition: given a system a conserved quantity is a realvalued continuous function E(x) that is constant on trajectories (de/dt=0), but nonconstant on every open set
39 Example I A conservative system cannot have any attracting fixed points If there were a fixed point x *, all points in its basin of attraction are bound to have the same energy E(x * ) (since energy would be constant on all trajectories leading to x * ), so there would be an open set with constant energy What kind of fixed points can there be?
40 Example II Particle of mass m=1 moving in a doublewell potential Fixed points: (0, 0), (±1,0)
41 Example II (0, 0) > Δ = 1 <0 > saddle point! (±1, 0) > τ=0, Δ = 2 > centers! Question: will the centers survive the nonlinearity? Answer: yes, because of energy conservation!
42 Example II
43 Example II 1) Near the centers there are small periodic orbits 2) There are also large periodic orbits encircling all fixed points 3) Solutions are periodic except for equilibria (fixed points) and the homoclinic orbits, which approach the origin when t > ±
44 Example II 1) Neutrally stable equilibria correspond to the particle at rest at the bottom of either one of the wells 2) Small closed orbits > small oscillations about equilibria 3) Large closed orbits > oscillations taking the particle back and forth over the hump 4) Saddle point? Homoclinic orbits?
45 Example II Sketch the graph of the energy function 1) Local minima of E project down to centers in the phase plane 2) Contours of slightly higher energy > small closed orbits 3) Cut at Evalue of local maximum > homoclinic orbits 4) Cuts at higher Evalues > large periodic orbits
46 Nonlinear centers Theorem (nonlinear centers for conservative systems): Consider the system, where and f is continuously differentiable. Suppose there exists a conserved quantity E(x) and an isolated fixed point x *. If x * is a local minimum of E, then all trajectories sufficiently close to x * are closed. Ideas behind the proof: 1) Since E is constant on trajectories, each trajectory is contained in some contour of E. 2) Near a local maximum (or minimum), contours are closed 3) The orbit is periodic, i.e. it does not stop at some point of the contour because x * is isolated, so there are no other fixed points on its close proximity
47 Reversible systems Many mechanical systems have timereversal symmetry, i.e. if one watches the motion backwards no anomaly is seen, the reversed motion is a possible future trajectory of the system (Ex. Pendulum) Any system is symmetric under time reversal! The acceleration does not change, the velocity changes sign!
48 Reversible systems Consequence: if [x(t), y(t)] is a solution, also [x(t), y(t)] is a solution!
49 Any system Reversible systems such that f(x,y) = f(x,y) and g(x,y)=g(x,y) is reversible! Reversible systems are different from conservative systems, but they share some properties Theorem (nonlinear centers for reversible systems): Suppose the origin x * = 0 is a linear center of a reversible system. Then, sufficiently close to the origin, all orbits are closed.
50 Reversible systems Ideas behind the proof: 1) Let us take a trajectory starting on the positive xaxis near the origin 2) Because of the influence of the linear center (if the system is close enough to it), the trajectory will bend and intersect the negative xaxis
51 Reversible systems Ideas behind the proof: 3) By using reversibility we can reflect the trajectory above the xaxis, obtaining a twin trajectory (we know that it is a solution of the equation of motion, it must be the only one) 4) The two trajectories form a closed orbit, as desired
52 Example I The system is reversible and the origin (0, 0) is a fixed point, which kind? Jacobian at the origin: τ = 0, Δ = 1 > linear center > nonlinear center (theorem) Other fixed points are (1, 1) and (1, 1) Δ = 2 < 0 > saddle points
53 Example I The twin saddle points are joined by a pair of trajectories (heteroclinic orbits or saddle connections) Homoclinic and heteroclinic orbits are more common in conservative and reversible systems
54 Example II Show that there is a homoclinic orbit in the halfplane x 0 The origin is a saddle point, since the Jacobian is which has Δ = 1 < 0 The eigenvectors are v 1 = (1, 1) and v 2 = (1, 1) The unstable manifold leaves the origin following v 1 = (1, 1)
55 Example II 1) Initially we are in the first quadrant (x > 0 and y > 0) 2) Velocity in the x direction is positive, in the y direction it is positive until the system passes x = 1 3) For x > 1 the velocity in the y direction becomes negative and the particle ends up hitting the xaxis 4) By reversibility there must be a twin trajectory with the same endpoints and arrows reversed 5) The two trajectories together form a homoclinic orbit
56 Nondimensionalization: Pendulum
57 Pendulum Fixed points: (θ *, ν * ) = (kπ, 0), where k is any integer Focus: (0, 0), (π, 0) (the other fixed points coincide with either of them) (0, 0) > τ = 0, Δ = 1 > 0 > linear center! (π, 0) > τ = 0, Δ = 1 < 0 > saddle point!
58 Pendulum The system is reversible The system is conservative The energy function has a local minimum at (0, 0) The origin is then a nonlinear center
59 Pendulum The eigenvectors at the saddle fixed point are v 1 = (1, 1) and v 2 = (1, 1)
60 Pendulum Let us add the energy contours for different values of E Portrait is periodic in the θdirection
61 Pendulum Physical interpretation: 1) Center > pendulum at rest straight down (minimum energy E = 1) 2) Small orbits about the center > small oscillations (librations) 3) If the energy increases, the amplitude of the oscillations increases 4) When E = 1 the motion asymptotically reaches the position of unstable equilibrium (pendulum straight up, heteroclinic orbit) 5) When E > 1 the pendulum whirls repeatedly over the top
62 Pendulum Cylindrical phase space Natural space for pendulum: one variable (θ) is periodic, the other (ν) is not Periodic whirling motions (E > 1) look periodic Saddle points indicate the same physical state Heteroclinic trajectories become homoclinic orbits
63 Pendulum Cylindrical phase space Plotting vertically the energy instead of the velocity: Utube
64 Pendulum Cylindrical phase space Orbits are sections at constant height/energy The two arms correspond to the two senses of rotations Homoclinic orbits lie at E = 1, borderline between librations and rotations
65 All trajectories continuously lose altitude, except for the fixed points Pendulum Damping Centers > stable spirals Saddle points > saddle points
66 Pendulum Damping Change of energy along trajectory: Consequence: E decreases monotonically along trajectories, except at fixed points (where )
67 Pendulum Damping Physics: pendulum rotates over the top with decreasing energy, until it cannot complete the rotation and makes damped oscillations about equilibrium, where it eventually stops
68 Index theory Global information about the phase portrait, as opposed to the local information provided by linearization Questions: 1) Must a closed trajectory always encircle a fixed point? 2) If so, what types of fixed points are permitted? 3) What types of fixed points can coalesce in bifurcations? 4) Trajectories near higherorder fixed points? 5) Possibility of closed orbits? Index of a closed curve C: integer that measures the winding of the vector field on C Similarity with electrostatics: from the behavior of electric field on a surface one may deduce the total amount of charges inside the surface; here one gets info on possible fixed points
69 Index theory Vector field, simple (= nonselfintersecting) closed curve C, not passing through fixed points of the system At each point of C the vector field makes a welldefined angle with the positive xaxis
70 Index theory As x moves counterclockwise around C, the angle φ changes continuously (the vector field is smooth) > when x comes back to the starting position φ has varied by a multiple of 2π [φ] C = net change in φ over one circuit Index of the closed curve C:
71 Example I What s the index? The vector field makes one complete rotation counterclockwise, so I C = + 1
72 Trick The index is the number of (counterclockwise) rotations that need to be done to go from 1 to 1 following the order of the numbers
73 Example II The vector field makes one complete rotation clockwise: I C =  1
74 Example III The curve C is the unit circle: What is I C? [φ] C = π + 2π π = 0 I C = 0
75 Properties of the index 1) If C can be continuously deformed into C without passing through a fixed point, I C = I C. Proof: the index cannot vary continuously, but only by integer values, so it cannot be altered by a continuous change of C 2) If C does not enclose any fixed points, I C = 0. Proof: by squeezing C until it becomes a very small circle the index does not change because of 1), and it equals zero because all vectors on the tiny circle are essentially identical (as the field varies smoothly in space)
76 Properties of the index 3) Under time reversal (t > t), the index is the same. Proof: the time reversal changes the signs of the velocity vectors, so the angles change from φ to φ + π, hence [φ] C stays the same 4) If C is a trajectory of the system, I C = + 1
77 Index of a point The index of an isolated fixed point x * is the index of the vector field on any closed curve encircling x * and no other fixed point By property 1), the value of the index is the same on any curve C, since it can be continuously deformed onto any other What is the index of a stable node? The vector field makes one complete rotation counterclockwise, so I = + 1 The value is the same for unstable nodes as well, as the situation would be the same, only with reversed arrows (property 3)
78 Index of a point What is the index of a saddle point? The vector field makes one complete rotation clockwise: I =  1 Spiral, centers, degenerate nodes and stars all have I = + 1, only saddle points have a different value
79 Index of a point Theorem: If a closed curve C surrounds n isolated fixed points, the index of the vector field on C equals the sum of the indices of the enclosed fixed points Proof: C can be deformed to the contour Γ of the figure; the contributions of the bridges cancel as each bridge is crossed in both directions so the net changes in the angle are equal and opposite
80 Index of a point Theorem: Any closed orbit in the phase plane must enclose fixed points whose indices sum to +1 Proof: If C is a closed orbit, I C = +1. From the previous theorem this is also the sum of the indices of the fixed points enclosed by C Consequence: any closed orbit encloses at least one fixed point (if there were none, the index on the curve would be 0, instead of + 1). If there is a unique fixed point, it cannot be a saddle (as in this case the index would be 1)
81 Example I Show that closed orbits are impossible for the rabbits versus sheep system The only orbits enclosing good fixed points (i.e. with index + 1) would have to cross the x/y axes, which contain trajectories of the system, and trajectories cannot cross!
82 Example II Show that the system has no closed orbits Solution: the system has no fixed points, so it cannot have closed orbits, since the latter have to enclose at least one fixed points
83 Exercises IV 1) Find and classify all fixed points, and sketch the nullclines and the phase portrait for 2) Find the fixed points, classify them, sketch the neighboring trajectories and try to fill in the phase portrait, for the systems
Physics: springmass system, planet motion, pendulum. Biology: ecology problem, neural conduction, epidemics
Applications of nonlinear ODE systems: Physics: springmass system, planet motion, pendulum Chemistry: mixing problems, chemical reactions Biology: ecology problem, neural conduction, epidemics Economy:
More informationIs chaos possible in 1d?  yes  no  I don t know. What is the long term behavior for the following system if x(0) = π/2?
Is chaos possible in 1d?  yes  no  I don t know What is the long term behavior for the following system if x(0) = π/2? In the insect outbreak problem, what kind of bifurcation occurs at fixed value
More informationProblem set 7 Math 207A, Fall 2011 Solutions
Problem set 7 Math 207A, Fall 2011 s 1. Classify the equilibrium (x, y) = (0, 0) of the system x t = x, y t = y + x 2. Is the equilibrium hyperbolic? Find an equation for the trajectories in (x, y) phase
More information7 Twodimensional bifurcations
7 Twodimensional bifurcations As in onedimensional systems: fixed points may be created, destroyed, or change stability as parameters are varied (change of topological equivalence ). In addition closed
More informationNonlinear Autonomous Systems of Differential
Chapter 4 Nonlinear Autonomous Systems of Differential Equations 4.0 The Phase Plane: Linear Systems 4.0.1 Introduction Consider a system of the form x = A(x), (4.0.1) where A is independent of t. Such
More information1. < 0: the eigenvalues are real and have opposite signs; the fixed point is a saddle point
Solving a Linear System τ = trace(a) = a + d = λ 1 + λ 2 λ 1,2 = τ± = det(a) = ad bc = λ 1 λ 2 Classification of Fixed Points τ 2 4 1. < 0: the eigenvalues are real and have opposite signs; the fixed point
More information1 The pendulum equation
Math 270 Honors ODE I Fall, 2008 Class notes # 5 A longer than usual homework assignment is at the end. The pendulum equation We now come to a particularly important example, the equation for an oscillating
More informationA plane autonomous system is a pair of simultaneous firstorder differential equations,
Chapter 11 PhasePlane Techniques 11.1 Plane Autonomous Systems A plane autonomous system is a pair of simultaneous firstorder differential equations, ẋ = f(x, y), ẏ = g(x, y). This system has an equilibrium
More informationMATH 215/255 Solutions to Additional Practice Problems April dy dt
. For the nonlinear system MATH 5/55 Solutions to Additional Practice Problems April 08 dx dt = x( x y, dy dt = y(.5 y x, x 0, y 0, (a Show that if x(0 > 0 and y(0 = 0, then the solution (x(t, y(t of the
More informationChapter 6 Nonlinear Systems and Phenomena. Friday, November 2, 12
Chapter 6 Nonlinear Systems and Phenomena 6.1 Stability and the Phase Plane We now move to nonlinear systems Begin with the firstorder system for x(t) d dt x = f(x,t), x(0) = x 0 In particular, consider
More informationStability of Dynamical systems
Stability of Dynamical systems Stability Isolated equilibria Classification of Isolated Equilibria Attractor and Repeller Almost linear systems Jacobian Matrix Stability Consider an autonomous system u
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 informationLECTURE 8: DYNAMICAL SYSTEMS 7
15382 COLLECTIVE INTELLIGENCE S18 LECTURE 8: DYNAMICAL SYSTEMS 7 INSTRUCTOR: GIANNI A. DI CARO GEOMETRIES IN THE PHASE SPACE Damped pendulum One cp in the region between two separatrix Separatrix Basin
More informationMath 232, Final Test, 20 March 2007
Math 232, Final Test, 20 March 2007 Name: Instructions. Do any five of the first six questions, and any five of the last six questions. Please do your best, and show all appropriate details in your solutions.
More information8.1 Bifurcations of Equilibria
1 81 Bifurcations of Equilibria Bifurcation theory studies qualitative changes in solutions as a parameter varies In general one could study the bifurcation theory of ODEs PDEs integrodifferential equations
More informationFundamentals of Dynamical Systems / DiscreteTime Models. Dr. Dylan McNamara people.uncw.edu/ mcnamarad
Fundamentals of Dynamical Systems / DiscreteTime Models Dr. Dylan McNamara people.uncw.edu/ mcnamarad Dynamical systems theory Considers how systems autonomously change along time Ranges from Newtonian
More information6.3. Nonlinear Systems of Equations
G. NAGY ODE November,.. Nonlinear Systems of Equations Section Objective(s): Part One: TwoDimensional Nonlinear Systems. ritical Points and Linearization. The HartmanGrobman Theorem. Part Two: ompeting
More informationPractice Problems for Final Exam
Math 1280 Spring 2016 Practice Problems for Final Exam Part 2 (Sections 6.6, 6.7, 6.8, and chapter 7) S o l u t i o n s 1. Show that the given system has a nonlinear center at the origin. ẋ = 9y 5y 5,
More informationOne Dimensional Dynamical Systems
16 CHAPTER 2 One Dimensional Dynamical Systems We begin by analyzing some dynamical systems with onedimensional phase spaces, and in particular their bifurcations. All equations in this Chapter are scalar
More informationOutline. Learning Objectives. References. Lecture 2: Secondorder Systems
Outline Lecture 2: Secondorder Systems! Techniques based on linear systems analysis! Phaseplane analysis! Example: Neanderthal / Early man competition! HartmanGrobman theorem  validity of linearizations!
More informationEE222  Spring 16  Lecture 2 Notes 1
EE222  Spring 16  Lecture 2 Notes 1 Murat Arcak January 21 2016 1 Licensed under a Creative Commons AttributionNonCommercialShareAlike 4.0 International License. Essentially Nonlinear Phenomena Continued
More informationMATH 415, WEEKS 7 & 8: Conservative and Hamiltonian Systems, Nonlinear Pendulum
MATH 415, WEEKS 7 & 8: Conservative and Hamiltonian Systems, Nonlinear Pendulum Reconsider the following example from last week: dx dt = x y dy dt = x2 y. We were able to determine many qualitative features
More information154 Chapter 9 Hints, Answers, and Solutions The particular trajectories are highlighted in the phase portraits below.
54 Chapter 9 Hints, Answers, and Solutions 9. The Phase Plane 9.. 4. The particular trajectories are highlighted in the phase portraits below... 3. 4. 9..5. Shown below is one possibility with x(t) and
More informationModels Involving Interactions between Predator and Prey Populations
Models Involving Interactions between Predator and Prey Populations Matthew Mitchell Georgia College and State University December 30, 2015 Abstract Predatorprey models are used to show the intricate
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 informationComplex Dynamic Systems: Qualitative vs Quantitative analysis
Complex Dynamic Systems: Qualitative vs Quantitative analysis Complex Dynamic Systems Chiara Mocenni Department of Information Engineering and Mathematics University of Siena (mocenni@diism.unisi.it) Dynamic
More information2DVolterraLotka Modeling For 2 Species
Majalat AlUlum AlInsaniya wat  Tatbiqiya 2DVolterraLotka Modeling For 2 Species Alhashmi Darah 1 University of Almergeb Department of Mathematics Faculty of Science Zliten Libya. Abstract The purpose
More information= F ( x; µ) (1) where x is a 2dimensional vector, µ is a parameter, and F :
1 Bifurcations Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University Tallahassee, Florida 32306 A bifurcation is a qualitative change
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 informationTWO DIMENSIONAL FLOWS. Lecture 5: Limit Cycles and Bifurcations
TWO DIMENSIONAL FLOWS Lecture 5: Limit Cycles and Bifurcations 5. Limit cycles A limit cycle is an isolated closed trajectory [ isolated means that neighbouring trajectories are not closed] Fig. 5.1.1
More informationMathematical Modeling I
Mathematical Modeling I Dr. Zachariah Sinkala Department of Mathematical Sciences Middle Tennessee State University Murfreesboro Tennessee 37132, USA November 5, 2011 1d systems To understand more complex
More informationMath 312 Lecture Notes Linear Twodimensional Systems of Differential Equations
Math 2 Lecture Notes Linear Twodimensional Systems of Differential Equations Warren Weckesser Department of Mathematics Colgate University February 2005 In these notes, we consider the linear system of
More informationNonlinear Autonomous Dynamical systems of two dimensions. Part A
Nonlinear Autonomous Dynamical systems of two dimensions Part A Nonlinear Autonomous Dynamical systems of two dimensions x f ( x, y), x(0) x vector field y g( xy, ), y(0) y F ( f, g) 0 0 f, g are continuous
More informationDynamical Systems and Chaos Part I: Theoretical Techniques. Lecture 4: Discrete systems + Chaos. Ilya Potapov Mathematics Department, TUT Room TD325
Dynamical Systems and Chaos Part I: Theoretical Techniques Lecture 4: Discrete systems + Chaos Ilya Potapov Mathematics Department, TUT Room TD325 Discrete maps x n+1 = f(x n ) Discrete time steps. x 0
More informationChapter 1, Section 1.2, Example 9 (page 13) and Exercise 29 (page 15). Use the Uniqueness Tool. Select the option ẋ = x
Use of Tools from Interactive Differential Equations with the texts Fundamentals of Differential Equations, 5th edition and Fundamentals of Differential Equations and Boundary Value Problems, 3rd edition
More informationarxiv: v1 [physics.classph] 5 Jan 2012
Damped bead on a rotating circular hoop  a bifurcation zoo Shovan Dutta Department of Electronics and Telecommunication Engineering, Jadavpur University, Calcutta 700 032, India. Subhankar Ray Department
More informationPhase portraits in two dimensions
Phase portraits in two dimensions 8.3, Spring, 999 It [ is convenient to represent the solutions to an autonomous system x = f( x) (where x x = ) by means of a phase portrait. The x, y plane is called
More informationNonLinear Dynamics Homework Solutions Week 6
NonLinear Dynamics Homework Solutions Week 6 Chris Small March 6, 2007 Please email me at smachr09@evergreen.edu with any questions or concerns reguarding these solutions. 6.8.3 Locate annd find the index
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 informationANSWERS Final Exam Math 250b, Section 2 (Professor J. M. Cushing), 15 May 2008 PART 1
ANSWERS Final Exam Math 50b, Section (Professor J. M. Cushing), 5 May 008 PART. (0 points) A bacterial population x grows exponentially according to the equation x 0 = rx, where r>0is the per unit rate
More information7.7 LOTKAVOLTERRA M ODELS
77 LOTKAVOLTERRA M ODELS sstems, such as the undamped pendulum, enabled us to draw the phase plane for this sstem and view it globall In that case we were able to understand the motions for all initial
More informationLinear Planar Systems Math 246, Spring 2009, Professor David Levermore We now consider linear systems of the form
Linear Planar Systems Math 246, Spring 2009, Professor David Levermore We now consider linear systems of the form d x x 1 = A, where A = dt y y a11 a 12 a 21 a 22 Here the entries of the coefficient matrix
More information4 SecondOrder Systems
4 SecondOrder Systems Secondorder autonomous systems occupy an important place in the study of nonlinear systems because solution trajectories can be represented in the plane. This allows for easy visualization
More informationMath 266: Phase Plane Portrait
Math 266: Phase Plane Portrait Long Jin Purdue, Spring 2018 Review: Phase line for an autonomous equation For a single autonomous equation y = f (y) we used a phase line to illustrate the equilibrium solutions
More informationCHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD
CHALMERS, GÖTEBORGS UNIVERSITET EXAM for DYNAMICAL SYSTEMS COURSE CODES: TIF 155, FIM770GU, PhD Time: Place: Teachers: Allowed material: Not allowed: January 08, 2018, at 08 30 12 30 Johanneberg Kristian
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 informationMATH 415, WEEK 11: Bifurcations in Multiple Dimensions, Hopf Bifurcation
MATH 415, WEEK 11: Bifurcations in Multiple Dimensions, Hopf Bifurcation 1 Bifurcations in Multiple Dimensions When we were considering onedimensional systems, we saw that subtle changes in parameter
More information11 Chaos in Continuous Dynamical Systems.
11 CHAOS IN CONTINUOUS DYNAMICAL SYSTEMS. 47 11 Chaos in Continuous Dynamical Systems. Let s consider a system of differential equations given by where x(t) : R R and f : R R. ẋ = f(x), The linearization
More informationCopyright (c) 2006 Warren Weckesser
2.2. PLANAR LINEAR SYSTEMS 3 2.2. Planar Linear Systems We consider the linear system of two first order differential equations or equivalently, = ax + by (2.7) dy = cx + dy [ d x x = A x, where x =, and
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 informationChapter 7. Nonlinear Systems. 7.1 Introduction
Nonlinear Systems Chapter 7 The scientist does not study nature because it is useful; he studies it because he delights in it, and he delights in it because it is beautiful.  Jules Henri Poincaré (18541912)
More informationPart II. Dynamical Systems. Year
Part II Year 2017 2016 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2017 34 Paper 1, Section II 30A Consider the dynamical system where β > 1 is a constant. ẋ = x + x 3 + βxy 2, ẏ = y + βx 2
More informationComputational Neuroscience. Session 42
Computational Neuroscience. Session 42 Dr. Marco A Roque Sol 06/21/2018 TwoDimensional TwoDimensional System In this section we will introduce methods of phase plane analysis of twodimensional systems.
More informationMATH 415, WEEKS 14 & 15: 1 Recurrence Relations / Difference Equations
MATH 415, WEEKS 14 & 15: Recurrence Relations / Difference Equations 1 Recurrence Relations / Difference Equations In many applications, the systems are updated in discrete jumps rather than continuous
More informationREUNotes08ODEs May 30, Chapter Two. Differential Equations
Chapter Two Differential Equations 4 CHAPTER 2 2.1 LINEAR EQUATIONS We consider ordinary differential equations of the following form dx = f(t, x). dt where x is either a real variable or a vector of several
More informationPhase Portraits of Nonlinear Differential Equations
ODE4net.nb Phase Portraits of Nonlinear Differential Equations Nonlinear Differential Equations: x' = f (x, y) (1) Consider the system where f and g are functions of two y' = g(x, y) () variables x and
More informationDynamical Systems in Neuroscience: Elementary Bifurcations
Dynamical Systems in Neuroscience: Elementary Bifurcations Foris Kuang May 2017 1 Contents 1 Introduction 3 2 Definitions 3 3 HodgkinHuxley Model 3 4 MorrisLecar Model 4 5 Stability 5 5.1 Linear ODE..............................................
More informationLotka Volterra PredatorPrey Model with a Predating Scavenger
Lotka Volterra PredatorPrey Model with a Predating Scavenger Monica Pescitelli Georgia College December 13, 2013 Abstract The classic Lotka Volterra equations are used to model the population dynamics
More informationLab 5: Nonlinear Systems
Lab 5: Nonlinear Systems Goals In this lab you will use the pplane6 program to study two nonlinear systems by direct numerical simulation. The first model, from population biology, displays interesting
More informationFIRSTORDER SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS III: Autonomous Planar Systems David Levermore Department of Mathematics University of Maryland
FIRSTORDER SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS III: Autonomous Planar Systems David Levermore Department of Mathematics University of Maryland 4 May 2012 Because the presentation of this material
More informationAPPPHYS217 Tuesday 25 May 2010
APPPHYS7 Tuesday 5 May Our aim today is to take a brief tour of some topics in nonlinear dynamics. Some good references include: [Perko] Lawrence Perko Differential Equations and Dynamical Systems (SpringerVerlag
More informationMath 273 (51)  Final
Name: Id #: Math 273 (5)  Final Autumn Quarter 26 Thursday, December 8, 266: to 8: Instructions: Prob. Points Score possible 25 2 25 3 25 TOTAL 75 Read each problem carefully. Write legibly. Show all
More informationB5.6 Nonlinear Systems
B5.6 Nonlinear Systems 4. Bifurcations Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Local bifurcations for vector fields 1.1 The problem 1.2 The extended centre
More informationChaotic motion. Phys 750 Lecture 9
Chaotic motion Phys 750 Lecture 9 Finitedifference equations Finite difference equation approximates a differential equation as an iterative map (x n+1,v n+1 )=M[(x n,v n )] Evolution from time t =0to
More informationStability Analysis for ODEs
Stability Analysis for ODEs Marc R Roussel September 13, 2005 1 Linear stability analysis Equilibria are not always stable Since stable and unstable equilibria play quite different roles in the dynamics
More informationDynamics of a Population Model Controlling the Spread of Plague in Prairie Dogs
Dynamics of a opulation Model Controlling the Spread of lague in rairie Dogs Catalin Georgescu The University of South Dakota Department of Mathematical Sciences 414 East Clark Street, Vermillion, SD USA
More informationCHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD
CHALMERS, GÖTEBORGS UNIVERSITET EXAM for DYNAMICAL SYSTEMS COURSE CODES: TIF 155, FIM770GU, PhD Time: Place: Teachers: Allowed material: Not allowed: January 14, 2019, at 08 30 12 30 Johanneberg Kristian
More informationCHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD
CHALMERS, GÖTEBORGS UNIVERSITET EXAM for DYNAMICAL SYSTEMS COURSE CODES: TIF 155, FIM770GU, PhD Time: Place: Teachers: Allowed material: Not allowed: April 06, 2018, at 14 00 18 00 Johanneberg Kristian
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 informationBasic Theory of Dynamical Systems
1 Basic Theory of Dynamical Systems Page 1 1.1 Introduction and Basic Examples Dynamical systems is concerned with both quantitative and qualitative properties of evolution equations, which are often ordinary
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 information3.5 Competition Models: Principle of Competitive Exclusion
94 3. Models for Interacting Populations different dimensional parameter changes. For example, doubling the carrying capacity K is exactly equivalent to halving the predator response parameter D. The dimensionless
More informationLogistic Map, Euler & RungeKutta Method and LotkaVolterra Equations
Logistic Map, Euler & RungeKutta Method and LotkaVolterra 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 informationA ProducerConsumer Model With Stoichiometry
A ProducerConsumer Model With Stoichiometry Plan B project toward the completion of the Master of Science degree in Mathematics at University of Minnesota Duluth Respectfully submitted by Laura Joan Zimmermann
More informationIntroduction to Dynamical Systems
Introduction to Dynamical Systems Autonomous Planar Systems Vector form of a Dynamical System Trajectories Trajectories Don t Cross Equilibria Population Biology RabbitFox System Trout System Trout System
More informationSection 9.3 Phase Plane Portraits (for Planar Systems)
Section 9.3 Phase Plane Portraits (for Planar Systems) Key Terms: Equilibrium point of planer system yꞌ = Ay o Equilibrium solution Exponential solutions o Halfline solutions Unstable solution Stable
More informationVector Field Topology. Ronald Peikert SciVis Vector Field Topology 81
Vector Field Topology Ronald Peikert SciVis 2007  Vector Field Topology 81 Vector fields as ODEs What are conditions for existence and uniqueness of streamlines? For the initial value problem i x ( t)
More informationNonlinear Control Lecture 2:Phase Plane Analysis
Nonlinear Control Lecture 2:Phase Plane Analysis Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Fall 2010 r. Farzaneh Abdollahi Nonlinear Control Lecture 2 1/53
More informationSolutions of a PTsymmetric Dimer with Constant Gainloss
Solutions of a PTsymmetric Dimer with Constant Gainloss G14DIS Mathematics 4th Year Dissertation Spring 2012/2013 School of Mathematical Sciences University of Nottingham John Pickton Supervisor: Dr
More informationCHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD
CHALMERS, GÖTEBORGS UNIVERSITET EXAM for DYNAMICAL SYSTEMS COURSE CODES: TIF 155, FIM770GU, PhD Time: Place: Teachers: Allowed material: Not allowed: August 22, 2018, at 08 30 12 30 Johanneberg Jan Meibohm,
More informationChapter 9 Global Nonlinear Techniques
Chapter 9 Global Nonlinear Techniques Consider nonlinear dynamical system 0 Nullcline X 0 = F (X) = B @ f 1 (X) f 2 (X). f n (X) x j nullcline = fx : f j (X) = 0g equilibrium solutions = intersection of
More informationLecture Notes for Math 251: ODE and PDE. Lecture 27: 7.8 Repeated Eigenvalues
Lecture Notes for Math 25: ODE and PDE. Lecture 27: 7.8 Repeated Eigenvalues Shawn D. Ryan Spring 22 Repeated Eigenvalues Last Time: We studied phase portraits and systems of differential equations with
More informationMathQuest: Differential Equations
MathQuest: Differential Equations Geometry of Systems 1. The differential equation d Y dt = A Y has two straight line solutions corresponding to [ ] [ ] 1 1 eigenvectors v 1 = and v 2 2 = that are shown
More informationAn Undergraduate s Guide to the HartmanGrobman and PoincaréBendixon Theorems
An Undergraduate s Guide to the HartmanGrobman and PoincaréBendixon Theorems Scott Zimmerman MATH181HM: Dynamical Systems Spring 2008 1 Introduction The HartmanGrobman and PoincaréBendixon Theorems
More informationNonlinear Dynamics. Moreno Marzolla Dip. di Informatica Scienza e Ingegneria (DISI) Università di Bologna.
Nonlinear Dynamics Moreno Marzolla Dip. di Informatica Scienza e Ingegneria (DISI) Università di Bologna http://www.moreno.marzolla.name/ 2 Introduction: Dynamics of Simple Maps 3 Dynamical systems A dynamical
More informationMATH 614 Dynamical Systems and Chaos Lecture 24: Bifurcation theory in higher dimensions. The Hopf bifurcation.
MATH 614 Dynamical Systems and Chaos Lecture 24: Bifurcation theory in higher dimensions. The Hopf bifurcation. Bifurcation theory The object of bifurcation theory is to study changes that maps undergo
More informationLecture 3. Dynamical Systems in Continuous Time
Lecture 3. Dynamical Systems in Continuous Time University of British Columbia, Vancouver YueXian Li November 2, 2017 1 3.1 Exponential growth and decay A Population With Generation Overlap Consider a
More informationThe Liapunov Method for Determining Stability (DRAFT)
44 The Liapunov Method for Determining Stability (DRAFT) 44.1 The Liapunov Method, Naively Developed In the last chapter, we discussed describing trajectories of a 2 2 autonomous system x = F(x) as level
More informationFinal Exam Review. Review of Systems of ODE. Differential Equations Lia Vas. 1. Find all the equilibrium points of the following systems.
Differential Equations Lia Vas Review of Systems of ODE Final Exam Review 1. Find all the equilibrium points of the following systems. (a) dx = x x xy (b) dx = x x xy = 0.5y y 0.5xy = 0.5y 0.5y 0.5xy.
More informationSolutions for B8b (Nonlinear Systems) Fake Past Exam (TT 10)
Solutions for B8b (Nonlinear Systems) Fake Past Exam (TT 10) Mason A. Porter 15/05/2010 1 Question 1 i. (6 points) Define a saddlenode bifurcation and show that the first order system dx dt = r x e x
More informationGerardo Zavala. Math 388. PredatorPrey Models
Gerardo Zavala Math 388 PredatorPrey Models Spring 2013 1 History In the 1920s A. J. Lotka developed a mathematical model for the interaction between two species. The mathematician Vito Volterra worked
More informationLecture 20/Lab 21: Systems of Nonlinear ODEs
Lecture 20/Lab 21: Systems of Nonlinear ODEs MAR514 Geoffrey Cowles Department of Fisheries Oceanography School for Marine Science and Technology University of MassachusettsDartmouth Coupled ODEs: Species
More informationQualitative Analysis of TumorImmune ODE System
of TumorImmune ODE System LG de Pillis and AE Radunskaya August 15, 2002 This work was supported in part by a grant from the WM Keck Foundation 00 QUALITATIVE ANALYSIS Overview 1 Simplified System of
More informationM.5 Modeling the Effect of Functional Responses
M.5 Modeling the Effect of Functional Responses The functional response is referred to the predation rate as a function of the number of prey per predator. It is recognized that as the number of prey increases,
More information+ i. cos(t) + 2 sin(t) + c 2.
MATH HOMEWORK #7 PART A SOLUTIONS Problem 7.6.. Consider the system x = 5 x. a Express the general solution of the given system of equations in terms of realvalued functions. b Draw a direction field,
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 informationSolutions to Final Exam Sample Problems, Math 246, Spring 2011
Solutions to Final Exam Sample Problems, Math 246, Spring 2 () Consider the differential equation dy dt = (9 y2 )y 2 (a) Identify its equilibrium (stationary) points and classify their stability (b) Sketch
More informationStable and Unstable Manifolds for Planar Dynamical Systems
Stable and Unstable Manifolds for Planar Dynamical Systems H.L. Smith Department of Mathematics Arizona State University Tempe, AZ 85287 184 October 7, 211 Key words and phrases: Saddle point, stable and
More informationPHY411 Lecture notes Part 5
PHY411 Lecture notes Part 5 Alice Quillen January 27, 2016 Contents 0.1 Introduction.................................... 1 1 Symbolic Dynamics 2 1.1 The Shift map.................................. 3 1.2
More informationChaotic motion. Phys 420/580 Lecture 10
Chaotic motion Phys 420/580 Lecture 10 Finitedifference equations Finite difference equation approximates a differential equation as an iterative map (x n+1,v n+1 )=M[(x n,v n )] Evolution from time t
More information