STABILITY ANALYSIS OF DYNAMIC SYSTEMS
|
|
- Pearl Hicks
- 5 years ago
- Views:
Transcription
1 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 dell Informazione (DEI) Università di Bologna claudio.melchiorri@unibo.it
2 C. Melchiorri (DEI) Automatic Control & System Theory Definition of stability Let us consider a generic dynamic systems (linear/non-linear, time varying, stationary) described by Assume that U, U f, X, Y are normed vector spaces (equipped with a norm). Given the initial time instant t, the initial state x(t ) and the input function u(.), let consider the reference motion We are interested to study two types of perturbations : 1) Variation of motion due to a perturbation of the initial state: Difference between motions ) Variation of motion due to a perturbation of the input: Difference between motions
3 C. Melchiorri (DEI) Automatic Control & System Theory 3 Local stability Initial state perturbation R 3 The motion x(.) corresponding to the input function u(.) is stable with respect to perturbations of the initial state if: if The motion x(.) is asimptotically stable with respect to perturbations of the initial state if: if
4 Local stability Input perturbation The motion x(.) corresponding to the input function u(.) is stable with respect to perturbations of the input function if: if Notes 1. These definitions of stability refer to a particular motion, defined by the three elements t, x(t ), u(.).. The stability definitions, besides to ϕ(t,t,x(t ),u(.)), can also be applied to the response function γ(t,t,x(t ),u(.)). 3. Stability can be referred to an equilibrum state x(t)=ϕ(t,t,x(t ),u(.)), that is a particular motion of the system. C. Melchiorri (DEI) Automatic Control & System Theory 4
5 C. Melchiorri (DEI) Automatic Control & System Theory 5 Local stability Notes 4. If the system is stationary, stability does not depend on t but on x(t ), u(.) only 5. A particular motion may correspond to different input functions, and its stability may depend on the applied input Example: x=, u(.) is an equilibrium state: u > Mg stable u < Mg unstable 6. The given definition of stability refers to a local stability property (for small perturbations).
6 C. Melchiorri (DEI) Automatic Control & System Theory 6 Stability in large Stability in large: initial state perturbations Let consider the asymptotic stability with respect to perturbations of the initial state; if δx 1 (t ) such that x(t )+δ x 1 (t ) X (t, x(t ), u( )), then X is called asymptotic stability domain with respect to t, x(t ), u( ). If X (t, x(t ), u( )) = X (i.e. the whole state space) then the motion x( ) corresponding to the input function u( ) is globally asymptotically stable. If this condiiton holds for all t and all u( ) the system is globally asymptotically stable.
7 C. Melchiorri (DEI) Automatic Control & System Theory 7 Stability in large Example: x= is an equilibrium point whose stability domain is X < x < X.
8 Stability in large Stability in large: input perturbations Let consider again the definition of stability with respect to variations of the input function. The system is bounded input - bounded state stable (b.i.b.s.) with respect to the motion t, x(t ), u( ) if two positive real numbers M u and M x exist, in general function of t, x(t ), u( ), such that for all δu( ) such that δu(t) < M u, t t. Similarly, by considering the response function γ(t,t,x(t ),u( )), it is possible to define the bounded input bounded output stability (b.i.b.o.). C. Melchiorri (DEI) Automatic Control & System Theory 8
9 C. Melchiorri (DEI) Automatic Control & System Theory 9 Stability of the equilibrium Stability of an equilibrium state Let consider an autonomous (homogeneous) system Def. 1: x e is an equilibrium point (state) if =f(x e,t), 8t t Let consider now system (1) and assume, without loss of generality, that the zero state is an equilibrium state (x e = ), i.e.
10 C. Melchiorri (DEI) Automatic Control & System Theory 1 Stability of the equilibrium Stability of an equilibrium state The assumption of considering the zero state as an equilibrium state is not restrictive. As a matter of fact, if the equilibrium is for x e (i.e. f(x e,t)=), it is possible to consider a new state z and a new function g(z, t) such that In which, obviously, we have
11 C. Melchiorri (DEI) Automatic Control & System Theory 11 Stability of the equilibrium Stability of an equilibrium state Let consider an autonomous system ẋ(t) =f(x(t)), (1) Note: A system is autonomous if it is without input and does not depend explicitly on time. Def. : The zero state of system (1) is stable in the sense of Lyapunov at time instant t if for all ε > there exists a parameter η > such that Def. 3: The zero state of system (1) is asymptotically stable in the sense of Lyapunov at time instant t if it is stable and if Note: These definitions apply to the equilibrium state, not to the system. A dynamic system can have several equilibrium points.
12 C. Melchiorri (DEI) Automatic Control & System Theory 1 Stability of the equilibrium Def. 4: The zero state of system (1) is said to be globally asymptotically stable or asymptotically stable at large in the sense of Lyapunov at time instant t if it is stable and if where X is the whole state space. In this case, since no other equilibrium state exists, the system (1) is said to be globally asymptotically stable.
13 C. Melchiorri (DEI) Automatic Control & System Theory 13 Stability of the equilibrium t x 1 3 ε t In summary: 1. Stable. Asymptotically stable 3. Unstable η x 1 If the previous definitions hold for all initial time instants t 1 [t, ], the stability property is said to be uniform in [t, ].
14 C. Melchiorri (DEI) Automatic Control & System Theory 14 Stability of the equilibrium Def. 5: The zero state of system (1) is said to be (locally) exponentially stable if there exist two real constants α, λ > such that x(t) x e α x() x e e λt t > whenever x() x e < δ. It is said to be globally exponentially stable if this property holds for any x R n. Note: Clearly, exponential stability implies asymptotic stability. The converse is, however, not true.
15 The integrator as a dynamic system Example: the integrator u(s) 1 s Integrator y(s) We want to study the time evolution of the state x(t) with initial condition x and null input. 1) Time evolution of the state x(t) due to perturbations of x. Given x x(t) = x Given a perturbation of the initial state δx x(t) = x + δx The perturbations on the time evolution of the state and of the output are bounded (globally stable in the sense of Lyapunov non asymptotically stable) ) Time evolution of the state x(t) due to perturbations of the input. Given x and u(t)= x(t) = x In case of a (constant) perturbation of the input δu x(t) = x + δu t The perturbations on the time evolution of the state and of the output are not bounded, x(t) for t (non b.i.b.s. stable and non b.i.b.o. stable) C. Melchiorri (DEI) Automatic Control & System Theory 15
16 C. Melchiorri (DEI) Automatic Control & System Theory 16 Analysis of stability for dynamic systems The study of the stability properties is of fundamental importance for the analysis of dynamic systems and the design of control laws. In the following, we discuss two methods applicable to both linear and non linear systems expressed in state form ẋ(t) =f(x, u, t), x(t )=x 1) First Lyapunov method (local linearization) ) Second Lyapunov method (direct method) Aleksandr Mikhailovich Lyapunov ( ) was a Russian mathematician, mechanician and physicist. Lyapunov is known for his development of the stability theory of a dynamical system, as well as for his many contributions to mathematical physics and probability theory.
17 C. Melchiorri (DEI) Automatic Control & System Theory 17 First Lyapunov method (local linearization) Let us consider the autonoumous system where x = is an equilibrium point. If it is possible to compute the derivative of the functions f( ), then in a proper neighbourhood of the origin it is possible to write, by using the Taylor series expansion: ẋ(t) =f(x(t)) = f() + f x x= x(t)+f (x(t)) where f (x(t)) are the residuals of the Taylor series expansion, assumed negligible. Since f() =, by defining then the homogeneous linear system ẋ(t) =Ax(t) (3) represents the local linearization or local approximation of the non linear system ().
18 C. Melchiorri (DEI) Automatic Control & System Theory 18 First Lyapunov method (local linearization) Theorem: (Lyapunov linearization method) 1) If the linearized system (3) is strictly stable (all the eigenvalues of A are with negative real part), the equilibrium point is asymptotically stable for system (). ) If the linearized system (3) is unstable (at least an eigenvalue of A has real positive part), the equilibrium point is unstable for system (). 3) If A has eigenvalues with null real part, the linearization does not give information on the stability of the considered equilibrium point.
19 C. Melchiorri (DEI) Automatic Control & System Theory 19 First Lyapunov method (local linearization) Example 1 f(x) x = 4 f 1 (x) x 1 f (x) x 1 f 1 (x) x f (x) x 3 5 = " sin x x 1 cos x 1 x 1 cos x x 1 cos x # Asimptotycally stable system
20 C. Melchiorri (DEI) Automatic Control & System Theory First Lyapunov method (local linearization) Example 1 Evolution of the original and of the linearized system (dash) 1.4 Original system Linearized system
21 C. Melchiorri (DEI) Automatic Control & System Theory 1 First Lyapunov method Example - Pendulum Gravity force Viscous friction Input If u(t) =, the system has equilibrium states: x e,1 = apple x e, = apple m g
22 First Lyapunov method Example - Pendulum f(x) x = 4 f 1 (x) x 1 f (x) x 1 f 1 (x) x f (x) x 3 5 = " 1 g L cos x 1 d L m # Then: STABLE UNSTABLE An eigenvalue is positive C. Melchiorri (DEI) Automatic Control & System Theory
23 First Lyapunov method Example - Pendulum f(x) x = 4 f 1 (x) x 1 f (x) x 1 f 1 (x) x f (x) x 3 5 = " 1 g L cos x 1 d L m # By considering d =, then:??????? UNSTABLE C. Melchiorri (DEI) Automatic Control & System Theory 3
24 C. Melchiorri (DEI) Automatic Control & System Theory 4 Lyapunov method Definitions Positive definite functions: A function V(x): R n -> R is positive definite in ϑ(,ρ), a spherical neighbourhood of the origin of the state space, if 1. V() = ;. V(x) > x Lyapunov functions: A function V(x): R n à R is a Lyapunov function in ϑ(,ρ) for system () if 1. V(x) is positive definite in ϑ(,ρ);. V(x) has continuous first-order partial derivatives with respect to x; 3.
25 C. Melchiorri (DEI) Automatic Control & System Theory 5 Lyapunov method Positive definite functions
26 C. Melchiorri (DEI) Automatic Control & System Theory 6 Lyapunov method Undefinite function Positive definite function X
27 Lyapunov method Lyapunov theorems 1. (simple stability): The zero state is stable in the sense of Lyapunov if there exists a Lyapunov function V in a neighbourhood of the origin ϑ(,ρ).. (asymptotic stability): The zero state is asymptotically stable (AS) in the sense of Lyapunov if there exists a Lyapunov function V in a neighbourhood of the origin ϑ(,ρ) such that dv(x)/dt < for all x ϑ(,ρ), x. 3. (global asymptotic stability): System () is globally asymptotically stable (GAS) in x = if there exists a Lyapunov function V defined in all the state space such that: (radially unlimited Lyapunov function). C. Melchiorri (DEI) Automatic Control & System Theory 7
28 Lyapunov method 4. (asymptotic stability domain): Let V(x) be a Lyapunov function and h a positive real number such that the (open) set is bounded and let dv(x)/dt < for all x D, x. Then, all trajectories starting from a point in D converge asymptotically to zero. V(x) t Time evolution of the Lyapunov function x 1 x() x Time evolution of the state of an autonomous system C. Melchiorri (DEI) Automatic Control & System Theory 8
29 C. Melchiorri (DEI) Automatic Control & System Theory 9 Lyapunov method NOTE: The function describes the variation in time of function V(x) when this is computed along a solution x(t) of the differential equation. NOTE: The previous theorems can be applied to non autonomous systems This can be made by defining a Lyapunov function V(x,t) (with time) such that and by properly modifying the theorems in order to consider the factor V/ t.
30 C. Melchiorri (DEI) Automatic Control & System Theory 3 Lyapunov method Chetaev theorem (instability theorem) If there exists a function V(x) positive definite, with positive definite derivative in an arbitrary small neighbourhood ϑ(, ρ) of the origin, that is then the zero state x = of the system is unstable.
31 C. Melchiorri (DEI) Automatic Control & System Theory 31 Lyapunov method Chetaev theorem (instability theorem -- another formulation) Consider an autonomous dynamical systems and assume that x = is an equilibrium point. Let V : D R have the following properties: (i) V() = (ii) x R n, arbitrarily close to x =, such that V(x ) > (iii) dv/dt > x U, where the set U is defined as follows: U = {x D : x ρ, and V (x) > }. Under these conditions, x = is unstable.
32 C. Melchiorri (DEI) Automatic Control & System Theory 3 Lyapunov method Example: mass with spring and dumper m x(t) k b m ẍ = kx b ẋ Function V(x) decreases since a dissipative element ( b ) is present, otherwise V(x) would be constant and therefore the energy would not be dissipated. The system evolves until velocity reaches the zero value, and then stops (in a stable equilibrium configuration). In this configuration, from the diff. equation of the system we obtain: = k x +, x = Therefore, indeed the system stops in the zero state x =.
33 C. Melchiorri (DEI) Automatic Control & System Theory 33 Lyapunov method x =[1, ] T 1.5 m = 5 K = 5 b = V (x) V (x)
34 C. Melchiorri (DEI) Automatic Control & System Theory 34 Lyapunov method x =[1, ] T 4 m = 5 K = 5 b = V (x) V (x)
35 C. Melchiorri (DEI) Automatic Control & System Theory 35 Lyapunov method The pendulum Gravity force Viscous friction Input By considering u(t) =, the system has equilibrium states: x e,1 = apple x e, = apple m g
36 C. Melchiorri (DEI) Automatic Control & System Theory 36 Lyapunov method The pendulum Gravity force Viscous friction Input Let define the function V(x) as the sum of the potential and kinetic energy m g Clearly, V(x) is positive definite
37 C. Melchiorri (DEI) Automatic Control & System Theory 37 Lyapunov method The pendulum x 1 : posizione angolare θ(t) x 1 : angular position θ(t) x x : velocità angolare : angular velocity ω(t)
38 C. Melchiorri (DEI) Automatic Control & System Theory 38 Lyapunov method The pendulum 3 x 1 Angular position θ (t) x : angular velocity ω (t)
39 C. Melchiorri (DEI) Automatic Control & System Theory 39 Lyapunov method The pendulum The derivative of V(x) is: If d = : The system is conservative: V(x) = c, a constant depending on the initial condition V(x) = c is a closed curve about x = If c is sufficiently small, then x(t) is constrained to remain close to zero x = is then a stable equilibrium. If d : The system dissipates energy: V(x) decreases and dv(x)/dt < V(x) decreases until the system is in motion (velocity non zero) x goes to x = is then an asymptotically stable equilibrium.
40 C. Melchiorri (DEI) Automatic Control & System Theory 4 Applications of the Lyapunov method First order systems Let consider the following first order system: with f(x) continuous, and the condition: ẋ + f(x) = (1) xf(x) >, x 6= () Let consider the Lyapunov function V = x, then for The origin x = is a global asymptotically stable (GAS) point. In fact V is unlimited for x
41 C. Melchiorri (DEI) Automatic Control & System Theory 41 Applications of the Lyapunov method First order systems - Example Given the system the function f(x)= -sin x + x verifies condition () since:
42 Applications of the Lyapunov method 3 First order systems - Example 1 1 Position Velocity ẋ Phase plane x x = Plot of V(x) = x C. Melchiorri (DEI) Automatic Control & System Theory 4
43 C. Melchiorri (DEI) Automatic Control & System Theory 43 Applications of the Lyapunov method Second order systems Let consider the following second order system: where f and g are continuous functions, and the conditions:
44 C. Melchiorri (DEI) Automatic Control & System Theory 44 Applications of the Lyapunov method Second order systems Consider the Lyapunov function Then: for for If the integral function is unlimited for x, V is unlimited and the origin is GAS.
45 C. Melchiorri (DEI) Automatic Control & System Theory 45 Applications of the Lyapunov method Second order systems - Example Equations such as (3) and (4) are common when mechanical or electrical systems are considered with elements able to accumulate energy (potential or kinetic) (mass-spring-damper or LRC). Let consider the circuit shown in the figure. The equation is of type (3) and verifies (4). If x(t)=i(t), then: The integral is radially unlimited The system is GAS The Lyapunov function can be considered as the sum of the kinetic and potential energy, the fact that it decreases is related to the dissipativity of the system.
46 C. Melchiorri (DEI) Automatic Control & System Theory 46 Applications of the Lyapunov method Second order systems - Example
47 C. Melchiorri (DEI) Automatic Control & System Theory 47 Limit cycles Limit cycles and invariant sets Non-linear systems may have: 1. Equilibrium points Ø stable Ø unstable. Limit cycles Ø stable Ø unstable (closed trajectory) Definition: Limit set A positive limit set S + is a set in the state space that attracts sufficiently close trajectories when t A negative limit set S + is a set in the state space that attracts sufficiently close trajectories when τ where τ, = t An invariant J is a set defined in the state space with the property that for each initial state belonging to J, the whole trajectory belongs to J for both positive and negative values of t.
48 C. Melchiorri (DEI) Automatic Control & System Theory 48 Limit cycles Limit cycles and invariant sets Let consider a non-linear autonomous system: Assume that 1) a positive definite function V 1 exists such that in the limited domain D 1 ={ x: V 1 (x) < h 1 } the following condition holds: ) a positive definite function V exists such that in the domain D ={ x: h < V (x) < h 3 } the following condition holds: D can be unlimited: in this case the term h 3 is not considered. It is possible to show that the closed domain contains a positive limit set.
49 C. Melchiorri (DEI) Automatic Control & System Theory 49 Limit cycles D 1 ={ x: V 1 (x) < h 1 } D ={ x: h < V (x) < h 3 } x D D 1 x 1 D 3 V h h 1 V 1 D D 3 D D 3 1 D
50 Limit cycles Let consider the system It is simple to verify that a limit cycle exists for r = 1 and that such limit cycle is globally asymptotically stable for the system, that is for any initial conditions the trajectories tend to this limit cycle..8.6 V(z) X X The limit cycle r = 1 is asymptotically stable N.B.: r > in any case C. Melchiorri (DEI) Automatic Control & System Theory 5-1.5
51 Limit cycles Let consider the system ẋ 1 = x ẋ = µ (1 x 1) x x 1 ẋ 3 = x 3 Van der Pool oscillator (λ < ) C. Melchiorri (DEI) Automatic Control & System Theory 51
52 Limit cycles Evolution of the system defined by In the state space the system has an evolution, with non zero initial conditions, as shown in the figure. NB: this system is linear, then this is NOT a limit cycle! C. Melchiorri (DEI) Automatic Control & System Theory
53 Lyapunov method Non stationary linear systems Let us consider the non stationary homogeneous system: As known, this system is stable in the sense of Lyapunov if, for all t and for all ε >, a parameter η > exists such that: The system is asymptotically stable in the sense of Lyapunov if it is stable and C. Melchiorri (DEI) Automatic Control & System Theory 53
54 C. Melchiorri (DEI) Automatic Control & System Theory 54 Lyapunov method Non stationary linear systems The following theorem relates the stability properties of system (1) to the state transition matrix φ(t,t ). Theorem: System (1) is stable in the sense of Lyapunov if and only if for all t a positive real number M exists such that: System (1) is asymptoticaly stable in the sense of Lyapunov if it is stable and if:
55 C. Melchiorri (DEI) Automatic Control & System Theory 55 Lyapunov method Non stationary linear systems Let us consider the non stationary linear system : The system is B.I.B.S. if for all t and for all ε > a parameter η > exists such that if u(t) < η, t t then Theorem: System (3) is B.I.B.S. stable if and only if:
56 C. Melchiorri (DEI) Automatic Control & System Theory 56 Lyapunov method Non stationary linear systems Similarly, system (3) is B.IB.O. if for all t and for all ε > a parameter η > exists such that if u(t) < η, t t then Theorem: System (3) is B.I.B.O. stable if and only if:
57 C. Melchiorri (DEI) Automatic Control & System Theory 57 Lyapunov method Continuous-time, time invariant linear systems Let us consider the following linear time invariant system: By considering V(x) = x T P x, then: where: If matrix M is positive definite, system (1) is globally asymptotically stable. The following theorem holds: Theorem: The Lyapunov matrix equation () admits a unique solution P (symmetric positive definite matrix) for each symmetric positive definite matrix M if and only if all the eigenvalues of A are with negative real part.
58 C. Melchiorri (DEI) Automatic Control & System Theory 58 Lyapunov method Discrete time, time invariant linear systems Let us consider the following linear time invariant system: By considering V(x) = x T P x, then: where: If matrix M is positive definite, system (3) is globally asymptotically stable. The following theorem holds: Theorem: The Lyapunov matrix equation (4) admits a unique solution P (symmetric positive definite matrix) for each symmetric positive definite matrix M if and only if all the eigenvalues of A have magnitude less than 1.
59 C. Melchiorri (DEI) Automatic Control & System Theory 59 AUTOMATIC CONTROL AND SYSTEM THEORY STABILITY ANALYSIS OF DYNAMIC SYSTEMS THE END
AC&ST AUTOMATIC CONTROL AND SYSTEM THEORY SYSTEMS AND MODELS. Claudio Melchiorri
C. Melchiorri (DEI) Automatic Control & System Theory 1 AUTOMATIC CONTROL AND SYSTEM THEORY SYSTEMS AND MODELS Claudio Melchiorri Dipartimento di Ingegneria dell Energia Elettrica e dell Informazione (DEI)
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 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 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 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 informationConsider a particle in 1D at position x(t), subject to a force F (x), so that mẍ = F (x). Define the kinetic energy to be.
Chapter 4 Energy and Stability 4.1 Energy in 1D Consider a particle in 1D at position x(t), subject to a force F (x), so that mẍ = F (x). Define the kinetic energy to be T = 1 2 mẋ2 and the potential energy
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 informationBIBO STABILITY AND ASYMPTOTIC STABILITY
BIBO STABILITY AND ASYMPTOTIC STABILITY FRANCESCO NORI Abstract. In this report with discuss the concepts of bounded-input boundedoutput stability (BIBO) and of Lyapunov stability. Examples are given to
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 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 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 informationP321(b), Assignement 1
P31(b), Assignement 1 1 Exercise 3.1 (Fetter and Walecka) a) The problem is that of a point mass rotating along a circle of radius a, rotating with a constant angular velocity Ω. Generally, 3 coordinates
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 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. 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 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 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 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 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 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 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 informationSymmetries 2 - Rotations in Space
Symmetries 2 - Rotations in Space This symmetry is about the isotropy of space, i.e. space is the same in all orientations. Thus, if we continuously rotated an entire system in space, we expect the system
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 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 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 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 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 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 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 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 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 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 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 informationNonlinear Systems and Control Lecture # 19 Perturbed Systems & Input-to-State Stability
p. 1/1 Nonlinear Systems and Control Lecture # 19 Perturbed Systems & Input-to-State Stability p. 2/1 Perturbed Systems: Nonvanishing Perturbation Nominal System: Perturbed System: ẋ = f(x), f(0) = 0 ẋ
More informationDepartment of Mathematics IIT Guwahati
Stability of Linear Systems in R 2 Department of Mathematics IIT Guwahati A system of first order differential equations is called autonomous if the system can be written in the form dx 1 dt = g 1(x 1,
More 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 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 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 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 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 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 informationCONTROL OF DIGITAL SYSTEMS
AUTOMATIC CONTROL AND SYSTEM THEORY CONTROL OF DIGITAL SYSTEMS Gianluca Palli Dipartimento di Ingegneria dell Energia Elettrica e dell Informazione (DEI) Università di Bologna Email: gianluca.palli@unibo.it
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 informationSystem Control Engineering 0
System Control Engineering 0 Koichi Hashimoto Graduate School of Information Sciences Text: Nonlinear Control Systems Analysis and Design, Wiley Author: Horacio J. Marquez Web: http://www.ic.is.tohoku.ac.jp/~koichi/system_control/
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 informationControl Systems. Internal Stability - LTI systems. L. Lanari
Control Systems Internal Stability - LTI systems L. Lanari definitions (AS) - A system S is said to be asymptotically stable if its state zeroinput response converges to the origin for any initial condition
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 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 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 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 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 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 informationChapter 14 Periodic Motion
Chapter 14 Periodic Motion 1 Describing Oscillation First, we want to describe the kinematical and dynamical quantities associated with Simple Harmonic Motion (SHM), for example, x, v x, a x, and F x.
More informationThe Pendulum. The purpose of this tab is to predict the motion of various pendulums and compare these predictions with experimental observations.
The Pendulum Introduction: The purpose of this tab is to predict the motion of various pendulums and compare these predictions with experimental observations. Equipment: Simple pendulum made from string
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 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 informationAn homotopy method for exact tracking of nonlinear nonminimum phase systems: the example of the spherical inverted pendulum
9 American Control Conference Hyatt Regency Riverfront, St. Louis, MO, USA June -, 9 FrA.5 An homotopy method for exact tracking of nonlinear nonminimum phase systems: the example of the spherical inverted
More informationAnalytical Dynamics: Lagrange s Equation and its Application A Brief Introduction
Analytical Dynamics: Lagrange s Equation and its Application A Brief Introduction D. S. Stutts, Ph.D. Associate Professor of Mechanical Engineering Missouri University of Science and Technology Rolla,
More informationOscillations. Simple Harmonic Motion of a Mass on a Spring The equation of motion for a mass m is attached to a spring of constant k is
Dr. Alain Brizard College Physics I (PY 10) Oscillations Textbook Reference: Chapter 14 sections 1-8. Simple Harmonic Motion of a Mass on a Spring The equation of motion for a mass m is attached to a spring
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 informationStability of Stochastic Differential Equations
Lyapunov stability theory for ODEs s Stability of Stochastic Differential Equations Part 1: Introduction Department of Mathematics and Statistics University of Strathclyde Glasgow, G1 1XH December 2010
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 informationLagging Pendulum. Martin Ga
02 Lagging Pendulum Martin Ga 1 Task A pendulum consists of a strong thread and a bob. When the pivot of the pendulum starts moving along a horizontal circumference, the bob starts tracing a circle which
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 informationLyapunov Theory for Discrete Time Systems
Università degli studi di Padova Dipartimento di Ingegneria dell Informazione Nicoletta Bof, Ruggero Carli, Luca Schenato Technical Report Lyapunov Theory for Discrete Time Systems This work contains a
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 informationECEN 420 LINEAR CONTROL SYSTEMS. Lecture 6 Mathematical Representation of Physical Systems II 1/67
1/67 ECEN 420 LINEAR CONTROL SYSTEMS Lecture 6 Mathematical Representation of Physical Systems II State Variable Models for Dynamic Systems u 1 u 2 u ṙ. Internal Variables x 1, x 2 x n y 1 y 2. y m Figure
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 informationNonlinear Control Lecture 2:Phase Plane Analysis
Nonlinear Control Lecture 2:Phase Plane Analysis Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Fall 2010 r. Farzaneh Abdollahi Nonlinear Control Lecture 2 1/53
More 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 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 informationthe EL equation for the x coordinate is easily seen to be (exercise)
Physics 6010, Fall 2016 Relevant Sections in Text: 1.3 1.6 Examples After all this formalism it is a good idea to spend some time developing a number of illustrative examples. These examples represent
More informationRotational motion problems
Rotational motion problems. (Massive pulley) Masses m and m 2 are connected by a string that runs over a pulley of radius R and moment of inertia I. Find the acceleration of the two masses, as well as
More informationDOMAIN OF ATTRACTION: ESTIMATES FOR NON-POLYNOMIAL SYSTEMS VIA LMIS. Graziano Chesi
DOMAIN OF ATTRACTION: ESTIMATES FOR NON-POLYNOMIAL SYSTEMS VIA LMIS Graziano Chesi Dipartimento di Ingegneria dell Informazione Università di Siena Email: chesi@dii.unisi.it Abstract: Estimating the Domain
More informationAutonomous Systems and Stability
LECTURE 8 Autonomous Systems and Stability An autonomous system is a system of ordinary differential equations of the form 1 1 ( 1 ) 2 2 ( 1 ). ( 1 ) or, in vector notation, x 0 F (x) That is to say, an
More 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 informationLecture 4. Alexey Boyarsky. October 6, 2015
Lecture 4 Alexey Boyarsky October 6, 2015 1 Conservation laws and symmetries 1.1 Ignorable Coordinates During the motion of a mechanical system, the 2s quantities q i and q i, (i = 1, 2,..., s) which specify
More informationTheoretical physics. Deterministic chaos in classical physics. Martin Scholtz
Theoretical physics Deterministic chaos in classical physics Martin Scholtz scholtzzz@gmail.com Fundamental physical theories and role of classical mechanics. Intuitive characteristics of chaos. Newton
More informationDiscrete and continuous dynamic systems
Discrete and continuous dynamic systems Bounded input bounded output (BIBO) and asymptotic stability Continuous and discrete time linear time-invariant systems Katalin Hangos University of Pannonia Faculty
More informationIn the presence of viscous damping, a more generalized form of the Lagrange s equation of motion can be written as
2 MODELING Once the control target is identified, which includes the state variable to be controlled (ex. speed, position, temperature, flow rate, etc), and once the system drives are identified (ex. force,
More informationJ07M.1 - Ball on a Turntable
Part I - Mechanics J07M.1 - Ball on a Turntable J07M.1 - Ball on a Turntable ẑ Ω A spherically symmetric ball of mass m, moment of inertia I about any axis through its center, and radius a, rolls without
More informationHigh-Gain Observers in Nonlinear Feedback Control
High-Gain Observers in Nonlinear Feedback Control Lecture # 1 Introduction & Stabilization High-Gain ObserversinNonlinear Feedback ControlLecture # 1Introduction & Stabilization p. 1/4 Brief History Linear
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 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 informationInput to state Stability
Input to state Stability Mini course, Universität Stuttgart, November 2004 Lars Grüne, Mathematisches Institut, Universität Bayreuth Part IV: Applications ISS Consider with solutions ϕ(t, x, w) ẋ(t) =
More informationChapter 8. Potential Energy and Conservation of Energy
Chapter 8 Potential Energy and Conservation of Energy 8.2 Conservative and non-conservative forces A system consists of two or more particles. A configuration of the system is just a specification of the
More informationL = 1 2 a(q) q2 V (q).
Physics 3550, Fall 2011 Motion near equilibrium - Small Oscillations Relevant Sections in Text: 5.1 5.6 Motion near equilibrium 1 degree of freedom One of the most important situations in physics is motion
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 informationCDS 101 Precourse Phase Plane Analysis and Stability
CDS 101 Precourse Phase Plane Analysis and Stability Melvin Leok Control and Dynamical Systems California Institute of Technology Pasadena, CA, 26 September, 2002. mleok@cds.caltech.edu http://www.cds.caltech.edu/
More informationAutomatic Control Systems theory overview (discrete time systems)
Automatic Control Systems theory overview (discrete time systems) Prof. Luca Bascetta (luca.bascetta@polimi.it) Politecnico di Milano Dipartimento di Elettronica, Informazione e Bioingegneria Motivations
More informationExam in Systems Engineering/Process Control
Department of AUTOMATIC CONTROL Exam in Systems Engineering/Process Control 27-6-2 Points and grading All answers must include a clear motivation. Answers may be given in English or Swedish. The total
More informationarxiv: v3 [math.ds] 22 Feb 2012
Stability of interconnected impulsive systems with and without time-delays using Lyapunov methods arxiv:1011.2865v3 [math.ds] 22 Feb 2012 Sergey Dashkovskiy a, Michael Kosmykov b, Andrii Mironchenko b,
More information28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod)
28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod) θ + ω 2 sin θ = 0. Indicate the stable equilibrium points as well as the unstable equilibrium points.
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 informationLinear Differential Equations. Problems
Chapter 1 Linear Differential Equations. Problems 1.1 Introduction 1.1.1 Show that the function ϕ : R R, given by the expression ϕ(t) = 2e 3t for all t R, is a solution of the Initial Value Problem x =
More informationwhere x 0 is arbitrary.
The forces internal to a system are of two types. Conservative forces, such as gravity; and dissipative forces such as friction. Internal forces arise from the natural dynamics of the system in contract
More informationControl Systems Design
ELEC4410 Control Systems Design Lecture 14: Controllability Julio H. Braslavsky julio@ee.newcastle.edu.au School of Electrical Engineering and Computer Science Lecture 14: Controllability p.1/23 Outline
More information1 Controllability and Observability
1 Controllability and Observability 1.1 Linear Time-Invariant (LTI) Systems State-space: Dimensions: Notation Transfer function: ẋ = Ax+Bu, x() = x, y = Cx+Du. x R n, u R m, y R p. Note that H(s) is always
More informationChapter One. Introduction
Chapter One Introduction A system is a combination of components or parts that is perceived as a single entity. The parts making up the system may be clearly or vaguely defined. These parts are related
More information