# Analysis of the Takens-Bogdanov bifurcation on m parameterized vector fields

Save this PDF as:

Size: px
Start display at page:

Download "Analysis of the Takens-Bogdanov bifurcation on m parameterized vector fields" ## Transcription

1 Analysis of the Takens-Bogdanov bifurcation on m parameterized vector fields Francisco Armando Carrillo Navarro, Fernando Verduzco G., Joaquín Delgado F. Programa de Doctorado en Ciencias (Matemáticas), Departamento de Matemáticas, Universidad de Sonora.

2 Contenido 1 Takens-Bogdanov bifurcation 2 Statement of the problem 3 Dynamics on the center manifold 4 Main theorem 5 An example

3 Takens-Bogdanov bifurcation We consider a planar dynamical system ż = f (z), with f (0) = 0 and ( ) 0 1 Df (0) = J (1) From normal form theory, there is a change of coordinates such that the original system can be reduced up to terms of second order to the form

4 ż = f 0 (z) = ( z 2 az bz 1z 2 ). (2) This is called the truncated normal form of the original system. Bogdanov  showed that the family ( ) z ż = 2 λ 1 + λ 2 z 1 + az1 2 + bz, (3) 1z 2 for ab 0 is a versal deformation of the truncated system (2).

5 Figure: Different dynamics of the Takens-Bogdanov bifurcation.

6 The point λ 1 = λ 2 = 0 separates two branches of the saddle-node bifurcation curve T +( ) = {(λ 1, λ 2 ) : λ 1 = 1 4 λ2 2, λ 2 > 0 (λ 2 < 0)}, the half-line H = {(λ 1, λ 2 ) : λ 1 = 0, λ 2 < 0} corresponds to the Hopf bifurcation that generates a stable limit cycle.

7 This cycle exists and remains hyperbolic between the line H and a smooth curve P = {(λ 1, λ 2 ) : λ 1 = 6 25 λ2 2+O( λ 2 3 ), λ 2 < 0}, where saddle homoclinic bifurcation occurs. When the cycle approaches the homoclinic orbit, its period tends to infinity.

8 Statement of the problem Consider the m-parameterized vector field ẋ = F (x, µ), (4) where x R n, µ R m, with m 2, and F C r (R n R m ), with r 2. Suppose that there exists (x 0, µ 0 ) R n R m such that H1) F (x 0, µ 0 ) = 0, and

9 H2) The spectrum of the linearization at the critical point satisfies σ (DF (x 0, µ 0 )) = {λ j C λ 1,2 = 0, Re(λ j ) 0, for j = 3,..., n }, the nonsemisimple case.

10 Our goal is to find sufficient conditions on the vector field F, such that the dynamics on the center manifold at x = x 0, is locally topologically equivalent to the versal deformation of the generic planar Takens-Bogdanov bifurcation

11 ż 1 = z 2 (5) ż 2 = β 1 + β 2 z 1 + az bz 1 z 2, where ab 0. (6)

12 Dynamics on the center manifold We use the center manifold theory to determine the dynamics on the m-parameterized two-dimensional center manifold at the equilibrium x = x 0 for µ µ 0.

13 Jordan form Let us consider the Taylor series around (x 0, µ 0 ), F (x, µ) = DF (x 0, µ 0 )(x x 0 ) + F µ (x 0, µ 0 )(µ µ 0 ) D2 F (x 0, µ 0 )(x x 0, x x 0 ) +F µx (x 0, µ 0 )(µ µ 0, x x 0 ) +. (7) From hypothesis H2) we have that,

14 The matrix A = DF (x 0, µ 0 ) R n n is similar to the matrix ( ) J0 0 J =, 0 J 1

15 where J 0 = ( ), and J 1 R (n 2) (n 2) is such that σ(j 1 ) = {λ j C Re(λ j ) 0, for j = 3,..., n }

16 That is, there exists a invertible matrix P = (p 1 p 2 P 0 ), such that P 1 AP = J, with P 1 = q T 1 q T 2 Q 0,

17 Change of coordinates Let us consider the change of coordinates and parameters y = P 1 (x x 0 ), and α = µ µ 0, (8) then (7) transforms into ẏ = Jy + P 1 F µ (x 0, µ 0 )α P 1 D 2 F (x 0, µ 0 )(Py, Py) +P 1 F µx (x 0, µ 0 )(α, Py) +. (9)

18 Let us define v 0 = ( p 1 p 2 ), and w 0 = ( q ( 1 q 2 ) T, then P = ( v 0 P 0 ) and P 1 w0 =. Now we define ( y1 y 2 ) Q 0 = y = P 1 (x x 0 ) = that is, y 1 = w 0 (x x 0 ) R 2, and y 2 = Q 0 (x x 0 ) R n 2. ( w0 (x x 0 ) Q 0 (x x 0 ) ),

19 then (9) transforms into ẏ 1 = J 0 y 1 + w 0 F µ (x 0, µ 0 )α w 0D 2 F (x 0, µ 0 )(v 0 y 1, v 0 y 1 ) + w 0 D 2 F (x 0, µ 0 )(v 0 y 1, P 0 y 2 ) w 0D 2 F (x 0, µ 0 )(P 0 y 2, P 0 y 2 ) +w 0 F µx (x 0, µ 0 )(α, v 0 y 1 ) + w 0 F µx (x 0, µ 0 )(α, P 0 y 2 ) +

20 ẏ 2 = J 1 y 2 + Q 0 F µ (x 0, µ 0 )α Q 0D 2 F (x 0, µ 0 )(v 0 y 1, v 0 y 1 ) + Q 0 D 2 F (x 0, µ 0 )(v 0 y 1, P 0 y 2 ) Q 0D 2 F (x 0, µ 0 )(P 0 y 2, P 0 y 2 ) +Q 0 F µx (x 0, µ 0 )(α, v 0 y 1 ) + Q 0 F µx (x 0, µ 0 )(α, P 0 y 2 ) +.

21 In order to simplify this system we will use the next Definition Given ν R n, ν = L = L 1. L n product by ν 1. ν n and L R n (r s),, where L i R r s, we define the ν L = n ν i L i. i=1

22 Therefore, the system can be written as the extended system ẏ 1 J 0 w 0 F µ (x 0, µ 0 ) 0 y 1 α = α ẏ 2 0 Q 0 F µ (x 0, µ 0 ) J 1 y 2 where F 1 (y 1, α, y 2 ) 0 F 2 (y 1, α, y 2 ) +, (10)

23 F 1 (y 1, α, y 2 ) = 1 [( w0 D 2 F (x 0, µ 0 ) ) (v 0, v 0 ) ] (y 1, y 1 ) 2 + [( w 0 D 2 F (x 0, µ 0 ) ) (v 0, P 0 ) ] (y 1, y 2 ) + 1 [( w0 D 2 F (x 0, µ 0 ) ) (P 0, P 0 ) ] (y 2, y 2 ) 2 + [(w 0 F µx (x 0, µ 0 )) v 0 ] (α, y 1 ) + [(w 0 F µx (x 0, µ 0 )) P 0 ] (α, y 2 ) +,

24 To calculate the m-parameterized local center manifold at the equilibrium point y = 0, we first consider the change coordinates where P = ξ α ζ I I m 0 0 K 2 I n 2 = P 1 y 1 α y 2 with K 2 = J 1 1 Q 0F µ (x 0, µ 0 ). and P 1 =, (11) I I m 0 0 K 2 I n 2,

25 Then, system (10) transforms into ξ J 0 R 0 0 ξ α = α + ζ 0 0 J 1 ζ where R 0 = w 0 F µ (x 0, µ 0 ), and f 1 (ξ, α, ζ) = F 1 (ξ, α, ζ K 2 α), f 2 (ξ, α, ζ) = F 2 (ξ, α, ζ K 2 α). f 1 (ξ, α, ζ) 0 f 2 (ξ, α, ζ) (12)

26 From the center manifold theory, system (12) has a center manifold ζ = h(ξ, α) = O( ξ, α 2 ), with h(0, 0) = Dh(0, 0) = 0, and the dynamics on the center manifold is given by ξ = J 0 ξ + R 0 α + f 1 (ξ, α, h(ξ, α))

27 where f 1 (ξ, α, h(ξ, α)) = F 1 (ξ, α, h(ξ, α) J 1 1 Q 0F µ (x 0, µ 0 )α). Observe that it is not necessary calculate the center manifold h(ξ, α) because it does not affect the quadratic terms in f 1.

28 Then, we have proved the next Lemma Let the nonlinear system ẋ = F (x, µ), satisfy the non-hyperbolicity conditions H1) H2) at the equilibrium point (x 0, µ 0 ). Then the dynamics on the m-parameterized two-dimensional center manifold at the equilibrium point x = x 0 for µ µ 0, is given by

29 ξ = J 0 ξ + R 0 α + R 1 (α, ξ) R 2(ξ, ξ) where +O( α 2 ) + O( ξ, α 3 ), (13) ξ = w 0 (x x 0 ), α = µ µ 0, R 0 = w 0 F µ (x 0, µ 0 ), (14)

30 and R 1 = (w 0 F µx (x 0, µ 0 )) v 0 (15) + ( w 0 D 2 F (x 0, µ 0 ) ) (A 0 F µ (x 0, µ 0 ), v 0 ), R 2 = ( w 0 D 2 F (x 0, µ 0 ) ) (v 0, v 0 ), (16) with A 0 = P 0 J 1 1 Q 0.

31 Main theorem Then our goal is to find a change of coordinates ξ = z + L 0 α + α T L 1 z zt L 2 z, where L 0 R 2 m, L 1 R 2 (m 2), L 2 R 2 (2 2), such that (13) is transformed into the versal deformation of the Takens-Bogdanov bifurcation (5).

32 Observe that ξ = (I + α T L 1 + z T L 2 )ż, and for z 0 we have that (I +α T L 1 +z T L 2 ) 1 = I α T L 1 z T L 2 +, then, ż = J 0 z + R 0 α + α T R1 z zt R2 z +,

33 where, R 0 = J 0 L 0 + R 0, R 1 = L 1 + R 1 L 1 J 0 + L T 0 R 2 R 0 T L 2 R 2 = L 2 + R 2 2L 2 J 0, with L i = (L i2, 0) T.

34 Lemma There exists L 0 such that R 0 α = β 1 e 2. Proof: Observe that R 0 = w 0 F µ (x 0, µ 0 ) = ( R T 01 = ( q T 1 F µ (x 0, µ 0 ) q T 2 F µ(x 0, µ 0 ) ), R T 02 )

35 Then, if we define L 02 = F T µ (x 0, µ 0 )q 1, and the results it follows. β 1 = q T 2 F µ (x 0, µ 0 )α,

36 Lemma There exists L 2 such that 1 2 z T R2 z = (az bz 1 z 2 )e 2. Proof: If we define D i = q i D 2 F (x 0, µ 0 ), and ( 1 L 21 = 2 pt 2 D 2p 2 + p1 T D ) 1p 2 l12 1 l12 1 l22 1, ( p T 1 L 22 = 1 D 1 p 1 2 pt 2 D ) 2p pt 2 D 2p 2 2l12 1 pt 2 D, 1p 2

37 where l12 1 and l 22 1 are free. The results follows with a = 1 2 pt 1 D 2 p 1, b = p1 T D 1 p 1 + p1 T D 2 p 2. (17)

38 Lemma There exist L 01 and L 12 such that Proof: If α T R1 z = β 2 z 1 e 2. L 01 = 1 b (( ) p T 1 D 1 p 2 + p2 T D 2 p 2 R01 + ( ) l l22 2 R02 R11 1 R12) 2, l12 1 = 1 ( 3 pt 2 D 1 + b ) 4a D 2 p 2, and L 12 = L 11 J 0 R 11 L T 0 R 21 + R T 0 L 21,

39 then, R 11 1 = R 11 2 = R 12 2 = 0, and ( ) α T 0 ) R1 z = α T = (α T R1 R1 12 z 12 z 1 e 2, 1 then, if we define β 2 = α T R1 12, the results follows.

40 Before establish the main theorem, let us define S 1 = F T µ (x 0, µ 0 )q 2, S 2 = R 12 1 (18) [ ] 2a = b (pt 1 D 1 p 2 + p2 T D 2 p 2 ) p1 T D 2 p 2 Fµ T (x 0, µ 0 )q 1 2a 2 R1i i + R 1 b 12. (19) i=1

41 To unfold the double-zero bifurcation from the center manifold, the transformation T : R m R 2, given by ( ) S T T (µ) = 1 (µ µ 0 ) S2 T(µ µ 0) must have rank two, that is, S 1 and S 2 must to be linearly independent.

42 Theorem Given the nonlinear system ẋ = F (x, µ), (20) where x R n, µ R m with m 2, such that, there exists (x 0, µ 0 ), that satisfies the conditions

43 H1) F (x 0, µ 0 ) = 0, H2) σ[df (x 0, µ 0 )] = {λ 1,2 = 0; Re(λ j ) 0, j = 3,..., n }, (non-hyperbolicity) H3) ab 0, (nondegeneracy) H4) S 1 and S 2 are linearly independent, (transversality)

44 Then, the dynamics on the center manifold of system (20) at x = x 0 and µ µ 0, which is given by (13), is locally topologically equivalent to the versal deformation of the Takens-Bogdanov bifurcation ż 1 = z 2 ż 2 = β 1 + β 2 z 1 + az1 2 + bz 1z 2, where β 1 = S T 1 (µ µ 0), and β 2 = S T 2 (µ µ 0).

45 An example Consider the following system of two differential equations: ẋ 1 = x 1 x 1x 2 1+αx 1 εx 2 1, ẋ 2 = γx 2 + x 1x 2 1+αx 1 δx 2 2, (21)

46 The equations model the dynamics of a predator-prey ecosystem. The variables x 1 and x 2 are (scaled) population numbers of prey and predator, respectively

47 while α, γ, ε, and δ are nonnegative parameters describing the behavior of isolated populations and their interaction.

48 Assume that ε = 0 is fixed. In order to the bifurcation diagram of the system with respect to the three remaining parameters (α, γ, δ) can exhibit the codim 2 Takens-Bogdanov bifurcation in planar systems, we consider that α (0, 1) and αγ < 1.

49 If we consider µ = (α, γ, δ) T, then (x 0, µ 0 ) = ( ( 1 ( α, α + 1 α 2, α ) T, 1 α α(1 + α), α (1 + α) 2 ) T ) is a family of equilibrium points whose linearization has a double-zero eigenvalue,

50 and p 1 = ( 1 α q 1 = ( 0 1 α ) ( ) 1 + α, p 2 =, 0 ) ( ), q 2 =, 1 1+α 1 α(1+α) are the right and left (generalized) eigenvectors, respectively, associated to the eigenvalue zero.

51 Besides, from (17), a = 1 2 pt 1 D 2 p 1 = α 2 (1 + α) 3, b = p T 1 D 1 p 1 + p T 1 D 2 p 2 = α2 (α 1) (α + 1) 2,

52 and, from (19), ( ) S1 T 1 = α 4 (1 + α), α α 2,, α 3 ( 2(α S2 T 2 2α 1) = α 2 (α + 1) 2 (α 1), α 2 3α 2, α 2 1 α 2 ) 2α 1. α(α 1)

53 Thus the conditions H1), H2), H3) and H4) are satisfied for all α (0, 1). Therefore, if we choose α = 1 2, the vector field (21) is locally topologically equivalent to ż 1 = z 2, ż 2 = β 1 + β 2 z z z 1z 2, where

54 with β 1 =< S 1, µ >, and β 2 =< S 2, µ > µ T = (α 1 2, γ 2 3, δ 2 9 ).

55 Figure: Surfaces of bifurcation: S sn surface of saddle-node bifurcation, S H surface of the Hopf bifurcation and S hom surface of homoclinic bifurcation.

56

57 Referencias Bazykin, A. Mathematical Biophysics of Interacting Populations. Nauka, Mmoskow. In Russian Bogdanov, R.I. Versal deformations of a singular point of a vector field on the plane in the case of zero eigenvalues. Funct. Anal. Appl. 9, Carrillo, F.A. and Verduzco, F. Control of the Planar Takens-Bogdanov Bifurcation with Applications. Acta Appl. Math. 105: Guckenheimer, J., Holmes, P. Nonlinear oscillations, dynamical systems, and bifurcations of vector fields. Applied Mathematical Sciences, Vol. 42. Springer-Verlag Kuznetsov, Y.A. Elements of Applied Bifurcation Theory. Applied Mathematical Sciences, Vol Springer. Second Edition Takens, F. Forced osccilations and bifurcations. Comm. Math. Inst. Rijkuniversiteit Utrecht, 3, S. Wiggins. Introduction to Applied Nonlinear Dynamical Systems and Chaos. Texts in Applied Mathematics, Vol. 2. Springer-Verlag. Second Edition

### Elements of Applied Bifurcation Theory Yuri A. Kuznetsov Elements of Applied Bifurcation Theory Third Edition With 251 Illustrations Springer Introduction to Dynamical Systems 1 1.1 Definition of a dynamical system 1 1.1.1 State space 1 1.1.2

More information

### Elements of Applied Bifurcation Theory Yuri A. Kuznetsov Elements of Applied Bifurcation Theory Second Edition With 251 Illustrations Springer Preface to the Second Edition Preface to the First Edition vii ix 1 Introduction to Dynamical Systems

More information

### BIFURCATION PHENOMENA Lecture 4: Bifurcations in n-dimensional ODEs BIFURCATION PHENOMENA Lecture 4: Bifurcations in n-dimensional ODEs Yuri A. Kuznetsov August, 2010 Contents 1. Solutions and orbits: equilibria cycles connecting orbits compact invariant manifolds strange

More information

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

More information

### Elements of Applied Bifurcation Theory Yuri A. Kuznetsov Elements of Applied Bifurcation Theory Third Edition With 251 Illustrations Springer Yuri A. Kuznetsov Department of Mathematics Utrecht University Budapestlaan 6 3584 CD Utrecht The

More information

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

More information

### On the Takens-Bogdanov Bifurcation in the Chua s Equation 1722 PAPER Special Section on Nonlinear Theory and Its Applications On the Takens-Bogdanov Bifurcation in the Chua s Equation Antonio ALGABA, Emilio FREIRE, Estanislao GAMERO, and Alejandro J. RODRÍGUEZ-LUIS,

More information

### Dynamics of Modified Leslie-Gower Predator-Prey Model with Predator Harvesting International Journal of Basic & Applied Sciences IJBAS-IJENS Vol:13 No:05 55 Dynamics of Modified Leslie-Gower Predator-Prey Model with Predator Harvesting K. Saleh Department of Mathematics, King Fahd

More information

### Homoclinic Orbits of Planar Maps: Asymptotics and Mel nikov Functions Department of Mathematics Mathematical Sciences Homoclinic Orbits of Planar Maps: Asymptotics and Mel nikov Functions A thesis submitted for the degree of Master of Science Author: Dirk van Kekem Project

More information

### UNIVERSIDADE DE SÃO PAULO 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

More information

### ZERO-HOPF BIFURCATION FOR A CLASS OF LORENZ-TYPE SYSTEMS This is a preprint of: Zero-Hopf bifurcation for a class of Lorenz-type systems, Jaume Llibre, Ernesto Pérez-Chavela, Discrete Contin. Dyn. Syst. Ser. B, vol. 19(6), 1731 1736, 214. DOI: [doi:1.3934/dcdsb.214.19.1731]

More information

### Hopf Bifurcation Control for Affine Systems opf Bifurcation Control for Affine Systems Ferno Verduzco Joaquin Alvarez Abstract In this paper we establish conditions to control the opf bifurcation of nonlinear systems with two uncontrollable modes

More information

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

More information

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

More information

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

### Shilnikov bifurcations in the Hopf-zero singularity Shilnikov bifurcations in the Hopf-zero singularity Geometry and Dynamics in interaction Inma Baldomá, Oriol Castejón, Santiago Ibáñez, Tere M-Seara Observatoire de Paris, 15-17 December 2017, Paris Tere

More information

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

More information

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

More information

### Application demonstration. BifTools. Maple Package for Bifurcation Analysis in Dynamical Systems Application demonstration BifTools Maple Package for Bifurcation Analysis in Dynamical Systems Introduction Milen Borisov, Neli Dimitrova Department of Biomathematics Institute of Mathematics and Informatics

More information

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

More information

### IMPROVED HOMOCLINIC PREDICTOR FOR BOGDANOV-TAKENS BIFURCATION International Journal of Bifurcation and Chaos c World Scientific Publishing Company IMPROVED HOMOCLINIC PREDICTOR FOR BOGDANOV-TAKENS BIFURCATION YU.A. KUZNETSOV, H.G.E. MEIJER Department of Applied Mathematics,

More information

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

More information

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

More information

### Canards at Folded Nodes Canards at Folded Nodes John Guckenheimer and Radu Haiduc Mathematics Department, Ithaca, NY 14853 For Yulij Il yashenko with admiration and affection on the occasion of his 60th birthday March 18, 2003

More information

### The Invariant Curve in a Planar System of Difference Equations Advances in Dynamical Systems and Applications ISSN 0973-5321, Volume 13, Number 1, pp. 59 71 2018 http://campus.mst.edu/adsa The Invariant Curve in a Planar System of Difference Equations Senada Kalabušić,

More information

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

More information

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

More information

### arxiv: v1 [math.ds] 20 Sep 2016 THE PITCHFORK BIFURCATION arxiv:1609.05996v1 [math.ds] 20 Sep 2016 Contents INDIKA RAJAPAKSE AND STEVE SMALE Abstract. We give development of a new theory of the Pitchfork bifurcation, which removes the

More information

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

More information

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

More information

### CANARDS AT FOLDED NODES MOSCOW MATHEMATICAL JOURNAL Volume 5, Number 1, January March 2005, Pages 91 103 CANARDS AT FOLDED NODES JOHN GUCKENHEIMER AND RADU HAIDUC For Yulij Ilyashenko with admiration and affection on the occasion

More information

### Period-doubling cascades of a Silnikov equation Period-doubling cascades of a Silnikov equation Keying Guan and Beiye Feng Science College, Beijing Jiaotong University, Email: keying.guan@gmail.com Institute of Applied Mathematics, Academia Sinica,

More information

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

More information

### UNIVERSIDADE DE SÃO PAULO 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

More information

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

More information

### Melnikov Method for Autonomous Hamiltonians Contemporary Mathematics Volume 00, 19xx Melnikov Method for Autonomous Hamiltonians Clark Robinson Abstract. This paper presents the method of applying the Melnikov method to autonomous Hamiltonian systems

More information

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

More information

### Homoclinic Bifurcation in an SIQR Model for Childhood Diseases Journal of Differential Equations 168, 150167 (2000) doi:10.1006jdeq.2000.3882, available online at http:www.idealibrary.com on Homoclinic Bifurcation in an SIQR Model for Childhood Diseases Lih-Ing Wu

More information

### Connecting Orbits with Bifurcating End Points Connecting Orbits with Bifurcating End Points Thorsten Hüls Fakultät für Mathematik Universität Bielefeld Postfach 100131, 33501 Bielefeld Germany huels@mathematik.uni-bielefeld.de June 24, 2004 Abstract

More information

### NBA Lecture 1. Simplest bifurcations in n-dimensional ODEs. Yu.A. Kuznetsov (Utrecht University, NL) March 14, 2011 NBA Lecture 1 Simplest bifurcations in n-dimensional ODEs Yu.A. Kuznetsov (Utrecht University, NL) March 14, 2011 Contents 1. Solutions and orbits: equilibria cycles connecting orbits other invariant sets

More information

### Abstract. We consider a generalised Gause predator-prey system with a generalised Holling. ax 2. We study the cases where b is positive or negative. BIFURCATION ANALYSIS OF A PREDATOR-PREY SYSTEM WIT GENERALISED OLLING TYPE III FUNCTIONAL RESPONSE YANN LAMONTAGNE, CAROLINE COUTU, AND CRISTIANE ROUSSEAU Abstract. We consider a generalised Gause predator-prey

More information

### UNIFORM SUBHARMONIC ORBITS FOR SITNIKOV PROBLEM Manuscript submitted to Website: http://aimsciences.org AIMS Journals Volume 00, Number 0, Xxxx XXXX pp. 000 000 UNIFORM SUBHARMONIC ORBITS FOR SITNIKOV PROBLEM CLARK ROBINSON Abstract. We highlight the

More information

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

More information

### Invariant manifolds of the Bonhoeffer-van der Pol oscillator Invariant manifolds of the Bonhoeffer-van der Pol oscillator R. Benítez 1, V. J. Bolós 2 1 Dpto. Matemáticas, Centro Universitario de Plasencia, Universidad de Extremadura. Avda. Virgen del Puerto 2. 10600,

More information

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

More information

### The Hopf-van der Pol System: Failure of a Homotopy Method DOI.7/s259--9-5 ORIGINAL RESEARCH The Hopf-van der Pol System: Failure of a Homotopy Method H. G. E. Meijer T. Kalmár-Nagy Foundation for Scientific Research and Technological Innovation 2 Abstract The

More information

### Example of a Blue Sky Catastrophe PUB:[SXG.TEMP]TRANS2913EL.PS 16-OCT-2001 11:08:53.21 SXG Page: 99 (1) Amer. Math. Soc. Transl. (2) Vol. 200, 2000 Example of a Blue Sky Catastrophe Nikolaĭ Gavrilov and Andrey Shilnikov To the memory of

More information

### M3-4-5 A16 Notes for Geometric Mechanics: Oct Nov 2011 M3-4-5 A16 Notes for Geometric Mechanics: Oct Nov 2011 Text for the course: Professor Darryl D Holm 25 October 2011 Imperial College London d.holm@ic.ac.uk http://www.ma.ic.ac.uk/~dholm/ Geometric Mechanics

More information

### Survey of strong normal-internal k : l resonances in quasi-periodically driven oscillators for l = 1, 2, 3. June, : WSPC - Proceedings Trim Size: in x in SPT-broer Survey of strong normal-internal k : l resonances in quasi-periodically driven oscillators for l =,,. H.W. BROER and R. VAN DIJK Institute for mathematics

More information

### DYNAMICAL SYSTEMS. I Clark: Robinson. Stability, Symbolic Dynamics, and Chaos. CRC Press Boca Raton Ann Arbor London Tokyo DYNAMICAL SYSTEMS Stability, Symbolic Dynamics, and Chaos I Clark: Robinson CRC Press Boca Raton Ann Arbor London Tokyo Contents Chapter I. Introduction 1 1.1 Population Growth Models, One Population 2

More information

### BIFURCATION DIAGRAM OF A CUBIC THREE-PARAMETER AUTONOMOUS SYSTEM Electronic Journal of Differential Equations, Vol. 2005(2005), No. 8, pp. 1 16. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) BIFURCATION

More information

### arxiv: v1 [nlin.cd] 20 Jul 2010 Invariant manifolds of the Bonhoeffer-van der Pol oscillator arxiv:1007.3375v1 [nlin.cd] 20 Jul 2010 R. Benítez 1, V. J. Bolós 2 1 Departamento de Matemáticas, Centro Universitario de Plasencia, Universidad

More information

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

More information

### Homoclinic saddle to saddle-focus transitions in 4D systems Faculty of Electrical Engineering, Mathematics & Computer Science Homoclinic saddle to saddle-focus transitions in 4D systems Manu Kalia M.Sc. Thesis July 2017 Assessment committee: Prof. Dr. S. A. van

More information

### Dynamical Systems. Pierre N.V. Tu. An Introduction with Applications in Economics and Biology Second Revised and Enlarged Edition. Pierre N.V. Tu Dynamical Systems An Introduction with Applications in Economics and Biology Second Revised and Enlarged Edition With 105 Figures Springer-Verlag Berlin Heidelberg New York London Paris

More information

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

More information

### PIECEWISE LINEAR DIFFERENTIAL SYSTEMS WITH TWO REAL SADDLES. 1. Introduction PIECEWISE LINEAR DIFFERENTIAL SYSTEMS WITH TWO REAL SADDLES JOAN C. ARTÉS1, JAUME LLIBRE 1, JOAO C. MEDRADO 2 AND MARCO A. TEIXEIRA 3 Abstract. In this paper we study piecewise linear differential systems

More information

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

More information

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

More information

### Towards a Global Theory of Singularly Perturbed Dynamical Systems John Guckenheimer Cornell University Towards a Global Theory of Singularly Perturbed Dynamical Systems John Guckenheimer Cornell University Dynamical systems with multiple time scales arise naturally in many domains. Models of neural systems

More information

### The Existence of Chaos in the Lorenz System The Existence of Chaos in the Lorenz System Sheldon E. Newhouse Mathematics Department Michigan State University E. Lansing, MI 48864 joint with M. Berz, K. Makino, A. Wittig Physics, MSU Y. Zou, Math,

More information

### HOPF BIFURCATION AND HETEROCLINIC CYCLES IN A CLASS OF D 2 EQUIVARIANT SYSTEMS HOPF BIFURCATION AND HETEROCLINIC CYCLES IN A CLASS OF D 2 EQUIVARIANT SYSTEMS ADRIAN C. MURZA Communicated by Vasile Brînzănescu In this paper, we analyze a generic dynamical system with D 2 constructed

More information

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

### The reduction to invariant cones for nonsmooth systems. Tassilo Küpper Mathematical Institute University of Cologne The reduction to invariant cones for nonsmooth systems Tassilo Küpper Mathematical Institute University of Cologne kuepper@math.uni-koeln.de International Workshop on Resonance Oscillations and Stability

More information

### Perturbation Theory of Dynamical Systems Perturbation Theory of Dynamical Systems arxiv:math/0111178v1 [math.ho] 15 Nov 2001 Nils Berglund Department of Mathematics ETH Zürich 8092 Zürich Switzerland Lecture Notes Summer Semester 2001 Version:

More information

### Hopf bifurcation in coupled cell networks with abelian symmetry 1 Hopf bifurcation in coupled cell networks with abelian symmetry 1 Ana Paula S. Dias e Rui C. Paiva CMUP Departamento de Matemática Universidade do Porto CMUP Instituto Politécnico de Leiria Abstract We

More information

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

More information

### LIE TRANSFORM CONSTRUCTION OF EXPONENTIAL NORMAL FORM OBSERVERS. Suba Thomas Harry G. Kwatny Bor-Chin Chang LIE TRANSFORM CONSTRUCTION OF EXPONENTIAL NORMAL FORM OBSERVERS Suba Thomas Harry G Kwatny Bor-Chin Chang Department of Mechanical Engineering and Mechanics, Drexel University, Philadelphia, PA 1914, USA

More information

### Improved Homoclinic Predictor for Bogdanov Takens Bifurcation International Journal of Bifurcation and Chaos, Vol. 24, No. 4 (214) 14557 (12 pages) c World Scientific Publishing Company DOI: 1.1142/S2181274145576 Improved Homoclinic Predictor for Bogdanov Takens

More information

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

More information

### THE ROSENZWEIG-MACARTHUR PREDATOR-PREY MODEL THE ROSENZWEIG-MACARTHUR PREDATOR-PREY MODEL HAL L. SMITH* SCHOOL OF MATHEMATICAL AND STATISTICAL SCIENCES ARIZONA STATE UNIVERSITY TEMPE, AZ, USA 8587 Abstract. This is intended as lecture notes for nd

More information

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

More information

### DIFFERENTIAL EQUATIONS, DYNAMICAL SYSTEMS, AND AN INTRODUCTION TO CHAOS DIFFERENTIAL EQUATIONS, DYNAMICAL SYSTEMS, AND AN INTRODUCTION TO CHAOS Morris W. Hirsch University of California, Berkeley Stephen Smale University of California, Berkeley Robert L. Devaney Boston University

More information

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

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

More information

### Computational Methods in Dynamical Systems and Advanced Examples and Advanced Examples Obverse and reverse of the same coin [head and tails] Jorge Galán Vioque and Emilio Freire Macías Universidad de Sevilla July 2015 Outline Lecture 1. Simulation vs Continuation. How

More information

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

More information

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

More information

### ONE-PARAMETER BIFURCATIONS IN PLANAR FILIPPOV SYSTEMS International Journal of Bifurcation and Chaos, Vol. 3, No. 8 (23) 257 288 c World Scientific Publishing Company ONE-PARAMETER BIFURCATIONS IN PLANAR FILIPPOV SYSTEMS YU. A. KUZNETSOV, S. RINALDI and A.

More information

### THE HARTMAN-GROBMAN THEOREM AND THE EQUIVALENCE OF LINEAR SYSTEMS THE HARTMAN-GROBMAN THEOREM AND THE EQUIVALENCE OF LINEAR SYSTEMS GUILLAUME LAJOIE Contents 1. Introduction 2 2. The Hartman-Grobman Theorem 2 2.1. Preliminaries 2 2.2. The discrete-time Case 4 2.3. The

More information

### Stability Analysis of Uzawa-Lucas Endogenous Growth Model Abstract: Stability Analysis of Uzawa-Lucas Endogenous Growth Model William A. Barnett* University of Kansas, Lawrence and Center for Financial Stability, NY City and Taniya Ghosh Indira Gandhi Institute

More information

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

More information

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

More information

### VERSAL DEFORMATIONS OF BILINEAR SYSTEMS UNDER OUTPUT-INJECTION EQUIVALENCE PHYSCON 2013 San Luis Potosí México August 26 August 29 2013 VERSAL DEFORMATIONS OF BILINEAR SYSTEMS UNDER OUTPUT-INJECTION EQUIVALENCE M Isabel García-Planas Departamento de Matemàtica Aplicada I Universitat

More information

### Here f and f2 are holomorphic functions of the state variables x; _x; y; _y and parameters 2 R k. These functions are \higher order" in the sense that Interactions of Andronov-Hopf and Bogdanov-Takens Bifurcations W. F. Langford and K. Zhan Department of Mathematics and Statistics University of Guelph Guelph, Ontario, Canada NG 2W June 5, 998 Abstract

More information

### Global bifurcation diagram for the Kerner-Konhäuser traffic flow model International Journal of Bifurcation and Chaos c World Scientific Publishing Company Global bifurcation diagram for the Kerner-Konhäuser traffic flow model Joaquín Delgado Mathematics Department. UAM Iztapalapa.

More information

### Stability Computations for Nilpotent Hopf Bifurcations in Coupled Cell Systems Stability Computations for Nilpotent Hopf Bifurcations in Coupled Cell Systems M. Golubitsky M. Krupa April 18, 006 Abstract Vanderbauwhede and van Gils, Krupa, and Langford studied unfoldings of bifurcations

More information

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

More information

### Global bifurcation diagram for the Kerner-Konhäuser traffic flow model International Journal of Bifurcation and Chaos c World Scientific Publishing Company Global bifurcation diagram for the Kerner-Konhäuser traffic flow model Joaquín Delgado Mathematics Department. UAM Iztapalapa.

More information

### Determining the Global Topology of Resonance Surfaces for Periodically Forced Oscillator Families Determining the Global Topology of Resonance Surfaces for Periodically Forced Oscillator Families Richard P. McGehee School of Mathematics University of Minnesota Minneapolis, MN 55455 Bruce B. Peckham

More information

### Preprint Preprint Preprint Preprint CADERNOS DE MATEMÁTICA 15, 179 193 October 2014) ARTIGO NÚMERO SMA#398 On the integrability and the zero Hopf bifurcation of a Chen Wang differential system Jaume Llibre Departament de Matemàtiques, Universitat

More information

### Spike-adding canard explosion of bursting oscillations Spike-adding canard explosion of bursting oscillations Paul Carter Mathematical Institute Leiden University Abstract This paper examines a spike-adding bifurcation phenomenon whereby small amplitude canard

More information

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

More information

### Defining Equations for Bifurcations and Singularities Defining Equations for Bifurcations and Singularities John Guckenheimer and Yongjian Xiang Mathematics Department, Ithaca, NY 14853 For Vladimir Arnold on the occasion of his 65th birthday July 1, 2002

More information

### THE CUSP HOPF BIFURCATION * International Journal of Bifurcation and Chaos, Vol. 17, No. 8 (2007) 2547 2570 c World Scientific Publishing Company THE CUSP HOPF BIFURCATION * J. HARLIM Courant Institute of Mathematical Sciences, New

More information

### PERIODIC SOLUTIONS WITH NONCONSTANT SIGN IN ABEL EQUATIONS OF THE SECOND KIND PERIODIC SOLUTIONS WITH NONCONSTANT SIGN IN ABEL EQUATIONS OF THE SECOND KIND JOSEP M. OLM, XAVIER ROS-OTON, AND TERE M. SEARA Abstract. The study of periodic solutions with constant sign in the Abel equation

More information

### An Undergraduate s Guide to the Hartman-Grobman and Poincaré-Bendixon Theorems An Undergraduate s Guide to the Hartman-Grobman and Poincaré-Bendixon Theorems Scott Zimmerman MATH181HM: Dynamical Systems Spring 2008 1 Introduction The Hartman-Grobman and Poincaré-Bendixon Theorems

More information