Exponential Stability of Conformable Fractional Differential Equations

Total Page:16

File Type:pdf, Size:1020Kb

Exponential Stability of Conformable Fractional Differential Equations 2017 International Conference on Electrical Engineering and Automation Control (ICEEAC 2017) ISBN: 978-1-60595-447-9 Exponential Stability of Conformable Fractional Differential Equations Zhong-Qian WANGa and Ying SUNb School of Mathematics, University of Jinan, Jinan, China [email protected], [email protected] Keywords: Fractional Differential Equation; Laplace Transform,; Conformable Fractional Calculus Abstract. In this paper, by utilizing Laplace transforms and Lyapunove functions, we investigate the global ultimate boundedness and exponential stability for some impulsive conformable fractional differential equations. As an application, an example is given to demonstrate the main results. Introduction Fractional differential equation is the theory of the differential and integral of any order, it is an extension of the integer order differential equation. Fractional calculus is nowadays one of the most intensively developing areas of mathematical analysis including several definitions of fractional operators like riemann-Liouville, Caputo and Grunwald-Letnikov, for details, see [1-4] and reference therein. Very recently, a new fractional derivative called the conformal be fractional derivative is given in [5, 6]. Some studies about theory and applications of the conformable fractional differential equations were conducted [7-10]. In this paper, we propose the application of Lyapunove's method to nonlinear fractional differential equations thanks to the definition of Laplace transforms of conformable fractional derivative, and give the global ultimate boundedness of such systems. The article is organized as follows. In Section 2 we recall the main concepts and tools necessary in the sequel. Our results are given in Section3. An example is given to demonstrate the main result in Section 4. Preliminaries Let RRn( n× n ) be the set of n( n× n) -dimensional real vectors (matrices). Let ||⋅ || denote the Euclidean norm in Rn , R =[0, ∞ ) , R=[,) t ∞ , and N ={1,2,3, …} . C RRR× n , denote the set of all + t0 0 ( t0 + ) piecewise continuous functions from RR× n to R . τ ()⋅ and τ ()⋅ denote the maximum and the t0 + M m minimum eigenvalue of the corresponding matrix, respectively. For convenience, we give some useful definitions and facts, which can be found in [5, 6]. Definition 2.1 Given a function f:[,) a∞ → R . Then the Conformable fractional derivative of order 0<α ≤ 1 is defined by 1−α a f(())() t+ε t − a − f t Tα f( t )= lim . (1) ε →0 ε Ta f t Ta f a= lim Ta f t If α () exists on (a,b) then α () t→ a+ α () . Definition 2.2 Let a∈ R,α ∈ (0,1] , and f :[ a,∞) → R be real valued function. Then the fractional Laplace transform of order α starting from a of f is defined by 520 a Lα { f( t)}(s) = F a ()s α (2) ()t−a α ∞ −s = eα f()() t t− aα −1 dt ∫a Lemma 2.1 Let a∈ R,α ∈ (0,1] and f:(,) a∞ → R be differentiable real valued function. Then a a Lα{( T α f)( t)}(s) = s Fα ( s) − f( a) . (3) a a Lemma 2.2 let f:(,) a∞ → R be a function such that Lα { f ( t )}( s )= Fα ( s) exists. Then a 1/α Fα (s) =L { f ( a + (α t ) )}( s ), (4) ∞ where L{g ( t )}( s )= e−st g ( t ) dt is the usual Laplace transform. ∫0 After a simple calculation, we see that t 1 L0 ()()1 s=, s > 0. α s 1 1 1+ t t0 α L 0 t s = +α α ,s > 0. α ()() 1 s 1+ s α α ()t−t0 t λ 1 L0 e α = . α s − λ Consider the following systems: Tα u( t) = Au( t) + g( t,, u( t)) u t u t+ u t − J t − u t − ∆()k =()()() k − k = k() k,, k (5) u t+ u ()0= 0 n× n where 0<α < 1, Tα u() t is the conformable fractional derivative of u() t of order α , A ∈ R is a n n n n constant matrix, g: Rt × R → R and JRRk : → , tk ( k∈ N ) satisfy t1< t 2 < < tk < and limtk = ∞ , 0 k→∞ + − with u( tk ) =lim u() tk + ε , u( tk ) =lim u() tk + ε . ε →0+ ε →0− n In this paper, we assume that for any given u0 ∈ R , system (5) owns at least one solution which is u t t u u t t u t− u t denoted by ( ,,0 0 ) (simply () ) and is left continuous at each k , that is ( k) = () k . For results on the existence and uniqueness of the solution of IFDEs, one can refer to [5-10]. Main Results Theorem 3.1 Suppose there exist a function P t, u∈ C R × Rn , R and several constants () ( t0 + ) d≥0, a1 > 0, a 2 > 0,ωk > 0 and ρ > 0 such that (I) for all (t, u) ∈ R × Rn , t0 2 P 2 a1 || g ||≤( t , u) ≤ a2 || g || ; (6) 521 (II) For all k∈ N , and u∈ Rn , + P (tk, u()() t k+ J k( t k, u t k)) ≤ γ kP ( t k,; u() t k ) (7) n (III) for all k∈ N, t ∈ ∆k−1( t k − 1 ,] t k , and u∈ R , ()(,)(,)Tα PP t u≤ −ρ t u + d (8) (IV) for all i∈ N, letγ= sup{ γ i} < 1 , δk−1 =t k − t k −1 . Then system (5) is global ultimately bounded, and the solution u (t) converges to the compact set P ={u ∈ Rn || u() t | ≤ G} , (9) ρ d − δ α with G=(1 − e α ). ρ− ργ Proof: According to (8), we can choose a nonnegative function N( t) such that (Tα PP )(,()) tut+ρ (,) tuNtdt + () = , ∈ ∆k , kN ∈ . (10) Laplace transform of (10) yields d sPP( s )− ( t+ , u ( t + ))+ ρPN ( s )+ ( s ) = , α k−1 k − 1 α α s t t k−1 N k−1 where PPα()s = L α { ( t, u() t )} and α(s) = L α { N( t)} . Writing Pα (s) in the form d P t+, u t + −N () s + ( k() k −1 ) α s P ()s = . (11) s + ρ The inverse Laplace transform of (11) gives P (t, u( t)) α α −ρβk−1 −ρβk−1 P t+ u t+ eα N t e α = k−1, () k − 1 − ()* () (12) α βk d −ρ −1 α + 1−e , k ∈ N , t ∈ ∆k , ρ where βk−1=t − t k − 1 , ∗ is the convolution operator. Using the no negativity, we obtain β α −ρ k−1 PPt u t t+ u t+ e α (),() ≤ ()k−1, () k − 1 α βk (13) d −ρ −1 +1 −eα , k ∈ N , t ∈ ∆ , ρ k ρ Taking k=1 in (\ref{eq13}) and set = ω , we get α 522 −ωβα PP 0 (t, u() t) ≤ ( t0, u() t 0 ) e d −ωβα (14) +1 −e0 , t ∈ ∆ ρ ( ) 1 Using (7) and (13), for t ∈ ∆2 we have −ω βα + δ α PP ()1 0 (,t u() t≤ γ1 ( t0, u() t 0 ) e dγ −ωδ α−ωβ αd −ωβ α (15) + 1 (1− e0 ) e1+ (1 − e 1 ) ρ ρ Repeating the above procedure, for t ∈ ∆k , we have P(t , u ( t )) k −ωδα − ωβ α i−1 k P ≤ ∏γ ie e( t0, u ( t 0 )) i=1 k −1 j α dγ k j α α (16) −ωδi+ k − j − −ωβk −ρδk− j + { γ i+ k − je e(1− e )} ∑ ∏ ρ j=1 i=1 d α αd α −ωδk−1 − ωβ k −ωβ k +γ k (1 − e ) e + (1 − e ). ρ ρ Setting δ= sup{ δk } and considering γ= sup{ γ i } , βk < δ and the no negativity of some terms, from (16) we have PP k− kωδ α (t, u() t) ≤ ( t0, u() t 0 )γ e k+1 α 1−γ d (17) + (1− e−ωδ ) 1− γ ρ This together with Condition (i), we have ‖‖u( t) k +1 a α α d 2‖‖2 −kωδ −ωδ (γ − 1) ≤ u0 γ k e+(1 − e ) , (18) a1 ρ( γ − 1) t∈ ∆k +1, k ≥ 0 which implies the conclusion holds for γ < 1 . Corollary 3.1 If d = 0 and the conditions (i)-(iv) of Theorem 3.2 hold, then the system (5) is globally exponetially stable with the exponential convergence rate α 1 eωδ λ = ln( ), (19) 2δ γ where δ= sup{} δk . k∈ N Proof. If d = 0 , then from (18) we have a α ‖ ‖‖‖2 2 −kωδ u ()t≤ u0 γ k e , t ∈ ∆k+1 , k ≥ 0 . (20) a1 523 For t∈ ∆k+1, k ≥ 0 , we have t− t0 ≤ tk+ 1 − t 0 ≤( k +1)δ . (21) Using above two inequalities (20) and (21), we have ‖u( t)‖ a α 2‖ ‖2 −ωδ k ≤ u0 ()γ e a1 t− t0 a α α 2‖‖2 −ωδ −1 −ωδ 2δ ≤ u0 ()()γ eγ e a1 ωδ α a e −t − t ⋅λ 2‖‖2 ()0 = u0 e a1 γ α 1 eωδ where λ = ln( ).Which ends the proof of Corollary 3.1. 2δ γ Lemma 3.2 Let y() t∈ Rn be a vector of differential function and B ∈ Rn× n a symmetric and positive definite matrix. Then T yT tBB y t yT t T y t ()α ( ()()) = 2 ( )()α ( ), (22) for all 0<α < 1,t ≥ t0. Proof. Since matrix B is symmetric, there exist an orthogonal matrix H ∈ Rn× n and a diagonal matrix Λ ∈ Rn× n such that yT ()() tB y t TT =y ()()t H Λ H y t (23) T =()HTT y() t Λ () H y() t . Let z( t) = HT y( t), from (23) we have T yT ()()()() tB y t= z t Λ z t .
Recommended publications
  • Exponential Stability for the Resonant D’Alembert Model of Celestial Mechanics
    DISCRETE AND CONTINUOUS Website: http://aimSciences.org DYNAMICAL SYSTEMS Volume 12, Number 4, April 2005 pp. 569{594 EXPONENTIAL STABILITY FOR THE RESONANT D'ALEMBERT MODEL OF CELESTIAL MECHANICS Luca Biasco and Luigi Chierchia Dipartimento di Matematica Universit`a\Roma Tre" Largo S. L. Murialdo 1, 00146 Roma (Italy) (Communicated by Carlangelo Liverani) Abstract. We consider the classical D'Alembert Hamiltonian model for a rotationally symmetric planet revolving on Keplerian ellipse around a ¯xed star in an almost exact \day/year" resonance and prove that, notwithstanding proper degeneracies, the system is stable for exponentially long times, provided the oblateness and the eccentricity are suitably small. 1. Introduction. Perturbative techniques are a basic tool for a deep understand- ing of the long time behavior of conservative Dynamical Systems. Such techniques, which include the so{called KAM and Nekhoroshev theories (compare [1] for gen- eralities), have, by now, reached a high degree of sophistication and have been applied to a great corpus of di®erent situations (including in¯nite dimensional sys- tems and PDE's). KAM and Nekhoroshev techniques work under suitable non{ degeneracy assumptions (e.g, invertibility of the frequency map for KAM or steep- ness for Nekhoroshev: see [1]); however such assumptions are strongly violated exactly in those typical examples of Celestial Mechanics, which { as well known { were the main motivation for the Dynamical System investigations of Poincar¶e, Birkho®, Siegel, Kolmogorov, Arnold, Moser,... In this paper we shall consider the day/year (or spin/orbit) resonant planetary D'Alembert model (see [10]) and will address the problem of the long time stability for such a model.
    [Show full text]
  • Adding Integral Action for Open-Loop Exponentially Stable Semigroups and Application to Boundary Control of PDE Systems A
    1 Adding integral action for open-loop exponentially stable semigroups and application to boundary control of PDE systems A. Terrand-Jeanne, V. Andrieu, V. Dos Santos Martins, C.-Z. Xu Abstract—The paper deals with output feedback stabilization of boundary regulation is considered for a general class of of exponentially stable systems by an integral controller. We hyperbolic PDE systems in Section III. The proof of the propose appropriate Lyapunov functionals to prove exponential theorem obtained in the context of hyperbolic systems is given stability of the closed-loop system. An example of parabolic PDE (partial differential equation) systems and an example of in Section IV. hyperbolic systems are worked out to show how exponentially This paper is an extended version of the paper presented in stabilizing integral controllers are designed. The proof is based [16]. Compare to this preliminary version all proofs are given on a novel Lyapunov functional construction which employs the and moreover more general classes of hyperbolic systems are forwarding techniques. considered. Notation: subscripts t, s, tt, ::: denote the first or second derivative w.r.t. the variable t or s. For an integer n, I I. INTRODUCTION dn is the identity matrix in Rn×n. Given an operator A over a ∗ The use of integral action to achieve output regulation and Hilbert space, A denotes the adjoint operator. Dn is the set cancel constant disturbances for infinite dimensional systems of diagonal matrix in Rn×n. has been initiated by S. Pojohlainen in [12]. It has been extended in a series of papers by the same author (see [13] II.
    [Show full text]
  • Exponential Stability of Distributed Parameter Systems Governed by Symmetric Hyperbolic Partial Differential Equations Using Lyapunov’S Second Method ∗
    ESAIM: COCV 15 (2009) 403–425 ESAIM: Control, Optimisation and Calculus of Variations DOI: 10.1051/cocv:2008033 www.esaim-cocv.org EXPONENTIAL STABILITY OF DISTRIBUTED PARAMETER SYSTEMS GOVERNED BY SYMMETRIC HYPERBOLIC PARTIAL DIFFERENTIAL EQUATIONS USING LYAPUNOV’S SECOND METHOD ∗ Abdoua Tchousso1, 2, Thibaut Besson1 and Cheng-Zhong Xu1 Abstract. In this paper we study asymptotic behaviour of distributed parameter systems governed by partial differential equations (abbreviated to PDE). We first review some recently developed results on the stability analysis of PDE systems by Lyapunov’s second method. On constructing Lyapunov functionals we prove next an asymptotic exponential stability result for a class of symmetric hyper- bolic PDE systems. Then we apply the result to establish exponential stability of various chemical engineering processes and, in particular, exponential stability of heat exchangers. Through concrete examples we show how Lyapunov’s second method may be extended to stability analysis of nonlin- ear hyperbolic PDE. Meanwhile we explain how the method is adapted to the framework of Banach spaces Lp,1<p≤∞. Mathematics Subject Classification. 37L15, 37L45, 93C20. Received November 2, 2005. Revised September 4, 2006 and December 7, 2007. Published online May 30, 2008. 1. Introduction Lyapunov’s second method, called also Lyapunov’s direct method, is one of efficient and classical methods for studying asymptotic behavior and stability of the dynamical systems governed by ordinary differential equations. It dates back to Lyapunov’s thesis published at the beginning of the twentieth century [22]. The applications of the method and its theoretical development have been marked by the works of LaSalle and Lefschetz in the past [19,35] and numerous successful achievements later in the control of mechanical, robotic and aerospace systems (cf.
    [Show full text]
  • Stability Analysis for a Class of Partial Differential Equations Via
    1 Stability Analysis for a Class of Partial terms of differential matrix inequalities. The proposed Differential Equations via Semidefinite method is then applied to study the stability of systems q Programming of inhomogeneous PDEs, using weighted norms (on Hilbert spaces) as LF candidates. H Giorgio Valmorbida, Mohamadreza Ahmadi, For the case of integrands that are polynomial in the Antonis Papachristodoulou independent variables as well, the differential matrix in- equality formulated by the proposed method becomes a polynomial matrix inequality. We then exploit the semidef- Abstract—This paper studies scalar integral inequal- ities in one-dimensional bounded domains with poly- inite programming (SDP) formulation [21] of optimiza- nomial integrands. We propose conditions to verify tion problems with linear objective functions and poly- the integral inequalities in terms of differential matrix nomial (sum-of-squares (SOS)) constraints [3] to obtain inequalities. These conditions allow for the verification a numerical solution to the integral inequalities. Several of the inequalities in subspaces defined by boundary analysis and feedback design problems have been studied values of the dependent variables. The results are applied to solve integral inequalities arising from the using polynomial optimization: the stability of time-delay Lyapunov stability analysis of partial differential equa- systems [23], synthesis control laws [24], [29], applied to tions. Examples illustrate the results. optimal control law design [15] and system analysis [11]. Index Terms—Distributed Parameter Systems, In the context of PDEs, preliminary results have been PDEs, Stability Analysis, Sum of Squares. presented in [19], where a stability test was formulated as a SOS programme (SOSP) and in [27] where quadratic integrands were studied.
    [Show full text]
  • Introduction to Lyapunov Stability Analysis
    Tools for Analysis of Dynamic Systems: Lyapunov’s Methods Stanisław H. Żak School of Electrical and Computer Engineering ECE 680 Fall 2013 1 A. M. Lyapunov’s (1857--1918) Thesis 2 Lyapunov’s Thesis 3 Lyapunov’s Thesis Translated 4 Some Details About Translation 5 Outline Notation using simple examples of dynamical system models Objective of analysis of a nonlinear system Equilibrium points Lyapunov functions Stability Barbalat’s lemma 6 A Spring-Mass Mechanical System x---displacement of the mass from the rest position 7 Modeling the Mass-Spring System Assume a linear mass, where k is the linear spring constant Apply Newton’s law to obtain Define state variables: x1=x and x2=dx/dt The model in state-space format: 8 Analysis of the Spring-Mass System Model The spring-mass system model is linear time-invariant (LTI) Representing the LTI system in standard state-space format 9 Modeling of the Simple Pendulum The simple pendulum 10 The Simple Pendulum Model Apply Newton’s second law J mglsin where J is the moment of inertia, 2 J ml Combining gives g sin l 11 State-Space Model of the Simple Pendulum Represent the second-order differential equation as an equivalent system of two first-order differential equations First define state variables, x1=θ and x2=dθ/dt Use the above to obtain state–space model (nonlinear, time invariant) 12 Objectives of Analysis of Nonlinear Systems Similar to the objectives pursued when investigating complex linear systems Not interested in detailed solutions, rather one seeks to characterize the system behavior---equilibrium points and their stability properties A device needed for nonlinear system analysis summarizing the system behavior, suppressing detail 13 Summarizing Function (D.G.
    [Show full text]
  • Small Denominators and Exponential Stability: from Poincare´ to the Present Time
    SMALL DENOMINATORS AND EXPONENTIAL STABILITY: FROM POINCARE´ TO THE PRESENT TIME ANTONIO GIORGILLI Dipartimento di Matematica dell’ Universit`adi Milano Via Saldini 50, 20133 Milano, Italy. Abstract. The classical problem of small denominators is revisited in its historical develop- ment, ending with recent results on exponential stability. Conference given at the Seminario Matematico e Fisico di Milano. Small denominators and exponential stability 1 1. Overview In 1959 J.E. Littlewood published two papers devoted to the stability of the Lagrangian equilateral equilibria of the problem of three bodies.[30][31] Le me quote the incipit of the first of those papers: “ The configuration is one of two point masses S and J, with mJ = µmS, and a body P of zero mass, the law of attraction being the inverse square; the equilateral triangle SJP , of side 1, rotates in equilibrium with unit angular velocity about the centre of gravity of S and J. If µ < µ0, where µ0 is the −2 1 smaller root of µ0(1 + µ0) = 27 , the configuration is stable in the sense that a small disturbance does not double in a few revolutions; there are real periods 2π/λ1, 2π/λ2 of ‘ normal oscillations ’. It is notorious that no non–trivial system whatever is known to be sta- ble (or bounded) over infinite time. It is possible, however, to ask a less far–reaching question: given that the initial disturbance in coordinates and velocities is of order ε, for how long a time, in terms of ε, can we say the disturbance will remain of order ε? We shall find that, for almost all values −1/2 −3/4 of µ < µ0, this time is as long as exp Aε log ε , where A depends only on µ; while not eternity, this is a considerable| | slice of it.
    [Show full text]
  • Existence and Exponential Stability for Neutral Stochastic Integro–Differential Equations with Impulses Driven by a Rosenblatt Process
    DISCRETE AND CONTINUOUS doi:10.3934/dcdsb.2019251 DYNAMICAL SYSTEMS SERIES B Volume 25, Number 2, February 2020 pp. 507{528 EXISTENCE AND EXPONENTIAL STABILITY FOR NEUTRAL STOCHASTIC INTEGRO{DIFFERENTIAL EQUATIONS WITH IMPULSES DRIVEN BY A ROSENBLATT PROCESS Tomas´ Caraballo Dpto. Ecuaciones Diferenciales y An´alisisNum´erico Facultad de Matem´aticas,Universidad de Sevilla C/ Tarfia s/n, 41012-Sevilla, Spain Carlos Ogouyandjou and Fulbert Kuessi Allognissode Universit´ed'Abomey-Calavi(UAC) Institut de Math´ematiqueset de Sciences Physiques(IMSP) 01 B.P. 613, Porto-Novo, R´epubliquedu B´enin Mamadou Abdoul Diop Universit´eGaston Berger de Saint-Louis, UFR SAT D´epartement de Math´ematiques B.P234, Saint-Louis, S´en´egal Dedicated to Prof. Dr. Juan J. Nieto on the occasion of his 60th birthday Abstract. The existence and uniqueness of mild solution of an impulsive sto- chastic system driven by a Rosenblatt process is analyzed in this work by using the Banach fixed point theorem and the theory of resolvent operator developed by R. Grimmer in [12]. Furthermore, the exponential stability in mean square for the mild solution to neutral stochastic integro–differential equations with Rosenblatt process is obtained by establishing an integral inequality. Finally, an example is exhibited to illustrate the abstract theory. 1. Introduction. Integro-differential equations are of great importance in the mod- eling of several physical phenomena. Their resolution can be done through the theory of resolvent operators (see Grimmer [12]). However, the resolvent operator does not satisfy semigroup properties. The study of the quantitative and qualita- tive properties of solutions to stochastic neutral differential equations like existence, uniqueness and stability, have been widely examined by many researchers by analyz- ing various mathematical models in different areas such as mechanics, electronics, control theory, engineering and economics, etc (see [3,5,7,8, 25]).
    [Show full text]
  • Nonoscillation, Maximum Principles, and Exponential Stability of Second Order Delay Differential Equations Without Damping Term Alexander Domoshnitsky*
    Domoshnitsky Journal of Inequalities and Applications 2014, 2014:361 http://www.journalofinequalitiesandapplications.com/content/2014/1/361 RESEARCH Open Access Nonoscillation, maximum principles, and exponential stability of second order delay differential equations without damping term Alexander Domoshnitsky* *Correspondence: [email protected] Ariel University, Ariel, Israel Abstract Delays, arising in nonoscillatory and stable ordinary differential equations, can induce oscillation and instability of their solutions. That is why the traditional direction in the study of nonoscillation and stability of delay equations is to establish a smallness of delay, allowing delay differential equations to preserve these convenient properties of ordinary differential equations with the same coefficients. In this paper, we find cases in which delays, arising in oscillatory and asymptotically unstable ordinary differential equations, induce nonoscillation and stability of delay equations. We demonstrate that, although the ordinary differential equation x(t)+c(t)x(t) = 0 can be oscillating and asymptotically unstable, the delay equation x(t)+a(t)x(t – h(t)) – b(t)x(t – g(t)) = 0, where c(t)=a(t)–b(t), can be nonoscillating and exponentially stable. Results on nonoscillation and exponential stability of delay differential equations are obtained. On the basis of these results on nonoscillation and stability, the new possibilities of non-invasive (non-evasive) control, which allow us to stabilize a motion of single mass point, are proposed. Stabilization of this
    [Show full text]
  • 4 Lyapunov Stability Theory
    4 Lyapunov Stability Theory In this section we review the tools of Lyapunov stability theory. These tools will be used in the next section to analyze the stability properties of a robot controller. We present a survey of the results that we shall need in the sequel, with no proofs. The interested reader should consult a standard text, such as Vidyasagar [?] or Khalil [?], for details. 4.1 Basic definitions Consider a dynamical system which satisfies n x˙ = f(x, t) x(t0)=x0 x ∈ R . (4.31) We will assume that f(x, t) satisfies the standard conditions for the exis- tence and uniqueness of solutions. Such conditions are, for instance, that f(x, t) is Lipschitz continuous with respect to x, uniformly in t, and piece- wise continuous in t.Apointx∗ ∈ Rn is an equilibrium point of (4.31) if f(x∗,t) ≡ 0. Intuitively and somewhat crudely speaking, we say an equi- librium point is locally stable if all solutions which start near x∗ (meaning that the initial conditions are in a neighborhood of x∗) remain near x∗ for all time. The equilibrium point x∗ is said to be locally asymptotically stable if x∗ is locally stable and, furthermore, all solutions starting near x∗ tend towards x∗ as t →∞. We say somewhat crude because the time-varying nature of equation (4.31) introduces all kinds of additional subtleties. Nonetheless, it is intuitive that a pendulum has a locally sta- ble equilibrium point when the pendulum is hanging straight down and an unstable equilibrium point when it is pointing straight up.
    [Show full text]
  • Module 06 Stability of Dynamical Systems
    Stability of CT Systems Stability of DT Linear Systems Lyapunov Stability, Nonlinear Systems Module 06 Stability of Dynamical Systems Ahmad F. Taha EE 5143: Linear Systems and Control Email: [email protected] Webpage: http://engineering.utsa.edu/ataha October 10, 2017 ©Ahmad F. Taha Module 06 — Stability of Dynamical Systems 1 / 24 Stability of CT Systems Stability of DT Linear Systems Lyapunov Stability, Nonlinear Systems Stability of CT LTV Systems The following CT LTI system without inputs n x˙ (t) = A(t)x(t), x(t) ∈ R has an equilibrium at xe = 0. Asymptotic Stability The above system is asymptotically stable at xe = 0 if its solution x(t) starting from any initial condition x(t0) satisfies x(t) → 0, as t → ∞ Exponential Stability The above system is exponentially stable at xe = 0 if its solution x(t) starting from any initial condition x(t0) satisfies −rt kx(t)k ≤ Ke kx(t0)k, ∀t ≥ t0 for some positive constants K and r. ©Ahmad F. Taha Module 06 — Stability of Dynamical Systems 2 / 24 Stability of CT Systems Stability of DT Linear Systems Lyapunov Stability, Nonlinear Systems Example 1 Consider the following TV LTI system from Homework 4: − 1 0 x˙ (t) = t+1 x(t) − 1 0 t+1 Recall that the solution to this system is 1 0 1 1 x(t) = φ(t, 0)x(0) = t+1 = t+1 − t 1 −1/5 − t − 1/5 t+1 t+1 Is this system asymptotically stable? Solution: it’s not, since the states do not go to zero for any initial conditions ©Ahmad F.
    [Show full text]
  • Exponential Mean-Square Stability of Numerical Solutions to Stochastic Differential Equations
    London Mathematical Society ISSN 1461–1570 EXPONENTIAL MEAN-SQUARE STABILITY OF NUMERICAL SOLUTIONS TO STOCHASTIC DIFFERENTIAL EQUATIONS DESMOND J. HIGHAM, XUERONG MAO and ANDREW M. STUART Abstract Positive results are proved here about the ability of numerical simula- tions to reproduce the exponential mean-square stability of stochastic differential equations (SDEs). The first set of results applies under finite-time convergence conditions on the numerical method. Under these conditions, the exponential mean-square stability of the SDE and that of the method (for sufficiently small step sizes) are shown to be equivalent, and the corresponding second-moment Lyapunov exponent bounds can be taken to be arbitrarily close. The required finite-time convergence conditions hold for the class of stochastic theta methods on globally Lipschitz problems. It is then shown that exponential mean-square stability for non-globally Lipschitz SDEs is not inherited, in general, by numerical methods. However, for a class of SDEs that satisfy a one-sided Lipschitz condition, positive results are obtained for two implicit methods. These results highlight the fact that for long-time simulation on nonlinear SDEs, the choice of numerical method can be crucial. 1. Introduction Suppose that we are required to find out whether a stochastic differential equation (SDE) is exponentially stable in mean square. In the absence of an appropriate Lyapunov function, we may carry out careful numerical simulations using a numerical method with a ‘small’ step size t. Two key questions then follow. (Q1) If the SDE is exponentially stable in mean square, will the numerical method be exponentially stable in mean square for sufficiently small t? (Q2) If the numerical method is exponentially stable in mean square for small t, can we infer that the underlying SDE is exponentially stable in mean square? These questions deal with an asymptotic (t →∞) property, and hence they cannot be answered directly by applying traditional finite-time convergence results.
    [Show full text]
  • Exponential Stabilization of Driftless Nonlinear Control Systems" Robert Thomas M'closkey I11
    TECHNICAL MEMORANDUM NO. CIT-CDS 95-003 January 1995 "Exponential Stabilization of Driftless Nonlinear Control Systems" Robert Thomas M'Closkey I11 Control and Dynamical Systems California Institute of Technology Pasadena, CA 91125 Exponential Stabilization of Driftless Nonlinear Control Systems Robert Thomas M7Closkey I11 Department of Mechanical Engineering California Institute of Technology Abstract This dissertation lays the foundation for practical exponential stabilization of driftless control systems. Driftless systems have the form, Such systems arise when modeling mechanical systems with nonholonomic con- straints. In engineering applications it is often required to maintain the mechan- ical system around a desired configuration. This task is treated as a stabilization problem where the desired configuration is made an asymptotically stable equi- librium point. The control design is carried out on an approximate system. The approximation process yields a nilpotent set of input vector fields which, in a special coordinate system, are homogeneous with respect to a non-standard dilation. Even though the approximation can be given a coordinate-free in- terpretation, the homogeneous structure is useful to exploit: the feedbacks are required to be homogeneous functions and thus preserve the homogeneous struc- ture in the closed-loop system. The stability achieved is called p-exponential stability. The closed-loop system is stable and the equilibrium point is exponen- tially attractive. This extended notion of exponential stability is required since the feedback, and hence the closed-loop system, is not Lipschitz. However, it is shown that the convergence rate of a Lipschitz closed-loop driftless system cannot be bounded by an exponential envelope. The synthesis methods generate feedbacks which are smooth on Rn \ (0).
    [Show full text]