# B5.6 Nonlinear Systems

Size: px
Start display at page:

## Transcription

1 B5.6 Nonlinear Systems 4. Bifurcations Alain Goriely 2018 Mathematical Institute, University of Oxford

2 Table of contents 1. Local bifurcations for vector fields 1.1 The problem 1.2 The extended centre manifold 1.3 The Hopf bifurcation 1.4 Sotomayor s theorem 1.5 Local bifurcation of maps 1

3 What we learned from Section 1&2. For a system of linear autonomous equations ẋ = Ax, the solutions live on invariant spaces that can be classified according to the eigenvalues of A. The stable (resp. unstable, centre) linear subspace is the span of eigenvectors whose eigenvalues have a negative (resp. positive, null) real part. For nonlinear systems, we define the notion of asymptotic sets (α and ω limit set), the notion of attracting set, and basin of attraction. We define two important notions of stability for a fixed point: (Lyapunov) stability (i.e. solutions remain close ) and exponential stability i.e. ( fixed point is stable AND all nearby solutions converge to the fixed point asymptotically for long time ). Lyapunov functions can be used to test stability. But, finding a Lyapunov function can be difficult. 2

4 What we learned from Section 3. For a system of autonomous equations ẋ = f(x), we are interested in the trajectories and asymptotic sets in phase space. At a fixed point, we can define local stable, unstable, and centre manifolds based on the corresponding linear subspaces of the linearised system. From the local stable and unstable manifolds, we can define the global stable and unstable manifolds by extending them to all times. If the unstable manifold is non-empty, the fixed point is unstable. If the unstable and centre manifolds are empty, the fixed point is asymptotically stable. If the unstable manifold is empty but the centre manifold is non-empty, we can study the dynamics on the centre manifold by centre manifold reduction. The same notions can be defined for iterative mappings and for periodic orbits. 3

5 4. Local bifurcations for vector fields

6 4. Local bifurcations for vector fields 4.1 The problem

7 The problem Consider the nonlinear, autonomous, first-order system of differential equations ẋ = f(x, µ) where x : E R n, µ : U R p (1) where µ is a vector of parameters. Problem: How does the dynamics change when the parameters are varied? What are the special values where qualitative changes occur? Can the possible changes be classified? Can they be obtained algorithmically? 4

8 An important example: the transcritical bifurcation Take p = n = 1 and consider the system ẋ = µx x 2 (2) 5

9 Main idea We consider the problem of determining bifurcations at fixed points. Take µ R (p=1). Then the fixed points are given by the solution of f(x, µ) = 0. (3) For values of µ for which a solution can be found, the solution defines a branch of equilibria x = x(µ). Along this branch, define the matrix D(µ) = D x f (x(µ),µ). (4) 6

10 Main idea D(µ) = D x f (x(µ),µ). (5) Assume that there is a value µ 0 for which D(µ) has only eigenvalues with non-zero real parts (i.e. the fixed point is hyperbolic for that value). Then, D(µ) is invertible and the local branch of equilibria can be continued locally. We increase (or decrease) µ up to a critical bifurcation value µ c where D(µ) is not invertible. At this point, the branch of equilibria is non-smooth. 7

11 Definitions For a system ẋ = f(x, µ) where x : E R n, µ : U R p the value µ c is a bifurcation value, if there is a qualitative change of behaviour in any neighborhood V of µ c (and for any small enough pertubations of the vector field). Note 1: The precise definition of bifurcation is related to the the notion of structural stability that implies that small enough perturbations of the vector field do not change locally the phase portrait. That is, there exists a homeomorphism mapping the trajectories of the original system to the trajectories of the perturbed system. A bifurcation occurs when the system fails to be structurally stable. Note 2: Here, as a starting point, we use the working definition for bifurcation at a fixed point: a bifurcation occurs at values where the number or stability of a fixed point change (for small enough perturbations of the vector field or parameters). 8

12 More definitions For a system ẋ = f(x, µ) where x : E R n, µ : U R p the bifurcation set is the loci in µ-space of bifurcation values. A bifurcation diagram is the loci in (x, µ)-space where f(x, µ) = 0 with at least one bifurcation value. 9

13 Another important example: the saddle-node bifurcation Consider the system ẋ = µ x 2 (6) 10

14 Yet another important example: the pitchfork bifurcation Consider the system ẋ = µx x 3 (7) 11

15 4. Local bifurcations for vector fields 4.2 The extended centre manifold

16 Bifurcation through reduction We are interested in problems for which a single parameter is responsible for the bifurcation (co-dimension 1 bifurcations). Therefore, we take p = 1 and we have a system of the form ż = F(z, λ), z R n, λ R (8) Assume that for λ = λ c, there is a fixed point z c and the matrix DF(z c ) has eigenvalues with zero real parts. We use the change of variables { z = zc + C z, µ = λ λ c (9) where C is chosen such that C 1 MC = [ A 0 0 B ]. (10) 12

17 Bifurcation through reduction The matrices A and B of respective dimension n c and n s + n u are s.t.: { Re(λ) = 0 λ Spec(A), Re(λ) 0 λ Spec(B). (11) After the change of variables, the new system in the variable z = (x, y) is {ẋ = Ax + f(x, y, µ), x R n c ẏ = By + g(x, y, µ) y R ns +nu (12) and (x, y) = (0, 0) is a fixed point for µ = 0. 13

18 Bifurcation through reduction MAIN IDEA: We extend the centre manifold to include the parameter. ẋ = Ax + f(x, y, µ), x R nc ns +nu ẏ = By + g(x, y, µ) y R (13) µ = 0 We can view this system as a dynamical system in the extended phase space in m = n s + n u + n c + 1 dimensions. Note 1: the vector y now denotes variables that are both in the stable and unstable manifolds. Note 2: The important realisation is that the centre manifold has now dimension n c + 1 and is parameterised by the vector (x, µ). 14

19 Bifurcation through reduction We can proceed as before and look for a center manifold of the form y = h(x, µ), (14) and obtain the extended centre manifold. Once this is known, we can write the dynamics on the extended centre manifold as: ẋ = Ax + f(x, h(x, µ), µ). (15) This equation captures the relevant part of the bifurcation. 15

20 Example Consider the system {ẋ = µ(x + y) (x + y) 2 ẏ = y µ(x + y) + (x + y) 2 (16) 16

21 (c). Dynamics on the centre manifold. Locally on the extended Centre Manifold µ =0 is trivial so it is the ẋ equation that is interesting: Example ẋ = µ(x + x 2 µx +...) (x 2 +2x(x 2 µx)+...) = µx x 2 + O(3) Substituting back into the equation for ẋ we get (to leading order) ẋ = µx x }{{} 2 +O(x 3 ). On the center manifold, the Standard dynamics Form for transcritical and bifurcation diagram are: The phase portrait for the reduced dynamics for x is shown in Figure 5.6 and the phase portrait for the original system is in Figure 5.7. µ<0 x µ<0 µ µ<0 Figure 5.6: Phase portraits on the (one-dimensional) centre manifoldandthebifurcation diagram. Remark. If the stable manifold is of higher dimension, then y 1 = h 1 (x, µ),y 2 = h 2 (x, µ) and we need to find h 1,h 2 using the same method. For example, for the system ẋ = µx yz, ẏ = y + x 2, ż = z + x 3. 17

22 Example In the full space (x, y) the dynamics for µ < 0 (left) and µ > 0 (right) is: MATH64041/44041: Applied Dynamical Systems y Figure 5.7: Full phase portraits of the dynamics in µ<0andµ>0. x Then a 1 =1,a 2 = a 3 = b 1 = b 2 = b 3 =0. Thatisy = x 2 +,butwehavetogoto cubic polynomials to find the stable manifold for z, whichgivesz = x 3 +.Therefore,the reduced dynamics on the stable manifold is ẋ = f µ (x) =µx x 5. 18

23 Generic bifurcations with one zero eigenvalue Note: If n c = 1, the equation on the center manifold is a single equation and the typical bifurcations are: 1. The saddle-node bifurcation (generic) ẋ = µ ax 2. (17) 2. The transcritical bifurcation (generic, assuming that there is always at least an equilibrium solution) ẋ = µx ax 2. (18) 3. The pitchfork bifurcation (generic, assuming that there is a symmetry in the system that enforces the existence of pair of fixed points after bifurcation) ẋ = µx ax 3. (19) 19

24 giving two new solutions in whichever sign of µ makes the right hand side positive. There are no other ways of balancing leading order terms (by posing x µ α )sothesearetheonly bifurcating solutions. Since Subcritical vs supercitical Consider the two pitchfork bifurcations x f(x, µ) =f xµµ f xxxx 2 +, (5.10) we see that the solution x ẋ = µx x 3 (20) fµµ 2f xµ µ is stable (locally) if f xµ µ<0andunstableiff xµ µ>0. So the sign of f ẋ = µx + x 3 xµ determines on which side of µ =0thisbranchisstable. (21) The stability of second set of solutions is determined by substituting (5.9) into (5.10) giving In the 2f first xµ µ case, and sothe thenon-trivial stability is the branch opposite is stable of the simple (supercritical branch described bifurcation) above. InThis the is second called case, a pitchfork the non-trivial bifurcation: branch if the non-trivial is unstable branch (subcritical is stablebif.) it is called a supercritical pitchfork bifurcation and if the non-trivial branch is unstable it is called a subcritical pitchfork bifurcation, as shown in Figure 5.8. Supercritical Pitchfork Bifurcation ẋ = µx x 3 Subcritical Pitchfork Bifurcation ẋ = µx + x 3 Figure 5.8: Two types of Pitchfork Bifurcation 20

25 4. Local bifurcations for vector fields 4.3 The Hopf bifurcation

26 The Hopf bifurcation There is another generic bifurcation with one parameter. It happens when the eigenvalues at the bifurcation are imaginary. Recall that we can bring a system to its canonical form ẋ = Ax + f(x, y, µ), x R nc ns +nu ẏ = By + g(x, y, µ) y R µ = 0 (22) We have studied the case where A is of dimension 1 and vanishes at the bifurcation. Next, we study the case where [ ] 0 ω A =. (23) ω 0 21

27 The Hopf bifurcation On the center manifold, at the bifurcation, the dynamics of the linear part (with x = (x, y)) is {ẋ = ωy ẏ = ωx. (24) To obtain the behaviour of the system close to the bifurcation (unfolding), we consider the generic perturbation around the linear system: Close to the bifurcation, the system unfolds to {ẋ = µx ωy + f (x, y, µ) ẏ = ωx + µy + g(x, y, µ). (25) 22

28 Example Consider the typical example of a Hopf bifurcation {ẋ = µx ωy x(x 2 + y 2 ) ẏ = ωx + µy y(x 2 + y 2 ). (26) In polar coordinates, it reads {ṙ = µr r 3 θ = ω. (27) 23

29 Example 24

30 Example 25

31 The Hopf bifurcation In polar coordinates, the general form of a Hopf bifurcation is {ṙ = dµr + ar 3 θ = ω + cµ + br 2, (28) where a, b, c, d are parameters that depend on the vector field at the bifurcation. The parameters c and d can be found from a linear analysis. Let λ(µ) be the eigenvalue such that λ(0) = iω, then { d = d dµ Re(λ(µ)), c = d dµ Im(λ(µ)). (29) A step-by-step method to compute the coefficients of a Hopf normal form for a planar system is given in the file hopfalgebra.pdf. 26

32 4. Local bifurcations for vector fields 4.4 Sotomayor s theorem

33 Sotomayor s theorem (1973) Consider the system ẋ = f(x, µ), x E R n, (30) and assume the general case where = x 0, is a fixed point at µ = µ 0 and there is a single eigenvalue that crosses 0 at µ = µ 0. We can obtain the form of the bifurcation without computing the extended centre manifold every time. 27

34 Sotomayor s theorem (1973) Define D = Df, (31) (x=x0,µ=µ 0) D has a single zero eigenvalue with the left and right eigenvectors Dv = 0, and wd = 0. (32) Define α = 1 v w w f µ (x=x0,µ=µ 0) β = 1 n ( ) 2 f i w i v j v k v w x j x k (x=x0,µ=µ 0) i,j,k=1 (33) (34) 28

35 Sotomayor s theorem (1973) Theorem 4.1 (Sotomayor) If α 0 β, then there exists a smooth curve of equilibrium point in R n R passing by (x 0, µ 0 ) and tangent to R n {µ 0 } so that locally: either, there are no fixed points for µ < µ 0 and two for µ > µ 0 ; or, there are no fixed points for µ > µ 0 and two for µ < µ 0 ; This theorem states that if α 0 β, there is a saddle-node bifurcation. We define the centre manifold variable y as the dynamics along v: x(t) = x 0 + y(t)v, (35) λ = µ µ 0, (36) and the dynamics on the extended centre manifold is given by ẏ = αλ + βy 2. (37) 29

36 Sotomayor s theorem (1973) What happens if α = 0? In this case, we define γ = 2 v w n i,j=1 w i v j ( 2 f i x j µ) (x=x0,µ=µ 0) (38) and in the variable y and λ defined by x(t) = x 0 + y(t)v, (39) λ = µ µ 0, (40) the dynamics on the extended centre manifold is a transcritical bifurcation given by ẏ = γλy + βy 2. (41) Question: what happens if β = 0? 30

37 Example Consider the bifurcation at the origin for system {ẋ = (1 + µ)x 4y + x 2 2xy ẏ = 2x 4µy y 2 x 2 (42) Compute Df(0, 0) = [ 1 + µ 4 2 4µ ]. (43) There is a bifurcation at µ = 1 with matrix [ ] 2 4 D =, (44) 2 4 and left/right kernels [ ] 2 [ v =, w = ], v w = 1. (45) 31

38 Example Since, µ f(0, 0) = (0, 0), we have α = 0 and compute γ = 12, β = 10, (46) Writing [x, y] T = z(t)[2, 1] T, the bifurcation is transcritical ż = 12z(µ 1) + 10z 2. (47) 32

39 4. Local bifurcations for vector fields 4.5 Local bifurcation of maps

40 Bifurcation of maps Consider the mapping x n+1 = G(x n ), (48) where x E R m. Assume that x 0 is a fixed point (G(x 0 ) = x 0 ). This fixed point is as. stable if λ < 1 for all λ Spec(DG(x 0 )). The fixed point is unstable if there λ Spec(DG(x 0 )) s.t. λ > 1. So bifurcation will occur when an eigenvalue is on the unit complex circle. 33

41 Bifurcation of maps We consider a mapping with one parameter µ x n+1 = G(x n, µ), (49) where x E R m. Assume that x 0 = x 0 (µ) is a fixed point (G(x 0, µ) = x 0 ). We are interested in the case where one of the eigenvalue crosses the unit disk. It gives three possibilities at the bifurcation: either λ = 1, λ = 1 or λ λ with λ = 1. 34

42 Bifurcation of maps: λ = 1 at bifurcation Case I: λ = 1 at bifurcation. This case is similar to the cases obtained for vector fields. Namely, we have the saddle-node bifurcation: The transcritical bifurcation: The pitchfork bifurcation: x x + µ x 2. (50) x x + µx x 2 (51) x x + µx x 3 (52) 35

43 Bifurcation of periodic orbit Example in polar coordinates: {ṙ = r ( µ (r 2 1) 2) θ = 1 (53) 36

44 Bifurcation of maps: λ = 1 at bifurcation Case II: λ = 1 at bifurcation. Period-doubling bifurcation x x µx + x 3 (54) 37

45 The logistic map An important example of period-doubling cascade leading to chaotic dynamics is the logistic map (See file LogisticMap.pdf) x µx(1 x) (55) 38

### B5.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

### Part 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

### 7 Two-dimensional bifurcations

7 Two-dimensional bifurcations As in one-dimensional systems: fixed points may be created, destroyed, or change stability as parameters are varied (change of topological equivalence ). In addition closed

### One Dimensional Dynamical Systems

16 CHAPTER 2 One Dimensional Dynamical Systems We begin by analyzing some dynamical systems with one-dimensional phase spaces, and in particular their bifurcations. All equations in this Chapter are scalar

### Problem 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

### 2.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

### Boulder School for Condensed Matter and Materials Physics. Laurette Tuckerman PMMH-ESPCI-CNRS

Boulder School for Condensed Matter and Materials Physics Laurette Tuckerman PMMH-ESPCI-CNRS laurette@pmmh.espci.fr Dynamical Systems: A Basic Primer 1 1 Basic bifurcations 1.1 Fixed points and linear

### Bifurcation of Fixed Points

Bifurcation of Fixed Points CDS140B Lecturer: Wang Sang Koon Winter, 2005 1 Introduction ẏ = g(y, λ). where y R n, λ R p. Suppose it has a fixed point at (y 0, λ 0 ), i.e., g(y 0, λ 0 ) = 0. Two Questions:

### Half 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

### EE222 - Spring 16 - Lecture 2 Notes 1

EE222 - Spring 16 - Lecture 2 Notes 1 Murat Arcak January 21 2016 1 Licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License. Essentially Nonlinear Phenomena Continued

### B5.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

### CHALMERS, 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,

### APPPHYS217 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 (Springer-Verlag

### 2 Lecture 2: Amplitude equations and Hopf bifurcations

Lecture : Amplitude equations and Hopf bifurcations This lecture completes the brief discussion of steady-state bifurcations by discussing vector fields that describe the dynamics near a bifurcation. From

### Summary of topics relevant for the final. p. 1

Summary of topics relevant for the final p. 1 Outline Scalar difference equations General theory of ODEs Linear ODEs Linear maps Analysis near fixed points (linearization) Bifurcations How to analyze a

### STABILITY. 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),

### Nonlinear 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

### Entrance 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

### Dynamical systems tutorial. Gregor Schöner, INI, RUB

Dynamical systems tutorial Gregor Schöner, INI, RUB Dynamical systems: Tutorial the word dynamics time-varying measures range of a quantity forces causing/accounting for movement => dynamical systems dynamical

### TWO 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

### 1. < 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

### CHALMERS, 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

### Mathematical 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

### 8.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 integro-differential equations

### 1. (i) Determine how many periodic orbits and equilibria can be born at the bifurcations of the zero equilibrium of the following system:

1. (i) Determine how many periodic orbits and equilibria can be born at the bifurcations of the zero equilibrium of the following system: ẋ = y x 2, ẏ = z + xy, ż = y z + x 2 xy + y 2 + z 2 x 4. (ii) Determine

### = F ( x; µ) (1) where x is a 2-dimensional 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

### CHALMERS, 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

### 10 Back to planar nonlinear systems

10 Back to planar nonlinear sstems 10.1 Near the equilibria Recall that I started talking about the Lotka Volterra model as a motivation to stud sstems of two first order autonomous equations of the form

### MATH 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

### BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs

BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs Yuri A. Kuznetsov August, 2010 Contents 1. Solutions and orbits. 2. Equilibria. 3. Periodic orbits and limit cycles. 4. Homoclinic orbits.

### Math 273 (51) - Final

Name: Id #: Math 273 (5) - Final Autumn Quarter 26 Thursday, December 8, 26-6: to 8: Instructions: Prob. Points Score possible 25 2 25 3 25 TOTAL 75 Read each problem carefully. Write legibly. Show all

### April 13, We now extend the structure of the horseshoe to more general kinds of invariant. x (v) λ n v.

April 3, 005 - Hyperbolic Sets We now extend the structure of the horseshoe to more general kinds of invariant sets. Let r, and let f D r (M) where M is a Riemannian manifold. A compact f invariant set

### 7 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

### Def. (a, b) is a critical point of the autonomous system. 1 Proper node (stable or unstable) 2 Improper node (stable or unstable)

Types of critical points Def. (a, b) is a critical point of the autonomous system Math 216 Differential Equations Kenneth Harris kaharri@umich.edu Department of Mathematics University of Michigan November

### MATH 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 one-dimensional systems, we saw that subtle changes in parameter

### 5.2.2 Planar Andronov-Hopf bifurcation

138 CHAPTER 5. LOCAL BIFURCATION THEORY 5.. Planar Andronov-Hopf bifurcation What happens if a planar system has an equilibrium x = x 0 at some parameter value α = α 0 with eigenvalues λ 1, = ±iω 0, ω

### Practice 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,

### Bifurcation Analysis of Non-linear Differential Equations

Bifurcation Analysis of Non-linear Differential Equations Caitlin McCann 0064570 Supervisor: Dr. Vasiev September 01 - May 013 Contents 1 Introduction 3 Definitions 4 3 Ordinary Differential Equations

### 11 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

### EN 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

### MCE693/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

### Chapter 2 Hopf Bifurcation and Normal Form Computation

Chapter 2 Hopf Bifurcation and Normal Form Computation In this chapter, we discuss the computation of normal forms. First we present a general approach which combines center manifold theory with computation

### Chapter 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é (1854-1912)

### Nonlinear 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

### Problem List MATH 5173 Spring, 2014

Problem List MATH 5173 Spring, 2014 The notation p/n means the problem with number n on page p of Perko. 1. 5/3 [Due Wednesday, January 15] 2. 6/5 and describe the relationship of the phase portraits [Due

UNIVERSIDADE DE SÃO PAULO Instituto de Ciências Matemáticas e de Computação ISSN 0103-2577 GLOBAL DYNAMICAL ASPECTS OF A GENERALIZED CHEN-WANG DIFFERENTIAL SYSTEM REGILENE D. S. OLIVEIRA CLAUDIA VALLS

### CHAPTER 6 HOPF-BIFURCATION IN A TWO DIMENSIONAL NONLINEAR DIFFERENTIAL EQUATION

CHAPTER 6 HOPF-BIFURCATION IN A TWO DIMENSIONAL NONLINEAR DIFFERENTIAL EQUATION [Discussion on this chapter is based on our paper entitled Hopf-Bifurcation Ina Two Dimensional Nonlinear Differential Equation,

### Lecture 5. Outline: Limit Cycles. Definition and examples How to rule out limit cycles. Poincare-Bendixson theorem Hopf bifurcations Poincare maps

Lecture 5 Outline: Limit Cycles Definition and examples How to rule out limit cycles Gradient systems Liapunov functions Dulacs criterion Poincare-Bendixson theorem Hopf bifurcations Poincare maps Limit

### Calculus 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,

### Dynamics and Bifurcations in Predator-Prey Models with Refuge, Dispersal and Threshold Harvesting

Dynamics and Bifurcations in Predator-Prey Models with Refuge, Dispersal and Threshold Harvesting August 2012 Overview Last Model ẋ = αx(1 x) a(1 m)xy 1+c(1 m)x H(x) ẏ = dy + b(1 m)xy 1+c(1 m)x (1) where

### WIDELY SEPARATED FREQUENCIES IN COUPLED OSCILLATORS WITH ENERGY-PRESERVING QUADRATIC NONLINEARITY

WIDELY SEPARATED FREQUENCIES IN COUPLED OSCILLATORS WITH ENERGY-PRESERVING QUADRATIC NONLINEARITY J.M. TUWANKOTTA Abstract. In this paper we present an analysis of a system of coupled oscillators suggested

### CHARACTERIZATION OF SADDLE-NODE EQUILIBRIUM POINTS ON THE STABILITY BOUNDARY OF NONLINEAR AUTONOMOUS DYNAMICAL SYSTEM

CHARACTERIZATION OF SADDLE-NODE EQUILIBRIUM POINTS ON THE STABILITY BOUNDARY OF NONLINEAR AUTONOMOUS DYNAMICAL SYSTEM Fabíolo Moraes Amaral, Josaphat Ricardo Ribeiro Gouveia Júnior, Luís Fernando Costa

### Chapter 1 Bifurcations and Chaos in Dynamical Systems

Chapter 1 Bifurcations and Chaos in Dynamical Systems Complex system theory deals with dynamical systems containing often a large number of variables. It extends dynamical system theory, which deals with

### Complex 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

### Characterization of the stability boundary of nonlinear autonomous dynamical systems in the presence of a saddle-node equilibrium point of type 0

Anais do CNMAC v.2 ISSN 1984-82X Characterization of the stability boundary of nonlinear autonomous dynamical systems in the presence of a saddle-node equilibrium point of type Fabíolo M. Amaral Departamento

### CENTER MANIFOLD AND NORMAL FORM THEORIES

3 rd Sperlonga Summer School on Mechanics and Engineering Sciences 3-7 September 013 SPERLONGA CENTER MANIFOLD AND NORMAL FORM THEORIES ANGELO LUONGO 1 THE CENTER MANIFOLD METHOD Existence of an invariant

### Andronov Hopf and Bautin bifurcation in a tritrophic food chain model with Holling functional response types IV and II

Electronic Journal of Qualitative Theory of Differential Equations 018 No 78 1 7; https://doiorg/10143/ejqtde018178 wwwmathu-szegedhu/ejqtde/ Andronov Hopf and Bautin bifurcation in a tritrophic food chain

### 4 Second-Order Systems

4 Second-Order Systems Second-order 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

### Stability 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

### Solutions to Dynamical Systems 2010 exam. Each question is worth 25 marks.

Solutions to Dynamical Systems exam Each question is worth marks [Unseen] Consider the following st order differential equation: dy dt Xy yy 4 a Find and classify all the fixed points of Hence draw the

### Invariant Manifolds of Dynamical Systems and an application to Space Exploration

Invariant Manifolds of Dynamical Systems and an application to Space Exploration Mateo Wirth January 13, 2014 1 Abstract In this paper we go over the basics of stable and unstable manifolds associated

### 2 Discrete growth models, logistic map (Murray, Chapter 2)

2 Discrete growth models, logistic map (Murray, Chapter 2) As argued in Lecture 1 the population of non-overlapping generations can be modelled as a discrete dynamical system. This is an example of an

### Nonlinear dynamics & chaos BECS

Nonlinear dynamics & chaos BECS-114.7151 Phase portraits Focus: nonlinear systems in two dimensions General form of a vector field on the phase plane: Vector notation: Phase portraits Solution x(t) describes

### 2 Qualitative theory of non-smooth dynamical systems

2 Qualitative theory of non-smooth dynamical systems In this chapter, we give an overview of the basic theory of both smooth and non-smooth dynamical systems, to be expanded upon in later chapters. In

### TWELVE LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM WITH Z 2 -SYMMETRY. P. Yu 1,2 and M. Han 1

COMMUNICATIONS ON Website: http://aimsciences.org PURE AND APPLIED ANALYSIS Volume 3, Number 3, September 2004 pp. 515 526 TWELVE LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM WITH Z 2 -SYMMETRY P. Yu 1,2

### Department of Mathematics IIT Guwahati

Stability of Linear Systems in R 2 Department of Mathematics IIT Guwahati A system of first order differential equations is called autonomous if the system can be written in the form dx 1 dt = g 1(x 1,

UNIVERSIDADE DE SÃO PAULO Instituto de Ciências Matemáticas e de Computação ISSN 010-577 GLOBAL DYNAMICAL ASPECTS OF A GENERALIZED SPROTT E DIFFERENTIAL SYSTEM REGILENE OLIVEIRA CLAUDIA VALLS N o 41 NOTAS

### (8.51) ẋ = A(λ)x + F(x, λ), where λ lr, the matrix A(λ) and function F(x, λ) are C k -functions with k 1,

2.8.7. Poincaré-Andronov-Hopf Bifurcation. In the previous section, we have given a rather detailed method for determining the periodic orbits of a two dimensional system which is the perturbation of a

### Math 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

### ME 680- Spring Geometrical Analysis of 1-D Dynamical Systems

ME 680- Spring 2014 Geometrical Analysis of 1-D Dynamical Systems 1 Geometrical Analysis of 1-D Dynamical Systems Logistic equation: n = rn(1 n) velocity function Equilibria or fied points : initial conditions

### DYNAMICAL SYSTEMS WITH A CODIMENSION-ONE INVARIANT MANIFOLD: THE UNFOLDINGS AND ITS BIFURCATIONS

International Journal of Bifurcation and Chaos c World Scientific Publishing Company DYNAMICAL SYSTEMS WITH A CODIMENSION-ONE INVARIANT MANIFOLD: THE UNFOLDINGS AND ITS BIFURCATIONS KIE VAN IVANKY SAPUTRA

### Stability lectures. Stability of Linear Systems. Stability of Linear Systems. Stability of Continuous Systems. EECE 571M/491M, Spring 2008 Lecture 5

EECE 571M/491M, Spring 2008 Lecture 5 Stability of Continuous Systems http://courses.ece.ubc.ca/491m moishi@ece.ubc.ca Dr. Meeko Oishi Electrical and Computer Engineering University of British Columbia,

### tutorial ii: One-parameter bifurcation analysis of equilibria with matcont

tutorial ii: One-parameter bifurcation analysis of equilibria with matcont Yu.A. Kuznetsov Department of Mathematics Utrecht University Budapestlaan 6 3508 TA, Utrecht February 13, 2018 1 This session

### Lotka Volterra Predator-Prey Model with a Predating Scavenger

Lotka Volterra Predator-Prey 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

### Lectures on Dynamical Systems. Anatoly Neishtadt

Lectures on Dynamical Systems Anatoly Neishtadt Lectures for Mathematics Access Grid Instruction and Collaboration (MAGIC) consortium, Loughborough University, 2007 Part 3 LECTURE 14 NORMAL FORMS Resonances

### Solution to Homework #4 Roy Malka

1. Show that the map: Solution to Homework #4 Roy Malka F µ : x n+1 = x n + µ x 2 n x n R (1) undergoes a saddle-node bifurcation at (x, µ) = (0, 0); show that for µ < 0 it has no fixed points whereas

### Applied 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,

### Stability 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,

### Hopf Bifurcation and Limit Cycle Analysis of the Rikitake System

ISSN 749-3889 (print), 749-3897 (online) International Journal of Nonlinear Science Vol.4(0) No.,pp.-5 Hopf Bifurcation and Limit Cycle Analysis of the Rikitake System Xuedi Wang, Tianyu Yang, Wei Xu Nonlinear

### Mathematical 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

### Problem Set Number 02, j/2.036j MIT (Fall 2018)

Problem Set Number 0, 18.385j/.036j MIT (Fall 018) Rodolfo R. Rosales (MIT, Math. Dept., room -337, Cambridge, MA 0139) September 6, 018 Due October 4, 018. Turn it in (by 3PM) at the Math. Problem Set

### The Higgins-Selkov oscillator

The Higgins-Selkov oscillator May 14, 2014 Here I analyse the long-time behaviour of the Higgins-Selkov oscillator. The system is ẋ = k 0 k 1 xy 2, (1 ẏ = k 1 xy 2 k 2 y. (2 The unknowns x and y, being

### Math 312 Lecture Notes Linearization

Math 3 Lecture Notes Linearization Warren Weckesser Department of Mathematics Colgate University 3 March 005 These notes discuss linearization, in which a linear system is used to approximate the behavior

### CDS 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/

### 2 Lyapunov Stability. x(0) x 0 < δ x(t) x 0 < ɛ

1 2 Lyapunov Stability Whereas I/O stability is concerned with the effect of inputs on outputs, Lyapunov stability deals with unforced systems: ẋ = f(x, t) (1) where x R n, t R +, and f : R n R + R n.

### Models 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 Predator-prey models are used to show the intricate

### Chapter 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 first-order system for x(t) d dt x = f(x,t), x(0) = x 0 In particular, consider

### Introduction to Applied Nonlinear Dynamical Systems and Chaos

Stephen Wiggins Introduction to Applied Nonlinear Dynamical Systems and Chaos Second Edition With 250 Figures 4jj Springer I Series Preface v L I Preface to the Second Edition vii Introduction 1 1 Equilibrium

### CONTROLLING IN BETWEEN THE LORENZ AND THE CHEN SYSTEMS

International Journal of Bifurcation and Chaos, Vol. 12, No. 6 (22) 1417 1422 c World Scientific Publishing Company CONTROLLING IN BETWEEN THE LORENZ AND THE CHEN SYSTEMS JINHU LÜ Institute of Systems

### A Study of the Van der Pol Equation

A Study of the Van der Pol Equation Kai Zhe Tan, s1465711 September 16, 2016 Abstract The Van der Pol equation is famous for modelling biological systems as well as being a good model to study its multiple

### A plane autonomous system is a pair of simultaneous first-order differential equations,

Chapter 11 Phase-Plane Techniques 11.1 Plane Autonomous Systems A plane autonomous system is a pair of simultaneous first-order differential equations, ẋ = f(x, y), ẏ = g(x, y). This system has an equilibrium

### Continuous 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

### Handout 2: Invariant Sets and Stability

Engineering Tripos Part IIB Nonlinear Systems and Control Module 4F2 1 Invariant Sets Handout 2: Invariant Sets and Stability Consider again the autonomous dynamical system ẋ = f(x), x() = x (1) with state

### Introduction to Bifurcation and Normal Form theories

Introduction to Bifurcation and Normal Form theories Romain Veltz / Olivier Faugeras October 9th 2013 ENS - Master MVA / Paris 6 - Master Maths-Bio (2013-2014) Outline 1 Invariant sets Limit cycles Stable/Unstable

### On a Codimension Three Bifurcation Arising in a Simple Dynamo Model

On a Codimension Three Bifurcation Arising in a Simple Dynamo Model Anne C. Skeldon a,1 and Irene M. Moroz b a Department of Mathematics, City University, Northampton Square, London EC1V 0HB, England b

### Numerical techniques: Deterministic Dynamical Systems

Numerical techniques: Deterministic Dynamical Systems Henk Dijkstra Institute for Marine and Atmospheric research Utrecht, Department of Physics and Astronomy, Utrecht, The Netherlands Transition behavior

### Chapter #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.

### Clearly the passage of an eigenvalue through to the positive real half plane leads to a qualitative change in the phase portrait, i.e.

Bifurcations We have already seen how the loss of stiffness in a linear oscillator leads to instability. In a practical situation the stiffness may not degrade in a linear fashion, and instability may