Example Chaotic Maps (that you can analyze)

Save this PDF as:
 WORD  PNG  TXT  JPG

Size: px
Start display at page:

Download "Example Chaotic Maps (that you can analyze)"

Transcription

1 Example Chaotic Maps (that you can analyze) Reading for this lecture: NDAC, Sections Lecture 7: Natural Computation & Self-Organization, Physics 256A (Winter 24); Jim Crutchfield Monday, January 27, 4

2 Example D Maps... Shift Map: x n [, ] x n+ = f(x n ) = 2x n (mod ) x n+ Fixed Point: x = Unstable: f (x )=2> Period-2 Orbit: {x } = {/3, 2/3} Unstable: (f 2 ) (x )=4> All periodic orbits unstable x n x n+2 x n Lecture 7: Natural Computation & Self-Organization, Physics 256A (Winter 24); Jim Crutchfield Monday, January 27, 4 2

3 Example D Maps... Shift Map: x n [, ] x n+ =2x n (mod ) x n+ x n Lecture 7: Natural Computation & Self-Organization, Physics 256A (Winter 24); Jim Crutchfield Monday, January 27, 4 3

4 Example D Maps... Shift Map: x n [, ] x n+ =2x n (mod ) x n+ Solvable! x n Lecture 7: Natural Computation & Self-Organization, Physics 256A (Winter 24); Jim Crutchfield Monday, January 27, 4 3

5 Example D Maps... Shift Map: x n [, ] x n+ =2x n (mod ) x n+ Solvable! x n =2 n x (mod ) x n Lecture 7: Natural Computation & Self-Organization, Physics 256A (Winter 24); Jim Crutchfield Monday, January 27, 4 3

6 Example D Maps... Shift Map: x n [, ] x n+ =2x n (mod ) x n+ Solvable! x n =2 n x (mod ) Chaotic mechanism: shift up least significant digits x n Lecture 7: Natural Computation & Self-Organization, Physics 256A (Winter 24); Jim Crutchfield Monday, January 27, 4 3

7 Example D Maps... Shift Map: x n [, ] x n+ =2x n (mod ) x n+ Solvable! x n =2 n x (mod ) Chaotic mechanism: shift up least significant digits x n x =... x =... Lecture 7: Natural Computation & Self-Organization, Physics 256A (Winter 24); Jim Crutchfield Monday, January 27, 4 3

8 Example D Maps... Shift Map: x n [, ] x n+ =2x n (mod ) x n+ Solvable! x n =2 n x (mod ) Chaotic mechanism: shift up least significant digits x n x =... x =... Lecture 7: Natural Computation & Self-Organization, Physics 256A (Winter 24); Jim Crutchfield Monday, January 27, 4 3

9 Example D Maps... Shift Map: x n [, ] x n+ =2x n (mod ) x n+ Solvable! x n =2 n x (mod ) Chaotic mechanism: shift up least significant digits x n x =... x =... Lecture 7: Natural Computation & Self-Organization, Physics 256A (Winter 24); Jim Crutchfield Monday, January 27, 4 3

10 Example D Maps... Lyapunov Characteristic Exponent for D Maps: x n+ = f(x n ) f(x) f (x ) 2 f (x ) 2 f (x )f (x )... δ N δ e λ N f(x + δ ) f(x) δ δ x x + δ or, Definition: LCE λ = lim N N log 2 δ N δ Lecture 7: Natural Computation & Self-Organization, Physics 256A (Winter 24); Jim Crutchfield Monday, January 27, 4 4

11 Example D Maps... Lyapunov Characteristic Exponent for D Maps... λ = lim N N log 2 f N (x + δ ) f N (x ) (x + δ ) x δ λ = lim N N log 2 (f N ) (x ) (f N ) (x )=f (x N )(f N ) (x )=f (x )f (x ) f (x N ) λ = lim N N N n= log 2 f (x n ) λ< stable λ> unstable Lecture 7: Natural Computation & Self-Organization, Physics 256A (Winter 24); Jim Crutchfield Monday, January 27, 4 5

12 Example D Maps... Back to Shift Map: Its LCE... x n+ = f(x n ) = 2x n (mod ) λ = lim N N N n= log 2 f (x n ) 2 x Independent of state: f (x) =2 x Amplification per step (or bits of resolution lost): λ = Lecture 7: Natural Computation & Self-Organization, Physics 256A (Winter 24); Jim Crutchfield Monday, January 27, 4 6

13 Example D Maps... Tent Map: x n [, ] { ax n, x n x n+ = 2 a( x n ), 2 <x n Slope: a [, 2] Height at max: a 2 x n+ x n Lecture 7: Natural Computation & Self-Organization, Physics 256A (Winter 24); Jim Crutchfield Monday, January 27, 4 7

14 Example D Maps... Tent Map Bifurcation Diagram: {x n }.5. a.5 2. Lecture 7: Natural Computation & Self-Organization, Physics 256A (Winter 24); Jim Crutchfield Monday, January 27, 4 8

15 Example D Maps... Tent Map... Stable fixed point: x =, a< a Unstable fixed points: {, +a }, a 2 All periodic orbits unstable: p> No periodic windows x n+2 x n Lecture 7: Natural Computation & Self-Organization, Physics 256A (Winter 24); Jim Crutchfield Monday, January 27, 4 9

16 Example D Maps... Tent Map LCE: λ = lim N N N n= log 2 f (x n ) Lecture 7: Natural Computation & Self-Organization, Physics 256A (Winter 24); Jim Crutchfield Monday, January 27, 4

17 Example D Maps... Tent Map LCE: λ = lim N N N n= log 2 f (x n ) λ = lim N N N n= log 2 ±a = log 2 a Lecture 7: Natural Computation & Self-Organization, Physics 256A (Winter 24); Jim Crutchfield Monday, January 27, 4

18 Example D Maps... Tent Map LCE: λ = lim N N N n= log 2 ±a = log 2 a.5 Chaotic -.5 Periodic: P Lecture 7: Natural Computation & Self-Organization, Physics 256A (Winter 24); Jim Crutchfield Monday, January 27, 4

19 Example D Maps... Logistic map: State space: Parameter (height): x n+ = rx n ( x n ) x n [, ] r [, 4] x n+ x n Lecture 7: Natural Computation & Self-Organization, Physics 256A (Winter 24); Jim Crutchfield Monday, January 27, 4 2

20 Example D Maps... Logistic map bifurcation diagram.... {x n }. 3. r 4. Lecture 7: Natural Computation & Self-Organization, Physics 256A (Winter 24); Jim Crutchfield Monday, January 27, 4 3

21 Example D Maps... Logistic map LCE: Local stability depends on state: λ = lim N λ = lim N N N N n= N n= log 2 f (x n ) log 2 r( 2x n ) f (x) =r( 2x) Period : x =, r f (x ) = r λ = log 2 r Period : x = r, r 3 λ = log r 2 2 r Superstable: f (x i )= λ r =2 Bifurcations: λ = Onset of chaos: λ = Lecture 7: Natural Computation & Self-Organization, Physics 256A (Winter 24); Jim Crutchfield Monday, January 27, 4 4

22 Example D Maps... LCE for D Maps... an aside on the Ergodic Theorem Rather than time average: λ = lim N N N n= log 2 f (x n ) Average over attractor s distribution: Pr(x), x Λ Invariant distribution: Pr(x) =f Pr(x) State-space averaged LCE: λ = Λ dx Pr(x) log 2 f (x) Lecture 7: Natural Computation & Self-Organization, Physics 256A (Winter 24); Jim Crutchfield Monday, January 27, 4 5

23 Example D Maps... Logistic map LCE: r =4 Invariant distribution: Pr(x) = π x( x) f(x) log 2 Pr(x) λ = dx log 2 4 8x π x( x) λ = bit per step x Lecture 7: Natural Computation & Self-Organization, Physics 256A (Winter 24); Jim Crutchfield Monday, January 27, 4 6

24 Example D Maps... LCE view of period-doubling route to chaos: r B r 2B B.5 r c λ = r 2 4 λ = Superstable Period-2 λ Superstable Period-4 r =3.83 Period-3 Window Lecture 7: Natural Computation & Self-Organization, Physics 256A (Winter 24); Jim Crutchfield Monday, January 27, 4 7

25 Example D Maps... Logistic map bifurcation diagram self-similarity. {x n }. 3. r {x n } r {x n } Lecture 7: Natural Computation & Self-Organization, Physics 256A (Winter 24); Jim Crutchfield r Monday, January 27, 4 8

26 Example D Maps... Logistic map bifurcation diagram self-similarity. {x n }. 3. r {x n } r {x n } Lecture 7: Natural Computation & Self-Organization, Physics 256A (Winter 24); Jim Crutchfield r Monday, January 27, 4 8

27 Example D Maps... Bifurcation Theory of D Maps... Scaling analysis of period-doubling cascade:. {x n } d d 3 x =/2 d r 2ss r r 2 r 4ss r 8ss Universal constants: Lecture 7: Natural Computation & Self-Organization, Physics 256A (Winter 24); Jim Crutchfield δ = lim n α = lim n r n = r n+ d n = d n+ Monday, January 27, 4 9

28 Example D Maps... Bifurcation Theory of D Maps... Renormalization group analysis of period-doubling: r r f 2 r (x) f r (x) αy y αx x > <, flip Lecture 7: Natural Computation & Self-Organization, Physics 256A (Winter 24); Jim Crutchfield Monday, January 27, 4 2

29 Example D Maps... Bifurcation Theory of D Maps... Renormalization group analysis of period-doubling... f(x, r ) αf 2 ( x α,r ) f 2 ( x α,r ) α 2 f 4 ( x α 2,r 2)... f(x, r ) α n f (2n) ( x α n,r n) Lecture 7: Natural Computation & Self-Organization, Physics 256A (Winter 24); Jim Crutchfield Monday, January 27, 4 2

30 Example D Maps... Bifurcation Theory of D Maps... Renormalization group analysis of period-doubling... Universal Map: g (x) = lim n αn f (2n) ( x α n,r n) for x x max Lecture 7: Natural Computation & Self-Organization, Physics 256A (Winter 24); Jim Crutchfield Monday, January 27, 4 22

31 Example D Maps... Bifurcation Theory of D Maps... Renormalization group analysis of period-doubling... r : f(x, r ) αf 2 ( x α,r ) for x x max Limiting functional equation: ( x max = ) ( ) g(x) =αg 2 x α Boundary conditions: g() = & g () = Lecture 7: Natural Computation & Self-Organization, Physics 256A (Winter 24); Jim Crutchfield Monday, January 27, 4 23

32 Example D Maps... Bifurcation Theory of D Maps... Renormalization group analysis of period-doubling... How to solve? g(x) =αg 2 ( x α ) with g() = & g () = Taylor expansion: g(x) =a + bx 2 + cx 4 + g(x) even Find: α = Parameter rescaling: (more work) Find: δ = Lecture 7: Natural Computation & Self-Organization, Physics 256A (Winter 24); Jim Crutchfield Monday, January 27, 4 24

33 Example D Maps... Reading for next lecture: Lecture Notes. Lecture 7: Natural Computation & Self-Organization, Physics 256A (Winter 24); Jim Crutchfield Monday, January 27, 4 25

Solutions to homework assignment #7 Math 119B UC Davis, Spring for 1 r 4. Furthermore, the derivative of the logistic map is. L r(x) = r(1 2x).

Solutions to homework assignment #7 Math 119B UC Davis, Spring for 1 r 4. Furthermore, the derivative of the logistic map is. L r(x) = r(1 2x). Solutions to homework assignment #7 Math 9B UC Davis, Spring 0. A fixed point x of an interval map T is called superstable if T (x ) = 0. Find the value of 0 < r 4 for which the logistic map L r has a

More information

One dimensional Maps

One 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 information

Renormalization Group Theory

Renormalization Group Theory Chapter 16 Renormalization Group Theory In the previous chapter a procedure was developed where hiher order 2 n cycles were related to lower order cycles throuh a functional composition and rescalin procedure.

More information

Chapter 4. Transition towards chaos. 4.1 One-dimensional maps

Chapter 4. Transition towards chaos. 4.1 One-dimensional maps Chapter 4 Transition towards chaos In this chapter we will study how successive bifurcations can lead to chaos when a parameter is tuned. It is not an extensive review : there exists a lot of different

More information

PHY411 Lecture notes Part 4

PHY411 Lecture notes Part 4 PHY411 Lecture notes Part 4 Alice Quillen February 1, 2016 Contents 0.1 Introduction.................................... 2 1 Bifurcations of one-dimensional dynamical systems 2 1.1 Saddle-node bifurcation.............................

More information

Introduction to Dynamical Systems Basic Concepts of Dynamics

Introduction to Dynamical Systems Basic Concepts of Dynamics Introduction to Dynamical Systems Basic Concepts of Dynamics A dynamical system: Has a notion of state, which contains all the information upon which the dynamical system acts. A simple set of deterministic

More information

Bifurcation of Fixed Points

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

Anatomy of a Bit. Reading for this lecture: CMR articles Yeung and Anatomy.

Anatomy of a Bit. Reading for this lecture: CMR articles Yeung and Anatomy. Anatomy of a Bit Reading for this lecture: CMR articles Yeung and Anatomy. 1 Information in a single measurement: Alternative rationale for entropy and entropy rate. k measurement outcomes Stored in log

More information

Lecture 7. Order Out of Chaos

Lecture 7. Order Out of Chaos Lecture 7 Order Out of Chaos Lyapunov Exponents: Recall from Last Time The Lyapunov exponent for maps is: < 0if f 0 (x ) < 1 =ln f 0 (x ) = > 0if f 0 (x ) > 1 Lyapunov Exponents: Numerical Estimation Lyapunov

More information

Logistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations

Logistic 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 information

CHAPTER 3 CHAOTIC MAPS BASED PSEUDO RANDOM NUMBER GENERATORS

CHAPTER 3 CHAOTIC MAPS BASED PSEUDO RANDOM NUMBER GENERATORS 24 CHAPTER 3 CHAOTIC MAPS BASED PSEUDO RANDOM NUMBER GENERATORS 3.1 INTRODUCTION Pseudo Random Number Generators (PRNGs) are widely used in many applications, such as numerical analysis, probabilistic

More information

Pascal (Yang Hui) triangles and power laws in the logistic map

Pascal (Yang Hui) triangles and power laws in the logistic map Journal of Physics: Conference Series PAPER OPEN ACCESS Pascal (Yang Hui) triangles and power laws in the logistic map To cite this article: Carlos Velarde and Alberto Robledo 2015 J. Phys.: Conf. Ser.

More information

Part II. Dynamical Systems. Year

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

One Dimensional Dynamical Systems

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

More information

DYNAMICAL SYSTEMS PROBLEMS. asgor/ (1) Which of the following maps are topologically transitive (minimal,

DYNAMICAL SYSTEMS PROBLEMS.  asgor/ (1) Which of the following maps are topologically transitive (minimal, DYNAMICAL SYSTEMS PROBLEMS http://www.math.uci.edu/ asgor/ (1) Which of the following maps are topologically transitive (minimal, topologically mixing)? identity map on a circle; irrational rotation of

More information

6.2 Brief review of fundamental concepts about chaotic systems

6.2 Brief review of fundamental concepts about chaotic systems 6.2 Brief review of fundamental concepts about chaotic systems Lorenz (1963) introduced a 3-variable model that is a prototypical example of chaos theory. These equations were derived as a simplification

More information

Nonlinear Dynamics and Chaos

Nonlinear Dynamics and Chaos Ian Eisenman eisenman@fas.harvard.edu Geological Museum 101, 6-6352 Nonlinear Dynamics and Chaos Review of some of the topics covered in homework problems, based on section notes. December, 2005 Contents

More information

Integration - Past Edexcel Exam Questions

Integration - Past Edexcel Exam Questions Integration - Past Edexcel Exam Questions 1. (a) Given that y = 5x 2 + 7x + 3, find i. - ii. - (b) ( 1 + 3 ) x 1 x dx. [4] 2. Question 2b - January 2005 2. The gradient of the curve C is given by The point

More information

arxiv:cond-mat/ v1 8 Jan 2004

arxiv:cond-mat/ v1 8 Jan 2004 Multifractality and nonextensivity at the edge of chaos of unimodal maps E. Mayoral and A. Robledo arxiv:cond-mat/0401128 v1 8 Jan 2004 Instituto de Física, Universidad Nacional Autónoma de México, Apartado

More information

The Existence of Chaos in the Lorenz System

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

2 Discrete Dynamical Systems: Maps

2 Discrete Dynamical Systems: Maps 19 2 Discrete Dynamical Systems: Maps 2.1 Introduction Many physical systems displaying chaotic behavior are accurately described by mathematical models derived from well-understood physical principles.

More information

QUANTIFYING THE FOLDING MECHANISM IN CHAOTIC DYNAMICS

QUANTIFYING THE FOLDING MECHANISM IN CHAOTIC DYNAMICS QUANTIFYING THE FOLDING MECHANISM IN CHAOTIC DYNAMICS V. BARAN, M. ZUS,2, A. BONASERA 3, A. PATURCA Faculty of Physics, University of Bucharest, 405 Atomiştilor, POB MG-, RO-07725, Bucharest-Măgurele,

More information

LECTURE 16 GAUSS QUADRATURE In general for Newton-Cotes (equispaced interpolation points/ data points/ integration points/ nodes).

LECTURE 16 GAUSS QUADRATURE In general for Newton-Cotes (equispaced interpolation points/ data points/ integration points/ nodes). CE 025 - Lecture 6 LECTURE 6 GAUSS QUADRATURE In general for ewton-cotes (equispaced interpolation points/ data points/ integration points/ nodes). x E x S fx dx hw' o f o + w' f + + w' f + E 84 f 0 f

More information

arxiv:chao-dyn/ v1 5 Mar 1996

arxiv:chao-dyn/ v1 5 Mar 1996 Turbulence in Globally Coupled Maps M. G. Cosenza and A. Parravano Centro de Astrofísica Teórica, Facultad de Ciencias, Universidad de Los Andes, A. Postal 26 La Hechicera, Mérida 5251, Venezuela (To appear,

More information

Integration of Rational Functions by Partial Fractions

Integration of Rational Functions by Partial Fractions Title Integration of Rational Functions by Partial Fractions MATH 1700 December 6, 2016 MATH 1700 Partial Fractions December 6, 2016 1 / 11 Readings Readings Readings: Section 7.4 MATH 1700 Partial Fractions

More information

OPER 627: Nonlinear Optimization Lecture 2: Math Background and Optimality Conditions

OPER 627: Nonlinear Optimization Lecture 2: Math Background and Optimality Conditions OPER 627: Nonlinear Optimization Lecture 2: Math Background and Optimality Conditions Department of Statistical Sciences and Operations Research Virginia Commonwealth University Aug 28, 2013 (Lecture 2)

More information

The logistic difference equation and the route to chaotic behaviour

The logistic difference equation and the route to chaotic behaviour The logistic difference equation and the route to chaotic behaviour Level 1 module in Modelling course in population and evolutionary biology (701-1418-00) Module author: Sebastian Bonhoeffer Course director:

More information

Chaotic motion. Phys 750 Lecture 9

Chaotic 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 information

RELAXATION AND TRANSIENTS IN A TIME-DEPENDENT LOGISTIC MAP

RELAXATION AND TRANSIENTS IN A TIME-DEPENDENT LOGISTIC MAP International Journal of Bifurcation and Chaos, Vol. 12, No. 7 (2002) 1667 1674 c World Scientific Publishing Company RELAATION AND TRANSIENTS IN A TIME-DEPENDENT LOGISTIC MAP EDSON D. LEONEL, J. KAMPHORST

More information

Dynamical Systems. Bernard Deconinck Department of Applied Mathematics University of Washington Campus Box Seattle, WA, 98195, USA

Dynamical Systems. Bernard Deconinck Department of Applied Mathematics University of Washington Campus Box Seattle, WA, 98195, USA Dynamical Systems Bernard Deconinck Department of Applied Mathematics University of Washington Campus Box 352420 Seattle, WA, 98195, USA June 4, 2009 i Prolegomenon These are the lecture notes for Amath

More information

MA 1128: Lecture 19 4/20/2018. Quadratic Formula Solving Equations with Graphs

MA 1128: Lecture 19 4/20/2018. Quadratic Formula Solving Equations with Graphs MA 1128: Lecture 19 4/20/2018 Quadratic Formula Solving Equations with Graphs 1 Completing-the-Square Formula One thing you may have noticed when you were completing the square was that you followed the

More information

Contents Dynamical Systems Stability of Dynamical Systems: Linear Approach

Contents Dynamical Systems Stability of Dynamical Systems: Linear Approach Contents 1 Dynamical Systems... 1 1.1 Introduction... 1 1.2 DynamicalSystems andmathematical Models... 1 1.3 Kinematic Interpretation of a System of Differential Equations... 3 1.4 Definition of a Dynamical

More information

MS 2001: Test 1 B Solutions

MS 2001: Test 1 B Solutions MS 2001: Test 1 B Solutions Name: Student Number: Answer all questions. Marks may be lost if necessary work is not clearly shown. Remarks by me in italics and would not be required in a test - J.P. Question

More information

Math123 Lecture 1. Dr. Robert C. Busby. Lecturer: Office: Korman 266 Phone :

Math123 Lecture 1. Dr. Robert C. Busby. Lecturer: Office: Korman 266 Phone : Lecturer: Math1 Lecture 1 Dr. Robert C. Busby Office: Korman 66 Phone : 15-895-1957 Email: rbusby@mcs.drexel.edu Course Web Site: http://www.mcs.drexel.edu/classes/calculus/math1_spring0/ (Links are case

More information

CALCULUS JIA-MING (FRANK) LIOU

CALCULUS JIA-MING (FRANK) LIOU CALCULUS JIA-MING (FRANK) LIOU Abstract. Contents. Power Series.. Polynomials and Formal Power Series.2. Radius of Convergence 2.3. Derivative and Antiderivative of Power Series 4.4. Power Series Expansion

More information

8 + 6) x 2 ) y = h(x)

8 + 6) x 2 ) y = h(x) . a. Horizontal shift 6 left and vertical shift up. Notice B' is ( 6, ) and B is (0, 0). b. h(x) = 0.5(x + 6) + (Enter in a grapher to check.) c. Use the graph. Notice A' to see h(x) crosses the x-axis

More information

The Research of Railway Coal Dispatched Volume Prediction Based on Chaos Theory

The Research of Railway Coal Dispatched Volume Prediction Based on Chaos Theory The Research of Railway Coal Dispatched Volume Prediction Based on Chaos Theory Hua-Wen Wu Fu-Zhang Wang Institute of Computing Technology, China Academy of Railway Sciences Beijing 00044, China, P.R.

More information

Chapter 1 Bifurcations and Chaos in Dynamical Systems

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

More information

Chapter 1 Mathematical Preliminaries and Error Analysis

Chapter 1 Mathematical Preliminaries and Error Analysis Chapter 1 Mathematical Preliminaries and Error Analysis Per-Olof Persson persson@berkeley.edu Department of Mathematics University of California, Berkeley Math 128A Numerical Analysis Limits and Continuity

More information

Lecture 32: Taylor Series and McLaurin series We saw last day that some functions are equal to a power series on part of their domain.

Lecture 32: Taylor Series and McLaurin series We saw last day that some functions are equal to a power series on part of their domain. Lecture 32: Taylor Series and McLaurin series We saw last day that some functions are equal to a power series on part of their domain. For example f(x) = 1 1 x = 1 + x + x2 + x 3 + = ln(1 + x) = x x2 2

More information

Linear Algebra Review

Linear Algebra Review January 29, 2013 Table of contents Metrics Metric Given a space X, then d : X X R + 0 and z in X if: d(x, y) = 0 is equivalent to x = y d(x, y) = d(y, x) d(x, y) d(x, z) + d(z, y) is a metric is for all

More information

f (x) = k=0 f (0) = k=0 k=0 a k k(0) k 1 = a 1 a 1 = f (0). a k k(k 1)x k 2, k=2 a k k(k 1)(0) k 2 = 2a 2 a 2 = f (0) 2 a k k(k 1)(k 2)x k 3, k=3

f (x) = k=0 f (0) = k=0 k=0 a k k(0) k 1 = a 1 a 1 = f (0). a k k(k 1)x k 2, k=2 a k k(k 1)(0) k 2 = 2a 2 a 2 = f (0) 2 a k k(k 1)(k 2)x k 3, k=3 1 M 13-Lecture Contents: 1) Taylor Polynomials 2) Taylor Series Centered at x a 3) Applications of Taylor Polynomials Taylor Series The previous section served as motivation and gave some useful expansion.

More information

Tropical Polynomials

Tropical Polynomials 1 Tropical Arithmetic Tropical Polynomials Los Angeles Math Circle, May 15, 2016 Bryant Mathews, Azusa Pacific University In tropical arithmetic, we define new addition and multiplication operations on

More information

Dynamical Behavior of Logistic Equation through Symbiosis Model

Dynamical Behavior of Logistic Equation through Symbiosis Model International Journal of Scientific and Research Publications, Volume 3, Issue 6, June 213 1 ISSN 225-3153 Dynamical Behavior of Logistic Equation through Symbiosis Model M.Sambasiavm Assistant Professor

More information

UNIVERSITY OF BRISTOL. Mock exam paper for examination for the Degrees of B.Sc. and M.Sci. (Level 3)

UNIVERSITY OF BRISTOL. Mock exam paper for examination for the Degrees of B.Sc. and M.Sci. (Level 3) UNIVERSITY OF BRISTOL Mock exam paper for examination for the Degrees of B.Sc. and M.Sci. (Level 3) DYNAMICAL SYSTEMS and ERGODIC THEORY MATH 36206 (Paper Code MATH-36206) 2 hours and 30 minutes This paper

More information

Chaotic motion. Phys 420/580 Lecture 10

Chaotic 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 information

The Laplace Transform. Background: Improper Integrals

The Laplace Transform. Background: Improper Integrals The Laplace Transform Background: Improper Integrals Recall: Definite Integral: a, b real numbers, a b; f continuous on [a, b] b a f(x) dx 1 Improper integrals: Type I Infinite interval of integration

More information

Problem Set Number 1, j/2.036j MIT (Fall 2014)

Problem Set Number 1, j/2.036j MIT (Fall 2014) Problem Set Number 1, 18.385j/2.036j MIT (Fall 2014) Rodolfo R. Rosales (MIT, Math. Dept., Cambridge, MA 02139) September 9, 2014 Due Fri., September 19, 2014. 1 Get equation from phase line portrait problem

More information

Basic Theory of Linear Differential Equations

Basic Theory of Linear Differential Equations Basic Theory of Linear Differential Equations Picard-Lindelöf Existence-Uniqueness Vector nth Order Theorem Second Order Linear Theorem Higher Order Linear Theorem Homogeneous Structure Recipe for Constant-Coefficient

More information

NUMERICALLY COMPUTING THE LYAPUNOV EXPONENTS OF MATRIX-VALUED COCYCLES

NUMERICALLY COMPUTING THE LYAPUNOV EXPONENTS OF MATRIX-VALUED COCYCLES NUMERICALLY COMPUTING THE LYAPUNOV EXPONENTS OF MATRIX-VALUED COCYCLES RODRIGO TREVIÑO This short note is based on a talk I gave at the student dynamical systems seminar about using your computer to figure

More information

Simple Iteration, cont d

Simple Iteration, cont d Jim Lambers MAT 772 Fall Semester 2010-11 Lecture 2 Notes These notes correspond to Section 1.2 in the text. Simple Iteration, cont d In general, nonlinear equations cannot be solved in a finite sequence

More information

Chapter 6: Ensemble Forecasting and Atmospheric Predictability. Introduction

Chapter 6: Ensemble Forecasting and Atmospheric Predictability. Introduction Chapter 6: Ensemble Forecasting and Atmospheric Predictability Introduction Deterministic Chaos (what!?) In 1951 Charney indicated that forecast skill would break down, but he attributed it to model errors

More information

Math 2: Algebra 2, Geometry and Statistics Ms. Sheppard-Brick Chapter 4 Test Review

Math 2: Algebra 2, Geometry and Statistics Ms. Sheppard-Brick Chapter 4 Test Review Chapter 4 Test Review Students will be able to (SWBAT): Write an explicit and a recursive function rule for a linear table of values. Write an explicit function rule for a quadratic table of values. Determine

More information

Oscillatory 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. 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 information

Nonlinear Systems and Control Lecture # 12 Converse Lyapunov Functions & Time Varying Systems. p. 1/1

Nonlinear Systems and Control Lecture # 12 Converse Lyapunov Functions & Time Varying Systems. p. 1/1 Nonlinear Systems and Control Lecture # 12 Converse Lyapunov Functions & Time Varying Systems p. 1/1 p. 2/1 Converse Lyapunov Theorem Exponential Stability Let x = 0 be an exponentially stable equilibrium

More information

Smooth Ergodic Theory and Nonuniformly Hyperbolic Dynamics

Smooth Ergodic Theory and Nonuniformly Hyperbolic Dynamics CHAPTER 2 Smooth Ergodic Theory and Nonuniformly Hyperbolic Dynamics Luis Barreira Departamento de Matemática, Instituto Superior Técnico, 1049-001 Lisboa, Portugal E-mail: barreira@math.ist.utl.pt url:

More information

Entrainment and Chaos in the Hodgkin-Huxley Oscillator

Entrainment and Chaos in the Hodgkin-Huxley Oscillator Entrainment and Chaos in the Hodgkin-Huxley Oscillator Kevin K. Lin http://www.cims.nyu.edu/ klin Courant Institute, New York University Mostly Biomath - 2005.4.5 p.1/42 Overview (1) Goal: Show that the

More information

CHAPTER 2 POLYNOMIALS KEY POINTS

CHAPTER 2 POLYNOMIALS KEY POINTS CHAPTER POLYNOMIALS KEY POINTS 1. Polynomials of degrees 1, and 3 are called linear, quadratic and cubic polynomials respectively.. A quadratic polynomial in x with real coefficient is of the form a x

More information

Math 12: Discrete Dynamical Systems Homework

Math 12: Discrete Dynamical Systems Homework Math 12: Discrete Dynamical Systems Homework Department of Mathematics, Harvey Mudd College Several of these problems will require computational software to help build our insight about discrete dynamical

More information

Role of multistability in the transition to chaotic phase synchronization

Role of multistability in the transition to chaotic phase synchronization Downloaded from orbit.dtu.dk on: Jan 24, 2018 Role of multistability in the transition to chaotic phase synchronization Postnov, D.E.; Vadivasova, T.E.; Sosnovtseva, Olga; Balanov, A.G.; Anishchenko, V.S.;

More information

Dynamical Systems in Biology

Dynamical Systems in Biology Dynamical Systems in Biology Hal Smith A R I Z O N A S T A T E U N I V E R S I T Y H.L. Smith (ASU) Dynamical Systems in Biology ASU, July 5, 2012 1 / 31 Outline 1 What s special about dynamical systems

More information

Dynamical behaviour of a controlled vibro-impact system

Dynamical behaviour of a controlled vibro-impact system Vol 17 No 7, July 2008 c 2008 Chin. Phys. Soc. 1674-1056/2008/17(07)/2446-05 Chinese Physics B and IOP Publishing Ltd Dynamical behaviour of a controlled vibro-impact system Wang Liang( ), Xu Wei( ), and

More information

Discrete (Difference) Equations

Discrete (Difference) Equations Discrete (Difference) Equations A discrete, or difference, equation expresses a relationship between the elements of a sequence, {y n }, where n N 0 {0, 1, 2,...}. For example, the trivial difference equation

More information

ENGI 9420 Lecture Notes 4 - Stability Analysis Page Stability Analysis for Non-linear Ordinary Differential Equations

ENGI 9420 Lecture Notes 4 - Stability Analysis Page Stability Analysis for Non-linear Ordinary Differential Equations ENGI 940 Lecture Notes 4 - Stability Analysis Page 4.01 4. Stability Analysis for Non-linear Ordinary Differential Equations A pair of simultaneous first order homogeneous linear ordinary differential

More information

1. State the informal definition of the limit of a function. 4. What are the values of the following limits: lim. 2x 1 if x 1 discontinuous?

1. State the informal definition of the limit of a function. 4. What are the values of the following limits: lim. 2x 1 if x 1 discontinuous? Review for Limits Key Concepts 1 State the informal definition of the limit of a function 2 State 9 limit theorems 3 State the squeeze theorem 4 What are the values of the following limits: lim sin x x

More information

Section Taylor and Maclaurin Series

Section Taylor and Maclaurin Series Section.0 Taylor and Maclaurin Series Ruipeng Shen Feb 5 Taylor and Maclaurin Series Main Goal: How to find a power series representation for a smooth function us assume that a smooth function has a power

More information

Solutions 2017 AB Exam

Solutions 2017 AB Exam 1. Solve for x : x 2 = 4 x. Solutions 2017 AB Exam Texas A&M High School Math Contest October 21, 2017 ANSWER: x = 3 Solution: x 2 = 4 x x 2 = 16 8x + x 2 x 2 9x + 18 = 0 (x 6)(x 3) = 0 x = 6, 3 but x

More information

Report due date. Please note: report has to be handed in by Monday, May 16 noon.

Report due date. Please note: report has to be handed in by Monday, May 16 noon. Lecture 23 18.86 Report due date Please note: report has to be handed in by Monday, May 16 noon. Course evaluation: Please submit your feedback (you should have gotten an email) From the weak form to finite

More information

Study on Proportional Synchronization of Hyperchaotic Circuit System

Study on Proportional Synchronization of Hyperchaotic Circuit System Commun. Theor. Phys. (Beijing, China) 43 (25) pp. 671 676 c International Academic Publishers Vol. 43, No. 4, April 15, 25 Study on Proportional Synchronization of Hyperchaotic Circuit System JIANG De-Ping,

More information

Bifurcation control and chaos in a linear impulsive system

Bifurcation control and chaos in a linear impulsive system Vol 8 No 2, December 2009 c 2009 Chin. Phys. Soc. 674-056/2009/82)/5235-07 Chinese Physics B and IOP Publishing Ltd Bifurcation control and chaos in a linear impulsive system Jiang Gui-Rong 蒋贵荣 ) a)b),

More information

PUTNAM PROBLEMS SEQUENCES, SERIES AND RECURRENCES. Notes

PUTNAM PROBLEMS SEQUENCES, SERIES AND RECURRENCES. Notes PUTNAM PROBLEMS SEQUENCES, SERIES AND RECURRENCES Notes. x n+ = ax n has the general solution x n = x a n. 2. x n+ = x n + b has the general solution x n = x + (n )b. 3. x n+ = ax n + b (with a ) can be

More information

Generalized projective synchronization of a class of chaotic (hyperchaotic) systems with uncertain parameters

Generalized projective synchronization of a class of chaotic (hyperchaotic) systems with uncertain parameters Vol 16 No 5, May 2007 c 2007 Chin. Phys. Soc. 1009-1963/2007/16(05)/1246-06 Chinese Physics and IOP Publishing Ltd Generalized projective synchronization of a class of chaotic (hyperchaotic) systems with

More information

NONLINEAR DYNAMICS PHYS 471 & PHYS 571

NONLINEAR DYNAMICS PHYS 471 & PHYS 571 NONLINEAR DYNAMICS PHYS 471 & PHYS 571 Prof. R. Gilmore 12-918 X-2779 robert.gilmore@drexel.edu Office hours: 14:00 Quarter: Winter, 2014-2015 Course Schedule: Tuesday, Thursday, 11:00-12:20 Room: 12-919

More information

Disturbance Decoupling Problem

Disturbance Decoupling Problem DISC Systems and Control Theory of Nonlinear Systems, 21 1 Disturbance Decoupling Problem Lecture 4 Nonlinear Dynamical Control Systems, Chapter 7 The disturbance decoupling problem is a typical example

More information

Math 113/113H Winter 2006 Departmental Final Exam

Math 113/113H Winter 2006 Departmental Final Exam Name KEY Instructor Section No. Student Number Math 3/3H Winter 26 Departmental Final Exam Instructions: The time limit is 3 hours. Problems -6 short-answer questions, each worth 2 points. Problems 7 through

More information

Chaotic Vibrations. An Introduction for Applied Scientists and Engineers

Chaotic Vibrations. An Introduction for Applied Scientists and Engineers Chaotic Vibrations An Introduction for Applied Scientists and Engineers FRANCIS C. MOON Theoretical and Applied Mechanics Cornell University Ithaca, New York A WILEY-INTERSCIENCE PUBLICATION JOHN WILEY

More information

WHAT IS A CHAOTIC ATTRACTOR?

WHAT IS A CHAOTIC ATTRACTOR? WHAT IS A CHAOTIC ATTRACTOR? CLARK ROBINSON Abstract. Devaney gave a mathematical definition of the term chaos, which had earlier been introduced by Yorke. We discuss issues involved in choosing the properties

More information

2. Metric Spaces. 2.1 Definitions etc.

2. Metric Spaces. 2.1 Definitions etc. 2. Metric Spaces 2.1 Definitions etc. The procedure in Section for regarding R as a topological space may be generalized to many other sets in which there is some kind of distance (formally, sets with

More information

DIFFERENT TYPES OF CRITICAL BEHAVIOR IN CONSERVATIVELY COUPLED HÉNON MAPS

DIFFERENT TYPES OF CRITICAL BEHAVIOR IN CONSERVATIVELY COUPLED HÉNON MAPS DIFFERENT TYPES OF CRITICAL BEHAVIOR IN CONSERVATIVELY COUPLED HÉNON MAPS D.V. Savin, *, A.P. Kuznetsov,, A.V. Savin and U. Feudel 3 Department of Nonlinear Processes, ernyshevsky Saratov State University,

More information

Dynamical analysis and circuit simulation of a new three-dimensional chaotic system

Dynamical analysis and circuit simulation of a new three-dimensional chaotic system Dynamical analysis and circuit simulation of a new three-dimensional chaotic system Wang Ai-Yuan( 王爱元 ) a)b) and Ling Zhi-Hao( 凌志浩 ) a) a) Department of Automation, East China University of Science and

More information

A simple electronic circuit to demonstrate bifurcation and chaos

A simple electronic circuit to demonstrate bifurcation and chaos A simple electronic circuit to demonstrate bifurcation and chaos P R Hobson and A N Lansbury Brunel University, Middlesex Chaos has generated much interest recently, and many of the important features

More information

Families of Functions, Taylor Polynomials, l Hopital s

Families of Functions, Taylor Polynomials, l Hopital s Unit #6 : Rule Families of Functions, Taylor Polynomials, l Hopital s Goals: To use first and second derivative information to describe functions. To be able to find general properties of families of functions.

More information

CHALLENGE! (0) = 5. Construct a polynomial with the following behavior at x = 0:

CHALLENGE! (0) = 5. Construct a polynomial with the following behavior at x = 0: TAYLOR SERIES Construct a polynomial with the following behavior at x = 0: CHALLENGE! P( x) = a + ax+ ax + ax + ax 2 3 4 0 1 2 3 4 P(0) = 1 P (0) = 2 P (0) = 3 P (0) = 4 P (4) (0) = 5 Sounds hard right?

More information

Numerical Methods for PDEs

Numerical Methods for PDEs Numerical Methods for PDEs Partial Differential Equations (Lecture 1, Week 1) Markus Schmuck Department of Mathematics and Maxwell Institute for Mathematical Sciences Heriot-Watt University, Edinburgh

More information

2 Qualitative theory of non-smooth dynamical systems

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

More information

1 Question related to polynomials

1 Question related to polynomials 07-08 MATH00J Lecture 6: Taylor Series Charles Li Warning: Skip the material involving the estimation of error term Reference: APEX Calculus This lecture introduced Taylor Polynomial and Taylor Series

More information

Deterministic chaos and diffusion in maps and billiards

Deterministic chaos and diffusion in maps and billiards Deterministic chaos and diffusion in maps and billiards Rainer Klages Queen Mary University of London, School of Mathematical Sciences Mathematics for the Fluid Earth Newton Institute, Cambridge, 14 November

More information

E209A: Analysis and Control of Nonlinear Systems Problem Set 6 Solutions

E209A: Analysis and Control of Nonlinear Systems Problem Set 6 Solutions E9A: Analysis and Control of Nonlinear Systems Problem Set 6 Solutions Michael Vitus Gabe Hoffmann Stanford University Winter 7 Problem 1 The governing equations are: ẋ 1 = x 1 + x 1 x ẋ = x + x 3 Using

More information

Thermalization in Quantum Systems

Thermalization in Quantum Systems Thermalization in Quantum Systems Jonas Larson Stockholm University and Universität zu Köln Dresden 18/4-2014 Motivation Long time evolution of closed quantum systems not fully understood. Cold atom system

More information

A Classical Approach to the Stark-Effect. Mridul Mehta Advisor: Prof. Enrique J. Galvez Physics Dept., Colgate University

A Classical Approach to the Stark-Effect. Mridul Mehta Advisor: Prof. Enrique J. Galvez Physics Dept., Colgate University A Classical Approach to the Stark-Effect Mridul Mehta Advisor: Prof. Enrique J. Galvez Physics Dept., Colgate University Abstract The state of an atom in the presence of an external electric field is known

More information

Strange-Nonchaotic-Attractor-Like Behaviors in Coupled Map Systems

Strange-Nonchaotic-Attractor-Like Behaviors in Coupled Map Systems Commun. Theor. Phys. (Beijing, China) 35 (2001) pp. 389 394 c International Academic Publishers Vol. 35, No. 4, April 15, 2001 Strange-Nonchaotic-Attractor-Like Behaviors in Coupled Map Systems SHI Peng-Liang,

More information

Period-doubling cascades of a Silnikov equation

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

APPPHYS217 Tuesday 25 May 2010

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

THE ESSENCE OF CHAOS. Carlos E. Puente. Department of Land, Air and Water Resources University of California, Davis

THE ESSENCE OF CHAOS. Carlos E. Puente. Department of Land, Air and Water Resources University of California, Davis THE ESSENCE OF CHAOS Carlos E. Puente Department of Land, Air and Water Resources University of California, Davis http://puente.lawr.ucdavis.edu 2 Outline Introduces the logistic map and its incredible

More information

AP Calculus Testbank (Chapter 9) (Mr. Surowski)

AP Calculus Testbank (Chapter 9) (Mr. Surowski) AP Calculus Testbank (Chapter 9) (Mr. Surowski) Part I. Multiple-Choice Questions n 1 1. The series will converge, provided that n 1+p + n + 1 (A) p > 1 (B) p > 2 (C) p >.5 (D) p 0 2. The series

More information

Equations with regular-singular points (Sect. 5.5).

Equations with regular-singular points (Sect. 5.5). Equations with regular-singular points (Sect. 5.5). Equations with regular-singular points. s: Equations with regular-singular points. Method to find solutions. : Method to find solutions. Recall: The

More information

PERIODICITY AND CHAOS ON A MODIFIED SAMUELSON MODEL

PERIODICITY AND CHAOS ON A MODIFIED SAMUELSON MODEL PERIODICITY AND CHAOS ON A MODIFIED SAMUELSON MODEL Jose S. Cánovas Departamento de Matemática Aplicada y Estadística. Universidad Politécnica de Cartagena e-mail:jose.canovas@upct.es Manuel Ruiz Marín

More information

Maps and differential equations

Maps 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 information

MATH 2554 (Calculus I)

MATH 2554 (Calculus I) MATH 2554 (Calculus I) Dr. Ashley K. University of Arkansas February 21, 2015 Table of Contents Week 6 1 Week 6: 16-20 February 3.5 Derivatives as Rates of Change 3.6 The Chain Rule 3.7 Implicit Differentiation

More information

Symbolic Dynamics of One-Dimensional Maps: Entropies, Finite Precision, and Noise

Symbolic Dynamics of One-Dimensional Maps: Entropies, Finite Precision, and Noise International Journal of Theoretical Physics, VoL 21, Nos. 6/7, 1982 Symbolic Dynamics of One-Dimensional Maps: Entropies, Finite Precision, and Noise J. P. Crutchfield and N. H. Packard Physics Board

More information