Next Article in Journal
Quantum Circuit Template Matching Optimization Method for Constrained Connectivity
Next Article in Special Issue
Model-Free Sliding Mode Enhanced Proportional, Integral, and Derivative (SMPID) Control
Previous Article in Journal
Solving Integral Equations via Fixed Point Results Involving Rational-Type Inequalities
Previous Article in Special Issue
A New Reverse Extended Hardy–Hilbert’s Inequality with Two Partial Sums and Parameters
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Robust Stability of Switched-Interval Positive Linear Systems with All Modes Unstable Using the Φ-Dependent Dwell Time Technique

School of Mathematics and Computer Science, Shanxi Normal University, Taiyuan 030032, China
*
Author to whom correspondence should be addressed.
Submission received: 25 May 2023 / Revised: 8 July 2023 / Accepted: 10 July 2023 / Published: 13 July 2023
(This article belongs to the Special Issue Advances in Analysis and Control of Systems with Uncertainties II)

Abstract

:
In this study, some stability and robust stability conditions for switched positive linear systems in which all subsystems are unstable in continuous time and discrete time were obtained using the Φ-dependent dwell time technique and the discretized co-positive Lyapunov functions approach, respectively. The co-positive Lyapunov functions constructed in this study are functions of time during the dwell time, and after that, they are independent of time. In addition, the above method was applied to switched-interval positive systems, and corresponding conclusions are presented. The Φ-dependent dwell time technique used in this paper is more effective than the dwell time and mode-dependent dwell time used in other studies. The results are verified with an illustrative example.

1. Introduction

Many dynamic systems in the real world are restricted to orthogonal variables, and such systems are usually referred to as positive systems in the literature. Switched positive linear systems (SPLSs), a class of positive systems, consist of a series of positive linear systems and switching signals that control the switching between them. SPLSs have received extensive attention due to their wide applications in control fields, such as mobile robot navigation [1], load forecasting [2], and fuel balancing [3]. In SPLSs, stability is the most basic property of the system. Some important results regarding SPLSs have been obtained in the literature [4,5,6,7,8,9,10,11], especially with respect to stability analysis. The linear co-positive Lyapunov function (LCLF) approach is a very effective approach in the stability analysis of SPLSs because it holds that the state of the systems is positive under natural constraints. The sufficient and necessary conditions for the existence of a common LCLF for SPLSs are mentioned in [12]. The stability of discrete-time SPLSs was analyzed in [13] by using the switched LCLF approach. The multiple LCLF for a given SPLS was established for the first time in [14]. Based on this, the stability analysis problem for SPLSs with average-dwell-time (ADT) switching was investigated and sufficient conditions for both the continuous-time and the discrete-time cases were given in [14].
Most practical switched systems have both stable and unstable subsystems due to disturbances, failures, and possibl changing environmental factors [15]. The problem of exponential stability for SPLSs consisting of both stable and unstable subsystems was investigated in [16]. Sufficient stability conditions were proposed for discrete-time switched delay positive systems with stable and unstable subsystems in [17]. The above conclusions all concern switched systems with at least one stable subsystem. However, these conclusions are generally not valid when all subsystems are unstable. As is well known, even if all subsystems are unstable, one may carefully switch between them to make the total system asymptotically stable. Using a discretized Lyapunov function approach, a sufficient condition ensuring the asymptotic stability of switched continuous-time systems with all modes unstable is proposed in [18]. The stability of switched linear systems with dwell time was studied by constructing a discretized Lyapunov function in [19]. Some conclusions have been reached [20,21] about SPLSs with uncertainty. In [20], the constructed Lyapunov functions were time-varying during the dwell time and time-invariant afterwards. There are currently few research conclusions about the stability of SPLSs with all modes unstable and uncertainty, which are significant for our research in this paper.
The main objective of this study was to establish the stability and robust stability criteria for a system where all subsystems are unstable in continuous-time and discrete-time cases. The research results in this paper include the following points. First of all, a group of switching signals that make SPLSs with all modes unstable asymptotically stable was found by using discretized co-positive Lyapunov functions and the Φ -dependent dwell time technique. Then, we extended the above method and conclusion to SPLSs with interval uncertainty. Finally, an illustrative example is given to verify that the Φ -dependent dwell time technique is more effective than using mode-dependent dwell time [22] or dwell time.
The structure of this paper is as follows. Section 2 gives a description of the system and the necessary definitions and lemmas. Section 3 presents the main conclusions with their proofs. An illustrative example is given in Section 4. Section 5 provides the conclusion of this paper.
The notations used in this paper are shown in Table 1.

2. Problem Formulation and Preliminaries

Consider the continuous-time switched linear system
x ˙ ( s ) = G ρ ( s ) x ( s ) , s R , s s 0 ,
and the discrete-time switched linear system
x ( s + 1 ) = G ρ ( s ) x ( s ) , s N , s s 0 ,
where x ( s ) R n is the state of the system, and ρ ( s ) is called a switching signal. ρ ( s ) is a piece-wise constant right-continuous function: [ 0 , ) M = { 1 , 2 , , m } , where m N + is the number of subsystems. G u R n × n and u M are known constant matrices of appropriate dimensions.
All subsystems studied in this paper are unstable. The switching instants have the following relationship: 0 s 0 < s 1 < < s r < s r + 1 < , where s 0 represents the initial time of system operation and s r represents the rth switching instant. In addition, we assume that x ( s r + ) = x ( s r ) , r = 0 , 1 , 2 , . When s [ s r , s r + 1 ) , we say the ρ ( s r ) th subsystem of the switched system is active. The length of time between adjacent switching instants is called the dwell time τ r = s r + 1 s r , r = 0 , 1 , 2 , . Let K = { 1 , 2 , , k } where k N and k m . Define the surjection operator: Φ : M K . Set Φ i = { u M Φ ( u ) = i } [23]. This work supposes that each family of Φ i subsystems has a dwell time, denoted as τ Φ i , r , r = 0 , 1 , 2 , , i K and called the Φ -dependent dwell time. If τ r is too large, then the total system will be unstable due to it running unstable subsystems for a long time. If τ r is too small, then overly fast switching also makes the system unstable. Thus, we limit the Φ -dependent dwell time to a range that ensures the asymptotic stability of the system; namely, τ Φ i , r [ τ Φ i , min , 2 τ a Φ i , max τ Φ i , min ] , i K , r = 0 , 1 , 2 , , where τ Φ i , min represents the minimum dwell time of a family of Φ i subsystems, τ a Φ i , max represents the maximum average dwell time of a family of Φ i subsystems, and 0 < τ Φ i , min τ a Φ i , max . D [ τ Φ i , min , 2 τ a Φ i , max τ Φ i , min ] is called the switching strategy set of the Φ -dependent dwell time.
The aim of this study was to analyze the problem of the stability and robust stability of SPLSs where all subsystems are unstable. Before that, we first provide some definitions and lemmas to be used.
Definition 1 
([24]). Systems (1) and (2) are said to be positive if x ( s ) 0 , s > s 0 , ∀ ρ ( s ) , x ( s 0 ) 0 .
Definition 2 
([12]). A matrix G is said to be a Metzler matrix if its non-diagonal elements are positive or zero.
Consider the continuous-time system
x ˙ ( s ) = G x ( s ) , s R , s s 0 ,
and the discrete-time system
x ( s + 1 ) = G x ( s ) , s N , s s 0 .
We note the following lemmas.
Lemma 1 
([19]). System (3) is positive ⟺ G is a Metzler matrix, and System (4) is positive ⟺ G 0 .
From Lemma 1, we can deduce that System (1) is positive with regard to ρ ( s ) G u R n × n and u M are Metzler matrices, and System (2) is positive with regard to ρ ( s ) G u 0 , u M .
Lemma 2 
([19]). Let Systems (3) and (4) be positive; then, Systems (3) and (4) is asymptotically stable ⟺ there exists a vector ι 0 such that G T ι 0 ( ( G I ) T ι 0 ).
The function F ( s ) = x T ( s ) ι is said to be a linear co-positive Lyapunov function for the Systems (3) and (4) if there exists a vector ι 0 such that G T ι 0 ( ( G I ) T ι 0 ).

3. Main Results

This section presents sufficient conditions for the stability and robust stability of SPLSs in which all subsystems are unstable in the continuous-time and discrete-time cases.

3.1. Continuous-Time Case

In this section, sufficient conditions for the stability and robust stability of continuous SPLSs and their proofs are given.
Theorem 1. 
Consider SPLS (1). Given scalars λ i > 0 , 0 < μ i < 1 , i K , 0 < τ Φ i , min τ a Φ i , max , if there exists a set of vectors ι u , h 0 , h = 0 , 1 , 2 , , H , u M such that h = 0 , 1 , 2 , , H , u , v M , i K ,
Λ u , h T + ι u , h T G u λ i ι u , h T 0 ,
Λ u , h T + ι u , h + 1 T G u λ i ι u , h + 1 T 0 ,
ι u , H T G u λ i ι u , H T 0 ,
ι v , 0 μ i ι u , H 0 ,
ln μ i + λ i τ a Φ i , max < 0 ,
where Λ u , h = H ( ι u , h + 1 ι u , h ) τ Φ i , min , Φ ( u ) = i , then System (1) is globally asymptotically stable under any switching law ρ ( s ) D [ τ Φ i , min , 2 τ a Φ i , max τ Φ i , min ] .
Analysis . The Lyapunov functions constructed in most previous articles that were continuous in the dwell-time interval are not applicable when all subsystems are unstable. Using another method, we consider the construction of a discretized Lyapunov function to break the commutativity of Lyapunov functions in adjacent dwell-time intervals. In this way, systems for which all subsystems are unstable can be stabilized by proper switching.
Proof. 
Step 1: Prove that System (1) is stable.
For the convenience of narration, assume that ρ ( s r ) = u , ρ ( s r + 1 ) = v . Divide the interval [ s r , s r + τ Φ i , min ) into H equal parts, with each interval represented as I r , h u = [ s r + ω h u , s r + ω h + 1 u ) , h = 0 , 1 , , H 1 . The length of the equal division is l u = τ Φ i , min H , where ω h u = h · l u , h = 0 , 1 , , H .
We use the linear interpolation formula to construct linear functions for each segment of the minimum dwell-time interval. When s I r , h u , let ι u ( s ) = ( 1 α ) ι u , h + α ι u , h + 1 , where α = s s r ω h u l u . Substituting the expressions of α and l u into the above equation yields ι ˙ u ( s ) = Λ u , h . When α = 0 , we denote ι u ( s r + ω h u ) as ι u , h .
When s [ s r + τ Φ i , min , s r + 1 ) , let ι u ( s ) remain the value of the left endpoint during this period; i.e., ι u ( s ) = ι u , H . Obviously, ι ˙ u ( s ) = 0 .
Since the system we are studying is positive, we can construct the following multiple co-positive Lyapunov functions:
F u ( s ) = ι u T ( s ) x ( s ) , u M .
Among them, the vector ι u ( s ) is as defined above.
When s I r , h u ,
F ˙ u ( s ) = Λ u , h T x ( s ) + ι u T ( s ) G u x ( s ) = Λ u , h T x ( s ) + [ ( 1 α ) ι u , h T + α ι u , h + 1 T ] G u x ( s ) = ( 1 α ) [ Λ u , h T + ι u , h T G u ] x ( s ) + α [ Λ u , h T + ι u , h + 1 T G u ] x ( s ) .
Combining Equations (5) and (6), it can be concluded that F ˙ u ( s ) < λ i F u ( s ) .
When s [ s r + τ Φ i , min , s r + 1 ) , F ˙ u ( s ) = ι u , H T G u x ( s ) . Combining this with Equation (7), it can be concluded that F ˙ u ( s ) < λ i F u ( s ) .
From the above analysis, it can be seen that
F ˙ u ( s ) < λ i F u ( s ) , s [ s r , s r + 1 ) .
According to Equation (8) and the assumption that the system state does not jump during switching, it can be concluded that
F v ( s r + 1 + ) μ i F u ( s r + 1 ) , u v , u , v M .
Combining Equations (12) and (13), we have
F ρ ( s r ) ( s ) < e λ Φ ( ρ ( s r ) ) ( s s r ) F ρ ( s r ) ( s r ) e λ Φ ( ρ ( s r ) ) ( s s r ) μ Φ ( ρ ( s r 1 ) ) F ρ ( s r 1 ) ( s r ) < e λ Φ ( ρ ( s r ) ) ( s s r ) μ Φ ( ρ ( s r 1 ) ) e λ Φ ( ρ ( s r 1 ) ) ( s r s r 1 ) F ρ ( s r 1 ) ( s r 1 ) < e λ Φ ( ρ ( s r ) ) ( s s r ) μ Φ ( ρ ( s r 1 ) ) μ Φ ( ρ ( s r 2 ) ) μ Φ ( ρ ( s 0 ) ) e λ Φ ( ρ ( s r 1 ) ) ( s r s r 1 ) e λ Φ ( ρ ( s 0 ) ) ( s 1 s 0 ) F ρ ( s 0 ) ( s 0 ) = e λ Φ ( ρ ( s r ) ) ( s s r ) i = 1 k μ i n i e i = 1 k λ i T i F ρ ( s 0 ) ( s 0 ) e λ Φ ( ρ ( s r ) ) · 2 τ a Φ i , m a x i = 1 k μ i n i e λ i n i τ a Φ i , max F ρ ( s 0 ) ( s 0 ) ,
where n i denotes the switching numbers of the i th family subsystems over the interval [ 0 , s ] , and T i denotes the total activated time of the i th family subsystems over the interval [ 0 , s ] . Let τ = max { τ a Φ i , m a x } , λ = max { λ i } , i K . From Equation (9), we can derive μ i e λ i τ a Φ i , max < 1 ; namely, μ i n i e λ i n i τ a Φ i , max < 1 , i K . Substituting into Equation (14), we obtain:
F ρ ( s r ) ( s ) < e 2 λ τ F ρ ( s 0 ) ( s 0 ) .
As we all know, the following inequalities are true: F ρ ( s 0 ) ( s 0 ) δ ¯ n x ( s 0 ) 2 , δ ̲ x ( s ) 2 F ρ ( s r ) ( s ) , where δ ¯ = max u M , h = 0 , 1 , , H { δ ¯ ( ι u , h ) } , δ ̲ = min u M , h = 0 , 1 , , H { δ ̲ ( ι u , h ) } . Substituting into Equation (15), we obtain:
x ( s ) 2 Θ x ( s 0 ) 2 ,
where Θ = δ ¯ n δ ̲ e 2 λ τ . Then, ϵ > 0 , we can choose x ( s 0 ) 2 < ξ ( ϵ ) = Θ 1 ϵ . This results in x ( s ) 2 < ϵ , and the stability of System (1) can be obtained.
Step 2: Prove that System (1) is globally asymptotically stable.
Combining Equations (12) and (13), it can be seen that F v ( s r + 1 ) < γ F u ( s r ) , 0 < γ < 1 , where γ = max ( μ i e 2 λ i τ a Φ i , max ) . Recursively, the following relationship can be obtained: F ρ ( s r ) ( s r ) < γ r F ρ ( s 0 ) ( s 0 ) . We can derive lim r F ρ ( s r ) ( s r ) = 0 , which implies lim r x ( s r ) = 0 .
Below, we use the proof by contradiction to obtain the desired conclusion. Assume the existence of x ( s f ) such that lim f x ( s f ) = c , where c is a positive constant. From the definition of the limit, it can be seen that there exists q > 0 such that x ( s f ) 2 > c a whenever f > q , where a ( < c ) is any positive constant.
From the stability of System (1), it can be seen that there exists ζ ( p ) > 0 such that x ( s ) 2 < p whenever x ( s 0 ) < ζ , where p is positive. Let p = c a ; then, x ( s ) 2 < c a . This contradicts x ( s f ) 2   > c a , so the assumption is not valid, Then, we have lim s x ( s ) = 0 .
Hence, SPLS (1) is globally asymptotically stable with regard to ρ ( s ) D [ τ Φ i , min , 2 τ a Φ i , max τ Φ i , min ] . □
Remark 1. 
The co-positive Lyapunov functions in Theorem 1 are time-varying during [ s r , s r + τ Φ i , min ) and time-invariant during [ s r + τ Φ i , min , s r + 1 ) , where r = 0 , 1 , . The discretized co-positive Lyapunov function divides the dwell-time interval into a finite number of small regions, and the vector function ι u ( s ) varies linearly with respect to s in each small region.
Remark 2. 
The Lyapunov functions F u ( s ) in Theorem 1 are allowed to increase while the unstable subsystem is working, and the increase rate is limited to F ˙ u ( s ) < λ i F u ( s ) . Moreover, this increment can be compensated for by the switching behavior; that is, the Lyapunov functions do not increase at the switching instants.
Remark 3. 
The switching strategy set can be obtained by the following procedure. Firstly, appropriate λ i , μ i , H , and τ Φ i , min are given according to the constraints of each parameter. Then, the values of vectors ι u , h are calculated by linear programming. We can then obtain the switching strategy set. Later, we will visually see through examples that the Φ-dependent dwell time technique is more flexible than the mode-dependent dwell time and the dwell time.
Remark 4. 
Unlike most stability results based on various dwell time strategies or their extended forms, Theorem 1 applies to the situation where all subsystems are unstable. The theorem is based on the fact that the energy increment caused by unstable subsystem may be compensated for by the suitable “stable” switching. This has two limitations: the computational complexity caused by group diversity and the blindness in determining the minimum dwell time. The conditions of the theorem are to some extent harsh, but it provides an effective switching design scheme for related research.
Due to the fact that real-world systems are always influenced by various factors, the following investigates the robust stability of continuous switched-interval positive systems where all subsystems are unstable.
Consider the continuous switched-interval positive system
x ˙ ( s ) = G ρ ( s ) x ( s ) ,
where G u [ G ̲ u , G ¯ u ] , u M , G ̲ u , and G ¯ u , u M are constant matrices with suitable dimensions, representing interval uncertainty. In addition, we assume that G ̲ u are Metzler matrices.
Theorem 2. 
Consider SPLS (17). Given scalars λ i > 0 , 0 < μ i < 1 , i K , and 0 < τ Φ i , min τ a Φ i , max , if there exists a set of vectors ι u , h 0 , h = 0 , 1 , 2 , , H , u M such that h = 0 , 1 , 2 , , H , u , v M , i K ,
Λ u , h T + ι u , h T G ¯ u λ i ι u , h T 0 ,
Λ u , h T + ι u , h + 1 T G ¯ u λ i ι u , h + 1 T 0 ,
ι u , H T G ¯ u λ i ι u , H T 0 ,
ι v , 0 μ i ι u , H 0 ,
ln μ i + λ i τ a Φ i , max < 0 ,
where Λ u , h = H ( ι u , h + 1 ι u , h ) τ Φ i , min , Φ ( u ) = i , then System (17) is globally asymptotically stable under any switching law ρ ( s ) D [ τ Φ i , min , 2 τ a Φ i , max τ Φ i , min ] .
Proof. 
Under the assumption that G ̲ u are Metzler matrices, according to Definition 2, it can be concluded that G u are also Metzler matrices. We still use the discretized co-positive Lyapunov function (Equation (10)).
When s I r , h u ,
F ˙ u ( s ) Λ u , h T x ( s ) + ι u T ( s ) G ¯ u x ( s ) = ( 1 α ) [ Λ u , h T + ι u , h T G ¯ u ] x ( s ) + α [ Λ u , h T + ι u , h + 1 T G ¯ u ] x ( s ) .
Combining Equations (18) and (19), it can be concluded that F ˙ u ( s ) < λ i F u ( s ) . Similarly, when s [ s r + τ Φ i , min , s r + 1 ) , there is the same formula. In summary, it can be concluded that F ˙ u ( s ) < λ i F u ( s ) , s [ s r , s r + 1 ) . The subsequent proof is similar to Theorem 1. Therefore, there will be no further explanation. □

3.2. Discrete-Time Case

This section provides the sufficient conditions for the stability and robust stability of discrete-time SPLSs.
Theorem 3. 
Consider SPLS (2). Given scalars λ i > 0 , 0 < μ i < 1 , i K , and 0 < τ Φ i , min τ a Φ i , max , if there exists a set of vectors ι u , h 0 , h = 0 , 1 , 2 , , H , u M such that h = 0 , 1 , 2 , , H , u , v M , i K ,
Λ u , h T G u + ι u , h T ( G u I ) λ i ι u , h T 0 ,
Λ u , h T G u + ι u , h + 1 T ( G u I ) λ i ι u , h + 1 T 0 ,
ι u , H T ( G u I ) λ i ι u , H T 0 ,
ι v , 0 μ i ι u , H 0 ,
ln μ i + τ a Φ i , max ln ( 1 + λ i ) < 0 ,
where Λ u , h = H ( ι u , h + 1 ι u , h ) τ Φ i , min , Φ ( u ) = i , then System (2) is globally asymptotically stable under any switching law ρ ( s ) D [ τ Φ i , min , 2 τ a Φ i , max τ Φ i , min ] .
Proof. 
We still use the linear interpolation formula to construct linear functions for each segment of the minimum dwell-time interval I r , h u = [ s r + ω h u , s r + ω h + 1 u ) , h = 0 , 1 , , H 1 , as shown in Theorem 1. We can easily obtain ι u ( s + 1 ) = Λ u , h + ι u ( s ) . Then,
Δ F u ( s ) = F u ( s + 1 ) F u ( s ) = ( ι u T ( s ) + Λ u , h T ) G u x ( s ) ι u T ( s ) x ( s ) = [ ( 1 α ) ι u , h T + α ι u , h + 1 T + Λ u , h T ] G u x ( s ) [ ( 1 α ) ι u , h T + α ι u , h + 1 T ] x ( s ) = ( 1 α ) [ ι u , h T ( G u I ) + Λ u , h T G u ] x ( s ) + α [ ι u , h + 1 T ( G u I ) + Λ u , h T G u ] x ( s ) .
Combining Equations (24) and (25), it can be concluded that F u ( s + 1 ) < ( 1 + λ i ) F u ( s ) .
Similarly, when s [ s r + τ Φ i , min , s r + 1 ) , take ι u ( s ) = ι u , H . Using Equation (26), we can derive
Δ F u ( s ) = ι u , H T G u x ( s ) ι u , H T x ( s ) = ι u , H T ( G u I ) x ( s ) < λ i F u ( s ) .
From the above analysis, it can be seen that
F u ( s + 1 ) < ( 1 + λ i ) F u ( s ) , s [ s r , s r + 1 ) .
According to Equation (27) and the assumption that the system state does not jump during switching, it can be concluded that
F v ( s r + 1 + ) μ i F u ( s r + 1 ) , u v , u , v M .
From Equations (31) and (32), we can obtain
F ρ ( s r ) ( s ) < ( 1 + λ Φ ( ρ ( s r ) ) ) s s r F ρ ( s r ) ( s r ) ( 1 + λ Φ ( ρ ( s r ) ) ) s s r μ Φ ( ρ ( s r 1 ) ) F ρ ( s r 1 ) ( s r ) < ( 1 + λ Φ ( ρ ( s r ) ) ) s s r μ Φ ( ρ ( s r 1 ) ) ( 1 + λ Φ ( ρ ( s r 1 ) ) ) s r s r 1 F ρ ( s r 1 ) ( s r 1 ) < ( 1 + λ Φ ( ρ ( s r ) ) ) s s r μ Φ ( ρ ( s r 1 ) ) μ Φ ( ρ ( s r 2 ) ) μ Φ ( ρ ( s 0 ) ) ( 1 + λ Φ ( ρ ( s r 1 ) ) ) s r s r 1 ( 1 + λ Φ ( ρ ( s 0 ) ) ) s 1 s 0 F ρ ( s 0 ) ( s 0 ) = ( 1 + λ Φ ( ρ ( s r ) ) ) s s r i = 1 k μ i n i ( 1 + λ i ) T i F ρ ( s 0 ) ( s 0 ) ( 1 + λ Φ ( ρ ( s r ) ) ) 2 τ a Φ i , max i = 1 k μ i n i ( 1 + λ i ) n i τ a Φ i , max F ρ ( s 0 ) ( s 0 ) ,
where n i denotes the switching numbers of the i th family subsystems over the interval [ 0 , s ] , and T i denotes the total activated time of the i th family subsystems over the interval [ 0 , s ] . Let Ψ = max { ( 1 + λ i ) 2 τ a Φ i , max } . From Equation (28), we can derive μ i ( 1 + λ i ) τ a Φ i , max < 1 . Substituting into Equation (33), we obtain:
F ρ ( s r ) ( s ) < Ψ F ρ ( s 0 ) ( s 0 ) .
As we all know, the following inequalities are true: F ρ ( s 0 ) ( s 0 ) δ ¯ n x ( s 0 ) 2 , δ ̲ x ( s ) 2 F ρ ( s r ) ( s ) , where δ ¯ = max u M , h = 0 , 1 , , H { δ ¯ ( ι u , h ) } , δ ̲ = min u M , h = 0 , 1 , , H { δ ̲ ( ι u , h ) } . Substituting into Equation (34), we obtain:
x ( s ) 2 Ω x ( s 0 ) 2 ,
where Ω = δ ¯ n δ ̲ Ψ . Then, ϵ > 0 , and we can choose x ( s 0 ) 2 < ξ ( ϵ ) = Ω 1 ϵ . This results in x ( s ) 2 < ϵ , and the stability of System (2) can be obtained.
The rest of the proof of Theorem 3 is similar to Theorem 1 and is therefore omitted. □
Below, we provide sufficient conditions for the case of interval uncertainty.
Consider the discrete switched-interval positive system
x ( s + 1 ) = G ρ ( s ) x ( s ) ,
where G u [ G ̲ u , G ¯ u ] , u M , G ̲ u , and G ¯ u , u M are the same as Equation (17). In addition, we assume that G ̲ u 0 .
Theorem 4. 
Consider Equation (36). Given scalars λ i > 0 , 0 < μ i < 1 , i K , and 0 < τ Φ i , min τ a Φ i , max , if there exists a set of vectors ι u , h 0 , h = 0 , 1 , 2 , , H , u M such that h = 0 , 1 , 2 , , H , u , v M , i K ,
Λ u , h T G ¯ u + ι u , h T ( G ¯ u I ) λ i ι u , h T 0 ,
Λ u , h T G ¯ u + ι u , h + 1 T ( G ¯ u I ) λ i ι u , h + 1 T 0 ,
ι u , H T ( G ¯ u I ) λ i ι u , H T 0 ,
ι v , 0 μ i ι u , H 0 ,
ln μ i + τ a Φ i , max ln ( 1 + λ i ) < 0 ,
where Λ u , h = H ( ι u , h + 1 ι u , h ) τ Φ i , min and Φ ( u ) = i ; then, System (36) is globally asymptotically stable under any switching law ρ ( s ) D [ τ Φ i , min , 2 τ a Φ i , max τ Φ i , min ] .
The proof of this theorem is basically the same as the previous proof and is omitted here.

4. Illustrative Example

In this section, we illustrate through a specific example that the Φ -dependent dwell time technique is more flexible than mode-dependent dwell time and dwell time.
Consider an epidemiological model with n population groups [25]. Every ı th group is divided into infectives and susceptibles. Let I ı ( s ) and S ı ( s ) , respectively, denote the number of infectives and susceptibles at time s. Suppose the total number I ı ( s ) + S ı ( s ) = N ı is constant at any time s 0 . By taking x ı ( s ) = I ı ( s ) / N ı , we have, for ı = 1 , 2 , , n :
x ˙ ı ( s ) = ( 1 x ı ( s ) ) ȷ = 1 N β ı ȷ N ȷ N ı x ȷ ( s ) ( η ı + ϖ ı ) x ı ( s ) ,
where β ı ȷ > 0 and η ı > 0 are some known constants, and ϖ ı > 0 is the death rate in the ı group. Suppose there are m different therapies to fight the epidemic, and the infection rate η ρ ( s ) is not constant but depends, at each time s, on ρ ( s ) { 1 , 2 , , m } , which orchestrates the different therapies for the different population groups. By linearizing the system around the disease-free point x = 0 , one can get the following system:
x ˙ ( s ) = G ρ ( s ) x ( s ) ;
i.e., Equation (1) with the system’s matrices given as follows:
G 1 = 2.75 0.05 0.045 0.15 , G 2 = 0.1 0.035 0.03 1.75 , G 3 = 0.2 0.06 0.05 3 .
The eigenvalues of G 1 are a 1 = 2.7508 and a 2 = 0.1508 , the eigenvalues of G 2 are a 1 = 0.1006 and a 2 = 1.7506 , and the eigenvalues of G 3 are a 1 = 0.2009 and a 2 = 3.0009 . Obviously, all three subsystems are unstable.
For the convenience of calculation, we take H = 1 , x ( 0 ) = [ 3 , 5 ] T . D i represents the switching strategy set of the i th family subsystem, and τ i represents the dwell time of the i th family subsystem. Due to the number of subsystems being three, the set K can be divided into the following three situations. According to the choice of the different Φ , one can get different switching strategies. When K = { 1 } , then Φ 1 = { 1 , 2 , 3 } , which corresponds to the case with the dwell time technique. When K = { 1 , 2 } , there are three types of grouping: (1) Φ 1 = { 1 , 2 } , Φ 2 = { 3 } ; (2) Φ 1 = { 1 , 3 } , Φ 2 = { 2 } ; and (3) Φ 1 = { 1 } , Φ 2 = { 2 , 3 } . When K = { 1 , 2 , 3 } , then Φ 1 = { 1 } , Φ 2 = { 2 } , and Φ 3 = { 3 } , which corresponds to the mode-dependent dwell time technique. In each group, select appropriate λ i > 0 , 0 < μ i < 1 , and τ Φ i , min , and use linear programming to solve the vector ι u , h that satisfies Equations (5)–(8) in Theorem 1. When Equations (5)–(8) have a solution, by bringing the values of λ i and μ i into Equation (9), we can obtain τ a Φ i , max . The switching strategy set D i can be obtained from τ Φ i , min and τ a Φ i , max . Select the appropriate switching τ i in D i . Based on the selected τ i , the state response (Figure 1, Figure 2, Figure 3, Figure 4 and Figure 5) can be obtained for different groups Φ . The following table provides the parameter selection and switching strategy set corresponding to the dwell time, Φ -dependent dwell time, and mode-dependent dwell time techniques, respectively. As can be seen from the state response figures, all three techniques can stabilize the system with appropriate switching. It can be seen from Table 2 that the τ Φ i , min of the dwell tine technique is 0.5 , while the τ Φ i , min can be 0.4 in the three cases where K = { 1 , 2 } . Similarly, the Φ -dependent dwell time technique has a broader range of switching strategy sets than the mode-dependent one. By comparing the switching signal sets and state response graphs of the different techniques, it can be seen that the conclusions obtained with the Φ -dependent dwell time cover the previous dwell time and mode-dependent dwell time. The technique used in this paper can make the system stable in the shortest time with the appropriate switching.

5. Conclusions

In this paper, we studied the stability and robust stability of SPLSs in which all subsystems are unstable by means of Φ -dependent dwell-time switching. By using the discretized co-positive Lyapunov functions, the sufficient conditions for the stability of SPLSs were obtained in the form of linear matrix inequalities. The stability and robust stability of SPLSs in continuous-time and discrete-time cases were studied in this paper, respectively. At the end of the paper, an illustrative example showed that the Φ -dependent dwell time technique is more effective than the mode-dependent dwell time and the dwell time.

Author Contributions

Conceptualization, Q.Y.; Methodology, Q.Y.; Software, X.J.; Validation, X.J.; Formal analysis, Q.Y.; Investigation, Q.Y.; Resources, X.J.; Data curation, X.J.; Writing—original draft, Q.Y. and X.J.; Writing—review & editing, Q.Y.; Visualization, X.J.; Supervision, Q.Y.; Project administration, X.J.; Funding acquisition, Q.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Fundamental Research Program of Shanxi Province (202103021224249) and the Fund Program for the Scientific Activities of Selected Returned Overseas Professionals in Shanxi Province (20220023).

Data Availability Statement

The data supporting the reported results are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare that they do not have any commercial or associative interests that would represent a conflict of interests in connection with the work.

References

  1. Jin, J.; Kim, Y.; Wee, S.; Lee, D.; Gans, N. A stable switched-system approach to collision-free wheeled mobile robot navigation. J. Intell. Robot. Syst. Theory Appl. 2017, 86, 599–616. [Google Scholar] [CrossRef]
  2. Maia, C.A.; Goncalves, M.M. Application of switched adaptive system to load forecasting. Electr. Power Syst. Res. 2008, 78, 721–727. [Google Scholar] [CrossRef]
  3. Nilsson, P.; Ozay, N.; Topcu, U.; Murray, R.M. Temporal logic control of switched affine systems with an application in fuel balancing. In Proceedings of the 2012 American Control Conference (ACC), Montreal, QC, Canada, 27–29 June 2012; pp. 5302–5309. [Google Scholar]
  4. Fornasini, E.; Valcher, M.E. Asymptotic stability and stabilizability of special classes of discrete-time positive switched systems. Linear Algebra Appl. 2013, 438, 1814–1831. [Google Scholar] [CrossRef]
  5. Blanchini, F.; Colaneri, P.; Valcher, M.E. Co-positive Lyapunov functions for the stabilization of positive switched systems. IEEE Trans. Autom. Control 2012, 57, 3038–3050. [Google Scholar] [CrossRef]
  6. Fornasini, E.; Valcher, M.E. On the stability of continuous-time positive switched systems. In Proceedings of the American Control Conference, Baltimore, MD, USA, 30 June–2 July 2010; pp. 225–6230. [Google Scholar]
  7. Gurvits, L.; Shorten, R.; Mason, O. On the stability of switched positive linear systems. IEEE Trans. Autom. Control 2007, 52, 1099–1103. [Google Scholar] [CrossRef]
  8. Liu, Z.; Wang, Y.Z.; Zhao, G.D. Exponential stability for positive switched linear systems with delays and impulses. In Proceedings of the 33rd Chinese Control Conference, Nanjing, China, 28–30 July 2014; pp. 2469–2474. [Google Scholar]
  9. Ding, X.Y.; Shu, L.; Wang, Z.H. On stability for switched linear positive systems. Math. Comput. Model. 2011, 53, 1044–1055. [Google Scholar] [CrossRef]
  10. Zhao, X.D.; Zhang, L.X.; Shi, P. Stability of a class of switched positive linear time-delay systems. Int. J. Robust Nonlinear Control 2013, 23, 578–589. [Google Scholar] [CrossRef]
  11. Zhao, X.D.; Zhang, L.X.; Shi, P.; Karimi, H.R. Stabilization of a class of slowly switched positive linear systems: State-feedback control. In Proceedings of the American Control Conference (ACC), Montreal, QC, Canada, 27–29 June 2012; pp. 5282–5286. [Google Scholar]
  12. Mason, O.; Shorten, R. On linear copositive Lyapunov functions and the stability of switched positive linear systems. IEEE Trans. Autom. Control 2007, 52, 1346–1349. [Google Scholar] [CrossRef]
  13. Liu, X.W. Stability analysis of switched positive systems: A switched linear copositive Lyapunov function methed. IEEE Trans. Circuits Syst. II Express Briefs 2009, 56, 414–418. [Google Scholar]
  14. Zhao, Z.D.; Zhang, L.X.; Shi, P.; Liu, M. Stability of switched positive linear systems with average dwell time switching. Automatica 2012, 48, 1132–1137. [Google Scholar] [CrossRef]
  15. Yang, H.; Cocquempot, V.; Jiang, B. On stabilization of switched nonlinear systems with unstable modes. Syst. Control Lett. 2009, 58, 703–708. [Google Scholar] [CrossRef]
  16. Lian, J.; Liu, J. New results on stability of switched positive systems: An average dwell-time approach. IET Control Theory Appl. 2013, 7, 1651–1658. [Google Scholar] [CrossRef]
  17. Shi, R.Q.; Tan, X.M.; Zhao, X.D.; Zheng, X.L. Stability and l1-gain analysis for switched delay positive systems with stable and unstable subsystems. Circuits Syst. Signal Process. 2015, 34, 1683–1696. [Google Scholar] [CrossRef]
  18. Xiang, W.M.; Xiao, J. Stabilization of switched continuous-time systems with all modes unstable via dwell time switching. Automatica 2014, 50, 940–945. [Google Scholar] [CrossRef]
  19. Allerhand, L.; Shaked, U. Robust stability and stabilization of linear switched systems with dwell time. IEEE Trans. Autom. Control 2011, 56, 381–386. [Google Scholar] [CrossRef]
  20. Liu, J.J.; Zhang, K.J.; Wei, H.K. Robust stability of positive switched systems with dwell time. Int. J. Syst. Sci. 2016, 47, 2553–2562. [Google Scholar] [CrossRef]
  21. Zhang, J.F.; Han, Z.Z.; Wu, H.; Huang, J. Robust stabilization of discrete-time positive switched systems with uncertainties and average dwell time switching. Circuits Syst. Signal Process. 2014, 33, 71–95. [Google Scholar] [CrossRef]
  22. Feng, S.X.; Wang, J.; Zhao, J. Stability and robust stability of switched positive linear systems with all modes unstable. IEEE/CAA J. Autom. Sin. 2019, 6, 170–179. [Google Scholar] [CrossRef]
  23. Yu, Q.; Zhai, G.S. Stability analysis of switched systems under Φ-dependent average dwell time approach. IEEE Access 2020, 8, 30655–30663. [Google Scholar] [CrossRef]
  24. Zhang, J.F.; Zhang, W.; Cai, X.S.; Han, Z.Z. Stability and stabilization of positive switched systems under asynchronous switching. In Proceedings of the 2014 International Conference on Mechatronics and Control (ICMC), Jinzhou, China, 3–5 July 2014; pp. 347–352. [Google Scholar]
  25. Rami, M.A.; Bokharaie, V.S.; Mason, O.; Wirth, F.R. Stability criteria for sis epidemiological models under switching policies. Discret. Contin. Dyn. Syst.—Ser. B 2013, 19, 2865–2887. [Google Scholar]
Figure 1. The state response of the system with the dwell time technique.
Figure 1. The state response of the system with the dwell time technique.
Axioms 12 00686 g001
Figure 2. The state response of the system with group one of the Φ -dependent dwell time technique.
Figure 2. The state response of the system with group one of the Φ -dependent dwell time technique.
Axioms 12 00686 g002
Figure 3. The state response of the system with group two of the Φ -dependent dwell time technique.
Figure 3. The state response of the system with group two of the Φ -dependent dwell time technique.
Axioms 12 00686 g003
Figure 4. The state response of the system with group three of the Φ -dependent dwell time technique.
Figure 4. The state response of the system with group three of the Φ -dependent dwell time technique.
Axioms 12 00686 g004
Figure 5. The state response of the system with the mode-dependent dwell time technique.
Figure 5. The state response of the system with the mode-dependent dwell time technique.
Axioms 12 00686 g005
Table 1. The notations used in this paper.
Table 1. The notations used in this paper.
R The set of real numbers
R n The set of n-dimensional real vectors
R n × n The space of n × n real matrices
N ( N + )The set of nonnegative (positive) integers
G T The transpose of a matrix G
ι 0 ( ι 0 ) Each component of vector ι is nonnegative (positive)
δ ̲ ( ι ) ( δ ¯ ( ι ) ) The minimal (maximal) components of vector ι
G 0 ( G 0 ) Each component of matrix G is nonnegative (positive)
x 1 1-norm of x ( t ) ; i.e., x 1 = i = 1 n x i
x 2 Euclidean vector norm of x ( t ) ; i.e., x 2 = ( i = 1 n x i 2 ) 1 2
If and only if
Table 2. Comparison of the three methods.
Table 2. Comparison of the three methods.
Φ -Dependent Dwell Time
TechniqueDwell Time K = { 1 } K = { 1 , 2 } Mode-Dependent Dwell Time K = { 1 , 2 , 3 }
Φ Φ 1 = { 1 , 2 , 3 } Φ 1 = { 1 , 2 } Φ 1 = { 1 , 3 } Φ 1 = { 1 } Φ 1 = { 1 }
Φ 2 = { 3 } Φ 2 = { 2 } Φ 2 = { 2 , 3 } Φ 2 = { 2 }
Φ 3 = { 3 }
λ λ 1 = 0.5 λ 1 = 0.5 λ 1 = 0.5 λ 1 = 0.5 λ 1 = 0.5
λ 2 = 0.55 λ 2 = 0.57 λ 2 = 0.6 λ 2 = 0.6
λ 3 = 0.55
μ 1 = 0.75 μ 1 = 0.7 μ 1 = 0.7 μ 1 = 0.7 μ 1 = 0.75
μ μ 2 = 0.65 μ 2 = 0.6 μ 2 = 0.68 μ 2 = 0.7
μ 3 = 0.65
ι 1 , 0 0.0100 0.0116 0.0100 0.0124 0.0100 0.0124 0.0100 0.0127 0.0100 0.0119
ι 1 , 1 0.0177 0.0133 0.0234 0.0143 0.0255 0.0143 0.0226 0.0143 0.0214 0.0133
ι 2 , 0 0.0132 0.0100 0.0164 0.0100 0.0178 0.0100 0.0158 0.0100 0.0161 0.0100
ι 2 , 1 0.0157 0.0133 0.0195 0.0143 0.0211 0.0167 0.0189 0.0147 0.0195 0.0143
ι 3 , 0 0.0118 0.0100 0.0137 0.0100 0.0126 0.0100 0.0128 0.0100 0.0137 0.0100
ι 3 , 1 0.0133 0.0154 0.0154 0.0190 0.0143 0.0177 0.0147 0.0187 0.0154 0.0183
D 1 [ 0.5 , 0.55 ] [ 0.5 , 0.7 ] [ 0.5 , 0.7 ] [ 0.4 , 0.7 ] [ 0.4 , 0.55 ]
D 2 [ 0.4 , 0.75 ] [ 0.4 , 0.85 ] [ 0.45 , 0.6 ] [ 0.45 , 0.5 ]
D 3 [ 0.4 , 0.53 ]
τ 1 0.5 0.6 0.6 0.55 0.4
τ 2 0.5 0.6 0.45 0.5 0.45
τ 3 0.5 0.5 0.6 0.5 0.5
State responseFigure 1Figure 2Figure 3Figure 4Figure 5
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yu, Q.; Jiang, X. Robust Stability of Switched-Interval Positive Linear Systems with All Modes Unstable Using the Φ-Dependent Dwell Time Technique. Axioms 2023, 12, 686. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms12070686

AMA Style

Yu Q, Jiang X. Robust Stability of Switched-Interval Positive Linear Systems with All Modes Unstable Using the Φ-Dependent Dwell Time Technique. Axioms. 2023; 12(7):686. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms12070686

Chicago/Turabian Style

Yu, Qiang, and Xiujuan Jiang. 2023. "Robust Stability of Switched-Interval Positive Linear Systems with All Modes Unstable Using the Φ-Dependent Dwell Time Technique" Axioms 12, no. 7: 686. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms12070686

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop