Lyapunov Stability Analysis: Open Loop
|
|
- Jonas Henry Fowler
- 5 years ago
- Views:
Transcription
1 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 by examining the system poles. Continuous-time systems are: marginally stable (MS) if : Poles in the left half of the s-plane, with any jω-axis poles nonrepeated. asymptotically stable (AS) if : Poles strictly in the open left half of the s-plane. Discrete-time systems are: marginally stable (MS) if : Poles in or on the unit circle of the z-plane, with any poles on the unit circle nonrepeated. asymptotically stable (AS) if : Poles strictly inside the unit circle of the z-plane. he Routh est gives an elegant method of determining stability of polynomials without actually finding the roots. Now we will show another method of determining stability of dynamical systems that does not involve the poles. It is based on energy concepts, and can be extended to time-varying and nonlinear systems. Alexsandr Lyapunov was a Russian scientist who in 189 published his work defining the Lyapunov function and using it in the stability analysis of dynamical systems. Lyapunov stability analysis is a general method that can be used for nonlinear systems. It is noteworthy that the first satellite Sputnik was launched by Russia in Nonlinear Systems Analysis was in large part responsible for the control systems on this spacecraft. In the western world at that time, control systems analysis was primarily being carried out in the frequency domain and mainly for linear systems. his was a serious limitation, which became more clear as we advanced into the space age after Aerospace systems are fundamentally nonlinear systems. 1
2 Aleksandr Mikhailovich Lyapunov (from Wikipedia) Born June 6, 1857 Yaroslavl, Imperial Russia Died November 3, 1918 (aged 61) Residence Nationality Fields Institutions Alma mater Known for Russia Russian Applied mathematics Saint Petersburg State University Russian Academy of Sciences Kharkov University Saint Petersburg State University Lyapunov function Lyapunov Analysis for Continuous-ime Systems Stability Analysis Using Energy Consider the continuous-time linear time invariant autonomous (i.e. no control input) system x = Ax, x(0)
3 n with state x( t) R. It is desired to find a way of studying stability that does not involve the system poles or determining the characteristic polynomial Λ () s = si A. Example 1. Stability Analysis Using Energy R u(t) i(t) C x(t) Electric Circuit A circuit is shown. One has u x i = cx = R so that a state variable description is given by 1 1 x = x+ u, x(0) RC RC with x(0) the initial voltage across the capacitor. We are focusing here on open-loop stability, so set u=0 to get the open-loop undriven (or autonomous) dynamics 1 x = x RC We know this is asymptotically stable since the pole is at s = 1/ RC. he natural mode is t/ RC e and the time constant is τ = RC. he state is given by t/ RC x() t = e x(0). o determine stability by another method not involving the poles, consider the energy stored in the capacitor 1 V( x) = Cx 1. Asymptotic Stability. Clearly, since V(x) is quadratic in x(t), it is positive definite, i.e. V( x) 0 V ( x) = 0if and only if x = 0 he energy is a function of time, Vt () = Vxt ( ()). Its derivative is given by 3
4 dv V V ( x) = = x dt x 1 1 V = Cxx = Cx x = x RC R where in the last line we substituted the system dynamics 1 x = x. RC Now note that the derivative of the energy, evaluated thus along the system dynamics, is negative definite, i.e. V ( x) 0 V ( x) = 0if and only if x = 0 hus we are faced with a positive energy (V>0) which is always decreasing ( V < 0 ) as long as the state x(t) is not equal to zero. Consequently, the energy V(t) decreases until it is equal to zero, i.e. until x(t)=0. herefore, the state must go to zero for all initial conditions x(0). his means that the system is asymptotically stable.. Exponential Stability. o find the rate of decrease of the energy, we perform some extra analysis to relate V(x) and V ( x). Note that x = V( x) C so that 1 V = x = V R RC his is a dynamical equation with solution given by / t RC Vt () = e V(0) 1 with initial energy given by Vt ( = 0) = Cx0. herefore the energy decreases exponentially with time constant τ V = RC /. Note that x() t = V() t C x(0) = V(0) C so that 4
5 t/ RC t/ RC xt () = V(0) e = x(0) e. C his exponentially decaying state response was found without determining the poles or solving the state equation. Note that the energy time constant is one half the state time constant since energy is a quadratic function of x(t). Definition of Continuous-ime Lyapunov Function We now formalize these concepts for general continuous-time (C) systems. n Let there be given x( t) R and a C autonomous (i.e. no input) nonlinear dynamic system x = f( x), x(0) Let f(x) be continuously differentiable so that there exist unique solutions, and suppose f(x)=0 so that x=0 is an equilibrium point, i.e. when x=0 one has x = f(0) = 0. his is always the case for linear systems x = Ax. Stability Definitions for Nonlinear Systems 1. he origin of the system is said to be stable in the sense of Lyapunov (SISL), if, for every ε > 0, there exists a δ = δ(ε) > 0 such that, if x(0) < δ then xt () < ε for t 0.. he origin of the system is said to be asymptotically stable (AS) if it is Lyapunov stable and if there exists δ > 0 such that if x(0) < δ then lim xt ( ) = he origin of the system is said to be exponentially stable (ES) if it is asymptotically stable and if there exist α,β,δ > 0 such that if x(0) < δ then x() t α x(0) e βt for t 0. For linear time-invariant systems, SISL is equivalent to marginal stability, AS is the usual asymptotic stability we have discussed before, and AS and ES are equivalent. SISL means that we can keep the state as small as desired for all times t 0 (i.e. xt () any choice of ε) by selecting the initial state small enough (i.e. x(0) < δ ( ε ) ). Definiteness of Functions and Matrices t < ε for Since Lyapunov techniques are extensions of energy concepts, the following definitions are germane. 5
6 n Definition: Definiteness of Real Functions. We say a real function x ( ) : R Ris: positive semidefinite, x ( ) 0 if x ( ) 0, x positive definite, (x)>0, if x ( ) 0, x ( x) = 0if and only if x = 0 negative semidefinite, x ( ) 0if x ( ) 0, x negative definite, (x)<0, if x ( ) 0, x ( x) = 0if and only if x = 0 Let (x) be the quadratic form x ( ) = xx with an n n real matrix. uadratic forms are very convenient for Lyapunov analysis. We define the following properties for square matrices. Definition: Definiteness of Square Matrices. We say a real square n n matrix is: positive semidefinite, 0 if x x 0, x positive definite, >0, if xx 0, x x x = 0if and only if x = 0 negative semidefinite, 0 if x x 0, x negative definite, <0, if xx 0, x x x = 0if and only if x = 0 If is a symmetric matrix, then the definiteness properties of (x) can be determined in terms of the eigenvalues of : 6
7 Fact. Definiteness of Real uadratic Functions. A real square symmetric matrix = is: positive semidefinite, if all the eigenvalues of are greater than or equal to zero positive definite, if all the eigenvalues of are greater than zero negative semidefinite, if all the eigenvalues of are less than or equal to zero negative definite, if all the eigenvalues of are less than zero he next result gives an alternative test for determining definiteness of matrices that is n n often easier than finding the eigenvalues. Let = [ q ] R. he leading minors of are q11 q1 q13 q11 q1 m1 = q11, m =, m3 = q1 q q3,, mn = q1 q q q q In terms of the leading minors, we have ij > 0 if m i > 0, all i < 0 if mi mi < 0, > 0, all odd i all even i. Lyapunov heorems for Continuous-ime Systems Now we formalize the results of Example 1. Definition: Continuous-time (C) Lyapunov Function V(x) is said to be a C :Lyapunov Function if: 1. V(x) is a continuous real-valued function. V( x ) > 0, i.e. V(x) is positive definite 3. V = dv 0, i.e. V ( x) is negative semi-definite dt Condition 3 must be verified along the system trajectories. hat is, compute dv/dt and substitute in the system dynamics x = f( x). 7
8 he following results were proven by A.M. Lyapunov. Lyapunov heorem. SISL. Suppose that for a given system there exists a Lyapunov function. hen the system is SISL. Lyapunov heorem. AS. Suppose that for a given system there exists a Lyapunov function which also satisfies the stronger third condition V = dv < 0, i.e. V ( x) is negative definite. hen dt the system is AS. As shown in Example 1, exponential stability can be verified by extra analysis to relate V(x) and V ( x) and show that the Lyapunov function decreases exponentially. hese theorems are powerful and elegant. However, they do not say how to find a Lyapunov function. his may be difficult for nonlinear systems. Often, energies can be used. he next definition is useful in finding Lyapunov functions for specific systems. Definition. Lyapunov Function Candidate (LFC). V(x) is said to be a C :Lyapunov Function Candidate (LFC) if: 1. V(x) is a continuous real-valued function. V( x ) > 0, i.e. V(x) is positive definite Having selected a reasonable LFC, one computes its time derivative (along the system trajectories) to determine whether condition 3. holds, i.e. whether the LFC is in fact a Lyapunov function that verifies the stability of the system. Lyapunov Functions for C Linear systems: C Lyapunov Equation. Consider the continuous-time (C) linear time invariant (LI) autonomous system x = Ax, x(0) n with state x( t) R. For LI systems it is easy to find a Lyapunov function. In fact, quadratic energy always works. herefore define the energy quadratic function 1 V( x) = x Px with constant kernel matrix P symmetric and positive definite, i.e. P= P > 0. Clearly 1. V(x) is a continuous real-valued function 8
9 . V( x ) > 0, i.e. V(x) is positive definite herefore, V(x) is a Lyapunov function candidate LFC. If we can prove that V = dv 0 along the system trajectories, then V(x) is a Lyapunov dt function and the system is SISL. Compute the Lyapunov derivative d 1 1 V = ( x Px) = ( x Px + x Px + x Px ) dt Now P is constant, and substituting the system dynamics x = Ax yields 1 1 V = ( ( Ax) Px + x PAx) = x ( A P + PA) x along the system trajectories. Now suppose that there exists a positive semidefinite matrix such that AP+ PA= hen 1 1 V = x ( A P+ PA) x= x x 0 hen, V(x) is a Lyapunov function and the system is SISL. We call AP+ PA= the continuous-time Lyapunov Equation. Note that this equation is the same if it is transposed. herefore, if is symmetric, so is the solution P. herefore, we see that if there exists a symmetric positive definite solution P to the C Lyapunov equation with a symmetric positive semidefinite matrix, then V(x) is a Lyapunov function and the system is SISL. In fact, this condition turns out to be necessary and sufficient for LI C systems, as formulated in the next result. heorem. SISL for LI C Systems. Let be a symmetric positive semidefinite matrix. hen the system x = Ax is SISL (e.g. marginally stable) if and only if the (symmetric) matrix P which solves the C Lyapunov equation AP+ PA= is positive definite. If = is in fact positive definite, we have the next stronger result. heorem. AS for LI C Systems. Let be a symmetric positive definite matrix. hen the system x = Ax is AS if and only if the (symmetric) matrix P which solves the C Lyapunov equation 9
10 AP+ PA= is positive definite. In fact, it turns out that can always be selected as the n n identity matrix to test a system for AS using this method. Example. Stability by Lyapunov Equation Solution for C LI Systems. Consider the LI C autonomous system 0 1 x = Ax = x 3 Let us test for stability using the Lyapunov equation. First, let us try for AS. Select =I. he Lyapunov equation becomes PA + A P = p1 p p1 p 1 0 p p + = p p where symmetry of P has been used. his yields 4p p1 3p p3 1 0 = p1 3p p3 p 6p Now the top left (e.g. (1,1) element) yields 4p = 1 p = 1/4 and then the bottom left (.) equation yields p 6p3 = 1 p 3 = 1/4 Finally, p1 3p p3 = 0 so that p 1 = 5/4 herefore, 5 1 p1 p 4 4 P = 1 1 p p = he leading minors of P are 10
11 p = > P = > herefore P>0. hus the system is AS. Notes: 1. In this case one could solve for the elements of P sequentially because of the easy structure of matrix A. Generally this is not possible. hen, write 4p p1 3p p3 1 0 = p1 3p p3 p 6p as p p 0 = 0 6 p 3 1 and invert the coefficient matrix to find [ p p p ] 1 3. If in some example one selects =I and the solution P of the Lyapunov equation is not positive definite, the system is not AS. However, it may still be SISL. o verify this, select a that is positive semidefinite and see if the resulting Lyapunov 1 0 solution P is positive definite. If so, the system is MS. For instance, one could try = or = 0 1. Lyapunov Analysis for Discrete-ime Systems Consider now the discrete-time autonomous nonlinear system x f( x ), x = k+ 1 k 0 Suppose f(x)=0 so that x=0 is an equilibrium point. his is always the case for linear D systems. Given a real function V( x k ) define the first difference as Δ V( xk) = V( xk+ 1) V( xk) Definition: Discrete-ime (D) Lyapunov Function V(x) is said to be a D :Lyapunov Function if: 11
12 1. V(x) is a continuous real-valued function. V( x ) > 0, i.e. V(x) is positive definite 3. Δ V( xk) = V( xk+ 1) V( xk) 0, i.e. the first difference Δ V( x k ) is negative semidefinite he following results were proven by A.M. Lyapunov. SISL and AS are defined for D systems in terms of time trajectories the same as for C systems. D Lyapunov heorem. SISL. Suppose that for a given system there exists a Lyapunov function. hen the system is SISL. D Lyapunov heorem. AS. Suppose that for a given system there exists a Lyapunov function which also satisfies the stronger third condition Δ V( xk) = V( xk+ 1) V( xk) < 0, i.e. Δ V( x k ) is negative definite. hen the system is AS. In the D case stability is referred to the unit circle in the z-plane. he next definition is useful in finding Lyapunov functions for specific systems. Definition. Lyapunov Function Candidate (LFC). V(x) is said to be a C :Lyapunov Function Candidate (LFC) if: 1. V(x) is a continuous real-valued function. V( x ) > 0, i.e. V(x) is positive definite Having selected a reasonable LFC, one computes its first difference (along the system trajectories) to determine whether condition 3. holds, i.e. whether the LFC is in fact a Lyapunov function that verifies the stability of the system. Lyapunov Functions for D Linear systems: D Lyapunov Equation. Consider the discrete-time (D) linear time invariant (LI) autonomous system x = k+ 1 Axk, x0 n with state R. xk For LI systems it is easy to find a Lyapunov function. In fact, quadratic energy always works. herefore define the energy quadratic function 1
13 1 V( xk) = xkpxk with constant kernel matrix P symmetric and positive definite, i.e. P= P > 0. Clearly 1. V(x) is a continuous real-valued function. V( x ) > 0, i.e. V(x) is positive definite herefore, V(x) is a Lyapunov function candidate LFC. If we can prove that Δ V( xk) = V( xk+ 1) V( xk) 0 along the system trajectories, then V(x) is a D Lyapunov function and the system is SISL. o study this, compute the Lyapunov function first difference Δ V( xk) = V( xk+ 1) V( xk) = xk+ 1Pxk+ 1 xkpxk Now substituting the system dynamics x = k+ 1 Ax yields k Δ V( xk) = xk ( A PA P) xk along the system trajectories. Now suppose that there exists a positive semidefinite matrix such that APA P= hen 1 1 Δ V( x ) = x ( A PA P) x = x x 0 k k k k k hen, V(x k ) is a Lyapunov function and the system is SISL. We call APA P= the discrete-time Lyapunov Equation. Note that this equation is the same if it is transposed. herefore, if is symmetric, so is the solution P. herefore, we see that if there exists a symmetric positive definite solution P to the D Lyapunov equation with a symmetric positive semidefinite matrix, then V(x k ) is a D Lyapunov function and the system is SISL. In fact, this condition turns out to be necessary and sufficient for LI C systems, as formulated in the next result. heorem. SISL for LI D Systems. Let be a symmetric positive semidefinite matrix. hen the system xk+ 1 = Axk is SISL (e.g. marginally stable) if and only if the (symmetric) matrix P which solves the D Lyapunov equation APA P= is positive definite. If = is in fact positive definite, we have the next stronger result. 13
14 heorem. AS for LI C Systems. Let be a symmetric positive definite matrix. hen the system xk+ 1 = Axk is AS if and only if the (symmetric) matrix P which solves the D Lyapunov equation APA P= is positive definite. In fact, it turns out that can always be selected as the n n identity matrix to test a system for AS using this method. Note that the solution properties of the C Lyapunov equation AP+ PA= refer to the locations of the poles of system matrix A with respect to the left half of the complex plane. On the other hand, the solution properties of the D Lyapunov equation APA P= refer to the locations of the poles of system matrix A with respect to the unit circle of the complex plane. 14
7.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 informationEE Control Systems LECTURE 9
Updated: Sunday, February, 999 EE - Control Systems LECTURE 9 Copyright FL Lewis 998 All rights reserved STABILITY OF LINEAR SYSTEMS We discuss the stability of input/output systems and of state-space
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 informationECE504: Lecture 8. D. Richard Brown III. Worcester Polytechnic Institute. 28-Oct-2008
ECE504: Lecture 8 D. Richard Brown III Worcester Polytechnic Institute 28-Oct-2008 Worcester Polytechnic Institute D. Richard Brown III 28-Oct-2008 1 / 30 Lecture 8 Major Topics ECE504: Lecture 8 We are
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 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 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 informationSolution of Linear State-space Systems
Solution of Linear State-space Systems Homogeneous (u=0) LTV systems first Theorem (Peano-Baker series) The unique solution to x(t) = (t, )x 0 where The matrix function is given by is called the state
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 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 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 informationLyapunov Stability Analysis for Feedback Control Design
Copyright F.L. Lewis 008 All rights reserved Updated: uesday, November, 008 Lyapuov Stability Aalysis for Feedbac Cotrol Desig Lyapuov heorems Lyapuov Aalysis allows oe to aalyze the stability of cotiuous-time
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 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 informationIntroduction to Nonlinear Control Lecture # 4 Passivity
p. 1/6 Introduction to Nonlinear Control Lecture # 4 Passivity È p. 2/6 Memoryless Functions ¹ y È Ý Ù È È È È u (b) µ power inflow = uy Resistor is passive if uy 0 p. 3/6 y y y u u u (a) (b) (c) Passive
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 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 informationModule 06 Stability of Dynamical Systems
Module 06 Stability of Dynamical Systems Ahmad F. Taha EE 5143: Linear Systems and Control Email: ahmad.taha@utsa.edu Webpage: http://engineering.utsa.edu/ataha October 10, 2017 Ahmad F. Taha Module 06
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 informationME Fall 2001, Fall 2002, Spring I/O Stability. Preliminaries: Vector and function norms
I/O Stability Preliminaries: Vector and function norms 1. Sup norms are used for vectors for simplicity: x = max i x i. Other norms are also okay 2. Induced matrix norms: let A R n n, (i stands for induced)
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 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 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 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 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 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 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 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 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 informationEE363 homework 2 solutions
EE363 Prof. S. Boyd EE363 homework 2 solutions. Derivative of matrix inverse. Suppose that X : R R n n, and that X(t is invertible. Show that ( d d dt X(t = X(t dt X(t X(t. Hint: differentiate X(tX(t =
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 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 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 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 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 informationEECE Adaptive Control
EECE 574 - Adaptive Control Model-Reference Adaptive Control - Part I Guy Dumont Department of Electrical and Computer Engineering University of British Columbia January 2010 Guy Dumont (UBC EECE) EECE
More informationNonlinear Control. Nonlinear Control Lecture # 6 Passivity and Input-Output Stability
Nonlinear Control Lecture # 6 Passivity and Input-Output Stability Passivity: Memoryless Functions y y y u u u (a) (b) (c) Passive Passive Not passive y = h(t,u), h [0, ] Vector case: y = h(t,u), h T =
More informationMULTIVARIABLE ZEROS OF STATE-SPACE SYSTEMS
Copyright F.L. Lewis All rights reserved Updated: Monday, September 9, 8 MULIVARIABLE ZEROS OF SAE-SPACE SYSEMS If a system has more than one input or output, it is called multi-input/multi-output (MIMO)
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 of Parameter Adaptation Algorithms. Big picture
ME5895, UConn, Fall 215 Prof. Xu Chen Big picture For ˆθ (k + 1) = ˆθ (k) + [correction term] we haven t talked about whether ˆθ(k) will converge to the true value θ if k. We haven t even talked about
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 informationControl engineering sample exam paper - Model answers
Question Control engineering sample exam paper - Model answers a) By a direct computation we obtain x() =, x(2) =, x(3) =, x(4) = = x(). This trajectory is sketched in Figure (left). Note that A 2 = I
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 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 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 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 informationLecture 2 and 3: Controllability of DT-LTI systems
1 Lecture 2 and 3: Controllability of DT-LTI systems Spring 2013 - EE 194, Advanced Control (Prof Khan) January 23 (Wed) and 28 (Mon), 2013 I LTI SYSTEMS Recall that continuous-time LTI systems can be
More informationEML5311 Lyapunov Stability & Robust Control Design
EML5311 Lyapunov Stability & Robust Control Design 1 Lyapunov Stability criterion In Robust control design of nonlinear uncertain systems, stability theory plays an important role in engineering systems.
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 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 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 informationTime-Invariant Linear Quadratic Regulators Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 2015
Time-Invariant Linear Quadratic Regulators Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 15 Asymptotic approach from time-varying to constant gains Elimination of cross weighting
More informationCALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems
CDS 101 1. For each of the following linear systems, determine whether the origin is asymptotically stable and, if so, plot the step response and frequency response for the system. If there are multiple
More informationSome of the different forms of a signal, obtained by transformations, are shown in the figure. jwt e z. jwt z e
Transform methods Some of the different forms of a signal, obtained by transformations, are shown in the figure. X(s) X(t) L - L F - F jw s s jw X(jw) X*(t) F - F X*(jw) jwt e z jwt z e X(nT) Z - Z X(z)
More informationTime-Invariant Linear Quadratic Regulators!
Time-Invariant Linear Quadratic Regulators Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 17 Asymptotic approach from time-varying to constant gains Elimination of cross weighting
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 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 informationOptimal control and estimation
Automatic Control 2 Optimal control and estimation Prof. Alberto Bemporad University of Trento Academic year 2010-2011 Prof. Alberto Bemporad (University of Trento) Automatic Control 2 Academic year 2010-2011
More informationCDS Solutions to Final Exam
CDS 22 - Solutions to Final Exam Instructor: Danielle C Tarraf Fall 27 Problem (a) We will compute the H 2 norm of G using state-space methods (see Section 26 in DFT) We begin by finding a minimal state-space
More informationRobot Manipulator Control. Hesheng Wang Dept. of Automation
Robot Manipulator Control Hesheng Wang Dept. of Automation Introduction Industrial robots work based on the teaching/playback scheme Operators teach the task procedure to a robot he robot plays back eecute
More informationStability and hybrid synchronization of a time-delay financial hyperchaotic system
ISSN 76-7659 England UK Journal of Information and Computing Science Vol. No. 5 pp. 89-98 Stability and hybrid synchronization of a time-delay financial hyperchaotic system Lingling Zhang Guoliang Cai
More informationEE 380. Linear Control Systems. Lecture 10
EE 380 Linear Control Systems Lecture 10 Professor Jeffrey Schiano Department of Electrical Engineering Lecture 10. 1 Lecture 10 Topics Stability Definitions Methods for Determining Stability Lecture 10.
More informationLinear Quadratic Regulator (LQR) Design I
Lecture 7 Linear Quadratic Regulator LQR) Design I Dr. Radhakant Padhi Asst. Proessor Dept. o Aerospace Engineering Indian Institute o Science - Bangalore LQR Design: Problem Objective o drive the state
More informationPerspective. ECE 3640 Lecture 11 State-Space Analysis. To learn about state-space analysis for continuous and discrete-time. Objective: systems
ECE 3640 Lecture State-Space Analysis Objective: systems To learn about state-space analysis for continuous and discrete-time Perspective Transfer functions provide only an input/output perspective of
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 informationLinear dynamical systems with inputs & outputs
EE263 Autumn 215 S. Boyd and S. Lall Linear dynamical systems with inputs & outputs inputs & outputs: interpretations transfer function impulse and step responses examples 1 Inputs & outputs recall continuous-time
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 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 informationVideo 8.1 Vijay Kumar. Property of University of Pennsylvania, Vijay Kumar
Video 8.1 Vijay Kumar 1 Definitions State State equations Equilibrium 2 Stability Stable Unstable Neutrally (Critically) Stable 3 Stability Translate the origin to x e x(t) =0 is stable (Lyapunov stable)
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 informationRobust Control 2 Controllability, Observability & Transfer Functions
Robust Control 2 Controllability, Observability & Transfer Functions Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University /26/24 Outline Reachable Controllability Distinguishable
More informationReview: control, feedback, etc. Today s topic: state-space models of systems; linearization
Plan of the Lecture Review: control, feedback, etc Today s topic: state-space models of systems; linearization Goal: a general framework that encompasses all examples of interest Once we have mastered
More information3 rd Tutorial on EG4321/EG7040 Nonlinear Control
3 rd Tutorial on EG4321/EG7040 Nonlinear Control Lyapunov Stability Dr Angeliki Lekka 1 1 Control Systems Research Group Department of Engineering, University of Leicester arch 9, 2017 Dr Angeliki Lekka
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 informationLinear Systems. Manfred Morari Melanie Zeilinger. Institut für Automatik, ETH Zürich Institute for Dynamic Systems and Control, ETH Zürich
Linear Systems Manfred Morari Melanie Zeilinger Institut für Automatik, ETH Zürich Institute for Dynamic Systems and Control, ETH Zürich Spring Semester 2016 Linear Systems M. Morari, M. Zeilinger - Spring
More information6.241 Dynamic Systems and Control
6.241 Dynamic Systems and Control Lecture 24: H2 Synthesis Emilio Frazzoli Aeronautics and Astronautics Massachusetts Institute of Technology May 4, 2011 E. Frazzoli (MIT) Lecture 24: H 2 Synthesis May
More informationw T 1 w T 2. w T n 0 if i j 1 if i = j
Lyapunov Operator Let A F n n be given, and define a linear operator L A : C n n C n n as L A (X) := A X + XA Suppose A is diagonalizable (what follows can be generalized even if this is not possible -
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 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 informationLIE TRANSFORM CONSTRUCTION OF EXPONENTIAL NORMAL FORM OBSERVERS. Suba Thomas Harry G. Kwatny Bor-Chin Chang
LIE TRANSFORM CONSTRUCTION OF EXPONENTIAL NORMAL FORM OBSERVERS Suba Thomas Harry G Kwatny Bor-Chin Chang Department of Mechanical Engineering and Mechanics, Drexel University, Philadelphia, PA 1914, USA
More informationRobust Control 5 Nominal Controller Design Continued
Robust Control 5 Nominal Controller Design Continued Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University 4/14/2003 Outline he LQR Problem A Generalization to LQR Min-Max
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 informationControl of industrial robots. Centralized control
Control of industrial robots Centralized control Prof. Paolo Rocco (paolo.rocco@polimi.it) Politecnico di Milano ipartimento di Elettronica, Informazione e Bioingegneria Introduction Centralized control
More informationSolution via Laplace transform and matrix exponential
EE263 Autumn 2015 S. Boyd and S. Lall Solution via Laplace transform and matrix exponential Laplace transform solving ẋ = Ax via Laplace transform state transition matrix matrix exponential qualitative
More informationApproximation Approach for Timing Jitter Characterization in Circuit Simulators
Approximation Approach for iming Jitter Characterization in Circuit Simulators MM.Gourary,S.G.Rusakov,S.L.Ulyanov,M.M.Zharov IPPM, Russian Academy of Sciences, Moscow, Russia K.K. Gullapalli, B. J. Mulvaney
More informationChap. 3. Controlled Systems, Controllability
Chap. 3. Controlled Systems, Controllability 1. Controllability of Linear Systems 1.1. Kalman s Criterion Consider the linear system ẋ = Ax + Bu where x R n : state vector and u R m : input vector. A :
More informationProblem Description The problem we consider is stabilization of a single-input multiple-state system with simultaneous magnitude and rate saturations,
SEMI-GLOBAL RESULTS ON STABILIZATION OF LINEAR SYSTEMS WITH INPUT RATE AND MAGNITUDE SATURATIONS Trygve Lauvdal and Thor I. Fossen y Norwegian University of Science and Technology, N-7 Trondheim, NORWAY.
More informationSome solutions of the written exam of January 27th, 2014
TEORIA DEI SISTEMI Systems Theory) Prof. C. Manes, Prof. A. Germani Some solutions of the written exam of January 7th, 0 Problem. Consider a feedback control system with unit feedback gain, with the following
More informationRobotics. Control Theory. Marc Toussaint U Stuttgart
Robotics Control Theory Topics in control theory, optimal control, HJB equation, infinite horizon case, Linear-Quadratic optimal control, Riccati equations (differential, algebraic, discrete-time), controllability,
More informationPhysics 116A Notes Fall 2004
Physics 116A Notes Fall 2004 David E. Pellett Draft v.0.9 Notes Copyright 2004 David E. Pellett unless stated otherwise. References: Text for course: Fundamentals of Electrical Engineering, second edition,
More informationFirst and Second Order Circuits. Claudio Talarico, Gonzaga University Spring 2015
First and Second Order Circuits Claudio Talarico, Gonzaga University Spring 2015 Capacitors and Inductors intuition: bucket of charge q = Cv i = C dv dt Resist change of voltage DC open circuit Store voltage
More informationState Regulator. Advanced Control. design of controllers using pole placement and LQ design rules
Advanced Control State Regulator Scope design of controllers using pole placement and LQ design rules Keywords pole placement, optimal control, LQ regulator, weighting matrixes Prerequisites Contact state
More informationLecture 8. Chapter 5: Input-Output Stability Chapter 6: Passivity Chapter 14: Passivity-Based Control. Eugenio Schuster.
Lecture 8 Chapter 5: Input-Output Stability Chapter 6: Passivity Chapter 14: Passivity-Based Control Eugenio Schuster schuster@lehigh.edu Mechanical Engineering and Mechanics Lehigh University Lecture
More informationOn the convergence of normal forms for analytic control systems
Chapter 1 On the convergence of normal forms for analytic control systems Wei Kang Department of Mathematics, Naval Postgraduate School Monterey, CA 93943 wkang@nps.navy.mil Arthur J. Krener Department
More informationInput-to-state stability and interconnected Systems
10th Elgersburg School Day 1 Input-to-state stability and interconnected Systems Sergey Dashkovskiy Universität Würzburg Elgersburg, March 5, 2018 1/20 Introduction Consider Solution: ẋ := dx dt = ax,
More informationLecture 12: Detailed balance and Eigenfunction methods
Miranda Holmes-Cerfon Applied Stochastic Analysis, Spring 2015 Lecture 12: Detailed balance and Eigenfunction methods Readings Recommended: Pavliotis [2014] 4.5-4.7 (eigenfunction methods and reversibility),
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 informationLinear Matrix Inequalities in Control
Linear Matrix Inequalities in Control Guido Herrmann University of Leicester Linear Matrix Inequalities in Control p. 1/43 Presentation 1. Introduction and some simple examples 2. Fundamental Properties
More informationIntermediate Differential Equations. Stability and Bifurcation II. John A. Burns
Intermediate Differential Equations Stability and Bifurcation II John A. Burns Center for Optimal Design And Control Interdisciplinary Center for Applied Mathematics Virginia Polytechnic Institute and
More information