MCE693/793: Analysis and Control of Nonlinear Systems
|
|
- Frank Atkinson
- 5 years ago
- Views:
Transcription
1 MCE693/793: Analysis and Control of Nonlinear Systems Lyapunov Stability - I Hanz Richter Mechanical Engineering Department Cleveland State University
2 Definition of Stability - Lyapunov Sense Lyapunov stability is related to the tendency of system trajectories to remain in prescribed regions surrounding an equilibrium point. This includes the possibility (but does not require) convergence to the equilibrium point, where the stability is asymptotic. The harmonic oscillator as an autonomous LTI system ẋ 1 = x 2 ẋ 2 = w 2 x 1 gives circular orbits passing through any initial condition (x 10,x 20 ). Each circle of radius x x2 20 determines a confinement region Ω, so that all initial conditions within Ω result in trajectories remaining in Ω. But no nonzero initial condition converges to the origin. The origin of the autonomous system is stable in the sense of Lyapunov, but not asymptotically stable. Warning: The forced undamped oscillator is unstable in the bounded input - bounded output sense. 2 / 33
3 Mathematical Definition of Lyapunov Stability Consider the autonomous nonlinear system ẋ = f(x) and let x eq be an equilibrium point. The system is stable about the equilibrium point if the following property holds: Fo any R > 0 there always exists r > 0 so that if x 0 x eq < r, then x(t) x eq < R for all t > 0. R defines a desired confinement region, while r defines the neighborhood of x eq where x 0 must belong so that x(t) does not exit the confinement region. R x eq r x 0 x(t) 3 / 33
4 Asymptotic Stability Instead of saying the nonlinear system is stable relative to equilibrium point x 0, we can just say x 0 is stable in context. The equilibrium point x 0 is asymptotically stable if it is stable and also there exists r > 0 such that the following implication holds: x(0) x 0 < r = x(t) x 0 as t That is, x 0 is the only positive limit point of all solutions with initial conditions in the ball B(x 0,r). This ball is a region of attraction of x 0. The maximal region of attraction does not have to be a ball. In general, regions of attraction have more complex shapes. 4 / 33
5 Exponential Stability The equilibrium point x 0 is exponentially stable if α,λ such that x(t) x 0 α x(0) x 0 e λt Checking whether a system is stable is in general impossible by using only the definition (a large number of initial conditions would have to be tested numerically, or the explicit solution could be used for very simple systems). Examine the example in S&L: ẋ = (1+sin 2 x)x Here exponential stability (and therefore asymptotic stability) can be established with the solution. 5 / 33
6 Local and Global (in the large) Stability The radius r involved in the definitions of stability and asymptotic stability provide a domain of attraction. When r is finite, stability is local. If r can be extended to infinity (if in the definition, r can always be chosen equal to R), stability is global. Immediate observations: If a system has multiple equilibrium points, can any one of them be globally stable? Asymptotic stability implies stability trivially, by definition, but: Why is stability included in the definition of asymptotic stability in the first place? Does the convergence condition alone imply stability? As an exercise, we show that the system given by Vinograd (S&L Figure 3.5) is unstable, but still satisfies the convergence property. 6 / 33
7 Lyapunov s Linearization Method This method allows us to determine the stability of the nonlinear system about the equilibrium point on the basis of the linearized system. Simply: If the eigenvalues of the A matrix in the linearized system have negative real parts, the nonlinear system is asymptotically stable about the equilibrium point. If at least one eigenvalue of the A matrix in the linearized system has positive real part, the nonlinear system is unstable about the equilibrium point. If at least one eigenvalue ofthe A matrix in the linearized system has zero real part, the test is inconclusive. The linear approximation is insufficient to determine stability. However, methods exist to include higher-order terms (center manifold technique). 7 / 33
8 Positive-Definite Functions and Matrices The function f : R n R is positive semidefinite if f(0) = 0 and f(x) 0 for all x R n. If f(x) > 0 for all nonzero x R n, then f is positive definite. Example: f(x 1,x 2 ) = x 2 1 +x 2 2 +sin 2 (x 1 ) is positive-definite. If f(x) 0 (resp. f(x) > 0) in a subset A R n, we say that f is positive-semidefinite (resp. positive definite) in A. When f(x) is a quadratic function defined through a symmetric matrix P as f(x) = x T Px then we say that P is positive-semidefinite (resp. positive definite) when the associated function is so. Finally, if f is positive-semidefinite (resp. positive definite), we say that f is negative-semidefinite (resp. negative definite). If f changes sign, it is sign-indefinite. 8 / 33
9 Testing for Positive-Definite Matrices Suppose f(x) = x x x2 3 x 1x 2 +8x 2 x 3 5x 1 x 3. How can we tell if f is sign-definite? The symmetric matrix associated with f is P = (1) The eigenvalue test can be applied. If all eigenvalues are nonnegative, we have a p.s.d. matrix. If they are all positive, we have a p.d. matrix. The same idea applies for n.s.d. and n.d. If the eigenvalues have mixed signs, the matrix is sign-indefinite. 9 / 33
10 Lyapunov s Direct (Second) Method For simplicity, suppose the origin is an equilibrium point. Suppose we find a function V : R n R which has continuous first partial derivatives and is positive-definite in a region surrounding the origin. We call this function a candidate Lyapunov function. If the function V is negative-semidefinite in the same region, then the origin is a stable equilibrium point. If V is negative-definite, the origin is asymptotically stable. We use Lyapunov s method with a quadratic function to prove that the van der Pol equation with µ = 2 is stable relative to the origin. Is the proven stability global? Lyapunov s direct method is only a sufficient condition. Hence the difficulties in finding Lyapunov functions 10 / 33
11 Global Asymptotic Stability Stronger conditions on V and its derivative are required for global asymptotic stability. Suppose V(x) has continuous first partial derivatives and satisfies: 1. V is positive definite in R n 2. V is negative-definite in R n 3. V(x) as x (V is radially unbounded) then the origin is globally asymptotically stable. 11 / 33
12 Examples The function V(x) = x is positive-definite, but cannot be used with the above results because it does not have a partial derivative at zero. However, results are available for piecewise-continuous Lyapunov functions. The system ẋ = x 3 is G.A.S. at zero, but the linearization method is inconclusive. Systems of the form ẋ = c(x), where c(x) is a continuous odd function with c(0) = 0 are always G.A.S. at x = 0. All stable linear systems ẋ = Ax admit a quadratic Lyapunov function V(x) = 1 2 xt Px. For linear stable systems and V(x) = x T Px, V = x T (A T P +PA)x. Admissible P s can be found by solving Lyapunov s equation: with any Q > 0. A T P +PA = Q 12 / 33
13 Constructive Controller Design with Lyapunov Functions Many techniques for control design are based on the idea of requiring a given Lyapunov function candidate to be negative (semi)definite as a way to obtain a control law. This is commonly done in robust and adaptive control, as well as the backstepping and sliding mode techniques. As an introductory example, consider the second-order control system ẋ 1 = x 2 ẋ 2 = x 3 1 +x 2 +u Examine system stability about the origin without control. Then use a quadratic Lyapunov function to find a stabilizing control law. Is the stability global or asymptotic? 13 / 33
14 Results for semi-definite V Frequently, Lyapunov s direct method goes only as far as showing stability, even when the system is actually asymptotically stable. This can happen with systems as simple as the linear mass-spring-damper with a quadratic V. Invariant set theorems, including the LaSalle-Krasovskii s theorem (known as LaSalle s invariance principle, ) and Barbalat s lemma (1959) are useful to analyze long-term system behavior when the Lyapunov function derivative is only negative semi-definite. The LaSalle invariance principle is used to prove asymptotic stability when the Lyapunov function derivative is only negative semi-definite. Barbǎlat s result is used very often in adaptive control (and in the proofs of the invariant set theorems). 14 / 33
15 Invariant Sets Let a nonlinear system be defined by ẋ = f(x) A set I R n is invariant if the following is true: x(t 0 ) I x(t) I t > t 0 Example: In the damped pendulum with states θ and θ, let E represent the total energy of the pendulum. For any l > 0, the set I = {(θ, θ) : E(θ, θ) l} is invariant. What is the shape of the boundary of I? 15 / 33
16 Local Invariant Set Theorem Refer to S&L, Th Examine Fig Let a nonlinear system be defined by ẋ = f(x) with f : R n R n continuous and f(0) = 0. Let V(x) be a scalar function with continuous first partial derivatives. For some l > 0 define sets Ω l and R by Ω l = {x R n : V(x) < l and V(x) 0} R = {x Ω l : V(x) = 0} Let M be the largest invariant set contained in R. Then every trajectory x(t) with x(0) Ω l is such that x(t) x M as t for some point x M M. X is the largest invariant set contained in Y means: if Z Y and Z is invariant, then Z X 16 / 33
17 LaSalle s Invariance Principle This is a special case of the local invariant set theorem. Suppose the conditions of the local invariant set theorem hold, and in addition, suppose V(x) is positive-definite in a set Ω containing the origin, and that V(x) 0 in Ω. Define R = {x Ω : V(x) = 0} If R = {0} (a singleton), then 0 is asymptotically stable. Moreover, let Ω l be the largest connected subset of Ω that can be defined by Ω l = {x R n : V(x) < l and V(x) 0} for some l > 0. Then Ω l is a domain of attraction of the origin. 17 / 33
18 LaSalle s Principle... Note that if R is invariant, then M = R. To use LaSalle s result, show that the only possible trajectory belonging to R is the zero trajectory: x(t) = 0 for t 0. Classical example: pendulum with viscous damping. We do this in class. Example 3.13 in S&L is very illustrative. We add details in class. 18 / 33
19 Global Invariance Theorem If Ω l is extended to R n and radial unboundedness is required, the invariance theorem can be applied to solutions starting from anywhere in R n. Consider the autonomous system ẋ = f(x), with f continuous. Let V(x) be a scalar function with continuous first partial derivatives and assume V(x) 0 x R n and V(x) as x. Define R = {x : V(x) = 0} and let M be the largest invariant set in R. Then every trajectory x(t) satisfies x(t) x M as t for some point x M M. 19 / 33
20 Barbǎlat s Lemma A function f : R n R is square integrable if: 0 f T (t)qf(t)dt for some matrix Q = Q T > 0. A function g : R n R is uniformly continuous if for any ǫ > 0 there is a δ > 0 such that: x y < δ g(x) g(y) < ǫ for all x,y R n. Note that δ can be a function of ǫ but not on the points x and y (that would be plain continuity for g; uniform continuity is a stronger property). Barbǎlat s lemma: If f(t) is square integrable and df(t) dt is uniformly continuous then df(t) dt 0 as t. 20 / 33
21 Barbǎlat s Lemma in Lyapunov Theory Let a nonlinear system be defined by ẋ = f(x,ψ(t)) where ψ(t) is some input function. Suppose the origin is an equilibrium point (f(0,ψ) = 0) for any ψ. Suppose we found a function V(x,t) which is lower-bounded in a set D containing the origin, and that V(x, t) 0 in D. If V(x, t) is uniformly continuous with respect to time, then V(x,t) 0 as t. Note that uniform continuity of V can be satisfied by verifying that V is bounded. 21 / 33
22 Variants of Barbǎlat s Lemma Some alternatives have surfaced over the years. For instance, Tao, G. A simple alternative to the Barbǎlat Lemma, IEEE Trans. Aut. Ctrl. V42,N5., If f is square integrable and has a bounded derivative, then f itself converges to zero asymptotically. 22 / 33
23 Lyapunov Function Construction - Krasovskii s Method When considering a system ẋ = f(x) suspected to be stable, Lyapunov theory provides a powerful method to verify stability. Note that f may be the closed-loop system resulting of applying a new control design method. Now we want to prove that the design is indeed stabilizing. Krasovskii s method involves using V(x) = f T (x)f(x) as a candidate Lyapunov function. Krasovskii s Theorem: Suppose f(0) = 0 and let A(x) be the Jacobian matrix function. If F(x) = A(x)+A T (x) is negative definite in some neighborhood Ω of 0, then 0 is asymptotically stable and V(x) = f T (x)f(x) is a Lyapunov function. If Ω = R n and V(x) is radially unbounded, 0 is globally asymptotically stable. Examine example Additional example in class. 23 / 33
24 Generalized Krasovskii s Method It is very likely that F(x) will not be negative definite, or that checking this might be really hard for high-order systems. Some flexibility is gained by extending the Lyapunov function candidate to a general quadratic function. Suppose f(0) = 0 and let A(x) be the Jacobian matrix function. Suppose P = P T > 0 and Q = Q T > 0 such that F(x) = A T (x)p +PA(x)+Q is negative semidefinite in some deleted neighborhood Ω of 0. Then 0 is asymptotically stable and V(x) = f T (x)pf(x) is a Lyapunov function. If Ω = R n and V(x) is radially unbounded, 0 is globally asymptotically stable. Although more flexibility has been afforded, the search for P and Q will be computationally intensive. If we shrink Ω to just the origin and obtain F < 0, what can we conclude? 24 / 33
25 Vector Calculus Refresher: Potential Functions This material is necessary for full understanding of the Variable Gradient Method for constructing Lyapunov functions. Refer to Hass, Weir and Thomas, Multivariable University Calculus: Early Transcendentals, 2nd Ed., Addison-Wesley, Definition: If F : R n R n is a vector field defined on a set D and Φ such that Φ = F on D, Φ is called a potential function for F. Examples: Φ(h) = mg(h 0 h) is a potential function for F = mg. Φ(x) = 1 2 kx2 is a potential function for F = kx. Φ(x 1,x 2 ) = sin(x 1 )x 2 +1 is a potential function for F = [x 2 cos(x 1 ) sin(x 1 )] In physics, a force field F is conservative if it is the gradient of some potential. 25 / 33
26 Curl Condition for the Existence of a Potential Function If a vector function F is given, how do we know if it has a potential function? Theorem: Let F(x) = [F 1 (x),f 2 (x),...f n (x)] be a vector field defined on a simply connected domain (a closed curve can be continuously shrunk to a point without leaving the domain). Assume the F i have continuous first partial derivatives. Then Φ such that Φ = F if and only if for all i,j = 1,2...n. F i x j = F j x i For n = 3, the condition is simply curl(f) = 0, hence the name. In fluid flow, we call a velocity field F satisfying this condition irrotational. 26 / 33
27 Conservative Fields and Path Independence Fundamental Theorem of Line Integrals: Let Φ be a scalar function defined on a domain D and let C be a curve joining two points x,y D R n. Then Φ(y) Φ(x) = Φ(r).dr In physics, we evaluate the work done by a conservative force between two points by evaluating the change in kinetic energy, instead of computing a line integral. C 27 / 33
28 Variable Gradient Method The idea is to propose the gradient of a Lyapunov function instead of the function directly. In doing this, we have the opportunity to enforce negative-definiteness of V, but we then have to check positive-definiteness of V (the method is often unsuccessful, but worth trying). 1. Start by assuming a simple form for V, for instance ( V) i = n j=1 a ijx j, where the a ij could be variable. 2. Impose the curl conditions 3. Impose conditions on a ij so that V 0 in a domain D surrounding the equilibrium point that is as large as possible. 4. Integrate to find V. Here we must perform the line integrals, but path independence allows us to use a convenient rectangular path. 5. Check if V is positive definite, at least in a subset of D. Follow S&L Example 3.20 carefully. 28 / 33
29 Example Schultz G. and Gibson, J.E., The variable gradient method for generating Liapunov functions, Trans. AIEE, Pt II 81 pp , Consider ẋ 1 = x 2 ẋ 2 = x 3 1 x 2 29 / 33
30 Regions of Attraction and Zubov s Construction Method Often, a controller design results in non-global asymptotic stability. For safe operation, it is important to obtain an estimate of the region of attraction, so the system is not allowed to operated away from its guaranteed stability region. Zubov s Theorem: Consider the autonomous system ẋ = f(x), with f(0) = 0 and let Ω be a domain containing the origin. Suppose there exist a function V : Ω R and φ : R n R such that: 1. V is positive definite and continously differentiable in Ω 2. φ is positive definite and continuous in R n 3. 0 V(x) < 1 for all x Ω and V(x) = 1 as x approaches the boundary of Ω (or x if Ω is unbounded). 4. V = (V 1)φ. Then the origin is asymptotically stable and Ω is the (exact) region of attraction. 30 / 33
31 Example (Hahn) ẋ 1 = x 1 +2x 2 1x 2 ẋ 2 = x 2 Take φ(x 1,x 2 ) = x 2 1 +x2 2 and compute the region of attraction. We must assume an arbitrary form for φ(x), usually a sum of squares. When an exact solution for V is not possible, we can obtain an estimate by approximating the nonlinear system. 31 / 33
32 Zubov s Approximation Method See Zubov, V.I., Methods of A.M. Lyapunov and Their Application, P. Noordhoff LTD, Groningen 1964 for full details (requires fairly advanced math to follow) For instance if we can assume ẋ = f(x) Ax+g(x) where A is Hurwitz and g(x) has second and higher-order terms, then we assume V = V 2 +V V j, where the subindex indicates quadratic, cubic, etc. homogeneous functions (sums of monomials of the same degree). Equate and match coefficients. V = T V.f = (V 1)φ 32 / 33
33 Zubov s Approximation... Zubov s result: for every approximation ˆV = V 2 +V V m, the boundary of the true region of attraction lies in the region defined by α ˆV(x) β, where: Notes: α = max ˆV=0ˆV β = min ˆV=0ˆV 1. The exact contour is obtained in the limit, with ˆV an infinite power series. 2. Convergence is not monotonic. Using 4 terms could be worse than using 3. But in the long run, more terms lead to a better approximation. Example: Damped pendulum ẍ + sin(x) + kẋ = 0. We know stability of the origin can t be global because other equilibrium points exist. We use Zubov s method to find the region of attraction of the origin. 33 / 33
MCE/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 informationGeorgia Institute of Technology Nonlinear Controls Theory Primer ME 6402
Georgia Institute of Technology Nonlinear Controls Theory Primer ME 640 Ajeya Karajgikar April 6, 011 Definition Stability (Lyapunov): The equilibrium state x = 0 is said to be stable if, for any R > 0,
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 informationNonlinear Control. Nonlinear Control Lecture # 2 Stability of Equilibrium Points
Nonlinear Control Lecture # 2 Stability of Equilibrium Points Basic Concepts ẋ = f(x) f is locally Lipschitz over a domain D R n Suppose x D is an equilibrium point; that is, f( x) = 0 Characterize and
More 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 informationNonlinear Systems and Control Lecture # 12 Converse Lyapunov Functions & Time Varying Systems. p. 1/1
Nonlinear Systems and Control Lecture # 12 Converse Lyapunov Functions & Time Varying Systems p. 1/1 p. 2/1 Converse Lyapunov Theorem Exponential Stability Let x = 0 be an exponentially stable equilibrium
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 informationNonlinear Control. Nonlinear Control Lecture # 3 Stability of Equilibrium Points
Nonlinear Control Lecture # 3 Stability of Equilibrium Points The Invariance Principle Definitions Let x(t) be a solution of ẋ = f(x) A point p is a positive limit point of x(t) if there is a sequence
More informationNonlinear Control. Nonlinear Control Lecture # 8 Time Varying and Perturbed Systems
Nonlinear Control Lecture # 8 Time Varying and Perturbed Systems Time-varying Systems ẋ = f(t,x) f(t,x) is piecewise continuous in t and locally Lipschitz in x for all t 0 and all x D, (0 D). The origin
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 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 Lyapunov theory of stability
M.Kawski, APM 581 Diff Equns Intro to Lyapunov theory. November 15, 29 1 1 Lyapunov theory of stability Introduction. Lyapunov s second (or direct) method provides tools for studying (asymptotic) stability
More informationNonlinear systems. Lyapunov stability theory. G. Ferrari Trecate
Nonlinear systems Lyapunov stability theory G. Ferrari Trecate Dipartimento di Ingegneria Industriale e dell Informazione Università degli Studi di Pavia Advanced automation and control Ferrari Trecate
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 informationStability in the sense of Lyapunov
CHAPTER 5 Stability in the sense of Lyapunov Stability is one of the most important properties characterizing a system s qualitative behavior. There are a number of stability concepts used in the study
More informationNonlinear Control Lecture 4: Stability Analysis I
Nonlinear Control Lecture 4: Stability Analysis I Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Fall 2010 Farzaneh Abdollahi Nonlinear Control Lecture 4 1/70
More informationDynamical Systems & Lyapunov Stability
Dynamical Systems & Lyapunov Stability Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University Outline Ordinary Differential Equations Existence & uniqueness Continuous dependence
More information2 Lyapunov Stability. x(0) x 0 < δ x(t) x 0 < ɛ
1 2 Lyapunov Stability Whereas I/O stability is concerned with the effect of inputs on outputs, Lyapunov stability deals with unforced systems: ẋ = f(x, t) (1) where x R n, t R +, and f : R n R + R n.
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 informationLyapunov Stability Theory
Lyapunov Stability Theory Peter Al Hokayem and Eduardo Gallestey March 16, 2015 1 Introduction In this lecture we consider the stability of equilibrium points of autonomous nonlinear systems, both in continuous
More informationNonlinear Control Lecture 5: Stability Analysis II
Nonlinear Control Lecture 5: Stability Analysis II Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Fall 2010 Farzaneh Abdollahi Nonlinear Control Lecture 5 1/41
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 information3. Fundamentals of Lyapunov Theory
Applied Nonlinear Control Nguyen an ien -.. Fundamentals of Lyapunov heory he objective of this chapter is to present Lyapunov stability theorem and illustrate its use in the analysis and the design of
More informationLyapunov stability ORDINARY DIFFERENTIAL EQUATIONS
Lyapunov stability ORDINARY DIFFERENTIAL EQUATIONS An ordinary differential equation is a mathematical model of a continuous state continuous time system: X = < n state space f: < n! < n vector field (assigns
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 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 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 informationChapter III. Stability of Linear Systems
1 Chapter III Stability of Linear Systems 1. Stability and state transition matrix 2. Time-varying (non-autonomous) systems 3. Time-invariant systems 1 STABILITY AND STATE TRANSITION MATRIX 2 In this chapter,
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 informationConverse Lyapunov theorem and Input-to-State Stability
Converse Lyapunov theorem and Input-to-State Stability April 6, 2014 1 Converse Lyapunov theorem In the previous lecture, we have discussed few examples of nonlinear control systems and stability concepts
More informationControl of Robotic Manipulators
Control of Robotic Manipulators Set Point Control Technique 1: Joint PD Control Joint torque Joint position error Joint velocity error Why 0? Equivalent to adding a virtual spring and damper to the joints
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 informationDO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED.
EE 533 Homeworks Spring 07 Updated: Saturday, April 08, 07 DO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED. Some homework assignments refer to the textbooks: Slotine
More informationLecture 9 Nonlinear Control Design
Lecture 9 Nonlinear Control Design Exact-linearization Lyapunov-based design Lab 2 Adaptive control Sliding modes control Literature: [Khalil, ch.s 13, 14.1,14.2] and [Glad-Ljung,ch.17] Course Outline
More informationIntroduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems
p. 1/5 Introduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems p. 2/5 Time-varying Systems ẋ = f(t, x) f(t, x) is piecewise continuous in t and locally Lipschitz in x for all t
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 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 informationHybrid Systems Course Lyapunov stability
Hybrid Systems Course Lyapunov stability OUTLINE Focus: stability of an equilibrium point continuous systems decribed by ordinary differential equations (brief review) hybrid automata OUTLINE Focus: stability
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 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 informationStability of Nonlinear Systems An Introduction
Stability of Nonlinear Systems An Introduction Michael Baldea Department of Chemical Engineering The University of Texas at Austin April 3, 2012 The Concept of Stability Consider the generic nonlinear
More 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 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 informationHybrid Systems - Lecture n. 3 Lyapunov stability
OUTLINE Focus: stability of equilibrium point Hybrid Systems - Lecture n. 3 Lyapunov stability Maria Prandini DEI - Politecnico di Milano E-mail: prandini@elet.polimi.it continuous systems decribed by
More informationSTABILITY ANALYSIS OF DYNAMIC SYSTEMS
C. Melchiorri (DEI) Automatic Control & System Theory 1 AUTOMATIC CONTROL AND SYSTEM THEORY STABILITY ANALYSIS OF DYNAMIC SYSTEMS Claudio Melchiorri Dipartimento di Ingegneria dell Energia Elettrica e
More informationLyapunov Stability Analysis: Open Loop
Copyright F.L. Lewis 008 All rights reserved Updated: hursday, August 8, 008 Lyapunov Stability Analysis: Open Loop We know that the stability of linear time-invariant (LI) dynamical systems can be determined
More informationNonlinear Control. Nonlinear Control Lecture # 8 Time Varying and Perturbed Systems
Nonlinear Control Lecture # 8 Time Varying and Perturbed Systems Time-varying Systems ẋ = f(t,x) f(t,x) is piecewise continuous in t and locally Lipschitz in x for all t 0 and all x D, (0 D). The origin
More informationModeling and Analysis of Dynamic Systems
Modeling and Analysis of Dynamic Systems Dr. Guillaume Ducard Fall 2017 Institute for Dynamic Systems and Control ETH Zurich, Switzerland G. Ducard c 1 / 57 Outline 1 Lecture 13: Linear System - Stability
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 information1.7. Stability and attractors. Consider the autonomous differential equation. (7.1) ẋ = f(x),
1.7. Stability and attractors. Consider the autonomous differential equation (7.1) ẋ = f(x), where f C r (lr d, lr d ), r 1. For notation, for any x lr d, c lr, we let B(x, c) = { ξ lr d : ξ x < c }. Suppose
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 informationUsing Lyapunov Theory I
Lecture 33 Stability Analysis of Nonlinear Systems Using Lyapunov heory I Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore Outline Motivation Definitions
More 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 informationHybrid Control and Switched Systems. Lecture #7 Stability and convergence of ODEs
Hybrid Control and Switched Systems Lecture #7 Stability and convergence of ODEs João P. Hespanha University of California at Santa Barbara Summary Lyapunov stability of ODEs epsilon-delta and beta-function
More informationDO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED.
EE 5 Homewors Fall 04 Updated: Sunday, October 6, 04 Some homewor assignments refer to the textboos: Slotine and Li, etc. For full credit, show all wor. Some problems require hand calculations. In those
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 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 informationConvergence Rate of Nonlinear Switched Systems
Convergence Rate of Nonlinear Switched Systems Philippe JOUAN and Saïd NACIRI arxiv:1511.01737v1 [math.oc] 5 Nov 2015 January 23, 2018 Abstract This paper is concerned with the convergence rate of the
More informationCONTROL SYSTEMS, ROBOTICS AND AUTOMATION - Vol. XII - Lyapunov Stability - Hassan K. Khalil
LYAPUNO STABILITY Hassan K. Khalil Department of Electrical and Computer Enigneering, Michigan State University, USA. Keywords: Asymptotic stability, Autonomous systems, Exponential stability, Global asymptotic
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 informationDelay-independent stability via a reset loop
Delay-independent stability via a reset loop S. Tarbouriech & L. Zaccarian (LAAS-CNRS) Joint work with F. Perez Rubio & A. Banos (Universidad de Murcia) L2S Paris, 20-22 November 2012 L2S Paris, 20-22
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 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 information7.1 Linear Systems Stability Consider the Continuous-Time (CT) Linear Time-Invariant (LTI) system
7 Stability 7.1 Linear Systems Stability Consider the Continuous-Time (CT) Linear Time-Invariant (LTI) system ẋ(t) = A x(t), x(0) = x 0, A R n n, x 0 R n. (14) The origin x = 0 is a globally asymptotically
More information1. Find the solution of the following uncontrolled linear system. 2 α 1 1
Appendix B Revision Problems 1. Find the solution of the following uncontrolled linear system 0 1 1 ẋ = x, x(0) =. 2 3 1 Class test, August 1998 2. Given the linear system described by 2 α 1 1 ẋ = x +
More informationNonlinear Autonomous Systems of Differential
Chapter 4 Nonlinear Autonomous Systems of Differential Equations 4.0 The Phase Plane: Linear Systems 4.0.1 Introduction Consider a system of the form x = A(x), (4.0.1) where A is independent of t. Such
More informationStability and Robustness Analysis of Nonlinear Systems via Contraction Metrics and SOS Programming
arxiv:math/0603313v1 [math.oc 13 Mar 2006 Stability and Robustness Analysis of Nonlinear Systems via Contraction Metrics and SOS Programming Erin M. Aylward 1 Pablo A. Parrilo 1 Jean-Jacques E. Slotine
More informationStabilization and Passivity-Based Control
DISC Systems and Control Theory of Nonlinear Systems, 2010 1 Stabilization and Passivity-Based Control Lecture 8 Nonlinear Dynamical Control Systems, Chapter 10, plus handout from R. Sepulchre, Constructive
More informationCDS 101/110a: Lecture 2.1 Dynamic Behavior
CDS 11/11a: Lecture.1 Dynamic Behavior Richard M. Murray 6 October 8 Goals: Learn to use phase portraits to visualize behavior of dynamical systems Understand different types of stability for an equilibrium
More informationENGI Duffing s Equation Page 4.65
ENGI 940 4. - Duffing s Equation Page 4.65 4. Duffing s Equation Among the simplest models of damped non-linear forced oscillations of a mechanical or electrical system with a cubic stiffness term is Duffing
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 informationHybrid Control and Switched Systems. Lecture #11 Stability of switched system: Arbitrary switching
Hybrid Control and Switched Systems Lecture #11 Stability of switched system: Arbitrary switching João P. Hespanha University of California at Santa Barbara Stability under arbitrary switching Instability
More informationNow I switch to nonlinear systems. In this chapter the main object of study will be
Chapter 4 Stability 4.1 Autonomous systems Now I switch to nonlinear systems. In this chapter the main object of study will be ẋ = f(x), x(t) X R k, f : X R k, (4.1) where f is supposed to be locally Lipschitz
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 informationChapter 13 Internal (Lyapunov) Stability 13.1 Introduction We have already seen some examples of both stable and unstable systems. The objective of th
Lectures on Dynamic Systems and Control Mohammed Dahleh Munther A. Dahleh George Verghese Department of Electrical Engineering and Computer Science Massachuasetts Institute of Technology 1 1 c Chapter
More informationCDS 101/110a: Lecture 2.1 Dynamic Behavior
CDS 11/11a: Lecture 2.1 Dynamic Behavior Richard M. Murray 6 October 28 Goals: Learn to use phase portraits to visualize behavior of dynamical systems Understand different types of stability for an equilibrium
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 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 informationModeling & Control of Hybrid Systems Chapter 4 Stability
Modeling & Control of Hybrid Systems Chapter 4 Stability Overview 1. Switched systems 2. Lyapunov theory for smooth and linear systems 3. Stability for any switching signal 4. Stability for given switching
More informationME 680- Spring Representation and Stability Concepts
ME 680- Spring 014 Representation and Stability Concepts 1 3. Representation and stability concepts 3.1 Continuous time systems: Consider systems of the form x F(x), x n (1) where F : U Vis a mapping U,V
More informationTTK4150 Nonlinear Control Systems Solution 6 Part 2
TTK4150 Nonlinear Control Systems Solution 6 Part 2 Department of Engineering Cybernetics Norwegian University of Science and Technology Fall 2003 Solution 1 Thesystemisgivenby φ = R (φ) ω and J 1 ω 1
More informationWe consider the problem of finding a polynomial that interpolates a given set of values:
Chapter 5 Interpolation 5. Polynomial Interpolation We consider the problem of finding a polynomial that interpolates a given set of values: x x 0 x... x n y y 0 y... y n where the x i are all distinct.
More informationPredictability: Does the Flap of a Butterfly s Wings in Brazil set off a Tornado
Chapter 4 Dynamic Behavior Predictability: Does the Flap of a Butterfly s Wings in Brazil set off a Tornado in Texas? Talk given by Edward Lorenz, December 972 meeting of the American Association for the
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 informationGLOBAL ANALYSIS OF PIECEWISE LINEAR SYSTEMS USING IMPACT MAPS AND QUADRATIC SURFACE LYAPUNOV FUNCTIONS
GLOBAL ANALYSIS OF PIECEWISE LINEAR SYSTEMS USING IMPACT MAPS AND QUADRATIC SURFACE LYAPUNOV FUNCTIONS Jorge M. Gonçalves, Alexandre Megretski y, Munther A. Dahleh y California Institute of Technology
More informationAutomatic Control 2. Nonlinear systems. Prof. Alberto Bemporad. University of Trento. Academic year
Automatic Control 2 Nonlinear systems Prof. Alberto Bemporad University of Trento Academic year 2010-2011 Prof. Alberto Bemporad (University of Trento) Automatic Control 2 Academic year 2010-2011 1 / 18
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 informationOutput Feedback and State Feedback. EL2620 Nonlinear Control. Nonlinear Observers. Nonlinear Controllers. ẋ = f(x,u), y = h(x)
Output Feedback and State Feedback EL2620 Nonlinear Control Lecture 10 Exact feedback linearization Input-output linearization Lyapunov-based control design methods ẋ = f(x,u) y = h(x) Output feedback:
More informationGlobal stabilization of feedforward systems with exponentially unstable Jacobian linearization
Global stabilization of feedforward systems with exponentially unstable Jacobian linearization F Grognard, R Sepulchre, G Bastin Center for Systems Engineering and Applied Mechanics Université catholique
More informationCDS Solutions to the Midterm Exam
CDS 22 - Solutions to the Midterm Exam Instructor: Danielle C. Tarraf November 6, 27 Problem (a) Recall that the H norm of a transfer function is time-delay invariant. Hence: ( ) Ĝ(s) = s + a = sup /2
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 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 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 informationLecture 9 Nonlinear Control Design. Course Outline. Exact linearization: example [one-link robot] Exact Feedback Linearization
Lecture 9 Nonlinear Control Design Course Outline Eact-linearization Lyapunov-based design Lab Adaptive control Sliding modes control Literature: [Khalil, ch.s 13, 14.1,14.] and [Glad-Ljung,ch.17] Lecture
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 informationESC794: Special Topics: Model Predictive Control
ESC794: Special Topics: Model Predictive Control Discrete-Time Systems Hanz Richter, Professor Mechanical Engineering Department Cleveland State University Discrete-Time vs. Sampled-Data Systems A continuous-time
More informationControl of Mobile Robots
Control of Mobile Robots Regulation and trajectory tracking Prof. Luca Bascetta (luca.bascetta@polimi.it) Politecnico di Milano Dipartimento di Elettronica, Informazione e Bioingegneria Organization and
More informationMCE/EEC 647/747: Robot Dynamics and Control. Lecture 12: Multivariable Control of Robotic Manipulators Part II
MCE/EEC 647/747: Robot Dynamics and Control Lecture 12: Multivariable Control of Robotic Manipulators Part II Reading: SHV Ch.8 Mechanical Engineering Hanz Richter, PhD MCE647 p.1/14 Robust vs. Adaptive
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 information