linear positive systems @xcite have been recently the subject of an increasing attention because of their natural ability to represent many real - world processes such as , among others , communication networks @xcite , biological networks @xcite , epidemiological networks @xcite , disease dynamics @xcite , etc . besides their applicative potential
, linear positive systems have been shown to exhibit a number of interesting theoretical properties of independent interest .
for instance , it is now well - known that linear copositive lyapunov functions can be used in order to formulate exact stability conditions taking the form of linear programs @xcite .
the design of structured and bounded state - feedback controllers @xcite and certain classes of static output feedback controllers @xcite are known to be convex and hence easily tractable .
the @xmath0-gains for @xmath1 can be exactly computed using convex programming and these gains are identical to the @xmath2-norm of the static matrix - gain of the system @xcite .
the famous kalman - yakubovich - popov lemma has been shown to admit a linear formulation in this setting @xcite .
robust analysis results also nicely extends and simplifies in this context , and often becomes necessary and sufficient criteria for stability @xcite .
their generalization to delay - systems with discrete - delays also led to the surprise that the system is stable if and only if the system with zero delay is stable @xcite .
design - wise , it has been shown in @xcite that the state - feedback controller that minimizes the @xmath3-gain of the closed - loop system is independent of the matrices driving the disturbance input to the state and the output of the system .
symmetrically , it is shown in @xcite that the observer - gain that minimizes the @xmath4-gain of map between the disturbance input and the observation error is independent of the output matrices of the error system .
extensions to deterministically @xcite or stochastically @xcite switched systems have also been considered .
we consider here the case of linear positive impulsive systems , a class of systems that seems to have been quite overlooked until now as only very few results can be found ; see e.g. @xcite .
such systems can be used to represent certain classes biochemical , population or epidemiological models having deterministic jumps in their dynamics .
they can also be used to represent processes that can be represented as linear positive switched systems ; see e.g. @xcite for some examples including epidemiology @xcite , traffic congestion models @xcite , etc .
impulsive systems are also well known to be able to exactly represent sampled - data systems as emphasized in @xcite .
the goal of this paper is hence the derivation of novel stability and stabilization conditions for linear positive impulsive systems using the concepts of arbitrary , constant , minimum , maximum and range dwell - times .
the concept of _ minimum dwell - time _ has been introduced by morse in @xcite in order to formulate stability conditions for general ( i.e. not necessarily positive ) switched systems .
the concept of _ average dwell - time _ has been proposed in @xcite in order to obtained less conservative conditions that by using minimum dwell - time conditions . since then , a large body of the literature has been focusing on these concepts as a way to efficiently characterize the stability of switched systems or , more generally , the stability of hybrid systems ; see e.g. @xcite . the notion of minimum dwell - time has been revisited in @xcite where novel sufficient lmi conditions derived from mode - dependent quadratic lyapunov functions were proposed .
based on a theoretical result proved in @xcite , these conditions were later extended and made necessary and sufficient in @xcite through the consideration of mode - dependent homogeneous lyapunov functions . unfortunately , these conditions were inapplicable to uncertain systems and to control design because of their complex nonlinear dependency in the matrices of the system .
this problem motivated the introduction of the so - called _ looped - functionals _
, a particular class of indefinite ( i.e. not necessarily positive definite ) functionals having the advantage of reformulating the complex conditions of @xcite into conditions being linear in the matrices of the system ; see e.g. @xcite , thereby extending the scope of the conditions to uncertain and nonlinear systems . yet
, these conditions were difficult to apply in the context of control design because of the presence of multiple products between decision matrices and the matrices of the system ; see e.g. @xcite . _
clock - dependent conditions _ have been shown to provide an essential framework for solving this latter problem as they produce stability conditions that are linear / convex in the matrices of the system and can be used for design purposes .
their computational complexity has also been shown to be much lower than that of looped - functionals @xcite . since then
, clock - dependent conditions have been used for the analysis and control of switched , impulsive , sampled - data and lpv systems ; see e.g. @xcite .
the first part of the paper is similar to the ones in @xcite where stability conditions are formulated in terms of the decrease of a lyapunov function of a given type . unlike in the previous references where quadratic lyapunov functions are involved ,
we exploit here the positivity of the system and consider a linear copositive lyapunov function @xcite .
the resulting conditions are stated in terms of finite - dimensional linear programs or semi - infinite dimensional linear programs , which are then relaxed into clock - dependent conditions using a lifting approach similar to that of @xcite .
since linear copositive lyapunov functions are used here , the clock - dependent conditions consist of infinite - dimensional linear programs .
this has to be contrasted with the fact that when quadratic lyapunov functions are used , clock - dependent conditions take the form of infinite - dimensional semidefinite programs , which may be harder to solve that their linear counterpart .
three possible ways for efficiently checking these conditions are then proposed .
the first one relies on a discretization approach which is largely inspired from @xcite and where the infinite - dimensional decision variable is assumed to be continuous and piecewise linear .
this method has also been considered , in turn , in @xcite . by doing so
, the infinite - dimensional program becomes finite - dimensional and can be solved using conventional algorithms such as interior point methods ; see e.g. @xcite .
the second method is based on handelman s theorem @xcite which characterizes the positivity of a given polynomial on a compact polytope by formulating it as a nonnegative linear combination of products of the ( affine ) basis functions that describe the polytope .
this result has been applied in various contexts @xcite and , notably , for characterizing the robust stability of uncertain linear positive systems in @xcite .
an important property of this approach is that the obtained characterization for the positivity of the polynomial can be exactly formulated as a finite - dimensional linear program , which can again be solved using well - known approaches .
finally , the last one is based on putinar s positivstellensatz @xcite which characterizes the positivity of a given polynomial on a compact semialgebraic set by formulating it as a weighted linear combination of the basis functions that describe the set and where the weights are sum of squares polynomials .
the resulting problem takes , in this case , the form of a finite - dimensional semidefinite program @xcite that can be solved using standard semidefinite programming solvers such as sedumi @xcite or sdpt3 @xcite used in conjunction with the package sostools @xcite .
it is notably emphasized that these relaxations are asymptotically exact meaning that when the discretization order , the number of products of basis functions or the degree of the sum of squares weights are sufficiently large , then the relaxed problem is feasible if the original one is .
several examples are considered in order to demonstrate the practicality of the relaxed conditions and to compare them in terms of number of variables and solving time .
the results are then extended to control design by considering the clock - dependent conditions and the dual impulsive system @xcite . by finally exploiting the possibility of formulating a switched system as an impulsive system
, we derive a number of stability conditions for linear positive switched systems .
notably , we recover stability conditions similar to those in @xcite which are the positive counterpart of those obtained in @xcite in the context of general linear switched systems while some other ones seem to be new .
novel stabilization conditions under arbitrary , minimum and mode - dependent range dwell - time conditions are then obtained and illustrated through few examples
. + * outline .
* the structure of the paper is as follows : in section [ sec : preliminary ] preliminary definitions and results are given .
section [ sec : stab ] is devoted to the dwell - time stability analysis of linear positive impulsive systems while section [ sec : stabz ] addresses their stabilization .
these results are then applied to switched systems in section [ sec : switched ] .
+ * notations . *
the cone of positive ( nonnegative ) vectors of dimension @xmath5 are denoted by @xmath6 ( @xmath7 ) .
the set of diagonal matrices of dimension @xmath5 is denoted by @xmath8 and that of diagonal matrices with positive diagonal elements by @xmath9 . for a set of matrices @xmath10 of compatible dimensions ,
we define @xmath11 to be the block diagonal matrix with the @xmath12 s as diagonal blocks and by @xmath13 the matrix consisting of the @xmath12 s stacked in column .
the @xmath5-dimensional vector of ones is denoted by @xmath14 .
let us consider here the following class of linear impulsive system : @xmath15 where @xmath16 are the state of the system and its initial condition , respectively .
the notation @xmath17 is a shorthand for @xmath18 , i.e. the trajectories are assumed to be left - continuous .
the sequence of impulse instants @xmath19 is assumed to be strictly increasing and to grow unboundedly .
the sequence hence does not admit any accumulation point , thereby excluding any zeno behavior for the system trajectories .
note that this assumption is not restrictive for the current paper as we only focus here on dwell - time results .
we define the dwell - times as the values @xmath20 , @xmath21 .    the following result states necessary and sufficient conditions for the impulsive system to be positive :    the following statements are equivalent :    a.   the system is positive ; i.e. for any @xmath22 , we have that @xmath23 for all @xmath24 .
b.   the matrix @xmath25 is metzler ( all its off - diagonal entries are nonnegative ) and @xmath26 is nonnegative .
the following result establishes conditions for which a metzler matrix is hurwitz stable :    [ prop : hurwitz ] let @xmath27 be a metzler matrix .
then , the following statements are equivalent :    a.   @xmath25 is hurwitz stable ; b.   @xmath25 is nonsingular and @xmath28 ; c.   there exists a @xmath29 such that @xmath30 ; d.   there exists a @xmath31 such that @xmath32 ;    similarly , the following result establishes conditions for which a nonnegative matrix is schur stable :    [ prop : schur ] let @xmath33 be a nonnegative matrix .
then , the following statements are equivalent :    a.   @xmath34 is schur stable ; b.   @xmath35 is hurwiz stable ; c.   there exists a @xmath29 such that @xmath36 ; d.   there exists a @xmath31 such that @xmath37 ;    the following generic stability result will be important in the sequel :    [ th : equivalence ] let the sequence @xmath38 for the system be given and assume that it satisfies the strict increase and unboundedness conditions .
then , the following statements are equivalent :    a.   [ item : equi : a ] the impulsive system is asymptotically stable . b.   [ item : equi : a2 ] the dual impulsive system @xmath39 is asymptotically stable . c.   [ item : equi : b ] the discrete - time system @xmath40 is asymptotically stable . d.   [ item : equi : c ] the discrete - time system @xmath41 is asymptotically stable .
e.   [ item : equi : d ] the discrete - time system @xmath42 is asymptotically stable .
f.   [ item : equi : e ] the discrete - time system @xmath43 is asymptotically stable .
* proof of @xmath44 . *
immediate .
noindent**proof of @xmath44 . * * to prove this implication , we need to consider two cases .
the first one is the case where the impulses are not persistent , i.e. there exists a @xmath45 such that @xmath46 for all @xmath47 and @xmath48 . in this case , the system is asymptotically stable if and only if @xmath25 is hurwitz stable .
if @xmath48 , then @xmath49 for all @xmath50 .
assuming that @xmath51 , this is equivalent to saying that @xmath25 is hurwitz stable .
let us now consider the case where the jumps are persistent , i.e. for any sequence @xmath38 we choose , there exists a @xmath52 such that @xmath53 for all @xmath54 .
note first that independently of the chosen sequence , we have that @xmath55 for some @xmath56 and where the norm is any vector norm .
this then implies that @xmath57 since the embedded discrete - time system is asymptotically stable , then @xmath58 as @xmath59 .
this then implies that @xmath60}||x(t_k+s)||\to0\ \textnormal{as}\ k\to\infty\ ] ] which proves the asymptotic stability of the system in the persistent jump case .
the proof is complete .
* proof of @xmath61 . *
note that while the system describes the evolution of the sequence @xmath62 , the system describes the evolution of the sequence @xmath63 ; i.e. @xmath64 .
using the fact that @xmath65 , then we get that @xmath66 hence , clearly , we have that @xmath67 as @xmath59 if and only if @xmath68 as @xmath59 .
the proof is complete .
* proof of @xmath61 .
* let assume first that the sequence @xmath38 is such that @xmath69 $ ] for some @xmath70 and define the following set of matrices @xmath71\right\}.\ ] ] clearly , this set is bounded since the interval @xmath72 $ ] is compact .
we know that the system is asymptotically stable if and only if the joint spectral radius condition @xmath73 holds ; see e.g. @xcite .
it is known that we have @xmath74 where the latter mean that we consider transposes of the matrices in @xmath75 ; see e.g. @xcite .
then the conclusion follows .
assuming now that we have no upper - bound for the dwell - time , we can assume that the matrix @xmath25 is hurwitz stable ( see the beginning of the proof that ( b ) implies ( a ) ) , which then implies that @xmath76 is compact and the same arguments as before apply .
the proof is complete .
* proof of @xmath61 .
* this follows from the same arguments .
* proof of @xmath61 .
* this follows from the same arguments as for the proof that @xmath61 .
the meaning of the above result is that we can choose the most convenient system to work with in order to derive stability and stabilization conditions for the original impulsive system .
this result is also important as it will be used to demonstrate that despite the systems can be all equivalent in terms of stability , the obtained stability conditions will not be necessarily equivalent .
this point is discussed for the case of positive switched systems in @xcite and the same discussion remains valid in the case of linear positive impulsive systems .
notably , the difference between the results obtained using the system and the  swapped " version will be emphasized in section [ sec : imp : examples ] and in section [ sec : sw : examples ] .
the interest for using the dual system and its  swapped " version will be emphasized in the sections related to stabilization ; i.e. section [ sec : stabz ] and section [ sec : switched ] .
note also that when the sequence @xmath38 is such that @xmath77 for all @xmath78 , @xmath79 , then all the discrete - time systems in theorem [ th : equivalence ] become time - invariant after time @xmath80 .
hence , the stability conditions reduce to standard stability conditions for linear time - invariant discrete - time systems and will all be equivalent with each other .
this result can be easily be understood through the fact that @xmath81 , where @xmath82 denotes the spectral radius .
we derive here several stability results for linear positive impulsive systems .
the case of arbitrary dwell - time is considered first in section [ sec : imp : arbitrary ] and is followed by stability results under constant dwell - time ( @xmath83 ) , minimum ( @xmath84 ) , maximum ( @xmath85 ) and range dwell - time ( @xmath86 $ ] ) in section [ sec : imp : constant ] , section [ sec : imp : minimum ] , section [ sec : imp : maximum ] and section [ sec : imp : range ] , respectively .
computational results are given in section [ sec : imp : computation ] together with some discussions regarding their conservatism .
comparative examples are finally given in section [ sec : imp : examples ] .
let us consider first the stability under arbitrary dwell - time ( @xmath87 ) , with the additional condition that the sequence @xmath88 grows unboundedly .
we then have the following result :    [ th : arbdt ] the following statements are equivalent :    a.   there exists a @xmath89 such that @xmath90 b.   we have that @xmath91    moreover , when one of the above statements holds , then the system is asymptotically stable under arbitrary dwell - time ; i.e. @xmath87 .
the first condition in is equivalent to saying that @xmath92 for some @xmath93 whereas the second one is equivalent to saying that @xmath94 for some @xmath95 and all @xmath54 .
these conditions , all together , imply that @xmath96 as @xmath97 regardless the impulse sequence @xmath38 . to prove the equivalence with ( ii ) , it is enough to remark that the conditions in statement ( i ) are identical to stability conditions for positive switched systems using a common copositive linear lyapunov functions @xcite .
the result follows from the latter paper .
[ rem : dual_arbdt][rem : arbdt2 ] dual conditions can be obtained by substituting the matrices of the dual system inside the conditions of theorem [ th : arbdt ] . in such a case ,
the conditions of statement ( a ) become @xmath98 whereas that of statement ( b ) become @xmath99    [ rem : arbdt1_flow ] in the case of persistent flowing , i.e. the flow will never stops , the conditions in theorem [ th : arbdt ] , ( a ) , can be relaxed to @xmath100 and those in remark [ rem : arbdt2 ] , ( a ) , to @xmath101    the example below illustrates the discussion below theorem [ th : equivalence ] about the non - equivalence between the stability conditions obtained from the impulsive system and its dual :    let us consider the system with matrices @xmath102 there is no vector @xmath29 such that the conditions hold since @xmath30 implies that @xmath103 while @xmath104 implies that @xmath105 , yielding then a contradiction . on the other hand , we can readily see that @xmath106 solves the conditions , emphasizing then the discrepancy between the conditions of theorem [ th : arbdt ] and remark [ rem : arbdt2 ] .
we consider now the case of constant dwell - times  the case where @xmath107 , for some @xmath108 and for all @xmath21  or , in other words , the case where jumps occur periodically .
we have the following result :    [ th : cstdt2 ] the following statements are equivalent :    a.   [ item : th1:0 ] the impulsive system is asymptotically stable under constant dwell - time @xmath109 .
b.   [ item : th1:1 ] the copositive linear form @xmath110 , @xmath29 is a discrete - time lyapunov function for the @xmath109-periodic impulsive system in the sense that the inequality @xmath111 holds for some @xmath93 , all @xmath112 and all @xmath21 . c.   [ item : th1:2 ] there exists a vector @xmath89 such that the inequality @xmath113<0\ ] ] holds or , equivalently ,
the matrix @xmath114 is schur stable .
d.   [ item : th1:3 ] there exists a vector @xmath89 such that the inequality @xmath115\lambda<0\ ] ] holds or , equivalently , the matrix @xmath116 is schur stable .
e.   [ item : th1:4 ] there exist a differentiable vector - valued function @xmath117\mapsto\mathbb{r}^n$ ] , @xmath118 , and a scalar @xmath119 such that the inequalities @xmath120 and @xmath121 hold for all @xmath122 $ ] . f.   [ item : th1:5 ] there exist a differentiable vector - valued function @xmath123\mapsto\mathbb{r}^n$ ] , @xmath124 , and a scalar @xmath119 such that the lmis @xmath125 and @xmath126 hold for all @xmath122 $ ] .
* proof of @xmath127 :* this follows from theorem [ th : equivalence ] and proposition [ prop : schur ] .
* proof of @xmath127 :* it is immediate to see that the left - hand side of is given by @xmath128^\t x(t_k).\ ] ] the equivalence between the two statements immediately follows .
* proof of @xmath127 :* this follows from proposition [ prop : hurwitz ] or proposition [ prop : schur ] .
* proof of @xmath127 :* the proof between the two statements follows from the definition that @xmath129 .
alternatively , similar arguments as in the proof that @xmath127 can also be considered .    the main advantages of the conditions of statements and in the above results are the following .
first of all , the conditions are affine in the matrices @xmath25 and @xmath26 of the system , allowing then for an immediate extension to uncertain matrices and to control design ( the latter requiring some additional steps that will be detailed in section [ sec : stabz ] ) .
a particularity of the approach is that we only need the vectors @xmath130 and @xmath131 to be positive while most of the lyapunov approaches would require the whole vector - valued functions @xmath132 and @xmath133 to be positive on their domain ( see e.g. @xcite ) , a constraint that is computationally way more complex .
the explanation for that in spite of the fact that the conditions in statements and statement look like continuous - time stability conditions , they are in fact lifted discrete - time stability condition for which only the behavior of the inequalities at the extremal points of the domain of the lifting variable ( i.e. the points @xmath134 and @xmath135 ) are actually meaningful .
the same discussion is made in @xcite in the context of looped - functionals and in @xcite in the context of clock - dependent conditions .
the main drawback of this procedure lies in the increase of the computational complexity as the conditions of statement and statement are infinite - dimensional linear programming conditions which are in general hard to solve .
note , however , that they should be less complex than the infinite - dimensional semidefinite programming conditions obtained in @xcite . verifying these conditions in an efficient way will be the topic of section [ sec : imp : computation ] .
[ sec : imp : minimum ]    let us consider now the minimum dwell - time case . in this setup ,
the dwell - times are assumed to satisfy the condition @xmath84 , @xmath21 .
even though we do not assume here the existence of an upper - bound for the dwell - times a priori , we will restrict ourselves to the case of persistent impulses ; i.e. for any chosen sequence @xmath136 , there exists an upper - bound for the dwell - times .
hence , we exclude the case that @xmath48 for some @xmath137 , which is not restrictive , as in such a case , the system would become a standard continuous - time linear positive system .
the following result can be seen as a  positive systems " counterpart of the results in @xcite and provides a stability condition in terms of the discrete - time system :    [ th : mindt2 ] the following statements are equivalent :    a.   the linear form @xmath110 , @xmath29 is a lyapunov function for the system in the sense that @xmath138 and @xmath139 hold for some @xmath140 , for all @xmath141 and for any sequence @xmath19 verifying @xmath84 , @xmath21 .
b.   there exists a @xmath89 such that the inequalities @xmath142 and @xmath143 hold for all @xmath144 . c.   there exists @xmath89 such that the conditions @xmath142 and @xmath145 hold .
d.   the pair @xmath146 admits a common linear copositive lyapunov function .
e.   we have that @xmath147 f.   there exist a matrix function @xmath117\mapsto\mathbb{r}^n$ ] , @xmath118 , and a scalar @xmath119 such that the inequalities @xmath148 @xmath149 and @xmath150 hold for all @xmath122 $ ] .
g.   there exist a matrix function @xmath123\mapsto\mathbb{r}^n$ ] , @xmath124 , and a scalar @xmath119 such that the inequalities @xmath151 @xmath152 and @xmath153 hold for all @xmath122 $ ] .
moreover , when one of the above statements holds , the impulsive system ( [ eq : mainsyst ] ) is asymptotically stable under minimum dwell - time @xmath154.@xmath155    the proof that ( a ) is equivalent to ( b ) follows from the same arguments as in the proof of theorem [ th : mindt2 ] .
the proof that ( b ) implies ( c ) is immediate . to prove the converse , we assume that @xmath156 , which implies that @xmath157 where we have used the fact that @xmath158 for all @xmath159
. note also that since @xmath160 , then we have that @xmath161 for all @xmath159 .
adding finally the latter inequality combined to implies that @xmath162 holds for all @xmath159 , which proves the implication .
the rest of the proof can be carried as for theorem [ th : mindt2 ] .
the above results interestingly connect the problem of finding a common linear copositive lyapunov functions for a linear positive switched system with two subsystems to the problem of establishing whether a linear positive impulsive system is stable under some minimum dwell - time constraint .
indeed , in the context of theorem [ th : mindt2 ] , the matrix of the first subsystem is @xmath25 and that of the second subsystem is @xmath163 . when the swapped system is considered ( see remark [ rem : swapped_minimum ] ) , the matrix of the first subsystem is also @xmath25 but that of the second subsystem is @xmath164 . in this regard , characterizing the existence of a linear copositive lyapunov function can be achieved using existing methods ; such as those developed in @xcite .
however , these conditions are difficult to consider when the matrix @xmath25 is uncertain or when design is the main goal .
the statements ( f ) and ( g ) provide alternative equivalent conditions that are affine in @xmath25 and @xmath26 and can then be hence considered for uncertain systems and for design purposes .
a last interesting point is that the results for arbitrary dwell - time can be retrieved by letting @xmath165 in the above theorem as also noticed in @xcite .
we consider now the maximum dwell - time case . in this setup , the dwell - times are assumed to satisfy a maximum dwell - time condition ; i.e. @xmath85 , @xmath21 .
the following result can be understood as being the  positive systems " version of the results in @xcite and is the maximum dwell - time counterpart of theorem [ th : mindt2 ] :    [ th : maxdt2 ] the following statements are equivalent :    a.   the linear form @xmath110 , @xmath29 is a lyapunov function for the system in the sense that @xmath166 and @xmath167 hold for some @xmath140 , all @xmath141 and any sequence @xmath168 .
b.   there exists @xmath89 such that the inequalities @xmath169 and @xmath170 hold for all @xmath144 . c.   there exists @xmath89 such that the conditions @xmath169 and @xmath171 hold .
d.   the pair @xmath172 admits a common linear copositive lyapunov function .
e.   we have that @xmath173 f.   there exist a matrix function @xmath117\mapsto\mathbb{r}^n$ ] , @xmath118 , and a scalar @xmath119 such that the inequalities @xmath174 @xmath149 and @xmath150 hold for all @xmath122 $ ] .
g.   there exist a matrix function @xmath123\mapsto\mathbb{r}^n$ ] , @xmath124 , and a scalar @xmath119 such that the inequalities @xmath175 @xmath152 and @xmath153 hold for all @xmath122 $ ] .
moreover , when one of the above statements holds , the impulsive system ( [ eq : mainsyst ] ) is asymptotically stable under maximum dwell - time @xmath154.@xmath155    the proof follows from the same lines as the proof of theorem [ th : mindt2 ] .
[ sec : imp : range ]    we address now the range dwell - time case , that is the case when @xmath86 $ ] , @xmath21 , for some prescribed bounds @xmath176 verifying @xmath70 .
we then have the following result :    [ th : rangedt ] the following statements are equivalent :    a.   [ item : th3:1 ] the copositive linear form @xmath110 , @xmath29 is a discrete - time lyapunov function for the impulsive system in the sense that the inequality @xmath177 holds for some @xmath93 , all @xmath112 , @xmath86 $ ] and all @xmath21 .
b.   [ item : th3:2 ] there exists a vector @xmath89 such that the inequality @xmath178<0\ ] ] holds for all @xmath179 $ ] . c.   [ item : th3:4 ]
there exist a differentiable vector - valued function @xmath117\mapsto\mathbb{r}^n$ ] , @xmath180 , and a scalar @xmath119 such that the inequalities @xmath181 and @xmath182 hold for all @xmath183 $ ] and @xmath179 $ ] .
it is readily seen that the lyapunov condition is equivalent to the inequality .
the proof of the equivalence between the statements ( b ) and ( c ) follows from exactly the same arguments as in the proof of theorem [ th : cstdt2 ] .
it is thus omitted .
it seems important to carefully address the problem of verifying the infinite - dimensional conditions of the previous results .
we detail here three methods for accurately and efficiently verifying them . only the relaxation of the conditions in theorem [ th : cstdt2 ] , are considered here for simplicity .
however , relaxed finite - dimensional results for the other conditions can be easily obtained using the same procedure .
this approach is based on the piecewise linear approximation of the function @xmath184 or @xmath185 in the results .
this method has been initially proposed in @xcite in order to relax infinite - dimensional lmi conditions arising in the analysis of switched systems into finite - dimensional lmi conditions .
below is the approximation of the conditions in theorem [ th : cstdt2 ] , statement :    [ prop : discretization ] the following statements are equivalent :    a.   for some @xmath186 , the function @xmath187,\ i=0,\ldots , d_g-1\ ] ] verifies the conditions of theorem [ th : cstdt2 ] , statement .
b.   there exist vectors @xmath188 , @xmath189 , such that the following conditions @xmath190 @xmath191 and @xmath192 hold for all @xmath193 .
the proof follows from direct substitutions .
this approach is based on putinar s positivstellensatz @xcite and is solved using semidefinite programming methods , as initially proposed in @xcite . before stating the main result , we need first to define some terminology .
a multivariate polynomial @xmath194 is said to be a sum - of - squares ( sos ) polynomial if it can be written as @xmath195 for some polynomials @xmath196 .
a polynomial vector @xmath197 is said to _ componentwise sum - of - squares _
( csos ) if each of its components is an sos polynomial .
checking whether a polynomial is sos can be exactly cast as a semidefinite program @xcite that can be easily solved using semidefinite programming solvers such as sedumi @xcite or sdpt3 @xcite .
the package sostools @xcite can be used to formulate the sos program in a convenient way .
this approach has been used for the analysis and control of impulsive systems @xcite , switched systems @xcite , pseudo - periodic systems with impulses @xcite and stochastic impulsive systems @xcite .
below is the approximation of the conditions in theorem [ th : cstdt2 ] , statement :    [ prop : sos ] let @xmath198 , @xmath119 and @xmath199 be given and assume that there exists a polynomial vectors @xmath200 and @xmath201 of degree @xmath202 such that    a.   @xmath203 is csos , b.   @xmath204 ( or is csos ) , c.   @xmath205 is csos and d.   @xmath206 ( or is csos ) .
then , the conditions of theorem [ th : cstdt2 ] , statement hold with the same @xmath132 and the system is asymptotically stable under constant dwell - time @xmath154 .    clearly , the condition in ( ii ) is equivalent to saying that @xmath207 and that in ( iv ) is equivalent to the condition .
we then consider the condition in ( iii ) , which implies that @xmath208 for all @xmath209 . since @xmath203 is csos
, then we have that @xmath210 for all @xmath209 and , therefore , this implies that @xmath211 for all @xmath212 $ ] which is exactly the condition .
the proof is complete .
this approach is based on handelman s theorem @xcite that states that if a polynomial is positive on a compact polytope of the form @xmath213 , then it can be expressed as a nonnegative linear combinations of products of @xmath214 .
this elegant result has been considered in @xcite for deriving finite - dimensional linear programs establishing the stability of linear uncertain positive systems or in @xcite for the construction of lyapunov functions and the verification of polynomials optimization problems .
below is the approximation of the conditions in theorem [ th : cstdt2 ] , statement :    [ prop : handelman ] let @xmath215 , @xmath119 and @xmath199 be given and assume that there exists a vector - valued polynomial @xmath200 of degree @xmath216 and parameters @xmath217 , @xmath218 such that    a.   @xmath219 for all @xmath220 , b.   @xmath221 , c.   the polynomial @xmath222 has nonnegative coefficients and d.   @xmath206 .
then , the conditions of theorem [ th : cstdt2 ] , statement hold with the same @xmath132 and the system is asymptotically stable under constant dwell - time @xmath154 .
the proof follows from the same line as for proposition [ prop : sos ] with the difference that we consider here handelman s theorem in place of putinar s positivstellensatz .
following similar arguments as in the proof of proposition [ prop : sos ] , if the conditions in statement ( i ) and ( iii ) are satisfied then we have that @xmath211 for all @xmath212 $ ] . the fact that the condition is stated in terms of the nonnegativity of the coefficients ( as opposed to the nonnegativity of the eigenvalues of the gram matrix for the sum of squares relaxation ) comes from the fact that handelman s theorem deals with an equality between the supposedly positive polynomial and a nonnegative sum of products of the basis functions of the polytope
however , we are not interested here in finding an equality but finding a positive lower bound for the polynomial that can be expressed as a nonnegative sum of products of basis functions .
the polynomial in statement ( iii ) can be rewritten in the form @xmath223 where @xmath223 .
using now the fact that @xmath224 , then we have that@xmath223 for all @xmath224 if @xmath225 .
it is also important to stress that the third condition could have also been formulated in terms of an sos condition instead of in terms of a condition on the coefficients on some polynomial .
however , the former would lead to semidefinite programming conditions as opposed to linear programming ones for the latter .
note also that the latter condition shares similarities with the polya s theorem @xcite .
we illustrate here the efficiency obtained results and compare them with existing ones .
comparisons between the proposed relaxation results are also made for completeness .
let us consider the system with the matrices
@xcite @xmath226 clearly , the matrix @xmath25 is hurwitz stable whereas the @xmath26 is schur anti - stable .
simple calculations show that the matrix @xmath227 is schur stable provided that @xmath228 . using the minimum dwell - time stability result of theorem [ th : mindt2 ] , ( a )
, we find the value 0.2311 for the estimated minimum dwell - time value .
as this value coincides with the lower - bounded obtained for the constant dwell - time case , then we can conclude that the estimated minimum dwell - time is exact .
this proves that the proposed method is way less conservative than the one proposed in @xcite where the value 1.0986 is found for the minimum dwell - time .
let us consider the system with the matrices @xmath229    when @xmath230 , the conditions of theorem [ th : mindt2 ] , ( a ) and of remark [ rem : swapped_minimum ] both yield the value @xmath231 for the minimum dwell - time .
this value was also found in @xcite using conditions based on quadratic lyapunov functions .
this value also coincides with the minimal value obtained in the constant dwell - time case , emphasizing then the exactness of the computed minimum dwell - time .    on the other hand , when @xmath232 , the conditions of remark [ rem : swapped_minimum ] yield the value 0.4290 while the value 0.3615 is obtained using the conditions of theorem [ th : mindt2 ] , ( a ) .
the latter value is also found when using the conditions based on the use of quadratic lyapunov functions stated in @xcite and coincides with the value obtained in the constant dwell - time case .
hence exactness of the estimate of the minimum dwell - time is again proved .
moreover , the discrepancy between the results obtained with the conditions in theorem [ th : mindt2 ] and in remark [ rem : swapped_minimum ] demonstrate the non - equivalence of the conditions and in the necessity of considering both in order to obtain less conservative result .
note that this phenomenon does not occur when considering quadratic lyapunov function ; see e.g. @xcite .
a comparison of the different computational methods described in section [ sec : imp : computation ] is performed and the results are summarized in table [ table : mindt1 ] where we can see that all the approaches are able to approach the value obtained using the condition of theorem [ th : mindt2 ] , ( b ) , when the discretization order @xmath233 or the polynomial degrees @xmath234 or @xmath216 increase .
however , this comes at the price of an increase of the computational complexity and the solving time .
.comparison of the numerical results obtained for the minimum dwell - time using the discretization approach ( section [ sec : discretization ] ) , the sum of squares approach ( section [ sec : sos ] ) and the handelman - based approach ( section [ sec : handelman ] ) . for each method
, we give the estimate for the minimum dwell - time @xmath154 , the number of primal / dual variables of the optimization problem and the overall solving time in seconds . for fairness ,
the solver sedumi is used for solving all the optimization problems . [
cols="^,^,^,^,^,^,^",options="header " , ]     let us consider the system with the matrices @xmath235 note that both subsystems are unstable .
the goal is then to design a controller of the form such that the closed - loop system is positive and stable under minimum dwell - time @xmath236 .
we consider theorem [ th : mindt_switched_2 ] and solve the conditions using the sos relaxation with @xmath237 ( degree of @xmath238 and @xmath239 is one while the additional polynomials are of degree 2 ) .
the obtained semidefinite program has 180 primal variables and 71 dual variables .
solving the conditions yields the controller matrices @xmath240 for simulation purposes , we generate a random sequence of dwell - times satisfying the minimum dwell - time constraint and we obtain the trajectories depicted in figure [ fig : mindt_switched_stabz ] where we can see observe the stabilizing effect of the controller .
let us consider the system with the matrices @xmath241 as in the previous examples , both subsystems are unstable and we seek for a controller of the form that makes the closed - loop system positive and stable under range dwell - time @xmath242 and @xmath243 .
we consider theorem [ th : range_switched_2 ] and solve the conditions using the sos relaxation with @xmath237 ( degree of @xmath238 and @xmath239 is one while the additional polynomials are of degree 2 ) . the semidefinite program to solve has 229 primal variables and 84 dual variables .
solving the conditions takes 1.5875 seconds and returns the controller matrices @xmath244 for simulation purposes , we generate a random sequence of dwell - times satisfying the range dwell - time constraint and we obtain the trajectories depicted in figure [ fig : rangedt_switched_stabz ] where we can see observe the stabilizing effect of the controller .
several stability and stabilization conditions for linear positive impulsive systems and linear positive switched systems have been obtained using the notion of dwell - times .
discrete - time stability conditions have been obtained and reformulated using a lifting approach in order to get conditions that can be applied for design purposes . because of their infinite - dimensional nature , several relaxation approaches have been proposed to make the lifted conditions finite - dimensional .
several examples have been given for illustration .
95 [ 1]#1 [ 1]`#1 ` urlstyle [ 1]doi : # 1        c.  briat , e.  a. yavuz , h.  hjalmarsson , k .- h .
johansson , u.  t. jnsson , g.  karlsson , and h.  sandberg .
the conservation of information , towards an axiomatized modular modeling approach to congestion control .
_ ieee transactions on networking _ , 23(3 ) , 2015 .
c.  briat and m.  khammash .
computer control of gene expression : robust setpoint tracking of protein mean and variance using integral feedback . in _
51st ieee conference on decision and control _ , pages 35823588 , maui , hawaii , usa , 2012 .
v.  jonsson , a.  rantzer , and r.  m. murray . a scalable formulation for engineering combination therapies for evolutionary dynamics of disease . in _ american control conference
_ , pages 27712778 , portland , usa , 2014
.    w.  m. haddad and v.  chellaboina .
stability and dissipativity theory for nonnegative dynamic systems : a unified analysis framework for biological and physiological systems .
_ nonlinear analysis : real world applications _ , 6:0 3565 , 2005 .      c.  briat .
robust stability and stabilization of uncertain linear positive systems via integral linear constraints - @xmath245- and @xmath246-gains characterizations .
_ international journal of robust and nonlinear control _ , 23(17):0 19321954 , 2013 .        c.  briat
robust stability analysis of uncertain linear positive systems via integral linear constraints - @xmath245- and @xmath246-gains characterizations . in _
50th ieee conference on decision and control _ , pages 31223129 , orlando ,
florida , usa , 2011 .
y.  ebihara , d.  peaucelle , and d.  arzelier .
@xmath245 gain analysis of linear positive systems and its applications . in _
50th conference on decision and control , orlando , florida , usa _ , pages 40294034 , 2011 .
m.  ait rami .
stability analysis and synthesis for linear positive systems with time - varying delays . in _ positive systems - proceedings of the 3rd multidisciplinary international symposium on positive systems : theory and applications ( posta 2009 ) _ , pages 205216 .
springer - verlag berlin heidelberg , 2009 .
y.  ebihara , d.  peaucelle , and d.  arzelier .
optimal @xmath249-controller synthesis for positive systems and its robustness properties . in _ american control conference _ , pages 59925997 , montreal , canada , 2012 .    c.  briat and m.  khammash .
interval peak - to - peak observers for continuous- and discrete - time systems with persistent inputs and delays . _ to appear in automatica _ , 2016 .
url http://arxiv.org/abs/1511.07888 .
a.  zappavigna , p.  colaneri , j.  c. geromel , and r.  h. middleton .
stabilization of continuous - time switched linear positive systems . in
_ american control conference _ , pages 32753280 , baltimore , maryland , usa , 2010 .
a.  zappavigna , p.  colaneri , j.  c. geromel , and r.  shorten .
dwell time analysis for continuous - time switched linear positive systems . in _
american control conference _ , pages 62566261 , baltimore , maryland , usa , 2010 .
f.  mazenc and m.  malisoff .
stability analysis for time - varying systems with delay using linear lyapunov functionals and a positive systems approach .
_ ieee transactions on automatic control _ , 61(3):0 771776 , 2016 .
n.  sivashankar and p.  p. khargonekar .
characterization of the @xmath251-induced norm for linear systems with jumps with applications to sampled - data systems .
_ siam journal on control and optimization _ , 32(4):0 11281150 , 1994 .
g.  chesi , p.  colaneri , j.  c. geromel , r.  middleton , and r.  shorten . a nonconservative lmi condition for stability of switched systems with guaranteed dwell - time .
_ ieee transactions on automatic control _ , 57(5):0 12971302 , 2012 .          c.  briat and a.  seuret
on the necessity of looped - functionals arising in the analysis of pseudo - periodic , sampled - data and hybrid systems .
_ international journal of control _ , 88((12):0 25632569 , 2015 .    c.  briat
theoretical and numerical comparisons of looped functionals and clock - dependent lyapunov functions - the case of periodic and pseudo - periodic systems with impulses . _ international journal of robust and nonlinear control _ , 26:0 22322255 , 2016 .          c.  briat .
convex lifted conditions for robust @xmath252-stability analysis and @xmath252-stabilization of linear discrete - time switched systems with minimum dwell - time constraint .
_ automatica _ , 50(3):0 976983 , 2014 .
a.  papachristodoulou , j.  anderson , g.  valmorbida , s.  prajna , p.  seiler , and p.  a. parrilo .
_ sostools : sum of squares optimization toolbox for matlab v3.00 _ , 2013 .
url http://www.cds.caltech.edu/sostools/sostools.pdf .
r.  kamyar , c.  murti , and m.  m. peet .
constructing piecewise - polynomial lyapunov functions for local stability of nonlinear systems using handelman s theorem . in _
53rd ieee conference on decision and control _
, pages 54815487 , los angeles , usa , 2014 .
g.  plya .
ber positive darstellung von polynomen vierteljschr , naturforsch . ges .
zrich , vol .
73 , pp . 141145 , 1928 .
volume  2 of _ collected papers ( ed .
r. p. boas ) _ , pages 309313 .
mit press , cambridge , ma , 1974 .              t.  tanaka and c.  langbort .
emma for internally positive systems and a tractable class of distributed h - infinity control problems . in _
american control conference _
, pages 62386243 , baltimore , maryland , usa , 2010 .    t.  tanaka and c.  langbort .
symmetric formulation of the s - procedure , kalman - yakubovich - popov lemma and their exact losslessness conditions .
_ ieee transactions on automatic control _ , 58(6):0 14861496 , 2013