ME 680 Spring Geometrical Analysis of 1D Dynamical Systems

 Merry Walters
 5 months ago
 Views:
Transcription
1 ME 680 Spring 2014 Geometrical Analysis of 1D Dynamical Systems 1
2 Geometrical Analysis of 1D 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 1D 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 1D 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 (HartmanGrobman 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 nonhyperbolic. According to the HartmanGrobman 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 1D 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 Saddlenode 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., HartmanGrobman 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 1D systems undergoing a saddlenode 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 supercritical saddle node bifurcation 20
21 f() = a + b 2 is the normal form of saddle  node bifurcation, i.e., all systems in 1D 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 selforganization 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 nonzero 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 nonlinear 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 saddlenode, 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 twoparameter 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 3D we can visualize the solutions set as follows: h h=0 catastrophic surface r=0 r 38
39 1D system on a circle: [overdamped 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 nonuniformly 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 saddlenode 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 onedimensional phase spaces, and in particular their bifurcations. All equations in this Chapter are scalar
More information8.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 integrodifferential equations
More information4 SecondOrder Systems
4 SecondOrder Systems Secondorder 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
More informationChapter 24 BIFURCATIONS
Chapter 24 BIFURCATIONS Abstract Keywords: Phase Portrait Fixed Point SaddleNode Bifurcation Diagram Codimension1 Hysteresis Hopf Bifurcation SNIC Page 1 24.1 Introduction In linear systems, responses
More informationAPPPHYS217 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 (SpringerVerlag
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 informationLecture 5. Outline: Limit Cycles. Definition and examples How to rule out limit cycles. PoincareBendixson 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 PoincareBendixson theorem Hopf bifurcations Poincare maps Limit
More informationNonlinear 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
More information1. < 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
More informationPHY411 Lecture notes Part 4
PHY411 Lecture notes Part 4 Alice Quillen February 1, 2016 Contents 0.1 Introduction.................................... 2 1 Bifurcations of onedimensional dynamical systems 2 1.1 Saddlenode bifurcation.............................
More informationCollege 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
More informationNonlinear 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
More informationAn Undergraduate s Guide to the HartmanGrobman and PoincaréBendixon Theorems
An Undergraduate s Guide to the HartmanGrobman and PoincaréBendixon Theorems Scott Zimmerman MATH181HM: Dynamical Systems Spring 2008 1 Introduction The HartmanGrobman and PoincaréBendixon Theorems
More informationPhysics 235 Chapter 4. Chapter 4 NonLinear Oscillations and Chaos
Chapter 4 NonLinear Oscillations and Chaos NonLinear Differential Equations Up to now we have considered differential equations with terms that are proportional to the acceleration, the velocity, and
More informationDynamical 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
More informationPart 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
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 informationStability 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,
More informationNonlinear 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
More informationENGI 9420 Lecture Notes 4  Stability Analysis Page Stability Analysis for Nonlinear Ordinary Differential Equations
ENGI 940 Lecture Notes 4  Stability Analysis Page 4.01 4. Stability Analysis for Nonlinear Ordinary Differential Equations A pair of simultaneous first order homogeneous linear ordinary differential
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 informationModelling in Biology
Modelling in Biology Dr GuyBart Stan Department of Bioengineering 17th October 2017 Dr GuyBart Stan (Dept. of Bioeng.) Modelling in Biology 17th October 2017 1 / 77 1 Introduction 2 Linear models of
More informationChapter 9 Global Nonlinear Techniques
Chapter 9 Global Nonlinear Techniques Consider nonlinear dynamical system 0 Nullcline X 0 = F (X) = B @ f 1 (X) f 2 (X). f n (X) x j nullcline = fx : f j (X) = 0g equilibrium solutions = intersection of
More information6.2 Brief review of fundamental concepts about chaotic systems
6.2 Brief review of fundamental concepts about chaotic systems Lorenz (1963) introduced a 3variable model that is a prototypical example of chaos theory. These equations were derived as a simplification
More informationarxiv: v1 [nlin.cd] 16 May 2017
A Tube Dynamics Perspective Governing Stability Transitions: An Example Based on Snapthrough Buckling Jun Zhong 1,, Lawrence N. Virgin,, Shane D. Ros, arxiv:175.593v1 [nlin.cd] 16 May 17 Abstract The
More informationChapter 1 Bifurcations and Chaos in Dynamical Systems
Chapter 1 Bifurcations and Chaos in Dynamical Systems Complex system theory deals with dynamical systems containing often a large number of variables. It extends dynamical system theory, which deals with
More informationDynamics of a massspringpendulum system with vastly different frequencies
Dynamics of a massspringpendulum system with vastly different frequencies Hiba Sheheitli, hs497@cornell.edu Richard H. Rand, rhr2@cornell.edu Cornell University, Ithaca, NY, USA Abstract. We investigate
More informationDIFFERENTIAL GEOMETRY APPLIED TO DYNAMICAL SYSTEMS
WORLD SCIENTIFIC SERIES ON NONLINEAR SCIENCE Series Editor: Leon O. Chua Series A Vol. 66 DIFFERENTIAL GEOMETRY APPLIED TO DYNAMICAL SYSTEMS JeanMarc Ginoux Université du Sud, France World Scientific
More informationMECH 401 Mechanical Design Applications
MECH 40 Mechanical Design Applications Dr. M. K. O Malley Master Notes Spring 008 Dr. D. M. McStravick Rice University Design Considerations Stress Deflection Strain Stiffness Stability Often the controlling
More information7 Planar systems of linear ODE
7 Planar systems of linear ODE Here I restrict my attention to a very special class of autonomous ODE: linear ODE with constant coefficients This is arguably the only class of ODE for which explicit solution
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 informationDynamics of Modified LeslieGower PredatorPrey Model with Predator Harvesting
International Journal of Basic & Applied Sciences IJBASIJENS Vol:13 No:05 55 Dynamics of Modified LeslieGower PredatorPrey Model with Predator Harvesting K. Saleh Department of Mathematics, King Fahd
More information7 Pendulum. Part II: More complicated situations
MATH 35, by T. Lakoba, University of Vermont 60 7 Pendulum. Part II: More complicated situations In this Lecture, we will pursue two main goals. First, we will take a glimpse at a method of Classical Mechanics
More informationNonlinear systems, chaos and control in Engineering
Nonlinear systems, chaos and control in Engineering Module 1 block 3 Onedimensional nonlinear systems Cristina Masoller Cristina.masoller@upc.edu http://www.fisica.edu.uy/~cris/ Schedule Flows on the
More informationMath 273 (51)  Final
Name: Id #: Math 273 (5)  Final Autumn Quarter 26 Thursday, December 8, 266: to 8: Instructions: Prob. Points Score possible 25 2 25 3 25 TOTAL 75 Read each problem carefully. Write legibly. Show all
More informationKinematics of fluid motion
Chapter 4 Kinematics of fluid motion 4.1 Elementary flow patterns Recall the discussion of flow patterns in Chapter 1. The equations for particle paths in a threedimensional, steady fluid flow are dx
More informationUniversity of Utrecht
University of Utrecht Bachelor thesis The bifurcation diagram of the second nontrivial normal form of an axially symmetric perturbation of the isotropic harmonic oscillator Author: T. Welker Supervisor:
More informationTangent spaces, normals and extrema
Chapter 3 Tangent spaces, normals and extrema If S is a surface in 3space, with a point a S where S looks smooth, i.e., without any fold or cusp or selfcrossing, we can intuitively define the tangent
More informationA ProducerConsumer Model With Stoichiometry
A ProducerConsumer Model With Stoichiometry Plan B project toward the completion of the Master of Science degree in Mathematics at University of Minnesota Duluth Respectfully submitted by Laura Joan Zimmermann
More informationProblem Set Number 1, j/2.036j MIT (Fall 2014)
Problem Set Number 1, 18.385j/2.036j MIT (Fall 2014) Rodolfo R. Rosales (MIT, Math. Dept., Cambridge, MA 02139) September 9, 2014 Due Fri., September 19, 2014. 1 Get equation from phase line portrait problem
More informationANSWERS Final Exam Math 250b, Section 2 (Professor J. M. Cushing), 15 May 2008 PART 1
ANSWERS Final Exam Math 50b, Section (Professor J. M. Cushing), 5 May 008 PART. (0 points) A bacterial population x grows exponentially according to the equation x 0 = rx, where r>0is the per unit rate
More informationIntroduction to the Phase Plane
Introduction to the Phase Plane June, 06 The Phase Line A single first order differential equation of the form = f(y) () makes no mention of t in the function f. Such a differential equation is called
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 informationDynamics of Decision Making in Animal Group. Motion
Dynamics of Decision Making in Animal Group Motion Benjamin Nabet, Naomi E Leonard, Iain D Couzin and Simon A Levin June, 007 Abstract We present a continuous model of a multiagent system motivated by
More informationTravelling waves. Chapter 8. 1 Introduction
Chapter 8 Travelling waves 1 Introduction One of the cornerstones in the study of both linear and nonlinear PDEs is the wave propagation. A wave is a recognizable signal which is transferred from one part
More informationMath 312 Lecture Notes Linear Twodimensional Systems of Differential Equations
Math 2 Lecture Notes Linear Twodimensional Systems of Differential Equations Warren Weckesser Department of Mathematics Colgate University February 2005 In these notes, we consider the linear system of
More information2DVolterraLotka Modeling For 2 Species
Majalat AlUlum AlInsaniya wat  Tatbiqiya 2DVolterraLotka Modeling For 2 Species Alhashmi Darah 1 University of Almergeb Department of Mathematics Faculty of Science Zliten Libya. Abstract The purpose
More informationTopic # /31 Feedback Control Systems. Analysis of Nonlinear Systems Lyapunov Stability Analysis
Topic # 16.30/31 Feedback Control Systems Analysis of Nonlinear Systems Lyapunov Stability Analysis Fall 010 16.30/31 Lyapunov Stability Analysis Very general method to prove (or disprove) stability of
More informationInClass Problems 3032: Moment of Inertia, Torque, and Pendulum: Solutions
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Physics Physics 8.01 TEAL Fall Term 004 InClass Problems 303: Moment of Inertia, Torque, and Pendulum: Solutions Problem 30 Moment of Inertia of a
More informationCHAPTER 5 KINEMATICS OF FLUID MOTION
CHAPTER 5 KINEMATICS OF FLUID MOTION 5. ELEMENTARY FLOW PATTERNS Recall the discussion of flow patterns in Chapter. The equations for particle paths in a threedimensional, steady fluid flow are dx 
More information2 Qualitative theory of nonsmooth dynamical systems
2 Qualitative theory of nonsmooth dynamical systems In this chapter, we give an overview of the basic theory of both smooth and nonsmooth dynamical systems, to be expanded upon in later chapters. In
More informationRigidBody Attitude Control USING ROTATION MATRICES FOR CONTINUOUS, SINGULARITYFREE CONTROL LAWS
RigidBody Attitude Control USING ROTATION MATRICES FOR CONTINUOUS, SINGULARITYFREE CONTROL LAWS 3 IEEE CONTROL SYSTEMS MAGAZINE» JUNE 33X//$. IEEE SHANNON MASH Rigidbody attitude control is motivated
More informationNonlinear Control Lecture 4: Stability Analysis I
Nonlinear Control Lecture 4: Stability Analysis I Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Fall 2010 Farzaneh Abdollahi Nonlinear Control Lecture 4 1/70
More informationSymmetry Properties of Confined Convective States
Symmetry Properties of Confined Convective States John Guckenheimer Cornell University 1 Introduction This paper is a commentary on the experimental observation observations of Bensimon et al. [1] of convection
More information2:1 Resonance in the delayed nonlinear Mathieu equation
Nonlinear Dyn 007) 50:31 35 DOI 10.1007/s110710069165 ORIGINAL ARTICLE :1 Resonance in the delayed nonlinear Mathieu equation Tina M. Morrison Richard H. Rand Received: 9 March 006 / Accepted: 19 September
More informationChapter 8 Equilibria in Nonlinear Systems
Chapter 8 Equilibria in Nonlinear Sstems Recall linearization for Nonlinear dnamical sstems in R n : X 0 = F (X) : if X 0 is an equilibrium, i.e., F (X 0 ) = 0; then its linearization is U 0 = AU; 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 informationSymmetries 2  Rotations in Space
Symmetries 2  Rotations in Space This symmetry is about the isotropy of space, i.e. space is the same in all orientations. Thus, if we continuously rotated an entire system in space, we expect the system
More informationUsing Lyapunov Theory I
Lecture 33 Stability Analysis of Nonlinear Systems Using Lyapunov heory I Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science  Bangalore Outline Motivation Definitions
More informationOscillations. Simple Harmonic Motion of a Mass on a Spring The equation of motion for a mass m is attached to a spring of constant k is
Dr. Alain Brizard College Physics I (PY 10) Oscillations Textbook Reference: Chapter 14 sections 18. Simple Harmonic Motion of a Mass on a Spring The equation of motion for a mass m is attached to a spring
More informationRepresent this system in terms of a block diagram consisting only of. g From Newton s law: 2 : θ sin θ 9 θ ` T
Exercise (Block diagram decomposition). Consider a system P that maps each input to the solutions of 9 4 ` 3 9 Represent this system in terms of a block diagram consisting only of integrator systems, represented
More informationQuestion 1: Spherical Pendulum
Question 1: Spherical Pendulum Consider a twodimensional pendulum of length l with mass M at its end. It is easiest to use spherical coordinates centered at the pivot since the magnitude of the position
More informationDef. (a, b) is a critical point of the autonomous system. 1 Proper node (stable or unstable) 2 Improper node (stable or unstable)
Types of critical points Def. (a, b) is a critical point of the autonomous system Math 216 Differential Equations Kenneth Harris kaharri@umich.edu Department of Mathematics University of Michigan November
More informationMultiple scale methods
Multiple scale methods G. Pedersen MEK 3100/4100, Spring 2006 March 13, 2006 1 Background Many physical problems involve more than one temporal or spatial scale. One important example is the boundary layer
More informationTHE GEOMETRY OF SLOW MANIFOLDS NEAR A FOLDED NODE
THE GEOMETRY OF SLOW MANIFOLDS NEAR A FOLDED NODE M DESROCHES, B KRAUSKOPF, HM OSINGA Abstract This paper is concerned with the geometry of slow manifolds of a dynamical system with two slow and one fast
More informationFor a rigid body that is constrained to rotate about a fixed axis, the gravitational torque about the axis is
Experiment 14 The Physical Pendulum The period of oscillation of a physical pendulum is found to a high degree of accuracy by two methods: theory and experiment. The values are then compared. Theory For
More informationMathematical Models of Biological Systems
Mathematical Models of Biological Systems CH924 November 25, 2012 Course outline Syllabus: Week 6: Stability of firstorder autonomous differential equation, genetic switch/clock, introduction to enzyme
More informationHarmonic Oscillator. MassSpring Oscillator Resonance The Pendulum. Physics 109 Experiment Number 12
Harmonic Oscillator MassSpring Oscillator Resonance The Pendulum Physics 109 Experiment Number 12 Outline Simple harmonic motion The vertical massspring system Driven oscillations and resonance The pendulum
More informationExact model reduction by a slow fast decomposition of nonlinear mechanical systems
Nonlinear Dyn 07 90:67 647 DOI 0.007/s070736859 ORIGINAL PAPER Exact model reduction by a slow fast decomposition of nonlinear mechanical systems George Haller Sten Ponsioen Received: 7 November 06
More informationPhysics: springmass system, planet motion, pendulum. Biology: ecology problem, neural conduction, epidemics
Applications of nonlinear ODE systems: Physics: springmass system, planet motion, pendulum Chemistry: mixing problems, chemical reactions Biology: ecology problem, neural conduction, epidemics Economy:
More informationSimple Iteration, cont d
Jim Lambers MAT 772 Fall Semester 201011 Lecture 2 Notes These notes correspond to Section 1.2 in the text. Simple Iteration, cont d In general, nonlinear equations cannot be solved in a finite sequence
More information1 / 5. Stability of 2 2 Systems of Differential Equations April 2014
/ 5 Stability of Systems of Differential Equations April 04 Alecos Papadopoulos Department of Economics Athens University of Economics and Business papadopalex@aueb.gr http://alecospapadopoulos.wordpress.com
More informationRotational Kinetic Energy
Lecture 17, Chapter 10: Rotational Energy and Angular Momentum 1 Rotational Kinetic Energy Consider a rigid body rotating with an angular velocity ω about an axis. Clearly every point in the rigid body
More informationt = g = 10 m/s 2 = 2 s T = 2π g
Annotated Answers to the 1984 AP Physics C Mechanics Multiple Choice 1. D. Torque is the rotational analogue of force; F net = ma corresponds to τ net = Iα. 2. C. The horizontal speed does not affect the
More informationTorsion Spring Oscillator with Dry Friction
Torsion Spring Oscillator with Dry Friction Manual Eugene Butikov Annotation. The manual includes a description of the simulated physical system and a summary of the relevant theoretical material for students
More informationNonlinear Marangoni convection in circular. and elliptical cylinders
Nonlinear rangoni convection in circular and elliptical cylinders P. Assemat 1, A. Bergeon 1 & E. Knobloch 2 1 IMFT UMR CNRS 5502  UPS UFR MIG, 31062 Toulouse Cedex, France 2 Department of Physics, University
More informationPiecewise Smooth Solutions to the BurgersHilbert Equation
Piecewise Smooth Solutions to the BurgersHilbert Equation Alberto Bressan and Tianyou Zhang Department of Mathematics, Penn State University, University Park, Pa 68, USA emails: bressan@mathpsuedu, zhang
More informationReport EProject Henriette Laabsch Toni Luhdo Steffen Mitzscherling Jens Paasche Thomas Pache
Potsdam, August 006 Report EProject Henriette Laabsch 7685 Toni Luhdo 7589 Steffen Mitzscherling 7540 Jens Paasche 7575 Thomas Pache 754 Introduction From 7 th February till 3 rd March, we had our laboratory
More information28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod)
28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod) θ + ω 2 sin θ = 0. Indicate the stable equilibrium points as well as the unstable equilibrium points.
More informationStudies in Applied Economics
SAE./No.32/April 2015 Studies in Applied Economics Bifurcation of MacroeconoMetric Models and robustness of dynamical inferences William A. Barnett and Guo Chen Johns Hopkins Institute for Applied Economics,
More informationStability of flow past a confined cylinder
Stability of flow past a confined cylinder Andrew Cliffe and Simon Tavener University of Nottingham and Colorado State University Stability of flow past a confined cylinder p. 1/60 Flow past a cylinder
More informationStructural model for the dynamic buckling of a column under constant rate compression
Structural model for the dynamic buckling of a column under constant rate compression Vitaly A. Kuzkin 1 Institute for Problems in Mechanical Engineering RAS Saint Petersburg Polytechnical University arxiv:1506.0047v
More informationInputtostate stability and interconnected Systems
10th Elgersburg School Day 1 Inputtostate stability and interconnected Systems Sergey Dashkovskiy Universität Würzburg Elgersburg, March 5, 2018 1/20 Introduction Consider Solution: ẋ := dx dt = ax,
More informationThermal Analysis. with SolidWorks Simulation 2013 SDC. Paul M. Kurowski. Better Textbooks. Lower Prices.
Thermal Analysis with SolidWorks Simulation 2013 Paul M. Kurowski SDC PUBLICATIONS Schroff Development Corporation Better Textbooks. Lower Prices. www.sdcpublications.com Visit the following websites to
More informationLecture 5: Rules of Differentiation. First Order Derivatives
Lecture 5: Rules of Differentiation First order derivatives Higher order derivatives Partial differentiation Higher order partials Differentials Derivatives of implicit functions Generalized implicit function
More informationLecture Notes. W. Tang. School of Mathematical & Statistical Sciences Arizona State University
Applied Math / Dynamical Systems Workshop Lecture Notes W. Tang School of Mathematical & Statistical Sciences Arizona State University May. 23 rd 25 th, 2011 Tang (ASU) AMDS Workshop Notes May. 23 rd
More informationTwoBody Problem. Central Potential. 1D Motion
TwoBody Problem. Central Potential. D Motion The simplest nontrivial dynamical problem is the problem of two particles. The equations of motion read. m r = F 2, () We already know that the center of
More informationCharacterization of the stability boundary of nonlinear autonomous dynamical systems in the presence of a saddlenode equilibrium point of type 0
Anais do CNMAC v.2 ISSN 198482X Characterization of the stability boundary of nonlinear autonomous dynamical systems in the presence of a saddlenode equilibrium point of type Fabíolo M. Amaral Departamento
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 informationSolving Differential Equations: First Steps
30 ORDINARY DIFFERENTIAL EQUATIONS 3 Solving Differential Equations Solving Differential Equations: First Steps Now we start answering the question which is the theme of this book given a differential
More informationL(q, q) = m 2 q2 V (q) 2 m + V (q)
Lecture 7 Phase Space, Part 1 MATHGA 71.1 Mechanics 1 Phase portraits 1.1 One dimensional system Consider the generic one dimensional case of a point mass m described by a generalized coordinate q and
More informationAnswer Key 1973 BC 1969 BC 24. A 14. A 24. C 25. A 26. C 27. C 28. D 29. C 30. D 31. C 13. C 12. D 12. E 3. A 32. B 27. E 34. C 14. D 25. B 26.
Answer Key 969 BC 97 BC. C. E. B. D 5. E 6. B 7. D 8. C 9. D. A. B. E. C. D 5. B 6. B 7. B 8. E 9. C. A. B. E. D. C 5. A 6. C 7. C 8. D 9. C. D. C. B. A. D 5. A 6. B 7. D 8. A 9. D. E. D. B. E. E 5. E.
More informationIn these chapter 2A notes write vectors in boldface to reduce the ambiguity of the notation.
1 2 Linear Systems In these chapter 2A notes write vectors in boldface to reduce the ambiguity of the notation 21 Matrix ODEs Let and is a scalar A linear function satisfies Linear superposition ) Linear
More informationT = 1 2 mv2. For a system of N particles, the kinetic energy of the system is obtained by adding up the individual kinetic energies: T = 2 m ivi 2.
Physics 3550, Fall 2011 Energy. Relevant Sections in Text: Chapter 4 Kinetic Energy. I probably don t need to impress upon you the importance of the notion of energy in analyzing/understanding physical
More informationNonlinear System Analysis
Nonlinear System Analysis Lyapunov Based Approach Lecture 4 Module 1 Dr. Laxmidhar Behera Department of Electrical Engineering, Indian Institute of Technology, Kanpur. January 4, 2003 Intelligent Control
More informationterm from the numerator yields 2
APPM 1350 Eam 2 Fall 2013 1. The following parts are not related: (a) (12 pts) Find y given: (i) y = (ii) y = sec( 2 1) tan() (iii) ( 2 + y 2 ) 2 = 2 2 2y 2 1 (b) (8 pts) Let f() be a function such that
More informationHoney Coiling  A Study on the Gravitational Regime of Liquid Rope Coiling
Honey Coiling  A Study on the Gravitational Regime of Liquid Rope Coiling Patrick Meister, MNG Rämibühl, patrickmeister@ymail.com 1 Introduction We report on the coiling motion a falling stream of viscous
More informationDownloaded from
Question 4.1: A circular coil of wire consisting of 100 turns, each of radius 8.0 cm carries a current of 0.40 A. What is the magnitude of the magnetic field B at the centre of the coil? Number of turns
More informationChapter 4. Transition towards chaos. 4.1 Onedimensional maps
Chapter 4 Transition towards chaos In this chapter we will study how successive bifurcations can lead to chaos when a parameter is tuned. It is not an extensive review : there exists a lot of different
More informationNumerical Methods. Root Finding
Numerical Methods Solving Non Linear 1Dimensional Equations Root Finding Given a real valued function f of one variable (say ), the idea is to find an such that: f() 0 1 Root Finding Eamples Find real
More information