LYAPUNOV EXPONENT AND DIMENSIONS OF THE ATTRACTOR FOR TWO DIMENSIONAL NEURAL MODEL
|
|
- Marilyn Greene
- 5 years ago
- Views:
Transcription
1 Volume 1, No. 7, July 2013 Journal of Global Research in Mathematical Archives RESEARCH PAPER Available online at LYAPUNOV EXPONENT AND DIMENSIONS OF THE ATTRACTOR FOR TWO DIMENSIONAL NEURAL MODEL T.K.Dutta 1, A. Kr. Jain 2, D.Bhattacharjee 3 1 Professor (Dr.) Tarini Kumar Dutta, Department of Mathematics, Gauhati University, Guwahati,Assam, Pin code , India. ( tkdutta2001@yahoo.co.in) 2 Anil Kumar Jain, Assistant Professor, Department of Mathematics, Barama College,Barama,Assam,Pincode ,India ( jainanil965@gmail.com) 3 Debasish Bhattacharjee, Assistant Professor, Department of Mathematics, B.Borooah College, Guwahati, Assam, Pin code , India. ( debabh2@gmail.com) Corresponding author: Anil Kumar Jain, Assistant Professor, Department of Mathematics, Barama College, Barama, Assam, Pin code , India. ( jainanil965@gmail.com) Abstract: In this paper a two dimensional non linear neural network model is considered and it is shown that chaotic attractor exists beyond accumulation point. To confirm the existence of chaotic attractor, Lyapunov exponent method is used. Further various fractal dimensions like Correlation dimension, Box-counting and Information dimension of the chaotic attractor were found to assess the geometry of the fractal set. Keywords: Lyapunov Exponent, Strange attractor, fractal dimension, Correlation dimension, Box-counting and Information dimension. 1. INTRODUCTION: In [3] we have seen that two dimensional discrete neural model follows the period doubling bifurcation route as the control parameter vary and reaches from order to chaos beyond some particular value(accumulation point) of the control parameter. It is also seen that with the help of Feigenbaum universal constant, the accumulation point is calculated and it is given as The Lyapunov exponents play a crucial role in the description of the behaviour of dynamical systems. They measure the average rate of divergence or convergence of orbits starting from nearby initial points. Therefore, they can be used to analyse the stability of limits sets and check sensitive dependence on initial conditions, i.e. the presence of chaotic attractors and we see that chaos begins beyond accumulation point [3,11] While the main purpose of Lyapunov exponents is to characterise the dynamical properties of trajectories on attractors, the fractal dimensions focus on the geometry of the chaotic attractor. Chaotic dynamical systems exhibit trajectories in their phase space that converges to a strange attractor. The fractal dimensions of this strange attractor counts the effective number of degree of freedom in the dynamical system and thus quantifies its complexity. [1, 5, 12]. JGRMA 2012, All Rights Reserved 50
2 Various fractal dimensions focus on the geometric structure of attractors in state space which occurs beyond accumulation point. If the system has atleast one positive Lyapunov exponent, we say that the system s behaviour is chaotic. If the attractor exhibits scaling with (in general) a noninteger dimension (scaling index), then we say the attractor is strange. In this paper we calculate maximum lyapunov exponent(mle) of the model and it is observed that chaos occurs beyond accumulation point where mle is positive.we also study various fractal dimensions like Correlation dimension, Box-counting and Information dimension of the chaotic attractor at the onset of chaos, i.e. at the accumulation point (i.e ) [3] which are found to be a non-integer dimension. 2. LYAPUNOV EXPONENTS [4, 11]: Lyapunov exponents provide a quantitative measure of the divergence or convergence of the nearby trajectories for a dynamical system. If we consider a small hypersphere of initial conditions in the phase space, for sufficiently short time scales, the effect will distort this set into a hyperellipsoid, stretched along some directions and contracted along others. The asymptotic rate of expansion of the largest axis, corresponding to the most unstable direction of the flow, is measured by the maximum Lyapunov exponent. [11] Fig 2.1: Describes the average expansions representing Lyapunov exponent[11]. Now if we take the Jacobian matrix of the two dimensional map and calculate it for the nth iterated point and the eigenvalues are calculated which in turn will give the eigenvectors. An eigen value having magnitude less than 1 indicates that the difference of closer points in the direction of corresponding eigen vector is converging i.e. decreasing and on the other hand an eigen value having magnitude greater than 1 indicate that the difference of closer points in the direction of corresponding eigen vector is diverging or increasing. Hence we can say that the eigen values evinces the nature of the magnetism of the nth iterative point. So the logarithm of the magnitude of the eigen values may be termed as the Lyapunov exponents. Let us consider a two dimensional difference equation, which is depending on a parameter a(say) with initial point (x 0,y 0 ).Let the evolved iterated points be (x 1,y 1 ),(x 2,y 2 ).(x n,y n ) at the parameter a=a 0. Suppose that the trajectory moves to the fixed point (x*,y*).since the fixed point is stable and doesn t move under iteration so after finite iterations the trajectory comes very close to the fixed point. Let the Jacobian matrix J * at the fixed point (x *, y * ) of the two dimensional map is JGRMA 2012, All Rights Reserved 51
3 Let the eigenvalue of the J * at the fixed point (x *,y * ) be.now for the nth iterated point (x n,y n ).The Jacobian matrix i.e. J n =J n-1.j n-2..j 0,as n is sufficiently large the Jacobian matrix will be of the form J=J*.J*.J* J n-1.j n-2..j 0, where J* is the Jacobian matrix at the fixed point (x*,y*). So for large iterative value, J* will govern and we may neglect the finite number of J n. Therefore for the large n, the eigen values of J are and then Same argument is valid for periodic points. Thus so long as fixed point or periodic point exists the eigen value will be less than 1 and hence Lyapunov exponent will be negative. At the bifurcation point one of the eigen values is -1 and hence maximum of the Lyapunov exponent will be 0. mle Fig2.2: Showing maximum Lyapunov exponent at the 1 st bifurcation point in y-axes and control parameter in x-axes. However if one of the eigenvalue decreases from positive value to negative and ultimately becomes -1 for the next bifurcation point and if the value is continuously decreasing then it must achieve 0 then the Lyapunov exponent will show negative infinity. However for the parameter value when the maximum Lyapunov exponent becomes positive chaos has started. 3. CALCULATION OF LYAPUNOV EXPONENT OF THE MAP: Here we take our two dimensional discrete model as F(x, y) = F ( ) where f(x, y) = JGRMA 2012, All Rights Reserved 52
4 and g(x, y) = where is the control parameter. The Lyapunov exponent for our two dimensional map is calculated by the eigen values of the limit of the following expression: Where J K is the jacobian given by Now to calculate the lyapunov exponent we first take an initial point and iterate it up to times so that we reach sufficiently nearer to the fixed points. Then we find (J 0.J 1...J n ), where n=100 say, and calculate the eigenvalues of that resultant matrix. Then is the Lyapunov exponent of our model. Fig3.1: Diagram of maximum Lyapunov exponent. Abscissa represents the control parameter while the ordinate represents maximum Lyapunov exponent. In the figure we plot the graph of the maximum lyapunov exponent against parameter values ( ) varies from 0.65 to 1.6. We have considered the initial point as (0.71, 0.71) for the estimate of lyapunov exponent with 100 iterations. Now what is the importance of the graph of the lyapunov exponent? The main importance of the graph is that one can easily identify chaotic region from the periodic one. We see, there are some points of the graph that hit the horizontal line and go to the negative side. These are perioddoubling bifurcation points where lyapunov exponents are almost zero. The points reside on negative side are periodic and the region which lies on the positive side of the parameter is chaotic. Interestingly, the area dominated by chaos occurs beyond the accumulation point Lyapunov exponent near bifurcation point and accumulation point: The following table shows maximum Lyapunov exponent at some parameter values. JGRMA 2012, All Rights Reserved 53
5 Table 2.2: Control parameter maximum Lyapunov exponent Fig 3.2 : Showing part of attractor of the model when initial values are taken as x=1.21 &y=1.21 JGRMA 2012, All Rights Reserved 54
6 Fig 3.3: Showing part of attractor of the model when initial values are taken as x=1.21 &y= GENERALISED CORRELATION DIMENSION [2, 4, 5, 6, 7, 8, 9, 12]: The term "chaotic" and "strange attractor" both has been used to describe the nonperiodic, random like motions in dynamical system. Whereas "chaotic" is meant to convey a loss of information or loss of predictability, the term "strange" is meant to describe the unfamiliar geometric structure on which the motion moves in phase space [9]. In this paper, we will describe a quantitative measure of the strangeness of the attractor. This measure is called the fractal dimension [9]. Since strange attractors are typically characterized by fractal dimensionality, for this purpose, we study different kind of dimensions for our strange attractor [4] Here, we try to determine the dimension of the attractor at the onset of chaos, i.e. at the accumulation point (i.e ) [3] of our two dimensional neural model and.so to study information and box counting dimension,we take the help of generalized dimension, i.e. Renyi s dimension which is defined as follows D q = N ( R) log Pi lim 1 i 1 R 0 q 1 log R q (4.1) where we divide the phase space to N ( R ) boxes of edge R. The probability to find a point in the ith box is denoted by P i and is defined as P i = N i / N The generalization, q, can be any real number; we usually use just an integer q. When we increase q, we give more weight to more populated boxes, and when we decrease q the less occupied boxes become dominant. If q = 0, it gives the box counting dimension For q = 1, by applying the L Hospital rule, equation (4.1) becomes JGRMA 2012, All Rights Reserved 55
7 dimension is where is referred to also as Shanon entropy of the system. The definition of generalised information Note that q = 2 gives the direct arithmetic average that defines the correlation dimension.for a finite set of points G q (N,R) can be approximate as Where H is the Heaviside function and means If, then and If, then. Also here we use the Euclidean norm i.e. if and then Moreover, and Thus, we can see that if we can find G q (N, R), we may get the value of G q (R) and hence we can obtain D q for a particular value of q.for our calculation purpose we calculate G q (30000,r) for different values of q and consider that as G q (r).however it will be justified due to the convergence nature of G q (N,R) as N becomes large. We consider the range of R as 10-7 to 10-1 and get the value of log G q ( R).Then we plot (log R, log G q ( R)) out of which scaling region is selected. In the scaling region we fit a straight line whose slope will give the dimension. 4.1 Box counting dimension [2, 4]: For q = 0 it gives the box counting dimension. We have done the calculation at the parameter with iterations. Calculated value of (log R, log G q (R)) in the scaling region are as follows: JGRMA 2012, All Rights Reserved 56
8 loggq R logr Fig.4.1a: log R vs. log G q (R) in the scaling region JGRMA 2012, All Rights Reserved 57
9 The slope of the above points when fitted with a straight line by least square method is with a mean deviation of The data is obtained from iterated points at the parameter : Information dimension [2, 4] If we take q very near to one i.e. q tends to 1 then that gives the information dimension. At q= at the parameter with iterations. Calculated value of (log R, log G q (R)) in the scaling region are as follows: JGRMA 2012, All Rights Reserved 58
10 loggq R logr Fig.4.2a: log R vs. log G q (R) in the scaling region The slope of the above points when fitted with a straight line by least square method is with a mean deviation of The data is obtained from iterated points at the parameter : Correlation dimension: [2, 4] G q (R) =G q (30000, R) is calculated. The parts of the plotted points (log R, log G q (R)) which follows equation (4.1) is taken. The slope of the fitted straight line in that scaling region is D q. JGRMA 2012, All Rights Reserved 59
11 Calculated value of (log R, log G q (R)) in the scaling region are as follows: JGRMA 2012, All Rights Reserved 60
12 loggq R T.K. Dutta et.al., Journal of Global Research in Mathematical Archives, July-2013, logr Fig4.3.a: log R vs. log G q (R) in the scaling region The slope of the above points when fitted with a straight line by least square method is with a mean deviation of The data is obtained from iterated points at the parameter Conclusion: In this paper we study the Lyapunov exponent of the map and in the figure 3.1, it is seen that at flip bifurcation points [3] maximum Lyapunov exponent are zero. It is also seen that negative Lyapunov exponent signatures stability whereas positive Lyapunov exponent indicates chaos. We also study box counting dimension ( ), Information dimension ( and Correlation dimension ( ) and we see that Correlation dimension D c = Information dimension D I = Box counting dimension D b = Clearly Moreover dimension of the attractors as shown in the figures 3.2 and 3.2 are interestingly found to be less than one. REFERENCES: [1]Crownover,R.M., Introduction to Fractals and Chaos, Jones and Bartlett,Sudbury,MA,1995 [2 ] Dutta, T.K, Jain.A.K and Bhattacharjee.D., Determination of various fractal dimensions in one dimensional models IJST,Vol-2,Issue-3,2013,ISSN JGRMA 2012, All Rights Reserved 61
13 [3]Dutta, T.K, Jain.A.K and Bhattacharjee.D., Period Doubling Scenario Exhibit on Artificial Neural Network Model JGRMA,Vol-1,Issue-6,2013,ISSN [4] Dutta, T.K, and Bhattacharjee.D, Bifurcation, Lyapunov Exponent and fractal Dimensions in a Non-linear Map [2010 AMS Classification: 37G15, 37G35, 37C45]. [5]Falconer, K., Fractal Geometry : Mathematical Foundations and Applications, Wiley, New York,2003. [6]Grassberger.P., Generalised Dimension of the Strange Attractors, Physics Letters,Vol-97A,No-6,1983. [7] Grassberger, P.,Procaccia,I., Characterization of strange Attractors, Physical Review Letters,Vol-50,No-5,1983. [8] Hilborn, R. C., Chaos and Nonlinear Dynamics: An Introduction for Scientists and Engineers., Oxford University Press, 1994 [9] Moon.F.C., Chaos and fractal dynamics,1992 [10] Ott.E., Strange attractors and chaotic motions of dynamical systems, Rev.Mod.Phys.,Vol- 53,No.4,Part 1,Oct [11] Sandri M., Numerical Calculation of Lyapunov Exponent The Mathematica Journal,Miller freeman Publications,1996, Vol 6,Issue 3. [12] Theiler.J., Estimating Fractal Dimension, J.Opt.Soc.Am.A pp ,vol-7,no.6, June JGRMA 2012, All Rights Reserved 62
Handling of Chaos in Two Dimensional Discrete Maps
Handling of Chaos in Two Dimensional Discrete Maps Anil Kumar Jain Assistant Professor, Department of Mathematics Barama College, Barama, Assam, Pincode-781346, India (Email: jainanil965@gmail.com) Abstract:
More informationANALYTICAL PROCESS FOR DETERMINATION OF HOPF BIFURCATIONS
Volume 1, No. 8, August 2013 Journal of Global Research in Mathematical Archives RESEARCH PAPER Available online at http://www.jgrma.info ANALYTICAL PROCESS FOR DETERMINATION OF HOPF BIFURCATIONS Tarini
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 informationExample Chaotic Maps (that you can analyze)
Example Chaotic Maps (that you can analyze) Reading for this lecture: NDAC, Sections.5-.7. Lecture 7: Natural Computation & Self-Organization, Physics 256A (Winter 24); Jim Crutchfield Monday, January
More informationCHAPTER 2 FEIGENBAUM UNIVERSALITY IN 1-DIMENSIONAL NONLINEAR ALGEBRAIC MAPS
CHAPTER 2 FEIGENBAUM UNIVERSALITY IN 1-DIMENSIONAL NONLINEAR ALGEBRAIC MAPS The chief aim of this chapter is to discuss the dynamical behaviour of some 1-dimensional discrete maps. In this chapter, we
More informationThe Hopf Bifurcation Theorem: Abstract. Introduction. Transversality condition; the eigenvalues cross the imginary axis with non-zero speed
Supercritical and Subcritical Hopf Bifurcations in Non Linear Maps Tarini Kumar Dutta, Department of Mathematics, Gauhati Universit Pramila Kumari Prajapati, Department of Mathematics, Gauhati Universit
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 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 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 information6.2 Their Derivatives
Exponential Functions and 6.2 Their Derivatives Copyright Cengage Learning. All rights reserved. Exponential Functions and Their Derivatives The function f(x) = 2 x is called an exponential function because
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 informationChapter 3. Gumowski-Mira Map. 3.1 Introduction
Chapter 3 Gumowski-Mira Map 3.1 Introduction Non linear recurrence relations model many real world systems and help in analysing their possible asymptotic behaviour as the parameters are varied [17]. Here
More informationCo-existence of Regular and Chaotic Motions in the Gaussian Map
EJTP 3, No. 13 (2006) 29 40 Electronic Journal of Theoretical Physics Co-existence of Regular and Chaotic Motions in the Gaussian Map Vinod Patidar Department of Physics, Banasthali Vidyapith Deemed University,
More informationSome Dynamical Behaviors In Lorenz Model
International Journal Of Computational Engineering Research (ijceronline.com) Vol. Issue. 7 Some Dynamical Behaviors In Lorenz Model Dr. Nabajyoti Das Assistant Professor, Department of Mathematics, Jawaharlal
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 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 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 informationONE DIMENSIONAL CHAOTIC DYNAMICAL SYSTEMS
Journal of Pure and Applied Mathematics: Advances and Applications Volume 0 Number 0 Pages 69-0 ONE DIMENSIONAL CHAOTIC DYNAMICAL SYSTEMS HENA RANI BISWAS Department of Mathematics University of Barisal
More information3 Polynomial and Rational Functions
3 Polynomial and Rational Functions 3.1 Polynomial Functions and their Graphs So far, we have learned how to graph polynomials of degree 0, 1, and. Degree 0 polynomial functions are things like f(x) =,
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 informationGenerating a Complex Form of Chaotic Pan System and its Behavior
Appl. Math. Inf. Sci. 9, No. 5, 2553-2557 (2015) 2553 Applied Mathematics & Information Sciences An International Journal http://dx.doi.org/10.12785/amis/090540 Generating a Complex Form of Chaotic Pan
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 informationUnit Ten Summary Introduction to Dynamical Systems and Chaos
Unit Ten Summary Introduction to Dynamical Systems Dynamical Systems A dynamical system is a system that evolves in time according to a well-defined, unchanging rule. The study of dynamical systems is
More informationIs the Hénon map chaotic
Is the Hénon map chaotic Zbigniew Galias Department of Electrical Engineering AGH University of Science and Technology, Poland, galias@agh.edu.pl International Workshop on Complex Networks and Applications
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 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: January 08, 2018, at 08 30 12 30 Johanneberg Kristian
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 informationA Novel Three Dimension Autonomous Chaotic System with a Quadratic Exponential Nonlinear Term
ETASR - Engineering, Technology & Applied Science Research Vol., o.,, 9-5 9 A Novel Three Dimension Autonomous Chaotic System with a Quadratic Exponential Nonlinear Term Fei Yu College of Information Science
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 informationCharacterizing Dynamics of a Physical System
Characterizing Dynamics of a Physical System How the dynamics of a physical system is investigated Abinash Chakraborty National Institute of Technology Rourkela, India Abstract In this paper, we shall
More informationA Two-dimensional Discrete Mapping with C Multifold Chaotic Attractors
EJTP 5, No. 17 (2008) 111 124 Electronic Journal of Theoretical Physics A Two-dimensional Discrete Mapping with C Multifold Chaotic Attractors Zeraoulia Elhadj a, J. C. Sprott b a Department of Mathematics,
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: January 14, 2019, at 08 30 12 30 Johanneberg Kristian
More informationBifurcations in the Quadratic Map
Chapter 14 Bifurcations in the Quadratic Map We will approach the study of the universal period doubling route to chaos by first investigating the details of the quadratic map. This investigation suggests
More informationChaos in adaptive expectation Cournot Puu duopoly model
Chaos in adaptive expectation Cournot Puu duopoly model Varun Pandit 1, Dr. Brahmadeo 2, and Dr. Praveen Kulshreshtha 3 1 Economics, HSS, IIT Kanpur(India) 2 Ex-Professor, Materials Science and 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 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 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 informationChaos in GDP. Abstract
Chaos in GDP R. Kříž Abstract This paper presents an analysis of GDP and finds chaos in GDP. I tried to find a nonlinear lower-dimensional discrete dynamic macroeconomic model that would characterize GDP.
More informationPersonal notes on renormalization
Personal notes on renormalization Pau Rabassa Sans April 17, 2009 1 Dynamics of the logistic map This will be a fast (and selective) review of the dynamics of the logistic map. Let us consider the logistic
More informationLesson 4: Non-fading Memory Nonlinearities
Lesson 4: Non-fading Memory Nonlinearities Nonlinear Signal Processing SS 2017 Christian Knoll Signal Processing and Speech Communication Laboratory Graz University of Technology June 22, 2017 NLSP SS
More informationAre numerical studies of long term dynamics conclusive: the case of the Hénon map
Journal of Physics: Conference Series PAPER OPEN ACCESS Are numerical studies of long term dynamics conclusive: the case of the Hénon map To cite this article: Zbigniew Galias 2016 J. Phys.: Conf. Ser.
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 information16 Period doubling route to chaos
16 Period doubling route to chaos We now study the routes or scenarios towards chaos. We ask: How does the transition from periodic to strange attractor occur? The question is analogous to the study of
More informationOne 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 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 informationDynamical 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 informationThe Pattern Recognition System Using the Fractal Dimension of Chaos Theory
Original Article International Journal of Fuzzy Logic and Intelligent Systems Vol. 15, No. 2, June 2015, pp. 121-125 http://dx.doi.org/10.5391/ijfis.2015.15.2.121 ISSN(Print) 1598-2645 ISSN(Online) 2093-744X
More informationChaos in multiplicative systems
Chaos in multiplicative systems Dorota Aniszewska 1 and Marek Rybaczuk 2 1 Institute of Materials Science and Applied Mechanics Wroclaw University of Technology 50-370 Wroclaw, Smoluchowskiego 25 (e-mail:
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 informationFORECASTING ECONOMIC GROWTH USING CHAOS THEORY
Article history: Received 22 April 2016; last revision 30 June 2016; accepted 12 September 2016 FORECASTING ECONOMIC GROWTH USING CHAOS THEORY Mihaela Simionescu Institute for Economic Forecasting of the
More informationA Lagrangian approach to the kinematic dynamo
1 A Lagrangian approach to the kinematic dynamo Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc/ 5 March 2001 with Allen
More informationTWO DIMENSIONAL FLOWS. Lecture 5: Limit Cycles and Bifurcations
TWO DIMENSIONAL FLOWS Lecture 5: Limit Cycles and Bifurcations 5. Limit cycles A limit cycle is an isolated closed trajectory [ isolated means that neighbouring trajectories are not closed] Fig. 5.1.1
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 informationAP * Calculus Review. Limits, Continuity, and the Definition of the Derivative
AP * Calculus Review Limits, Continuity, and the Definition of the Derivative Teacher Packet Advanced Placement and AP are registered trademark of the College Entrance Examination Board. The College Board
More informationNonlinear dynamics in the Cournot duopoly game with heterogeneous players
arxiv:nlin/0210035v1 [nlin.cd] 16 Oct 2002 Nonlinear dynamics in the Cournot duopoly game with heterogeneous players H. N. Agiza and A. A. Elsadany Department of Mathematics, Faculty of Science Mansoura
More informationChapter 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 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 informationRené Thomas Université de Bruxelles. Frontier diagrams: a global view of the structure of phase space.
René Thomas Université de Bruxelles Frontier diagrams: a global view of the structure of phase space. We have the tools to identify and characterise steady states and trajectories. But WHY several steady
More informationTHE theory of dynamical systems is an extremely
Attractors: Nonstrange to Chaotic Robert L. V. Taylor The College of Wooster rtaylor@wooster.edu Advised by Dr. John David The College of Wooster jdavid@wooster.edu Abstract The theory of chaotic dynamical
More informationChapter 2 Chaos theory and its relationship to complexity
Chapter 2 Chaos theory and its relationship to complexity David Kernick This chapter introduces chaos theory and the concept of non-linearity. It highlights the importance of reiteration and the system
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 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 informationChapitre 4. Transition to chaos. 4.1 One-dimensional maps
Chapitre 4 Transition to 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 manners
More informationComparative Analysis of Convergence of Various Numerical Methods
Journal of Computer and Mathematical Sciences, Vol.6(6),290-297, June 2015 (An International Research Journal), www.compmath-journal.org ISSN 0976-5727 (Print) ISSN 2319-8133 (Online) Comparative Analysis
More informationEdward Lorenz. Professor of Meteorology at the Massachusetts Institute of Technology
The Lorenz system Edward Lorenz Professor of Meteorology at the Massachusetts Institute of Technology In 1963 derived a three dimensional system in efforts to model long range predictions for the weather
More informationExample of a Blue Sky Catastrophe
PUB:[SXG.TEMP]TRANS2913EL.PS 16-OCT-2001 11:08:53.21 SXG Page: 99 (1) Amer. Math. Soc. Transl. (2) Vol. 200, 2000 Example of a Blue Sky Catastrophe Nikolaĭ Gavrilov and Andrey Shilnikov To the memory of
More 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 informationChaos Theory. Namit Anand Y Integrated M.Sc.( ) Under the guidance of. Prof. S.C. Phatak. Center for Excellence in Basic Sciences
Chaos Theory Namit Anand Y1111033 Integrated M.Sc.(2011-2016) Under the guidance of Prof. S.C. Phatak Center for Excellence in Basic Sciences University of Mumbai 1 Contents 1 Abstract 3 1.1 Basic Definitions
More informationControl and synchronization of Julia sets of the complex dissipative standard system
Nonlinear Analysis: Modelling and Control, Vol. 21, No. 4, 465 476 ISSN 1392-5113 http://dx.doi.org/10.15388/na.2016.4.3 Control and synchronization of Julia sets of the complex dissipative standard system
More informationarxiv: v1 [nlin.cd] 22 Feb 2011
Generalising the logistic map through the q-product arxiv:1102.4609v1 [nlin.cd] 22 Feb 2011 R W S Pessoa, E P Borges Escola Politécnica, Universidade Federal da Bahia, Rua Aristides Novis 2, Salvador,
More information2 Discrete growth models, logistic map (Murray, Chapter 2)
2 Discrete growth models, logistic map (Murray, Chapter 2) As argued in Lecture 1 the population of non-overlapping generations can be modelled as a discrete dynamical system. This is an example of an
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 informationDeterministic Chaos Lab
Deterministic Chaos Lab John Widloski, Robert Hovden, Philip Mathew School of Physics, Georgia Institute of Technology, Atlanta, GA 30332 I. DETERMINISTIC CHAOS LAB This laboratory consists of three major
More informationFundamentals of Dynamical Systems / Discrete-Time Models. Dr. Dylan McNamara people.uncw.edu/ mcnamarad
Fundamentals of Dynamical Systems / Discrete-Time Models Dr. Dylan McNamara people.uncw.edu/ mcnamarad Dynamical systems theory Considers how systems autonomously change along time Ranges from Newtonian
More informationHamiltonian Chaos and the standard map
Hamiltonian Chaos and the standard map Outline: What happens for small perturbation? Questions of long time stability? Poincare section and twist maps. Area preserving mappings. Standard map as time sections
More informationChapter 6 Nonlinear Systems and Phenomena. Friday, November 2, 12
Chapter 6 Nonlinear Systems and Phenomena 6.1 Stability and the Phase Plane We now move to nonlinear systems Begin with the first-order system for x(t) d dt x = f(x,t), x(0) = x 0 In particular, consider
More informationSolutions 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 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 informationCollege Algebra. Chapter 5 Review Created by: Lauren Atkinson. Math Coordinator, Mary Stangler Center for Academic Success
College Algebra Chapter 5 Review Created by: Lauren Atkinson Math Coordinator, Mary Stangler Center for Academic Success Note: This review is composed of questions from the chapter review at the end of
More informationChaotifying 2-D piecewise linear maps via a piecewise linear controller function
Chaotifying 2-D piecewise linear maps via a piecewise linear controller function Zeraoulia Elhadj 1,J.C.Sprott 2 1 Department of Mathematics, University of Tébéssa, (12000), Algeria. E-mail: zeraoulia@mail.univ-tebessa.dz
More informationCHAOS -SOME BASIC CONCEPTS
CHAOS -SOME BASIC CONCEPTS Anders Ekberg INTRODUCTION This report is my exam of the "Chaos-part" of the course STOCHASTIC VIBRATIONS. I m by no means any expert in the area and may well have misunderstood
More informationLecture 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 informationIntroduction 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 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 informationCharacterization of the natural measure by unstable periodic orbits in nonhyperbolic chaotic systems
PHYSICAL REVIEW E VOLUME 56, NUMBER 6 DECEMBER 1997 Characterization of the natural measure by unstable periodic orbits in nonhyperbolic chaotic systems Ying-Cheng Lai * Department of Physics and Astronomy
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 informationMATH 215/255 Solutions to Additional Practice Problems April dy dt
. For the nonlinear system MATH 5/55 Solutions to Additional Practice Problems April 08 dx dt = x( x y, dy dt = y(.5 y x, x 0, y 0, (a Show that if x(0 > 0 and y(0 = 0, then the solution (x(t, y(t of the
More informationLIMITS AND DERIVATIVES
2 LIMITS AND DERIVATIVES LIMITS AND DERIVATIVES 2.2 The Limit of a Function In this section, we will learn: About limits in general and about numerical and graphical methods for computing them. THE LIMIT
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 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 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 informationA Novel Hyperchaotic System and Its Control
1371371371371378 Journal of Uncertain Systems Vol.3, No., pp.137-144, 009 Online at: www.jus.org.uk A Novel Hyperchaotic System and Its Control Jiang Xu, Gouliang Cai, Song Zheng School of Mathematics
More informationWORD: EXAMPLE(S): COUNTEREXAMPLE(S): EXAMPLE(S): COUNTEREXAMPLE(S): WORD: EXAMPLE(S): COUNTEREXAMPLE(S): EXAMPLE(S): COUNTEREXAMPLE(S): WORD:
Bivariate Data DEFINITION: In statistics, data sets using two variables. Scatter Plot DEFINITION: a bivariate graph with points plotted to show a possible relationship between the two sets of data. Positive
More informationn-dimensional LCE code
n-dimensional LCE code Dale L Peterson Department of Mechanical and Aeronautical Engineering University of California, Davis dlpeterson@ucdavisedu June 10, 2007 Abstract The Lyapunov characteristic exponents
More informationPHY411 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 informationChaos in the Dynamics of the Family of Mappings f c (x) = x 2 x + c
IOSR Journal of Mathematics (IOSR-JM) e-issn: 78-578, p-issn: 319-765X. Volume 10, Issue 4 Ver. IV (Jul-Aug. 014), PP 108-116 Chaos in the Dynamics of the Family of Mappings f c (x) = x x + c Mr. Kulkarni
More informationThe phenomenon: complex motion, unusual geometry
Part I The phenomenon: complex motion, unusual geometry Chapter 1 Chaotic motion 1.1 What is chaos? Certain long-lasting, sustained motion repeats itself exactly, periodically. Examples from everyday life
More informationOn the Dynamics of a n-d Piecewise Linear Map
EJTP 4, No. 14 2007 1 8 Electronic Journal of Theoretical Physics On the Dynamics of a n-d Piecewise Linear Map Zeraoulia Elhadj Department of Mathematics, University of Tébéssa, 12000, Algeria. Received
More informationPeriod-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