5th IEEE Conference on Decision Control European Control Conference (CDC-ECC) Orlo, FL, USA, December -5, A LaSalle version of Matrosov theorem Alessro Astolfi Laurent Praly Abstract A weak version of Matrosov theorem, in the spirit of LaSalle invariance principle, is established. The result is clarified by means of two examples. The case of nested- Matrosov theorem is also discussed. I. INTRODUCTION The property of stability of equilibria ( motions) of dynamical systems has been studied extensively during the last century. Several fundamental results characterizations have been derived, see e.g. [], [], [3], [4] for a classical account [5], [6] for a modern discussion. In the present paper we study a variation of Matrosov theorem with the aim to establish stability properties, or otherwise, of sets rather than of equilibria. In particular, unlike the classical Matrosov theorem which allows to prove asymptotic stability of an equilibrium provided a linear combination of some functions is positive definite, we study the case in which such a linear combination is only positive semi-definite. In this respect, our result stays at Matrosov theorem as LaSalle invariance principle stays at Lyapunov theorem for asymptotic stability. Unlike LaSalle invariance principle, however, the conclusions that can be drawn from our weak version of Matrosov theorem are not necessarily of the LaSalle type, since additional technical assumptions are required to prove convergence results. For convenience, we recall (the simplest version of) Matrosov theorem, for time-invariant nonlinear systems, as proved in [7], see also [8] for a proof based on the use of a Lyapunov function [9] for some generalizations. Theorem (Matrosov): Consider a nonlinear, timeinvariant, system described by the equation ẋ = f(x), () where x(t) R n, f : R n R n is locally Lipschitz continuous. Assume x = is an equilibrium point, i.e. f() =. Consider the following conditions. (M) There exists a differentiable, positive definite radially unbounded function V : R n R + such that V. () A. Astolfi is with the Dept. of Electrical Electronic Engineering, Imperial College London, London, SW7 AZ, UK DISP, University of Roma Tor Vergata, Via del Politecnico, 33 Rome, Italy, E-mail: a.astolfi@ic.ac.uk. This work is partly supported by the EPSRC Programme Grant Control For Energy Sustainability EP/G66477 by the MIUR under PRIN Project Advanced Methods for Feedback Control of Uncertain Nonlinear Systems. L. Praly is with MINES ParisTech, CAS, Mathématiques et Systèmes, 35, rue St-Honoré, 7735 Fontainebleau CEDEX, France, E-mail: Laurent.Praly@mines-paristech.fr. (M) (M3) There exists two differentiable functionsv : R n R V : R n R, two continuous, positive semi-definite functions N : R n R + N : R n R + a continuous function ψ : R R, satisfying ψ () =, such that V N, (3) V N +ψ (N ). (4) There exists a positive definite function ω : R n R such that the following hold. N (x)+n (x) ω(x). (5) (M), (M) (M3) imply that the equilibrium x = of system () is globally asymptotically stable. (M) (M3) imply that all bounded trajectories of system () converge to the equilibrium x =. Remark : The functions V V are called auxiliary functions. Often, Matrosov theorem is used in the simplified version in which V = V. The problem that we address in this paper is to study what happens replacing the condition (5) with the following condition. (M3) The function is positive semi-definite. N (x)+n (x) (6) One might expect, borrowing from LaSalle invariance principle, that the use of the new, weaker, condition on N +N would allow to prove convergence of the trajectories of the system to the largest invariant set contained in Ω N = {x R n N (x) +N (x) = }. This is unfortunately not true in general, but may be true under additional technical assumptions (which clarify the role of stability in the proof of such convergence result). The rest of the paper is organized as follows. In Section II we present the main results of the paper, which are then proved using a technical result given in Section III. Section IV contains two examples which highlight the significance of the new stability theorem. Finally, Section V provides some concluding remarks. II. MAIN RESULTS The first result of this section is a weak version of Theorem, whereas the second result is a weak version of the general result of [8], in which several auxiliary functions are used. 978--684-799-3//$6. IEEE 3
Theorem : Consider the nonlinear, time-invariant, system () the following conditions. (M4) The largest invariant set contained in the set (M5) Ω N,N = {x R n N (x) = N (x) = } (7) is stable. The set Ω N = {n R + x : N (x) = n,n (x) = } (8) is a singleton, i.e. N takes a single value on the set {x R n N (x) = }. the following hold. ) (M), (M) (M3) imply that all trajectories of the system () are such that inf [N (x(t))+n (x(t))] =. (9) ) (M), (M), (M3) one of (M4) (M5) imply that all trajectories of the system () are such that [N (x(t))+n (x(t))] =. () 3) (M) (M3) imply that all bounded trajectories of the system () satisfy condition (9). 4) (M), (M3) one of (M4) (M5) imply that all bounded trajectories of the system () satisfy condition (). Theorem can be generalized to the case in which several auxiliary functions are used, thus providing a weak version of the Matrosov Theorem proved in [8], as detailed in the following statement, in which only the case exploiting condition (M4) is considered. The proof of the statement, relying on Lemma, on arguments similar to those in the proof of Theorem is omitted. Theorem 3: Consider the nonlinear, time-invariant, system () the following conditions. (M) v There exists k differentiable functions V : R n R,, V k : R n R, k continuous, positive semi-definite functions N : R n R +,, N k : R n R + k continuous functions ψ : R R,, ψ k : R R, satisfying ψ i () =, such that V N, V N +ψ (N ),. V k N k +ψ k (N,N,,N k ), V k N k +ψ k (N,N,,N k ). () (M3) v The function N (x)+n (x)+ +N k (x) () is positive semi-definite. (M4) v The largest invariant set contained in the set Ω N,,N k = {x R n N (x) = = N k (x) = } (3) is stable. the following hold. ) (M), (M) v (M3) v imply that all trajectories of the system () are such that inf [N (x(t)+n (x(t))+ +N k (x(t))] =. (4) ) (M), (M) v, (M3) v (M4) v imply that all trajectories of the system () are such that [N (x(t)+n (x(t))+ +N k (x(t))] =. (5) 3) (M) v (M3) v imply that all bounded trajectories of the system () satisfy condition (4). 4) (M) v, (M3) v (M4) v imply that all bounded trajectories of the system () satisfy condition (5). III. TWO LEMMAS The following Lemmas are intended to be used for studying solutions of ordinary differential equations which are known to exist on [,+ ), taking values in a compact set. This explains why boundedness of various functions can be assumed. Lemma : Let i {,}. Let a i : R + [,+ ) be bounded differentiable functions b i : R + R + be bounded continuous functions. Suppose there exist continuous positive definite functions α i : R + R + a continuous function β : R + R +, satisfying β () =, such that ȧ α (b ), ȧ α (b )+β (b ). α (b (s))ds < + inf [b (t)+b (t)] =. (6) Remark : As will be illustrated by the examples in Section IV, the hypotheses of Lemma, which imply inf [b (t)+b (t)] =, do not guarantee, in general, the condition [b (t)+b (t)] =. The extension of Lemma to k differential inequalities is given in the following statement. Lemma : Let k 3, i {,...,k} j {,...,k}. Suppose there exist bounded differentiable functions a i : R + [,+ ), bounded continuous functions b i : R + R +, continuous positive definite functions α i : R + R + 3
continuous functions β j : R j + R +, satisfying β j () =, such that ȧ α (b ), ȧ α (b )+β (b ),. ȧ k α k (b k )+β k (b,b,,b k ), ȧ k α k (b k )+β k (b,b,,b k ). α (b (s))ds < + inf [b (t)+b (t)+ +b k (t)] =. IV. EXAMPLES In this section we present two examples which allow to clarify the results of the paper to highlight that the condition () cannot be improved without additional assumption. Consider the 3-dimensional system ẋ = (x + +xp 3 )x, ẋ = (x + +xp 3 )x, (7) ẋ 3 = x q 3, where x + = max{x,}, p is a positive even integer q is a positive odd integer. The set of equilibrium points is given by {(x,x,x 3 ) : x + = x 3 = }. A. The case p =, q = 3 Let V (x,x,x 3 ) = V (x,x,x 3 ) = x +x +x 3 V = V = x 4 3 V (x,x,x 3 ) = x. V = (x + +x 3)x = ( x + +x3) x+ + V x, (8) wherex = min{x,}. As a result, conditions (M), (M) (M3) hold hence x 4 ( 3(τ)dτ < +, inf x 3 (t)+x + (t) 3) =. (9) We now study if it is possible to obtain stronger asymptotic properties. For, observe that the first equality in (8) implies that the origin is globally stable that all the solutions converge to the largest invariant set contained in the intersection of a level set of V with the set {(x,x,x 3 ) : x 3 = }. In such a set, we have V = x 3 +, 3 one could be tempted to claim that the solutions converge to the largest invariant set contained also in the set {(x,x,x 3 ) : x 3 = x + = }. This is not the case, since, in general, a generic solution does not inherit the properties of the solutions in the invariant set. To clarify this statement re-write the system using polar coordinates (θ,ρ) in the (x,x )-plane, i.e. note that ρ(t) = ρ(), As a result ρ =, θ = (ρ cos(θ) + +x 3 ), ẋ 3 = x 3 3, θ(t) θ() ( x 3 (t) = exp x 3(s)ds x 3 (s)ds )x 3 (). θ(t) [θ() log(x 3 ())] + log(x 3 (t)), since x 3 (t) =, θ(t) tends to modulo π, i.e. θ(t) does not converge. In addition, since ρ(t) is constant, the vector (x (t),x (t)) has a constant modulus does not stop turning around the origin, which implies that but also that x 3(t) =, inf x +(t) =, supx + (t) = x (). This last equation shows that the asymptotic property expressed by the second of equations (9) cannot be improved. Figure shows the phase portrait of the trajectories of the system with initial condition x() = [ ], whereas Figure shows the time histories of the states x, x x 3. Note that the time axis is in log-scale. Figure highlights that all trajectories with initial condition off the (x,x )- plane have an oscillatory behavior with a period that tends to infinity. Note that trajectories with initial conditions such that x 3 () = converge to the set {(x,x ) x +x = x () +x (),x }, i.e. to a semi-circle centered at the origin, the size of which depends upon the initial conditions. This set is not stable, hence condition (M4) does not hold. Remark 3: Theω-it set of the trajectories of the system starting off the (x,x )-plane is, as detailed in [], a chain recurrent set, which strictly contains the ω-it set of the trajectories of the system starting in the (x,x )- plane, consistently with the results in [] on asymptotically autonomous semiflows. Remark 4: As a consequence of the discussion in this section, the (x,x )-subsystem of system (7), with p = q = 3, x 3 regarded as an input, does not possess the converging-input converging-state property, see []. This
.8.6.8.6 4.4.4 6.. x x logx 3 8...4.4.6.8.6.8 4 4 6 8 log 4 6 8 log 6 4 6 8 log Fig.. Time histories of the states of the system (7) with p =, q = 3 x() = [ ]. Note that the states x x undergo fast transients. These perturbations are such that x 3..8.6.4..5 x.5.5 x Fig.. The trajectory of the system (7), with p =, q = 3 x() = [ ]..5 u (t) +u (t) = x 3 (t) 4 (x () +x () ) x 3 () 4 (x () +x () ), [u (s) +u (s) ]ds = (x () +x () ) x 3 (s) 4 ds (x () +x () ) x 3(), hence, selecting x 3 () = o( ε), the constraints on (u,u ) are satisfied. The consequence is that the ω-it set of any solution of the (x,x )-subsystem with x 3 =, which is composed of equilibrium points located on the portion of a circle centered at the origin having the x -coordinate non-positive, can be transformed into the full circle by such small perturbations. does not contradict the result in [], which highlights (among other things, similarly to what done in this paper) the importance of asymptotic stability (of an equilibrium, or of a set) to establish asymptotic properties of solutions. Remark 5: The result of the above discussion has to be interpreted on the basis of the results in []. Therein it is shown that the ω-it set of any solution of the sub-system ẋ = x +x, ẋ = x +x is chain recurrent. This says approximately that any point of this set can be made a point of an homoclinic orbit when solutions are taken as it, for ε, of solutions, on [,+ ), of the perturbed system ẋ = x + x +u (t), ẋ = x + x + u (t), where (u,u ) has L L -norm smaller than ε. For system (7) the perturbations are given by u (t) = x 3 (t) x (t), u (t) = x 3 (t) x (t). B. The cases p =, q = p = 4, q = 3 Consider system (7) with p =, q = or p = 4, q = 3. In this case we have V = (x + +x 3)x = ( x + +x p 3) x+ V x, yielding V +k V ( x + +x p 3) x+, () for all (x,x,x 3 ) such that x > k. Recall now that all trajectories are bounded, i.e. remain in a compact set. Hence it is possible to select k (as a function of the initial conditions) such that the inequality () holds for all t. By LaSalle invariance principle we conclude that all solutions converge to the largest invariant set contained in a level set of V in the set {(x,x,x 3 ) : x + = x 3 = }. This convergence result is illustrated in Figure 3, which shows a trajectory with initial state x() = [ ]. Similarly to the case in Section IV-A, the x x states show an oscillatory behavior but, unlike the case in Section IV-A, x (t) x (t) converge to a point such that x + =. 33
.8.6.4. x x x 3..4.6.8 4 6 8.8.6.4...4.6.8 4 6 8.8.6.4. 4 6 8 Fig. 3. Time histories of the states of the system (7) with p =, q = x() = [ ]. Note that the states x x undergo fast transients. V. CONCLUSIONS A LaSalle version of Matrosov theorem has been stated proved. It has been shown that stability plays a crucial role in the study of the asymptotic behavior of trajectories that a naive application of ideas borrowed from LaSalle invariance principle may yield wrong conclusions. The theory has been illustrated by means of two simple examples. REFERENCES [] J. La Salle S. Lefschetz, Stability by Lyapunov s Direct Method. Academic Press, New York, 96. [] W. Hahn, Stability of Motion. Springer Verlag, New York, 967. [3] J. La Salle, The Stability of Dynamical Systems. SIAM, Philadelphia, 976, with Appendix: Limiting Equations Stability of Nonautonomous Ordinary Differential Equations, By Z. Artstein. [4] A. Lyapunov, The General Problem of the Stability of Motion. Taylor & Francis, London, 99. [5] A. Bacciotti, Local Stabilizability of Nonlinear Control Systems. Series on Advances in Mathematics for Applied Sciences, World Scientific, 99. [6] H. Khalil, Nonlinear Systems, 3rd ed. Prentice-Hall, Upper Saddle River,. [7] V. Matrosov, On the stability of motion, J. Appl. Math. Mech., vol. 6, pp. 337 353, 96. [8] F. Mazenc D. Nesic, Lyapunov functions for time-varying systems satisfying generalized conditions of Matrosov theorem, Math. Control Signals Systems, vol. 9, pp. 5 8, 7. [9] A. Loria, E. Panteley, D. Popovic, A. Teel, A nested Matrosov theorem persistency of excitation for uniform convergence in stable non-autonomous systems, IEEE Trans. Autom. Control, vol. 5, pp. 83 98, 5. [] K. Mischaikow, H. Smith, H. Thieme, Asymptotically autonomous semiflows: chain recurrence lyapunov functions, Trans. Amer. Math. Soc., vol. 347, pp. 669 685, 995. [] E. Sontag, A remark on the converging-input converging-state property, IEEE Trans. Autom. Control, vol. 48, no., pp. 33 34, 3. [] C. Conley, Isolated invariant sets the Morse index, ser. CBMS Regional Conference series in Mathematics. American Mathematical Society, 978, no. 38. 34