MA371 the Qualitative Theory of Ordinary Differential Equations

Total Page:16

File Type:pdf, Size:1020Kb

MA371 the Qualitative Theory of Ordinary Differential Equations MA371 The Qualitative Theory of Ordinary Differential Equations As Lectured by Professor C. Sparrow Typeset by D. Kitson Assisted by K. Crooks October 2009 to Jan 2010 Contents 1 Introduction 2 2 One Dimensional ODES (X = R1) 5 2.1 Flows, Existence and Uniqueness . 5 2.2 Orbits . 7 2.3 Stability . 9 3 Two Dimensional Flows (X = R2) 11 3.1 Hamiltonian Flows and Hamiltonians . 11 3.2 Liaponov Functions . 13 3.3 Fixed Points and Nearby Behaviour . 15 3.3.1 The Linear Part . 15 3.3.2 The Non-Linear Part . 16 3.4 Stable and Unstable Manifolds . 16 3.5 Stable Manifold Theorem (SMT) . 18 3.6 Hartman-Grobman Theorem . 20 3.7 Non-Hyperbolic Fixed points . 20 4 Periodic Orbits 23 4.1 The Poincar´e-BendixsonTheorem . 23 4.2 Dulac's Criterion or the Negative Divergence Test . 24 4.3 Poincar´eIndex in R2 ................................... 26 4.4 Stability of Periodic Orbits . 27 4.5 Van de Pol Oscillators . 28 5 Bifurcations 34 5.1 Bifurcations in R1 .................................... 34 5.1.1 Summary . 37 5.2 Bifurcations in R2 and the Central Manifold Theorem . 37 5.3 Hopf Bifurcations . 40 5.4 Co-Dimension 1 Bifurcations of Periodic Orbits . 41 5.5 Global Bifurcations . 42 5.5.1 Saddle Node on a cycle . 43 5.5.2 Homoclinic Orbits . 44 Warning: I cannot guarantee that these notes are accurate or complete. Note: Definitions will be in Bold Text. Chapter 1 Introduction We will study Ordinary Differential Equations, or ODEs, of the formx _ = f(x). The State Variables x vary in the the State Space (or phase space) X. In this course X will be finite dimensional (Partial Differential Equations live in infinite dimensional spaces). X is a metric space, usually Rn. x(t) is the State of the system x 2 X at time t. t 2 R is the Independent Variable. f : X ! X is a vector field. f(x) gives the speed and direction of the motion at x 2 X. The graph of (x(t); t) in X is the Phase Portrait. The graph of x(t) is called the Solution, the trajectory, the orbit of the system or the flow.1 Figure 1.1: A phase portrait. x_ = f(x) gives an autonomous DE, i.e. one not depending on t. This means that given any point (x1; : : : ; xn) 2 X, f(x1; : : : ; xn) will take the same value 8 t 2 R. This course is concerned with the study of first order autonomous ODEs. Qualitative Theory helps one to understand the local and global behaviour of an ODE without actually having to find explicit solutions to them; thus the subject is of great importance as many ODEs have no explicit solutions. t Example 1.0.1: Take X = R,_x = x has solution x(t) = x0e , with x(0) = x0. Looking at the phase portrait, one can see that 0 is an unstable fixed point, which is also a global repeller and all orbits not starting at 0 go to ±∞. This is the qualitative description of the the behaviour of the system. Thus, if we know f or its graph, we don't need a solution to understand the behaviour of the system. This works for lots of \nice" fs in R, but this becomes very difficult in Rn≥3. Example 1.0.2: The below cubic-like function f displays the behaviour one would expect as the ODEx _ = f(x), and is similar to the type of thing studied in first year ODEs. 2 1 Example 1.0.3: Considerx _ = f(x) = Sign(x)x . f(x) = 0 at x = 0. For x > 0, x(t) = −1 . x0 −t The solution goes to infinity as t ! 1 . This is called Finite Time Blowup (FTB), i.e. x ! 1 as x0 t ! T < 1. This happens but we do not worry about it. It doesn't change the phase portrait. However, we do worry about the existence and uniqueness of solutions. 1From now on the vector x will not be underlined unless we have something like x = (x; y). 2 Figure 1.2: Two \nice" functions giving the same result, but only one has an explicit solution. Figure 1.3: Example 1.0.2: a cubic-like function giving two \unstable" fixed points and one \stable" fixed point, as per MA133 Differential Equations. 1 Example 1.0.4: Considerx _ = jxj 2 Sign(x). We find that for any T ≥ 0, we have a solution x = 0 t 2 [0;T ] x(t) = 1 2 x = 4 (t − T ) t 2 [T; 1) so the solution is not unique - what about 0 and T ? Example 1.0.5: With X = R2. Consider 1-dimensional Simple Harmonic Motion (SHM), given byx ¨ = −x. Write y =x _ to gety _ = −x,_x = y, which is a first order ODE in R2. We can write this in the form x_ 0 1 x = y_ −1 0 y We can find solutions of the form x(t) = r0 cos(t + φ); y(t) = r0 sin(t + φ) The phase portrait shows that all solutions are periodic orbits around a fixed point (0; 0), which is structurally unstable. Aside: One method of solving these is a change of coordinates. In X = R2 wherex _ = Ax with complex eigenvalues for A, write x = r cos(θ), y = r sin(θ). Then try to get to the formr _ = 0, θ_ = 1. Example 1.0.6: In X = R2, x_ 2 0 x = y_ 0 −1 y gives x ! ±∞ as t ! 1 and y ! 0 as t ! 1. So (0; 0) is a fixed point. but suppose we had some non-linear, i.e. higher order terms as well? Then the equations would look like x_ 2 0 x = + h: o: t: y_ 0 −1 y This usually can't be solved, but the qualitative question is \what does the phase portrait look like near (0; 0) now?". We will show that it is the same under certain conditions. 3 Figure 1.4: Example 1.0.5: the phase portrait of Simple Harmonic Motion. Example 1.0.7: With X = R3. x_ = 10(y − x) y_ = 28x − y − xz 8 z_ = z − xy 3 The Lorenz equations have only two quadratic but terms, but the phase portrait is chaotic with a strange attractor: Figure 1.5: The well-known Lorenz attractor. We are interested in behaviour (such as periodicity, going to a limit or infinity, etc) which is invariant under \nice" changes of coordinates. By a \nice" change of coordinates h with y = h(x) we mean a diffeomorphism, i.e. a differentiable bijection with differentiable inverse. We are also interested in properties invariant under parameterisation of time, e.g. new time s = α(t), where α : R ! R. Figure 1.6: These systems have been subjected to a diffeomorphism, but their phase portraits display the same sort of behaviors. What does it mean to say that there does not exist explicit solutions? Considery _ = y2 − t - this has no explicit solution, but we could define a new function L(y0; t) as being the solution, but obviously this is not helpful. Similarly the pendulum solution θ¨ = − sin(θ) can be solved using elliptic functions, but these are not \normal". 4 Chapter 2 1 One Dimensional ODES (X = R ) 2.1 Flows, Existence and Uniqueness It is more natural to write \solutions" as functions of starting position and time. Definition: φ : X × R ! R is called the Flow; φ(x0; t) = x(t) the place one gets to starting at x0 at t = 0 and solving until time t. We need φ to have some obvious properties: (I) φ(x0; 0) = x0 (II) φ(φ(x0; s); t) = φ(x0; s + t) A function φ satisfying these properties is a candidate for the solution of the ODE.1 To solve x_ = f(x), then it should satisfy d (φ(x; t)) = f(φ(x; t)) dt Equivalently, integrate both sides to get Z t φ(x0; t) = x0 + f (φ(x0; s)) ds 0 So, given f, does φ exist and is it unique? Definition: A function f : X ! Y metric spaces is Lipschitz with Lipschitz constant L if kf(x) − f(y)k ≤ Lkx − yk for some L 2 R, and for all x; y 2 X. Example 2.1.1: If f : R ! R, f is Lipschitz with constant L, the function remains in the region bounded by lines of slope ±L through the point (x; f(x)). Figure 2.1: Example 2.1.1: a Lipschitz function. f cannot escape the lines of slope ±L. 1If it does so for s; t ≥ 0 then it is called a semi flow in some books. 5 Note: Lipschitz ) Continuous, but Lipschitz ; Differentiable. Definition: f is Locally Lipschitz at the point x∗ with constant L if kf(x∗) − f(y)k ≤ Lkx∗ − yk for all y in a neighborhood of x∗. 1 Example 2.1.2: f(x) = jxj 2 is not locally Lipschitz at 0 because L would have to be ±∞. Outline of what we will do next: (1) Theorem: f Lipschitz ) 9! φ solvingx _ = f(x). (2) Worry: How long does φ exist for? The theorem is only local. (3) Proof: An algorithm to find a sequence of functions converging to φ. How to solvex _ = f(x)?2 We shall try to approximate the solution with a sequence of functions fuig: u0(x0; t) = x0 Z t u1(x0; t) = x0 + f (u0(x0; s)) ds 0 .
Recommended publications
  • Dynamics Forced by Homoclinic Orbits
    DYNAMICS FORCED BY HOMOCLINIC ORBITS VALENT´IN MENDOZA Abstract. The complexity of a dynamical system exhibiting a homoclinic orbit is given by the orbits that it forces. In this work we present a method, based in pruning theory, to determine the dynamical core of a homoclinic orbit of a Smale diffeomorphism on the 2-disk. Due to Cantwell and Conlon, this set is uniquely determined in the isotopy class of the orbit, up a topological conjugacy, so it contains the dynamics forced by the homoclinic orbit. Moreover we apply the method for finding the orbits forced by certain infinite families of homoclinic horseshoe orbits and propose its generalization to an arbitrary Smale map. 1. Introduction Since the Poincar´e'sdiscovery of homoclinic orbits, it is known that dynamical systems with one of these orbits have a very complex behaviour. Such a feature was explained by Smale in terms of his celebrated horseshoe map [40]; more precisely, if a surface diffeomorphism f has a homoclinic point then, there exists an invariant set Λ where a power of f is conjugated to the shift σ defined on the compact space Σ2 = f0; 1gZ of symbol sequences. To understand how complex is a diffeomorphism having a periodic or homoclinic orbit, we need the notion of forcing. First we will define it for periodic orbits. 1.1. Braid types of periodic orbits. Let P and Q be two periodic orbits of homeomorphisms f and g of the closed disk D2, respectively. We say that (P; f) and (Q; g) are equivalent if there is an orientation-preserving homeomorphism h : D2 ! D2 with h(P ) = Q such that f is isotopic to h−1 ◦ g ◦ h relative to P .
    [Show full text]
  • An Example of Bifurcation to Homoclinic Orbits
    JOURNAL OF DIFFERENTIAL EQUATIONS 37, 351-373 (1980) An Example of Bifurcation to Homoclinic Orbits SHUI-NEE CHOW* Department of Mathematics, Michigan State University, East Lansing, Michigan 48823 AND JACK K. HALE+ AND JOHN MALLET-PARET Lefschetz Center for Dynamical Systems, Division of Applied Mathematics, Brown University, Providence, Rhode Island 02912 Received October 19, 1979; revised January 23, 1980 Consider the equation f - x + x 2 = -A12 + &f(t) where f(t + 1) = f(t) and h = (h, , A,) is small. For h = 0, there is a homoclinic orbit r through zero. For X # 0 and small, there can be “strange” attractors near r. The pur- pose of this paper is to determine the curves in /\-space of bifurcation to “strange” attractors and to relate this to hyperbolic subharmonic bifurcations. I. INTRODUCTION At the present time, there is considerable interest in the theory of “strange” attractors of mappings and the manner in which strange attractors arise through successive bifurcations of periodic orbits (see [l, 5-121). The simplest “strange” attractors arise near homoclinic points. For a hyperbolic fixed point p of a diffeomorphism, a point 4 # p is a homo- clinic point of p if 4 is in the intersection of the stable and unstable manifolds of p. The point 4 is called a transverse homoclinic point of p if the intersection of the stable and unstable manifolds of p at 4 is transversal. It was known to Poincare that the existence of one transverse homoclinic point implies there must be infinitely many transverse homoclinic points. Birkhoff showed that every * This research was supported in part by the Netional Science Foundation, under MCS76-06739.
    [Show full text]
  • Attractors and Orbit-Flip Homoclinic Orbits for Star Flows
    PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 141, Number 8, August 2013, Pages 2783–2791 S 0002-9939(2013)11535-2 Article electronically published on April 12, 2013 ATTRACTORS AND ORBIT-FLIP HOMOCLINIC ORBITS FOR STAR FLOWS C. A. MORALES (Communicated by Bryna Kra) Abstract. We study star flows on closed 3-manifolds and prove that they either have a finite number of attractors or can be C1 approximated by vector fields with orbit-flip homoclinic orbits. 1. Introduction The notion of attractor deserves a fundamental place in the modern theory of dynamical systems. This assertion, supported by the nowadays classical theory of turbulence [27], is enlightened by the recent Palis conjecture [24] about the abun- dance of dynamical systems with finitely many attractors absorbing most positive trajectories. If confirmed, such a conjecture means the understanding of a great part of dynamical systems in the sense of their long-term behaviour. Here we attack a problem which goes straight to the Palis conjecture: The fini- tude of the number of attractors for a given dynamical system. Such a problem has been solved positively under certain circunstances. For instance, we have the work by Lopes [16], who, based upon early works by Ma˜n´e [18] and extending pre- vious ones by Liao [15] and Pliss [25], studied the structure of the C1 structural stable diffeomorphisms and proved the finitude of attractors for such diffeomor- phisms. His work was largely extended by Ma˜n´e himself in the celebrated solution of the C1 stability conjecture [17]. On the other hand, the Japanesse researchers S.
    [Show full text]
  • Codimension-Two Homoclinic Bifurcations Underlying Spike Adding in the Hindmarsh-Rose Burster Arxiv:1109.5689V1 [Math.DS] 26 S
    Codimension-two homoclinic bifurcations underlying spike adding in the Hindmarsh-Rose burster Daniele Linaro1, Alan Champneys2, Mathieu Desroches2, Marco Storace1 1Biophysical and Electronic Engineering Department, University of Genoa, Genova, Italy 2Department of Engineering Mathematics, University of Bristol, Bristol, UK July 2, 2018 Abstract The well-studied Hindmarsh-Rose model of neural action potential is revisited from the point of view of global bifurcation analysis. This slow-fast system of three paremeterised differential equations is arguably the simplest reduction of Hodgkin- Huxley models capable of exhibiting all qualitatively important distinct kinds of spiking and bursting behaviour. First, keeping the singular perturbation parameter fixed, a comprehensive two-parameter bifurcation diagram is computed by brute force. Of particular concern is the parameter regime where lobe-shaped regions of irregular bursting undergo a transition to stripe-shaped regions of periodic bursting. The boundary of each stripe represents a fold bifurcation that causes a smooth spike-adding transition where the number of spikes in each burst is increased by one. Next, numerical continuation studies reveal that the global structure is organised by various curves of homoclinic bifurcations. In particular the lobe to stripe transition is organised by a sequence of codimension- two orbit- and inclination-flip points that occur along each homoclinic branch. Each branch undergoes a sharp turning point and hence approximately has a double- cover of the same curve in parameter space. The sharp turn is explained in terms of the interaction between a two-dimensional unstable manifold and a one-dimensional slow manifold in the singular limit. Finally, a new local analysis is undertaken us- ing approximate Poincar´emaps to show that the turning point on each homoclinic branch in turn induces an inclination flip that gives birth to the fold curve that or- ganises the spike-adding transition.
    [Show full text]
  • Homoclinic and Heteroclinic Bifurcations in Vector Fields
    Homoclinic and Heteroclinic Bifurcations in Vector Fields Ale Jan Homburg Bj¨ornSandstede Korteweg{de Vries Institute for Mathematics Division of Applied Mathematics University of Amsterdam Brown University 1018 TV Amsterdam, The Netherlands Providence, RI 02906, USA June 17, 2010 Abstract An overview of homoclinic and heteroclinic bifurcation theory for autonomous vector fields is given. Specifically, homoclinic and heteroclinic bifurcations of codimension one and two in generic, equivariant, reversible, and conservative systems are reviewed, and results pertaining to the existence of multi-round homoclinic and periodic orbits and of complicated dynamics such as suspended horseshoes and attractors are stated. Bifurcations of homoclinic orbits from equilibria in local bifurcations are also considered. The main analytic and geometric techniques such as Lin's method, Shil'nikov variables and homoclinic center manifolds for analyzing these bifurcations are discussed. Finally, a few related topics, such as topological moduli, numerical algorithms, variational methods, and extensions to singularly perturbed and infinite- dimensional systems, are reviewed briefly. 1 Contents 1 Introduction 4 2 Homoclinic and heteroclinic orbits, and their geometry5 2.1 Homoclinic orbits to hyperbolic equilibria.............................6 2.2 Homoclinic orbits to nonhyperbolic equilibria........................... 10 2.3 Heteroclinic cycles with hyperbolic equilibria........................... 12 3 Analytical and geometric approaches 13 3.1 Normal forms and linearizability.................................
    [Show full text]
  • Homoclinic Tangles for Noninvertible Maps
    Nonlinear Analysis 41 (2000) 259–276 www.elsevier.nl/locate/na Homoclinic tangles for noninvertible maps Evelyn Sander Department of Mathematical Sciences, George Mason University, MS 3F2, 4400 University Drive, Fairfax, VA 22030, USA Received 2 March 1998; accepted 17 July 1998 Keywords: Homoclinic tangles; Noninvertible maps; Transverse homoclinic point 1. Introduction Not all naturally occurring iterated systems have the property that time is reversible. In which case, the dynamics must be modelled by a noninvertible map rather than by a di eomorphism. While one-dimensional theory describes dynamics of noninvertible maps, higher-dimensional theory has focussed on di eomorphisms, though there are many higher-dimensional examples of natural systems modelled by noninvertible maps. For example, noninvertible maps occur in the study of population dynamics [2], time one maps of delay equations [25], control theory algorithms [1, 5, 6], neural networks [16], and iterated di erence methods [14]. In order to use dynamical systems theory in applications for which maps are noninvertible, it is important to know how the theory di ers from the di eomorphism case. One of the principal indications of chaotic dynamics in deterministic systems is the existence of transverse homoclinic orbits. This paper examines the existence of homoclinic tangles and the associated dynamic behavior for noninvertible maps. In particular, we give examples of maps which have transverse intersection of stable and unstable manifolds but no horseshoe. Furthermore, we show that this is a codimension two phenomenon. Papers by Steinlein and Walther [24,25] and Lerman and Shil’nikov [13] both estab- lish conditions for a transverse homoclinic point of a noninvertible map to be contained in a horseshoe.
    [Show full text]
  • Homoclinic Orbits in Hamiltonian Systems
    CORE Metadata, citation and similar papers at core.ac.uk Provided by Elsevier - Publisher Connector JOURNAL OF DIFFERENTIAL EQUATIONS 21,431-438 (1976) Homoclinic Orbits in Hamiltonian Systems ROBERTL. DEVANEY Department of Mathematics, Northwestern University, Evanston, Illinois 60201 Received January 25, 1975 The object of this paper is to study the orbit structure of a Hamiltonian system in a neighborhood of a trajectory which is doubly asymptotic to an equilibrium solution, i.e., an orbit which lies in the intersection of the stable and unstable manifolds of a critical point. Such an orbit is called a homoclinic orbit. For diffeomorphisms, the analogous situation is fairly well understood. By a theorem of Smale [7], in every neighborhood of a transversal homoclinic point of a periodic point, there is a compact invariant set on which some iterate of the diffeomorphism is topologically conjugate to the Bernoulli shift on N symbols. Via Poincare maps on local transversal sections, there is thus an analogous result for hyperbolic closed orbits of vector fields. For critical points of vector fields, however, the situation is somewhat different. In the first place, by the Kupka-Smale theorem [4], the existence of orbits doubly asymptotic to an equilibrium point is not generic. In fact, the set of vector fields which admit homoclinic orbits at critical points is of the first Baire category in the set of all smooth vector fields. This follows since the sum of the dimensions of the stable and unstable manifolds at the critical point is at most equal to the dimension of the mani- fold itself.
    [Show full text]
  • Heteroclinic Connections Between Periodic Orbits and Resonance Transitions in Celestial Mechanics
    Heteroclinic Connections between Periodic Orbits and Resonance Transitions in Celestial Mechanics Wang Sang Koon Control and Dynamical Systems and JPL Caltech 107-81, Pasadena, CA 91125 [email protected] Martin W. Lo Navigation and Flight Mechanics Jet Propulsion Laboratory M/S: 301-142 4800 Oak Grove Drive Pasadena, CA 91109-8099 [email protected] Jerrold E. Marsden Control and Dynamical Systems Caltech 107-81, Pasadena, CA 91125 [email protected] Shane D. Ross Control and Dynamical Systems and JPL Caltech 107-81, Pasadena, CA 91125 [email protected] October, 1997. This version: April 12, 2000 Web site for the color version of this article: http://www.cds.caltech.edu/˜marsden/ Abstract This paper applies dynamical systems techniques to the problem of heteroclinic con- nections and resonance transitions in the planar circular restricted three-body problem. These related phenomena have been of concern for some time in topics such as the cap- ture of comets and asteroids and with the design of trajectories for space missions such as the Genesis Discovery Mission. The main new technical result in this paper is the numerical demonstration of the existence of a heteroclinic connection between pairs of periodic orbits, one around the libration point L1 and the other around L2, with the two periodic orbits having the same energy. This result is applied to the resonance transition problem and to the explicit numerical construction of interesting orbits with prescribed itineraries. The point of view developed in this paper is that the invariant manifold structures associated to L1 and L2 as well as the aforementioned heteroclinic connection are fundamental tools that can aid in understanding dynamical channels throughout the solar system as well as transport between the “interior” and “exterior” Hill’s regions and other resonant phenomena.
    [Show full text]
  • Dynamical Systems Theory
    Dynamical Systems Theory Bjorn¨ Birnir Center for Complex and Nonlinear Dynamics and Department of Mathematics University of California Santa Barbara1 1 c 2008, Bjorn¨ Birnir. All rights reserved. 2 Contents 1 Introduction 9 1.1 The 3 Body Problem . .9 1.2 Nonlinear Dynamical Systems Theory . 11 1.3 The Nonlinear Pendulum . 11 1.4 The Homoclinic Tangle . 18 2 Existence, Uniqueness and Invariance 25 2.1 The Picard Existence Theorem . 25 2.2 Global Solutions . 35 2.3 Lyapunov Stability . 39 2.4 Absorbing Sets, Omega-Limit Sets and Attractors . 42 3 The Geometry of Flows 51 3.1 Vector Fields and Flows . 51 3.2 The Tangent Space . 58 3.3 Flow Equivalance . 60 4 Invariant Manifolds 65 5 Chaotic Dynamics 75 5.1 Maps and Diffeomorphisms . 75 5.2 Classification of Flows and Maps . 81 5.3 Horseshoe Maps and Symbolic Dynamics . 84 5.4 The Smale-Birkhoff Homoclinic Theorem . 95 5.5 The Melnikov Method . 96 5.6 Transient Dynamics . 99 6 Center Manifolds 103 3 4 CONTENTS 7 Bifurcation Theory 109 7.1 Codimension One Bifurcations . 110 7.1.1 The Saddle-Node Bifurcation . 111 7.1.2 A Transcritical Bifurcation . 113 7.1.3 A Pitchfork Bifurcation . 115 7.2 The Poincare´ Map . 118 7.3 The Period Doubling Bifurcation . 119 7.4 The Hopf Bifurcation . 121 8 The Period Doubling Cascade 123 8.1 The Quadradic Map . 123 8.2 Scaling Behaviour . 130 8.2.1 The Singularly Supported Strange Attractor . 138 A The Homoclinic Orbits of the Pendulum 141 List of Figures 1.1 The rotation of the Sun and Jupiter in a plane around a common center of mass and the motion of a comet perpendicular to the plane.
    [Show full text]
  • Saddle-Center and Periodic Orbit: Dynamics Near Symmetric Heteroclinic Connection
    Saddle-center and periodic orbit: dynamics near symmetric heteroclinic connection L.M. Lermana;b, K.N. Trifonovb aNational Research University Higher School of Economics, bScientific and Educational Mathematical Center “Mathematics of Future Technologies” Lobachevsky State Research University of Nizhny Novgorod, Russia Abstract An analytic reversible Hamiltonian system with two degrees of freedom is studied in a neighborhood of its symmetric heteroclinic connection made up of a symmetric saddle-center, a symmetric orientable saddle periodic orbit lying in the same level of a Hamiltonian and two non-symmetric heteroclinic orbits permuted by the involution. This is a codimension one structure and therefore it can be met generally in one-parameter families of reversible Hamiltonian systems. There exist two possible types of such connections in dependence on how the involution acts near the equilibrium. We prove a series of theorems which show a chaotic behavior of the system and those in its unfoldings, in particular, the existence of countable sets of transverse homoclinic orbits to the saddle periodic orbit in the critical level, transverse heteroclinic connections involving a pair of saddle periodic orbits, families of elliptic periodic orbits, homoclinic tangencies, families of homoclinic orbits to saddle-centers in the unfolding, etc. As a byproduct, we get a criterion of the existence of homoclinic orbits to a saddle-center. The orbit structure of a non-integrable Hamiltonian system with two or more degrees of freedom is very complicated and usually it is impossible, except for some specific model situations, to describe its structure more or less completely. By this reason, a fruitful way to understand the orbit behavior in some parts of the phase space is to detect some simple invariant subsets (usually containing a finite number of orbits) whose neighborhoods can be understood from the viewpoint of their orbit structure.
    [Show full text]
  • Transversal Heteroclinic and Homoclinic Orbits in Singular Perturbation Problems*
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by Elsevier - Publisher Connector JOURNAL OF DIFPERENTIAL EQUATIONS 92, 252-281 (1991) Transversal Heteroclinic and Homoclinic Orbits in Singular Perturbation Problems* P. SZMOLYAN Insritiir ftir Angewandie und Numerische Mathematik, Technische Universitiii Wien, Vienna, Austria Communicated by Jack K. Hale Received December 4, 1989; revised March 13, 1990 In this paper we give a geometric construction of heteroclinic and homoclinic orbits for singularly perturbed differential equations. By using methods from invariant manifold theory we show that transversal intersection of stable and unstable manifolds of the reduced problem implies the existence of transversal heteroclinic or homoclinic orbits of the singularly perturbed problem. We derive analytical conditions for transversality. We show how these results can be used to prove the existence of heteroclinic and homoclinic orbits in singularly perturbed problems which depend on additional parameters. We describe a configuration which implies transversal intersection of the stable and unstable manifolds of periodic orbits and the associated chaotic dynamics. 0 1991 Academic Press, Inc. 1. INTROOUCTI~N We consider singularly perturbed systems of differential equations “t =f(x, Y, E) (1.1) d=gb,Y,E) with E E ( - sO,so), so > 0 small, and (x, y) E A, an open subset A? c R” + k. We assumethat (A g) E c’(.M x ( -so, sO)+ R” + “), with r suffkiently large. For E z 0 these equations define a smooth dynamical system on A. We consider x and y as functions of the variable t. Problems of the form (1.1) arise frequently in applications.
    [Show full text]
  • Chapter 23 Predicting Chaos
    Chapter 23 Predicting Chaos We have discussed methods for diagnosing chaos, but what about predicting the existence of chaos in a dynamical system. This is a much harder problem, and it seems that the best that can be done is to predict the existence of “chaotic sets” that are not attractors but may perhaps “organize” the chaotic attractor. These methods rely on the understanding of the complicated structure that emerges in phase space when a homoclinic or heteroclinic orbit is perturbed. To introduce these ideas we first introduce a definition of chaos that is useful for sets that are not attractors, motivating the idea using the simplest of chaotic maps, the shift map. 23.1 The Shift Map and Symbolic Dynamics The map xn+1 = 2xn mod 1 (23.1) displays many of the features of chaos in an analytically accessible way. The tent map at a = 2 is equivalent to the shift map. To understand the dynamics write the initial point in binary representation x0 = 0.a1a2 ... equivalent to X −ν x0 = av2 (23.2) with each ai either 0 or 1. The action of the map is simply to shift the “binal” point to the right and discard the first digit x1 = f(x0) = 0.a2a3 ... (23.3) 1 CHAPTER 23. PREDICTING CHAOS 2 Many properties are now evident: 1. Sensitivity to initial conditions: two initial conditions differing by of order 2−m lead to orbits that differ by O(1) after m iterations. 2. The answer to the question whether the mth iteration from a random initial 1 condition is greater or less than 2 is as random as a coin toss.
    [Show full text]