Stability Analysis ofNonlinear SystemsUsing Lyapunov TheoryBy: Nafees Ahmed
Outline Motivation Definitions Lyapunov Stability Theorems Analysis of LTI System Stability Instability Theorem Examples
References Dr. Radhakant Padhi, AE Dept., IISc-Bangalore (NPTEL) H. K. Khalil: Nonlinear Systems, Prentice Hall, 1996. H. J. Marquez: Nonlinear Control Systems analysis andDesign, Wiley, 2003. J-J. E. Slotine and W. Li: Applied Nonlinear Control,Prentice Hall, 1991. Control system, principles and design by M. Gopal, McGraw Hill
Techniques of Nonlinear ControlSystems Analysis and Design Phase plane analysis: Up to 2nd order or maxi 3rd ordersystem (graphical method) Differential geometry (Feedback linearization) Lyapunov theory Intelligent techniques: Neural networks, Fuzzy logic,Genetic algorithm etc. Describing functions Optimization theory (variational optimization, dynamicprogramming etc.)
Motivation Eigenvalue analysis concept does not hold good for nonlinearsystems. Nonlinear systems can have multiple equilibrium points andlimit cycles. Stability behaviour of nonlinear systems need not be alwaysglobal (unlike linear systems). So we seek stability near theequilibrium point. Stability of non linear system depends on both initial value andits input (Unlike liner system). Stability of linear system isindependent of initial conditions. Need of a systematic approach that can be exploited forcontrol design as well.
Idea Lyapunovβs theory is based on the simpleconcept that the energy stored in a stablesystem canβt increase with time.
DefinitionsNote: Above system is an autonomous (i/p, u 0) Here Lyapunov stability is considered only for autonomous system (Itcan also extended to non autonomous system) We can have multiple equilibrium points We are interested in finding the stability at these equilibrium points Rn n dimensions (ie x1,x2 n 2 two dimensions )
DefinitionsOpen Set: Let set A be a subset of R then the set A is open if every point in Ahas a neighborhood lying in the set. Or open set means boundary lines are notincluded. Mathematically
Definitions Open set: A set π΄ βπ is called as open, if for each π₯ π΄ there exist an π 0 such thatthe interval π₯ π, π₯ π is contained in A. Such an interval is often called asπ -neighborhood of x or simply neighborhood of x.
Definitions
1. Starting with a small ball of radius Ξ΄(Ξ΅) from initialcondition Xo a system will move anywhere around the ballbut will not leave the ball of radius Ξ΅2. Ball Ξ΄(Ξ΅) is a function of Ξ΅.3. Size of Ξ΄(Ξ΅) may be larger then ball of radius Ξ΅* X0* XeΡΞ΄(Ρ)ΡΞ΄(Ρ)Ξ΄(Ρ)
Definitions
Convergent system: Starting from any initialcondition Xo, system may go anywhere but finallyconverges to equilibrium point Xe* X0* Xe
DefinitionsNote: System will never leave the Ξ΅ bound and finally will converge toequilibrium point Xe.
DefinitionsConversion :π π(π ππ )π π ππ π π ππ π π ππ 0 π π
DefinitionsA scalar function V : D R is said to be Positive definite function: if following condition aresatisfied(domain D excluding 0) Positive semi definite function: Negative define function: (i) condition same, (ii) Negative semi define function: (i) condition same, (ii) Note:1. Output of function V(x) is a scalar value, hence V(x) is scalar function .2. Negative define (semi definite) if βV(x) is definite ( semi definite)
Note:Condition (i) & (ii) V(X) positive definiteCondition (iii) π(π) Negative semi definite
What about V(X) There is no general method for selection of V(X). Some time select V(X) such that its properties are similar toenergy i.e.ππ π½ πΏ πΏπ» πΏ πΆπ π½ πΏ π²ππππππ π¬πππππ π·ππππππππ π¬πππππ πΆπ π½ πΏ πππ πππ etc How to calculate π½(πΏ)π»π» π½ π½π½ πΏ πΏ π(πΏ) π π
Note:Condition (i) & (ii) V(X) positive definiteCondition (iii) π(π) Negative definite
Radially Unbounded ? The more and more you go away from the equilibrium point, V(X) willincrease more and more.
Note: Global Subset D R
NOTE Here, pendulum with friction should beasymptotically stable as it comes to anequilibrium point finally due to friction ( π½(πΏ)should be negative definite not negative semidefinite nsdf) But we are not able to prove this. Becausex2 when x2 0, π½(πΏ) will always be βVe But when x2 0 There are multiple equilibrium pointson x1 line. Negative definite means the movement I go awayfrom the zero I should get βve valuex1
Example: Consider the system described by the equationsππ ππππ ππ πππ Solution:Chooseπ½ π πππ πππWhich satisfies following two conditions that is it ispositive definiteπ½ π π&π½ π ππ½(π) πππ ππ πππ ππ πππ ππ πππ ππ πππ πππππ½(π) π nsdf (similar to pendulum with friction)So system is stable, we canβt say asymptotically stable
Analysis of LTI system using LyapunovstabilityNote:π π΄π π π π΄ππ π π π΄π
Analysis of LTI system using Lyapunovstability
Analysis of LTI system using Lyapunovstability .
Step to solveAnalysis of LTI system using Lyapunovstability .
Example: Analysis of LTI system usingLyapunov stability Determine the stability of the system described by the following equation π₯ π΄π₯Withπ΄ 11 2 4 Solution:π΄π π ππ΄ π πΌ 1 21 π11 4 π12π12π11 π22π12π12 1 2 1 0 π22 1 40 1 Note here we took p12 p21 because Matrix P will be real symmetricmatrix
-2p11 2p12 -1 -2p11-5p12 p22 0 -4p12-8p22 -1 Solving above three equations π π11π12π12 π2223607 607 601160 which is seen to be positive definite. Hence this system is asymptoticallystable
Till now ? All were Lyapunov Directmethods There are some indirectmethods also
In rough way In rough way instability theorem state that if V(X) positive definite tπ‘ππ§ π½(πΏ) should also be positive definite
Thanks?
Techniques of Nonlinear Control Systems Analysis and Design Phase plane analysis: Up to 2nd order or maxi 3rd order system (graphical method) Differential geometry (Feedback linearization) Lyapunov theory Intelligent techniques: Neural networks, Fuzzy logic, Genetic algorithm etc. Describing functions Optimization theory (variational optimization, dynamic
Nonlinear Finite Element Analysis Procedures Nam-Ho Kim Goals What is a nonlinear problem? How is a nonlinear problem different from a linear one? What types of nonlinearity exist? How to understand stresses and strains How to formulate nonlinear problems How to solve nonlinear problems
Outline Nonlinear Control ProblemsSpecify the Desired Behavior Some Issues in Nonlinear ControlAvailable Methods for Nonlinear Control I For linear systems I When is stabilized by FB, the origin of closed loop system is g.a.s I For nonlinear systems I When is stabilized via linearization the origin of closed loop system isa.s I If RoA is unknown, FB provideslocal stabilization
function verifying maximal regions of orbital stability via iterated of sum-of-squares programs. The construction of the transverse dynamics is novel, and valid for a broad class of nonlinear hybrid systems. Keywords: Veri cation, Stability Analysis, Periodic Motion 1. INTRODUCTION Nonlinear dynamical systems exhibiting oscillating solu-
Third-order nonlinear effectThird-order nonlinear effect In media possessing centrosymmetry, the second-order nonlinear term is absent since the polarization must reverse exactly when the electric field is reversed. The dominant nonlinearity is then of third order, 3 PE 303 Ξ΅Ο The third-order nonlinear material is called a Kerr medium. P 3 E
Introduction to Nonlinear Optics 1 1.2. Descriptions of Nonlinear Optical Processes 4 1.3. Formal Deο¬nition of the Nonlinear Susceptibility 17 1.4. Nonlinear Susceptibility of a Classical Anharmonic . Rabi Oscillations and Dressed Atomic States 301 6.6. Optical Wave Mixing in Two-Level Systems 313 Problems 326 References 327 7. Processes .
oriented nonlinear analysis proceduresβ based on the so-called βpushover analysisβ. All pushover analysis procedures can be considered as approximate extensions of the response spectrum method to the nonlinear response analysis with varying degrees of sophistication. For example, βNonlinear Static ProcedureβNSPβ (ATC, 1996; FEMA, 2000) may be looked upon as a βsingle-mode .
Nonlinear analysis for improved designs Nature is nonlinear. Using Marc, accurately capture the inherent nonlinear behavior of your designs to improve product quality, reduce your testing costs, and improve reliability incorporating the true . Perform global-local analysis to better capture local behavior Dynamic analysis
eigenvalue buckling analysis, nonlinear stress analysis, and graphical post-processing. In this paper a brief description of CALEB version 1.4 and of its main features is presented. INTRODUCTION CALEB is a nonlinear finite element program for geometric and material nonlinear analysis of offshore platforms and general framed structures.