arxiv:6.58v [cs.sy] 7 Apr 7 Observer design for piecewise smooth and switched systems via contraction theory Davide Fiore Marco Coraggio Mario di Bernardo, Department of Electrical Engineering and Information Technology, University of Naples Federico II, Via Claudio, 85 Naples, Italy Department of Engineering Mathematics, University of Bristol, BS8 TR Bristol, U.K. Abstract: The aim of this paper is to present the application of an approach to study contraction theory recently developed for piecewise smooth and switched systems. The approach that can be used to analyze incremental stability properties of so-called Filippov systems (or variable structure systems) is based on the use of regularization, a procedure to make the vector field of interest differentiable before analyzing its properties. We show that by using this extension of contraction theory to nondifferentiable vector fields, it is possible to design observers for a large class of piecewise smooth systems using not only Euclidean norms, as also done in previous literature, but also non-euclidean norms. This allows greater flexibility in the design and encompasses the case of both piecewise-linear and piecewise-smooth (nonlinear) systems. The theoretical methodology is illustrated via a set of representative examples. Keywords: contraction theory, observer design, incremental stability, discontinuous control, regularization. INTRODUCTION The problem of designing state observers for nondifferentiable systems is the subject of current research. For example, the design of observers for Lipschitz continuous nonlinear systems was investigated in (Rajamani, 998; Zemouche and Boutayeb, 3), while in (Arcak and Kokotovic, ; Brogliato and Heemels, 9) design approaches based on passivity theory were proposed for Lur e-type systems. Also, in (Juloski et al., 7; Doris et al., 8) sufficient conditions were presented to ensure stability of the estimation error for state observers of bimodal piecewise linear (PWL) systems (both continuous and discontinuous on the switching surface). The analysis was conducted analyzing the quadri-modal estimation error dynamics based on quadratic Lyapunov functions and LMIs. Related results were presented in (van de Wouw and Pavlov, 8) for the case of piecewise affine (PWA) systems. Therein, using theoretical results developed in (Pavlov et al., 7), sufficient conditions guaranteeing exponential stability of the estimation error were given in terms of a set of appropriate LMIs. More recently, the state estimation problem was investigated in (Heemels et al., ) for linear complementarity systems and in (Forni et al., 3) for hybrid systems with impacts. Contraction theory (Lohmiller and Slotine, 998; Russo et al., ; Jouffroy, 5; Forni and Sepulchre, 4; Aminzare and Sontag, 4) is a powerful analysis tool providing sufficient conditions for incremental stability (Angeli, ) of a dynamical system. Namely, if the system vector field is contracting in a set of interest, any two of its Corresponding author Mario di Bernardo. E-mail address: mario.dibernardo@unina.it; tel. +39-8-768399; fax +39-8- 768386. trajectories will converge towards each other in that set, a property that can be effectively exploited to design state observers and solve tracking control problems as discussed, for instance, in (Lohmiller and Slotine, 998; van de Wouw and Pavlov, 8; Bonnabel et al., ; Dinh et al., 3; Manchester and Slotine, 4b,a; di Bernardo and Fiore, 6). More specifically, incremental exponential stability over a given forward invariant set is guaranteed if some matrix measure, say, of the system Jacobian matrix is uniformly negative in that set for all time. The original results on contraction analysis were presented for continuously differentiable vector fields limiting their application to observer design for this class of dynamical systems. Recently, extensions have been presented in the literature for applying contraction and convergence analysis to different classes of nondifferentiable and discontinuous vector fields (Lohmiller and Slotine, ; Pavlov et al., 7; di Bernardo et al., 4; Lu and di Bernardo, 6; di Bernardo and Liuzza, 3; di Bernardo and Fiore, 4; Fiore et al., 6). In this paper we propose a methodology to design state observers for nondifferentiable bimodal vector fields, which stems from the results presented in (Fiore et al., 6) on extending contraction analysis to Filippov systems. Specifically, we derive conditions on the observer dynamics for the estimation error to converge exponentially to zero. These conditions, when particularized to the case of PWA systems, generalize those presented in (van de Wouw and Pavlov, 8) to the case of non-euclidean norms. In what follows, after reviewing some key results on contraction analysis of switched systems, we present our procedure for state observer design complementing the theoretical derivations with some illustrative examples.
. CONTRACTION ANALYSIS OF SWITCHED SYSTEMS. Incremental Stability and Contraction Theory Let U R n be an open set. Consider the system of ordinary differential equations ẋ = f(t, x), () where f is a continuously differentiable vector field defined for t [, ) and x U, that is f C (R + U, R n ). We denote by ψ(t, t, x ) the value of the solution x(t) at time t of the differential equation () with initial value x(t ) = x. We say that a set C R n is forward invariant for system () if x C implies ψ(t, t, x ) C for all t t. A nonlinear dynamical system () is contracting if it forgets initial conditions or temporary state perturbations exponentially fast, implying convergence of system trajectories towards each other and consequently towards a steady-state solution which is determined only by the input (entrainment property). The following theorem summarize the basic results of contraction theory (Russo et al., ; Lohmiller and Slotine, 998). Theorem. Let C U be a forward invariant K-reachable set. The continuously differentiable vector field () is said to be contracting on C if there exists some norm in C, with associated matrix measure (see Appendix A), such that, for some constant c > (the contraction rate), ( ) f (t, x) c, x C, t t. () x Then, for every two solutions x(t) = ψ(t, t, x ) and y(t) = ψ(t, t, y ) with initial conditions x, y C we have that x(t) y(t) K e c(t t) x y, t t, (3) that is, system () is incrementally exponentially stable (IES) in C. In this paper we analyze contraction properties of dynamical systems based on norms and matrix measures (Lohmiller and Slotine, 998; Russo et al., ). Other more general definitions exist in the literature, for example results based on Riemannian metrics (Lohmiller and Slotine, 998) and Finsler-Lyapunov functions (Forni and Sepulchre, 4). The relationships between these three definitions and the definition of convergence (Pavlov et al., 4) were investigated in (Forni and Sepulchre, 4).. Switched systems Switched (or bimodal) Filippov systems are dynamical systems ẋ = f(x) where f(x) is a piecewise continuous vector field having a codimension-one submanifold Σ as its discontinuity set (Filippov, 988; Utkin, 99). The submanifold Σ is called the switching manifold and is defined as the zero set of a smooth function h : U R n R, that is Σ := {x U : h(x) = }, where R is a regular value of h, i.e. h(x), x Σ. Σ divides U in two disjoint regions, S + := {x U : h(x) > } and S := {x U : h(x) < } (see Fig. ). Hence, a bimodal Filippov system can be defined as { f + (x), if x S + ẋ = f (x), if x S, (4) Fig.. Regions of state space: switching manifold Σ, S +, S (hatched zone) and S ε (grey zone). where f +, f C (U, R n ). We assume that solutions of system (4) are defined in the sense of Filippov (and therefore admitting sliding motions on Σ) and they have the property of right-uniqueness in U (Filippov, 988, pag. 6). Condition () was previously presented as a sufficient condition for a dynamical system to be incrementally exponentially stable, but it cannot be directly applied to system (4) because its vector field is not continuously differentiable. In recent work reported in (Fiore et al., 6), sufficient conditions were derived for convergence of any two trajectories of a Filippov system towards each other. Instead of directly analyzing the Filippov vector field on Σ, the analysis is conducted on its regularization, say f ε (x), defined as f ε (x) = + ϕ ε (h(x)) f + (x) + ϕ ε (h(x)) f (x), where ϕ ε C (R, R) is the so-called transition function. In this new system the switching manifold Σ is replaced by a boundary layer S ε (Fig. ) of width ε, defined as S ε := {x U : ε < h(x) < ε}, and, more importantly, f ε is continuously differentiable in U, so that condition () can be applied to it. Finally, contraction properties of Filippov systems (4) are recovered taking the limit for ε and considering the following Lemma. Lemma. Denoting by x ε (t) a solution to the regularized system and by x(t) a solution to the switched system with the same initial conditions x, then x ε (t) x(t) = O(ε), uniformly for all t t and for all x U. For further details see (Sotomayor and Teixeira, 996; Utkin, 99; Fiore et al., 6). The resulting sufficient conditions for a bimodal Filippov system to be incrementally exponentially stable in a certain set are stated in the following theorem (see (Fiore et al., 6) for a complete proof and further details). Theorem 3. The bimodal switched system (4) is incrementally exponentially stable in a K-reachable set C U with convergence rate c := min {c, c } if there exists some norm in C, with associated matrix measure, such that, for some positive constants c, c, ( ) f + x (x) c, x S +, ( ) f x (x) c, x S, ] ) f + (x) f (x) h(x) =, x Σ. In the above relations S + and S represent the closures of the sets S + and S, respectively.
3. STATE OBSERVER DESIGN 3. Problem formulation Consider the bimodal switched system ẋ = { f + (x) + u(t), h(x) > f (x) + u(t), h(x) <, (5) y = g(x), (6) where x R n, y R p, u R n are the state, output and the input of the system, respectively, and f +, f, g are continuously differentiable vector fields. As an observer for the system (5)-(6), we propose a bimodal Luenberger-like switched observer of the form ˆx = { f + (ˆx) + L + (y ŷ) + u(t), h(ˆx) > f (ˆx) + L (y ŷ) + u(t), h(ˆx) <, (7) ŷ = g(ˆx), (8) where ˆx(t) R n is the estimated state and L +, L R n p are observer gain matrices to be selected appropriately. We are interested in deriving conditions on the observer gain matrices L + and L that guarantee exponential convergence to of the estimation error e(t) := x(t) ˆx(t) for all x(t) : R + R n satisfying (5)-(6) for any given continuous function u(t) : R + R n. Note that in what follows we will not require system (5)-(6) to be contracting, i.e. Theorem 3 must not necessarily hold for this system. Instead, contraction theory will be exploited to analyze convergence of the system describing the dynamics of the estimation error. 3. Main results Theorem 4. The state estimation error e(t) converges exponentially to zero, that is, for some c >, e(t) K e c(t t) x(t ), t t, (9) if there exists some matrix measure, such that, for some positive constants c, c, ( ) f + g (ˆx) L+ x x (ˆx) c, ˆx : h(ˆx) >, () ( ) f g (ˆx) L x x (ˆx) c, ˆx : h(ˆx) <, () ] ) f(ˆx) + L(y ŷ) h(ˆx) =, ˆx : h(ˆx) =, () where f(ˆx) = f + (ˆx) f (ˆx) and L = L + L. Moreover, the convergence rate c can be estimated as min{c, c }. Proof. Conditions ()-() come from the application of Theorem 3 to the dynamics of the state observer (7)-(8) by rewriting them as ˆx = { f + (ˆx) + η + (t), h(ˆx) > f (ˆx) + η (t), h(ˆx) <, where f ± (ˆx) = f ± (ˆx) L ± g(ˆx) depends only on ˆx, and η ± (t) = L ± g(x(t)) + u(t) is a function of t. Hence, if such conditions are satisfied, then the state observer is contracting; this in turn implies that, for two generic solutions ˆx (t) and ˆx (t), (3) holds, i.e. ˆx (t) ˆx (t) K e c(t t) ˆx (t ) ˆx (t ), t t. Now, notice that a solution x(t) of system (5) is a particular solution of the observer (7) because (5) and (7) have the same structure, except for the correction term g(x) g(ˆx), which is null when considering x(t) as a solution of the observer. Then, we can replace ˆx (t) with x(t), rename ˆx (t) as the general solution ˆx(t), and write e(t) = x(t) ˆx(t) K e c(t t) x(t ), for all t t, where ˆx(t ) = as usual in observer design. Hence, the exponential convergence to zero of the estimation error is proved. Remark. Alternatively, the theorem can be proved considering the regularized dynamics of both system (5) and observer (7). Denoting by x ε (t) a solution to the regularized switched system (5), and by ˆx ε (t) a solution to the regularized observer (7), we have e(t) = x(t) ˆx(t) x(t) x ε (t) + x ε (t) ˆx ε (t) + ˆx ε (t) ˆx(t). The first and the third terms are the error between a solution to the discontinuous system and a solution to its regularized counterpart; hence, from Lemma we know that x(t) x ε (t) = O(ε), and ˆx(t) ˆx ε (t) = O(ε). Furthermore, similarly to what done in (Fiore et al., 6), it can be shown that conditions ()-() imply incremental stability of the trajectories of the regularized observer, thus ˆx ε, (t) ˆx ε, (t) K e c(t t) ˆx ε, (t ) ˆx ε, (t ), t t. The theorem is finally proved by taking the limit for ε + and taking the same last step as that in the proof of Theorem 4. Remark. If one of the two modes, f + or f, of the observed system (5) is already contracting, the corresponding observer gain matrix, L + or L, in (7) can be set to zero to simplify the design problem. The drawback is that the convergence rate of the estimation error will depend on that of the contracting mode that cannot be altered if this choice is made. Remark 3. In the presence of bounded disturbances or uncertainties on the models, contraction properties of the vector fields guarantee boundedness of the estimation error (a more detailed analysis is not the aim of the current paper; the interested reader can refer to (Lohmiller and Slotine, 998)). 4. EXAMPLES Here we present some examples to illustrate the use of Theorem 4 for the design of observers for switched systems. All simulations presented in this section have been computed using the numerical solver in (Piiroinen and Kuznetsov, 8). Example Consider a nonlinear bimodal switched system as in (5)-(6) with [ ] [ ] f + 9x 3x (x)= 8, f 4x 9x + 3x (x)= + 8, 4x
a) 4 b) - -4.5.5.5 6 4.5 3.5.5.5.5 Fig.. Panel a: Time evolution of the states x (t) (solid line) and ˆx (t) (dashed line) of Example, with initial conditions x = [3 3] T, ˆx = [ ] T. Panel b: Norm of the estimation error e(t). The dashed line represents the analytical estimate (9) with c = 4 and K =. Parameters: L + = [ ] T and L = [ ] T. and h(x) = x, y = g(x) = x. According to Theorem 4, a state observer as in (7)-(8) with L + = [l + l + ]T and L = [l l ]T for this system has the property that its estimation error converges exponentially to zero if there exist choices of the gain matrices L + and L so that all three conditions ()-() are satisfied. To find L + and L, it is first necessary to select a specific matrix measure; here we use the measure, associated to the so-called l -norm (see Appendix A). Therefore, conditions () and () translate respectively to 9 6ˆx l + ˆx l + ˆx 4 9 + 6ˆx l ˆx l ˆx 4 ]) <, with ˆx >, ]) <, with ˆx <. Selecting for simplicity l + = l =, the above inequalities are satisfied if max{ 9 6ˆx l + ˆx ; 4} <, with ˆx >, max{ 9 + 6ˆx l ˆx ; 4} <, with ˆx <. This is true if l + > 3 and l < 3. Next, from the the third condition (), we have 6ˆx 36 + (l + ] ) l )(x ˆx ) [ ] =, with ˆx =, which is verified if max{ 36 + (l + l )x ; } =, i.e. if 36 + (l + l )x <, which holds for all x if l + < l. Therefore, to satisfy all three conditions of Theorem 4, it is possible for example to select L + = [ ] T and L = [ ] T. The resulting state observer is contracting and its estimation error satisfies (9) with convergence rate c = 4. In Fig. (a) we show numerical simulations of the evolution of the states x and ˆx when an input u(t) = [ ] T sin(πt) of period T = is applied to the system. In Fig. (b) the evolution of the l -norm of the state estimation error e(t) is reported, confirming the analytical estimate (9). Example the form where A = Consider a piecewise affine (PWA) system of ẋ = [ ], b = { A x + b + Bu, if h T x > A x + b + Bu, if h T x <, (3) y = c T x, (4) [ ] [ ] [, A 3 =, b 3 =, 4] and B = [ ] T, h = [ ] T, c = [ ] T. A state observer as in (7)-(8) for this system has the structure ˆx = { Aˆx + b + L + (y ŷ) + Bu, if h Tˆx > A ˆx + b + L (y ŷ) + Bu, if h Tˆx <, (5) ŷ = c Tˆx, (6) where, for the sake of simplicity, we choose L + = L = L. Again we decide to proceed using the matrix measure induced by the l -norm. ( In this case, conditions () and () yield respectively A ( Lc T) = max{ l + l ; l + l }, and A Lc T) = max{ l + l ; 3 l + l }. It is easy to verify that choosing l = l = both measures are equal to. Condition () is verified independently of L. Hence, the designed observer (5) is contracting and the estimation error converges exponentially to zero with rate c =. In Fig. 3(a) we show numerical simulations of the evolution of the states x and ˆx when an input u(t) = 4 sin(πt) of period T = is applied to the system. In Fig. 3(b) the evolution is reported of the l -norm of the state estimation error e(t). Note that faster convergence can be obtained by choosing higher values of l and l fulfilling conditions (5)-(6). For example choosing L = [.5 ] T we obtain a convergence rate c =.5, as shown in Fig. 3(c). Example 3 Consider now a harmonic oscillator affected by Coulomb friction, described by the equations ẋ = x, ẋ = ω n x ω n Q x F f m sgn(x ) + F d m sin(ω dt), (7) y = x, (8) where x R is the position of the oscillator, x R is its velocity, ω n is its natural frequency, Q is said Q factor and is inversely proportional to the damping, m is the mass of the oscillator, F d is the amplitude of the driving force, ω d is the driving frequency and F f is the amplitude of the dry friction force which is modeled through the sign function as in (Csernák and Stépán, 6). The proposed observer for system (7)-(8) has the form ˆx = ˆx + l (x ˆx ) ˆx = ω nˆx ω n Q ˆx F f m sgn(ˆx ), ŷ = ˆx. +l (x ˆx ) + F d m sin(ω dt) Note that system (7) may be viewed as a PWA system (3) where
a) b) c).4.5. -.5 -. 3 4.6.45.3.5 3 4.5.375.5.5 3 4 Fig. 3. Panel a: Time evolution of the states x (t) (solid line) and ˆx (t) (dashed line) of Example, with initial conditions x = [.3.3] T, ˆx = [ ] T. Panel b: Norm of the estimation error e(t). The dashed line represents the analytical estimate (9) with c = and K =. Parameters: L + = L = [ ] T. Panel c: Norm of the estimation error using observer gain L = [.5 ] T. A = A = [ ], ω d ω d /Q B = [ /m] T, b = [ F f /m] T, b = [ F f /m] T, h = [ ] T, and excited by an input u(t) = F d sin(ω d t). Using the measure induced by the uniform norm (see Appendix A), conditions () and () of Theorem 4, combined, yield ]) l <, with ˆx ω n l ω n /Q, which in turn is equivalent to max { l + ; ω n /Q + ω n l } <, with ˆx. Therefore l and l must be chosen so that l >, and ω n ( + /Q) < l < ω n ( /Q). Furthermore, condition () is verified if max{; F f /m} =, which always holds because F f, m >. Numerical simulations reported in Fig. 4(a)-(b) confirm the theoretical predictions, showing that the estimation error converges to zero. In practice, the exact value of the parameter F f is not known. This implies bounded convergence of the estimation error, as stated in Remark 3. 5. CONCLUSIONS We presented an approach based on contraction for the design of state observers for a large class of nonlinear a) b).5 -.5-5 5 5 3.75.5.5 5 5 5 3 Fig. 4. Panel a: time evolution of the states x (t) (solid line) and ˆx (t) (dashed line) of Example 3, with initial conditions x = [ ] T, ˆx = [ ] T. Panel b: Norm of the estimation error e(t). The dashed line represents the analytical estimate (9) with c =. and K =. Parameters: ω n = rad/s, Q =, m = kg, F d = N, ω d = π rad/s, F f =. N, l =., l =. switched systems including those exhibiting sliding motion, such as the friction oscillator. The design methodology is based on the analysis of incremental exponential stability based on the extension of contraction theory to switched bimodal Filippov systems derived in (Fiore et al., 6). The conditions were formulated in terms of matrix measures of the Jacobians of the observer dynamics and of an additional condition on the vector fields on the discontinuity set. The theoretical results were illustrated through simple but representative examples demonstrating the effectiveness of the proposed methodology. Future work will be aimed at extending the approach to a wider class of switched systems, investigating constructive methods to design both metrics and observer gains, and reformulating the design procedure as a convex optimization problem to compute them numerically. REFERENCES Aminzare, Z. and Sontag, E.D. (4). Contraction methods for nonlinear systems: A brief introduction and some open problems. In Proc. of IEEE Conference on Decision and Control, 3835 3847. Angeli, D. (). A Lyapunov approach to incremental stability properties. IEEE Transactions on Automatic Control, 47(3), 4 4. Arcak, M. and Kokotovic, P. (). Observer-based control of systems with slope-restricted nonlinearities. IEEE Transactions on Automatic Control, 46(7), 46 5. Bonnabel, S., Astolfi, A., and Sepulchre, R. (). Contraction and observer design on cones. In Proc. of IEEE Conference on Decidion and Control, 747 75. Brogliato, B. and Heemels, W.P.M.H. (9). Observer design for Lur e systems with multivalued mappings: a passivity approach. IEEE Transactions on Automatic Control, 54(8), 996.
Csernák, G. and Stépán, G. (6). On the periodic response of a harmonically excited dry friction oscillator. Journal of Sound and Vibration, 95(3), 649 658. Dahlquist, G. (958). Stability and error bounds in the numerical integration of ordinary differential equations. PhD Dissertation; reprinted in Transactions of the Royal Institute of Technology, (3). Desoer, C. and Haneda, H. (97). The measure of a matrix as a tool to analyze computer algorithms for circuit analysis. IEEE Transactions on Circuit Theory, 9(5), 48 486. di Bernardo, M. and Fiore, D. (4). Incremental stability of bimodal Filippov systems in R n. In Proc. of IEEE Conference on Decision and Control, 4679 4684. di Bernardo, M. and Fiore, D. (6). Switching control for incremental stabilization of nonlinear systems via contraction theory. In Proc. of European Control Conference, 54 59. di Bernardo, M. and Liuzza, D. (3). Incremental stability of planar Filippov systems. In Proc. of European Control Conference, 376 37. di Bernardo, M., Liuzza, D., and Russo, G. (4). Contraction analysis for a class of nondifferentiable systems with applications to stability and network synchronization. SIAM Journal on Control and Optimization, 5(5), 33 37. Dinh, T.N., Bonnabel, S., and Sepulchre, R. (3). Contraction-based design of positive observers. In Proc. of IEEE Conference on Decision and Control, 6574 6579. Doris, A., Juloski, A.L., Mihajlovic, N., Heemels, W., van de Wouw, N., and Nijmeijer, H. (8). Observer designs for experimental non-smooth and discontinuous systems. IEEE Transactions on Control Systems Technology, 6(6), 33 33. Filippov, A.F. (988). Differential equations with discontinuous righthand sides. Kluwer. Fiore, D., Hogan, S.J., and di Bernardo, M. (6). Contraction analysis of switched systems via regularization. Automatica, 73, 79 88. Forni, F. and Sepulchre, R. (4). A differential Lyapunov framework for contraction analysis. IEEE Transactions on Automatic Control, 59(3), 64 68. Forni, F., Teel, A.R., and Zaccarian, L. (3). Follow the bouncing ball: Global results on tracking and state estimation with impacts. IEEE Transactions on Automatic Control, 58(6), 47 485. Heemels, W.P.M.H., Camlibel, M.K., Schumacher, J.M., and Brogliato, B. (). Observer-based control of linear complementarity systems. International Journal of Robust and Nonlinear Control, (), 93 8. Jouffroy, J. (5). Some ancestors of contraction analysis. In Proc. of IEEE Conference on Decision and Control and European Control Conference, 545 5455. Juloski, A.L., Heemels, W.P.M.H., and Weiland, S. (7). Observer design for a class of piecewise linear systems. International Journal of Robust and Nonlinear Control, 7(5), 387 44. Lohmiller, W. and Slotine, J.J.E. (998). On contraction analysis for non-linear systems. Automatica, 34(6), 683 696. Lohmiller, W. and Slotine, J.J.E. (). Nonlinear process control using contraction theory. AIChE Journal, 46(3), 588 596. Lu, W. and di Bernardo, M. (6). Contraction and incremental stability of switched Carathéodory systems using multiple norms. Automatica, 7, 8. Manchester, I.R. and Slotine, J.J.E. (4a). Outputfeedback control of nonlinear systems using control contraction metrics and convex optimization. In Proc. of Australian Control Conference, 5. Manchester, I.R. and Slotine, J.J.E. (4b). Control contraction metrics and universal stabilizability. In Proc. of IFAC World Congress, volume 9, 83 88. Pavlov, A., Pogromsky, A., van de Wouw, N., Nijmeijer, H., and Rooda, K. (4). Convergent dynamics, a tribute to Boris Pavlovich Demidovich. Systems & Control Letters, 5(3), 57 6. Pavlov, A., Pogromsky, A., van de Wouw, N., and Nijmeijer, H. (7). On convergence properties of piecewise affine systems. International Journal of Control, 8(8), 33 47. Piiroinen, P.T. and Kuznetsov, Y.A. (8). An eventdriven method to simulate Filippov systems with accurate computing of sliding motions. ACM Transactions on Mathematical Software, 34(3), 3: 3:4. Rajamani, R. (998). Observers for Lipschitz nonlinear systems. IEEE Transactions on Automatic Control, 43(3), 397 4. Russo, G., di Bernardo, M., and Sontag, E.D. (). Global entrainment of transcriptional systems to periodic inputs. PLoS Comput Biol, 6(4), e739. Sotomayor, J. and Teixeira, M.A. (996). Regularization of discontinuous vector fields. In Proc. of International Conference on Differential Equations, Lisbon, 7 3. Utkin, V.I. (99). Sliding modes in control and optimization. Springer. van de Wouw, N. and Pavlov, A. (8). Tracking and synchronisation for a class of PWA systems. Automatica, 44(), 99 95. Vidyasagar, M. (). Nonlinear systems analysis. SIAM, nd edition. Zemouche, A. and Boutayeb, M. (3). On LMI conditions to design observers for Lipschitz nonlinear systems. Automatica, 49(), 585 59. Appendix A. MATRIX MEASURES The matrix measure (Dahlquist, 958; Vidyasagar, ) associated to a matrix A R n n is the function ( ) : R n n R defined as (A) = lim h + ( I + ha ) /h. See (Vidyasagar, ; Desoer and Haneda, 97) for a list of properties of this measure. The most common used matrix measures are those associated to the l -norm, the Euclidean norm and the uniform norm, and they are defined as follows: (A) = max j [ a jj + ] a ij, i j ) (A) = λ max ( A + A T (A) = max i, ]. [ a ii + a ij j i