MCE693/793: Analysis and Control of Nonlinear Systems
|
|
- Catherine Greene
- 5 years ago
- Views:
Transcription
1 MCE693/793: Analysis and Control of Nonlinear Systems Systems of Differential Equations Phase Plane Analysis Hanz Richter Mechanical Engineering Department Cleveland State University
2 Systems of Nonlinear Differential Equations We often work with systems in the general form: ẋ = f(t,x) where solutions x(t) take values in R n and f(t,x) is a function defined on a subset of R R n. Whether the above systems admits one, none or multiple solutions given an initial point x(t 0 ) = x 0, is a fundamental question. Can we find a solution to How about ẋ = sign (x)? ẋ = x 1/3? (from Khalil, Sect. 2.2). We need a few mathematical concepts to answer existence and uniqueness of solutions. 2 / 36
3 Norms in Euclidean Space The n dimensional Euclidean space (over the real field) is formed by all vectors x = [x 1,x 2,...x n ] T where x i R for all i = 1,2...n. The norm x of a vector is a real-valued function satisfying four axioms for all x,y R n : 1. x 0 2. x = 0 x = 0 3. αx = α x for all α R 4. x+y x + y We are familiar with the Euclidean norm since childhood: x = n i=1 x2 i. But other norms are commonly used. Show that x = max i { x i } is a valid norm. How to extend the idea of norm to matrices will be discussed later 3 / 36
4 Continuity and Uniform Continuity of Functions We consider functions in Euclidean space, f : R n R m. We look at two equivalent definitions: 1. f is continuous at a point x if the sequence f(x k ) converges to f(x) for any sequence x k which converges to x. Succintly, f(x k ) x as x k x. 2. Given any ǫ > 0, we can find δ(ǫ,x) > 0 such that x y < δ(x,ǫ) f(x) f(y) < ǫ When δ(x,ǫ) is actually independent of x in the second definition, we speak of uniform continuity. Important: Uniform continuity obviously implies continuity, but the reverse is also true when f : S R m, where S R n is compact (closed and bounded). 4 / 36
5 Continuity in Some Variables and Piecewise Continuity The function f(x,y) = x sign (y) yx 2 is continuous in x but not in y. Specifically, continuity in a subset of variables is evaluated by regarding the remaining variables as constants. The function round(x) for x R is not continuous, but it is piecewise-continuous. That is, its domain can be partitioned into regions where the function is continuous. In this case, round(x) is continuous on any interval of the form: (k+0.5,k+1.5), k = 0,±1,±2,± / 36
6 Lipschitz Condition, Local Existence and Uniqueness Theorem (Khalil 2.2) Let f(t,x) be piecewise continuous in t and suppose the Lipschitz condition f(t,x) f(t,y) L x y holds for some constant L > 0, for all x,y in some closed ball of radius r centered at x 0 and for all t in [t,t 0 ]. Then the differential equation ẋ = f(t,x) with initial condition x(t 0 ) = x 0 has a unique solution in an interval [t 0,t 0 +δ], for some δ > 0. Example: for f(x) = x 2, use the above Theorem to analyze existence and uniqueness. Can the solution be continued indefinitely for any initial conditions? 6 / 36
7 Global Existence and Uniqueness The previous example shows that solutions may be short-lived unless we place additional conditions. The Lipschitz condition will have to hold not only locally, but globally in x and in an extended set for t. Boundedness of f(t,x 0 ) will also be required. Theorem (Khalil 2.3) Let f(t,x) be piecewise continuous in t and suppose the conditions f(t,x) f(t,y) L x y f(t,x 0 ) h hold for all x,y R n, all t [t 0,t 1 ] and some h > 0. Then the differential equation ẋ = f(t,x) with initial condition x(t 0 ) = x 0 has a unique solution in the interval [t 0,t 1 ] 7 / 36
8 Constant Solutions (Equilibrium Points) Equilibrium points are solutions of the form x(t) = c, for some constant c, valid for all t. They reflect operating setpoints in engineering systems and must be studied carefully. Nonlinear systems may have empty, finite or infinite equilibrium sets. The last class may consist of a continuum of equilibria or countless isolated points, or a mix. Finding and characterizing the cardinality of equilibrium sets involves the solution of the nonlinear equations f(t,x) = 0 We point out that most of the time we use autonomous systems (where t does not appear in f), but it is possible to find equilibria in non-autonomous systems. For example: ẋ = f(t,x) = cos(t)x has x = 0 as the single equilibrium point, regardless of t. But the behavior of the solution away from zero is influenced by the cos(t) factor. 8 / 36
9 Examples: Equilibrium Finding Whenever an analytical solution to f(x) = 0 is not possible, we must use numerical tools. Matlab s fsolve is useful. Find all equilibrium points for the tunnel diode example in Khalil (1.2) and verify the solution provided there using fsolve. Find all equilibrium points for the second-order system ẋ = round(y) ẏ = sin(x) and sketch them in the phase plane (x-y plane) 9 / 36
10 Control Systems and Equilibrium A control system (a system with an exogenous input) includes an additional argument in f: ẋ = f(t,x,u) where u may be specified as a function u(t,x) with values in R m. Although is may be possible to find constant solutions when u is non constant (take f(x,u) = cos(u)x, then x = 0 is an equilibrium point regardless of how u varies), the accepted definition of equilibrium points for control systems involves setting u to some constant u 0. It is clear then that the equilibrium solutions depend on u 0. For equilibrium-finding purposes, u 0 is simply treated as a parameter of f. 10 / 36
11 Equilibrium Point Stability We know that equilibrium solutions are stationary, that is, the system remains at its initial condition x 0 indefinitely whenever x 0 is an equilibrium point. More importantly, we want to determine whether the system has the tendency to return tox 0 if a small initial condition offset is used, or if it moves away fromx 0. Stable Unstable 11 / 36
12 Equilibrium Point Stability In engineering we normally prefer systems that regulate themselves back to their original operating point when perturbed. When systems don t have that natural tendency (example: magnetic levitation), we must add feedback controls to change system behavior from unstable to stable. A precise definition of stability is pending. With Lyapunov s linearization method, we determine the stability of an equilibrium point of a nonlinear system by examining the stability of the linearized system about the origin. Except for one case (zero eigenvalues), the stability of the linearized system is the same as the local stability of the nonlinear system. 12 / 36
13 Review of Small-Signal (Jacobian) Linearization Let x 0 be an equilibrium point for ẋ = f(x,u) for u = u 0, and consider f sufficiently smooth to possess a Taylor series expansion at (x 0,u 0 ): f(x,u) f(x 0,u 0 )+A(x x 0 )+B(u u 0 )+H.O.T. where A = J x (x 0,u 0 ) and B = J u (x 0,u 0 ) are the Jacobians of f(x,u) relative to x and u evaluated at the equilibrium point. Define δ x = x x 0 and δ u = u u 0. Neglecting the higher-order terms we obtain a linear dynamic description of δ x with δ u as an input: δ x = Aδ x +Bδ u Note that the above is a good approximation provided: 1. δ x and δ u are small (solutions remain close to equilibrium values). Othewise the H.O.T. introduce significant differences. 2. f(x 0,u 0 ) = 0, that is, we linearize relative to an equilibrium point. It is a popular mistake to linearize at non-equilibrium points (trajectories) using the above formulas. Linearization methods relative to trajectories are available. 13 / 36
14 Equilibria in Linear Systems Consider the n-state linear system ẋ = Ax As we saw before, the equilibrium set is either {0} (when A is nonsingular) or a linear subspace of R n (when A is singular). We also know that the system is asymptotically stable whenever all eigenvalues of A have negative real parts. This implies non-singularity and uniqueness of equilibrium. If at least one eigenvalue is zero, singularity follows and a dense equilibrium set arises. Second-order systems are highly relevant and allow graphical study through phase-plane analysis. In linear systems, we can completely characterize solution behavior by examining the eigenvalues and eigenvectors of A. In the second-order case, graphical representations are very useful. 14 / 36
15 2nd-Order Linear Systems: Distinct real eigenvalues When the two eigenvalues are real, distinct and nonzero, two linearly-independent real eigenvectors v 1 and v 2 can be found. Then A is diagonalizable by the transformation x = Tz, with T = [v 1 v 2 ]. In z (modal) coordinates, the system becomes decoupled: ż 1 = λ 1 z 1 ż 2 = λ 2 z 2 The solutions are single exponential functions, and satisfy z 2 = cz λ 2/λ 1 1 for some constant c dependent on initial conditions. The signs of λ 1 and λ 2 determine stability and the shape of the phase portrait. 15 / 36
16 Stable and Unstable Nodes in Modal Phase Space If both eigenvalues are negative, we have a stable node. If at least one is positive, we have an unstable node. In modal space, eigenvectors are (0,1) and (1,0). Recall that an initial condition belonging to the span of an eigenvector generates a trajectory along the eigenvector. 16 / 36
17 Stable and Unstable Nodes in Original Phase Space In phase space, eigenvectors are v 1 (slow) and v 2 (fast). All solutions are linear combinations of eigensolutions. Again, an initial condition belonging to the span of an eigenvector generates a trajectory along the eigenvector. 17 / 36
18 Saddle Point If the eigenvalues have opposite signs, the equilibrium point is a saddle. Examine solution behaviors near the stable and unstable eigenvectors. 18 / 36
19 Straight-Line Solutions If the eigenvalues are real and one of them is zero, there are straight line solutions. To see this, note that there is still a set of linearly-independent eigenvectors which diagonalize the system via x = Mz, with M = [v 1 v 2 ]. In these modal coordinates: ż 1 = 0 ż 2 = λ 2 z 2 The fact that phase trajectories are straight lines is obvious from these equations. Note that A is singular, so the equilibrium set is dense (the null space of A). 19 / 36
20 Straight-Line Solutions... If the eigenvalues are real and both of them are zero, A could be the zero matrix. In this case, there are no trajectories and the entire phase plane is made up of equilibrium points. But there are matrices whose nullspace has dimension 1 (nonzero A with λ 1 = λ 2 = 0). In this case, diagonalization leads to ż 1 = z 2 ż 2 = 0 Again, trajectories in modal space are horizontal lines. The direction of the trajectories depends on the initial condition for z / 36
21 Complex Eigenvalues: Stable Focus Points When the eigenvalues have the form a±bi with a < 0 and b 0, trajectories spiral into the equilibrium point. Clocwkise or counterclockwise? 21 / 36
22 Complex Eigenvalues: Focus Points When the eigenvalues have the form a±bi with a > 0 and b 0, trajectories spiral out of the equilibrium point. Clockwise or counterclockwise? 22 / 36
23 Imaginary Axis Eigenvalue(s): Center Points When the eigenvalues have the form bi with b 0, trajectories are periodic. The phase plane is partitioned by periodic orbits, that is, any two periodic orbits have empty intersection and the union of all periodic orbits is the entire plane. We will not call these orbits limit cycles, because there are no initial conditions outside an orbit which converge to it. 23 / 36
24 Nonlinear 2nd-Order Systems - Isoclines Method For nonlinear second-order systems, phase portraits can be generated by: 1. Computer simulation and graphics 2. Analytical solution and elimination of time from solutions 3. Method of isoclines 4. Other methods Methods 2 and 3 can be used to obtain insightful sketches and induce learning. We will focus on the method of isoclines. Given ẋ 1 = f 1 (x 1,x 2 ) ẋ 2 = f 2 (x 1,x 2 ) we can show that the slope of phase trajectories is given by dx 2 dx 1 = f 2(x 1,x 2 ) f 1 (x 1,x 2 ) 24 / 36
25 Isoclines Method To use the isoclines method, find and plot all equilibrium points. Then determine the locus of phase space points having particular values of slope that might be informative, such as,±1 and 0. From dx 2 = f 2(x 1,x 2 ) dx 1 f 1 (x 1,x 2 ) note that the slope at equilibrium points is undefined. Example: Let f 1 (x 1,x 2 ) = x 2 and f 2 (x 1,x 2 ) = x 1 +x 2. Sketch the phase portrait. Verify against known characteristics of linear systems. 25 / 36
26 Analytical Method In some cases, explicit solutions x 1 (t) and x 2 (t) can be found, given an initial point (x 10,x 20 ). If t can be solved for as a function of x 1 (or x 2 ), direct substitution into x 2 (t) (or x 1 (t)) results in an implicit relationship describing the phase trajectory passing through (x 10,x 20 ). Example: Use the analytical method for the previous example. 26 / 36
27 Guidelines for Computer-Generated Phase Portraits 1. Take advantage of symmetry (also valid for hand sketching, see S&L ) 2. Negative-time solutions: start close to the equilibrium point and compute backwards. Define τ = t and see what happens to the system relative to this new time. See Khalil, Use Matlab s ode45 or stiff methods as appropriate and compute solutions for short times. 4. See S&L 2.3 to estimate time spans needed for given trajectory lengths. 27 / 36
28 Limit Cycles - Background Definitions Consider the autonomous system ẋ = f(x) and let x = φ(x 0,t) be a solution with x(0) = x 0, assumed valid for t (, ). Define: The positive semiorbit through x 0 : γ + (x 0 ) = {φ(x 0,t) : 0 t < } The negative semiorbit through x 0 : γ (x 0 ) = {φ(x 0,t) : < t 0} 28 / 36
29 Limit Cycles - Background Definitions Definition of limit points and sets: If there is a sequence of times {t n } with t n as n such that the sequence {φ(t n,x)} p, we call p a positive limit point of the solution. If there is a sequence of times {t n } with t n as n such that the sequence {φ(t n,x)} p, we call p a negative limit point of the solution. The set of all positive limit points of φ(t,x 0 ) is the positive limit set of the solution. The set of all negative limit points of φ(t,x 0 ) is the negative limit set of the solution. 29 / 36
30 Limit Cycles - Background Definitions Definitions of periodic solutions and closed orbits: The solution φ(t,x 0 ) is nontrivial periodic if x 0 is not an equilibrium point and there is a constant T > 0 such that φ(t+t,x 0 ) = φ(x 0 ) for all t 0. The period is the smallest T that meets the above definition. If φ(t,x 0 ) is nontrivial periodic, the set {x R n : x = φ(t,x 0 ) for some t} is the periodic orbit or closed orbit associated with the solution. 30 / 36
31 Limit Cycles - Definition (Khalil 7.1) A limit cycle is a closed orbit γ such that γ is the positive limit set of some positive semiorbit γ + (x 0 ) or the negative limit set of some negative semiorbit γ (x 0 ), for some x 0 / γ. According to this definition, are linear harmonic oscillator solutions limit cycles? How about nonlinear pendulum oscillations? 31 / 36
32 The van der Pol Equation and its Limit Cycle The van der Pol equation ẍ µ(1 x 2 )ẋ+x = 0 was found by the eponymous Dutch researcher when investigating oscillations in circuits with vacuum tubes. We can interpret the equation as a mass-spring-damper system with nonlinear damping (which can have positive, negative or zero coefficient according to x). 3 van der Pol equation trajectories, µ= ẍ µ(1 x 2 )ẋ+x = 0 ẋ x 32 / 36
33 van der Pol Equation We use the van der Pol equation (µ 0) for an overview of equilibrium points and linearization.we also use it to demonstrate the isoclines method and computation of time from trajectories. 1. The origin is the only equilibrium point 2. The eigenvalues of the linearized system matrix A always have positive real parts for µ > Depending on µ, the origin may be a node or a focus, but always unstable for µ > Find values of µ to obtain each of the above 2 cases. 33 / 36
34 Limit Cycles in 2nd-Order Systems: 3 Theorems Knowing that the van der Pol equation is unstable at its equilibrium point, can we detect the presence of a limit cycle? We will consider a second-order system in the form of Eq. 2.1 in S&L: ẋ 1 = f 1 (x 1,x 2 ) ẋ 2 = f 2 (x 1,x 2 ) Poincaré s Theorem (index theorem): Let N be the total number of centers, foci and nodes enclosed by a limit cycle. Let S be the total number of enclosed saddle points. If a limit cycle exists, then N = S +1. The above implies that there can t be a limit cycle that doesn t contain an equilibrium point inside. The theorem provides a necessary condition for the existence of a limit cycle, not a sufficient condition. Make sure you have your logic straight. 34 / 36
35 Poincaré-Bendixson Theorem Statement from S&L: If a trajectory remains in a finite region of the phase plane for all t 0, then exactly one of the following holds: 1. The trajectory converges to an equilibrium point 2. The trajectory converges to a limit cycle 3. The trajectory is itself a limit cycle Examples 35 / 36
36 Poincaré-Bendixson and Bendixson Theorems Poincaré-Bendixson: Statement from Khalil: Let γ + be a bounded positive semiorbit and let L + be its positive limit set. If L + contains no equilibrium points, then it is a periodic orbit. Bendixson (S&L): If a limit cycle exists in a region Ω of the phase plane, then f 1 vanish and change sign somewhere in Ω. x 1 + f 2 x 2 must See the interesting proof with Stokes theorem in S&L. Examine the applicability of the three theorems to the van der Pol equation. Again, we have only a necessary condition. In general, limit cycle detection is a very hard problem. But when they don t exist, ruling them out is immediate using the above theorems. 36 / 36
Nonlinear Control Lecture 2:Phase Plane Analysis
Nonlinear Control Lecture 2:Phase Plane Analysis Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Fall 2010 r. Farzaneh Abdollahi Nonlinear Control Lecture 2 1/53
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 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 # 1 Introduction. Nonlinear Control
Nonlinear Control Lecture # 1 Introduction Nonlinear State Model ẋ 1 = f 1 (t,x 1,...,x n,u 1,...,u m ) ẋ 2 = f 2 (t,x 1,...,x n,u 1,...,u m ).. ẋ n = f n (t,x 1,...,x n,u 1,...,u m ) ẋ i denotes the derivative
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 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 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 informationLecture 4. Chapter 4: Lyapunov Stability. Eugenio Schuster. Mechanical Engineering and Mechanics Lehigh University.
Lecture 4 Chapter 4: Lyapunov Stability Eugenio Schuster schuster@lehigh.edu Mechanical Engineering and Mechanics Lehigh University Lecture 4 p. 1/86 Autonomous Systems Consider the autonomous system ẋ
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 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 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 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 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 informationMCE693/793: Analysis and Control of Nonlinear Systems
MCE693/793: Analysis and Control of Nonlinear Systems Lyapunov Stability - I Hanz Richter Mechanical Engineering Department Cleveland State University Definition of Stability - Lyapunov Sense Lyapunov
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 informationExamples include: (a) the Lorenz system for climate and weather modeling (b) the Hodgkin-Huxley system for neuron modeling
1 Introduction Many natural processes can be viewed as dynamical systems, where the system is represented by a set of state variables and its evolution governed by a set of differential equations. Examples
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 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 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 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 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 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 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 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 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 informationPoincaré Map, Floquet Theory, and Stability of Periodic Orbits
Poincaré Map, Floquet Theory, and Stability of Periodic Orbits CDS140A Lecturer: W.S. Koon Fall, 2006 1 Poincaré Maps Definition (Poincaré Map): Consider ẋ = f(x) with periodic solution x(t). Construct
More 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 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 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 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 informationMCE693/793: Analysis and Control of Nonlinear Systems
MCE693/793: Analysis and Control of Nonlinear Systems Input-Output and Input-State Linearization Zero Dynamics of Nonlinear Systems Hanz Richter Mechanical Engineering Department Cleveland State University
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 informationThe goal of this chapter is to study linear systems of ordinary differential equations: dt,..., dx ) T
1 1 Linear Systems The goal of this chapter is to study linear systems of ordinary differential equations: ẋ = Ax, x(0) = x 0, (1) where x R n, A is an n n matrix and ẋ = dx ( dt = dx1 dt,..., dx ) T n.
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 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 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 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 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 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 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 informationApplied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.
Printed Name: Signature: Applied Math Qualifying Exam 11 October 2014 Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. 2 Part 1 (1) Let Ω be an open subset of R
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 informationEE363 homework 7 solutions
EE363 Prof. S. Boyd EE363 homework 7 solutions 1. Gain margin for a linear quadratic regulator. Let K be the optimal state feedback gain for the LQR problem with system ẋ = Ax + Bu, state cost matrix Q,
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 informationLMI Methods in Optimal and Robust Control
LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 15: Nonlinear Systems and Lyapunov Functions Overview Our next goal is to extend LMI s and optimization to nonlinear
More informationMA 527 first midterm review problems Hopefully final version as of October 2nd
MA 57 first midterm review problems Hopefully final version as of October nd The first midterm will be on Wednesday, October 4th, from 8 to 9 pm, in MTHW 0. It will cover all the material from the classes
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 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 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 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 informationAdditional Homework Problems
Additional Homework Problems These problems supplement the ones assigned from the text. Use complete sentences whenever appropriate. Use mathematical terms appropriately. 1. What is the order of a differential
More information2.3. VECTOR SPACES 25
2.3. VECTOR SPACES 25 2.3 Vector Spaces MATH 294 FALL 982 PRELIM # 3a 2.3. Let C[, ] denote the space of continuous functions defined on the interval [,] (i.e. f(x) is a member of C[, ] if f(x) is continuous
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 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 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 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 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 information(x k ) sequence in F, lim x k = x x F. If F : R n R is a function, level sets and sublevel sets of F are any sets of the form (respectively);
STABILITY OF EQUILIBRIA AND LIAPUNOV FUNCTIONS. By topological properties in general we mean qualitative geometric properties (of subsets of R n or of functions in R n ), that is, those that don t depend
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 08, 2018, at 08 30 12 30 Johanneberg Kristian
More informationMath 302 Outcome Statements Winter 2013
Math 302 Outcome Statements Winter 2013 1 Rectangular Space Coordinates; Vectors in the Three-Dimensional Space (a) Cartesian coordinates of a point (b) sphere (c) symmetry about a point, a line, and a
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 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 informationCharacterization of the stability boundary of nonlinear autonomous dynamical systems in the presence of a saddle-node equilibrium point of type 0
Anais do CNMAC v.2 ISSN 1984-82X Characterization of the stability boundary of nonlinear autonomous dynamical systems in the presence of a saddle-node equilibrium point of type Fabíolo M. Amaral Departamento
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 information27. Topological classification of complex linear foliations
27. Topological classification of complex linear foliations 545 H. Find the expression of the corresponding element [Γ ε ] H 1 (L ε, Z) through [Γ 1 ε], [Γ 2 ε], [δ ε ]. Problem 26.24. Prove that for any
More 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 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 informationAutonomous system = system without inputs
Autonomous system = system without inputs State space representation B(A,C) = {y there is x, such that σx = Ax, y = Cx } x is the state, n := dim(x) is the state dimension, y is the output Polynomial representation
More informationStability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games
Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,
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 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 informationConstruction of Lyapunov functions by validated computation
Construction of Lyapunov functions by validated computation Nobito Yamamoto 1, Kaname Matsue 2, and Tomohiro Hiwaki 1 1 The University of Electro-Communications, Tokyo, Japan yamamoto@im.uec.ac.jp 2 The
More informationELEC 3035, Lecture 3: Autonomous systems Ivan Markovsky
ELEC 3035, Lecture 3: Autonomous systems Ivan Markovsky Equilibrium points and linearization Eigenvalue decomposition and modal form State transition matrix and matrix exponential Stability ELEC 3035 (Part
More informationODE Homework 1. Due Wed. 19 August 2009; At the beginning of the class
ODE Homework Due Wed. 9 August 2009; At the beginning of the class. (a) Solve Lẏ + Ry = E sin(ωt) with y(0) = k () L, R, E, ω are positive constants. (b) What is the limit of the solution as ω 0? (c) Is
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: August 22, 2018, at 08 30 12 30 Johanneberg Jan Meibohm,
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 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 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 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 informationA Study of the Van der Pol Equation
A Study of the Van der Pol Equation Kai Zhe Tan, s1465711 September 16, 2016 Abstract The Van der Pol equation is famous for modelling biological systems as well as being a good model to study its multiple
More informationNonlinear System Analysis
Nonlinear System Analysis Lyapunov Based Approach Lecture 4 Module 1 Dr. Laxmidhar Behera Department of Electrical Engineering, Indian Institute of Technology, Kanpur. January 4, 2003 Intelligent Control
More informationApplied Dynamical Systems
Applied Dynamical Systems Recommended Reading: (1) Morris W. Hirsch, Stephen Smale, and Robert L. Devaney. Differential equations, dynamical systems, and an introduction to chaos. Elsevier/Academic Press,
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 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 informationLecture Notes 6: Dynamic Equations Part C: Linear Difference Equation Systems
University of Warwick, EC9A0 Maths for Economists Peter J. Hammond 1 of 45 Lecture Notes 6: Dynamic Equations Part C: Linear Difference Equation Systems Peter J. Hammond latest revision 2017 September
More information5.2.2 Planar Andronov-Hopf bifurcation
138 CHAPTER 5. LOCAL BIFURCATION THEORY 5.. Planar Andronov-Hopf bifurcation What happens if a planar system has an equilibrium x = x 0 at some parameter value α = α 0 with eigenvalues λ 1, = ±iω 0, ω
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 informationLecture 1: Systems of linear equations and their solutions
Lecture 1: Systems of linear equations and their solutions Course overview Topics to be covered this semester: Systems of linear equations and Gaussian elimination: Solving linear equations and applications
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 informationDIAGONALIZABLE LINEAR SYSTEMS AND STABILITY. 1. Algebraic facts. We first recall two descriptions of matrix multiplication.
DIAGONALIZABLE LINEAR SYSTEMS AND STABILITY. 1. Algebraic facts. We first recall two descriptions of matrix multiplication. Let A be n n, P be n r, given by its columns: P = [v 1 v 2... v r ], where the
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 informationMATH 44041/64041 Applied Dynamical Systems, 2018
MATH 4441/6441 Applied Dynamical Systems, 218 Answers to Assessed Coursework 1 1 Method based on series expansion of the matrix exponential Let the coefficient matrix of α the linear system be A, then
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 informationNonconstant Coefficients
Chapter 7 Nonconstant Coefficients We return to second-order linear ODEs, but with nonconstant coefficients. That is, we consider (7.1) y + p(t)y + q(t)y = 0, with not both p(t) and q(t) constant. The
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 informationNonlinear Systems Theory
Nonlinear Systems Theory Matthew M. Peet Arizona State University Lecture 2: Nonlinear Systems Theory Overview Our next goal is to extend LMI s and optimization to nonlinear systems analysis. Today we
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 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 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 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 information