E209A: Analysis and Control of Nonlinear Systems Problem Set 3 Solutions
|
|
- Phoebe Bell
- 5 years ago
- Views:
Transcription
1 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. Consider the closed contour shown in Figure. For the contour, which is positively oriented in the counter-clockwise direction, we find that I f D) = dθ f x) = ) π The only equilibrium point in D is an unstable node. Therefore, If x i ) = ) J which is consistent with the equation ). Figure : Problem a - contour
2 Figure 3: Problem b Part b Figure 4: Problem b - contour This phase portrait is impossible. Consider the closed contour shown in Figure 4. For the contour, which is positively oriented in the counter-clockwise direction, we find that I f D) = dθ f x) = 3) π Inside the contour, there are two equilibria, which are both foci. Therefore, If x i ) = 4) J which is inconsistent with equation 3). Figure shows the correct phase portrait in which the right middle arrow has been reversed. Figure : Problem b - fixed Part c This phase portrait is impossible. Consider the closed contour shown in Figure 7. For the contour, which is positively oriented in the counter-clockwise direction, we find that I f D) = dθ f x) = ) π Inside the contour, there are 8 equilibria of which 4 are saddles. Therefore, If x i ) = 0 6) J
3 Figure 6: Problem c Figure 7: Problem c - contour which is inconsistent with Eqn. ). To correct the phase portrait, we need to add an equilibrium point inside the middle circle. This equilibrium can be either an unstable/stable node, unstable/stable focus or a center. Problem The governing equations are x = x ) x = x + x x x 7) The equilibrium of the system is 0, 0). To analyze the system, lets define the function V x) = x + x ) which is proportional to the squared distance from a point x, x ) to the origin. Now we take the derivative of V x) w.r.t. time along the system s trajectory [ )] V x) = x ẋ + x ẋ = x x + x x + x x x V x) = x x x) 9) The purpose of V x) is to show how the system s trajectories are evolving w.r.t. the contours of V x). If V x) < 0 then the trajectories of the system are moving closer to 0, 0). If V x) > 0 then the trajectories of the system are moving away from 0, 0). 8) Let s define a region such that x + x 0) 3
4 Let s first analyze the outer boundary in which x + x =, plugging this into Eqn. 9) and simplifying we obtain V x) = x 4 0 ) which shows that trajectories either enter or orbit the region. Now let s analyze the inner boundary in which x + x =, plugging this into Eqn. 9) and simplifying we obtain V x) = x x 0 ) which shows that trajectories either enter or orbit the region. Therefore we have shown that the annulus is invariant. By Poincare -Bendixson s theorem, there must be either a closed orbit or an equilibrium point inside, but the only equilibrium point is at 0, 0) which is outside the annulus. Therefore, the annulus must contain a closed orbit. Problem 3 Part i) The governing equations are ẋ = x + x ẋ = g x ) + ax, where a is a constant, a 3) Applying Bendixson s theorem, divf) = f x + f x = + a 4) which is neither zero nor changes sign in all of R so there are no closed orbits or limit cycles. Part ii) The governing equations are ẋ = x x ẋ = x ) Solving for the equiibria results in the x -axis being an equilibrium set. We know that a closed orbit cannot cross the x -axis, therefore it must be either be in the upper or lower halves of the plane. Since there are no equilibrium points other than the x -axis, we can conclude that there are no closed orbits as a result from index theory that a closed orbit must enclose an equilibrium point. Problem 4 The governing equation is Letting x = x and x = ẋ results in ẍ + x x) + ẋ = 0 6) ẋ = x ẋ = x x ) x 7) 4
5 Part a) Solving for the equilibria yields 0, 0) and, 0). To determine the stability of the equilibrium points, we will need to linearize the state model about them. The Jacobian is: [ ] 0 Df = 8) x Substituting in for the equilibrium point at 0, 0): 0 Df 0,0) = The eigenvalues of Eqn. 9 are λ =, which is a stable node. 9) Substituting in for the equilibrium point at, 0) yields 0 Df,0) = 0) The eigenvalues of Eqn. 0 are λ = ± which is a saddle. Part b) For this part of the problem we need to show that the region in Figure 8 is positively invariant. Figure 8: Problem 4 - Invariant Region On x = 0, the governing equations are for x [0, ], ẋ 0. ẋ = 0 ẋ = x x ) ) On x =, the governing equations are for x [0, ], ẋ 0 and ẋ 0. ẋ = x ẋ = x )
6 On x = x, the inward pointing normal on this line with respect to K is [ ] T. Therefore, the magnitude of f along the vector [ ] T is [ ] x [ ] = x x x ) x x x ) x = x 3) for x = x, x 0. Since the magnitude is greater than zero, we know that the trajectories are pointing in toward K along the line x = x. Therefore, the region K is positively invariant. Part c) By Poincare -Bendixson, all trajectories inside K converge to either a closed orbit or an equilibrium point. There are no closed orbits in K since ẋ < 0 in K, or by index theory - if there were a closed orbit it would surround an equilibrium point and there are no equilibrium points in the interior of K. Therefore, all trajectories not at, 0) converge to 0, 0). Problem Model The governing equations of the system are ẋ = ax bxy ẏ = dy + cxy 4) The equilibria of the system are: 0, 0) and d c, ) a b. To determine the stability of the equilibrium points, we will need to linearize the state model about them. The Jacobian is: a by bx Df = ) cy d + cx For 0, 0): Df 0,0) = [ a 0 0 d which has eigenvalues of λ = a, d and therefore is a saddle point. ] 6) For d c, a b ) : 0 Df d c, a b ) = db c ac b 0 which has eigenvalues of λ = ± adj. By the Hartman-Grobman theorem, we can t trust the linearization since the eigenvalues are on the jω-axis. As you will see from the simulation, Figure 9, the equilibrium point at d c, a b ) is indeed a center. The intuition behind this is that any perturbations from steady state will cause cyclic variations. Consider when sardines are low and sharks are high, the sharks will start to die off because their food source is low. Since the sharks are dying off, the sardines will start to increase in population because there are less predators around. With a large sardine population and low shark population, the sharks will flourish and increase in number which in turn decreases the sardine population. Now this cyclic trend will repeat. 7) 6
7 Figure 9: Problem - Model 7
8 Model The governing equations of the system are The equilibrium points of the system are: 0, 0) ẋ = a by λx)x ẏ = d + cx µy)y ) a λ, 0) bd+aµ bc+λµ, ac dλ bc+λµ To determine the stability of the equilibrium points, we will need to linearize the state model about them. The Jacobian is: a by λx bx Df = 30) cy d + cx µy For 0, 0): Df = [ a 0 0 d which has eigenvalues of λ = a, d and therefore is a saddle point. For a λ, 0) : which has eigenvalues of λ = a, ac λ ] 8) 9) 3) a ba Df = λ 0 d + ca 3) λ d. If ac > λd then it is a saddle. If ac λd then it is a stable node. ) For bd+aµ bc+λµ, ac dλ bc+λµ : If ac < λd, then the equilibrium point will be negative which can t happen. If ac λd, then the linearization is: Plugging in x e, y e ) = The characteristic equation is: which simplifies to Defining Since x e > 0, y e > 0, and λ, µ, b > 0, a bye λx Df xe,y e) = e bx e cy e d + cx e µy e ) bd+aµ bc+λµ, ac dλ bc+λµ and simplifying the expression we obtain: λxe bx Df xe,y e) = e cy e µy e Therefore, the equilibrium will always be stable, and if 33) 34) s + λx e ) s + µy e ) + bcx e y e = 0 3) s + λx e + µy e ) s + λµ + bc) x e y e = 0 36) α 4β > 0 α 4β < 0 α 4β = 0 β = λx e + µy e α = λµ + bc) x e y e 37) β > 0 α > 0 38) stable node stable focus improper stable node 39) 8
9 Figure 0 depicts an example solution. The plot depicts that there always exists a stable population of sharks and sardines. No matter what the initial populations are, they always reach the stable equilibrium. Figure 0: Problem - Model Problem 6 The governing equations are ẋ = a x 4xx +x ) 40) ẋ = bx x +x Part a) Solve for the equilibrium points. x from Eqn. 4), x = 0 or =. +x If x = 0, then from Eqn. 4) a = 0 which is a contradiction. x If =, then from Eqn. 4) a x +x 4x = 0 which yields x = a. Therefore there is a unique equilibrium at a a ) ), + a x 4x x + x = 0 4) bx x ) + x = 0 4) To determine the stability of the equilibrium points, we will need to linearize the state model about them. The Jacobian is: 4x + 8x +x x 4x +x Df = ) +x ) b x + bx +x x bx 44) +x ) +x 43) 9
10 simplifying yields, Df = + x a Evaluating the jacobian at, + ) ) a is Df = [ + ) a The characteristic equation of Eqn. 46 is a ) s a b This if we can choose a value for b such that Solving for this value of b yields, + 3x 4x bx bx 3 a + 3 a ) 4 a b a ) b a ] 4) 46) ) s + b a a ) ) + = 0 47) ) ) + b a < 0 then the equilibrium point is unstable. b < 3 a a 48) Part b) Now we need to show the invariance of the region M = { x, x ) x 0, x 0, x a, x + a }. For x = 0, ẋ = a > 0 For x = 0, ẋ = bx > 0 if x > 0 and at 0, 0), ẋ > 0 For x = a, ẋ < 0 for x > 0 and at a, 0), ẋ > 0 For x = + a, ẋ < 0 if 0 < x < a Therefore the region M is invariant. Since M is a closed and bounded region, and for b < 3 a a contains only one unstable equilibrium point there must be a closed orbit inside M by the Poincare -Bendixson theorem. Problem 7 Part a) The governing equations are where σ, λ, and b are positive constants. ẋ = σx x ) ẋ = + λ x 3 )x x ẋ 3 = x x bx 3 49) To solve for the equilibrium points, first set ẋ = 0 which yields x = x. Plugging that into ẋ = 0, yields x 3 = λ. Finally, plugging that into ẋ 3 = 0 gives x = x = ± bλ. Also, by inspection 0, 0, 0) is an equilibrium point. Therefore the three equilibrium points are: 0, 0, 0) bλ, bλ, λ) bλ, bλ, λ) 0) 0
11 To determine the stability of the equilibrium points, we will need to linearize the state model about them. The Jacobian is σ σ 0 Df = + λ x 3 x ) x x b For 0, 0, 0): The characteristic equation of Eqn. is: Df 0,0,0) = + λ 0 σ σ b ) s 3 + b + σ + ) s + σ b λ σ + b) s bλ σ = 0 3) Since the constant coefficient is negative we can conclude that the equilibrium point at 0, 0, 0) is always unstable, which comes from the necessary condition for Routh stability. For bλ, bλ, λ): Df bλ, bλ,λ) = s + σ σ 0 s + b λ b λ b λ s + b 4) The characteristic equation of Eqn. 4 is: s 3 + b + σ + ) s + bλ + σ b + b) s + bλ σ = 0 ) To determine the stability we can construct the Routh array. bλ + σ + ) b + σ + bλσ b + σ + )bλ + bσ + b) bλσ b + σ + 6) bλσ In order for the system to be stable, all of the coefficients in the first column have to be positive. Considering the third coefficient from the top: After Simplifying: Now factor the σ terms: Finally, solving for λ: b + σ + )bλ + bσ + b) bλσ b + σ + > 0 7) b + σ)λ + σ + + b)σ + b + ) > 0 8) b + σ)λ + σ + )σ + b + ) > 0 9) λ > σ + )σ + b + ) σ b ) 60) Also, since we know that λ is positive, σ b ) > 0 which simplifies to σ > b + 6)
12 Figure : Problem 7 Part b) Figure shows the projection on the x, x ) plane of a trajectory of the Lorentz equations with σ = 0, λ = 4, and b = with initial condition of,, ). Part c) Differentiating the function V x, x, x 3 ) = x + x + x 3 a) ) : Now plugging in the state equations and simplifying gives us: From this point there are two ways to solve it. Method V = x ẋ + x ẋ + x 3 a)ẋ 3 6) V = σx x bx 3 x 3 a) 63) In this method, we will try to show that V αv + β 64) where α and β are positive constants. Completing the square for x 3 in Eqn. 63) results in: V = σx + x + b x 3 ax 3 + a ) + ) bx 3 ba) V = σx + x + b x 3 a) + bx 3 + 6) ba Since we are trying to bound V we can drop the bx 3 since it will always be positive and will only make V smaller, therefore Eqn. 6 becomes: V σx + x + b ) x 3 a) + ba 66)
13 From this we can easily verify that β = ba 67) In order to determine α, we need to compare the coefficients in Eqn. 66 and V x). Through this comparison we find that α = minσ,, b) 68) Since Eqn. 64 holds, we can then find the region in which V is always less than zero. αv + β < 0 V > β α 69) Therefore, V < 0 when V > β α. No trajectory can leave the region given by V β α and any trajectory starting outside of this region must either enter it or asymptotically approach its boundary. Method This method is more of a geometric argument. By completing the square for x 3 in Eqn. 63 we obtain, V = σx + x + b x 3 a ) ) + 4 ba 70) Realizing that the right hand side of Eqn. 70 represents an ellipse in which we will now consider the region bounded by: which is an ellipse center at 0, 0, a ). On this surface, V = 0. σx + x + b x 3 a ) = 4 ba 7) Outside the region, σx + x + b x 3 a ) > 4 ba, we can show that V < 0. Therefore, we can conclude that all trajectories are bounded within the region: σx + x + b x 3 a ) 4 ba 7) 3
Outline. 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 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 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 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 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 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 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 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 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 information2D-Volterra-Lotka Modeling For 2 Species
Majalat Al-Ulum Al-Insaniya wat - Tatbiqiya 2D-Volterra-Lotka Modeling For 2 Species Alhashmi Darah 1 University of Almergeb Department of Mathematics Faculty of Science Zliten Libya. Abstract The purpose
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 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 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 informationLyapunov functions and stability problems
Lyapunov functions and stability problems Gunnar Söderbacka, Workshop Ghana, 29.5-10.5, 2013 1 Introduction In these notes we explain the power of Lyapunov functions in determining stability of equilibria
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 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 informationPhase Portraits of Nonlinear Differential Equations
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
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 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 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 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 informationNon-Linear Dynamics Homework Solutions Week 6
Non-Linear Dynamics Homework Solutions Week 6 Chris Small March 6, 2007 Please email me at smachr09@evergreen.edu with any questions or concerns reguarding these solutions. 6.8.3 Locate annd find the index
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 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 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 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 informationTWO DIMENSIONAL FLOWS. Lecture 5: Limit Cycles and Bifurcations
TWO DIMENSIONAL FLOWS Lecture 5: Limit Cycles and Bifurcations 5. Limit cycles A limit cycle is an isolated closed trajectory [ isolated means that neighbouring trajectories are not closed] Fig. 5.1.1
More 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 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 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 informationSection 5. Graphing Systems
Section 5. Graphing Systems 5A. The Phase Plane 5A-1. Find the critical points of each of the following non-linear autonomous systems. x = x 2 y 2 x = 1 x + y a) b) y = x xy y = y + 2x 2 5A-2. Write each
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 informationThe Big Picture. Discuss Examples of unpredictability. Odds, Stanisław Lem, The New Yorker (1974) Chaos, Scientific American (1986)
The Big Picture Discuss Examples of unpredictability Odds, Stanisław Lem, The New Yorker (1974) Chaos, Scientific American (1986) Lecture 2: Natural Computation & Self-Organization, Physics 256A (Winter
More information1. (i) Determine how many periodic orbits and equilibria can be born at the bifurcations of the zero equilibrium of the following system:
1. (i) Determine how many periodic orbits and equilibria can be born at the bifurcations of the zero equilibrium of the following system: ẋ = y x 2, ẏ = z + xy, ż = y z + x 2 xy + y 2 + z 2 x 4. (ii) Determine
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 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 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 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 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 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 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 informationContinuous time population models
Continuous time population models Jaap van der Meer jaap.van.der.meer@nioz.nl Abstract Many simple theoretical population models in continuous time relate the rate of change of the size of two populations
More informationDynamical Systems: Ecological Modeling
Dynamical Systems: Ecological Modeling G Söderbacka Abstract Ecological modeling is becoming increasingly more important for modern engineers. The mathematical language of dynamical systems has been applied
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 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 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 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 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 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 informationEG4321/EG7040. Nonlinear Control. Dr. Matt Turner
EG4321/EG7040 Nonlinear Control Dr. Matt Turner EG4321/EG7040 [An introduction to] Nonlinear Control Dr. Matt Turner EG4321/EG7040 [An introduction to] Nonlinear [System Analysis] and Control Dr. Matt
More information10 Back to planar nonlinear systems
10 Back to planar nonlinear sstems 10.1 Near the equilibria Recall that I started talking about the Lotka Volterra model as a motivation to stud sstems of two first order autonomous equations of the form
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 informationStability Implications of Bendixson s Criterion
Wilfrid Laurier University Scholars Commons @ Laurier Mathematics Faculty Publications Mathematics 1998 Stability Implications of Bendixson s Criterion C. Connell McCluskey Wilfrid Laurier University,
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 informationMath 3301 Homework Set Points ( ) ( ) I ll leave it to you to verify that the eigenvalues and eigenvectors for this matrix are, ( ) ( ) ( ) ( )
#7. ( pts) I ll leave it to you to verify that the eigenvalues and eigenvectors for this matrix are, λ 5 λ 7 t t ce The general solution is then : 5 7 c c c x( 0) c c 9 9 c+ c c t 5t 7 e + e A sketch of
More informationIntroduction to Applied Nonlinear Dynamical Systems and Chaos
Stephen Wiggins Introduction to Applied Nonlinear Dynamical Systems and Chaos Second Edition With 250 Figures 4jj Springer I Series Preface v L I Preface to the Second Edition vii Introduction 1 1 Equilibrium
More 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 informationGerardo Zavala. Math 388. Predator-Prey Models
Gerardo Zavala Math 388 Predator-Prey Models Spring 2013 1 History In the 1920s A. J. Lotka developed a mathematical model for the interaction between two species. The mathematician Vito Volterra worked
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 informationModelling and Mathematical Methods in Process and Chemical Engineering
Modelling and Mathematical Methods in Process and Chemical Engineering Solution Series 3 1. Population dynamics: Gendercide The system admits two steady states The Jacobi matrix is ẋ = (1 p)xy k 1 x ẏ
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 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 informationENGI 9420 Lecture Notes 4 - Stability Analysis Page Stability Analysis for Non-linear Ordinary Differential Equations
ENGI 940 Lecture Notes 4 - Stability Analysis Page 4.01 4. Stability Analysis for Non-linear Ordinary Differential Equations A pair of simultaneous first order homogeneous linear ordinary differential
More informationLinearization problem. The simplest example
Linear Systems Lecture 3 1 problem Consider a non-linear time-invariant system of the form ( ẋ(t f x(t u(t y(t g ( x(t u(t (1 such that x R n u R m y R p and Slide 1 A: f(xu f(xu g(xu and g(xu exist and
More informationContinuous Threshold Policy Harvesting in Predator-Prey Models
Continuous Threshold Policy Harvesting in Predator-Prey Models Jon Bohn and Kaitlin Speer Department of Mathematics, University of Wisconsin - Madison Department of Mathematics, Baylor University July
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 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 informationME 680- Spring Geometrical Analysis of 1-D Dynamical Systems
ME 680- Spring 2014 Geometrical Analysis of 1-D Dynamical Systems 1 Geometrical Analysis of 1-D Dynamical Systems Logistic equation: n = rn(1 n) velocity function Equilibria or fied points : initial conditions
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 informationMath 216 First Midterm 19 October, 2017
Math 6 First Midterm 9 October, 7 This sample exam is provided to serve as one component of your studying for this exam in this course. Please note that it is not guaranteed to cover the material that
More informationSection 5.4 (Systems of Linear Differential Equation); 9.5 Eigenvalues and Eigenvectors, cont d
Section 5.4 (Systems of Linear Differential Equation); 9.5 Eigenvalues and Eigenvectors, cont d July 6, 2009 Today s Session Today s Session A Summary of This Session: Today s Session A Summary of This
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 informationB5.6 Nonlinear Systems
B5.6 Nonlinear Systems 1. Linear systems Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Linear systems 1.1 Differential Equations 1.2 Linear flows 1.3 Linear maps
More informationA Novel Three Dimension Autonomous Chaotic System with a Quadratic Exponential Nonlinear Term
ETASR - Engineering, Technology & Applied Science Research Vol., o.,, 9-5 9 A Novel Three Dimension Autonomous Chaotic System with a Quadratic Exponential Nonlinear Term Fei Yu College of Information Science
More informationPhase portraits in two dimensions
Phase portraits in two dimensions 8.3, Spring, 999 It [ is convenient to represent the solutions to an autonomous system x = f( x) (where x x = ) by means of a phase portrait. The x, y plane is called
More informationLecture 5. Outline: Limit Cycles. Definition and examples How to rule out limit cycles. Poincare-Bendixson theorem Hopf bifurcations Poincare maps
Lecture 5 Outline: Limit Cycles Definition and examples How to rule out limit cycles Gradient systems Liapunov functions Dulacs criterion Poincare-Bendixson theorem Hopf bifurcations Poincare maps Limit
More informationThe Liapunov Method for Determining Stability (DRAFT)
44 The Liapunov Method for Determining Stability (DRAFT) 44.1 The Liapunov Method, Naively Developed In the last chapter, we discussed describing trajectories of a 2 2 autonomous system x = F(x) as level
More informationEntrance Exam, Differential Equations April, (Solve exactly 6 out of the 8 problems) y + 2y + y cos(x 2 y) = 0, y(0) = 2, y (0) = 4.
Entrance Exam, Differential Equations April, 7 (Solve exactly 6 out of the 8 problems). Consider the following initial value problem: { y + y + y cos(x y) =, y() = y. Find all the values y such that the
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 information3.5 Competition Models: Principle of Competitive Exclusion
94 3. Models for Interacting Populations different dimensional parameter changes. For example, doubling the carrying capacity K is exactly equivalent to halving the predator response parameter D. The dimensionless
More informationCHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD
CHALMERS, GÖTEBORGS UNIVERSITET EXAM for DYNAMICAL SYSTEMS COURSE CODES: TIF 155, FIM770GU, PhD Time: Place: Teachers: Allowed material: Not allowed: January 14, 2019, at 08 30 12 30 Johanneberg Kristian
More informationTHE SEPARATRIX FOR A SECOND ORDER ORDINARY DIFFERENTIAL EQUATION OR A 2 2 SYSTEM OF FIRST ORDER ODE WHICH ALLOWS A PHASE PLANE QUANTITATIVE ANALYSIS
THE SEPARATRIX FOR A SECOND ORDER ORDINARY DIFFERENTIAL EQUATION OR A SYSTEM OF FIRST ORDER ODE WHICH ALLOWS A PHASE PLANE QUANTITATIVE ANALYSIS Maria P. Skhosana and Stephan V. Joubert, Tshwane University
More information14 Periodic phenomena in nature and limit cycles
14 Periodic phenomena in nature and limit cycles 14.1 Periodic phenomena in nature As it was discussed while I talked about the Lotka Volterra model, a great deal of natural phenomena show periodic behavior.
More 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 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 informationProf. Krstic Nonlinear Systems MAE281A Homework set 1 Linearization & phase portrait
Prof. Krstic Nonlinear Systems MAE28A Homework set Linearization & phase portrait. For each of the following systems, find all equilibrium points and determine the type of each isolated equilibrium. Use
More informationProblem set 6 Math 207A, Fall 2011 Solutions. 1. A two-dimensional gradient system has the form
Problem set 6 Math 207A, Fall 2011 s 1 A two-dimensional gradient sstem has the form x t = W (x,, x t = W (x, where W (x, is a given function (a If W is a quadratic function W (x, = 1 2 ax2 + bx + 1 2
More informationBIFURCATION PHENOMENA Lecture 4: Bifurcations in n-dimensional ODEs
BIFURCATION PHENOMENA Lecture 4: Bifurcations in n-dimensional ODEs Yuri A. Kuznetsov August, 2010 Contents 1. Solutions and orbits: equilibria cycles connecting orbits compact invariant manifolds strange
More informationẋ = 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 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 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 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 informationNon-Linear Models. Non-Linear Models Cont d
Focus on more sophistiated interaction models between systems. These lead to non-linear, rather than linear, DEs; often not soluble exactly in analytical form so use Phase-Plane Analysis. This is a method
More informationStable Coexistence of a Predator-Prey model with a Scavenger
ANTON DE KOM UNIVERSITEIT VAN SURINAME INSTITUTE FOR GRADUATE STUDIES AND RESEARCH Stable Coexistence of a Predator-Prey model with a Scavenger Thesis in partial fulfillment of the requirements for the
More informationDynamical Systems and Chaos Part I: Theoretical Techniques. Lecture 4: Discrete systems + Chaos. Ilya Potapov Mathematics Department, TUT Room TD325
Dynamical Systems and Chaos Part I: Theoretical Techniques Lecture 4: Discrete systems + Chaos Ilya Potapov Mathematics Department, TUT Room TD325 Discrete maps x n+1 = f(x n ) Discrete time steps. x 0
More informationInteractions. Yuan Gao. Spring Applied Mathematics University of Washington
Interactions Yuan Gao Applied Mathematics University of Washington yuangao@uw.edu Spring 2015 1 / 27 Nonlinear System Consider the following coupled ODEs: dx = f (x, y). dt dy = g(x, y). dt In general,
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 informationDynamics of a Population Model Controlling the Spread of Plague in Prairie Dogs
Dynamics of a opulation Model Controlling the Spread of lague in rairie Dogs Catalin Georgescu The University of South Dakota Department of Mathematical Sciences 414 East Clark Street, Vermillion, SD USA
More informationHandout 2: Invariant Sets and Stability
Engineering Tripos Part IIB Nonlinear Systems and Control Module 4F2 1 Invariant Sets Handout 2: Invariant Sets and Stability Consider again the autonomous dynamical system ẋ = f(x), x() = x (1) with state
More informationMath 5490 November 5, 2014
Math 549 November 5, 214 Topics in Applied Mathematics: Introduction to the Mathematics of Climate Mondays and Wednesdays 2:3 3:45 http://www.math.umn.edu/~mcgehee/teaching/math549-214-2fall/ Streaming
More information