ENGI 9420 Lecture Notes 4  Stability Analysis Page Stability Analysis for Nonlinear Ordinary Differential Equations


 Phyllis McKinney
 1 years ago
 Views:
Transcription
1 ENGI 940 Lecture Notes 4  Stability Analysis Page Stability Analysis for Nonlinear Ordinary Differential Equations A pair of simultaneous first order homogeneous linear ordinary differential equations for two functions x(t), y(t) of one independent variable t, x = = ax + by dy y = = cx + d y x a b x may be represented by the matrix equation = y c d y A single second order linear homogeneous ordinary differential equation for x(t) with constant coefficients, d x dy + p + qx = 0 + py + qx = 0 may be rewritten as a linked pair of first order homogeneous ordinary differential equations, by introducing a second dependent variable: dy d d x = y = = dy = qx py x 0 1 x and may also be represented in matrix form = y q p y The general solution for (x, y) in either case can be displayed graphically as a set of contour curves (or level curves) in a phase space. Sections in this chapter: 4.01 Motion of a Pendulum 4.0 Stability of Stationary Points 4.03 Linear Approximation to a System of NonLinear ODEs (1) 4.04 Reminder of Linear Ordinary Differential Equations 4.05 Stability Analysis for a Linear System 4.06 Linear Approximation to a System of NonLinear ODEs () 4.07 Limit Cycles 4.08 Van der Pol s Equation 4.09 Theorem for Limit Cycles 4.10 Lyapunov Functions [for reference only  not examinable] 4.11 Duffing s Equation 4.1 More Examples 4.13 Liénard s Theorem
2 ENGI Pendulum Page Motion of a Pendulum Consider a pendulum, moving under its own weight, without friction. The pendulum bob has mass m, the shaft has length L and negligible mass, and the angle of the shaft with the vertical is x. The tension along the shaft is T. The acceleration due to gravity is g ( 9.81 m s ). Resolving forces radially (centripetal force) Resolving forces transverse to the pendulum: T + mg cos x = ml x mg sin x = ml x g (1) L x + k sin x = 0, where k = The Maclaurin series expansion of sin x is: 3 5 x x sin x = x + 3! 5! Provided the oscillations of the pendulum are small, (x << 1), sin x x and the ordinary differential equation governing the motion of the pendulum is, to a good approximation, d x + kx= 0 () (which is the ODE of simple harmonic motion)
3 ENGI Pendulum Page 4.03 Let the angular velocity of the pendulum be v Then, using the chain rule of differentiation, dv dv dv x = v = = = v = x =. The ODE () becomes dv v k x 0 v dv k x 0 + = + = v x c v dv + k x = 0 + k = v + k x = c (3) If at time t = 0 the pendulum is passing through its equilibrium position with angular speed v o, then the initial conditions are x v v ( 0) = 0, = ( 0) = o t = 0 Substituting the initial conditions into (), vo + k 0 = c, which leads to a complete solution for v as an implicit function of x, v + kx = v o (4) A plot of this solution for various choices of v o generates a family of concentric ellipses. Recall that this solution is valid only for small displacements x. The xv plane is called the phase plane.
4 ENGI Pendulum Page 4.04 Returning to the more general case g x + k sin x = 0, where k = (1) L and again using dv dv dv x = v = = = v the ODE can be rewritten as dv v k sin x 0 v dv k sin x 0 + = + = v c v dv + k sin x = 0 k cos x = 1 1 ml v + ( ml k cos x) = ml c (5) KE.. PE.. Total energy However, the kinetic energy is 1 ( ) m Lv. In the absence of friction, the sum of kinetic and potential energy is constant, so that the potential energy of the pendulum must be ml k cos x (= mg L cos x, which makes sense upon examining the diagram on page 4.0). Each value of total energy E = 1 ml c generates an orbit (or energy curve). The relationship between total energy and initial angular velocity is obtained from substituting the initial conditions (x = 0 and v = v o when t = 0) into (5): vo c k = c = vo k v k cos x = vo k (6) vo v = k k cos x = k ( 1 cos x) 0 vo v 0 vo v v vo The maximum angular speed of the pendulum, v o, occurs when x = 0. Also, using cos 1 sin sin 1 cos θ = θ θ = θ, v o v x = 1 cos x = sin k vo v x = sin 4k
5 ENGI Pendulum Page 4.05 Recall that the angular velocity is just v =. Differentiating the complete solutions (6): v k cosx = vo k implicitly with respect to time, we obtain dv dv v + k sin x = 0 = k sin x This expression can also be derived directly from the ODE x + k sin x = 0 (1). d x dv v 0 and k sinx 0 When 0 < x < π and v > 0, = > = = <. Therefore, in the phase plane, as x increases from the starting point (0, v o ), v decreases in the first quadrant, until the maximum value of x (label that maximum value of x as M). Tracking back into the second quadrant, before (0, v o ), d x dv π < x < 0 and v > 0 = v > 0 and = = k sinx > 0. The orbit increases from x = M to a maximum at (0, v o ), then decreases until x = +M. This tracks the motion of the pendulum on its complete swing from left to right. By symmetry, the swing in the opposite direction should generate a mirror image in the x axis of the phase plane, to complete the orbit.
6 ENGI Pendulum Page 4.06 (6) v k cosx = v k o x v v k ( 1 cosx) v k sin = o = o x v = vo ksin Three cases arise: (7) v o < k : v will decrease to zero: x x vo o v= 0 0 = v ksin sin = ± k In the first quadrant of the phase plane, the orbit will move right and down to an intercept M vo on the x axis at (M, 0), where sin = + and 0 < M < π. Extending to the k other three quadrants, the orbits resemble ellipses, centred on the origin. v o = k : v will just barely decrease to zero: x x v= 0 0 = ( k) ksin = ( k) cos M = ± π The pendulum swings all the way to the upsidedown position and comes to rest there, before either swinging back or continuing on in the same direction. v o > k : The pendulum will never come to rest, reaching a nonzero minimum speed as it passes through the upsidedown position: x v sin o = ± has no real solution for x when v o > k. k If it is swinging anticlockwise, then the orbit stays in the first and second quadrants. If it is swinging clockwise, then the orbit stays in the third and fourth quadrants. These orbits are not closed and extend beyond the range π < x < π. π vmin = vo ksin = vo 4 k We can then generate the full set of orbits in the phase plane for the general pendulum problem.
7 ENGI Pendulum Page 4.07 As time progresses, one moves along an orbit to the right above the x axis, but to the left below the x axis v =. The relationship (7) between angular velocity v and angle x is itself a first order nonlinear ordinary differential equation for x as a function of the time t: o sin x v k vo ksin x = = ± ± = t + C x vo 4k sin For the case of closed orbits ( v o < k), the time to complete one orbit (the period T of the pendulum) can be shown to be π / 4 dθ M vo g T =, where b= sin = and k = k 0 1 b sin θ k L This is a complete elliptic integral of the first kind, which has no analytic solution in terms of finite combinations of algebraic functions, (except for special choices of v o and k). As v o k, the period T diverges to infinity it takes forever for the zero energy pendulum to reach the upsidedown position.
8 ENGI Stability Page Stability of Stationary Points Consider the (generally nonlinear) system of simultaneous first order ordinary differential equations dy = x = P( x, y), = y = Q( x, y) (1) (, ) (, ) dy dy y Q x y Using the chain rule, = = x =. P x y This can be integrated with respect to x to obtain a solution for y as an implicit function of x, provided P(x, y) 0. At points where P(x, y) = 0 but Q(x, y) 0, one may P( x, y) integrate = instead. dy Q x, y ( ) Points on the phase plane where P(x, y) = Q(x, y) = 0 are singular points. A unique slope does not exist at such points. Alternative names for singular points are equilibrium points or stationary points (because both x and y do not [instantaneously] change with time there) or critical points or fixed points. A singular point is stable (and is called an "attractor") if the response to a small disturbance remains small for all time. Stable singular point: all paths starting inside the inner circle stay closer than the outer circle forever. Unstable singular point. (or source )
9 ENGI Stability Page 4.09 Consider the system x = Pxy,, y = Qxy,, x0 = x*, y0 = y*, P0,0 = Q0,0 = 0 () ( ) ( ) ( ) ( ) ( ) ( ) which has a stationary point at the origin. Let x(t; x*), y(t; y*) be the complete solution to this system. The stationary point at the origin is stable if and only if, for every ε > 0 (however small), there exists a δ(ε) such that whenever the point (x*, y*) = (x(0; x*), y(0; y*)) is closer than δ to the origin, the point (x(t; x*), y(t; y*)) remains closer than ε to the origin for all time, or ( x* ) ( y* ) δ x ( t; x* ) y ( t; y* ) + < + < ε t A stationary point is asymptotically stable (also known as a sink ) if it is stable and any disturbance ultimately vanishes: lim x ( t; x* ) + y ( t; y* ) = 0. t
10 ENGI Stability Page 4.10 Here are three types of stationary points with nearby orbits: Stable centre (but not asymptotically stable) Unstable saddle point [all saddle points are unstable] Asymptotically stable focus (or spiral sink)
11 ENGI Linear Approximation (1) Page Linear Approximation to a System of NonLinear ODEs (1) The Taylor series of any function f (x, y) about the point (x 0, y 0 ) is f f f ( x, y) = f ( x0, y0) + ( x x0) + ( y y0) + x ( x, y ) y ( x, y ) ( ) ( )( ) ( ) f x x0 f x x0 y y0 f y y + + x ( 0, 0)! x y ( 0, 0)! y ( 0, 0)! x y x y x y provided that the series converges to f (x, y). 0 + (1) This allows us to create a linear approximation to the nonlinear system x = Pxy,, y = Qxy,, x0 = x*, y0 = y*, P0,0 = Q0,0 = 0. () ( ) ( ) ( ) ( ) ( ) ( ) P P P( x, y) = P( 0,0) + x + y + P1 ( x, y) (3) x y (, ) ( 0,0) ( 0,0) P1 x y where lim = 0, (because P 1 (x, y) is at least second order in x, y) ( xy, ) ( 0,0) x + y and similarly for Q(x, y), so that the system becomes x = ax + by + P x, y 1 ( ) (, ) y = cx + d y + Q x y where a, b, c, d are all constants. In the neighbourhood of the singular point (0, 0), this system can be modelled by the linear system x = ax + by (5) y = cx + d y P P Q Q where a =, b =, c =, d =. x y x y ( 0,0) ( 0,0) ( 0,0) ( 0,0) 1 (4)
12 ENGI ODEs Page Reminder of Linear Ordinary Differential Equations To find the general solution of the homogeneous second order linear ODE d y dy + p + q y = 0, with constant real coefficients p and q, form the auxiliary equation or characteristic equation λ + pλ + q = 0 and evaluate the discriminant D = p ± D p 4q and the roots λ1, λ =. Three cases arise. D > 0 : The characteristic equation has a pair of distinct real roots λ 1, λ. The general solution is y = λ1 Ae x + λx Be. D = 0: The characteristic equation has a pair of equal real roots λ. x The general solution is y = ( Ax + B) e λ. D < 0 : The characteristic equation has a complex conjugate pair of roots λ 1, λ = a ± bj. The general solution is y ax = e ( Acosbx + Bsin bx), where A, B are arbitrary constants of integration.
13 ENGI ODEs Page 4.13 To find the general solution of the system of simultaneous first order linear ODEs = ax + by dy = cx + d y λ, t λt x t y t = αe, βe into the ODE, to obtain substitute the trial solution () () ( ) ( ) e λt αλ = aα e λt + bβ e λt βλe λt = cα e λt + dβ e λt or, in matrix form, a λ b α 0 = c d λ β 0 α = β = 0 is a solution (the trivial solution) for any choice of a, b, c, d and λ. Nontrivial solutions exist when the determinant of the matrix of coefficients is zero: a λ b = 0 ( a λ)( d λ) bc = 0 c d λ ( a d ) λ ( ad bc) λ + + = which is the characteristic equation of the system. The solutions to the characteristic equation are the eigenvalues of the coefficient matrix a b A = c d and, for each eigenvalue λ, a nonzero vector α that satisfies the equation β a λ b α 0 = c d λ β 0 is an eigenvector for that eigenvalue. The general solution to the system of ODEs is a linear combination of the solutions arising from each eigenvalue: ( x() t, y() t ) ( c1 1e λ 1 1 ) 1t c e λt, c e λ1t λ = α + α β + c β e t unless the eigenvalues are equal, in which case the general solution is t t x t, y t c c t e λ λ = α + α, cβ + c β t e ( ) ( () ()) ( 1 1 ) ( 1 1 ) (where, in this case, (α 1, β 1 ) is not necessarily an eigenvector). 0
14 ENGI Stability Analysis (Linear) Page Stability Analysis for a Linear System In the case where (0, 0) is the only critical point of the system = ax + by dy = cx + d y it follows that the characteristic equation λ ( a + d ) λ + ( ad bc) = has only nonzero roots and that det A = ad bc 0. 0 Proof: If λ = 0 then at least one eigenvalue of the coefficient matrix A is zero, from which it follows immediately that a 0 b = det A = ad bc = 0 c d 0 Both roots nonzero ad bc 0. If (0, 0) is the only critical point of the system, then no other choice of (x, y) satisfies both equations ax + by = 0 ( ad bc) x = 0 cx + d y = 0 ( ad bc) y = 0 from which it follows immediately that ad bc = det A 0. If the roots are both nonzero and (x, y) is a critical point of the system, then ax + by = 0 ( ad bc) x = 0 cx + d y = 0 ( ad bc) y = 0 But λ1, λ 0 ad bc 0 (0, 0) is the only solution to this pair of simultaneous linear equations. Therefore (0, 0) is the only critical point of the system if and only if both roots of the characteristic equation are nonzero.
ENGI 9420 Lecture Notes 4  Stability Analysis Page Stability Analysis for Nonlinear Ordinary Differential Equations
ENGI 940 Lecture Notes 4  Stability Analysis Page 4.0 4. Stability Analysis for Nonlinear Ordinary Differential Equations A pair of simultaneous first order homogeneous linear ordinary differential equations
More informationENGI Linear Approximation (2) Page Linear Approximation to a System of NonLinear ODEs (2)
ENGI 940 4.06  Linear Approximation () Page 4. 4.06 Linear Approximation to a System of NonLinear ODEs () From sections 4.0 and 4.0, the nonlinear system dx dy = x = P( x, y), = y = Q( x, y) () with
More informationENGI 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 nonlinear forced oscillations of a mechanical or electrical system with a cubic stiffness term is Duffing
More informationNonlinear 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
More informationSTABILITY. Phase portraits and local stability
MAS271 Methods for differential equations Dr. R. Jain STABILITY Phase portraits and local stability We are interested in system of ordinary differential equations of the form ẋ = f(x, y), ẏ = g(x, y),
More informationConsider 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
More informationConsider an ideal pendulum as shown below. l θ is the angular acceleration θ is the angular velocity
1 Second Order Ordinary Differential Equations 1.1 The harmonic oscillator Consider an ideal pendulum as shown below. θ l Fr mg l θ is the angular acceleration θ is the angular velocity A point mass m
More information154 Chapter 9 Hints, Answers, and Solutions The particular trajectories are highlighted in the phase portraits below.
54 Chapter 9 Hints, Answers, and Solutions 9. The Phase Plane 9.. 4. The particular trajectories are highlighted in the phase portraits below... 3. 4. 9..5. Shown below is one possibility with x(t) and
More informationENGI 9420 Lecture Notes 1  ODEs Page 1.01
ENGI 940 Lecture Notes  ODEs Page.0. Ordinary Differential Equations An equation involving a function of one independent variable and the derivative(s) of that function is an ordinary differential equation
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 informationMath 266: Phase Plane Portrait
Math 266: Phase Plane Portrait Long Jin Purdue, Spring 2018 Review: Phase line for an autonomous equation For a single autonomous equation y = f (y) we used a phase line to illustrate the equilibrium solutions
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 informationA plane autonomous system is a pair of simultaneous firstorder differential equations,
Chapter 11 PhasePlane Techniques 11.1 Plane Autonomous Systems A plane autonomous system is a pair of simultaneous firstorder differential equations, ẋ = f(x, y), ẏ = g(x, y). This system has an equilibrium
More informationStability of Nonlinear Systems An Introduction
Stability of Nonlinear Systems An Introduction Michael Baldea Department of Chemical Engineering The University of Texas at Austin April 3, 2012 The Concept of Stability Consider the generic nonlinear
More informationLECTURE 8: DYNAMICAL SYSTEMS 7
15382 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 informationComplex Dynamic Systems: Qualitative vs Quantitative analysis
Complex Dynamic Systems: Qualitative vs Quantitative analysis Complex Dynamic Systems Chiara Mocenni Department of Information Engineering and Mathematics University of Siena (mocenni@diism.unisi.it) Dynamic
More informationAutonomous Systems and Stability
LECTURE 8 Autonomous Systems and Stability An autonomous system is a system of ordinary differential equations of the form 1 1 ( 1 ) 2 2 ( 1 ). ( 1 ) or, in vector notation, x 0 F (x) That is to say, an
More informationChapter 6 Nonlinear Systems and Phenomena. Friday, November 2, 12
Chapter 6 Nonlinear Systems and Phenomena 6.1 Stability and the Phase Plane We now move to nonlinear systems Begin with the firstorder system for x(t) d dt x = f(x,t), x(0) = x 0 In particular, consider
More informationCopyright (c) 2006 Warren Weckesser
2.2. PLANAR LINEAR SYSTEMS 3 2.2. Planar Linear Systems We consider the linear system of two first order differential equations or equivalently, = ax + by (2.7) dy = cx + dy [ d x x = A x, where x =, and
More informationMatrix Solutions to Linear Systems of ODEs
Matrix Solutions to Linear Systems of ODEs James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University November 3, 216 Outline 1 Symmetric Systems of
More informationAdvanced Eng. Mathematics
Koya University Faculty of Engineering Petroleum Engineering Department Advanced Eng. Mathematics Lecture 6 Prepared by: Haval Hawez Email: haval.hawez@koyauniversity.org 1 Second Order Linear Ordinary
More informationCDS 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/
More informationMath 4B Notes. Written by Victoria Kala SH 6432u Office Hours: T 12:45 1:45pm Last updated 7/24/2016
Math 4B Notes Written by Victoria Kala vtkala@math.ucsb.edu SH 6432u Office Hours: T 2:45 :45pm Last updated 7/24/206 Classification of Differential Equations The order of a differential equation is the
More informationSection 9.3 Phase Plane Portraits (for Planar Systems)
Section 9.3 Phase Plane Portraits (for Planar Systems) Key Terms: Equilibrium point of planer system yꞌ = Ay o Equilibrium solution Exponential solutions o Halfline solutions Unstable solution Stable
More informationPower series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) a n z n. n=0
Lecture 22 Power series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) Recall a few facts about power series: a n z n This series in z is centered at z 0. Here z can
More informationODE. Philippe Rukimbira. Department of Mathematics Florida International University PR (FIU) MAP / 92
ODE Philippe Rukimbira Department of Mathematics Florida International University PR (FIU) MAP 2302 1 / 92 4.4 The method of Variation of parameters 1. Second order differential equations (Normalized,
More informationNonlinear differential equations  phase plane analysis
Nonlinear differential equations  phase plane analysis We consider the general first order differential equation for y(x Revision Q(x, y f(x, y dx P (x, y. ( Curves in the (x, yplane which satisfy this
More informationReview: control, feedback, etc. Today s topic: statespace models of systems; linearization
Plan of the Lecture Review: control, feedback, etc Today s topic: statespace models of systems; linearization Goal: a general framework that encompasses all examples of interest Once we have mastered
More informationExam 4 SCORE. MA 114 Exam 4 Spring Section and/or TA:
Exam 4 Name: Section and/or TA: Last Four Digits of Student ID: Do not remove this answer page you will return the whole exam. You will be allowed two hours to complete this test. No books or notes may
More informationMATH 415, WEEKS 7 & 8: Conservative and Hamiltonian Systems, Nonlinear Pendulum
MATH 415, WEEKS 7 & 8: Conservative and Hamiltonian Systems, Nonlinear Pendulum Reconsider the following example from last week: dx dt = x y dy dt = x2 y. We were able to determine many qualitative features
More informationSolutions to Dynamical Systems 2010 exam. Each question is worth 25 marks.
Solutions to Dynamical Systems exam Each question is worth marks [Unseen] Consider the following st order differential equation: dy dt Xy yy 4 a Find and classify all the fixed points of Hence draw the
More informationChapter #4 EEE8086EEE8115. 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.
More informationPolytechnic Institute of NYU MA 2132 Final Practice Answers Fall 2012
Polytechnic Institute of NYU MA Final Practice Answers Fall Studying from past or sample exams is NOT recommended. If you do, it should be only AFTER you know how to do all of the homework and worksheet
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 informationSection 5. Graphing Systems
Section 5. Graphing Systems 5A. The Phase Plane 5A1. Find the critical points of each of the following nonlinear autonomous systems. x = x 2 y 2 x = 1 x + y a) b) y = x xy y = y + 2x 2 5A2. Write each
More informationLecture 9: Eigenvalues and Eigenvectors in Classical Mechanics (See Section 3.12 in Boas)
Lecture 9: Eigenvalues and Eigenvectors in Classical Mechanics (See Section 3 in Boas) As suggested in Lecture 8 the formalism of eigenvalues/eigenvectors has many applications in physics, especially in
More informationEvenNumbered Homework Solutions
6 EvenNumbered Homework Solutions Suppose that the matric B has λ = + 5i as an eigenvalue with eigenvector Y 0 = solution to dy = BY Using Euler s formula, we can write the complexvalued solution Y
More informationGrade 11/12 Math Circles Conics & Applications The Mathematics of Orbits Dr. Shahla Aliakbari November 18, 2015
Faculty of Mathematics Waterloo, Ontario N2L 3G1 Centre for Education in Mathematics and Computing Grade 11/12 Math Circles Conics & Applications The Mathematics of Orbits Dr. Shahla Aliakbari November
More informationChapter 5 Oscillatory Motion
Chapter 5 Oscillatory Motion Simple Harmonic Motion An object moves with simple harmonic motion whenever its acceleration is proportional to its displacement from some equilibrium position and is oppositely
More informationDynamical 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
More informationENGI 4430 PDEs  d Alembert Solutions Page 11.01
ENGI 4430 PDEs  d Alembert Solutions Page 11.01 11. Partial Differential Equations Partial differential equations (PDEs) are equations involving functions of more than one variable and their partial derivatives
More informationChapter 6. Second order differential equations
Chapter 6. Second order differential equations A second order differential equation is of the form y = f(t, y, y ) where y = y(t). We shall often think of t as parametrizing time, y position. In this case
More informationFirst order differential equations
First order differential equations Samy Tindel Purdue University Differential equations and linear algebra  MA 262 Taken from Differential equations and linear algebra by Goode and Annin Samy T. First
More informationWhy are Discrete Maps Sufficient?
Why are Discrete Maps Sufficient? Why do dynamical systems specialists study maps of the form x n+ 1 = f ( xn), (time is discrete) when much of the world around us evolves continuously, and is thus well
More information7.3 Singular points and the method of Frobenius
284 CHAPTER 7. POWER SERIES METHODS 7.3 Singular points and the method of Frobenius Note: or.5 lectures, 8.4 and 8.5 in [EP], 5.4 5.7 in [BD] While behaviour of ODEs at singular points is more complicated,
More informationfor changing independent variables. Most simply for a function f(x) the Legendre transformation f(x) B(s) takes the form B(s) = xs f(x) with s = df
Physics 106a, Caltech 1 November, 2018 Lecture 10: Hamiltonian Mechanics I The Hamiltonian In the Hamiltonian formulation of dynamics each second order ODE given by the Euler Lagrange equation in terms
More information2. Secondorder Linear Ordinary Differential Equations
Advanced Engineering Mathematics 2. Secondorder Linear ODEs 1 2. Secondorder Linear Ordinary Differential Equations 2.1 Homogeneous linear ODEs 2.2 Homogeneous linear ODEs with constant coefficients
More informationLinear Planar Systems Math 246, Spring 2009, Professor David Levermore We now consider linear systems of the form
Linear Planar Systems Math 246, Spring 2009, Professor David Levermore We now consider linear systems of the form d x x 1 = A, where A = dt y y a11 a 12 a 21 a 22 Here the entries of the coefficient matrix
More informationDepartment of Mathematics IIT Guwahati
Stability of Linear Systems in R 2 Department of Mathematics IIT Guwahati A system of first order differential equations is called autonomous if the system can be written in the form dx 1 dt = g 1(x 1,
More informationLecture Notes for PHY 405 Classical Mechanics
Lecture Notes for PHY 405 Classical Mechanics From Thorton & Marion s Classical Mechanics Prepared by Dr. Joseph M. Hahn Saint Mary s University Department of Astronomy & Physics September 1, 2005 Chapter
More informationSTABILITY ANALYSIS OF DYNAMIC SYSTEMS
C. Melchiorri (DEI) Automatic Control & System Theory 1 AUTOMATIC CONTROL AND SYSTEM THEORY STABILITY ANALYSIS OF DYNAMIC SYSTEMS Claudio Melchiorri Dipartimento di Ingegneria dell Energia Elettrica e
More informationDo not turn over until you are told to do so by the Invigilator.
UNIVERSITY OF EAST ANGLIA School of Mathematics Main Series UG Examination 2013 14 CALCULUS AND MULTIVARIABLE CALCULUS MTHA4005Y Time allowed: 3 Hours Attempt QUESTIONS 1 and 2, and THREE other questions.
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 information3 2 6 Solve the initial value problem u ( t) 3. a If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b True or false and why 1. if A is
More informationMB4018 Differential equations
MB4018 Differential equations Part II http://www.staff.ul.ie/natalia/mb4018.html Prof. Natalia Kopteva Spring 2015 MB4018 (Spring 2015) Differential equations Part II 0 / 69 Section 1 SecondOrder Linear
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 information{ } is an asymptotic sequence.
AMS B Perturbation Methods Lecture 3 Copyright by Hongyun Wang, UCSC Recap Iterative method for finding asymptotic series requirement on the iteration formula to make it work Singular perturbation use
More informationNORMAL MODES, WAVE MOTION AND THE WAVE EQUATION. Professor G.G.Ross. Oxford University Hilary Term 2009
NORMAL MODES, WAVE MOTION AND THE WAVE EQUATION Professor G.G.Ross Oxford University Hilary Term 009 This course of twelve lectures covers material for the paper CP4: Differential Equations, Waves and
More informationMCE693/793: Analysis and Control of Nonlinear Systems
MCE693/793: Analysis and Control of Nonlinear Systems Systems of Differential Equations Phase Plane Analysis Hanz Richter Mechanical Engineering Department Cleveland State University Systems of Nonlinear
More informationdx n a 1(x) dy
HIGHER ORDER DIFFERENTIAL EQUATIONS Theory of linear equations Initialvalue and boundaryvalue problem nthorder initial value problem is Solve: a n (x) dn y dx n + a n 1(x) dn 1 y dx n 1 +... + a 1(x)
More informationOne 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 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 informationSolution of Additional Exercises for Chapter 4
1 1. (1) Try V (x) = 1 (x 1 + x ). Solution of Additional Exercises for Chapter 4 V (x) = x 1 ( x 1 + x ) x = x 1 x + x 1 x In the neighborhood of the origin, the term (x 1 + x ) dominates. Hence, the
More informationODE Final exam  Solutions
ODE Final exam  Solutions May 3, 018 1 Computational questions (30 For all the following ODE s with given initial condition, find the expression of the solution as a function of the time variable t You
More informationAutonomous 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
More informationFunctions, Graphs, Equations and Inequalities
CAEM DPP Learning Outcomes per Module Module Functions, Graphs, Equations and Inequalities Learning Outcomes 1. Functions, inverse functions and composite functions 1.1. concepts of function, domain and
More informationAPPLIED MATHEMATICS ADVANCED LEVEL
APPLIED MATHEMATICS ADVANCED LEVEL INTRODUCTION This syllabus serves to examine candidates knowledge and skills in introductory mathematical and statistical methods, and their applications. For applications
More informationLinear Systems of ODE: Nullclines, Eigenvector lines and trajectories
Linear Systems of ODE: Nullclines, Eigenvector lines and trajectories James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University October 6, 2013 Outline
More informationNonlinear Control. Nonlinear Control Lecture # 2 Stability of Equilibrium Points
Nonlinear Control Lecture # 2 Stability of Equilibrium Points Basic Concepts ẋ = f(x) f is locally Lipschitz over a domain D R n Suppose x D is an equilibrium point; that is, f( x) = 0 Characterize and
More informationCh 10.1: Two Point Boundary Value Problems
Ch 10.1: Two Point Boundary Value Problems In many important physical problems there are two or more independent variables, so the corresponding mathematical models involve partial differential equations.
More information1.11 Some HigherOrder Differential Equations
page 99. Some HigherOrder Differential Equations 99. Some HigherOrder Differential Equations So far we have developed analytical techniques only for solving special types of firstorder differential equations.
More information6 Second Order Linear Differential Equations
6 Second Order Linear Differential Equations A differential equation for an unknown function y = f(x) that depends on a variable x is any equation that ties together functions of x with y and its derivatives.
More informationTEST CODE: MMA (Objective type) 2015 SYLLABUS
TEST CODE: MMA (Objective type) 2015 SYLLABUS Analytical Reasoning Algebra Arithmetic, geometric and harmonic progression. Continued fractions. Elementary combinatorics: Permutations and combinations,
More informationwhere x 0 is arbitrary.
The forces internal to a system are of two types. Conservative forces, such as gravity; and dissipative forces such as friction. Internal forces arise from the natural dynamics of the system in contract
More informationNonDimensionalization of Differential Equations. September 9, 2014
NonDimensionalization of Differential Equations September 9, 2014 1 Nondimensionalization The most elementary idea of dimensional analysis is this: in any equation, the dimensions of the right hand side
More informationNonlinear dynamics & chaos BECS
Nonlinear dynamics & chaos BECS114.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
More informationSection 5.2 Series Solution Near Ordinary Point
DE Section 5.2 Series Solution Near Ordinary Point Page 1 of 5 Section 5.2 Series Solution Near Ordinary Point We are interested in second order homogeneous linear differential equations with variable
More informationHello 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
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 information!T = 2# T = 2! " The velocity and acceleration of the object are found by taking the first and second derivative of the position:
A pendulum swinging back and forth or a mass oscillating on a spring are two examples of (SHM.) SHM occurs any time the position of an object as a function of time can be represented by a sine wave. We
More informationProblem Set Number 5, j/2.036j MIT (Fall 2014)
Problem Set Number 5, 18.385j/2.036j MIT (Fall 2014) Rodolfo R. Rosales (MIT, Math. Dept.,Cambridge, MA 02139) Due Fri., October 24, 2014. October 17, 2014 1 Large µ limit for Liénard system #03 Statement:
More informationENGI Second Order Linear ODEs Page Second Order Linear Ordinary Differential Equations
ENGI 344  Second Order Linear ODEs age 01. Second Order Linear Ordinary Differential Equations The general second order linear ordinary differential equation is of the form d y dy x Q x y Rx dx dx Of
More informationMAT 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
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 informationSimple Harmonic Motion Concept Questions
Simple Harmonic Motion Concept Questions Question 1 Which of the following functions x(t) has a second derivative which is proportional to the negative of the function d x! " x? dt 1 1. x( t ) = at. x(
More informationElementary ODE Review
Elementary ODE Review First Order ODEs First Order Equations Ordinary differential equations of the fm y F(x, y) () are called first der dinary differential equations. There are a variety of techniques
More information2.2 The derivative as a Function
2.2 The derivative as a Function Recall: The derivative of a function f at a fixed number a: f a f a+h f(a) = lim h 0 h Definition (Derivative of f) For any number x, the derivative of f is f x f x+h f(x)
More informationM2A2 Problem Sheet 3  Hamiltonian Mechanics
MA Problem Sheet 3  Hamiltonian Mechanics. The particle in a cone. A particle slides under gravity, inside a smooth circular cone with a vertical axis, z = k x + y. Write down its Lagrangian in a) Cartesian,
More informationControl Systems. Internal Stability  LTI systems. L. Lanari
Control Systems Internal Stability  LTI systems L. Lanari outline LTI systems: definitions conditions South stability criterion equilibrium points Nonlinear systems: equilibrium points examples stable
More informationPHYSICS 1 Simple Harmonic Motion
Advanced Placement PHYSICS 1 Simple Harmonic Motion Student 014015 What I Absolutely Have to Know to Survive the AP* Exam Whenever the acceleration of an object is proportional to its displacement and
More informationMath Review for Exam Answer each of the following questions as either True or False. Circle the correct answer.
Math 22  Review for Exam 3. Answer each of the following questions as either True or False. Circle the correct answer. (a) True/False: If a n > 0 and a n 0, the series a n converges. Soln: False: Let
More informationB. Differential Equations A differential equation is an equation of the form
B Differential Equations A differential equation is an equation of the form ( n) F[ t; x(, xʹ (, x ʹ ʹ (, x ( ; α] = 0 dx d x ( n) d x where x ʹ ( =, x ʹ ʹ ( =,, x ( = n A differential equation describes
More information6.11 MECHANICS 3, M3 (4763) A2
6.11 MECHANICS 3, M3 (4763) A Objectives To build on the work in Mechanics 1 and Mechanics, further extending the range of mechanics concepts which students are able to use in modelling situations. The
More informationFIRSTORDER SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS III: Autonomous Planar Systems David Levermore Department of Mathematics University of Maryland
FIRSTORDER SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS III: Autonomous Planar Systems David Levermore Department of Mathematics University of Maryland 4 May 2012 Because the presentation of this material
More informationDifferential Equations 2280 Sample Midterm Exam 3 with Solutions Exam Date: 24 April 2015 at 12:50pm
Differential Equations 228 Sample Midterm Exam 3 with Solutions Exam Date: 24 April 25 at 2:5pm Instructions: This inclass exam is 5 minutes. No calculators, notes, tables or books. No answer check is
More informationTHE NCUK INTERNATIONAL FOUNDATION YEAR (IFY) Further Mathematics
IFYFM00 Further Maths THE NCUK INTERNATIONAL FOUNDATION YEAR (IFY) Further Mathematics Examination Session Summer 009 Time Allowed hours 0 minutes (Including 0 minutes reading time) INSTRUCTIONS TO STUDENTS
More information1 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
More informationThe HigginsSelkov oscillator
The HigginsSelkov oscillator May 14, 2014 Here I analyse the longtime behaviour of the HigginsSelkov oscillator. The system is ẋ = k 0 k 1 xy 2, (1 ẏ = k 1 xy 2 k 2 y. (2 The unknowns x and y, being
More informationA matrix is a rectangular array of. objects arranged in rows and columns. The objects are called the entries. is called the size of the matrix, and
Section 5.5. Matrices and Vectors A matrix is a rectangular array of objects arranged in rows and columns. The objects are called the entries. A matrix with m rows and n columns is called an m n matrix.
More informationYou may use a calculator, but you must show all your work in order to receive credit.
Math 2410010/015 Exam II April 7 th, 2017 Name: Instructions: Key Answer each question to the best of your ability. All answers must be written clearly. Be sure to erase or cross out any work that you
More information