Phase Portraits of Nonlinear Differential Equations
|
|
- Annabelle Hicks
- 5 years ago
- Views:
Transcription
1 ODE4-net.nb Phase Portraits of Nonlinear Differential Equations Nonlinear Differential Equations: x' = f (x, y) (1) Consider the system where f and g are functions of two y' = g(x, y) () variables x and y. Assume that f and g are not both linear. The functions x and y are functions of time t and the derivatives are also taken x' (t) = f (x(t), y(t)) (1) with respect to the variable t. y' (t) = g(x(t), y(t)) (). x' = ax + by (1) A linear system with invertible coefficient matrix has only y' = cx + dy () one equilibrium point (critical point, steady state solution) at (x,y)=(,). Analyzing the stability of this equilibrium solutions gives complete information about the Phase Portrait of the ODE (the trajectories). Unlike a linear system, a nonlinear system could have none, one, two, three, or any number of critical points. Like a linear system, however, the critical points are found by setting x = y =, and solve the resulting system f(x,y)= and g(x,y)=. Since there might be multiple critical points present on the phase portrait, each trajectory could be influenced by more than one critical point. This results in a much more chaotic appearance of the phase portrait. Consequently, the type and stability of each critical point need to be determined locally (in a small neighborhood on the phase plane around the critical point in question) on a caseby-case basis. The trajectories will appear very differently than those of the linear systems, except very near the critical points. General Strategy : We approximate the nonlinear system by a linear system near the equilibirum points. Suppose that (a, b) is an equilibrium point, that is f(a, b) = g(a, b) =. x' = f (x, y) f (a.b) + ' x f (a, b) (x ( a) + ' y f (a, b) (y ( b) = ' x f (a, b) (x ( a) + ' y f (a, b) (y ( b) y' = g(x, y) g(a, b) + ' x g(a, b) (x ( a) + ' y g(a, b) (y ( b) = ' x g(a, b) (x ( a) + ' y g(a, b) (y ( b) The critical point can be translated to via the coordinate change x := x - a and y := y - b and still retains all of its properties. The approximated system becomes a linear system: x' = ' x f (a, b) x + ' y f (a, b) y (1) with coefficient matrix y' = ' x g(a, b) x + ' y g(a, b) y () J (a,b) = x f(a, b) y f(a, b) x g(a, b) y g(a, b). The matrix J (x,y) = x f(x, y) y f(x, y) is called the Jacobian x g(x, y) y g(x, y) matrix of the system. It contains the first order partial derivatives of f and g evaluated at the point (x,y). We can use the eigenvalues of the matrix Jacobian matrix to decide the type and stability of the critical point (a,b). Let " and # be the two eigenvalues. Assume that " and # are different from. Eigenvalues Linear System NonlinearSyste λ, μ!! Real λ > μ > Nodal Source (Unstable) Nodal Source (Unstab λ < μ < Nodal Sink (Stable) Nodal Sink (Stable) λ > > μ Saddle point (Unstable) Saddle point (Unstable λ = μ > Degenerate Source or Nodal Source (Unstable) dependingof the geometric multiplicity of # Source (Degenerate, Nodal, Spiral Sou dependingon the nonlinearterms and the λ = μ < Degenerate Sink or Nodal Sink (Stable) dependingof the geometric multiplicity of # (1) () Sink (Degenerate, Nodal, Spiral Sink) (S dependingon the nonlinearterms and the
2 ODE4-net.nb 3 4 ODE4-net.nb λ, μ!! Complex Re (λ) > Spiral Source (Unstable) Spiral Source (Unstab Re (λ) < Spiral Sink (Stable) Spiral Sink (Stable) Re (λ) = Center (Stable) Center, Spiral Sink, Spiral Source (Stability cannot be determined based o Check also the Bifurcation Diagram on Wikipedia. Hyperbolic critical point Definition: The equilibrium is said to be hyperbolic if all eigenvalues of the Jacobian matrix have non - zero real parts. Hyperbolic equilibria are robust : Small perturbations do not change qualitatively the phase portrait near the equilibria. Moreover, local phase portrait of a hyperbolic equilibrium of a nonlinear system is equivalent to that of its linearization. This statement has a mathematically precise form known as the Hartman - Grobman Theorem. Examples of nonlinear differential equations: Example 1. Solve the system x' = x(1 ( y) (1) y' = y(x ( 1) () Usually an explicit solution cannot be found, but you can still ask Mathematica to find a numeric solution using the NDSolve command. DSolve[{x'[t]. x[t] (1 / y[t]), y'[t]. y[t] (x[t] / 1), x[]. 1, y[]. 1}, {x[t], y[t]}, t] Solve::ifun: Inverse functionsare beingused by Solve, so some solutionsmay notbe found; use Reduce forcompletesolutioninformation. ) Solve::ifun: Inverse functionsare beingused by Solve, so some solutionsmay notbe found; use Reduce forcompletesolutioninformation. ) Solve::ifun: Inverse functionsare beingused by Solve, so some solutionsmay notbe found; use Reduce forcompletesolutioninformation. ) General::stop: Furtheroutputof Solve::ifunwillbe suppressedduringthiscalculation. ) DSolve::bvnul: For some branchesof thegeneralsolution, thegivenboundaryconditionslead toan emptysolution. ) {} NDSolve[{x'[t]. x[t] (1 / y[t]), y'[t]. y[t] (x[t] / 1), x[]., y[]. 1}, {x[t], y[t]}, {t,, 1}] x[t] & InterpolatingFunction y[t] & InterpolatingFunction Domain: {{., 1.}} Domain: {{., 1.}} [t], [t] The most reliable information here comes from the Vector Field Plot. 4 # #4 StreamPlot[{x (1 / y), y (x / 1)}, {x, /5, 5}, {y, /5, 5}, StreamScale 1 Automatic, StreamPoints 1 {{{{4, 1}, Red}, Automatic}}] #4 # 4 What is the nature of the equilibrium solutions?
3 ODE4-net.nb 5 6 ODE4-net.nb What is the nature of the equilibrium solutions? 1. Set x' = and y' =. In this case, it is easy to solve the equation by hand: x(1- y)= and y(x-1)= implies (x,y)=(,) or (x,y)=(1,1). For more complicated cases, we will use Solve[..] or NSolve[..] EqPoints = Solve[{x (1 / y)., y (x / 1). }, {x, y}] {{x &, y & }, {x & 1, y & 1}}. Next we find the Jacobian matrix of the system. In this case, it is easy to compute it by hand: J (x,y) = x (x (1 / y)) y (x (1 / y)) 1 / y /x = x (y (x / 1)) y (y (x / 1)) y x / 1 In more complicated cases, we will ask Mathematica to compute the Jacobian matrix for us, as follows : The commands D[f[x], {x}] or f'[x] give the dervative of f with respect to x. If f is a function of several variables, D[f[x, y], {x}] or y (f[x, y]) gives the partial derivative of f with respect to x. Try them out! 3. We evaluate the Jacobian matrix at the two equilibrium points and find its eigenvalues: B = J[f] 5. Flatten[EqPoints[[1]]]; Print["The Jacobian Matrix at the equilibrium point (,) is ", MatrixForm[B]] Print["The eigenvalues of B are ", Eigenvalues[B]] The Jacobian Matrix at the equilibrium point (,) is 1 *1 The eigenvalues of B are {*1, 1} T = J[f] 5. Flatten[EqPoints[[]]]; Print["The Jacobian Matrix at the equilibrium point (1,1) is ", MatrixForm[ T]] Print["The eigenvalues are ", Eigenvalues[T]] The Jacobian Matrix at the equilibrium point (1,1) is *1 1 The eigenvalues are {, *} y x y + y^ D[y^, {y}] The following code computes the Jacobian matrix: Clear[f, J, x, y] f[x_, y_] := {x (1 / y), y (x / 1)} J[f_] := Module[{A = Table[,, ]}, A[[1, 1]] = D[f[x, y][[1]], {x}]; A[[1, ]] = D[f[x, y][[1]], {y}]; A[[, 1]] = D[f[x, y][[]], {x}]; A[[, ]] = D[f[x, y][[]], {y}]; Return[A]] J[f] {{1 * y, *x}, {y, *1 + x}}
4 ODE4-net.nb 7 8 ODE4-net.nb By analyzing the Real and Imaginary Part of the Eigenvalues, we conclude that the point (, ) is a saddle (the Jacobian matrix has two real eigenvalues, one smaller than, one greater than ) whereas the point (1, 1) is a center (the Jacobian matrix has two complex conjugate eigenvalues with Real part equal to, and the Phase Portrait depicts closed circular trajectories). This gives us information about the shape of trajectories in a vecinity of the equilibrium solutions. Example : 1 StreamPlot[{y / x (x^ + y^), /x / y (x^ + y^) }, {x, /, }, {y, /, }, StreamPoints 1 {{{{1, }, Red}, Automatic}}] x' = y ( x(x^ + y^) (1) Consider the system. Study the nature of the y' = (x ( y(x^ + y^) () critical point (x,y)=(,). Clear[h] h[x_, y_] := {y / x (x^ + y^), /x / y (x^ + y^)} NSolve[{y / x (x^ + y^)., /x / y (x^ + y^). }, {x, y}] MatrixForm[J[h]] Eigenvalues[J[h] 5. {x 1, y 1 }] {{x &., y &.}} *3 x * y 1 * x y *1 * x y *x * 3 y {, *} #1 # # #1 1
5 ODE4-net.nb 9 1 ODE4-net.nb sol = NDSolve[{x'[t] == y[t] / x[t] (x[t]^ + y[t]^), y'[t] == /x[t] / y[t] (x[t]^ + y[t]^), x[]. 1, y[]. }, {x[t], y[t]}, {t,, 1}] ParametricPlot[{x[t], y[t]} 5. Flatten[sol], {t,, 1}, PlotPoints 1 1] StreamPlot[{y + x (x^ + y^), /x + y (x^ + y^) }, {x, /, }, {y, /, }] x[t] & InterpolatingFunction Domain:., [t], 1 y[t] & InterpolatingFunction Domain:., [t].3. #1.1 # #.4 #.3 #. # # #1 1 #.1 #. #.3 x(y / x (x^ + y^)) + y(/x / y (x^ + y^)) *4 x * 4 y Clear[h] h[x_, y_] := {y + x (x^ + y^), /x + y (x^ + y^)} MatrixForm[J[h]] Eigenvalues[J[h] 5. {x 1, y 1 }] 3 x + y 1 + x y *1 + x y x + 3 y {, *}
6 ODE4-net.nb 11 1 ODE4-net.nb Conclusions : The linearized system has a center at the origin. The nonlinear system does not have a center. Example 3 Consider the system and study their nature. x' = x + y ( x(x^ + y^) (1). Find all critical points y' = (x + y ( y(x^ + y^) () Clear[h] h[x_, y_] := {x + y / x (x^ + y^), /x + y / y (x^ + y^)} NSolve[{y / x (x^ + y^)., /x / y (x^ + y^). }, {x, y}] MatrixForm[J[h]] Eigenvalues[J[h] 5. {x 1, y 1 }] StreamPlot[h[x, y], {x, /, }, {y, /, }] {{x &., y &.}} 1 * 3 x * y 1 * x y *1 * x y 1 * x * 3 y {1 +, 1 * } Conclusions : The origin is an spiral source (unstable spiral point) for both the linear system and the nonlinear system. One may therefore think that all trajectories different from should spiral out to infinity. This is clearly not the case, as far away from the origin, the arrows are pointing inward. Moreover, there is a closed trajectory (a circle) and all trajectories different from (,) spiral towards this circle (limit cycle). You can find the equation of this limit cycle by switching to polar coordinates x=rcos', y=rsin'. The literature on the existence and properties of limits cycles is very rich: Poincare - Bendixon Theorem Suppose f and g have continuous partial derivatives. a) Every closed trajectory must enclose at least one critical point. If it encloses exactly one critical point, then this cannot be a saddle point. b) If ( x f + ( y g has the same sign in a disk D, then there can be no closed trajectory of the system lying entirely inside D. c) Suppose that U is a connected open subset that contains no critical point of the system. Let U denote the set U together with its boundary. If the vector field in U points towards the interior of U, then U contains a closed trajectory. d) Suppose that U is a connected open subset that contains no critical point of the system. Let U denote the set U together with its boundary. If there exists some trajectory that never leaves the set U (is trapped inside U), then U must necessarily contain a closed trajectory. 1 Exercise : x' = y (1) Consider the system. Find all critical points (x,y) of the y' = (y ( sin(x) () system which belong to the region -$ < x,y < $. Study the nature of these critical points (type, stability). Use only Mathematica commands, even when computing the solutions! #1 # # #1 1 Conclusions : The origin is an spiral source (unstable spiral point) for both the linear system and the nonlinear system. One may therefore think that all trajectories different from should spiral out to infinity. This is clearly not the case, as far away from the origin, the arrows are pointing inward. Moreover, there is a closed trajectory (a circle) and all trajectories different from (,) spiral towards this circle (limit cycle). You can find the equation of this limit cycle by switching to polar coordinates x=rcos', y=rsin'.
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
More information2.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
More informationProblem 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
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 first-order system for x(t) d dt x = f(x,t), x(0) = x 0 In particular, consider
More informationEE222 - Spring 16 - Lecture 2 Notes 1
EE222 - Spring 16 - Lecture 2 Notes 1 Murat Arcak January 21 2016 1 Licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License. Essentially Nonlinear Phenomena Continued
More informationPractice Problems for Final Exam
Math 1280 Spring 2016 Practice Problems for Final Exam Part 2 (Sections 6.6, 6.7, 6.8, and chapter 7) S o l u t i o n s 1. Show that the given system has a nonlinear center at the origin. ẋ = 9y 5y 5,
More information6.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
More informationClassification of Phase Portraits at Equilibria for u (t) = f( u(t))
Classification of Phase Portraits at Equilibria for u t = f ut Transfer of Local Linearized Phase Portrait Transfer of Local Linearized Stability How to Classify Linear Equilibria Justification of the
More informationNonlinear 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
More informationMathematical Modeling I
Mathematical Modeling I Dr. Zachariah Sinkala Department of Mathematical Sciences Middle Tennessee State University Murfreesboro Tennessee 37132, USA November 5, 2011 1d systems To understand more complex
More informationLecture 38. Almost Linear Systems
Math 245 - Mathematics of Physics and Engineering I Lecture 38. Almost Linear Systems April 20, 2012 Konstantin Zuev (USC) Math 245, Lecture 38 April 20, 2012 1 / 11 Agenda Stability Properties of Linear
More informationE209A: 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.
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 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 integro-differential equations
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, y-plane which satisfy this
More informationENGI Linear Approximation (2) Page Linear Approximation to a System of Non-Linear ODEs (2)
ENGI 940 4.06 - Linear Approximation () Page 4. 4.06 Linear Approximation to a System of Non-Linear ODEs () From sections 4.0 and 4.0, the non-linear system dx dy = x = P( x, y), = y = Q( x, y) () with
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 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 information= 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
More informationMath 312 Lecture Notes Linearization
Math 3 Lecture Notes Linearization Warren Weckesser Department of Mathematics Colgate University 3 March 005 These notes discuss linearization, in which a linear system is used to approximate the behavior
More informationMATH 415, WEEKS 7 & 8: Conservative and Hamiltonian Systems, Non-linear Pendulum
MATH 415, WEEKS 7 & 8: Conservative and Hamiltonian Systems, Non-linear 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 informationOutline. Learning Objectives. References. Lecture 2: Second-order Systems
Outline Lecture 2: Second-order Systems! Techniques based on linear systems analysis! Phase-plane analysis! Example: Neanderthal / Early man competition! Hartman-Grobman theorem -- validity of linearizations!
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 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 informationA 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
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 Half-line solutions Unstable solution Stable
More informationAn 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
More informationSolutions to Math 53 Math 53 Practice Final
Solutions to Math 5 Math 5 Practice Final 20 points Consider the initial value problem y t 4yt = te t with y 0 = and y0 = 0 a 8 points Find the Laplace transform of the solution of this IVP b 8 points
More informationMath 273 (51) - Final
Name: Id #: Math 273 (5) - Final Autumn Quarter 26 Thursday, December 8, 26-6: to 8: Instructions: Prob. Points Score possible 25 2 25 3 25 TOTAL 75 Read each problem carefully. Write legibly. Show all
More informationB5.6 Nonlinear Systems
B5.6 Nonlinear Systems 4. Bifurcations Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Local bifurcations for vector fields 1.1 The problem 1.2 The extended centre
More informationAPPM 2360: Final Exam 10:30am 1:00pm, May 6, 2015.
APPM 23: Final Exam :3am :pm, May, 25. ON THE FRONT OF YOUR BLUEBOOK write: ) your name, 2) your student ID number, 3) lecture section, 4) your instructor s name, and 5) a grading table for eight questions.
More informationDynamics and Bifurcations in Predator-Prey Models with Refuge, Dispersal and Threshold Harvesting
Dynamics and Bifurcations in Predator-Prey Models with Refuge, Dispersal and Threshold Harvesting August 2012 Overview Last Model ẋ = αx(1 x) a(1 m)xy 1+c(1 m)x H(x) ẏ = dy + b(1 m)xy 1+c(1 m)x (1) where
More informationFIRST-ORDER SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS III: Autonomous Planar Systems David Levermore Department of Mathematics University of Maryland
FIRST-ORDER 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 informationChapter #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.
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 2009 Farzaneh Abdollahi Nonlinear Control Lecture 2 1/68
More informationMath 1270 Honors ODE I Fall, 2008 Class notes # 14. x 0 = F (x; y) y 0 = G (x; y) u 0 = au + bv = cu + dv
Math 1270 Honors ODE I Fall, 2008 Class notes # 1 We have learned how to study nonlinear systems x 0 = F (x; y) y 0 = G (x; y) (1) by linearizing around equilibrium points. If (x 0 ; y 0 ) is an equilibrium
More informationMATH 215/255 Solutions to Additional Practice Problems April dy dt
. For the nonlinear system MATH 5/55 Solutions to Additional Practice Problems April 08 dx dt = x( x y, dy dt = y(.5 y x, x 0, y 0, (a Show that if x(0 > 0 and y(0 = 0, then the solution (x(t, y(t of the
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 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 (Springer-Verlag
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 informationPhysics: spring-mass system, planet motion, pendulum. Biology: ecology problem, neural conduction, epidemics
Applications of nonlinear ODE systems: Physics: spring-mass system, planet motion, pendulum Chemistry: mixing problems, chemical reactions Biology: ecology problem, neural conduction, epidemics Economy:
More informationMATH 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
More informationMath 232, Final Test, 20 March 2007
Math 232, Final Test, 20 March 2007 Name: Instructions. Do any five of the first six questions, and any five of the last six questions. Please do your best, and show all appropriate details in your solutions.
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 informationMathQuest: Differential Equations
MathQuest: Differential Equations Geometry of Systems 1. The differential equation d Y dt = A Y has two straight line solutions corresponding to [ ] [ ] 1 1 eigenvectors v 1 = and v 2 2 = that are shown
More informationOne 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
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 informationThere is a more global concept that is related to this circle of ideas that we discuss somewhat informally. Namely, a region R R n with a (smooth)
82 Introduction Liapunov Functions Besides the Liapunov spectral theorem, there is another basic method of proving stability that is a generalization of the energy method we have seen in the introductory
More 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 informationE209A: Analysis and Control of Nonlinear Systems Problem Set 6 Solutions
E9A: Analysis and Control of Nonlinear Systems Problem Set 6 Solutions Michael Vitus Gabe Hoffmann Stanford University Winter 7 Problem 1 The governing equations are: ẋ 1 = x 1 + x 1 x ẋ = x + x 3 Using
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 informationCALCULUS MATH*2080 SAMPLE FINAL EXAM
CALCULUS MATH*28 SAMPLE FINAL EXAM Sample Final Exam Page of 2 Prof. R.Gentry Print Your Name Student No. SIGNATURE Mark This exam is worth 45% of your final grade. In Part I In Part II In part III In
More information11 Chaos in Continuous Dynamical Systems.
11 CHAOS IN CONTINUOUS DYNAMICAL SYSTEMS. 47 11 Chaos in Continuous Dynamical Systems. Let s consider a system of differential equations given by where x(t) : R R and f : R R. ẋ = f(x), The linearization
More informationLECTURE 8: DYNAMICAL SYSTEMS 7
15-382 COLLECTIVE INTELLIGENCE S18 LECTURE 8: DYNAMICAL SYSTEMS 7 INSTRUCTOR: GIANNI A. DI CARO GEOMETRIES IN THE PHASE SPACE Damped pendulum One cp in the region between two separatrix Separatrix Basin
More informationB5.6 Nonlinear Systems
B5.6 Nonlinear Systems 5. Global Bifurcations, Homoclinic chaos, Melnikov s method Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Motivation 1.1 The problem 1.2 A
More informationCalculus for the Life Sciences II Assignment 6 solutions. f(x, y) = 3π 3 cos 2x + 2 sin 3y
Calculus for the Life Sciences II Assignment 6 solutions Find the tangent plane to the graph of the function at the point (0, π f(x, y = 3π 3 cos 2x + 2 sin 3y Solution: The tangent plane of f at a point
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 information4 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
More informationA review of stability and dynamical behaviors of differential equations:
A review of stability and dynamical behaviors of differential equations: scalar ODE: u t = f(u), system of ODEs: u t = f(u, v), v t = g(u, v), reaction-diffusion equation: u t = D u + f(u), x Ω, with boundary
More information. For each initial condition y(0) = y 0, there exists a. unique solution. In fact, given any point (x, y), there is a unique curve through this point,
1.2. Direction Fields: Graphical Representation of the ODE and its Solution Section Objective(s): Constructing Direction Fields. Interpreting Direction Fields. Definition 1.2.1. A first order ODE of the
More informationSample Solutions of Assignment 10 for MAT3270B
Sample Solutions of Assignment 1 for MAT327B 1. For the following ODEs, (a) determine all critical points; (b) find the corresponding linear system near each critical point; (c) find the eigenvalues of
More informationODE, part 2. Dynamical systems, differential equations
ODE, part 2 Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2011 Dynamical systems, differential equations Consider a system of n first order equations du dt = f(u, t),
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 Dynamical Systems
Introduction to Dynamical Systems Autonomous Planar Systems Vector form of a Dynamical System Trajectories Trajectories Don t Cross Equilibria Population Biology Rabbit-Fox System Trout System Trout System
More informationEN Nonlinear Control and Planning in Robotics Lecture 3: Stability February 4, 2015
EN530.678 Nonlinear Control and Planning in Robotics Lecture 3: Stability February 4, 2015 Prof: Marin Kobilarov 0.1 Model prerequisites Consider ẋ = f(t, x). We will make the following basic assumptions
More informationStable and Unstable Manifolds for Planar Dynamical Systems
Stable and Unstable Manifolds for Planar Dynamical Systems H.L. Smith Department of Mathematics Arizona State University Tempe, AZ 85287 184 October 7, 211 Key words and phrases: Saddle point, stable 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 informationSolutions to Final Exam Sample Problems, Math 246, Spring 2011
Solutions to Final Exam Sample Problems, Math 246, Spring 2 () Consider the differential equation dy dt = (9 y2 )y 2 (a) Identify its equilibrium (stationary) points and classify their stability (b) Sketch
More informationQuestion: Total. Points:
MATH 308 May 23, 2011 Final Exam Name: ID: Question: 1 2 3 4 5 6 7 8 9 Total Points: 0 20 20 20 20 20 20 20 20 160 Score: There are 9 problems on 9 pages in this exam (not counting the cover sheet). Make
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 informationẋ = f(x, y), ẏ = g(x, y), (x, y) D, can only have periodic solutions if (f,g) changes sign in D or if (f,g)=0in D.
4 Periodic Solutions We have shown that in the case of an autonomous equation the periodic solutions correspond with closed orbits in phase-space. Autonomous two-dimensional systems with phase-space R
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 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 informationComputational Neuroscience. Session 4-2
Computational Neuroscience. Session 4-2 Dr. Marco A Roque Sol 06/21/2018 Two-Dimensional Two-Dimensional System In this section we will introduce methods of phase plane analysis of two-dimensional systems.
More informationChapter 4: First-order differential equations. Similarity and Transport Phenomena in Fluid Dynamics Christophe Ancey
Chapter 4: First-order differential equations Similarity and Transport Phenomena in Fluid Dynamics Christophe Ancey Chapter 4: First-order differential equations Phase portrait Singular point Separatrix
More informationSolutions for B8b (Nonlinear Systems) Fake Past Exam (TT 10)
Solutions for B8b (Nonlinear Systems) Fake Past Exam (TT 10) Mason A. Porter 15/05/2010 1 Question 1 i. (6 points) Define a saddle-node bifurcation and show that the first order system dx dt = r x e x
More informationMath 331 Homework Assignment Chapter 7 Page 1 of 9
Math Homework Assignment Chapter 7 Page of 9 Instructions: Please make sure to demonstrate every step in your calculations. Return your answers including this homework sheet back to the instructor as a
More informationHalf of Final Exam Name: Practice Problems October 28, 2014
Math 54. Treibergs Half of Final Exam Name: Practice Problems October 28, 24 Half of the final will be over material since the last midterm exam, such as the practice problems given here. The other half
More informationMathematical Foundations of Neuroscience - Lecture 7. Bifurcations II.
Mathematical Foundations of Neuroscience - Lecture 7. Bifurcations II. Filip Piękniewski Faculty of Mathematics and Computer Science, Nicolaus Copernicus University, Toruń, Poland Winter 2009/2010 Filip
More informationLocal Phase Portrait of Nonlinear Systems Near Equilibria
Local Phase Portrait of Nonlinear Sstems Near Equilibria [1] Consider 1 = 6 1 1 3 1, = 3 1. ( ) (a) Find all equilibrium solutions of the sstem ( ). (b) For each equilibrium point, give the linear approimating
More informationCalculus and Differential Equations II
MATH 250 B Second order autonomous linear systems We are mostly interested with 2 2 first order autonomous systems of the form { x = a x + b y y = c x + d y where x and y are functions of t and a, b, c,
More informationBIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs
BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs Yuri A. Kuznetsov August, 2010 Contents 1. Solutions and orbits. 2. Equilibria. 3. Periodic orbits and limit cycles. 4. Homoclinic orbits.
More informationFundamentals 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
More informationPhase Plane Analysis
Phase Plane Analysis Phase plane analysis is one of the most important techniques for studying the behavior of nonlinear systems, since there is usually no analytical solution for a nonlinear system. Background
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 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 informationUnderstand the existence and uniqueness theorems and what they tell you about solutions to initial value problems.
Review Outline To review for the final, look over the following outline and look at problems from the book and on the old exam s and exam reviews to find problems about each of the following topics.. Basics
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 informationSection 3.5 The Implicit Function Theorem
Section 3.5 The Implicit Function Theorem THEOREM 11 (Special Implicit Function Theorem): Suppose that F : R n+1 R has continuous partial derivatives. Denoting points in R n+1 by (x, z), where x R n and
More informationLotka Volterra Predator-Prey Model with a Predating Scavenger
Lotka Volterra Predator-Prey Model with a Predating Scavenger Monica Pescitelli Georgia College December 13, 2013 Abstract The classic Lotka Volterra equations are used to model the population dynamics
More informationChapter 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)
More informationProblem List MATH 5173 Spring, 2014
Problem List MATH 5173 Spring, 2014 The notation p/n means the problem with number n on page p of Perko. 1. 5/3 [Due Wednesday, January 15] 2. 6/5 and describe the relationship of the phase portraits [Due
More informationMCE/EEC 647/747: Robot Dynamics and Control. Lecture 8: Basic Lyapunov Stability Theory
MCE/EEC 647/747: Robot Dynamics and Control Lecture 8: Basic Lyapunov Stability Theory Reading: SHV Appendix Mechanical Engineering Hanz Richter, PhD MCE503 p.1/17 Stability in the sense of Lyapunov A
More informationHigher Order Singularities in Piecewise Linear Vector Fields
Higher Order Singularities in Piecewise Linear Vector Fields Xavier Tricoche, Gerik Scheuermann, Hans Hagen Computer Science Department, University of Kaiserslautern, Germany Summary. Piecewise linear
More informationThe Derivative. Appendix B. B.1 The Derivative of f. Mappings from IR to IR
Appendix B The Derivative B.1 The Derivative of f In this chapter, we give a short summary of the derivative. Specifically, we want to compare/contrast how the derivative appears for functions whose domain
More informationMath 312 Lecture Notes Linear Two-dimensional Systems of Differential Equations
Math 2 Lecture Notes Linear Two-dimensional Systems of Differential Equations Warren Weckesser Department of Mathematics Colgate University February 2005 In these notes, we consider the linear system of
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 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 informationThird In-Class Exam Solutions Math 246, Professor David Levermore Thursday, 3 December 2009 (1) [6] Given that 2 is an eigenvalue of the matrix
Third In-Class Exam Solutions Math 26, Professor David Levermore Thursday, December 2009 ) [6] Given that 2 is an eigenvalue of the matrix A 2, 0 find all the eigenvectors of A associated with 2. Solution.
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 in-class exam is 5 minutes. No calculators, notes, tables or books. No answer check is
More information