Complicated behavior of dynamical systems. Mathematical methods and computer experiments.

Size: px
Start display at page:

Download "Complicated behavior of dynamical systems. Mathematical methods and computer experiments."

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

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

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

More information

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

WHAT IS A CHAOTIC ATTRACTOR?

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

More information

B5.6 Nonlinear Systems

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

Lyapunov Stability Theory

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

SIMPLE CHAOTIC FLOWS WITH ONE STABLE EQUILIBRIUM

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

LECTURE 8: DYNAMICAL SYSTEMS 7

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

BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs

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

Poincaré Map, Floquet Theory, and Stability of Periodic Orbits

Poincaré 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 information

Stabilization of Hyperbolic Chaos by the Pyragas Method

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

Stability of Nonlinear Systems An Introduction

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

MCE693/793: Analysis and Control of Nonlinear Systems

MCE693/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 information

Recent new examples of hidden attractors

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

Example of a Blue Sky Catastrophe

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

Nonlinear Analysis and Design of Phase-Locked Loops (PLL)

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

There 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)

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

2 Discrete growth models, logistic map (Murray, Chapter 2)

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

Introduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems

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

B5.6 Nonlinear Systems

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

CHALMERS, 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 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 information

Edward Lorenz. Professor of Meteorology at the Massachusetts Institute of Technology

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

EN Nonlinear Control and Planning in Robotics Lecture 3: Stability February 4, 2015

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

Handout 2: Invariant Sets and Stability

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

LMI Methods in Optimal and Robust Control

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

TWO DIMENSIONAL FLOWS. Lecture 5: Limit Cycles and Bifurcations

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

Lecture 4. Chapter 4: Lyapunov Stability. Eugenio Schuster. Mechanical Engineering and Mechanics Lehigh University.

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

3 Stability and Lyapunov Functions

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

STABILITY. Phase portraits and local stability

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

Introduction to Applied Nonlinear Dynamical Systems and Chaos

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

EE222 - Spring 16 - Lecture 2 Notes 1

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

PHY411 Lecture notes Part 5

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

Hidden oscillations in dynamical systems

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

Solutions for B8b (Nonlinear Systems) Fake Past Exam (TT 10)

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

LYAPUNOV EXPONENTS AND STABILITY FOR THE STOCHASTIC DUFFING-VAN DER POL OSCILLATOR

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

April 13, We now extend the structure of the horseshoe to more general kinds of invariant. x (v) λ n v.

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

Introduction to Dynamical Systems Basic Concepts of Dynamics

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

More information

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

Lectures on Periodic Orbits

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

ASTATISM IN NONLINEAR CONTROL SYSTEMS WITH APPLICATION TO ROBOTICS

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

1 Lyapunov theory of stability

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

Passivity-based Stabilization of Non-Compact Sets

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

Analytical method for computation of phase-detector characteristic Leonov G.A., Kuznetsov N.V., Yuldashev M.V., Yuldashev R.V.

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

Topological Bifurcations of Knotted Tori in Quasiperiodically Driven Oscillators

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

Input to state Stability

Input 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,

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

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

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

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

11 Chaos in Continuous Dynamical Systems.

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

Nonlinear Control. Nonlinear Control Lecture # 8 Time Varying and Perturbed Systems

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

Thermodynamics of Chaotic systems by C. Beck and F Schlögl (1993) LecturesonGeometryandDynamicalSystemsbyY.PesinandV.Clemenhaga

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

THREE DIMENSIONAL SYSTEMS. Lecture 6: The Lorenz Equations

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

Mathematical Modeling I

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

Applied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.

Applied 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é

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

14 Periodic phenomena in nature and limit cycles

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

One dimensional Maps

One dimensional Maps Chapter 4 One dimensional Maps The ordinary differential equation studied in chapters 1-3 provide a close link to actual physical systems it is easy to believe these equations provide at least an approximate

More information

High-Gain Observers in Nonlinear Feedback Control. Lecture # 2 Separation Principle

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

K. Pyragas* Semiconductor Physics Institute, LT-2600 Vilnius, Lithuania Received 19 March 1998

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

The nonsmooth pitchfork bifurcation. Glendinning, Paul. MIMS EPrint: Manchester Institute for Mathematical Sciences School of Mathematics

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

Hilbert s 16. problem

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

In essence, Dynamical Systems is a science which studies differential equations. A differential equation here is the equation

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

A Model of Evolutionary Dynamics with Quasiperiodic Forcing

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

Convergence Rate of Nonlinear Switched Systems

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

2 Qualitative theory of non-smooth dynamical systems

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

More information

Practice Problems for Final Exam

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

Method of Generation of Chaos Map in the Centre Manifold

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

27. Topological classification of complex linear foliations

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

Synchronization, 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 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 information

Available online at ScienceDirect. Procedia IUTAM 19 (2016 ) IUTAM Symposium Analytical Methods in Nonlinear Dynamics

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

5.2.2 Planar Andronov-Hopf bifurcation

5.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 (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 information

High-Gain Observers in Nonlinear Feedback Control

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

BIFURCATION PHENOMENA Lecture 4: Bifurcations in n-dimensional ODEs

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

Summary of topics relevant for the final. p. 1

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

Nonlinear Dynamical Systems Lecture - 01

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

Nonlinear Systems Theory

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

AVERAGING AND RECONSTRUCTION IN HAMILTONIAN SYSTEMS

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

Nonlinear Control. Nonlinear Control Lecture # 3 Stability of Equilibrium Points

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

Course 212: Academic Year Section 1: Metric Spaces

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

Computers and Mathematics with Applications. Chaos suppression via periodic pulses in a class of piece-wise continuous systems

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

Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games

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

Stability of Dynamical systems

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

Lecture 6. Lorenz equations and Malkus' waterwheel Some properties of the Lorenz Eq.'s Lorenz Map Towards definitions of:

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

ZERO-HOPF BIFURCATION FOR A CLASS OF LORENZ-TYPE SYSTEMS

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

Converse Lyapunov theorem and Input-to-State Stability

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

Higher Order Averaging : periodic solutions, linear systems and an application

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

Prof. Krstic Nonlinear Systems MAE281A Homework set 1 Linearization & phase portrait

Prof. 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é

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

CDS Solutions to the Midterm Exam

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

Nonlinear Control Lecture 1: Introduction

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

Simple approach to the creation of a strange nonchaotic attractor in any chaotic system

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

An introduction to Mathematical Theory of Control

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

Nonlinear Control Lecture # 1 Introduction. Nonlinear Control

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

Analytic Exact Upper Bound for the Lyapunov Dimension of the Shimizu Morioka System

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

Stability of Parameter Adaptation Algorithms. Big picture

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

Chapter 23. Predicting Chaos The Shift Map and Symbolic Dynamics

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

CONTROLLING IN BETWEEN THE LORENZ AND THE CHEN SYSTEMS

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

A Study of the Van der Pol Equation

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

Department of Mathematics IIT Guwahati

Department 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