The Contopoulos and Hénon-Heiles model. Numerical integration. Numerical computation of the Poincaré section.
|
|
- Juliet Ryan
- 5 years ago
- Views:
Transcription
1 ½ ÈÇÁÆ Ê Ë ÌÁÇÆ The Contopoulos and Hénon-Heiles model. Numerical integration. Numerical computation of the Poincaré section.
2 Do or do not, there is no try, Yoda Operating system: GNU/Linux (I don t care about other OS, do you?). Programming language: C. Graphic library: Mizar Giorgilli (1979), fair enough! I will give you the code as is (GPLv3)... enjoy! REMARK: you are supposed to code! Be a man, write your own programs!, Carles Simó. A mathematical formula should never be owned by anybody! Mathematics belong to God., D.E. Knuth. Poincaré section 2
3 Yes, you need to code! Keep in mind: Real Programmers aren t afraid to use GOTOs. Real Programmers don t need comments the code is obvious. Real Programmers don t comment their code. If it was hard to write, it should be hard to understand. The computers in the Space Shuttle were programmed by Real Programmers. No Real Programmer works 9 to 5. (Unless it s the ones at night.) References: Real Programmers Don t Write Specs, Tom Van Vleck (1982). Real Programmers Don t Use Pascal, Ed Post, Tektronix, Wilsonville OR USA. (1983) Poincaré section 3
4 Tools needed C compiler (GCC). Make utility (GNU Make). Unix shell (Bash). Mizar. All the material for the exercises is on-line at the page ravello2014.html. I strongly suggest to use the live GNU/Linux Mint distribution, see info.txt. Poincaré section 4
5 1.1 The rediscovery of chaos: Contopoulos (1960), Hénon and Heiles (1964) H(q,p) = ω 0 2 (p2 0 +q2 0 )+ ω 1 2 (p2 1 +q2 1 )+q2 0 q q3 1. Equations of motion q = H p, ṗ = H q. Problem: study the dynamics of the system. Tool: numerical integrations and Poincaré sections. The idea is to transform the continuous flow induced of an autonomous differential equation in R n into a map Π : Σ Σ, with Σ R n a (n 1)-dimensional surface. Poincaré section 5
6 1.1.1 How to compute the Poincaré section H(q,p) = ω 0 2 (p2 0 +q2 0 )+ ω 1 2 (p2 1 +q2 1 )+q2 0 q q3 1, ω 0,ω 1 > 0. We need a surface transversal to the flow (q 0 = 0). Fix the initial energy E 0 and a point (q 1,p 1 ), and compute p 1 > 0 such that (0,q 1,p 0,p 1 ) Σ. The Poincaré section is the intersection of the orbit with the transversal surface. If E 0 satisfies the Poincaré section is bounded in the region { } ω E 0 E c, E c 3 = min ω2 0ω 1 8, ω3 1 6 ω 1 2 (p2 1 +q 2 1) 1 3 q3 1 E 0., We need an accurate numerical integrator. Random numbers should not be generated with a method chosen at random., D.E. Knuth. (it applies also to integration schemes...) Poincaré section 6
7 1.1.2 Runge-Kutta vs. Leap-Frog (order 2) Runge-Kutta (RK2): ẋ = f(x), dx n+ 1 2 = τ 2 f(x n), x n+1 = x n +τf(x n +dx n+ 1 2 ). Leapfrog (LF2): ẍ = f(x), x n+1 = x n +τẋ n+ 1 2, ẋ n+ 3 2 = ẋ n τf(x n+1). Question: How to choose the good one? Answer: The Hamiltonian flow is symplectic, the leapfrog is symplectic... wanna bet? An algorithm must be seen to be believed., D.E. Knuth. Poincaré section 7
8 1.1.3 Error against time-step RK2 ( ) and LF2 ( ) have the same accuracy. For τ too small, the accumulation of round-off errors is dominant Error against integration time RK2 ( ) shows a linear drift in the error. LF2 ( ) keeps its accuracy! Poincaré section 8
9 1.1.5 Open the box... (error against time-step) The file hh rk2 errstep.c contains the program that produces the data file needed to reproduce the plot error against time-step (RK2). Performs the following operations: read from the file hh rk2 errstep.inp some parameters; compute some rough Poincaré sections; compute the maximum relative error for different time-steps; store the results in hh rk2 errstep.dat The file hh lf2 errstep.c contains the program that produces the data file needed to reproduce the plot error against time-step (LF2). Performs the following operations: read from the file hh lf2 errstep.inp some parameters; compute some rough Poincaré sections; compute the maximum relative error for different time-steps; store the results in hh lf2 errstep.dat The file hh comp errstep.c permits to reproduce the plot error against time-step (RK2 vs. LF2). Poincaré section 9
10 1.1.6 Open the box... (error against integration time) The file hh rk2 errtot.c contains the program that produces the data file needed to reproduce the plot error against integration time (RK2). Performs the following operations: read from the file hh rk2 errtot.inp some parameters; compute some rough Poincaré sections; compute the maximum relative error for different integration times; store the results in hh rk2 errtot.dat The file hh lf2 errtot.c contains the program that produces the data file needed to reproduce the plot error against integration time (LF2). Performs the following operations: read from the file hh lf2 errtot.inp some parameters; compute some rough Poincaré sections; compute the maximum relative error for different integration times; store the results in hh lf2 errtot.dat The file hh comp errtot.c permits to reproduce the plot error against integration time (RK2 vs. LF2). Poincaré section 10
11 1.1.7 Poincaré sections for the Hénon-Heiles (resonant case) H(q,p) = ω 0 2 (p2 0 +q 2 0)+ ω 1 2 (p2 1 +q 2 1)+q 2 0q q3 1, ω 0 = ω 1 = 1. Energies: E = 1/100 (upper left); E = 1/12 (upper right); E = 1/8 (lower left); E = 1/6 (lower right). For low energy the dynamics appears to be completely ordered. The separation into two main regions is due to the 1:1 resonance. At the center of the main regions there are fixed points that correspond to periodic orbits; the surrounding curves represent beat phenomena due to resonance. For increasing energy a chaotic sea apparently filled by an orbit appears. Higher order resonances create ordered islands inside the chaotic sea. Poincaré section 11
12 1.1.8 Poincaré sections for the Hénon-Heiles (non-resonant case) H(q,p) = ω 0 2 (p2 0 +q2 0 )+ ω 1 2 (p2 1 +q2 1 )+q2 1 q q3 2, ω 0 = 1, ω 1 =. 2 E = E = E = E = E = E = Poincaré section 12
13 1.1.9 Poincaré sections for the Hénon-Heiles (frequencies with different sign) H(q,p) = ω 0 2 (p2 0 +q2 0 )+ ω 1 2 (p2 1 +q2 1 )+q2 1 q q3 2, ω 0 = 1, ω 1 =. 2 E = The energy surface is not compact. Many orbits escape to infinity. Successive enlargements of small regions (see the scale on the axes). The phenomenon of boxes into boxes : resonant islands and separatrices creating homoclinic intersections show up at every enlargement. Poincaré section 13
14 Open the box... (Poincaré sections) The file hh section nox.c allows to interactively set all the parameter of the Hénon-Heiles model and to compute the Poincaré section. The program store the points of the Poincaré section in the file hh section nox points.dat and the corresponding relative errors in hh section nox errors.dat. The accuracy of the Poincaré section can be improved changing a parameter in the code. The data file can be easily visualized using Gnuplot ( The file hh section X.c allows to interactively compute the Poincaré section for the Hénon-Heiles model. If the file hh section X.dat is found, it is supposed to be valid and the previous state is loaded. All the operations are stored in hh section X 2.dat. The idea is that after some runs, the good points are stored in a file that can be used as a skelton to produce nice plots. The program uses the Mizar graphics library in interactive mode via Tektronix Graphics (xterm). Active keys: E: change the energy F: change the frequencies R: reset Z: zoom (select a rectangle) B: change the boundary C + n: continue the previous orbit (n sections) <other> + n: compute n sections Poincaré section 14
15 1.2 Proposed exercises Write your own program to compute the Poincaré sections. If the graphics is a problem just save the date in a text file and use Gnuplot! Reproduce the figures reported in this chapter. Open the box and have a look at the code. Modify the programs that compute the Poincaré sections so as to read the input parameters from file. Improve the integration scheme (e.g., RK4, LF4, S BAB ). Science is what we understand well enough to explain to a computer. Art is everything else we do., D.E. Knuth. Poincaré section 15
16 Poincaré section 16
Painting chaos: OFLI 2 TT
Monografías de la Real Academia de Ciencias de Zaragoza. 28: 85 94, (26). Painting chaos: OFLI 2 TT R. Barrio Grupo de Mecánica Espacial. Dpto. Matemática Aplicada Universidad de Zaragoza. 59 Zaragoza.
More informationOn the numerical integration of variational equations
On the numerical integration of variational equations E. Gerlach 1, S. Eggl 2 and H. Skokos 3 1 Lohrmann Observatory, Technical University Dresden 2 Institute for Astronomy, University of Vienna 3 Max
More informationOn the numerical integration of variational equations
On the numerical integration of variational equations Haris Skokos Max Planck Institute for the Physics of Complex Systems Dresden, Germany E-mail: hskokos@pks.mpg.de, URL: http://www.pks.mpg.de/~hskokos/
More informationInfinity Unit 2: Chaos! Dynamical Systems
Infinity Unit 2: Chaos! Dynamical Systems Iterating Linear Functions These questions are about iterating f(x) = mx + b. Seed: x 1. Orbit: x 1, x 2, x 3, For each question, give examples and a symbolic
More informationLecture 1: A Preliminary to Nonlinear Dynamics and Chaos
Lecture 1: A Preliminary to Nonlinear Dynamics and Chaos Autonomous Systems A set of coupled autonomous 1st-order ODEs. Here "autonomous" means that the right hand side of the equations does not explicitly
More informationOne dimensional Maps
Chapter 4 One dimensional Maps The ordinary differential equation studied in chapters 1-3 provide a close link to actual physical systems it is easy to believe these equations provide at least an approximate
More informationvii Contents 7.5 Mathematica Commands in Text Format 7.6 Exercises
Preface 0. A Tutorial Introduction to Mathematica 0.1 A Quick Tour of Mathematica 0.2 Tutorial 1: The Basics (One Hour) 0.3 Tutorial 2: Plots and Differential Equations (One Hour) 0.4 Mathematica Programs
More information11 Chaos in Continuous Dynamical Systems.
11 CHAOS IN CONTINUOUS DYNAMICAL SYSTEMS. 47 11 Chaos in Continuous Dynamical Systems. Let s consider a system of differential equations given by where x(t) : R R and f : R R. ẋ = f(x), The linearization
More informationOn the numerical integration of variational equations
On the numerical integration of variational equations Haris Skokos Max Planck Institute for the Physics of Complex Systems Dresden, Germany E-mail: hskokos@pks.mpg.de, URL: http://www.pks.mpg.de/~hskokos/
More informationHamiltonian Dynamics
Hamiltonian Dynamics CDS 140b Joris Vankerschaver jv@caltech.edu CDS Feb. 10, 2009 Joris Vankerschaver (CDS) Hamiltonian Dynamics Feb. 10, 2009 1 / 31 Outline 1. Introductory concepts; 2. Poisson brackets;
More informationNumerical integration of variational equations
Numerical integration of variational equations Haris Skokos Max Planck Institute for the Physics of Complex Systems Dresden, Germany E-mail: hskokos@pks.mpg.de, URL: http://www.pks.mpg.de/~hskokos/ Enrico
More informationDynamics of Periodically Perturbed Homoclinic Solutions
Dynamics of Periodically Perturbed Homoclinic Solutions Qiudong Wang University of Arizona Part I Structure of Homoclinic Tangles Unperturbed equation Let (x,y) R 2 be the phase variables. dx dt = αx+f(x,y),
More informationChaos in open Hamiltonian systems
Chaos in open Hamiltonian systems Tamás Kovács 5th Austrian Hungarian Workshop in Vienna April 9. 2010 Tamás Kovács (MPI PKS, Dresden) Chaos in open Hamiltonian systems April 9. 2010 1 / 13 What s going
More informationB5.6 Nonlinear Systems
B5.6 Nonlinear Systems 5. Global Bifurcations, Homoclinic chaos, Melnikov s method Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Motivation 1.1 The problem 1.2 A
More informationProjects in Geometry for High School Students
Projects in Geometry for High School Students Goal: Our goal in more detail will be expressed on the next page. Our journey will force us to understand plane and three-dimensional geometry. We will take
More informationarxiv: v4 [nlin.cd] 12 Sep 2011
International Journal of Bifurcation and Chaos c World Scientific Publishing Company Chains of rotational tori and filamentary structures close to high multiplicity periodic orbits in a 3D galactic potential
More informationHÉNON HEILES HAMILTONIAN CHAOS IN 2 D
ABSTRACT HÉNON HEILES HAMILTONIAN CHAOS IN D MODELING CHAOS & COMPLEXITY 008 YOUVAL DAR PHYSICS DAR@PHYSICS.UCDAVIS.EDU Chaos in two degrees of freedom, demonstrated b using the Hénon Heiles Hamiltonian
More informationChaos in the Hénon-Heiles system
Chaos in the Hénon-Heiles system University of Karlstad Christian Emanuelsson Analytical Mechanics FYGC04 Abstract This paper briefly describes how the Hénon-Helies system exhibits chaos. First some subjects
More informationLogistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations
Logistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations S. Y. Ha and J. Park Department of Mathematical Sciences Seoul National University Sep 23, 2013 Contents 1 Logistic Map 2 Euler and
More informationChapter 23. Predicting Chaos The Shift Map and Symbolic Dynamics
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
More informationAMADEU DELSHAMS AND RAFAEL RAMíREZ-ROS
POINCARÉ-MELNIKOV-ARNOLD METHOD FOR TWIST MAPS AMADEU DELSHAMS AND RAFAEL RAMíREZ-ROS 1. Introduction A general theory for perturbations of an integrable planar map with a separatrix to a hyperbolic fixed
More informationMaps and differential equations
Maps and differential equations Marc R. Roussel November 8, 2005 Maps are algebraic rules for computing the next state of dynamical systems in discrete time. Differential equations and maps have a number
More informationChaos in the Planar Two-Body Coulomb Problem with a Uniform Magnetic Field
Annual Review of Chaos Theory, Bifurcations and Dynamical Systems Vol. 3, (2013) 23-33, www.arctbds.com. Copyright (c) 2013 (ARCTBDS). ISSN 2253 0371. All Rights Reserved. Chaos in the Planar Two-Body
More informationDerivation of Euler's Method - Numerical Methods for Solving Differential Equations
Derivation of Euler's Method - Numerical Methods for Solving Differential Equations Let s start with a general first order Initial Value Problem dx = f(x, y) y(x 0) = y 0 (1) where f(x, y) is a known function
More informationControl of a Chaotic Double Pendulum Using the OGY Method. Abstract
Control of a Chaotic Double Pendulum Using the OGY Method Matthew Cammack School of Physics Georgia Institute of Technology, Atlanta, GA 3033-0430, U.S.A (Dated: December 11, 003) Abstract This paper discusses
More informationSynchronization, Chaos, and the Dynamics of Coupled Oscillators. Supplemental 1. Winter Zachary Adams Undergraduate in Mathematics and Biology
Synchronization, Chaos, and the Dynamics of Coupled Oscillators Supplemental 1 Winter 2017 Zachary Adams Undergraduate in Mathematics and Biology Outline: The shift map is discussed, and a rigorous proof
More information1.7. Stability and attractors. Consider the autonomous differential equation. (7.1) ẋ = f(x),
1.7. Stability and attractors. Consider the autonomous differential equation (7.1) ẋ = f(x), where f C r (lr d, lr d ), r 1. For notation, for any x lr d, c lr, we let B(x, c) = { ξ lr d : ξ x < c }. Suppose
More informationAccel Alg E. L. E. Notes Solving Quadratic Equations. Warm-up
Accel Alg E. L. E. Notes Solving Quadratic Equations Warm-up Solve for x. Factor. 1. 12x 36 = 0 2. x 2 8x Factor. Factor. 3. 2x 2 + 5x 7 4. x 2 121 Solving Quadratic Equations Methods: (1. By Inspection)
More informationChaotic transport through the solar system
The Interplanetary Superhighway Chaotic transport through the solar system Richard Taylor rtaylor@tru.ca TRU Math Seminar, April 12, 2006 p. 1 The N -Body Problem N masses interact via mutual gravitational
More informationCDS 101 Precourse Phase Plane Analysis and Stability
CDS 101 Precourse Phase Plane Analysis and Stability Melvin Leok Control and Dynamical Systems California Institute of Technology Pasadena, CA, 26 September, 2002. mleok@cds.caltech.edu http://www.cds.caltech.edu/
More informationLectures 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 informationThe integral test and estimates of sums
The integral test Suppose f is a continuous, positive, decreasing function on [, ) and let a n = f (n). Then the series n= a n is convergent if and only if the improper integral f (x)dx is convergent.
More informationPhysics 750: Exercise 2 Tuesday, August 29, 2017
Physics 750: Exercise 2 Tuesday, August 29, 207 Use the curl command to download from the class website everything you ll need for this exercise $ WEBPATH=http://wwwphyolemissedu/~kbeach/ $ curl $WEBPATH/courses/fall207/phys750/src/exercise2tgz
More informationDynamical Systems: Lecture 1 Naima Hammoud
Dynamical Systems: Lecture 1 Naima Hammoud Feb 21, 2017 What is dynamics? Dynamics is the study of systems that evolve in time What is dynamics? Dynamics is the study of systems that evolve in time a system
More informationTheFourierTransformAndItsApplications-Lecture28
TheFourierTransformAndItsApplications-Lecture28 Instructor (Brad Osgood):All right. Let me remind you of the exam information as I said last time. I also sent out an announcement to the class this morning
More informationPhysics 115/242 The leapfrog method and other symplectic algorithms for integrating Newton s laws of motion
Physics 115/242 The leapfrog method and other symplectic algorithms for integrating Newton s laws of motion Peter Young (Dated: April 14, 2009) I. INTRODUCTION One frequently obtains detailed dynamical
More informationHandout 2 (Correction of Handout 1 plus continued discussion/hw) Comments and Homework in Chapter 1
22M:132 Fall 07 J. Simon Handout 2 (Correction of Handout 1 plus continued discussion/hw) Comments and Homework in Chapter 1 Chapter 1 contains material on sets, functions, relations, and cardinality that
More informationA short introduction with a view toward examples. Short presentation for the students of: Dynamical Systems and Complexity Summer School Volos, 2017
A short introduction with a view toward examples Center of Research and Applications of Nonlinear (CRANS) Department of Mathematics University of Patras Greece sanastassiou@gmail.com Short presentation
More informationDynamical Systems with Applications
Stephen Lynch Dynamical Systems with Applications using MATLAB Birkhauser Boston Basel Berlin Preface xi 0 A Tutorial Introduction to MATLAB and the Symbolic Math Toolbox 1 0.1 Tutorial One: The Basics
More informationAstronomy 8824: Numerical Methods Notes 2 Ordinary Differential Equations
Astronomy 8824: Numerical Methods Notes 2 Ordinary Differential Equations Reading: Numerical Recipes, chapter on Integration of Ordinary Differential Equations (which is ch. 15, 16, or 17 depending on
More informationMaths III - Numerical Methods
Maths III - Numerical Methods Matt Probert matt.probert@york.ac.uk 4 Solving Differential Equations 4.1 Introduction Many physical problems can be expressed as differential s, but solving them is not always
More informationChaotic motion. Phys 420/580 Lecture 10
Chaotic motion Phys 420/580 Lecture 10 Finite-difference equations Finite difference equation approximates a differential equation as an iterative map (x n+1,v n+1 )=M[(x n,v n )] Evolution from time t
More informationSolutions for B8b (Nonlinear Systems) Fake Past Exam (TT 10)
Solutions for B8b (Nonlinear Systems) Fake Past Exam (TT 10) Mason A. Porter 15/05/2010 1 Question 1 i. (6 points) Define a saddle-node bifurcation and show that the first order system dx dt = r x e x
More informationHandout 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
More informationIn Arnold s Mathematical Methods of Classical Mechanics (1), it
Near strongly resonant periodic orbits in a Hamiltonian system Vassili Gelfreich* Mathematics Institute, University of Warwick, Coventry CV4 7AL, United Kingdom Communicated by John N. Mather, Princeton
More informationAM 205 Final Project The N-Body Problem
AM 205 Final Project The N-Body Problem Leah Birch Elizabeth Finn Karen Yu December 14, 2012 Abstract The N-Body Problem can be solved using a variety of numeric integrators. Newton s Law of Universal
More informationx x2 2 + x3 3 x4 3. Use the divided-difference method to find a polynomial of least degree that fits the values shown: (b)
Numerical Methods - PROBLEMS. The Taylor series, about the origin, for log( + x) is x x2 2 + x3 3 x4 4 + Find an upper bound on the magnitude of the truncation error on the interval x.5 when log( + x)
More informationNewton s Cooling Model in Matlab and the Cooling Project!
Newton s Cooling Model in Matlab and the Cooling Project! James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University March 10, 2014 Outline Your Newton
More informationChaotic motion. Phys 750 Lecture 9
Chaotic motion Phys 750 Lecture 9 Finite-difference equations Finite difference equation approximates a differential equation as an iterative map (x n+1,v n+1 )=M[(x n,v n )] Evolution from time t =0to
More informationPH 120 Project # 2: Pendulum and chaos
PH 120 Project # 2: Pendulum and chaos Due: Friday, January 16, 2004 In PH109, you studied a simple pendulum, which is an effectively massless rod of length l that is fixed at one end with a small mass
More informationA 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
More informationTwo hours UNIVERSITY OF MANCHESTER SCHOOL OF COMPUTER SCIENCE. Date: Monday 22nd May 2017 Time: 09:45-11:45
COMP 14112 Two hours UNIVERSITY OF MANCHESTER SCHOOL OF COMPUTER SCIENCE Fundamentals of Artificial Intelligence Date: Monday 22nd May 2017 Time: 09:45-11:45 Answer One Question from Section A and One
More informationLocalization of Compact Invariant Sets of Nonlinear Systems
Localization of Compact Invariant of Nonlinear Systems ALEXANDER P. KRISHCHENKO Bauman Moscow State Technical University Department of Mathematical Modeling 2-aja Baumanskaja ul., 5, 105005 Moscow RUSSIA
More informationDynamical Systems with Applications using Mathematica
Stephen Lynch Dynamical Systems with Applications using Mathematica Birkhäuser Boston Basel Berlin Contents Preface xi 0 A Tutorial Introduction to Mathematica 1 0.1 A Quick Tour of Mathematica 2 0.2 Tutorial
More informationSymplectic maps. James D. Meiss. March 4, 2008
Symplectic maps James D. Meiss March 4, 2008 First used mathematically by Hermann Weyl, the term symplectic arises from a Greek word that means twining or plaiting together. This is apt, as symplectic
More informationExponentially small splitting of separatrices of the pendulum: two different examples. Marcel Guardia, Carme Olivé, Tere M-Seara
Exponentially small splitting of separatrices of the pendulum: two different examples Marcel Guardia, Carme Olivé, Tere M-Seara 1 A fast periodic perturbation of the pendulum We consider a non-autonomous
More informationPhysics 115/242 The leapfrog method and other symplectic algorithms for integrating Newton s laws of motion
Physics 115/242 The leapfrog method and other symplectic algorithms for integrating Newton s laws of motion Peter Young (Dated: April 17, 23) I. INTRODUCTION One frequently obtains detailed dynamical information
More information2 Discrete Dynamical Systems (DDS)
2 Discrete Dynamical Systems (DDS) 2.1 Basics A Discrete Dynamical System (DDS) models the change (or dynamics) of single or multiple populations or quantities in which the change occurs deterministically
More informationNonlinear Dynamical Systems Lecture - 01
Nonlinear Dynamical Systems Lecture - 01 Alexandre Nolasco de Carvalho August 08, 2017 Presentation Course contents Aims and purpose of the course Bibliography Motivation To explain what is a dynamical
More informationThank you for your purchase!
TM Thank you for your purchase! Please be sure to save a copy of this document to your local computer. This activity is copyrighted by the AIMS Education Foundation. All rights reserved. No part of this
More informationIdentifying Safe Zones for Planetary Satellite Orbiters
AIAA/AAS Astrodynamics Specialist Conference and Exhibit 16-19 August 2004, Providence, Rhode Island AIAA 2004-4862 Identifying Safe Zones for Planetary Satellite Orbiters M.E. Paskowitz and D.J. Scheeres
More informationWhen a function is defined by a fraction, the denominator of that fraction cannot be equal to zero
As stated in the previous lesson, when changing from a function to its inverse the inputs and outputs of the original function are switched, because we take the original function and solve for x. This
More information... it may happen that small differences in the initial conditions produce very great ones in the final phenomena. Henri Poincaré
Chapter 2 Dynamical Systems... it may happen that small differences in the initial conditions produce very great ones in the final phenomena. Henri Poincaré One of the exciting new fields to arise out
More informationGeometric control and dynamical systems
Université de Nice - Sophia Antipolis & Institut Universitaire de France 9th AIMS International Conference on Dynamical Systems, Differential Equations and Applications Control of an inverted pendulum
More informationOscillatory Motion. Simple pendulum: linear Hooke s Law restoring force for small angular deviations. small angle approximation. Oscillatory solution
Oscillatory Motion Simple pendulum: linear Hooke s Law restoring force for small angular deviations d 2 θ dt 2 = g l θ small angle approximation θ l Oscillatory solution θ(t) =θ 0 sin(ωt + φ) F with characteristic
More informationFinal Physics of Schwarzschild
Physics 4 Lecture 32 Final Physics of Schwarzschild Lecture 32 Physics 4 Classical Mechanics II November 6th, 27 We have studied a lot of properties of the Schwarzschild metric I want to finish with a
More informationLogarithmic, Exponential, and Other Transcendental Functions. Copyright Cengage Learning. All rights reserved.
5 Logarithmic, Exponential, and Other Transcendental Functions Copyright Cengage Learning. All rights reserved. 5.5 Bases Other Than e and Applications Copyright Cengage Learning. All rights reserved.
More informationSolving the Time-Dependent Schrödinger Equation: Instructor Notes a. Abstract
Solving the Time-Dependent Schrödinger Equation: Instructor Notes a David G. Robertson Department of Physics, Otterbein University, Westerville, OH 43081 (Dated: October 10, 2011) Abstract Some notes on
More informationExamples of linear systems and explanation of the term linear. is also a solution to this equation.
. Linear systems Examples of linear systems and explanation of the term linear. () ax b () a x + a x +... + a x b n n Illustration by another example: The equation x x + 5x 7 has one solution as x 4, x
More informationOscillatory Motion. Simple pendulum: linear Hooke s Law restoring force for small angular deviations. Oscillatory solution
Oscillatory Motion Simple pendulum: linear Hooke s Law restoring force for small angular deviations d 2 θ dt 2 = g l θ θ l Oscillatory solution θ(t) =θ 0 sin(ωt + φ) F with characteristic angular frequency
More informationCHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD
CHALMERS, GÖTEBORGS UNIVERSITET EXAM for DYNAMICAL SYSTEMS COURSE CODES: TIF 155, FIM770GU, PhD Time: Place: Teachers: Allowed material: Not allowed: August 22, 2018, at 08 30 12 30 Johanneberg Jan Meibohm,
More informationMultiple Representations: Equations to Tables and Graphs Transcript
Algebra l Teacher: It s good to see you again. Last time we talked about multiple representations. If we could, I would like to continue and discuss the subtle differences of multiple representations between
More informationPhysics 115/242 Comparison of methods for integrating the simple harmonic oscillator.
Physics 115/4 Comparison of methods for integrating the simple harmonic oscillator. Peter Young I. THE SIMPLE HARMONIC OSCILLATOR The energy (sometimes called the Hamiltonian ) of the simple harmonic oscillator
More informationPhysics 110. Electricity and Magnetism. Professor Dine. Spring, Handout: Vectors and Tensors: Everything You Need to Know
Physics 110. Electricity and Magnetism. Professor Dine Spring, 2008. Handout: Vectors and Tensors: Everything You Need to Know What makes E&M hard, more than anything else, is the problem that the electric
More informationGeometric methods for orbit integration. PiTP 2009 Scott Tremaine
Geometric methods for orbit integration PiTP 2009 Scott Tremaine spacecraft trajectories Cassini-Huygens trajectory around Saturn, 2004-2008 Planetary orbits lines = current orbits of the four inner planets
More informationUNIFORM 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 informationA GUI FOR EVOLVE ZAMS
A GUI FOR EVOLVE ZAMS D. R. Schlegel Computer Science Department Here the early work on a new user interface for the Evolve ZAMS stellar evolution code is presented. The initial goal of this project is
More informationQuantitative global phase space analysis of APM. Stability and Instability in Mechanical Systems: Applications and Numerical Tools
Quantitative global phase space analysis of APM Workshop on Stability and Instability in Mechanical Systems: Applications and Numerical Tools Carles Simó & Arturo Vieiro carles@maia.ub.es vieiro@maia.ub.es
More informationLeap Frog Solar System
Leap Frog Solar System David-Alexander Robinson Sch. 08332461 20th November 2011 Contents 1 Introduction & Theory 2 1.1 The Leap Frog Integrator......................... 2 1.2 Class.....................................
More informationExponential Integrators
Exponential Integrators John C. Bowman (University of Alberta) May 22, 2007 www.math.ualberta.ca/ bowman/talks 1 Exponential Integrators Outline Exponential Euler History Generalizations Stationary Green
More information2.2 Separable Equations
2.2 Separable Equations Definition A first-order differential equation that can be written in the form Is said to be separable. Note: the variables of a separable equation can be written as Examples Solve
More informationApplications of Advanced Mathematics (C4) Paper B: Comprehension TUESDAY 22 JANUARY 2008 Time:Upto1hour
ADVANCED GCE 4754/01B MATHEMATICS (MEI) Applications of Advanced Mathematics (C4) Paper B: Comprehension TUESDAY JANUARY 008 Afternoon Time:Upto1hour Additional materials: Rough paper MEI Examination Formulae
More informationf(x) = lim x 0 + x = lim f(x) =
Infinite Limits Having discussed in detail its as x ±, we would like to discuss in more detail its where f(x) ±. Once again we would like to emphasize that ± are not numbers, so if we write f(x) = we are
More informationThe Euler Method for the Initial Value Problem
The Euler Method for the Initial Value Problem http://people.sc.fsu.edu/ jburkardt/isc/week10 lecture 18.pdf... ISC3313: Introduction to Scientific Computing with C++ Summer Semester 2011... John Burkardt
More information2 Problem Set 2 Graphical Analysis
2 PROBLEM SET 2 GRAPHICAL ANALYSIS 2 Problem Set 2 Graphical Analysis 1. Use graphical analysis to describe all orbits of the functions below. Also draw their phase portraits. (a) F(x) = 2x There is only
More informationSOLVING ODE s NUMERICALLY WHILE PRESERVING ALL FIRST INTEGRALS
SOLVING ODE s NUMERICALLY WHILE PRESERVING ALL FIRST INTEGRALS G.R.W. QUISPEL 1,2 and H.W. CAPEL 3 Abstract. Using Discrete Gradient Methods (Quispel & Turner, J. Phys. A29 (1996) L341-L349) we construct
More informationJane and Joe are measuring the circumference of a dime with a string. Jane' s result is: 55 mm Joe's result is: 58 mm
A LESSON ON ABSOLUTE VALUE Jane and Joe are measuring the circumference of a dime with a string. Jane' s result is: 55 mm Joe's result is: 58 mm Tom knows the true length of the circumference: 56 mm. He
More informationAngle estimation using gyros and accelerometers
Angle estimation using gyros and accelerometers This version: January 23, 2018 Name: LERTEKNIK REG P-number: Date: AU T O MA RO TI C C O N T L Passed: LINKÖPING Chapter 1 Introduction The purpose of this
More informationLarge Fluctuations in Chaotic Systems
Large Fluctuations in in Chaotic Systems Igor Khovanov Physics Department, Lancaster University V.S. Anishchenko, Saratov State University N.A. Khovanova, D.G. Luchinsky, P.V.E. McClintock, Lancaster University
More informationFinding numerically Newhouse sinks near a homoclinic tangency and investigation of their chaotic transients. Takayuki Yamaguchi
Hokkaido Mathematical Journal Vol. 44 (2015) p. 277 312 Finding numerically Newhouse sinks near a homoclinic tangency and investigation of their chaotic transients Takayuki Yamaguchi (Received March 13,
More informationIntroductory Analysis I Fall 2014 Homework #5 Solutions
Introductory Analysis I Fall 2014 Homework #5 Solutions 6. Let M be a metric space, let C D M. Now we can think of C as a subset of the metric space M or as a subspace of the metric space D (D being a
More information4.4 Computing π, ln 2 and e
252 4.4 Computing π, ln 2 and e The approximations π 3.1415927, ln 2 0.69314718, e 2.7182818 can be obtained by numerical methods applied to the following initial value problems: (1) y = 4, 1 + x2 y(0)
More informationAlignment indices: a new, simple method for determining the ordered or chaotic nature of orbits
INSTITUTE OF PHYSICSPUBLISHING JOURNAL OFPHYSICSA: MATHEMATICAL AND GENERAL J. Phys. A: Math. Gen. 34 (2001) 10029 10043 PII: S0305-4470(01)25097-5 Alignment indices: a new, simple method for determining
More informationIntroduction - Motivation. Many phenomena (physical, chemical, biological, etc.) are model by differential equations. f f(x + h) f(x) (x) = lim
Introduction - Motivation Many phenomena (physical, chemical, biological, etc.) are model by differential equations. Recall the definition of the derivative of f(x) f f(x + h) f(x) (x) = lim. h 0 h Its
More informationMain topics for the First Midterm Exam
Main topics for the First Midterm Exam The final will cover Sections.-.0, 2.-2.5, and 4.. This is roughly the material from first three homeworks and three quizzes, in addition to the lecture on Monday,
More informationLecture 11 : Overview
Lecture 11 : Overview Error in Assignment 3 : In Eq. 1, Hamiltonian should be H = p2 r 2m + p2 ϕ 2mr + (p z ea z ) 2 2 2m + eφ (1) Error in lecture 10, slide 7, Eq. (21). Should be S(q, α, t) m Q = β =
More informationxt+1 = 1 ax 2 t + y t y t+1 = bx t (1)
Exercise 2.2: Hénon map In Numerical study of quadratic area-preserving mappings (Commun. Math. Phys. 50, 69-77, 1976), the French astronomer Michel Hénon proposed the following map as a model of the Poincaré
More informationUnbounded energy growth in Hamiltonian systems with a slowly varying parameter
Unbounded energy growth in Hamiltonian systems with a slowly varying parameter Vassili Gelfreich 1 Dmitry Turaev 2 1 Mathematics Institute, University of Warwick, Coventry CV4 7AL, United Kingdom Fax:
More informationQUADRATICS 3.2 Breaking Symmetry: Factoring
QUADRATICS 3. Breaking Symmetry: Factoring James Tanton SETTING THE SCENE Recall that we started our story of symmetry with a rectangle of area 36. Then you would say that the picture is wrong and that
More informationExponential Functions Dr. Laura J. Pyzdrowski
1 Names: (4 communication points) About this Laboratory An exponential function is an example of a function that is not an algebraic combination of polynomials. Such functions are called trancendental
More information