arxiv: v1 [math.ds] 17 Mar 2017

Transcription

1 To appear in Dynamics of Continuous, Discrete and Impulsive Systems journal Dwell time for switched systems with multiple equilibria on a finite time-interval Oleg Makarenkov and Anthony Phung arxiv: v [math.ds] 7 Mar 207 Department of Mathematical Sciences University of Texas at Dallas, Richardson, TX 75080, USA Abstract. We describe the behavior of solutions of switched systems with multiple globally exponentially stable equilibria. We introduce an ideal attractor and show that the solutions of the switched system stay in any given -inflation of the ideal attractor if the frequency of switchings is slower than a suitable dwell time T. In addition, we give conditions to ensure that the -inflation is a global attractor. Finally, we investigate the effect of the increase of the number of switchings on the total time that the solutions need to go from one region to another. Keywords. Switched system, dwell-time, global exponential stability, ideal attractor. AMS (MOS) subject classification: 93C30; 34D23 Introduction Dwell time is the lower bound on the time between successive switchings of the switched system ẋ = f u(t) (x), u(t) is a piecewise constant function, x R n, () which ensures a required dynamic behavior under the assumption that each of the subsystems ẋ = f u (x), u R, x R n, (2) possess a globally stable equilibrium x u. When all the equilibria {x u(t) } t t0 coincide, the dwell time T > 0 which gives global exponential stability of the common equilibrium x 0 is computed e.g. in Liberzon [4, 3.2.]. Specifically, the result of [4, 3.2.] gives a formula for T which makes x 0 globally exponentially stable for any piecewise constant function u(t) whose discontinuities t, t 2,... verify t i t i T. (3) The case where the equilibria are distinct is covered in Alpcan-Basar [], who offered a dwell time T that ensures global exponential stability of a

2 2 Oleg Makarenkov and Anthony Phung suitable set A {x u(t) } t t0 for any u(t) whose discontinuities verify (3). The problem of stability of switched systems with multiple equilibria appears e.g. in differential games, load balancing, agreement and robotic navigation (see [, 5] and references therein). A deeper analysis of the dynamics of switched systems with multiple equilibria was recently carried out in Xu et al [6], who gave a sharp formula for the attractor A in the case of quasi-linear switched systems (). Assuming that u(t) is periodic and denoting by t X u (t, x) the solution of (2) with the initial condition X u (0, x) = x, the paper [6] investigated the asymptotic attractivity of A = t t 0, τ t 0 {X u(τ) (t, x u(τ) )}. The motivation for our paper comes from the problem of planning the motion of a 3-D walking robot, where turn left, walk straight and turn right correspond to u(t) =, u(t) = 0 and u(t) = respectively, see Gregg et al [3]. It is not the asymptotic attractivity of A which is of importance for the robot turning maneuver but rather an appropriate attractivity of A during the time of the maneuver. The goal of this paper is to provide a dwell time which can ensure the required attractivity. The paper is organized as follows. In the next section of the paper, we prove our main result (Theorem 2.). Given > 0, Theorem 2. provides a dwell time T > 0 such that the solutions of () with the initial conditions in the -neighborhood B (A) of A never leave B (A) in the future. Theorem 2. can be viewed as a version of [6, Theorem ] for fully nonlinear systems. In section 3, we compute (Theorem 3.) a dwell time to ensure that the attractor B (A) is reached asymptotically from any initial condition. The proof of Theorem 3. follows the ideas of Alpcan-Basar []. However, we offer weaker conditions where the Lyapunov functions of subsystems (2) are not supposed to respect any uniform estimates. A particular case study where the Lyapunov functions of subsystems (2) are shifts of one another is addressed in section 4. In this section, we consider a switched system which switches between two subsystems u = u and u = u 2 and analyze the solutions of the switched system with the initial conditions in B (A). Let x and x 2 be the equilibria of subsystems u = u and u = u 2 respectively. The result of section 4 (Theorem 4.) clarifies whether or not the solutions from the neighborhood of x reach the neighborhood of x 2 faster if the switching signal is amended in such a way that an additional switching occurs between u = u and u = u 2. In other words, section 4 investigates whether or not adding more discrete events is alone capable of making the dynamics inside B (A) faster. Examples 2. and 4. illustrate the conclusions of Theorems 2. and 4..

3 Dwell time for switched systems on a finite time-interval 3 2 The local trapping region Let x u be the unique equilibrium of (2). We assume that for any u, system (2) admits a global Lyapunov function V u such that α u ( x x u ) V u (x) β u ( x x u ), x R n, (4) (V u ) (x)f u (x) k u V u (x), x R n, (5) where α, β are strictly monotonically increasing functions with α u (0) = β u (0), and k u > 0. Introduce the following trapping regions Nu = {x : V u (x) }, L u,u 2 (t) = {X u2 (t, x)} x Nu (6) and define the dwell time that the solutions need to go from N u to N u 2 as Theorem 2.. Assume that Tu,u 2 = ln ( k u2 β u2 xu2 x u + αu ). (7) () (A) f u C (R n, R n ) for any u R, (A2) for any u R, system (2) admits an equilibrium x u whose Lyapunov function V u satisfies (4)-(5), where α, β C 0 (R, R) are strictly increasing functions, α(0) = β(0) = 0, k u > 0, (A3) u : [t 0, ) R is a piecewise constant function. Let {t, t 2,...} = {t i } i I be a finite or infinite increasing sequence of points of discontinuity of u and u i = u(t i + 0). If t i t i T u i,u i, then for any solution x of () with x(t i ) Nu i, i I, one has x(t) L u i,u i (t), t i t t i, i I, (8) x(t i ) Nu i, i I. (9)

4 4 Oleg Makarenkov and Anthony Phung x ui {x:v ui (x)= {x:v ui (x)= ( ) u i x ui ( ) i u {x:v ui (x)= Figure : Illustration of the proof of Theorem 2.. The two gray rings are the estimates for the sets {x : V ui (x) = } and {x : V ui (x) = δ} given by condition (4). Proof. We only have to prove (9) because the validity of (8) follows directly from the definition of L u i,u i. Let x(t i ) Nu i. Our goal is to show that x(t i ) Nu i. Given > 0, define δ > 0 as ( ) δ = β ui x ui x ui + αu i (). By construction (see Fig. ), Nu δ i Nu i, and so x(t i ) Nu δ i. Introduce By (5) we have v(t) = V ui (x(t)). v(t) k i v(t), t i t t i, v(t i ) δ. By the comparison lemma (see e.g. [2, Lemma 6.4]), it holds that where p(t) is the solution of At the same time, v(t) p(t), ṗ(t) = k i p(t), t i t t i, p(t i ) = δ. p(t i ) = e ku i (ti ti ) δ e ku i T u i,u i δ = ( β ui xui x ui + αu )δ =. i ()

5 Dwell time for switched systems on a finite time-interval 5 Therefore, V ui (x(t i )), which completes the proof. Theorem 2. suggests the following definition of the -inflation A of the ideal attractor A of (). Given a function u : [t 0, ) R and the respective increasing sequence (t, t 2,...) = {t i } i I, let u i = u(t i ) and L u A (t) = i,u i (t), t i t < t i, i I, L u max(i),u max(i) (t), t t max(i), if I is finite. Corollary 2.. Let the assumptions (A)-(A3) of Theorem 2. hold. If t i t i sup Tu i,u i =: Tloc, i I, i I then, for any solution x of () with the initial condition x(t 0 ) A (t 0 ), one has x(t) A (t), t t 0. Note, sup Tu i,u i i I [t 0, ). is finite when t u(t) takes a finite number of values on Example 2.. To illustrate Theorem 2., we consider the following switched system (slightly modified from Example 2 in []) ( ) ( ) u ẋ = x +, (0) whose unique equilibrium is given by x u = ( u 2 u + ). Introduce the three discrete states u, u 2, and u 3 as and consider u =, u 2 = 0, u 3 =, = 5. If the Lyapunov function V u (x) is selected as then formulas (6) and (7) yield V u (x) = x x u 2, N u = { x : x x u }, T loc.426.

6 6 Oleg Makarenkov and Anthony Phung Therefore, for the control input u, t [0, T ), u2, t [T, 2T ), u(t) = u3, t 2T, T =.43, () and for any solution x of (0), Theorem 2. ensures the following: if x(0) Nu, then x(t ) Nu2 and x(2t ) Nu x2 x2 Figure 2(left) documents the sharpness of the dwell time T. Indeed, the figure shows that if the initial condition x(0) deviates to the outside of Nu just a little bit, then the dwell time T is no longer sufficient to get x(t ) Nu2 (though we still have x(2t ) Nu3 for this solution) x x Figure 2: Left: Solutions of switched system () with initial conditions (blue dots) inside Nu5, on the boundary of Nu5, and outside Nu5, and for the control input u(t) given by (). Right: The solution of switched system () with the initial condition (, )T for the 4T -periodic control input (2). To demonstrate that trapping regions Nu, Nu2, Nu3 (and thus the -inflated attractor A, see Corollary 2.) provide a rather sharp estimate for the location of the attractor of (0), we extend the input u(t) to [0, 4T ] as u, t [0, T ), u2, t [T, 2T ), u(t) = (2) u3, t [2T, 3T ), u2, t [3T, 4T ), and then continue it to the entire [0, ) by 4T -periodicity. The respective solution x of (0) with the initial condition x(0) = (, )T is plotted in Fig. 2 (right). The drawing shows that the switching points of the solution x are very close to the boundaries of the trapping regions Nu, Nu2, Nu3, i.e. there is only a little window to reduce the size of those regions.

7 Dwell time for switched systems on a finite time-interval 7 3 Global attractivity of the local trapping region Theorem 3.. Let the assumptions (A)-(A3) of Theorem 2. hold and I be infinite. Fix > 0 and suppose that there exists constants µ i () such that V ui+ (x) V ui (x) Finally, assume that µ 0 ()... µ i ()e µ i (), x R n \ N u i, i N {0}. ti+ t 0 k u(s) ds 0, as i. Then, x( ˆT ) N u i for some ˆT > 0 and some i N. Proof. Let W (t) = e ku(t)t V u(t) (x(t)), where t [t 0, ). Then, for t [t i, t i+ ), W (t) = k ui W (t) + e ku i t d dt V u i (x(t)) k ui W (t) k ui e ku i t V ui (x(t)) = 0, which means that W is decreasing on [t i, t i+ ). In particular, On the other hand, Therefore, W (t + i ) W (t i+ ). W (t + i+ ) W (t i+ ) = ekui+ ti+ Vu i+(x(t i+ )) eku i+ ti+ µ i () =: µ i. e ku i ti+ V ui (x(t i+ )) e ku i ti+ µ i W (t + i ) µ iw (t i+ ) W (t+ i+ ) (3) Replacing i by i and combining with (3), one gets µ i µ i W (t + i ) W (t + i+ ). Continuing this process for i 2, i 3, etc. we obtain µ 0... µ i W (t + 0 ) W (t+ i+ ). Applying e ku i+ ti+ yields or, equivalently, µ 0... µ i e ku i+ ti+ W (t + 0 ) V u i+ (x(t i+ )), V u0 (x(t 0 ))e i j=0 ku j (tj+ tj) µ 0 ()... µ i () V ui+ (x(t i+ )). The left-hand-side approaches 0 as i by the assumption of the theorem. Therefore, V ui+ (x(t i+ )) 0 as i. The proof is complete.

8 8 Oleg Makarenkov and Anthony Phung Corollary 3.. (Alpcan-Basar []) Let the assumptions (A)-(A3) of Theorem 2. hold and I be infinite. Fix > 0 and suppose that there exists a constant µ() > such that V u(t) (x) V u(τ) (x) µ(), x Rn \N u(τ), τ, t t 0. Finally, assume that k = inf t t 0 k u(t) > 0 and consider T glob satisfying If Tglob > ln(µ()). k t i t i T glob, i N, then, x( ˆT ) N u i for some ˆT > 0 and some i N. Proof. Let γ > 0 be such that Tglob as given by Theorem 3.. Then e i ln(µ() + γ) =. Let µ 0 (),... µ i () be k j=0 ku j (tj+ tj) µ 0 ()... µ i () e kit glob µ() i = ( ) i = (µ() + γ) i µ() µ() i = 0 as i. µ() + γ Corollary 3.2. Let the conditions of Corollary 3. hold. Let T loc and T glob be those given by Corollaries 2. and 3.. If t i t i max { T loc, T glob}, i N, then, for any solution x of (), there exists ˆT > t 0 such that x(t) A (t), t ˆT. 4 Dependence of the dwell time on the number of discrete states Suppose that u(t) switches from u 0 to u at t = t 0. According to Theorem 2., it takes at most time Tu 0,u (see formula (7)) for a trajectory x of () to go from Nu 0 to Nu. The next theorem shows that adding more discrete states between u 0 and u makes the travel time from Nu 0 to Nu longer. Theorem 4.. Let the assumptions (A)-(A2) of Theorem 2. hold and suppose α u =: α, β u =: β, k u =: k don t depend on u. Fix d > 0 and r > 0. Then there exists 0 > 0 such that for any T u 0,u < T u 0,v + T v,u, (0, 0 ), x u0 d, x u d, x u0 x v r, x v x u r. (4)

9 Dwell time for switched systems on a finite time-interval 9 Proof. By formula (7) one has Tu0,u Tu0,v Tv,u + = ln k β (kxu0 xu k + α ()) + ln + k β (kxu0 xv k + α ()) + ln = k β (kxv xu k + α ()) K = ln /k, (5) where K=!/k /k β kxv xu k + α (). β kxu0 xv k + α () β (kxu0 xu k + α ()) x2 x2 Observe that there exists K0 > 0 such that K K0 for any functions u0, u, v that verify (4) as long as d > 0 and r > 0 stay fixed. Therefore, it is possible to choose 0 > 0 (which depends on just d > 0 and r > 0) to satisfy K//k > for all (0, 0 ). The proof is complete x x Figure 3: Solutions of switched system () with the initial condition in Nu5 for the control inputs u e(t) (Left) and u(t) (Right) of Example 4.. Example 4.. In order to illustrate Theorem 4., we refer to Example 2. again. Figure 3 shows the graphs of the solutions x of (0) for two control inputs u, t [0, Tu,u2 ), u, t [0, Tu,u3 ), u2, t [Tu,u2, Tu,u2 + Tu2,u3 ), u e(t) = u(t) = u3, t Tu,u3, u3, t Tu,u2 + Tu2,u3, over the time interval [0, Tu,u3 ]. The plotting documents that Tu,u2 +Tu2,u3 turns out to be a longer time compared to Tu,u3.

10 0 Oleg Makarenkov and Anthony Phung 5 Conclusion In this paper we considered a switched system of differential equations under the assumption that the time between two successive switchings is greater than a certain number T called dwell time. We proved (Theorem 2.) that a suitable choice of the dwell time makes the solution stay within a required neighborhood A of a so-called ideal attractor. We further proved that the solutions reach A asymptotically if the initial conditions don t belong to A. By doing that we obtained a new integral condition (Theorem 3.) for global stability which didn t seem to appear in the literature before. Finally, we addressed a case study where the Lyapunov functions of different subsystems are just shifts of one another. Here we used the dwell time formulas from Theorem 2. to estimate the time that the trajectories need to go from the neighborhood of an equilibrium of one subsystem to the neighborhood of an equilibrium of another subsystem (i.e. we considered a switched system with two discrete states). We proved (Theorem 4.) that adding more discrete states makes this travel time longer. Examples 2. and 4. show that our theoretical conclusions agree with numeric simulations. 6 Acknowledgements The first author is partially supported by NSF Grant CMMI References [] T. Alpcan, T. Basar, A stability result for switched systems with multiple equilibria, Dyn. Contin. Discrete Impuls. Syst. Ser. A Math. Anal. 7, (200), no. 6, [2] Amann, Herbert Ordinary differential equations. An introduction to nonlinear analysis. Translated from the German by Gerhard Metzen. de Gruyter Studies in Mathematics, 3. Walter de Gruyter & Co., Berlin, 990. [3] R. D. Gregg, A. K. Tilton, S. Candido, T. Bretl, M. W. Spong, Control and Planning of 3-D Dynamic Walking With Asymptotically Stable Gait Primitives, IEEE Transactions on Robotics, 28 (202), issue 6, [4] D. Liberzon, Switching in systems and control. Systems & Control: Foundations & Applications. Birkhauser Boston, Inc., Boston, MA, [5] S. Mastellone, D. M. Stipanovic, M. W. Spong, Stability and Convergence for Systems with Switching Equilibria, Proceedings of the 46th IEEE Conference on Decision and Control (2007), [6] H. Xu, Y. Zhang, J. Yang, G. Zhou, L. Caccetta, Practical exponential set stabilization for switched nonlinear systems with multiple subsystem equilibria, J. Global Optim. 65 (206), no., journal@monotone.uwaterloo.ca journal/