# ME 680- Spring Geometrical Analysis of 1-D Dynamical Systems

Size: px
Start display at page:

Download "ME 680- Spring Geometrical Analysis of 1-D Dynamical Systems"

## Transcription

1 ME 680- Spring 2014 Geometrical Analysis of 1-D Dynamical Systems 1

2 Geometrical Analysis of 1-D Dynamical Systems Logistic equation: n = rn(1 n) velocity function Equilibria or fied points : initial conditions n where you start and stay without evolving for all time. They correspond to zeros of the velocity function: f(n) Phase diagram n =0 n =1 n The length of the arrows magnitude of the velocity (function) at that point. 2

3 α limit set of a point (initial condition) n 0 : It is defined as the set of limit points of the trajectory started at n 0, for t -. Thus, α(n 0 ) = n lim φ (n 0, t) = n t ω limit set of a point n 0 is the set ω(n 0 ) = n lim φ (n 0, t) = n t + Eistence of a potential function: Consider n = f(n) (Gradient Dynamical System) Let there be a function V(n) such that f(n) = V n Eample: for the Logistic equation f(n) = r n (1 n), i.e., V(n) = r n r n 3 3 3

4 Then, note that the equilibrium points for the system (a Gradient Dynamical System) are at the local etrema of the potential function. This is where the similarity with mechanical systems with potential energy functions ends!! Considering the Logistic equation: V(n) = rn2 2 + rn3 3 the plot of the potential function, and the equilibrium points are as follows: V(n) n=0 n=1 n 4

5 Observations: Oscillatory behavior is not possible in 1-D autonomous Systems Trajectories approach the equilibrium point n =1, but never reach it in finite time. Invariant subspaces are regions I R in phase space where if n 0 I, then φ t, n 0 I for all negative and positive flow times (- < t < ). For the Logistic Equation, the invariant subspaces are: I 1 =, 0, I 2 = 0, I 3 = 0, 1, I 4 = 1, I 5 = {1, + } Thus, the state space is decomposed into: R = I 1 I 2 I 3 I 4 I 5 α and ω limit sets of any initial condition n 0 : 5

6 If n 0 I, the α limit set is 0 ω limit set is If n 0 I 2, then α and ω limit sets are the same If n 0 I 3, then α (n 0 ) = {0}, and ω(n 0 ) = {1} so on. Stability of Equilibria/Fied Points An equilibrium point of = f(), say =, is stable if ε > 0, δ (ε) > 0 such that for any initial condition 0, with 0 < δ, φ (t, 0 ) < ε for all t, 0 < t < +. Otherwise, it is unstable. 6

7 δ 0 This definition of stability is very difficult to use directly to deduce stability of an equilibrium point. One needs to a priori know the solution for every given initial condition starting inside the region of size δ. Thus, one really needs to find other criteria that can be used to characterize stability without solving the differential equation. If in addition, φ(t, 0 ) 0 as t, then is an asymptotically stable equilibrium. 7

8 Eamples: 1. Logistic Equation f(n)=rn(1 n) n =0 (unstable) n =1 n (asymptotically stable) 2. Quadratic System f(n)=an 2 3. Cubic System f(n) = an 3 n =0 (unstable) n n =0 (stable) n Observe: Any isolated stable equilibrium in 1-D autonomous systems has to be asymptotically stable. 8

9 Linearization about equilibrium points Let be a fied point of = f(), i.e. f( ) = 0 To linearize about =, introduce a perturbation: Let = Then + = f( + ) df or f( ) d 0 0 (Taylor series epansion for small X ) df d 0 This is the linearized equation about = 9

10 Eample: Logistic Equation: n = r n(1 n) The equilibrium points are n = 1 and n = 0. Let us linearize the system about n = 1 Then n = n + n = 1 + n and n = r(1 + n) ( n) = nr rn 2 df r n r n dn n 0 This is the linearized system near n = 1. Note that n = 1 is linearly stable. We can make the connection between linear stability (i.e. stability of equilibrium for the linearized system) and nonlinear stability if (only if) df dn n0 0 (Hartman-Grobman theorem) 10

11 Closing Remarks on Linearization = f() ; f( ) = 0 is an equilibrium There are a few ways to linearize the system. (i): df f( ) ( ) d (Taylor series epansion) df f( ) ( ) d Let. Then df d linearized system around an Equilibrium 11

12 (ii): Let = +. Then f( ) f() ˆ dfˆ f(0) ˆ d 0 dfˆ or d 0 General solution of dfˆ d dfˆ (0)ep( t) d eigenvalue if eigenvalue < 0, = is asymptotically stable if eigenvalue > 0, = is unstable is 12

13 If eigenvalue is 0, the equilibrium is called hyperbolic. Otherwise, it is called non-hyperbolic. According to the Hartman-Grobman theorem, if is a hyperbolic equilibrium, stability conclusions drawn from linearized equation (linear stability) hold also for the nonlinear model (nonlinear stability) if df d 0,, then we have to look at higher order terms in the Taylor series to judge stability. 13

14 Bifurcations of equilibria in 1-D Interesting dynamics can occur as system (or control) parameters vary: Equilibria can suddenly change in number or stability type. E: Consider the eample of a cantilever beam with a mass on top, with the mass being a control parameter: For mg < P cr (1 equilibrium) For mg > P cr (3 equilibriums) mg g g 14

15 Saddle-node bifurcation (fold, or turning point, blue sky bifurcation) A prototypical eample: = r + 2 here r is some control parameter The velocity functions for three distinct cases are as follows: r r < 0 two equil r = 0 one equil r > 0 no equil 15

16 We can present these results in a diagram of equilibrium solutions as a function of the parameter r. X unstable r =0 r stable This is a bifurcation diagram. (r = 0, = 0) is the bifurcation point. This is called a subcritical saddle - node bifurcation. 16

17 Supercritical saddle node bifurcation: Consider the system: 2 r f() X r =0 stable unstable r Linear Stability Analysis r > 0 : the equilibrium points are df The function derivative is d r For df r, we get 2 r ( asympt stable) d 17

18 df For r, we get 2 r (unstable) d df Note that both equilibria are hyperbolic 0 d At r = 0, however, theorem fails! df d 0 Consider the velocity function at r = 0: 0 i.e., Hartman-Grobman The equilibrium at =0 is actually unstable! r = 0 one equilibrium 18

19 Another eample Consider a system governed by : r e f() We can determine the equilibria and find their stability via linearization: r < r r = r r > r 1 e - 1 e - 1 e - r (r-) r (r-) (r-) What is the critical value of r? At critical value, and r must satisfy f( ) = 0 r - - e - = 0 as well as r df d 0 1 e 0 0 r 1 19

20 Brief Introduction to Normal Forms In a sense, f() = r ± 2 are prototypical of all 1-D systems undergoing a saddle-node bifurcation. Consider the system just studied: f (,r) r e with r 1, 0. Near the critical point, r r r 1 r and for small r r and, write 2 2 higher order 1 r (1 / 2 ) r / 2 terms f(,r) increasing r for r=r same form as that of super-critical saddle node bifurcation 20

21 f() = a + b 2 is the normal form of saddle - node bifurcation, i.e., all systems in 1-D undergoing this bifurcation must locally possess this form. Transcritical Bifurcation The normal form for this bifurcation is = r 2 (similar to n = r n (1 n), the logistic equation). Consider the velocity function for different parameter values: r < 0 two equil r = 0 one equil r > 0 two equil 21

22 We can display the results in the form of a bifurcation diagram: =r =0 r=0 r This is called a transcritical bifurcation Eample: Lasers. See notes. 22

23 Eample: Laser threshold pump partially reflecting mirror Active material laser light At low energy levels each atom oscillates acting as a little antenna, but all atoms oscillate independently and emit randomly phased photons. At a threshold pumping level, all the atoms oscillate in phase producing laser! This is due to self-organization out of cooperative interaction of atoms. (Ref: Haken 1983, Strogatz s book) 23

24 Let n(t) - no. of photons Then, n = gain loss = G n N k n (escape or leakage thru endface) gain coeff > 0 no. of ecited atoms Note that k > 0, a rate constant Here τ = 1 k = typical life time of a photon in the laser Note however that N(t) = N o α n (because atoms after radiation of a photon, are not in an ecited state), i.e., n Gn(N n) kn or 2 o n (GNo k)n Gn 24

25 The corresponding bifurcation diagram is: lamp n laser =r n =0 N 0 =k/g N 0 No physical meaning 25

26 Pitchfork bifurcation Eamples: We have already seen the eample of buckling of a column as a function of the aial load: g mg g fluid Another eample is that of the onset of convection in a toroidal thermosyphan heating coil 26

27 The normal form for pitchfork bifurcation is: = r 3 = f( The behavior can be understood in terms of the velocity functions as follows: f() f() f() r < 0 r = 0 X r > 0 stable The bifurcation diagram is then: Supercritical pitchfork stable r =0 unstable stable r 27

28 Linear stability analysis Consider = r 3 = f() The equilibria are at 0, r df d 0 0 r 3(0) r is stable when r < 0 is unstable r > 0 what about when r = 0? The linear analysis fails!! For the non-zero equilibria: df d eigenvalue is negative if r > 0 2 i.e., these bifurcating equilbiria are asymp. stable. r 2 r 3( r ) 2r 28

29 Subcritical pitchfork The normal form is r, with equilibria 0, r The resulting bifurcation diagram is: 3 unstable X stable r =0 unstable r unstable 29

30 Usually, the unstable behavior is stabilized by higher order non-linear terms, e.g., = r The resulting bifurcation diagram can be shown to be: unstable X stable r =0 unstable r unstable subcritical pitchfork, r = r P supercritical saddle-node, r = r S 30

31 Connection between simple bifurcations and the implicit function theorem Consider the system f(,r) Let f( 0, r 0 ) = 0 i. e., ( 0, r 0 ) be an equilibrium. Let f be continuously differentiable w.r.t. and r in some open region in the (, r) plane containing ( 0, r 0 ). Then if df 0 in a small neighborhood of ( 0, r 0 ), d (,r ) 0 0 we must have: f(, r) = 0 has a unique solution =(r) such that f((r),r)=0 furthermore, (r) is also continuously differentiable. No bifurcations arise so long as df 0 d (,r )

32 The figure below illustrates the idea through two points along a solution curve. At ( 1,r 1 ), the derivative df/d does not vanish, where as at ( 2,r 2 ), the derivative df/d vanishes. df/d0 ( 2,r 2 ) ( 1,r 1 ) df/d=0 r 32

33 Imperfect Bifurcations & Catastrophes Consider the buckling eample. If the load does not coincide with the ais of the column, what happens? g mg g mg symmetric loading asymmetric loading Real physical systems have imperfections and mathematical imposition of reflection symmetry is an idealization. Do the bifurcation diagrams change significantly if imperfections or perturbations are added to the model (velocity function)? This is related to the concept of structural stability or robustness of models. 33

34 Consider the two-parameter normal form: h r 3 we now have the two parameters, h and r. Note that it is a perturbation of the normal form for pitchfork bifurcation (h=0) y h y r 3 y h, h h (r) h c h (r) c y r r 0, only one equilibrium possible for any h 3 h h (r) c 3 equilibria in t his region h h (r) c r > 0, one or three equilibria possible 34

35 h increasing Imperfect bifurcations and catastrophes (cont d). Consider = h + r 3 = f( (h = 0 normal form of pitchfork) y() = r 3 with h We look for intersections of y h 3 y r y h, h h c(r) h h (r) c y r r 0, only one equilibrium possible for any h 3 r C / 3 h r > 0, one or three equilibria possible h (r) c 35

36 dy 2 r For critical point 0 r c c 3c 0 c d 3 Furthermore, 3 2r c r h c c rc c c 0 hc 3 3 h h=0 only one equil soln cusp r=0 r C equil solns equil solns r 36

37 1. 0 r C r 2. r C 0 r 3. 0 r C r 37

38 4. 0 h=0 h Alternately, in 3-D we can visualize the solutions set as follows: h h=0 catastrophic surface r=0 r 38

39 1-D system on a circle: [over-damped pendulum (acted by a constant torque )] The equation of motion is: O θ l g 2 ml b mglsin m e Let ml 2 θ e R be negligible (imagine pendulum in a vat of molasses) The resulting equation is b mglsin b or mgl sin mgl 39

40 mgl Let t ; b mgl (ratio of appl. torque to ma. gravitational torque) d Then sin where d We say that [0, 2 ] i.e., the phase space is a circle θ = θ θ =0 Consider the system: sin 40

41 if γ > 1, pendulum goes around the circle albeit non-uniformly If γ = 1, θ = π 2 is an equilibrium θ =/2 θ = θ =0 1, If 1, there are two equilibria : 1 and 2 which have opposite stability characteristics θ = θ 2 Clearly, there is a saddle-node bifurcation at γ = 1. θ 1 θ =0 41

### One Dimensional Dynamical Systems

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

### ONE DIMENSIONAL FLOWS. Lecture 3: Bifurcations

ONE DIMENSIONAL FLOWS Lecture 3: Bifurcations 3. Bifurcations Here we show that, although the dynamics of one-dimensional systems is very limited [all solutions either settle down to a steady equilibrium

### = F ( x; µ) (1) where x is a 2-dimensional vector, µ is a parameter, and F :

1 Bifurcations Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University Tallahassee, Florida 32306 A bifurcation is a qualitative change

### 7 Two-dimensional bifurcations

7 Two-dimensional bifurcations As in one-dimensional systems: fixed points may be created, destroyed, or change stability as parameters are varied (change of topological equivalence ). In addition closed

### 10 Back to planar nonlinear systems

10 Back to planar nonlinear sstems 10.1 Near the equilibria Recall that I started talking about the Lotka Volterra model as a motivation to stud sstems of two first order autonomous equations of the form

### Autonomous Mechanical systems of 1 degree of freedom

Autonomous Mechanical systems of degree of freedom = f (, ), = y = y = f (, y) = f (, ) : position y : velocity f : force function Solution : = ( t; c, c), y = y( t; c, c ), ci = ci(, y) Integral of motion

### 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

### Clearly the passage of an eigenvalue through to the positive real half plane leads to a qualitative change in the phase portrait, i.e.

Bifurcations We have already seen how the loss of stiffness in a linear oscillator leads to instability. In a practical situation the stiffness may not degrade in a linear fashion, and instability may

### Problem Set Number 02, j/2.036j MIT (Fall 2018)

Problem Set Number 0, 18.385j/.036j MIT (Fall 018) Rodolfo R. Rosales (MIT, Math. Dept., room -337, Cambridge, MA 0139) September 6, 018 Due October 4, 018. Turn it in (by 3PM) at the Math. Problem Set

### 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

### Solution to Homework #5 Roy Malka 1. Questions 2,3,4 of Homework #5 of M. Cross class. dv (x) dx

Solution to Homework #5 Roy Malka. Questions 2,,4 of Homework #5 of M. Cross class. Bifurcations: (M. Cross) Consider the bifurcations of the stationary solutions of a particle undergoing damped one dimensional

### 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),

### 4 Second-Order Systems

4 Second-Order Systems Second-order autonomous systems occupy an important place in the study of nonlinear systems because solution trajectories can be represented in the plane. This allows for easy visualization

### 8.1 Bifurcations of Equilibria

1 81 Bifurcations of Equilibria Bifurcation theory studies qualitative changes in solutions as a parameter varies In general one could study the bifurcation theory of ODEs PDEs integro-differential equations

### Lecture 3. Dynamical Systems in Continuous Time

Lecture 3. Dynamical Systems in Continuous Time University of British Columbia, Vancouver Yue-Xian Li November 2, 2017 1 3.1 Exponential growth and decay A Population With Generation Overlap Consider a

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

Problem Set Number 2, 18.385j/2.036j MIT (Fall 2014) Rodolfo R. Rosales (MIT, Math. Dept.,Cambridge, MA 02139) Due Mon., September 29, 2014. 1 Inverse function problem #01. Statement: Inverse function

### Chapter 24 BIFURCATIONS

Chapter 24 BIFURCATIONS Abstract Keywords: Phase Portrait Fixed Point Saddle-Node Bifurcation Diagram Codimension-1 Hysteresis Hopf Bifurcation SNIC Page 1 24.1 Introduction In linear systems, responses

### Chaos. Lendert Gelens. KU Leuven - Vrije Universiteit Brussel Nonlinear dynamics course - VUB

Chaos Lendert Gelens KU Leuven - Vrije Universiteit Brussel www.gelenslab.org Nonlinear dynamics course - VUB Examples of chaotic systems: the double pendulum? θ 1 θ θ 2 Examples of chaotic systems: the

### 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

### Nonlinear dynamics & chaos BECS

Nonlinear dynamics & chaos BECS-114.7151 Phase portraits Focus: nonlinear systems in two dimensions General form of a vector field on the phase plane: Vector notation: Phase portraits Solution x(t) describes

### APPPHYS217 Tuesday 25 May 2010

APPPHYS7 Tuesday 5 May Our aim today is to take a brief tour of some topics in nonlinear dynamics. Some good references include: [Perko] Lawrence Perko Differential Equations and Dynamical Systems (Springer-Verlag

### Consider a particle in 1D at position x(t), subject to a force F (x), so that mẍ = F (x). Define the kinetic energy to be.

Chapter 4 Energy and Stability 4.1 Energy in 1D Consider a particle in 1D at position x(t), subject to a force F (x), so that mẍ = F (x). Define the kinetic energy to be T = 1 2 mẋ2 and the potential energy

### Nonlinear Oscillations and Chaos

CHAPTER 4 Nonlinear Oscillations and Chaos 4-. l l = l + d s d d l l = l + d m θ m (a) (b) (c) The unetended length of each spring is, as shown in (a). In order to attach the mass m, each spring must be

### 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

### Chapter 6 Nonlinear Systems and Phenomena. Friday, November 2, 12

Chapter 6 Nonlinear Systems and Phenomena 6.1 Stability and the Phase Plane We now move to nonlinear systems Begin with the first-order system for x(t) d dt x = f(x,t), x(0) = x 0 In particular, consider

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

### Lecture 5. Outline: Limit Cycles. Definition and examples How to rule out limit cycles. Poincare-Bendixson theorem Hopf bifurcations Poincare maps

Lecture 5 Outline: Limit Cycles Definition and examples How to rule out limit cycles Gradient systems Liapunov functions Dulacs criterion Poincare-Bendixson theorem Hopf bifurcations Poincare maps Limit

### 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

### A plane autonomous system is a pair of simultaneous first-order differential equations,

Chapter 11 Phase-Plane Techniques 11.1 Plane Autonomous Systems A plane autonomous system is a pair of simultaneous first-order differential equations, ẋ = f(x, y), ẏ = g(x, y). This system has an equilibrium

### MATH 614 Dynamical Systems and Chaos Lecture 24: Bifurcation theory in higher dimensions. The Hopf bifurcation.

MATH 614 Dynamical Systems and Chaos Lecture 24: Bifurcation theory in higher dimensions. The Hopf bifurcation. Bifurcation theory The object of bifurcation theory is to study changes that maps undergo

### Chapter 7. Nonlinear Systems. 7.1 Introduction

Nonlinear Systems Chapter 7 The scientist does not study nature because it is useful; he studies it because he delights in it, and he delights in it because it is beautiful. - Jules Henri Poincaré (1854-1912)

### 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

### NONLINEAR DYNAMICS AND CHAOS. Numerical integration. Stability analysis

LECTURE 3: FLOWS NONLINEAR DYNAMICS AND CHAOS Patrick E McSharr Sstems Analsis, Modelling & Prediction Group www.eng.o.ac.uk/samp patrick@mcsharr.net Tel: +44 83 74 Numerical integration Stabilit analsis

### 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

### Lecture 3 : Bifurcation Analysis

Lecture 3 : Bifurcation Analysis D. Sumpter & S.C. Nicolis October - December 2008 D. Sumpter & S.C. Nicolis General settings 4 basic bifurcations (as long as there is only one unstable mode!) steady state

### 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

### 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: April 06, 2018, at 14 00 18 00 Johanneberg Kristian

### Stability Analysis of a Single-Degree-of-Freedom Mechanical Model with Distinct Critical Points: I. Bifurcation Theory Approach

World Journal of Mechanics, 013, 3, 6-81 doi:10.436/wjm.013.31005 Published Online February 013 (http://www.scirp.org/journal/wjm) Stability Analysis of a Single-Degree-of-Freedom Mechanical Model with

### Dynamical Systems in Neuroscience: Elementary Bifurcations

Dynamical Systems in Neuroscience: Elementary Bifurcations Foris Kuang May 2017 1 Contents 1 Introduction 3 2 Definitions 3 3 Hodgkin-Huxley Model 3 4 Morris-Lecar Model 4 5 Stability 5 5.1 Linear ODE..............................................

### CHAPTER 6 HOPF-BIFURCATION IN A TWO DIMENSIONAL NONLINEAR DIFFERENTIAL EQUATION

CHAPTER 6 HOPF-BIFURCATION IN A TWO DIMENSIONAL NONLINEAR DIFFERENTIAL EQUATION [Discussion on this chapter is based on our paper entitled Hopf-Bifurcation Ina Two Dimensional Nonlinear Differential Equation,

### arxiv: v1 [physics.class-ph] 5 Jan 2012

Damped bead on a rotating circular hoop - a bifurcation zoo Shovan Dutta Department of Electronics and Telecommunication Engineering, Jadavpur University, Calcutta 700 032, India. Subhankar Ray Department

### Physics 235 Chapter 4. Chapter 4 Non-Linear Oscillations and Chaos

Chapter 4 Non-Linear Oscillations and Chaos Non-Linear Differential Equations Up to now we have considered differential equations with terms that are proportional to the acceleration, the velocity, and

### Nonlinear Autonomous Systems of Differential

Chapter 4 Nonlinear Autonomous Systems of Differential Equations 4.0 The Phase Plane: Linear Systems 4.0.1 Introduction Consider a system of the form x = A(x), (4.0.1) where A is independent of t. Such

### Nonlinear Autonomous Dynamical systems of two dimensions. Part A

Nonlinear Autonomous Dynamical systems of two dimensions Part A Nonlinear Autonomous Dynamical systems of two dimensions x f ( x, y), x(0) x vector field y g( xy, ), y(0) y F ( f, g) 0 0 f, g are continuous

### 1. < 0: the eigenvalues are real and have opposite signs; the fixed point is a saddle point

Solving a Linear System τ = trace(a) = a + d = λ 1 + λ 2 λ 1,2 = τ± = det(a) = ad bc = λ 1 λ 2 Classification of Fixed Points τ 2 4 1. < 0: the eigenvalues are real and have opposite signs; the fixed point

### 2 Lecture 2: Amplitude equations and Hopf bifurcations

Lecture : Amplitude equations and Hopf bifurcations This lecture completes the brief discussion of steady-state bifurcations by discussing vector fields that describe the dynamics near a bifurcation. From

### Fundamentals of Dynamical Systems / Discrete-Time Models. Dr. Dylan McNamara people.uncw.edu/ mcnamarad

Fundamentals of Dynamical Systems / Discrete-Time Models Dr. Dylan McNamara people.uncw.edu/ mcnamarad Dynamical systems theory Considers how systems autonomously change along time Ranges from Newtonian

### PHY411 Lecture notes Part 4

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

### Boulder School for Condensed Matter and Materials Physics. Laurette Tuckerman PMMH-ESPCI-CNRS

Boulder School for Condensed Matter and Materials Physics Laurette Tuckerman PMMH-ESPCI-CNRS laurette@pmmh.espci.fr Dynamical Systems: A Basic Primer 1 1 Basic bifurcations 1.1 Fixed points and linear

### Achieve asymptotic stability using Lyapunov's second method

IOSR Journal of Mathematics (IOSR-JM) e-issn: 78-578, p-issn: 39-765X Volume 3, Issue Ver I (Jan - Feb 07), PP 7-77 wwwiosrjournalsorg Achieve asymptotic stability using Lyapunov's second method Runak

### 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

### Bifurcations. Chapter Types of Bifurcations Saddle-node bifurcation

Chapter 2 Bifurcations 2.1 Types of Bifurcations 2.1.1 Saddle-node bifurcation We remarked above how f (u) is in general nonzero when f(u) itself vanishes, since two equations in a single unknown is an

### 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

### 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

### 8 Autonomous first order ODE

8 Autonomous first order ODE There are different ways to approach differential equations. Prior to this lecture we mostly dealt with analytical methods, i.e., with methods that require a formula as a final

### 1 The pendulum equation

Math 270 Honors ODE I Fall, 2008 Class notes # 5 A longer than usual homework assignment is at the end. The pendulum equation We now come to a particularly important example, the equation for an oscillating

### 2.10 Saddles, Nodes, Foci and Centers

2.10 Saddles, Nodes, Foci and Centers In Section 1.5, a linear system (1 where x R 2 was said to have a saddle, node, focus or center at the origin if its phase portrait was linearly equivalent to one

### Phase Synchronization

Phase Synchronization Lecture by: Zhibin Guo Notes by: Xiang Fan May 10, 2016 1 Introduction For any mode or fluctuation, we always have where S(x, t) is phase. If a mode amplitude satisfies ϕ k = ϕ k

### 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

### 8 Ecosystem stability

8 Ecosystem stability References: May [47], Strogatz [48]. In these lectures we consider models of populations, with an emphasis on the conditions for stability and instability. 8.1 Dynamics of a single

### MOL 410/510: Introduction to Biological Dynamics Fall 2012 Problem Set #4, Nonlinear Dynamical Systems (due 10/19/2012) 6 MUST DO Questions, 1

MOL 410/510: Introduction to Biological Dynamics Fall 2012 Problem Set #4, Nonlinear Dynamical Systems (due 10/19/2012) 6 MUST DO Questions, 1 OPTIONAL question 1. Below, several phase portraits are shown.

### Oscillations. Phys101 Lectures 28, 29. Key points: Simple Harmonic Motion (SHM) SHM Related to Uniform Circular Motion The Simple Pendulum

Phys101 Lectures 8, 9 Oscillations Key points: Simple Harmonic Motion (SHM) SHM Related to Uniform Circular Motion The Simple Pendulum Ref: 11-1,,3,4. Page 1 Oscillations of a Spring If an object oscillates

### Continuous-Time Dynamical Systems: Sketching solution curves and bifurcations

Math 360 Winter 07 Section 0 HW#4 Continuous-Time Dynamical Systems: Sketching solution curves and bifurcations (Due in one week! Thursday Nov 9, 07) Materials contained in this assignment are of great

### 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

### Nonlinear Control Lecture 2:Phase Plane Analysis

Nonlinear Control Lecture 2:Phase Plane Analysis Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Fall 2010 r. Farzaneh Abdollahi Nonlinear Control Lecture 2 1/53

### Chapter #4 EEE8086-EEE8115. Robust and Adaptive Control Systems

Chapter #4 Robust and Adaptive Control Systems Nonlinear Dynamics.... Linear Combination.... Equilibrium points... 3 3. Linearisation... 5 4. Limit cycles... 3 5. Bifurcations... 4 6. Stability... 6 7.

### 2 One-dimensional models in discrete time

2 One-dimensional models in discrete time So far, we have assumed that demographic events happen continuously over time and can thus be written as rates. For many biological species with overlapping generations

### College of Natural Science Department of Mathematics

College of Natural Science Department of Mathematics Graduate Project Report on Difference Equations and Bifurcation Analysis of Discrete Dynamical System Submitted in partial fulfilment of the requirements

### 6.3. Nonlinear Systems of Equations

G. NAGY ODE November,.. Nonlinear Systems of Equations Section Objective(s): Part One: Two-Dimensional Nonlinear Systems. ritical Points and Linearization. The Hartman-Grobman Theorem. Part Two: ompeting

### Dynamical Systems & Lyapunov Stability

Dynamical Systems & Lyapunov Stability Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University Outline Ordinary Differential Equations Existence & uniqueness Continuous dependence

### 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 ẋ

### An Undergraduate s Guide to the Hartman-Grobman and Poincaré-Bendixon Theorems

An Undergraduate s Guide to the Hartman-Grobman and Poincaré-Bendixon Theorems Scott Zimmerman MATH181HM: Dynamical Systems Spring 2008 1 Introduction The Hartman-Grobman and Poincaré-Bendixon Theorems

### CDS 101 Precourse Phase Plane Analysis and Stability

CDS 101 Precourse Phase Plane Analysis and Stability Melvin Leok Control and Dynamical Systems California Institute of Technology Pasadena, CA, 26 September, 2002. mleok@cds.caltech.edu http://www.cds.caltech.edu/

### PHYSICS 110A : CLASSICAL MECHANICS

PHYSICS 110A : CLASSICAL MECHANICS 1. Introduction to Dynamics motion of a mechanical system equations of motion : Newton s second law ordinary differential equations (ODEs) dynamical systems simple 2.

### Stability lectures. Stability of Linear Systems. Stability of Linear Systems. Stability of Continuous Systems. EECE 571M/491M, Spring 2008 Lecture 5

EECE 571M/491M, Spring 2008 Lecture 5 Stability of Continuous Systems http://courses.ece.ubc.ca/491m moishi@ece.ubc.ca Dr. Meeko Oishi Electrical and Computer Engineering University of British Columbia,

### 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

### Recurrent Gait of Anthropomorphic, Bipedal Walkers

Recurrent Gait of Anthropomorphic, Bipedal Walkers by Colleen E. Shannon Thesis submitted to the faculty of the Virginia Polytechnic Institute and State University in partial fulfillment of the requirements

### Nonlinear FEM. Critical Points. NFEM Ch 5 Slide 1

5 Critical Points NFEM Ch 5 Slide Assumptions for this Chapter System is conservative: total residual is the gradient of a total potential energy function r(u,λ) = (u,λ) u Consequence: the tangent stiffness

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

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

### CDS 101/110a: Lecture 2.1 Dynamic Behavior

CDS 11/11a: Lecture.1 Dynamic Behavior Richard M. Murray 6 October 8 Goals: Learn to use phase portraits to visualize behavior of dynamical systems Understand different types of stability for an equilibrium

### Autonomous systems. Ordinary differential equations which do not contain the independent variable explicitly are said to be autonomous.

Autonomous equations Autonomous systems Ordinary differential equations which do not contain the independent variable explicitly are said to be autonomous. i f i(x 1, x 2,..., x n ) for i 1,..., n As you

### E209A: Analysis and Control of Nonlinear Systems Problem Set 3 Solutions

E09A: Analysis and Control of Nonlinear Systems Problem Set 3 Solutions Michael Vitus Stanford University Winter 007 Problem : Planar phase portraits. Part a Figure : Problem a This phase portrait is correct.

### 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

### MATH 415, WEEK 11: Bifurcations in Multiple Dimensions, Hopf Bifurcation

MATH 415, WEEK 11: Bifurcations in Multiple Dimensions, Hopf Bifurcation 1 Bifurcations in Multiple Dimensions When we were considering one-dimensional systems, we saw that subtle changes in parameter

### 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

### Problem set 7 Math 207A, Fall 2011 Solutions

Problem set 7 Math 207A, Fall 2011 s 1. Classify the equilibrium (x, y) = (0, 0) of the system x t = x, y t = y + x 2. Is the equilibrium hyperbolic? Find an equation for the trajectories in (x, y)- phase

### CDS 101/110a: Lecture 2.1 Dynamic Behavior

CDS 11/11a: Lecture 2.1 Dynamic Behavior Richard M. Murray 6 October 28 Goals: Learn to use phase portraits to visualize behavior of dynamical systems Understand different types of stability for an equilibrium

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

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

### Math 345 Intro to Math Biology Lecture 7: Models of System of Nonlinear Difference Equations

Math 345 Intro to Math Biology Lecture 7: Models of System of Nonlinear Difference Equations Junping Shi College of William and Mary, USA Equilibrium Model: x n+1 = f (x n ), here f is a nonlinear function

### Part II. Dynamical Systems. Year

Part II Year 2017 2016 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2017 34 Paper 1, Section II 30A Consider the dynamical system where β > 1 is a constant. ẋ = x + x 3 + βxy 2, ẏ = y + βx 2

### of the Schnakenberg model

Pulse motion in the semi-strong limit of the Schnakenberg model Newton Institute 2005 Jens Rademacher, Weierstraß Institut Berlin joint work with Michael Ward (UBC) Angelfish 2, 6, 12 months old [Kondo,

### 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

### MAT 22B - Lecture Notes

MAT 22B - Lecture Notes 4 September 205 Solving Systems of ODE Last time we talked a bit about how systems of ODE arise and why they are nice for visualization. Now we'll talk about the basics of how to

### Predictability: Does the Flap of a Butterfly s Wings in Brazil set off a Tornado

Chapter 4 Dynamic Behavior Predictability: Does the Flap of a Butterfly s Wings in Brazil set off a Tornado in Texas? Talk given by Edward Lorenz, December 972 meeting of the American Association for the

### Hello everyone, Best, Josh

Hello everyone, As promised, the chart mentioned in class about what kind of critical points you get with different types of eigenvalues are included on the following pages (The pages are an ecerpt from

### ENGI Duffing s Equation Page 4.65

ENGI 940 4. - Duffing s Equation Page 4.65 4. Duffing s Equation Among the simplest models of damped non-linear forced oscillations of a mechanical or electrical system with a cubic stiffness term is Duffing

### arxiv: v2 [physics.optics] 15 Dec 2016

Bifurcation analysis of the Yamada model for a pulsing semiconductor laser with saturable absorber and delayed optical feedback Abstract Soizic Terrien 1, Bernd Krauskopf 1, Neil G.R. Broderick 2 arxiv:1610.06794v2

### Problem Set 5 Solutions

Problem Set 5 Solutions Dorian Abbot APM47 0/8/04. (a) We consider π θ π. The pendulum pointin down corresponds to θ=0 and the pendulum pointin up corresponds to θ π. Define ν θ. The system can be rewritten: