Phase-Space Reconstruction. Gerrit Ansmann
|
|
- Brooke Rodgers
- 5 years ago
- Views:
Transcription
1 Phase-Space Reconstruction Gerrit Ansmann
2 Reprise: The Need for Non-Linear Methods. Lorenz oscillator x = 1(y x), y = x(28 z) y, z = xy 8z 3 Autoregressive process measured with non-linearity: y t =.8y t 1 + ε t x t = (1 + y t ) 2 2
3 Reprise: The Need for Non-Linear Methods. Lorenz oscillator x = 1(y x), y = x(28 z) y, z = xy 8z 3 Autoregressive process measured with non-linearity: y t =.8y t 1 + ε t x t = (1 + y t ) x t 1 1 x t t t 2
4 Reprise: The Need for Non-Linear Methods. Lorenz oscillator x = 1(y x), y = x(28 z) y, z = xy 8z 3 Skewness: v.4, Kurtosis: γ.71 Autoregressive process measured with non-linearity: y t =.8y t 1 + ε t x t = (1 + y t ) 2 Skewness: v 2.6 Kurtosis: γ 9.7 PDF (a.u.) PDF (a.u.) x t x t 2
5 Reprise: The Need for Non-Linear Methods. Lorenz oscillator x = 1(y x), y = x(28 z) y, z = xy 8z 3 Autoregressive process measured with non-linearity: y t =.8y t 1 + ε t x t = (1 + y t ) 2 spectral density (a.u) spectral density (a.u) period length period length 2
6 Reprise: The Need for Non-Linear Methods. Lorenz oscillator x = 1(y x), y = x(28 z) y, z = xy 8z 3 Autoregressive process measured with non-linearity: y t =.8y t 1 + ε t x t = (1 + y t ) R(x,τ) R(x,τ) τ τ 2
7 Reprise: The Need for Non-Linear Methods When faced with chaotic time series, linear methods fail to detect the dynamics / structure do not tell much about the dynamics cannot distinguish chaos from noise Structure can be seen in attractors. 3
8 Reminder: Attractors States of the dynamics for t. Type of dynamics can be deduced from topology of attractor: point fixed-point dynamics limit cycle periodic dynamics torus quasiperiodic dynamics strange attractor chaos Attractor reflects further central properties of dynamics. 4
9 Attractor of the Lorenz System adapted from commons.wikimedia.org/wiki/file:lorenz_system_r28_s1_b png 5
10 The Need for Phase-Space Reconstruction Directly observing the phase space/attractor requires access to all dynamical variables. But: Often, only one dynamical variable is available. The dimension of phase space is unknown. 6
11 The Need for Phase-Space Reconstruction Directly observing the phase space/attractor requires access to all dynamical variables. But: Often, only one dynamical variable is available. The dimension of phase space is unknown. Can we obtain from one time series a set that preserves important properties of the attractor? 6
12 The Need for Phase-Space Reconstruction Directly observing the phase space/attractor requires access to all dynamical variables. But: Often, only one dynamical variable is available. The dimension of phase space is unknown. Can we obtain from one time series a set that preserves important properties of the attractor? Yes, but it s a bit of work. 6
13 Good and Bad Reconstructions actual attractor perfect reconstruction perfect reconstruction preserves structure, practically bad preserves structure, practically bad 2D: structure not preserved 3D: structure may be preserved (phase as real number) structure not preserved (forth and back on the same line) structure not preserved 7
14 Reconstructions in the Language of Topology original attractor an n-manifold A R d measurement and reconstruction a map: Φ A R m structure-preserving reconstruction topology-preserving map an embedding 8
15 Reconstructions in the Language of Topology original attractor an n-manifold A R d measurement and reconstruction a map: Φ A R m structure-preserving reconstruction topology-preserving map an embedding. embedding. A map Φ A R m is called an embedding, if:. Φ has full rank. Φ is a diffeomorphism: Φ is differentiable Φ 1 exists and is differentiable 8
16 Whitney s Embedding Theorems. Strong Whitney embedding theorem..for m = 2n, there exists a Φ A R m that is an embedding. Problem: We do not know Φ. 9
17 Whitney s Embedding Theorems. Strong Whitney embedding theorem..for m = 2n, there exists a Φ A R m that is an embedding.. Weak Whitney embedding theorem. For m 2n + 1, almost every continuously differentiable (C 1 ) map. Φ A R m is an embedding. Problems: Often, we do not have m independent observables. (Redundant observables are one of the reasons for almost every.) We do not know n. 9
18 Delay Embedding: Idea Time series: x 1, x 2,, x N x t, x t, x t, are not fully redundant. Approximate: x t = x t +1 x t x t = x t +2 2x t +1 + x t x t, x t +1, are not fully redundant. inverse Taylor expansion. 1
19 Takens Delay-Embedding Theorem.. Takens theorem. Let A = {y, y 1,, y N } with the index indicating time.. y t is the state of the system at time t. 11
20 Takens Delay-Embedding Theorem.. Takens theorem. Let A = {y, y 1,, y N } with the index indicating time. Let h A R be the measurement function.. y t is the state of the system at time t. h(y ), h(y 1),, h(y N ) is the time series. 11
21 Takens Delay-Embedding Theorem.. Takens theorem. Let A = {y, y 1,, y N } with the index indicating time. Let h A R be the measurement function. If m 2n + 1,. Φ hτ(y t ) = (h(y t ), h(y t τ),, h(y t (m 1)τ)) is an embedding for almost all dynamics, τ, and h. y t is the state of the system at time t. h(y ), h(y 1),, h(y N ) is the time series. m is called embedding dimension. τ is called embedding delay. 11
22 Takens Phase-Space Reconstruction Once more, for application: Univariate time series x: x 1, x 2,, x N Consider m-dimensional states (on the attractor): x t x t τ x t 2τ x t (m 1)τ For a proper embedding dimension m and embedding delay τ, these states make up a topologically equivalent reconstruction of the attractor. 12
23 Example: Sine Wave / Harmonic Oscillator Blackboard. 13
24 Example: Lorenz Attractor original reconstruction (τ =.8)
25 U(t) EEG (Wachzustand) t x(t) x(t+τ)
26 U(t) EEG (Epilepsiepatient) t x(t) x(t+τ)
27 U(t) EEG (epileptischer Anfall) t x(t) x(t+τ)
28 Theorems of Sauer, Yorke, and Casdagli Problem: Strange Attractors are only a subset of some n-manifold.. Fractal embedding theorems (Sauer, Yorke, and Casdagli). In the weak Whitney embedding theorem and Takens theorem, one can replace:. n-manifold A compact invariant set A n the box-counting dimension d b of A almost all (open and dense subset) almost all (probability 1) Advantage: d b < n 15
29 Dynamical Invariants Important characteristics of the dynamics are invariant under the embedding transformation: Lyapunov exponents dimensions entropy 16
30 Impact of the Delay: Sine Wave / Harmonic Oscillator Blackboard. 17
31 Impact of the Delay: Lorenz Attractor 2 τ =
32 Impact of the Delay: Lorenz Attractor 2 τ =
33 Impact of the Delay: Lorenz Attractor 2 τ =
34 Impact of the Delay: Lorenz Attractor 2 τ =
35 Impact of the Delay: Lorenz Attractor 2 τ =
36 Impact of the Delay: Lorenz Attractor 2 τ =
37 Impact of the Delay: Lorenz Attractor 2 τ =
38 Impact of the Delay: Lorenz Attractor 2 τ =
39 Impact of the Delay: Lorenz Attractor 2 τ =
40 Impact of the Delay: Lorenz Attractor 2 τ =
41 Impact of the Delay: Lorenz Attractor 2 τ =
42 Impact of the Delay: Lorenz Attractor 2 τ =
43 Impact of the Delay: Lorenz Attractor 2 τ =
44 Impact of the Delay: Lorenz Attractor 2 τ =
45 Impact of the Delay: Lorenz Attractor 2 τ =
46 Impact of the Delay: Lorenz Attractor 2 τ =
47 Impact of the Delay: Lorenz Attractor 2 τ =
48 Impact of the Delay: Lorenz Attractor 2 τ =
49 Impact of the Delay: Lorenz Attractor 2 τ =
50 Impact of the Delay: Lorenz Attractor 2 τ =
51 Impact of the Delay: Lorenz Attractor 2 τ =
52 Choosing the Time Delay Requirement for an embedding: x t, x t τ, x t 2τ,, x t (m 1)τ not fully redundant. Aforementioned theorems: Almost every τ yields an embedding: Requirement for a good embedding: Minimal redundance of x t, x t τ, x t 2τ,, x t (m 1)τ (to unfold the attractor) Reasonably small τ (to avoid folding the attractor onto itself) Compare to: linear independence vs. orthogonality. 19
53 Zeros of the Autocorrelation Function Idea: R(x,τ) = x t and x t τ linearly independent on average. Choose τ as the first zero of the autocorrelation function. 2
54 First Minimum of the Mutual Information Mutual information I m (τ): Common information (redundancy) between x t and x t τ. To what extent does x t τ help to predict x t? Measure of independence / redundancy based on information theory (non-linear). (More on this in a few weeks) Choose τ as the first minimum of the mutual-information function. A. Fraser and H. Swinney, Independent coordinates for strange attractors from mutual information, Phys. Rev. A 33 (1986) 21
55 Autocorrelation vs. Mutual Information from A. Fraser and H. Swinney, Independent coordinates for strange attractors from mutual information, Phys. Rev. A 33 (1986) 22
56 Methods to Find the Embedding Delay zeros of the autocorrelation function minima of the mutual information many more No method is perfect or commonly agreed upon. Practically: Try at least two methods. Judge by further analysis. Alternative for m 3: Visually inspect the attractor. Keep (m 1)τ (time span in an embedded vector) constant. 23
57 Choosing the Embedding Dimension Aforementioned theorems: m = 2n + 1 or m = 2d b + 1 suffices. But: A priori, we do know neither n nor d b. This is the worst case a smaller m may also suffice: The surface of a torus (n = 2) can be embedded with m = 3. The Lorenz attractor (n = 3, d b 2) can be embedded with m = 3. Choosing m overly high may hamper further analyses (impact of noise, computational complexity). Need other ways to determine a good m. 24
58 Linear Dependence Idea: m is higher than necessary attractor only covers a subspace of reconstruction space (e.g, circle in m = 3) Check whether embedded vectors have full rank. Problems: Noise acts in all directions. Assumes linear dependence dependence may be non-linear. (e.g., bent circle in m = 3) 25
59 Assymptotic Invariants m too small wrong dynamical invariants (in general) m sufficient correct dynamical invariants Increase m until dynamical invariants converge. Difficulties: Criterion for convergence under real conditions (noise, finite data, ). Wrong (in general) is unpredictable. 26
60 False Nearest Neighbours Idea m too small trajectories intersect points close in reconstruction space that aren t close in actual phase space (false nearest neigbours) Raise m until false nearest neighbours vanish. M. B. Kennel, R. Brown, and H. D. I. Abarbanel, Determining minimum embedding dimension using a geometrical construction, Phys. Rev. A 45 (1992) 27
61 False Nearest Neighbours Practically Practically: Choose threshold ε for nearest neighbours. NN(m): Number of pairs of points in m-dimensional reconstruction space that are closer than ε. NN(m + 1) < NN(m) at least NN(m) NN(m + 1) false nearest neighbours in the m-dimensional reconstruction. Difficulties: Number of true nearest neighbours large and fluctuating (noise). 28
62 Methods to Determine the Embedding Dimension assymptotic invariants false nearest neighbours many more No method is perfect or commonly agreed upon. Practically: Try at least two methods. Judge by further analysis. Alternative for m 3: Visually inspect the attractor. If possible: embedding theorems as a sanity check. Keep (m 1)τ (time span in an embedded vector) constant. 29
63 Summary Delay embedding allows to reconstruct attractor from single observable. Parameters m and τ have to be carefully chosen. Reconstructed phase space may be used for: understanding prediction modelling Characteristics preserved by reconstruction: Lyapunov exponents dimensions entropy 3
64 Extensions multivariate time series different embedding delays for each component state-dependent embedding delays 31
Delay Coordinate Embedding
Chapter 7 Delay Coordinate Embedding Up to this point, we have known our state space explicitly. But what if we do not know it? How can we then study the dynamics is phase space? A typical case is when
More informationEstimating Lyapunov Exponents from Time Series. Gerrit Ansmann
Estimating Lyapunov Exponents from Time Series Gerrit Ansmann Example: Lorenz system. Two identical Lorenz systems with initial conditions; one is slightly perturbed (10 14 ) at t = 30: 20 x 2 1 0 20 20
More informationDynamical systems, information and time series
Dynamical systems, information and time series Stefano Marmi Scuola Normale Superiore http://homepage.sns.it/marmi/ Lecture 4- European University Institute October 27, 2009 Lecture 1: An introduction
More informationDynamics of partial discharges involved in electrical tree growth in insulation and its relation with the fractal dimension
Dynamics of partial discharges involved in electrical tree growth in insulation and its relation with the fractal dimension Daniela Contreras Departamento de Matemática Universidad Técnica Federico Santa
More informationCharacterizing chaotic time series
Characterizing chaotic time series Jianbo Gao PMB InTelliGence, LLC, West Lafayette, IN 47906 Mechanical and Materials Engineering, Wright State University jbgao.pmb@gmail.com http://www.gao.ece.ufl.edu/
More informationComplex system approach to geospace and climate studies. Tatjana Živković
Complex system approach to geospace and climate studies Tatjana Živković 30.11.2011 Outline of a talk Importance of complex system approach Phase space reconstruction Recurrence plot analysis Test for
More informationNo. 6 Determining the input dimension of a To model a nonlinear time series with the widely used feed-forward neural network means to fit the a
Vol 12 No 6, June 2003 cfl 2003 Chin. Phys. Soc. 1009-1963/2003/12(06)/0594-05 Chinese Physics and IOP Publishing Ltd Determining the input dimension of a neural network for nonlinear time series prediction
More informationNONLINEAR TIME SERIES ANALYSIS, WITH APPLICATIONS TO MEDICINE
NONLINEAR TIME SERIES ANALYSIS, WITH APPLICATIONS TO MEDICINE José María Amigó Centro de Investigación Operativa, Universidad Miguel Hernández, Elche (Spain) J.M. Amigó (CIO) Nonlinear time series analysis
More informationTHREE DIMENSIONAL SYSTEMS. Lecture 6: The Lorenz Equations
THREE DIMENSIONAL SYSTEMS Lecture 6: The Lorenz Equations 6. The Lorenz (1963) Equations The Lorenz equations were originally derived by Saltzman (1962) as a minimalist model of thermal convection in a
More informationDetection of Nonlinearity and Stochastic Nature in Time Series by Delay Vector Variance Method
International Journal of Engineering & Technology IJET-IJENS Vol:10 No:02 11 Detection of Nonlinearity and Stochastic Nature in Time Series by Delay Vector Variance Method Imtiaz Ahmed Abstract-- This
More informationAnalysis of Electroencephologram Data Using Time-Delay Embeddings to Reconstruct Phase Space
Dynamics at the Horsetooth Volume 1, 2009. Analysis of Electroencephologram Data Using Time-Delay Embeddings to Reconstruct Phase Space Department of Mathematics Colorado State University Report submitted
More informationSIMULATED CHAOS IN BULLWHIP EFFECT
Journal of Management, Marketing and Logistics (JMML), ISSN: 2148-6670 Year: 2015 Volume: 2 Issue: 1 SIMULATED CHAOS IN BULLWHIP EFFECT DOI: 10.17261/Pressacademia.2015111603 Tunay Aslan¹ ¹Sakarya University,
More informationChaos, Complexity, and Inference (36-462)
Chaos, Complexity, and Inference (36-462) Lecture 4 Cosma Shalizi 22 January 2009 Reconstruction Inferring the attractor from a time series; powerful in a weird way Using the reconstructed attractor to
More informationRevista Economica 65:6 (2013)
INDICATIONS OF CHAOTIC BEHAVIOUR IN USD/EUR EXCHANGE RATE CIOBANU Dumitru 1, VASILESCU Maria 2 1 Faculty of Economics and Business Administration, University of Craiova, Craiova, Romania 2 Faculty of Economics
More informationRemoving the Noise from Chaos Plus Noise
Removing the Noise from Chaos Plus Noise Steven P. Lalley Department of Statistics University of Chicago November 5, 2 Abstract The problem of extracting a signal x n generated by a dynamical system from
More informationINTRODUCTION TO CHAOS THEORY T.R.RAMAMOHAN C-MMACS BANGALORE
INTRODUCTION TO CHAOS THEORY BY T.R.RAMAMOHAN C-MMACS BANGALORE -560037 SOME INTERESTING QUOTATIONS * PERHAPS THE NEXT GREAT ERA OF UNDERSTANDING WILL BE DETERMINING THE QUALITATIVE CONTENT OF EQUATIONS;
More informationNonlinear Prediction for Top and Bottom Values of Time Series
Vol. 2 No. 1 123 132 (Feb. 2009) 1 1 Nonlinear Prediction for Top and Bottom Values of Time Series Tomoya Suzuki 1 and Masaki Ota 1 To predict time-series data depending on temporal trends, as stock price
More informationA benchmark test A benchmark test of accuracy and precision in estimating dynamical systems characteristics from a time series
2 A benchmark test A benchmark test of accuracy and precision in estimating dynamical systems characteristics from a time series, SM Rispens, M Pijnappels, JH van Dieën, KS van Schooten, PJ Beek, A Daffertshofer
More informationarxiv: v1 [nlin.ao] 21 Sep 2018
Using reservoir computers to distinguish chaotic signals T L Carroll US Naval Research Lab, Washington, DC 20375 (Dated: October 11, 2018) Several recent papers have shown that reservoir computers are
More informationMechanisms of Chaos: Stable Instability
Mechanisms of Chaos: Stable Instability Reading for this lecture: NDAC, Sec. 2.-2.3, 9.3, and.5. Unpredictability: Orbit complicated: difficult to follow Repeatedly convergent and divergent Net amplification
More informationMajid Sodagar, 1 Patrick Chang, 1 Edward Coyler, 1 and John Parke 1 School of Physics, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
Experimental Characterization of Chua s Circuit Majid Sodagar, 1 Patrick Chang, 1 Edward Coyler, 1 and John Parke 1 School of Physics, Georgia Institute of Technology, Atlanta, Georgia 30332, USA (Dated:
More informationInvestigation of Chaotic Nature of Sunspot Data by Nonlinear Analysis Techniques
American-Eurasian Journal of Scientific Research 10 (5): 272-277, 2015 ISSN 1818-6785 IDOSI Publications, 2015 DOI: 10.5829/idosi.aejsr.2015.10.5.1150 Investigation of Chaotic Nature of Sunspot Data by
More information1 Random walks and data
Inference, Models and Simulation for Complex Systems CSCI 7-1 Lecture 7 15 September 11 Prof. Aaron Clauset 1 Random walks and data Supposeyou have some time-series data x 1,x,x 3,...,x T and you want
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 informationWHAT 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 informationMathematical Foundations of Neuroscience - Lecture 7. Bifurcations II.
Mathematical Foundations of Neuroscience - Lecture 7. Bifurcations II. Filip Piękniewski Faculty of Mathematics and Computer Science, Nicolaus Copernicus University, Toruń, Poland Winter 2009/2010 Filip
More informationModeling and Predicting Chaotic Time Series
Chapter 14 Modeling and Predicting Chaotic Time Series To understand the behavior of a dynamical system in terms of some meaningful parameters we seek the appropriate mathematical model that captures the
More informationSimple approach to the creation of a strange nonchaotic attractor in any chaotic system
PHYSICAL REVIEW E VOLUME 59, NUMBER 5 MAY 1999 Simple approach to the creation of a strange nonchaotic attractor in any chaotic system J. W. Shuai 1, * and K. W. Wong 2, 1 Department of Biomedical Engineering,
More informationPHY411 Lecture notes Part 5
PHY411 Lecture notes Part 5 Alice Quillen January 27, 2016 Contents 0.1 Introduction.................................... 1 1 Symbolic Dynamics 2 1.1 The Shift map.................................. 3 1.2
More informationMultifractal Models for Solar Wind Turbulence
Multifractal Models for Solar Wind Turbulence Wiesław M. Macek Faculty of Mathematics and Natural Sciences. College of Sciences, Cardinal Stefan Wyszyński University, Dewajtis 5, 01-815 Warsaw, Poland;
More informationEnhanced sensitivity of persistent events to weak forcing in dynamical and stochastic systems: Implications for climate change. Khatiwala, et.al.
Enhanced sensitivity of persistent events to weak forcing in dynamical and stochastic systems: Implications for climate change Questions What are the characteristics of the unforced Lorenz system? What
More informationChaos and Liapunov exponents
PHYS347 INTRODUCTION TO NONLINEAR PHYSICS - 2/22 Chaos and Liapunov exponents Definition of chaos In the lectures we followed Strogatz and defined chaos as aperiodic long-term behaviour in a deterministic
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 informationIntroduction Knot Theory Nonlinear Dynamics Topology in Chaos Open Questions Summary. Topology in Chaos
Introduction Knot Theory Nonlinear Dynamics Open Questions Summary A tangled tale about knot, link, template, and strange attractor Centre for Chaos & Complex Networks City University of Hong Kong Email:
More informationThe Behaviour of a Mobile Robot Is Chaotic
AISB Journal 1(4), c SSAISB, 2003 The Behaviour of a Mobile Robot Is Chaotic Ulrich Nehmzow and Keith Walker Department of Computer Science, University of Essex, Colchester CO4 3SQ Department of Physics
More informationCHAPTER 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 informationForecasting Financial Indices: The Baltic Dry Indices
International Journal of Maritie, Trade & Econoic Issues pp. 109-130 Volue I, Issue (1), 2013 Forecasting Financial Indices: The Baltic Dry Indices Eleftherios I. Thalassinos 1, Mike P. Hanias 2, Panayiotis
More informationSmooth Local Subspace Projection for Nonlinear Noise Reduction
University of Rhode Island DigitalCommons@URI Mechanical, Industrial & Systems Engineering Faculty Publications Mechanical, Industrial & Systems Engineering 23 Smooth Local Subspace Projection for Nonlinear
More informationStabilization of Hyperbolic Chaos by the Pyragas Method
Journal of Mathematics and System Science 4 (014) 755-76 D DAVID PUBLISHING Stabilization of Hyperbolic Chaos by the Pyragas Method Sergey Belyakin, Arsen Dzanoev, Sergey Kuznetsov Physics Faculty, Moscow
More informationMore Details Fixed point of mapping is point that maps into itself, i.e., x n+1 = x n.
More Details Fixed point of mapping is point that maps into itself, i.e., x n+1 = x n. If there are points which, after many iterations of map then fixed point called an attractor. fixed point, If λ
More informationApril 13, We now extend the structure of the horseshoe to more general kinds of invariant. x (v) λ n v.
April 3, 005 - Hyperbolic Sets We now extend the structure of the horseshoe to more general kinds of invariant sets. Let r, and let f D r (M) where M is a Riemannian manifold. A compact f invariant set
More informationDeterminism, Complexity, and Predictability in Computer Performance
Determinism, Complexity, and Predictability in Computer Performance arxiv:1305.508v1 [nlin.cd] 23 May 2013 Joshua Garland, Ryan G. James and Elizabeth Bradley Dept. of Computer Science University of Colorado,
More informationHow much information is contained in a recurrence plot?
Physics Letters A 330 (2004) 343 349 www.elsevier.com/locate/pla How much information is contained in a recurrence plot? Marco Thiel, M. Carmen Romano, Jürgen Kurths University of Potsdam, Nonlinear Dynamics,
More informationDynamical Systems and Chaos Part I: Theoretical Techniques. Lecture 4: Discrete systems + Chaos. Ilya Potapov Mathematics Department, TUT Room TD325
Dynamical Systems and Chaos Part I: Theoretical Techniques Lecture 4: Discrete systems + Chaos Ilya Potapov Mathematics Department, TUT Room TD325 Discrete maps x n+1 = f(x n ) Discrete time steps. x 0
More informationChaos. Dr. Dylan McNamara people.uncw.edu/mcnamarad
Chaos Dr. Dylan McNamara people.uncw.edu/mcnamarad Discovery of chaos Discovered in early 1960 s by Edward N. Lorenz (in a 3-D continuous-time model) Popularized in 1976 by Sir Robert M. May as an example
More informationComplex Dynamics of Microprocessor Performances During Program Execution
Complex Dynamics of Microprocessor Performances During Program Execution Regularity, Chaos, and Others Hugues BERRY, Daniel GRACIA PÉREZ, Olivier TEMAM Alchemy, INRIA, Orsay, France www-rocq.inria.fr/
More informationphysics/ v2 21 Jun 1999
TIME SERIES FORECASTING: A NONLINEAR DYNAMICS APPROACH Stefano Sello Termo-Fluid Dynamics Research Center Enel Research Via Andrea Pisano, 120 56122 PISA - ITALY e-mail: sello@pte.enel.it physics/9906035
More informationPeriodic Sinks and Observable Chaos
Periodic Sinks and Observable Chaos Systems of Study: Let M = S 1 R. T a,b,l : M M is a three-parameter family of maps defined by where θ S 1, r R. θ 1 = a+θ +Lsin2πθ +r r 1 = br +blsin2πθ Outline of Contents:
More information6.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 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 informationHowever, in actual topology a distance function is not used to define open sets.
Chapter 10 Dimension Theory We are used to the notion of dimension from vector spaces: dimension of a vector space V is the minimum number of independent bases vectors needed to span V. Therefore, a point
More informationPHONEME CLASSIFICATION OVER THE RECONSTRUCTED PHASE SPACE USING PRINCIPAL COMPONENT ANALYSIS
PHONEME CLASSIFICATION OVER THE RECONSTRUCTED PHASE SPACE USING PRINCIPAL COMPONENT ANALYSIS Jinjin Ye jinjin.ye@mu.edu Michael T. Johnson mike.johnson@mu.edu Richard J. Povinelli richard.povinelli@mu.edu
More informationFoliations of Three Dimensional Manifolds
Foliations of Three Dimensional Manifolds M. H. Vartanian December 17, 2007 Abstract The theory of foliations began with a question by H. Hopf in the 1930 s: Does there exist on S 3 a completely integrable
More informationfy (X(g)) Y (f)x(g) gy (X(f)) Y (g)x(f)) = fx(y (g)) + gx(y (f)) fy (X(g)) gy (X(f))
1. Basic algebra of vector fields Let V be a finite dimensional vector space over R. Recall that V = {L : V R} is defined to be the set of all linear maps to R. V is isomorphic to V, but there is no canonical
More informationTakens embedding theorem for infinite-dimensional dynamical systems
Takens embedding theorem for infinite-dimensional dynamical systems James C. Robinson Mathematics Institute, University of Warwick, Coventry, CV4 7AL, U.K. E-mail: jcr@maths.warwick.ac.uk Abstract. Takens
More informationAuthor(s) Okamura, Noriaki; Okamoto, Hisao; T.
Title A New Classification of Strange Att Mutual Information Author(s) Okamura, Noriaki; Okamoto, Hisao; T Chaos Memorial Symposium in Asuka : Citation dedicated to professor Yoshisuke Ue 60th birthday,
More informationTesting for Chaos in Type-I ELM Dynamics on JET with the ILW. Fabio Pisano
Testing for Chaos in Type-I ELM Dynamics on JET with the ILW Fabio Pisano ACKNOWLEDGMENTS B. Cannas 1, A. Fanni 1, A. Murari 2, F. Pisano 1 and JET Contributors* EUROfusion Consortium, JET, Culham Science
More informationThe 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 informationA topological delay embedding theorem for infinite-dimensional dynamical systems
INSTITUTE OF PHYSICS PUBLISHING Nonlinearity 18 (2005) 2135 2143 NONLINEARITY doi:10.1088/0951-7715/18/5/013 A topological delay embedding theorem for infinite-dimensional dynamical systems 1. Introduction
More informationInformation Dynamics Foundations and Applications
Gustavo Deco Bernd Schürmann Information Dynamics Foundations and Applications With 89 Illustrations Springer PREFACE vii CHAPTER 1 Introduction 1 CHAPTER 2 Dynamical Systems: An Overview 7 2.1 Deterministic
More informationarxiv: v1 [nlin.cd] 7 Mar 2019
PRE Impact of Embedding View on Cross Mapping Convergence Robert Martin and Justin Koo In-Space Propulsion Branch, Air Force Research Laboratory Daniel Eckhardt Postdoctoral Fellow arxiv:93.369v [nlin.cd]
More informationIntroduction to Classical Chaos
Introduction to Classical Chaos WeiHan Hsiao a a Department of Physics, The University of Chicago E-mail: weihanhsiao@uchicago.edu ABSTRACT: This note is a contribution to Kadanoff Center for Theoretical
More informationNonlinear dynamics, delay times, and embedding windows
Physica D 127 (1999) 48 60 Nonlinear dynamics, delay times, and embedding windows H.S. Kim 1,a, R. Eykholt b,, J.D. Salas c a Department of Civil Engineering, Colorado State University, Fort Collins, CO
More informationInformation Mining for Friction Torque of Rolling Bearing for Space Applications Using Chaotic Theory
Research Journal of Applied Sciences, Engineering and Technology 5(22): 5223229, 213 ISSN: 24-7459; e-issn: 24-7467 Maxwell Scientific Organization, 213 Submitted: October 9, 212 Accepted: December 3,
More informationMultivariable Calculus
2 Multivariable Calculus 2.1 Limits and Continuity Problem 2.1.1 (Fa94) Let the function f : R n R n satisfy the following two conditions: (i) f (K ) is compact whenever K is a compact subset of R n. (ii)
More informationFocus was on solving matrix inversion problems Now we look at other properties of matrices Useful when A represents a transformations.
Previously Focus was on solving matrix inversion problems Now we look at other properties of matrices Useful when A represents a transformations y = Ax Or A simply represents data Notion of eigenvectors,
More informationSmooth Structure. lies on the boundary, then it is determined up to the identifications it 1 2
132 3. Smooth Structure lies on the boundary, then it is determined up to the identifications 1 2 + it 1 2 + it on the vertical boundary and z 1/z on the circular part. Notice that since z z + 1 and z
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 informationPredictability and Chaotic Nature of Daily Streamflow
34 th IAHR World Congress - Balance and Uncertainty 26 June - 1 July 211, Brisbane, Australia 33 rd Hydrology & Water Resources Symposium 1 th Hydraulics Conference Predictability and Chaotic Nature of
More informationExperimental Characterization of Nonlinear Dynamics from Chua s Circuit
Experimental Characterization of Nonlinear Dynamics from Chua s Circuit John Parker*, 1 Majid Sodagar, 1 Patrick Chang, 1 and Edward Coyle 1 School of Physics, Georgia Institute of Technology, Atlanta,
More informationUnique equilibrium states for geodesic flows in nonpositive curvature
Unique equilibrium states for geodesic flows in nonpositive curvature Todd Fisher Department of Mathematics Brigham Young University Fractal Geometry, Hyperbolic Dynamics and Thermodynamical Formalism
More informationUSING DELAY-COORDINATE EMBEDDING AS A DATA MINING TECHNIQUE: PART ONE
James Popoff, Ph.D. jamespopoffphd@gmail.com 561-371-0349 USIG DELAY-COORDIATE EMBEDDIG AS A DATA MIIG TECHIQUE: PART OE Data mining is an example of the theory of counting in applied form. The theory
More informationA MINIMAL 2-D QUADRATIC MAP WITH QUASI-PERIODIC ROUTE TO CHAOS
International Journal of Bifurcation and Chaos, Vol. 18, No. 5 (2008) 1567 1577 c World Scientific Publishing Company A MINIMAL 2-D QUADRATIC MAP WITH QUASI-PERIODIC ROUTE TO CHAOS ZERAOULIA ELHADJ Department
More informationApplication of Physics Model in prediction of the Hellas Euro election results
Journal of Engineering Science and Technology Review 2 (1) (2009) 104-111 Research Article JOURNAL OF Engineering Science and Technology Review www.jestr.org Application of Physics Model in prediction
More informationProblems in hyperbolic dynamics
Problems in hyperbolic dynamics Current Trends in Dynamical Systems and the Mathematical Legacy of Rufus Bowen Vancouver july 31st august 4th 2017 Notes by Y. Coudène, S. Crovisier and T. Fisher 1 Zeta
More informationTHE EULER CHARACTERISTIC OF A LIE GROUP
THE EULER CHARACTERISTIC OF A LIE GROUP JAY TAYLOR 1 Examples of Lie Groups The following is adapted from [2] We begin with the basic definition and some core examples Definition A Lie group is a smooth
More informationAttractor of a Shallow Water Equations Model
Thai Journal of Mathematics Volume 5(2007) Number 2 : 299 307 www.math.science.cmu.ac.th/thaijournal Attractor of a Shallow Water Equations Model S. Sornsanam and D. Sukawat Abstract : In this research,
More informationWhat is Chaos? Implications of Chaos 4/12/2010
Joseph Engler Adaptive Systems Rockwell Collins, Inc & Intelligent Systems Laboratory The University of Iowa When we see irregularity we cling to randomness and disorder for explanations. Why should this
More informationNumerics of Stochastic Processes II ( ) Illia Horenko
Numerics of Stochastic Processes II (14.11.2008) Illia Horenko Research Group Computational Time Series Analysis Institute of Mathematics Freie Universität Berlin (FU) DFG Research Center MATHEON Mathematics
More informationFAULT DETECTION IN DYNAMIC SYSTEMS USING THE LARGEST LYAPUNOV EXPONENT
FAULT DETECTION IN DYNAMIC SYSTEMS USING THE LARGEST LYAPUNOV EXPONENT A Thesis by YIFU SUN Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the requirements
More informationOn the choice of coarse variables for. dynamics. Workshop on Mathematics of Model Reduction University of Leicester, Aug.
On the choice of coarse variables for dynamics Amit Acharya Workshop on Mathematics of Model Reduction University of Leicester, Aug. 28-30, 2007 Acknowledgment: US DOE Locally Invariant Manifolds System
More informationThe dynamics of human gait
INSTITUTE OF PHYSICS PUBLISHING Eur. J. Phys. 26 (2005) 525 534 EUROPEAN JOURNAL OF PHYSICS doi:10.1088/0143-0807/26/3/017 The dynamics of human gait Matjaž Perc Department of Physics, Faculty of Education,
More informationDefending against Internet worms using a phase space method from chaos theory
Defending against Internet worms using a phase space method from chaos theory Jing Hu and Jianbo Gao Department of Electrical & Computer Engineering University of Florida Nageswara S. Rao Computer Science
More informationEssential hyperbolicity versus homoclinic bifurcations. Global dynamics beyond uniform hyperbolicity, Beijing 2009 Sylvain Crovisier - Enrique Pujals
Essential hyperbolicity versus homoclinic bifurcations Global dynamics beyond uniform hyperbolicity, Beijing 2009 Sylvain Crovisier - Enrique Pujals Generic dynamics Consider: M: compact boundaryless manifold,
More informationNonlinear dynamics of the Great Salt Lake: system identification and prediction
Climate Dynamics (1996) 12: 287 297 q Springer-Verlag 1996 Nonlinear dynamics of the Great Salt Lake: system identification and prediction Henry D. I. Abarbanel 1, Upmanu Lall 2 1 Department of Physics
More informationThe Big Picture. Python environment? Discuss Examples of unpredictability. Chaos, Scientific American (1986)
The Big Picture Python environment? Discuss Examples of unpredictability Email homework to me: chaos@cse.ucdavis.edu Chaos, Scientific American (1986) Odds, Stanislaw Lem, The New Yorker (1974) 1 Nonlinear
More informationDetermination of fractal dimensions of solar radio bursts
Astron. Astrophys. 357, 337 350 (2000) ASTRONOMY AND ASTROPHYSICS Determination of fractal dimensions of solar radio bursts A. Veronig 1, M. Messerotti 2, and A. Hanslmeier 1 1 Institute of Astronomy,
More informationLecture 6. Lorenz equations and Malkus' waterwheel Some properties of the Lorenz Eq.'s Lorenz Map Towards definitions of:
Lecture 6 Chaos Lorenz equations and Malkus' waterwheel Some properties of the Lorenz Eq.'s Lorenz Map Towards definitions of: Chaos, Attractors and strange attractors Transient chaos Lorenz Equations
More informationy 2 . = x 1y 1 + x 2 y x + + x n y n 2 7 = 1(2) + 3(7) 5(4) = 3. x x = x x x2 n.
6.. Length, Angle, and Orthogonality In this section, we discuss the defintion of length and angle for vectors and define what it means for two vectors to be orthogonal. Then, we see that linear systems
More informationDimension estimation for autonomous nonlinear systems
Dimension estimation for autonomous nonlinear systems Alberto Padoan and Alessandro Astolfi Abstract The problem of estimating the dimension of the state-space of an autonomous nonlinear systems is considered
More informationTwo Decades of Search for Chaos in Brain.
Two Decades of Search for Chaos in Brain. A. Krakovská Inst. of Measurement Science, Slovak Academy of Sciences, Bratislava, Slovak Republic, Email: krakovska@savba.sk Abstract. A short review of applications
More informationStatistical Pattern Recognition
Statistical Pattern Recognition Feature Extraction Hamid R. Rabiee Jafar Muhammadi, Alireza Ghasemi, Payam Siyari Spring 2014 http://ce.sharif.edu/courses/92-93/2/ce725-2/ Agenda Dimensionality Reduction
More informationFRACTAL GEOMETRY, DYNAMICAL SYSTEMS AND CHAOS
FRACTAL GEOMETRY, DYNAMICAL SYSTEMS AND CHAOS MIDTERM SOLUTIONS. Let f : R R be the map on the line generated by the function f(x) = x 3. Find all the fixed points of f and determine the type of their
More informationInterdependence and Predictability of Human Mobility and Social Interactions
Interdependence and Predictability of Human Mobility and Social Interactions Manlio De Domenico School of Computer Science University of Birmingham United Kingdom m.dedomenico@cs.bham.ac.uk Antonio Lima
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 informationBROUWER FIXED POINT THEOREM. Contents 1. Introduction 1 2. Preliminaries 1 3. Brouwer fixed point theorem 3 Acknowledgments 8 References 8
BROUWER FIXED POINT THEOREM DANIELE CARATELLI Abstract. This paper aims at proving the Brouwer fixed point theorem for smooth maps. The theorem states that any continuous (smooth in our proof) function
More informationNON-LINEAR DYNAMICS TOOLS FOR THE MOTION ANALYSIS AND CONDITION MONITORING OF ROBOT JOINTS
Mechanical Systems and Signal Processing (2001) 15(6), 1141}1164 doi:10.1006/mssp.2000.1394, available online at http://www.idealibrary.com on NON-LINEAR DYNAMICS TOOLS FOR THE MOTION ANALYSIS AND CONDITION
More informationImplicit Functions, Curves and Surfaces
Chapter 11 Implicit Functions, Curves and Surfaces 11.1 Implicit Function Theorem Motivation. In many problems, objects or quantities of interest can only be described indirectly or implicitly. It is then
More informationUnderstanding the dynamics of snowpack in Washington State - part II: complexity of
Understanding the dynamics of snowpack in Washington State - part II: complexity of snowpack Nian She and Deniel Basketfield Seattle Public Utilities, Seattle, Washington Abstract. In part I of this study
More informationEngineering Jump Start - Summer 2012 Mathematics Worksheet #4. 1. What is a function? Is it just a fancy name for an algebraic expression, such as
Function Topics Concept questions: 1. What is a function? Is it just a fancy name for an algebraic expression, such as x 2 + sin x + 5? 2. Why is it technically wrong to refer to the function f(x), for
More information