Introduction. - The Second Lyapunov Method. - The First Lyapunov Method

Similar documents
The equation of motion of a dynamical system is given by a set of differential equations. That is (1)

Dynamic Systems on Graphs

Note 10. Modeling and Simulation of Dynamic Systems

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Convexity preserving interpolation by splines of arbitrary degree

ALGORITHM FOR THE CALCULATION OF THE TWO VARIABLES CUBIC SPLINE FUNCTION

U.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017

Numerical Heat and Mass Transfer

Digital Signal Processing

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Non-negative Matrices and Distributed Control

STAT 309: MATHEMATICAL COMPUTATIONS I FALL 2018 LECTURE 16

On the symmetric character of the thermal conductivity tensor

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

Control of Uncertain Bilinear Systems using Linear Controllers: Stability Region Estimation and Controller Design

Slobodan Lakić. Communicated by R. Van Keer

10-801: Advanced Optimization and Randomized Methods Lecture 2: Convex functions (Jan 15, 2014)

Neuro-Adaptive Design - I:

Uniqueness of Weak Solutions to the 3D Ginzburg- Landau Model for Superconductivity

Stability analysis for class of switched nonlinear systems

Continuous Time Markov Chain

Improved delay-dependent stability criteria for discrete-time stochastic neural networks with time-varying delays

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS

Perfect Competition and the Nash Bargaining Solution

Appendix B. The Finite Difference Scheme

Perron Vectors of an Irreducible Nonnegative Interval Matrix

Variable Structure Control ~ Basics

The Order Relation and Trace Inequalities for. Hermitian Operators

APPENDIX A Some Linear Algebra

Group Analysis of Ordinary Differential Equations of the Order n>2

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

Decentralized Adaptive Regulation for Nonlinear Systems with iiss Inverse Dynamics

On Finite Rank Perturbation of Diagonalizable Operators

Physics 181. Particle Systems

Lecture 12: Discrete Laplacian

On the non-singularity of the thermal conductivity tensor and its consequences

On a direct solver for linear least squares problems

A p-adic PERRON-FROBENIUS THEOREM

Module 2. Random Processes. Version 2 ECE IIT, Kharagpur

ECEN 667 Power System Stability Lecture 21: Modal Analysis

Army Ants Tunneling for Classical Simulations

Maximum Likelihood Estimation of Binary Dependent Variables Models: Probit and Logit. 1. General Formulation of Binary Dependent Variables Models

Stability Analysis and Anti-windup Design of Switched Systems with Actuator Saturation

COMPLEX NUMBERS AND QUADRATIC EQUATIONS

An efficient algorithm for multivariate Maclaurin Newton transformation

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2

Projective change between two Special (α, β)- Finsler Metrics

CONDITIONS FOR INVARIANT SUBMANIFOLD OF A MANIFOLD WITH THE (ϕ, ξ, η, G)-STRUCTURE. Jovanka Nikić

THE GUARANTEED COST CONTROL FOR UNCERTAIN LARGE SCALE INTERCONNECTED SYSTEMS

Feb 14: Spatial analysis of data fields

CALCULUS CLASSROOM CAPSULES

Limit Cycle Generation for Multi-Modal and 2-Dimensional Piecewise Affine Control Systems

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction

Least Squares Fitting of Data

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Infinitely Split Nash Equilibrium Problems in Repeated Games

MEEM 3700 Mechanical Vibrations

Global Sensitivity. Tuesday 20 th February, 2018

1 Derivation of Point-to-Plane Minimization

Modeling of Dynamic Systems

PARTICIPATION FACTOR IN MODAL ANALYSIS OF POWER SYSTEMS STABILITY

The Second Anti-Mathima on Game Theory

EXPONENTIAL ERGODICITY FOR SINGLE-BIRTH PROCESSES

INTERNATIONAL JOURNAL OF ELECTRICAL ENGINEERING & TECHNOLOGY (IJEET)

Inexact Newton Methods for Inverse Eigenvalue Problems

Christian Aebi Collège Calvin, Geneva, Switzerland

SEMI-GLOBAL EXPONENTIAL STABILIZATION OF LINEAR SYSTEMS SUBJECT TO \INPUT SATURATION" VIA LINEAR FEEDBACKS. Zongli Lin & Ali Saberi

Randić Energy and Randić Estrada Index of a Graph

Lecture 21: Numerical methods for pricing American type derivatives

Constant Best-Response Functions: Interpreting Cournot

Neuro-Adaptive Design II:

Stability and Analytical Approximation of Limit. Cycles in Hopf Bifurcations of Four-dimensional. Economic Models

MEM 255 Introduction to Control Systems Review: Basics of Linear Algebra

A PROCEDURE FOR SIMULATING THE NONLINEAR CONDUCTION HEAT TRANSFER IN A BODY WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY.

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed

LocalReliableControlforLinearSystems with Saturating Actuators

A new Approach for Solving Linear Ordinary Differential Equations

SELECTED SOLUTIONS, SECTION (Weak duality) Prove that the primal and dual values p and d defined by equations (4.3.2) and (4.3.3) satisfy p d.

Nonlinear Oscillations in a three-dimensional Competition with Inhibition Responses in a Bio-Reactor

ON THE BURGERS EQUATION WITH A STOCHASTIC STEPPING STONE NOISY TERM

Visualization of the Economic Impact of Process Uncertainty in Multivariable Control Design

6) Derivatives, gradients and Hessian matrices

Lecture 6/7 (February 10/12, 2014) DIRAC EQUATION. The non-relativistic Schrödinger equation was obtained by noting that the Hamiltonian 2

Maximum Likelihood Estimation of Binary Dependent Variables Models: Probit and Logit. 1. General Formulation of Binary Dependent Variables Models

Lecture 5.8 Flux Vector Splitting

A RELAXED SUFFICIENT CONDITION FOR ROBUST STABILITY OF AFFINE SYSTEMS

The Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method

Three views of mechanics

Modelli Clamfim Equazioni differenziali 7 ottobre 2013

CALCU THIRD EDITION. Monfy J. Strauss. Texas Tech University. Gerald L Bradley. Claremont McKenna College. Karl J. Smith. Santa Rosa Junior College

Fundamental loop-current method using virtual voltage sources technique for special cases

Mathematical Preparations

The Two-scale Finite Element Errors Analysis for One Class of Thermoelastic Problem in Periodic Composites

6.854J / J Advanced Algorithms Fall 2008

Linear Approximation with Regularization and Moving Least Squares

CHEMICAL REACTIONS AND DIFFUSION

Georgia Tech PHYS 6124 Mathematical Methods of Physics I

Transcription:

Stablty Analyss A. Khak Sedgh Control Systems Group Faculty of Electrcal and Computer Engneerng K. N. Toos Unversty of Technology February 2009 1

Introducton Stablty s the most promnent characterstc of dynamcal systems Dynamcal Systems Stablty from: NLTV to LTV and LTI Stablty of LTI dynamcal systems: The Egenvalue concept Methods of examnng LTI dynamcal System Stablty Lyapunov : Stablty analyss from NLTV to LTV and LTI Two man Lyapunov Methods: - The Frst Lyapunov Method - The Second Lyapunov Method 2

Stablty Defntons x() t = f [ x(), t u(),] t t Equlbrum Pont Concept 3

Defnton Stable Equlbrum Pont n the sense of Lyapunov Defnton Asymptotcally Stable Equlbrum Pont 4

Defnton Doman of Attracton Defnton Globally Asymptotcal Stable Equlbrum Pont Defnton Unstable Equlbrum Pont Defnton Internal Stablty Defnton BIBO Stable Defnton T-Stable or Totally Stable 5

LTI System Stablty ( ) x () t = Ax t + Bu() t y() t = Cx() t Autonomous Systems wth Dstnct Egenvalues Yelds: x t ( ) n = µ = 1 Theorem A System s Stable n Lyapunov Sense f and only f: - All Egenvalues have Non-Postve Real Parts - Complex Egenvalues are Smple e λ t e 6

Theorem An LTI system s Asymptotcally Stable f and only f t s Globally Asymptotcally Stable. Theorem An LTI system s Asymptotcally Stable f and only f t s Egenvalues have Strctly Negatve Real Parts. Theorem An LTI system s BIBO (Bounden Input-Bounded Output) Stable f and only f All transfer Functon Poles have Strctly Negatve Real Parts. 7

Theorem If an LTI system s Controllable and Observable, then the followng terms are Equvalent - System s Totally Stable. - Systems Zero State Response s BIBO Stable. - Transfer Functon Poles have Negatve Real Parts. - Egenvalues of the State Matrx have Negatve Real Parts. 8

Lyapunov s Frst Method Stablty Analyss of Nonlnear Systems n Operatng Ponts After Lnearzaton Theorem System x( t) = f [ x( t ),0] and x e x () t = Ax () t Where, A = Jacoban of f at x e That s, State Matrx Asymptotc Stablty yelds Nonlnear System Asymptotc Stablty n a specfc operatng pont. 9

Lyapunov s Second Method A General Stablty Analyss of Dynamcal Systems from an Internal System Pont of vew. Lyapunov Theory s Common Practce n Analyss and Desgn of Dynamcal Systems. Lyapunov Canddate Functon. Exstence of Lyapunov Functon. Unqueness of Lyapunov Functon? 10

Mathematcal Prelmnares Quadratc Forms V ( x) n n = = 1 j= 1 a x x j j x = x1 x n (,, ) [ ] V ( x) = x x x 1 2 T x Ax = = x, Ax n aj R a11 a12 a1 n x1 a21 a22 a 2n x 2 a a a x n1 n2 nn n Symmetrc Matrx 11

Defnton Postve Defnte Scalar Functon V ( x) n vcnty of S Defnton Postve Sem Defnte Scalar functon V ( x) n vcnty of S Defnton Negatve Defnte Scalar Functon V ( x) n vcnty of S Defnton Negatve Sem Defnte Scalar functon V ( x) n vcnty of S Defnton Indefnte Scalar functon V ( x) n vcnty of S Scalar Functon Sgn Determnaton Methods. 12

Lyapunov s Second Method Classcal Mechancal Theory Prncpal: Oscllatng Systems wth No Stmulatng Input are Stable f ther Total Energy s Contnuously Decreasng. Lyapunov Theory Based on Energy Functon Generalzed Energy Functon or Lyapunov Functon Lyapunov Canddate Functon Propertes of Lyapunov Canddate Functon x = ( x1,, x n ) V ( x) ( ) dv x dt 13

Lyapunov Functon Illustraton 14

Theorem System s Asymptotcally Stable n Localty of Equlbrum Pont n orgn f there exsts a Scalar Functon such that: V ( x) > 0 for x 0 V (0) = 0 dv ( x ) dt < 0 for x 0 15

What f Lyapunov Functon s Dervatve s Negatve Sem Defnte? Do System Trajectores Zero Lyapunov Functon s Dervatve? Revsed Asymptotc Stablty Condton V ( x) 0, V ( x ) = 0 x f ( x( t)) 16

Theorem A system s Asymptotcally Stable n vcnty of Orgn f t has only One Equlbrum Pont, Furthermore there exsts a Scalar Functon so that: Functon s Contnuous on entre State Space and ts Partal Dervatves are also Contnuous V ( x) > 0 for x 0 V (0) = 0 V ( x) for V x ( x) 0, V ( x ) = 0 x f ( x( t)) 17

Theorem (System Instablty) A system s unstable n the vcnty of orgn f there exsts a scalar functon that: V ( x) 0, V (0) = 0, V ( x) Contnuous n S wth contnuous partal dervatves V ( x) > 0, V (0) = 0 18

Lyapunov Stablty Analyss for Lnear Tme Invarant Systems Lnear Tme Invarant System: x() t = Ax() t Necessary and Suffcent condtons for Lnear Tme Invarant Stablty based on Egenvalues and Characterstc Equaton Lyapunov Algebrac Method for Stablty Analyss of Lnear Tme Invarant Systems 19

General Structure T x() t = Ax() t and V ( x) = x Px T T V ( x) = x Px + x P x T T = ( Ax) Px + x P( Ax) T T = x ( A P + PA) x = T x Qx Symmetrc Matrx Symmetrc Matrx T A P + PA = Q Lyapunov Equaton 20

Stablty Evaluaton Usng Lyapunov s Second Method Step 1: PD or PSD Matrx Selecton Step 2: Solvng the Lyapunov Equaton Step 3: Sgn Determnaton of the Matrx Step 4: Stablty Evaluaton from the Sgn of the Matrx Q A T P + PA = Q P Matrx n Step 1 PD or PSD? 21

Theorem The followng system x () t = Ax() t s asymptotcally stable f and only f Q PD Matrx P PD Matrx Proof See Ref [1] PSD Case? 22