Complicated behavior of dynamical systems. Mathematical methods and computer experiments.
|
|
- Ross Blankenship
- 5 years ago
- Views:
Transcription
1 Complicated behavior of dynamical systems. Mathematical methods and computer experiments. Kuznetsov N.V. 1, Leonov G.A. 1, and Seledzhi S.M. 1 St.Petersburg State University Universitetsky pr St.Petersburg, Russia ( leonov@math.spbu.ru) Abstract. It is shown that the most suitable definition of chaos in dynamical systems is instability by Zhukovsky. There are considered computer experiments and analytical proofs are corresponding to this definition. Also it is shown that the Filter Hypothesis is untrue. Keywords: Zhukovsky stability, Harmonic linearization. 1 Introduction Complicated behavior in dynamical systems are attracting attentions of various researchers (Moon[1], Shuster[2], Kolesnikov[3], Magnitskii and Sidorov[4], Leonov[5]). In this work we consider some approaches for investigation such complicated behaviors of trajectories as strange attractors, complicated periodic solutions, flattering of trajectories. Here it is shown that from among the classical notions of instability for studying strange attractors in the continuous case, the most adequate ones are Zhukovsky instability. In order to investigate Zhukovsky stability, a new research tool a moving Poincaré section is introduced (Leonov[5]). The complicated behavior of periodic solutions and filter hypothesis are considered. Difficulties of numerical investigations of limit cycles of two dimensional autonomous systems are considered. 2 Zhukovsky stability It is well known (Leonov[5]) that the defining property of strange attractors is the sensitivity of their trajectories with respect to the initial data. The notion of Zhukovsky instability is adequate to the sensitivity of trajectories with respect to the initial data for continuous dynamical systems. In order to investigate Zhukovsky stability, a new research tool a moving Poincaré section is introduced. With the help of this tool, extensions of the widelyknown theorems of Andronov Witt, Demidovich, Borg, and Poincaré are carried out.
2 2 Kuznetsov, Leonov and Seledzhi Consider the dynamical systems, generated by the differential equations dx dt = f(x), t R1, x R n. (1) We introduce now the definition of Zhukovsky stability. For this purpose we need to consider the following set of homeomorphisms Hom = {τ( ) τ : [0, + ) [0, + ), τ(0) = 0}. The functions τ(t) from the set Hom play the role of the reparametrization of time for the trajectories of system (1). Definition. The trajectory x(t, x 0 ) of system (1) is said to be Zhukovsky stable if for any number ε > 0 there exists the number δ(ε) > 0 such that for any vector y 0, satisfying the inequality x 0 y 0 δ(ε), the function τ( ) Hom can be found such that the following inequality x(t, x 0 ) x(τ(t), y 0 ) ε, t 0 is valid. If, in addition, for the certain number δ 0 > 0 and any y 0 from the ball {y x 0 y δ 0 } the function τ( ) Hom can be found such that the relation holds lim t + x(t, x 0) x(τ(t), y 0 ) = 0, then we shall say that the trajectory x(t, x 0 ) is asymptotically stable by Zhukovsky. This means that the stability by Zhukovsky is a stability by Lyapunov for the suitable reparametrization of each of perturbed trajectories. Consider here relationships between the basic stability notions. There are well known examples of periodic trajectories of continuous systems that happen to be Lyapunov unstable but Poincaré stable. Example. Consider the Duffing equation ẍ + x + x 3 = 0 It is easily seen (Leonov[10]) that all solutions of are periodic and the period of the solutions is depend on initial date x 0. Here the trajectories of the equation remain wholly within small neighborhoods of each other for small changes in initial data; hence they are Poincaré stable. Since all trajectories considered are periodic, they are also Zhukovsky stable. However, in this case Lyapunov stability is lacking. Example. Consider the linearized equations of two decoupled pendula: ẋ 1 = y 1, ẏ 1 = ω1x 2 1, ẋ 2 = y 2, ẏ 2 = ω2 2x 2. (2) with ω 1 /ω 2 irrational. Change the flow of trajectories on the tori as follows. Cut the toroidal surface along a certain segment of the fixed trajectory from
3 Complicated behavior of dynamical systems 3 Fig.1. Heteroclinic trajectories with hole H. the point z 1 to the point z 2 (Fig. 1). Then the surface is stretched diffeomorphically along the torus so that a cut is mapped into the circle with the fixed points z 1 and z 2. Denote by H the interior of the circle. Outside the hole H, after the diffeomorphic stretching, the disposition of trajectories on the torus is the same. Outside the hole H, the trajectories are everywhere dense on torus. They are therefore, as before, asymptotically Poincaré stable. But it is possible here to choose two trajectories in a small neighborhood of z 1 that these trajectories are repelled from each other; hence, the trajectory leaving z 0 is Zhukovsky unstable. Thus, a trajectory can be asymptotically Poincaré stable and Zhukovsky unstable. Hence, from among the classical notions of instability for studying strange attractors in the continuous case, the most adequate ones are Zhukovsky instability. The reparametrization of trajectories permits us to introduce interesting and important tool of investigations, namely the moving Poincare section. The classical Poincare section is the transversal (n 1)-dimensional surface S in the phase space R n, which possesses a recurring property. We force now the Poincare section to move along the trajectory x(t, x 0 ) Fig. 2.
4 4 Kuznetsov, Leonov and Seledzhi We assume further that the vector-function f(x) is twice continuously differentiable and the trajectory x(t, x 0 ), whose the Zhukovsky stability (or instability) will be considered, is wholly situated in the certain bounded domain Ω R n for t 0. Suppose also that f(x) 0, x Ω. Here Ω is a closure of the domain Ω. Under these assumptions there exist the positive numbers δ and ε such that the following relation f(y) f(x) δ, y S(x, ε), x Ω is satisfied. Here S(x, ε) = {y (y x) f(x) = 0, x y < ε}. Definition. The set S(x(t, x 0 ), ε) is called a moving Poincare section. The classical Poincare section allow us to clear up the behavior of trajectories using the information at discrete time of their crossings with the Poincare section. The reparametrization of trajectories makes it possible to organize the motion of trajectories so that at time t all trajectories are situated on the same moving Poincare section S(x(t, x 0 ), ε): x(ϕ(t), y 0 ) S(x(t, x 0 ), ε). (3) Here ϕ(t) is a reparametrization of the trajectory x(t, y 0 ), y 0 S(x 0, ε). This consideration has, of course, a local property and it is only possible for the values t such that x(ϕ(t), y 0 ) x(t, x 0 ) < ε. Conside system (1) in the form where dz dt = A(x(t, x 0))z + g(t, z), z f(x(t, x 0 )) = 0, (4) A(x) = f [ ( ) ] f(x)f(x) f f (x) x f(x) 2 x (x) + x (x). Here the following relations g(t, z) f(x(t, x 0 )) = 0, g(t, z) = O( z 2 ) are valid. Thus, for system (4) we have the system of the first approximation dv dt = A(x(t, x 0))v, f(x(t, x 0 )) v = 0. (5) It differs from the usual system of the first approximation dw dt = f x (x(t, x 0))w (6) in that we introduce here the projector v = (I f(x(t, x 0)f(x(t, x 0 )) ) f(x(t, x 0 ) 2 w. (7) Theorem. Let Λ is the higher characteristic exponent of system (5), Γ is a coefficient of irregularity. Then if for system (5) the inequality Λ + Γ < 0,
5 Complicated behavior of dynamical systems 5 is satisfied, then the trajectory x(t, x 0 ) is asymptotically Zhukovsky stable. It possible to show that system (6) has the one null characteristic exponent λ 1. Then denote by λ 2 λ n the rest of characteristic exponents of system (6). From relations (7) it follows that the characteristic exponents of system (5) are not greater than the corresponding characteristic exponents of system (6). From this we have the inequality γ Γ. (8) Here γ is the coefficient of irregularity of system (6). Besides we have λ 2 Λ. (9) Theorem 4 and inequalities (8) and (9) give the following Theorem 5. If for system (6) the following inequality λ 2 + γ < 0 (10) is satisfied, then the trajectory x(t, x 0 ) is asymptotically Zhukovsky stable. This result is the generalization of the well-known Andronov Witt theorem. Theorem (Andronov and Witt[6]) 6. If the trajectory x(t, x 0 ) is periodic, differs from equilibria and for system (6) the inequality λ 2 < 0 is satisfied, then the trajectory x(t, x 0 ) is asymptotically orbitally stable (asymptotically Poincare stable). Theorem 6 is a corollary of Theorem 5 since system (6) with the periodic matrix f x (x(t, x 0)) is regular (Demidovich[7]). Recall that for the periodic trajectories the asymptotic stability by Zhukovsky and by Poincare are equivalent. The theorem of Demidovich is also a corollary of Theorem 5. Theorem 7 (Demidovich[8]). If system (6) is regular (i.e. γ = 0) and the inequality λ 2 < 0 is satisfied, then the trajectory x(t, x 0 ) is asymptotically orbitally stable. 3 The harmonic linearization method describing function method Consider a system dx dt = Px + qϕ(r x), (11) where P is a constant n n-matrix, q and r are constant n-vectors, ϕ(σ) is a piecewise-continuous function, and is the operation of transposition. When applied the harmonic linearization method (Khalil[9]) to this system, standard assumptions are the existence of a pair of purely imaginary eigen
6 6 Kuznetsov, Leonov and Seledzhi values ±iω 0 (ω 0 > 0) of the matrix P and a negativeness of the rest of eigen values. By nonsingular linear transformation, under the above assumptions system (1) can be reduced to the form ẋ 1 = ω 0 x 2 + b 1 ϕ(x 1 + c x 3 ) ẋ 2 = ω 0 x 1 + b 2 ϕ(x 1 + c x 3 ) ẋ 3 = Ax 3 + bϕ(x 1 + c x 3 ). (12) Here A is a constant (n 2) (n 2)-matrix, all eigen values of which have negative real parts, b and c are (n 2)-dimensional vectors, b 1 and b 2 are certain numbers. Combined application of the harmonic linearization method, the classical method of small parameter, and numerical methods permit us to compute periodic oscillations of certain multistage procedure, where at the first step the harmonic linearization method is applied. In the basic, noncritical, case we suppose that the relation ϕ(σ) = εψ(σ), where ε is a small parameter, is satisfied. In the sequel, without loss of generality, we can assume that for A there exists the number α > 0 such that x 3 (A + A )x 3 2α x 3 2, x 3 R n 2. We introduce the function K(a) = 2π/ω 0 0 ψ(cos(ω 0 t)a)cos(ω 0 t)dt. Theorem. If the conditions K(a) = 0, da < 0, are satisfied, then for sufficiently small ε > 0 system (11) has T-periodic solution such that b 1 dk(a) r x(t) = a cos(ω 0 t) + O(ε), T = 2π ω 0 + O(ε). This periodic solution is stable in the sense that there exists its certain ε- neighbourhood such that all solutions with the initial data from this ε-neighbourhood remain in it in increasing time t. The described here standard, basic, method of harmonic linearization turns out too rough to locate periodic oscillations in nonlinear systems, satisfying the generalized Routh-Hurwitz conditions. The extension of this Theorem in the spirit of classical research of critical cases in the theory of motion stability makes it possible to obtain effective estimates for periodic oscillations in systems, satisfying the generalized Routh-Hurwitz conditions. The described periodic solution can be considered as certain support (basic) periodic oscillations and system (11) with the considered above nonlinearities as generating start system in the algorithms of seeking the periodic solutions of another system, namely dx dt = P 0x + qf(r x), (13)
7 Complicated behavior of dynamical systems 7 In this case, we can organize a finite sequence of the functions ϕ j (σ) j = 1,...,m, such that the graphs of each pair ϕ j and ϕ j+1 are close to each other and ϕ 1 (σ) = ϕ(σ), ϕ m (σ) = f(σ). Then, for the system dx dt = Px + qϕ j(r x) (14) with ϕ 1 (σ) = ϕ(σ) and small ε we take the periodic solution g 1 (t), described in either Theorem 1 or Theorem 2. Two cases occur: either all points of this periodic solution are situated in a domain of attraction of the stable periodic solution g 2 (t) of system (14) with j = 2, or in passing from system (14) with j = 1 to system (14) with j = 2 we have a bifurcation of stability loss and a vanishing of periodic solution. In the first case, we can numerically find g 2 (t) when the trajectory of system (14) with j = 2 begins at the initial point x(0) = g 1 (0). Starting from the point g 1 (0), after transient process the computational procedure outputs into the periodic solution g 2 (t) and calculates it. For this purpose the interval [0, τ], on which the computation occurs, must be sufficiently large. After computation of g 2 (t) it is possible to go to the following system (14) with j = 3 and to organize a similar procedure for computation of the periodic solution g 3 (t) when a trajectory, which in increasing t tends to the periodic trajectory g 3 (t), starts from the initial point x(0) = g 2 (τ). Proceeding then this procedure for sequential computing g j (t) and making use of trajectories of system (14) with the initial data x(0) = g j 1 (τ), we arrive by numerical computation of periodic solution of system (13) or observe, at a certain step, a bifurcation of stability loss and a vanishing of periodic solution. We give two examples of applying this procedure. Fig. 3. Periodic solution of system for ε = 1 Example (Leonov[10], Leonov[11]). Consider system (12) with the function ϕ(σ) = εψ(σ), ψ(σ) = k 1 σ + k 3 σ 3. Then we have K(a) = (k 1 a k 3a 3 ) π ω 0. It follows that a can be determined from the equation K(a) = 0
8 8 Kuznetsov, Leonov and Seledzhi in such way a = a 1 = 4k1 3k 3, and a stability condition takes the form b 1 k 1 < 0. Let be k 1 = 3, k 3 = 4, ω 0 = 1, b 1 = 1, b 2 = 1, A = 1, c = 1, b = 1. Then a 1 = 1. Using the classical harmonic linearization method (Khalil[9]), we obtain that for any ε > 0, system (12) has a periodic solution and in this case σ(t) = r x(t) cost. By Theorem 1 for small ε > 0, system (12) has a periodic solution of the form x 1 (t) = cost + O(ε), x 2 (t) = sint + O(ε), x 3 (t) = O(ε). Further, using the above computational procedure, we obtain a periodic solution of system (12) for ϕ j (σ) = ε j ϕ(σ), ε 1 = 0, 1, ε 2 = 0, 3, ε 3 = 0, 6, ε 4 = 0, 7, ε 5 = 0, 9, ε 6 = 1. In Fig. 3 it is shown a projection of periodic solution thus computed on the plane {x 1, x 2 } when ε = 1. For these periodic solution the graph σ(t) = x 1 (t) + x 3 (t) is represented. In this case for ε = 1, the output σ(t) is substantially not harmonic and the filter hypothesis is untrue. Therefore, here it is impossible, in principle, to justify a standard harmonic linearization method, founding on the filter hypothesis. References 1.F. Moon. Chaotic Vibrations. John Wiley. New York, H. G. Schuster. Deterministic Chaos. Physik Verlag, Weinheim, A. Kolesnikov. Synergetics Methods of Complex Systems Control: System Synthesis Theory. Moscow: URSS, N. A. Magnitskii, S. V. Sidorov. New Methods for Chaotic Dynamics, World Scientific, G. A. Leonov. Strange attractors and classical stability theory. Saint-Petersburg State University Press, A. A. Andronov, A. A. Witt. On Lyapunov Stability. Journal of Experimental and Theoretical Physics. Vol. 3, N 3, B.P. Demidovich. Lectures on Mathematical Theory of Stability. Nauka, Moscow, B. P. Demidovich. Orbital Stability of the Solutions Autonomous Systems. I, II// Differential Equations, 1968, vol. 4, p , N 8, p (Translated from Differenzialnya Uravneniya). 9.H. K. Khalil. Nonlinear Systems. Prentice Hall, G. A. Leonov. Method of Harmonic Linearization. Automation and Remote Control, G. A. Leonov. Aizerman Problem. Automation and Remote Control, [in print]
Mathematical 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 informationNonlinear 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 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 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 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 informationLyapunov Stability Theory
Lyapunov Stability Theory Peter Al Hokayem and Eduardo Gallestey March 16, 2015 1 Introduction In this lecture we consider the stability of equilibrium points of autonomous nonlinear systems, both in continuous
More informationSIMPLE CHAOTIC FLOWS WITH ONE STABLE EQUILIBRIUM
International Journal of Bifurcation and Chaos, Vol. 23, No. 11 (2013) 1350188 (7 pages) c World Scientific Publishing Company DOI: 10.1142/S0218127413501885 SIMPLE CHAOTIC FLOWS WITH ONE STABLE EQUILIBRIUM
More informationLECTURE 8: DYNAMICAL SYSTEMS 7
15-382 COLLECTIVE INTELLIGENCE S18 LECTURE 8: DYNAMICAL SYSTEMS 7 INSTRUCTOR: GIANNI A. DI CARO GEOMETRIES IN THE PHASE SPACE Damped pendulum One cp in the region between two separatrix Separatrix Basin
More informationBIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs
BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs Yuri A. Kuznetsov August, 2010 Contents 1. Solutions and orbits. 2. Equilibria. 3. Periodic orbits and limit cycles. 4. Homoclinic orbits.
More informationPoincaré Map, Floquet Theory, and Stability of Periodic Orbits
Poincaré Map, Floquet Theory, and Stability of Periodic Orbits CDS140A Lecturer: W.S. Koon Fall, 2006 1 Poincaré Maps Definition (Poincaré Map): Consider ẋ = f(x) with periodic solution x(t). Construct
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 informationStability of Nonlinear Systems An Introduction
Stability of Nonlinear Systems An Introduction Michael Baldea Department of Chemical Engineering The University of Texas at Austin April 3, 2012 The Concept of Stability Consider the generic nonlinear
More informationMCE693/793: Analysis and Control of Nonlinear Systems
MCE693/793: Analysis and Control of Nonlinear Systems Systems of Differential Equations Phase Plane Analysis Hanz Richter Mechanical Engineering Department Cleveland State University Systems of Nonlinear
More informationRecent new examples of hidden attractors
Eur. Phys. J. Special Topics 224, 1469 1476 (2015) EDP Sciences, Springer-Verlag 2015 DOI: 10.1140/epjst/e2015-02472-1 THE EUROPEAN PHYSICAL JOURNAL SPECIAL TOPICS Review Recent new examples of hidden
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 informationNonlinear Analysis and Design of Phase-Locked Loops (PLL)
Nonlinear Analysis and Design of Phase-Locked Loops (PLL) Nikolay V. Kuznetsov Gennady A. Leonov Svetlana M. Seledzhi Department of Applied Cybernetics Faculty of Mathematics and Mechanics Saint-Petersburg
More informationThere is a more global concept that is related to this circle of ideas that we discuss somewhat informally. Namely, a region R R n with a (smooth)
82 Introduction Liapunov Functions Besides the Liapunov spectral theorem, there is another basic method of proving stability that is a generalization of the energy method we have seen in the introductory
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 informationIntroduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems
p. 1/5 Introduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems p. 2/5 Time-varying Systems ẋ = f(t, x) f(t, x) is piecewise continuous in t and locally Lipschitz in x for all t
More informationB5.6 Nonlinear Systems
B5.6 Nonlinear Systems 4. Bifurcations Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Local bifurcations for vector fields 1.1 The problem 1.2 The extended centre
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 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 informationEN Nonlinear Control and Planning in Robotics Lecture 3: Stability February 4, 2015
EN530.678 Nonlinear Control and Planning in Robotics Lecture 3: Stability February 4, 2015 Prof: Marin Kobilarov 0.1 Model prerequisites Consider ẋ = f(t, x). We will make the following basic assumptions
More informationHandout 2: Invariant Sets and Stability
Engineering Tripos Part IIB Nonlinear Systems and Control Module 4F2 1 Invariant Sets Handout 2: Invariant Sets and Stability Consider again the autonomous dynamical system ẋ = f(x), x() = x (1) with state
More informationLMI Methods in Optimal and Robust Control
LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 15: Nonlinear Systems and Lyapunov Functions Overview Our next goal is to extend LMI s and optimization to nonlinear
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 informationLecture 4. Chapter 4: Lyapunov Stability. Eugenio Schuster. Mechanical Engineering and Mechanics Lehigh University.
Lecture 4 Chapter 4: Lyapunov Stability Eugenio Schuster schuster@lehigh.edu Mechanical Engineering and Mechanics Lehigh University Lecture 4 p. 1/86 Autonomous Systems Consider the autonomous system ẋ
More information3 Stability and Lyapunov Functions
CDS140a Nonlinear Systems: Local Theory 02/01/2011 3 Stability and Lyapunov Functions 3.1 Lyapunov Stability Denition: An equilibrium point x 0 of (1) is stable if for all ɛ > 0, there exists a δ > 0 such
More informationSTABILITY. Phase portraits and local stability
MAS271 Methods for differential equations Dr. R. Jain STABILITY Phase portraits and local stability We are interested in system of ordinary differential equations of the form ẋ = f(x, y), ẏ = g(x, y),
More informationIntroduction to Applied Nonlinear Dynamical Systems and Chaos
Stephen Wiggins Introduction to Applied Nonlinear Dynamical Systems and Chaos Second Edition With 250 Figures 4jj Springer I Series Preface v L I Preface to the Second Edition vii Introduction 1 1 Equilibrium
More informationEE222 - Spring 16 - Lecture 2 Notes 1
EE222 - Spring 16 - Lecture 2 Notes 1 Murat Arcak January 21 2016 1 Licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License. Essentially Nonlinear Phenomena Continued
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 informationHidden oscillations in dynamical systems
Hidden oscillations in dnamical sems G.A. LEONOV a, N.V. KUZNETSOV b,a, S.M. SELEDZHI a a St.Petersburg State Universit, Universitetsk pr. 28, St.Petersburg, 19854, RUSSIA b Universit of Jväsklä, P.O.
More informationSolutions for B8b (Nonlinear Systems) Fake Past Exam (TT 10)
Solutions for B8b (Nonlinear Systems) Fake Past Exam (TT 10) Mason A. Porter 15/05/2010 1 Question 1 i. (6 points) Define a saddle-node bifurcation and show that the first order system dx dt = r x e x
More informationLYAPUNOV EXPONENTS AND STABILITY FOR THE STOCHASTIC DUFFING-VAN DER POL OSCILLATOR
LYAPUNOV EXPONENTS AND STABILITY FOR THE STOCHASTIC DUFFING-VAN DER POL OSCILLATOR Peter H. Baxendale Department of Mathematics University of Southern California Los Angeles, CA 90089-3 USA baxendal@math.usc.edu
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 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 informationA short introduction with a view toward examples. Short presentation for the students of: Dynamical Systems and Complexity Summer School Volos, 2017
A short introduction with a view toward examples Center of Research and Applications of Nonlinear (CRANS) Department of Mathematics University of Patras Greece sanastassiou@gmail.com Short presentation
More informationLectures on Periodic Orbits
Lectures on Periodic Orbits 11 February 2009 Most of the contents of these notes can found in any typical text on dynamical systems, most notably Strogatz [1994], Perko [2001] and Verhulst [1996]. Complete
More informationASTATISM IN NONLINEAR CONTROL SYSTEMS WITH APPLICATION TO ROBOTICS
dx dt DIFFERENTIAL EQUATIONS AND CONTROL PROCESSES N 1, 1997 Electronic Journal, reg. N P23275 at 07.03.97 http://www.neva.ru/journal e-mail: diff@osipenko.stu.neva.ru Control problems in nonlinear systems
More information1 Lyapunov theory of stability
M.Kawski, APM 581 Diff Equns Intro to Lyapunov theory. November 15, 29 1 1 Lyapunov theory of stability Introduction. Lyapunov s second (or direct) method provides tools for studying (asymptotic) stability
More informationPassivity-based Stabilization of Non-Compact Sets
Passivity-based Stabilization of Non-Compact Sets Mohamed I. El-Hawwary and Manfredi Maggiore Abstract We investigate the stabilization of closed sets for passive nonlinear systems which are contained
More informationAnalytical method for computation of phase-detector characteristic Leonov G.A., Kuznetsov N.V., Yuldashev M.V., Yuldashev R.V.
IEEE Transactions on Circuits and Systems Part II: Express Briefs, 59(0), 0, pp. 633-637 (doi:0.09/tcsii.0.336) Analytical method for computation of phase-detector characteristic Leonov G.A., Kuznetsov
More informationTopological Bifurcations of Knotted Tori in Quasiperiodically Driven Oscillators
Topological Bifurcations of Knotted Tori in Quasiperiodically Driven Oscillators Brian Spears with Andrew Szeri and Michael Hutchings University of California at Berkeley Caltech CDS Seminar October 24,
More informationInput to state Stability
Input to state Stability Mini course, Universität Stuttgart, November 2004 Lars Grüne, Mathematisches Institut, Universität Bayreuth Part IV: Applications ISS Consider with solutions ϕ(t, x, w) ẋ(t) =
More information(8.51) ẋ = A(λ)x + F(x, λ), where λ lr, the matrix A(λ) and function F(x, λ) are C k -functions with k 1,
2.8.7. Poincaré-Andronov-Hopf Bifurcation. In the previous section, we have given a rather detailed method for determining the periodic orbits of a two dimensional system which is the perturbation of a
More informationNBA Lecture 1. Simplest bifurcations in n-dimensional ODEs. Yu.A. Kuznetsov (Utrecht University, NL) March 14, 2011
NBA Lecture 1 Simplest bifurcations in n-dimensional ODEs Yu.A. Kuznetsov (Utrecht University, NL) March 14, 2011 Contents 1. Solutions and orbits: equilibria cycles connecting orbits other invariant sets
More informationProblem Set Number 5, j/2.036j MIT (Fall 2014)
Problem Set Number 5, 18.385j/2.036j MIT (Fall 2014) Rodolfo R. Rosales (MIT, Math. Dept.,Cambridge, MA 02139) Due Fri., October 24, 2014. October 17, 2014 1 Large µ limit for Liénard system #03 Statement:
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 informationNonlinear Control. Nonlinear Control Lecture # 8 Time Varying and Perturbed Systems
Nonlinear Control Lecture # 8 Time Varying and Perturbed Systems Time-varying Systems ẋ = f(t,x) f(t,x) is piecewise continuous in t and locally Lipschitz in x for all t 0 and all x D, (0 D). The origin
More informationThermodynamics of Chaotic systems by C. Beck and F Schlögl (1993) LecturesonGeometryandDynamicalSystemsbyY.PesinandV.Clemenhaga
ME597B/MATH597G/PHYS597C Spring, 2015 Chapter 01: Dynamical Systems Source Thermodynamics of Chaotic systems by C. Beck and F Schlögl (1993) LecturesonGeometryandDynamicalSystemsbyY.PesinandV.Clemenhaga
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 informationMathematical Modeling I
Mathematical Modeling I Dr. Zachariah Sinkala Department of Mathematical Sciences Middle Tennessee State University Murfreesboro Tennessee 37132, USA November 5, 2011 1d systems To understand more complex
More informationApplied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.
Printed Name: Signature: Applied Math Qualifying Exam 11 October 2014 Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. 2 Part 1 (1) Let Ω be an open subset of R
More information... it may happen that small differences in the initial conditions produce very great ones in the final phenomena. Henri Poincaré
Chapter 2 Dynamical Systems... it may happen that small differences in the initial conditions produce very great ones in the final phenomena. Henri Poincaré One of the exciting new fields to arise out
More information14 Periodic phenomena in nature and limit cycles
14 Periodic phenomena in nature and limit cycles 14.1 Periodic phenomena in nature As it was discussed while I talked about the Lotka Volterra model, a great deal of natural phenomena show periodic behavior.
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 informationHigh-Gain Observers in Nonlinear Feedback Control. Lecture # 2 Separation Principle
High-Gain Observers in Nonlinear Feedback Control Lecture # 2 Separation Principle High-Gain ObserversinNonlinear Feedback ControlLecture # 2Separation Principle p. 1/4 The Class of Systems ẋ = Ax + Bφ(x,
More informationK. Pyragas* Semiconductor Physics Institute, LT-2600 Vilnius, Lithuania Received 19 March 1998
PHYSICAL REVIEW E VOLUME 58, NUMBER 3 SEPTEMBER 998 Synchronization of coupled time-delay systems: Analytical estimations K. Pyragas* Semiconductor Physics Institute, LT-26 Vilnius, Lithuania Received
More informationThe nonsmooth pitchfork bifurcation. Glendinning, Paul. MIMS EPrint: Manchester Institute for Mathematical Sciences School of Mathematics
The nonsmooth pitchfork bifurcation Glendinning, Paul 2004 MIMS EPrint: 2006.89 Manchester Institute for Mathematical Sciences School of Mathematics The University of Manchester Reports available from:
More informationHilbert s 16. problem
Hilbert s 16. problem Aarhus University Aarhus 24 th March, 2017 Hilbert s Problem Does there exist (minimal) natural numbers H n, for n = 2, 3,..., such that any planar autonomous system of ODE s dx dt
More informationIn essence, Dynamical Systems is a science which studies differential equations. A differential equation here is the equation
Lecture I In essence, Dynamical Systems is a science which studies differential equations. A differential equation here is the equation ẋ(t) = f(x(t), t) where f is a given function, and x(t) is an unknown
More informationA Model of Evolutionary Dynamics with Quasiperiodic Forcing
paper presented at Society for Experimental Mechanics (SEM) IMAC XXXIII Conference on Structural Dynamics February 2-5, 205, Orlando FL A Model of Evolutionary Dynamics with Quasiperiodic Forcing Elizabeth
More informationConvergence Rate of Nonlinear Switched Systems
Convergence Rate of Nonlinear Switched Systems Philippe JOUAN and Saïd NACIRI arxiv:1511.01737v1 [math.oc] 5 Nov 2015 January 23, 2018 Abstract This paper is concerned with the convergence rate of the
More information2 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 informationPractice Problems for Final Exam
Math 1280 Spring 2016 Practice Problems for Final Exam Part 2 (Sections 6.6, 6.7, 6.8, and chapter 7) S o l u t i o n s 1. Show that the given system has a nonlinear center at the origin. ẋ = 9y 5y 5,
More informationMethod of Generation of Chaos Map in the Centre Manifold
Advanced Studies in Theoretical Physics Vol. 9, 2015, no. 16, 795-800 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/astp.2015.51097 Method of Generation of Chaos Map in the Centre Manifold Evgeny
More information27. Topological classification of complex linear foliations
27. Topological classification of complex linear foliations 545 H. Find the expression of the corresponding element [Γ ε ] H 1 (L ε, Z) through [Γ 1 ε], [Γ 2 ε], [δ ε ]. Problem 26.24. Prove that for any
More informationSynchronization, Chaos, and the Dynamics of Coupled Oscillators. Supplemental 1. Winter Zachary Adams Undergraduate in Mathematics and Biology
Synchronization, Chaos, and the Dynamics of Coupled Oscillators Supplemental 1 Winter 2017 Zachary Adams Undergraduate in Mathematics and Biology Outline: The shift map is discussed, and a rigorous proof
More informationAvailable online at ScienceDirect. Procedia IUTAM 19 (2016 ) IUTAM Symposium Analytical Methods in Nonlinear Dynamics
Available online at www.sciencedirect.com ScienceDirect Procedia IUTAM 19 (2016 ) 11 18 IUTAM Symposium Analytical Methods in Nonlinear Dynamics A model of evolutionary dynamics with quasiperiodic forcing
More information5.2.2 Planar Andronov-Hopf bifurcation
138 CHAPTER 5. LOCAL BIFURCATION THEORY 5.. Planar Andronov-Hopf bifurcation What happens if a planar system has an equilibrium x = x 0 at some parameter value α = α 0 with eigenvalues λ 1, = ±iω 0, ω
More information(Non-)Existence of periodic orbits in dynamical systems
(Non-)Existence of periodic orbits in dynamical systems Konstantin Athanassopoulos Department of Mathematics and Applied Mathematics University of Crete June 3, 2014 onstantin Athanassopoulos (Univ. of
More informationHigh-Gain Observers in Nonlinear Feedback Control
High-Gain Observers in Nonlinear Feedback Control Lecture # 1 Introduction & Stabilization High-Gain ObserversinNonlinear Feedback ControlLecture # 1Introduction & Stabilization p. 1/4 Brief History Linear
More informationBIFURCATION PHENOMENA Lecture 4: Bifurcations in n-dimensional ODEs
BIFURCATION PHENOMENA Lecture 4: Bifurcations in n-dimensional ODEs Yuri A. Kuznetsov August, 2010 Contents 1. Solutions and orbits: equilibria cycles connecting orbits compact invariant manifolds strange
More informationSummary of topics relevant for the final. p. 1
Summary of topics relevant for the final p. 1 Outline Scalar difference equations General theory of ODEs Linear ODEs Linear maps Analysis near fixed points (linearization) Bifurcations How to analyze a
More informationNonlinear Dynamical Systems Lecture - 01
Nonlinear Dynamical Systems Lecture - 01 Alexandre Nolasco de Carvalho August 08, 2017 Presentation Course contents Aims and purpose of the course Bibliography Motivation To explain what is a dynamical
More informationNonlinear Systems Theory
Nonlinear Systems Theory Matthew M. Peet Arizona State University Lecture 2: Nonlinear Systems Theory Overview Our next goal is to extend LMI s and optimization to nonlinear systems analysis. Today we
More informationAVERAGING AND RECONSTRUCTION IN HAMILTONIAN SYSTEMS
AVERAGING AND RECONSTRUCTION IN HAMILTONIAN SYSTEMS Kenneth R. Meyer 1 Jesús F. Palacián 2 Patricia Yanguas 2 1 Department of Mathematical Sciences University of Cincinnati, Cincinnati, Ohio (USA) 2 Departamento
More informationNonlinear Control. Nonlinear Control Lecture # 3 Stability of Equilibrium Points
Nonlinear Control Lecture # 3 Stability of Equilibrium Points The Invariance Principle Definitions Let x(t) be a solution of ẋ = f(x) A point p is a positive limit point of x(t) if there is a sequence
More informationCourse 212: Academic Year Section 1: Metric Spaces
Course 212: Academic Year 1991-2 Section 1: Metric Spaces D. R. Wilkins Contents 1 Metric Spaces 3 1.1 Distance Functions and Metric Spaces............. 3 1.2 Convergence and Continuity in Metric Spaces.........
More informationComputers and Mathematics with Applications. Chaos suppression via periodic pulses in a class of piece-wise continuous systems
Computers and Mathematics with Applications 64 (2012) 849 855 Contents lists available at SciVerse ScienceDirect Computers and Mathematics with Applications journal homepage: www.elsevier.com/locate/camwa
More informationStability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games
Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,
More informationStability of Dynamical systems
Stability of Dynamical systems Stability Isolated equilibria Classification of Isolated Equilibria Attractor and Repeller Almost linear systems Jacobian Matrix Stability Consider an autonomous system u
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 informationZERO-HOPF BIFURCATION FOR A CLASS OF LORENZ-TYPE SYSTEMS
This is a preprint of: Zero-Hopf bifurcation for a class of Lorenz-type systems, Jaume Llibre, Ernesto Pérez-Chavela, Discrete Contin. Dyn. Syst. Ser. B, vol. 19(6), 1731 1736, 214. DOI: [doi:1.3934/dcdsb.214.19.1731]
More informationConverse Lyapunov theorem and Input-to-State Stability
Converse Lyapunov theorem and Input-to-State Stability April 6, 2014 1 Converse Lyapunov theorem In the previous lecture, we have discussed few examples of nonlinear control systems and stability concepts
More informationHigher Order Averaging : periodic solutions, linear systems and an application
Higher Order Averaging : periodic solutions, linear systems and an application Hartono and A.H.P. van der Burgh Faculty of Information Technology and Systems, Department of Applied Mathematical Analysis,
More informationProf. Krstic Nonlinear Systems MAE281A Homework set 1 Linearization & phase portrait
Prof. Krstic Nonlinear Systems MAE28A Homework set Linearization & phase portrait. For each of the following systems, find all equilibrium points and determine the type of each isolated equilibrium. Use
More information... le souci du beau nous conduit aux mêmes choix que celui de l utile. Henri Poincaré
There is nothing more practical than a good theory. James C. Maxwell... le souci du beau nous conduit aux mêmes choix que celui de l utile. Henri Poincaré Preface Many phenomena in science and technology
More informationCDS Solutions to the Midterm Exam
CDS 22 - Solutions to the Midterm Exam Instructor: Danielle C. Tarraf November 6, 27 Problem (a) Recall that the H norm of a transfer function is time-delay invariant. Hence: ( ) Ĝ(s) = s + a = sup /2
More informationNonlinear Control Lecture 1: Introduction
Nonlinear Control Lecture 1: Introduction Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Fall 2011 Farzaneh Abdollahi Nonlinear Control Lecture 1 1/15 Motivation
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 informationAn introduction to Mathematical Theory of Control
An introduction to Mathematical Theory of Control Vasile Staicu University of Aveiro UNICA, May 2018 Vasile Staicu (University of Aveiro) An introduction to Mathematical Theory of Control UNICA, May 2018
More informationNonlinear Control Lecture # 1 Introduction. Nonlinear Control
Nonlinear Control Lecture # 1 Introduction Nonlinear State Model ẋ 1 = f 1 (t,x 1,...,x n,u 1,...,u m ) ẋ 2 = f 2 (t,x 1,...,x n,u 1,...,u m ).. ẋ n = f n (t,x 1,...,x n,u 1,...,u m ) ẋ i denotes the derivative
More informationAnalytic Exact Upper Bound for the Lyapunov Dimension of the Shimizu Morioka System
Entropy 25, 7, 5-56; doi:.339/e775 OPEN ACCESS entropy ISSN 99-43 www.mdpi.com/journal/entropy Article Analytic Exact Upper Bound for the Lyapunov Dimension of the Shimizu Morioka System Gennady A. Leonov,
More informationStability of Parameter Adaptation Algorithms. Big picture
ME5895, UConn, Fall 215 Prof. Xu Chen Big picture For ˆθ (k + 1) = ˆθ (k) + [correction term] we haven t talked about whether ˆθ(k) will converge to the true value θ if k. We haven t even talked about
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 informationCONTROLLING IN BETWEEN THE LORENZ AND THE CHEN SYSTEMS
International Journal of Bifurcation and Chaos, Vol. 12, No. 6 (22) 1417 1422 c World Scientific Publishing Company CONTROLLING IN BETWEEN THE LORENZ AND THE CHEN SYSTEMS JINHU LÜ Institute of Systems
More informationA Study of the Van der Pol Equation
A Study of the Van der Pol Equation Kai Zhe Tan, s1465711 September 16, 2016 Abstract The Van der Pol equation is famous for modelling biological systems as well as being a good model to study its multiple
More informationDepartment of Mathematics IIT Guwahati
Stability of Linear Systems in R 2 Department of Mathematics IIT Guwahati A system of first order differential equations is called autonomous if the system can be written in the form dx 1 dt = g 1(x 1,
More information