# Stability Analysis for a Type of Multiswitching System with Parallel Structure

Authors: Zhang Yan; Liu Yongqiang; Liu Yang
Journal: Mathematical Problems in Engineering
 (2018)

---

## Abstract

This paper proposes a new type of multiswitching system and a subsystems-group as a basic switching unit that obeys the law. Unlike traditional switching systems, the system selects multiple subsystems instead of one on each time interval. Thus, a framework of parallel structure organizes the subsystems as a group. A multiswitched system is widely used in engineering for modelling and control; this system reflects the actual industrial dynamical process. Thus, the stability of the system is studied. Assuming that these continuous and discrete-time subsystems are Hurwitz and Schur stable, the subsystems-groups matrices commute each other based on the subsystems matrices pairwise commutative. Then, the multiswitched system is exponentially stable under arbitrary switching, and there exists a common Lyapunov function for these subsystems. The main result is extended to a parallel-like structure; therefore, some stability results are gained under some reasonable assumption. At last, a numeral example is given to illustrate the structure and the stability of this system.


---

## Body

## 1. Introduction

A switched system is a particular kind of hybrid dynamical system that consist of a finite number of subsystems and a switching signal that switches between these subsystems in an orderly manner. The set of subsystems can be continuous-time and discrete-time; the mathematical models are differential equations and difference equations. According to the dynamics behaviours, it can be classified as continuous-time switched system, discrete-time switched system, or mixed switched system (both them) [1], linear or nonlinear and so on [2]. Actually, many fields of complex industrial processes include typical switching systems, which reflect the hybrid dynamical characteristics that are widely used in many fields [3–8], such as in aircraft control, the automotive industry, communication systems, electrical engineering, chemical processes, and so on.Multiswitched systems are a novel type of system, which can be composed of a set of continuous-time subsystems-groups(1)x.t=Ecixt,ci=c1,…,cNc(2)Eci=Ac1⋯0⋯0⋮⋱⋮0Acp0⋮⋱⋮0⋯0⋯AcM,cp=c1,…,cMand a set of discrete-time subsystems-groups(3)xk+1=Edjxk,dj=d1,…,dNd(4)Edj=Ad1⋯0⋯0⋮⋱⋮0Adq0⋮⋱⋮0⋯0⋯AdM,dq=d1,…,dMwherex(t),x(k) are the states,t is the time scalar, and τ is positive integer of the sample period.Eci andEdj are nonnull constant matrices with appropriate dimensions, which are combined by some constant matricesAcp andAdq, respectively. Subsystems-groups (1) and (2) are composed of continuous-time subsystems x.(t)=Acpx(t) with cp=c1,…,cM and discrete-time subsystems x(k+1)=Adqx(k) with dq=d1,…,dM, respectively. The system switches between matrices Ec1,Ec2,…,EcNc,Ed1,Ed2,…,EdNd belonging to the setE;E≜Ec,Ed, Ec≜Ec1,Ec2,…,EcNc, and Ed≜Ed1,Ed2,…,EdNd. Ac1,Ac2,…,AcM belong to the set Ac≜Ac1,Ac2,…,AcM, Ad1,Ad2,…,AdM belong to the set Ad≜Ad1,Ad2,…,AdM, and A≜Ad,Ac. A and E are finite sets.The major difference between a multiswitched system and a traditional switched system is that the former selects any amount of subsystems, and the latter selects only one each time. There is a structure for organizing the subsystems as a subsystems-group on each interval that includes mutual independence between subsystems. Many different structures can be considered, such as parallel, tandem, and other mixed modes. Equations (2) and (4) show the parallel structure of the matrices. However, the traditional switching system has only one matrix without any structure.A multiswitched system is widely used in engineering cases, such as in the chilled water system in central air conditioning. The dynamical pumping process is a kind of typical multiswitch system, and each different pump as a subsystem with switching and variable frequency behaviours. For variable flow control, there are several different pumps working together as a subsystem-group at each time period. According to the water pipe network, the pumps are parallel in structure. Thus, this physical model illustrates the proposed system. In this paper, only parallel and parallel-like structures are researched.Modelling [9], performance analysis [10] and optimal control [11, 12] are the key issues of switched systems. Obviously, stability analysis is a very important research branch that has attracted the attention of researchers globally. Stability is the most basic property and the primary issue to be solved by the control system. The stability problems of switched systems are reduced to the following three basic issues [13]: finding conditions to guarantee asymptotical stability for any switching signal; identifying some switching signals for asymptotical stability; and constructing a switching signal to make the system asymptotically stable. Accordingly, some excellent theoretical methods are proposed for solving those problems, such as common Lyapunov functions [14, 15], multiple Lyapunov functions [16–18], dwell time and average dwell time [19, 20], and piecewise quadratic Lyapunov functions [16].It is worth noting that stability under arbitrary switching is a fundamental in the design and analysis of switched systems [21]. It is known that all the subsystems that are exponentially stable are not sufficient to guarantee stability under arbitrary switching, expect for some reasonable assumptions, such as the state matrices of subsystems commute pairwise [2, 14] (e.g., AiAj=AjAi, for all i,j∈I, i≠j; AciAdj=AdjAci, for all i,j∈I, c and d denote continuous-time and discrete-time subsystems, respectively); the state matrices of subsystems are symmetric [22] (e.g., Ai=AiT, for all i∈I); the state matrices of subsystems are normal [23] (e.g., AciTAci=AciAciT, for all i∈I and AdjTAdj=AdjAdjT, for all j∈I,c andd denote the same substances above). Thus, some researchers present a sufficient condition in terms of the Lie algebra. If this Lie algebra is solvable, the exponential stability of the switched system for arbitrary switching can be gained [20]. A more general result is that the semisimple subalgebraS is a compact Lie algebra by considered Levi decomposition, and the exponentially stability can be ensured [24]. Previous research shows that the pairwise commutation of the vector fields is a sufficient condition for the stability [25]. By using Hurwitz stable matrix pencils σα[A1,A2] and σα[A1,A2-1], a common quadratic Lyapunov function exists that can guarantee stability under arbitrary switching [26]. Based on previous research, the main emphasis is placed on the conditions of stability analysis and a common Lyapunov function is constructed under an assumption.In this paper, a type of linear multiswitched system with parallel structure is put forward for the first time. The composition of the system is expounded, and the structure and form of the subsystems-group are described. From the switching law and the physical meaning, we illustrate the difference between the system and the traditional switching system. A case of chilled water system pumps in central air conditioning shows the engineering significance of the system. Next, stability under arbitrary switching for this type of system is studied. The inference of properties of state matrices of the subsystems-group is derived from the subsystems. Then, based on the conditions of Hurwitz stable and Schur stable for continuous-time subsystems and discrete-time subsystems, respectively, the stability of subsystems-groups can be obtained. The matricesEc and Ed of the subsystems-groups commute pairwise can be ensured based on the assumptions of the state matrices Ac and Ad [2, 14]. Thus, the stability under arbitrary switching is guaranteed, and a common Lyapunov function is given for all subsystems-groups and subsystems. Finally, the result is extended to a parallel-like structure, which reduces conservativeness in stability analysis.The body of this paper is organized as follows: in Section2, the system description, the example of the multiswitching system, and the preliminaries are presented. In Section 3, the main stability results in the continuous-time system and the mixed system are expressed; two examples are also given to illustrate the results. In Section 4, the studies extended to parallel-like structure is presented. Section 5 concludes the paper.


## 2. System Description and Mathematical Preliminaries

### 2.1. System Description and Examples

Throughout, the following notation is adopted:R and C denote the fields of real and complex number, respectively; Rn denotes then-dimensional real Euclidean space; Rn×n denotes the space of n×n matrices with real entries;xi andxj denote theith andjth components of the vectorx in Rn, respectively; acpp′ and adqq′ denote the entry in the (p,p′) and (q,q′) position of the matricesA orE in Rn×n, respectively.There are three examples to show the multiswitched system and the parallel structure.Example 1.
Consider a setA=A1,A2,A3 of a switched system, which are the constant matrices of 3 subsystems, respectively. Assume in the classical switched system that there are 3 subsystems. The system switches between A1, A2, and A3. Assume in the multiswitched system that there are 7 subsystems-groups. The system switches betweenE,(5)E=E1,E2,E3,E4,E5,E6,E7=A100000000,0000A20000,00000000A3,A1000A20000,A10000000A3,0000A2000A3,A1000A2000A3Here,Nc=CcM1+⋯+CcMcM and Nd=CdM1+⋯+CdMdM.Remark 2.
A continuous-time subsystems-group composed of a subsystem or all the subsystems. Obviously, there are subsystems, but the value is null. It is same with a discrete-time subsystems-group.Remark 3.
Continuous-time subsystems cannot mix with discrete-time subsystems to be a subsystems-group. Continuous-time subsystems-groups and discrete-time subsystems-groups must be distinguished.Remark 4.
WhenEci is a singular matrix, zero rows or zero columns remain with the aim of uniform description. In Example 1, E1=A1, but the uniform description is kept to show E1=A1000A2000A3=A100000000, with A2=0, A3=0. Thus, allEci values look the same, ci=c1,…,cNc. The same is true of Edj.Remark 5.
Equations (2) and (4) are standard parallel structures. A parallel-like structure will be considered in Section 4, which has a coupling phenomenon based on a parallel structure between some subsystems. Assumptions are made for reducing conservativeness to obtain ideal stability.Example 6.
In a chilled water system of a central air conditioning system, there are three pumps driving the chilled water from the evaporator to the air conditioning unit (see Figure1). All the pumps can be switched ON/OFF and the variable frequency obeying a range of 50%-100% rated frequency. For setting an energy-saving control strategy with variable water volume technology, the whole volume cannot be less than 50% of the rated volume. Assuming in a time interval that the cold load changes, the switching, and pump water volume are shown in Figures 2 and 3. The parameters and symbols are shown in Table 2.
It is assumed that the cold load in a certain area requires air conditioning to meet the cold demand within 5 minutes to achieve a balanced state. Then, in the first 5 minutes, 3 pumps work at the rated frequency. After reaching the preset value, the cold load is stable with small fluctuations. In this period, pumpA andB exit; only pumpC works in the form of frequency conversion, maintaining the driving chilled water and transferring cooling capacity. In the last 3 minutes, due to the cold load increase, a pump is not enough to meet the demand, and the system puts pumpA into work with pumpC together with frequency conversion.
From the three time intervals ([0,5min], [5min,12min], and [12min,15min]), we find that there are three, one, and two subsystems working, respectively. According to the framework in this paper, the three pumps denote different continuous-time subsystems, and the work combinations of the pumps denote the different continuous-time subsystems-groups. Assuming the pumpA,B, andC denote the subsystem X1, X2, and X3, respectively; and the combinationA/B/C,C andA/C we just said denote the subsystems-group Y1, Y2, and Y3, respectively.Figure 1
Chilled water system of central air conditioning.Figure 2
Total flow volume.Figure 3
Water volumes of pumps.Example 7.
In the same system with Example6, but the pumps are divided into one variable frequency pump and two switchable fixed frequency pumps. Thus, the former is a continuous-time subsystem, and the latter are discrete-time subsystems. According to the energy-saving control strategy, assuming pumpC is always working with variable frequency, and pumpsA andB are switched ON or OFF with fixed frequency. Thus, the chilled water system is composed of continuous and discrete-time subsystems. The parameters and symbols are shown in Table 3.There are two simulation figures (Figures4 and 5) that illustrate the switching and working situations of the system.Figure 4
Total flow volume.Figure 5
Pump water volume.In the first 5 minutes, 3 pumps work together; this seems a subsystems-group with one continuous and two discrete-time subsystems. In the interval of from 5 min to 12 min, only pumpC works, which is the only one continuous-time subsystem that is a subsystems-group. In the last three minutes, pumpA works with pumpC; therefore, the subsystems-group has one continuous and one discrete-time subsystem.


### 2.2. Important Theories for Stability Analysis

We analyse the stability of the multiswitched system based on the structural property ofE-matrices, which cannot exist without the basic component,A-matrices. All the inferences of multiswitched system are derived from a classical switched system. Thus, some important theorems and propositions of the classical switched system should be quoted first.Theorem 8 (see [14]).
Consider a classical continuous-time switched systemx.=Acx(t) with Ac≜Ac1,Ac2,…,AcM, where the matrices Acp are asymptotically stable and commute pairwise. Then, we get the following:(1)
The system is exponentially stable under any arbitrary switching between the elements ofAc.(2)
There exists a common Lyapunov function for all the subsystems. For any positive symmetric definite matrixPc0, let Pc1,Pc2,…,PcM be the unique symmetric positive definite solutions to the Lyapunov equations(6)AcpTPcp+PcpAcp=-Pcp-1,cp=c1,c2,…,cM
The functionV(x)=xTPcMx is a common Lyapunov function for each of the individual system x.=Acx(t) and hence a Lyapunov function for the switching system.Remark 9.
All of the proposed descriptions and notations are identical to the given multisystem.Theorem 10 (see [1]).
Consider a classical mixed switched system composed of continuous-time subsystemsx.=Acx(t) with Ac≜Ac1,Ac2,…,AcM and discrete-time subsystems x(k+1)=Adx(k) with Ad≜Ad1,Ad2,…,AdM. Let Acp be Hurwitz stable and Adq be Schur stable to commute pairwise. Then, we get the following:(1)
The system is exponentially stable for any arbitrary switching between the elements ofA.(2)
There exists a common Lyapunov function for all the subsystems. For any positive definite matrixPd0, let Pd1,Pd2,…,PdM be the unique positive definite solutions to the Lyapunov equations(7)Ad1TPd1Ad1+Pd1=-Pd0(8)AdqTPdqAdq+Pdq=-Pdq-1
letPc1,Pc2,…,PcM be the unique positive definite solutions to the Lyapunov equations(9)AcpTPcp+PcpAcp=-Pcp-1(10)Ac1TPc1+Pc1Ad1=-PdM
The functionV(x)=xTPcMx is a common Lyapunov function for each of the individual system x.=Acx(t), x(k+1)=Adx(k) and hence a Lyapunov function for the switching system.Remark 11.
Same as Remark9.Clearly, in Theorems8 and 10, as a necessary condition, all the subsystems that are asymptotically stable still cannot ensure a stable switched system under arbitrary switching. Then, the sufficient condition is needed that is pairwise commutative, which guarantees the stability of the switched system. In other words, if there exists a common Lyapunov function for all the subsystems, the stability is ensured under arbitrary switching, and the solutions are gained.In the following section, matricesE is analysed, and some lemmas and theories are derived. The exponential stability of the system under any arbitrary switching is discussed, and a common Lyapunov function can be obtained.


## 2.1. System Description and Examples

Throughout, the following notation is adopted:R and C denote the fields of real and complex number, respectively; Rn denotes then-dimensional real Euclidean space; Rn×n denotes the space of n×n matrices with real entries;xi andxj denote theith andjth components of the vectorx in Rn, respectively; acpp′ and adqq′ denote the entry in the (p,p′) and (q,q′) position of the matricesA orE in Rn×n, respectively.There are three examples to show the multiswitched system and the parallel structure.Example 1.
Consider a setA=A1,A2,A3 of a switched system, which are the constant matrices of 3 subsystems, respectively. Assume in the classical switched system that there are 3 subsystems. The system switches between A1, A2, and A3. Assume in the multiswitched system that there are 7 subsystems-groups. The system switches betweenE,(5)E=E1,E2,E3,E4,E5,E6,E7=A100000000,0000A20000,00000000A3,A1000A20000,A10000000A3,0000A2000A3,A1000A2000A3Here,Nc=CcM1+⋯+CcMcM and Nd=CdM1+⋯+CdMdM.Remark 2.
A continuous-time subsystems-group composed of a subsystem or all the subsystems. Obviously, there are subsystems, but the value is null. It is same with a discrete-time subsystems-group.Remark 3.
Continuous-time subsystems cannot mix with discrete-time subsystems to be a subsystems-group. Continuous-time subsystems-groups and discrete-time subsystems-groups must be distinguished.Remark 4.
WhenEci is a singular matrix, zero rows or zero columns remain with the aim of uniform description. In Example 1, E1=A1, but the uniform description is kept to show E1=A1000A2000A3=A100000000, with A2=0, A3=0. Thus, allEci values look the same, ci=c1,…,cNc. The same is true of Edj.Remark 5.
Equations (2) and (4) are standard parallel structures. A parallel-like structure will be considered in Section 4, which has a coupling phenomenon based on a parallel structure between some subsystems. Assumptions are made for reducing conservativeness to obtain ideal stability.Example 6.
In a chilled water system of a central air conditioning system, there are three pumps driving the chilled water from the evaporator to the air conditioning unit (see Figure1). All the pumps can be switched ON/OFF and the variable frequency obeying a range of 50%-100% rated frequency. For setting an energy-saving control strategy with variable water volume technology, the whole volume cannot be less than 50% of the rated volume. Assuming in a time interval that the cold load changes, the switching, and pump water volume are shown in Figures 2 and 3. The parameters and symbols are shown in Table 2.
It is assumed that the cold load in a certain area requires air conditioning to meet the cold demand within 5 minutes to achieve a balanced state. Then, in the first 5 minutes, 3 pumps work at the rated frequency. After reaching the preset value, the cold load is stable with small fluctuations. In this period, pumpA andB exit; only pumpC works in the form of frequency conversion, maintaining the driving chilled water and transferring cooling capacity. In the last 3 minutes, due to the cold load increase, a pump is not enough to meet the demand, and the system puts pumpA into work with pumpC together with frequency conversion.
From the three time intervals ([0,5min], [5min,12min], and [12min,15min]), we find that there are three, one, and two subsystems working, respectively. According to the framework in this paper, the three pumps denote different continuous-time subsystems, and the work combinations of the pumps denote the different continuous-time subsystems-groups. Assuming the pumpA,B, andC denote the subsystem X1, X2, and X3, respectively; and the combinationA/B/C,C andA/C we just said denote the subsystems-group Y1, Y2, and Y3, respectively.Figure 1
Chilled water system of central air conditioning.Figure 2
Total flow volume.Figure 3
Water volumes of pumps.Example 7.
In the same system with Example6, but the pumps are divided into one variable frequency pump and two switchable fixed frequency pumps. Thus, the former is a continuous-time subsystem, and the latter are discrete-time subsystems. According to the energy-saving control strategy, assuming pumpC is always working with variable frequency, and pumpsA andB are switched ON or OFF with fixed frequency. Thus, the chilled water system is composed of continuous and discrete-time subsystems. The parameters and symbols are shown in Table 3.There are two simulation figures (Figures4 and 5) that illustrate the switching and working situations of the system.Figure 4
Total flow volume.Figure 5
Pump water volume.In the first 5 minutes, 3 pumps work together; this seems a subsystems-group with one continuous and two discrete-time subsystems. In the interval of from 5 min to 12 min, only pumpC works, which is the only one continuous-time subsystem that is a subsystems-group. In the last three minutes, pumpA works with pumpC; therefore, the subsystems-group has one continuous and one discrete-time subsystem.


## 2.2. Important Theories for Stability Analysis

We analyse the stability of the multiswitched system based on the structural property ofE-matrices, which cannot exist without the basic component,A-matrices. All the inferences of multiswitched system are derived from a classical switched system. Thus, some important theorems and propositions of the classical switched system should be quoted first.Theorem 8 (see [14]).
Consider a classical continuous-time switched systemx.=Acx(t) with Ac≜Ac1,Ac2,…,AcM, where the matrices Acp are asymptotically stable and commute pairwise. Then, we get the following:(1)
The system is exponentially stable under any arbitrary switching between the elements ofAc.(2)
There exists a common Lyapunov function for all the subsystems. For any positive symmetric definite matrixPc0, let Pc1,Pc2,…,PcM be the unique symmetric positive definite solutions to the Lyapunov equations(6)AcpTPcp+PcpAcp=-Pcp-1,cp=c1,c2,…,cM
The functionV(x)=xTPcMx is a common Lyapunov function for each of the individual system x.=Acx(t) and hence a Lyapunov function for the switching system.Remark 9.
All of the proposed descriptions and notations are identical to the given multisystem.Theorem 10 (see [1]).
Consider a classical mixed switched system composed of continuous-time subsystemsx.=Acx(t) with Ac≜Ac1,Ac2,…,AcM and discrete-time subsystems x(k+1)=Adx(k) with Ad≜Ad1,Ad2,…,AdM. Let Acp be Hurwitz stable and Adq be Schur stable to commute pairwise. Then, we get the following:(1)
The system is exponentially stable for any arbitrary switching between the elements ofA.(2)
There exists a common Lyapunov function for all the subsystems. For any positive definite matrixPd0, let Pd1,Pd2,…,PdM be the unique positive definite solutions to the Lyapunov equations(7)Ad1TPd1Ad1+Pd1=-Pd0(8)AdqTPdqAdq+Pdq=-Pdq-1
letPc1,Pc2,…,PcM be the unique positive definite solutions to the Lyapunov equations(9)AcpTPcp+PcpAcp=-Pcp-1(10)Ac1TPc1+Pc1Ad1=-PdM
The functionV(x)=xTPcMx is a common Lyapunov function for each of the individual system x.=Acx(t), x(k+1)=Adx(k) and hence a Lyapunov function for the switching system.Remark 11.
Same as Remark9.Clearly, in Theorems8 and 10, as a necessary condition, all the subsystems that are asymptotically stable still cannot ensure a stable switched system under arbitrary switching. Then, the sufficient condition is needed that is pairwise commutative, which guarantees the stability of the switched system. In other words, if there exists a common Lyapunov function for all the subsystems, the stability is ensured under arbitrary switching, and the solutions are gained.In the following section, matricesE is analysed, and some lemmas and theories are derived. The exponential stability of the system under any arbitrary switching is discussed, and a common Lyapunov function can be obtained.


## 3. Main Results

In the framework of proposed multiswitched system, and based on Theorems8 and 10, we have the following lemmas.Lemma 12.
If matricesAc andAd are pairwise commutative, matricesEc andEd are also pairwise commutative.Proof.
ConsiderAciAci′=Aci′Aci of continuous-time subsystems, i≠i′.
We get(11)EciEci′=Eci′Eci=Ac1⋯0⋯0⋮⋱⋮0Acp0⋮⋱⋮0⋯0⋯AcMAc1⋯0⋯0⋮⋱⋮0Acp0⋮⋱⋮0⋯0⋯AcM
ConsiderAcAd=AdAc of all subsystems; we get(12)EcEd=Ac1⋯0⋯0⋮⋱⋮0Acp0⋮⋱⋮0⋯0⋯AcMAd1⋯0⋯0⋮⋱⋮0Adq0⋮⋱⋮0⋯0⋯AdM=Ac1Ad1⋯0⋯0⋮⋱⋮0AcpAdp0⋮⋱⋮0⋯0⋯AcMAdM=Ad1Ac1⋯0⋯0⋮⋱⋮0AdqAcp0⋮⋱⋮0⋯0⋯AdMAcM=EdEc
E d i E d i ′ = E d i ′ E d i is similar with (11).
Therefore, the matricesEc andEd are pairwise commutative. This completes the proof for Lemma 12.Lemma 13.
If matricesAc andAd are Hurwitz stable and Schur stable, respectively, matricesEc andEd are also Hurwitz stable and Schur stable, respectively.Proof.
According to the Hurwitz stable criterion, all the order principal minors of matricesEmust be positive.
The structure of (2) is parallel, and matricesEare diagonal. Thus, all the order principal minors are also diagonal. Acp, cp=c1,…,cM are Hurwitz stable, and all the order principal minors are positive; thus, it is proven that matricesE are Hurwitz stable.
E d is Schur stable, which is similar with the Hurwitz stable criterion.Remark 14.
Unnecessary zero rows or zero columns can be omitted.Based on Theorem8 and Lemmas 12 and 13, we have the following results.Proposition 15.
Assumex.=Acx(t) are Hurwitz stable andAcp are pairwise commutative. Then, we get the following:(1)
The continuous-time multiswitched system is exponentially stable for any arbitrary switching between the elements ofEc.(2)
There exists a common Lyapunov function for all the subsystems-groups and subsystems. For any positive symmetric definite matrixPc0, let Pc1,Pc2,…,PcM be the unique symmetric positive definite solutions to the Lyapunov equations(13)Ac1TPi+PiAc1+Ac2TPi+PiAc2+⋯+AcMTPi+PiAcM=-Pi-1,cp=c1,c2,…,cM;i=1,2,…,Nc
The functionV(x)=xTPNcx is a common Lyapunov function for each of the individual systems x.=Acx(t) and is thus a Lyapunov function for the switching system.Proof.
Based on Lemmas12 and 13,Eci are pairwise commutative, andEc are Hurwitz stable. Combining Theorem 8, we get(14)EciTPi+PiEci=-Pi-1,i=1,…,NcWe substituteEci in (2) into (14) to obtain(15)Ac1⋯0⋯0⋮⋱⋮0Acp0⋮⋱⋮0⋯0⋯AcMTPi+PiAc1⋯0⋯0⋮⋱⋮0Acp0⋮⋱⋮0⋯0⋯AcM=-Pi-1
Equation (15) equals (14), which implies that the system is exponentially stable for arbitrary switching; (15) transforms into (13), which implies the solutions to the Lyapunov equations.Example 16.
Consider a setAc=Ac1,Ac2,Ac3,Ac4,Ac5 of a switched system, which are the constant matrices of 5 continuous-time subsystems. The multiswitched system switches between subsystems-groupsEc. Here, we select only three subsystems-groups as an example.
LetEc=E1,E2,E3, E1=Ac1, E2=Ac200Ac3, and E3=Ac3000Ac4000Ac5.
Then, for any positive symmetric definite matrixPc0, there are Pc1, Pc2,andPc3 as the unique symmetric positive definite solutions to the Lyapunov equations(16a)Ec1TP1+P1Ec1=-Pc0(16b)Ec2TP2+P2Ec2=-Pc1(16c)Ec3TP3+P3Ec3=-Pc2
Equations (16a), (16b), and (16c) can be rewritten as(17a)Ac1TP1+P1Ac1=-Pc0(17b)Ac2TP2+P2Ac2+Ac3TP2+P2Ac3=-Pc1(17c)Ac3TP3+P3A34+Ac4TP3+P3Ac4+Ac5TP3+P3Ac5=-Pc2
The functionV(x)=xTP3x is a common Lyapunov function for each of the individual system x.=Acx(t) with cp=c1,c2,c3,c4,and c5.Remark 17.
The above theorem and proof can be extended to the case of discrete-time multiswitched systemx(k+1)=Edjx(k) with j=1,…,Nd. Assume the matrices x(k+1)=Adqx(k) with dq=d1,…,dM are Schur stable and commute pairwise. Then, the discrete-time multiswitched system is exponentially stable for arbitrary switching between the elements ofEd. The solution (13) can be modified as(18)Ad1TPjAd1+Ad2TPjAd2+⋯=AdMTPjAdM=Pj-Pj-1,dq=d1,d2,…,dM;j=1,…,Ndand the common Lyapunov function is V(x)=xTPNdx.Proposition 18.
Assumex.=Acx(t) are Hurwitz stable and x(k+1)=Adx(k) are Schur stable; Acp and Adq are pairwise commutative. Then, we get the following:(1)
The multiswitched system is exponentially stable for any arbitrary switching between the elements ofE.(2)
There exists a common Lyapunov function for all the subsystems. For any positive symmetric definite matrixPd0, let Pd1,Pd2,…,PdNd be the unique positive definite solutions to the Lyapunov equations(19a)Ad1TPd1Ad1+⋯+AdMTPd1AdM-Pd1=-Pd0,dq=d1,d2,…,dM(19b)Ad1TPdjAd1+⋯+AdMTPdjAdM-Pdj=-Pdj-1,dq=d1,d2,…,dM;j=1,…,NdLet Pc1,Pc2,…,PcN be the unique positive definite solutions to the Lyapunov equations(19c)Ac1TPc1+Pc1Ad1+⋯+AcMTPc1+Pc1AdM=-PdNd,cp=c1,c2,…,cM;dq=d1,d2,…,dM(19d)Ac1TPci+PciAc1+⋯+AcMTPci+PciAcM=-Pci-1,cp=c1,c2,…,cM;i=1,2,…,Nc
The functionV(x)=xTPcNcx is a common Lyapunov function for each of the individual system x.=Acx(t), and x(k+1)=Adqx(k); hence a Lyapunov function is used for the switching system.Proof.
Based on Theorem10 as well as Lemmas 12 and 13, we get(20a)Ed1TPd1Ed1-Pd1=-Pd0(20b)EdjTPdjEdj-Pdj=-Pdj-1(20c)Ec1TPc1+Pc1Ed1=-PdNd(20d)EciTPci+PciEci=-Pci-1We substitute Eci and Edj in (2) into (19a), (19b), (19c), and (19d) to obtain(21a)Ad1⋯0⋯0⋮⋱⋮0Adq0⋮⋱⋮0⋯0⋯AdMTPd1Ad1⋯0⋯0⋮⋱⋮0Adq0⋮⋱⋮0⋯0⋯AdM-Pd1=-Pd0(21b)Ad1⋯0⋯0⋮⋱⋮0Adq0⋮⋱⋮0⋯0⋯AdMTPdjAd1⋯0⋯0⋮⋱⋮0Adq0⋮⋱⋮0⋯0⋯AdM-Pdj=-Pdj-1(21c)Ac1⋯0⋯0⋮⋱⋮0Acp0⋮⋱⋮0⋯0⋯AcMTPc1+Pc1Ad1⋯0⋯0⋮⋱⋮0Adq0⋮⋱⋮0⋯0⋯AdM=-PdNd(21d)Ac1⋯0⋯0⋮⋱⋮0Acp0⋮⋱⋮0⋯0⋯AcMcpTPci+PciAc1⋯0⋯0⋮⋱⋮0Acp0⋮⋱⋮0⋯0⋯AcM=-Pci-1 Equations (21a), (21b), (21c), and (21d) equal (20a), (20b), (20c), and (20d), which implies that the system is exponentially stable for any arbitrary switching; and (21a), (21b), (21c), and (21d) transform into (19a), (19b), (19c), and (19d), which implies the solutions to the Lyapunov equations.Example 19.
Consider setsAc=Ac1,Ac2,Ac3 and Ad=Ad1,Ad2 of a switched system, which are the constant matrices of 3 continuous-time subsystems and 2 discrete-time subsystems, respectively. The multiswitched system switches between in subsystems-groupsE. Here, we select only five subsystems-groups as an example.
LetEc=Ec1,Ec2,Ec3, Ec1=Ac100Ac3, Ec2=Ac200Ac3, and Ec3=Ac1000Ac2000Ac3; let Ed=Ed1,Ed2, Ed1=Ad2, and Ed2=Ad100Ad2.
Then, for any positive definite matrixPd0, there are Pd1, Pd2, Pc1, Pc2,andPc3 as the unique symmetric positive definite solutions to the Lyapunov equations(22a)Ed1TPd1Ed1-Pd1=-Pd0(22b)Ed2TPd2Ed2-Pd2=-Pd1(22c)Ec1TPc1+Pc1Ed1=-Pd2(22d)Ec2TPc2+Pc2Ec2=-Pc1(22e)Ec3TPc3+Pc3Ec3=-Pc2
Equations (22a), (22b), (22c), (22d), and (22e) can be rewritten as(23a)Ad2TPd1Ad2-Pd1=-Pd0(23b)Ad1TPd2Ad1+Ad2TPd2Ad2-Pd2=-Pd1(23c)Ac1TPc1+Ac2TPc1+Pc1Ad2=-Pd2(23d)Ac2TPc2+Pc2Ac2+Ac3TPc2+Pc2Ac3=-Pc1(23e)Ac1TPc3+Pc3Ac1+Ac2TPc3+Pc3Ac2+Ac3TPc3+Pc3Ac3=-Pc2
The functionV(x)=xTPc3x is a common Lyapunov function for each of the individual systems x.=Acx(t) with cp=c1,c2,c3 and x(k+1)=Adx(k) with dq=d1,d2, respectively.


## 4. Results in Parallel-Like Structures

Lemma 20.
If the structures ofEc are not standard parallel and contain coupled components based on a parallel framework (i.e., the subsystems are not independent of each other), then the matrices should be modified as(24)Eci=Ac1⋯Ao1,p⋯Ao1,cM⋮⋱⋮Aop,1AcpAop,cM⋮⋱⋮AocM,1⋯AocM,p⋯AcMAo(p,p′) are matrices of the coupled components, p=c1,…,cM, and p≠p′. Clearly, there are another nonzero elements except in the main diagonal line.Remark 21.
Considering the coupled components,Ed can be modified as(25)Edj=Ad1⋯Ao1,q⋯Ao1,dM⋮⋱⋮Aoq,1AdqAoq,dM⋮⋱⋮AodM,1⋯AodM,q⋯AdM
A o ( q , q ′ ) are matrices of the coupled components, q=q1,…,qM, and q≠q′.Remark 22.
The structures proposed in Lemma20 and Remark 21 are named as parallel-like. Only the coupled component is the distinction of structure between parallel and parallel-like; in other words, the parallel-like structure is minor alteration on the parallel structure.Proposition 23.
Under the assumption of Lemma20, the framework ofEc in Lemma 20 can be classified as the following 3 types: upper triangular, lower triangular, and other nonregular modes. Then, Proposition 15 can be modified as follows:
(a) If the framework ofEc is upper triangular, Proposition 15 is true, but the solutions is(26a)Ac1TPi+PiAc1⋯PiAo1,p⋯PiAo1,cM⋮⋱⋮Ao1,pTPiAcpTPi+PiAcpPiAop,cM⋮⋱⋮Ao1,cMTPi⋯Aop,cMTP⋯AcMTPi+PiAcM=-Pi-1where cp=c1,c2,…,cM; i=1,2,…,Nc. The function V(x)=xTPNcx is a common Lyapunov function for each of the individual system x.=Acx(t).
(b) Lower triangular is same as upper triangular. The solution is(26b)Ac1TPi+PiAc1⋯Aop,1TPi⋯AocM,1TPi⋮⋱⋮PiAo1,pAcpTPi+PiAcpAocM,pTPi⋮⋱⋮PiAo1,cM⋯PiAop,cM⋯AcMTPi+PiAcM=-Pi-1(c) Other nonregular modes must satisfy Hurwitz stable, only this, which have ideal stability. The solution is(26c)Ac1TPi+PiAc1⋯Aop,1TPi+PiAo1,p⋯AocM,1TPi+PiAo1,cM⋮⋱⋮Ao1,pTPi+PiAo1,pAcpTPi+PiAcpAocM,pTPi+PiAop,cM⋮⋱⋮Ao1,cMTP+PiAo1,cM⋯Aop,cMTP+PiAop,cM⋯AcMTPi+PiAcM=-Pi-1Proof.
The condition ofEc that is pairwise commutative is determined as follows:(27)EciEci′=Eci′Eci=Ac1⋯Ao1,p⋯Ao1,cM⋮⋱⋮Aop,1AcpAop,cM⋮⋱⋮AocM,1⋯AocM,p⋯AcMAc1⋯Ao1,p⋯Ao1,cM⋮⋱⋮Aop,1AcpAop,cM⋮⋱⋮AocM,1⋯AocM,p⋯AcM
The Hurwitz stability ofEc can be ensured by the framework of upper triangular, in which all the order principal minors of matrices are positive, ifAci is Hurwitz stable. It satisfies the Hurwitz stable criterion, which is the same as the framework of the lower triangular.
Thus, the two frameworks of a continuous-time multiswitched system are exponentially stable for any arbitrary switching between the elements ofEc. There exists a common Lyapunov function for all the subsystems-groups and subsystems. In the framework of the upper triangular, we modify (14) as follows:(28)
A
c
1
⋯
A
o
1
,
p
⋯
A
o
1
,
c
M
⋮
⋱
⋮
0
0
A
o
p
,
c
M
⋮
⋱
⋮
0
⋯
0
⋯
A
c
M
T
P
i
+
P
i
A
c
1
⋯
A
o
1
,
p
⋯
A
o
1
,
c
M
⋮
⋱
⋮
0
A
c
p
A
o
p
,
c
M
⋮
⋱
⋮
0
⋯
0
⋯
A
c
M
=
-
P
i
-
1
Obviously, (28) transforms into (26a), which implies the solutions to the Lyapunov equations.
It is similar with upper triangular, in the framework of the lower triangular, we modify (14) as(29)Ac1⋯0⋯0⋮⋱⋮Aop,1Acp0⋮⋱⋮AocM,1⋯AocM,p⋯AcMTPi+PiAc1⋯0⋯0⋮⋱⋮Aop,1Acp0⋮⋱⋮AocM,1⋯AocM,p⋯AcM=-Pi-1
Of course, (29) transforms into (26b), which implies the solutions to the Lyapunov equations.
However, in other nonregular modes, the stability cannot be guaranteed. The framework is nonregular, so theEc Hurwitz needs to be stable. Then, the system stability under arbitrary switching is ensured, and a common Lyapunov function can be gained as (26c).Remark 24.
The above theorem and proof can be extended to the structure (see (25)) of the discrete-time multiswitched system x(k+1)=Edjx(k) with j=1,…,Nd. Remark 17 can be modified as follows:
(a) If the framework ofEd is an upper triangular, Remark 17 is true; however, the solutions are(30a)Ad1⋯Ao1,q⋯Ao1,dM⋮⋱⋮0AdqAoq,dM⋮⋱⋮0⋯0⋯AdMTPiAd1⋯Ao1,q⋯Ao1,dM⋮⋱⋮0AdqAoq,dM⋮⋱⋮0⋯0⋯AdM=Pi-Pi-1 where dq=d1,d2,…,dM; j=1,2,…,Nd. The function V(x)=xTPNdx is a common Lyapunov function for each of the individual systems x.=Adx(t).
(b) Lower triangular is same as upper triangular. The solution is(30b)Ad1⋯0⋯0⋮⋱⋮Aoq,1Adq0⋮⋱⋮AodM,1⋯AodM,q⋯AdMTPjAd1⋯0⋯0⋮⋱⋮Aoq,1Adq0⋮⋱⋮AodM,1⋯AodM,q⋯AdM=Pj-Pj-1(c) Other nonregular modes must satisfy Schur stable, which has ideal stability. The solution is(30c)Ad1⋯Ao1,q⋯Ao1,dM⋮⋱⋮Aoq,1AdqAoq,dM⋮⋱⋮AodM,1⋯AodM,q⋯AdMTPjAd1⋯Ao1,q⋯Ao1,dM⋮⋱⋮Aoq,1AdqAoq,dM⋮⋱⋮AodM,1⋯AodM,q⋯AdM=Pj-Pj-1
If the above theorems and remarks extend to mixed switched system, we find it difficult to get the condition of pairwise commutative, whatever the triangular and other frameworks inE. In some special situations, it satisfies the assuming condition EcEd=EdEc. Then, the matricesEc andEd should better be in upper or lower triangular framework, to ensure Hurwitz stable and Schur stable, respectively. Otherwise, the stable condition becomes strictly to require matricesE and be unconcerned with matricesA.
If all the conditions are satisfied, the goal of stability under arbitrary switching can be gained, and the functionV(x)=xTPcNcx is a common Lyapunov function for each of the individual system x.=Acx(t), and x(k+1)=Adqx(k). For any positive symmetric definite matrix Pd0, thus, Pd1,Pd2,…,PdNd are the unique positive definite solutions to the Lyapunov equations(31a)Ad1⋯Ao1,q⋯Ao1,dM⋮⋱⋮Aoq,1AdqAoq,dM⋮⋱⋮AodM,1⋯AodM,q⋯AdMTPd1Ad1⋯Ao1,q⋯Ao1,dM⋮⋱⋮Aoq,1AdqAoq,dM⋮⋱⋮AodM,1⋯AodM,q⋯AdM-Pd1=-Pd0(31b)Ad1⋯Ao1,q⋯Ao1,dM⋮⋱⋮Aoq,1AdqAoq,dM⋮⋱⋮AodM,1⋯AodM,q⋯AdMTPdjAd1⋯Ao1,q⋯Ao1,dM⋮⋱⋮Aoq,1AdqAoq,dM⋮⋱⋮AodM,1⋯AodM,q⋯AdM-Pdj=-Pdj-1(31c)Ac1⋯Ao1,p⋯Ao1,cM⋮⋱⋮Aop,1AcpAop,cM⋮⋱⋮AocM,1⋯AocM,p⋯AcMTPc1+Pc1Ad1⋯Ao1,q⋯Ao1,dM⋮⋱⋮Aoq,1AdqAoq,dM⋮⋱⋮AodM,1⋯AodM,q⋯AdM=-PdNd(31d)Ad1⋯Ao1,q⋯Ao1,dM⋮⋱⋮Aoq,1AdqAoq,dM⋮⋱⋮AodM,1⋯AodM,q⋯AdMcpTPci+PciAd1⋯Ao1,q⋯Ao1,dM⋮⋱⋮Aoq,1AdqAoq,dM⋮⋱⋮AodM,1⋯AodM,q⋯AdM=-Pci-1


## 5. Numerical Example

In this section, an engineering application of central air conditioning is introduced as a numerical example, which shows the framework of multiswitched system with parallel structure. In addition, the engineering environment is illustrated and thermal parameters are given in Table4. Simulation results are presented to illustrate the characteristics of the system and the situation of stability by different control strategies. The mathematical model is given in [27], which shows the thermal balance of a test room affected by factors, such as the structure and materials of a building, outdoor weather parameters, indoor lighting, radiating equipment, and number of occupants (see Figure 6). The cooling capacity is transferred from chilled water system to air system via air conditioning units by the measures of constant air volume and variable water volume/constant temperature difference.Figure 6
Thermal balance of air system.The thermal balance equation is(32)Camadθdτ=-γεCwΔθqw-1-RrqsaCaθ+1-RrqsaCaθout+Qrd+Qqr-∑KjAjθ+∑KjAjθjwhere θ is the real-time indoor temperature. On the left side of the equation, Camadθ/dτ means the time differential of the heat capacity of a room. On the left side of the equation, -γεCwΔθqw means the cooling capacity for chilled water system; -(1-Rr)qsaCaθ and (1-Rr)qsaCaθout represent the cooling capacity from return and fresh air systems, respectively; Qrd denotes random heat of occupants and equipment; Qqr means latent heat inside; -∑KjAjθ+∑KjAjθj is the heat transfer from building structure. The description of the symbols is presented in Table 4.A midsize conference room (length: 10 m, width: 6 m, and height: 3 m) is simulated, we use two different control strategies (the strategies 1 and 2 mentioned in Examples6 and 7, respectively) to adjust cooling capacity and illustrate system stability (corresponding Propositions 15 and 18, respectively) in the framework of multiswitching system with parallel structure, and use strategy 3 to reflect the unstable situation. The above three control strategies and the two types of pumps are shown in Table 1. The range of variable volume is 50%-100% and the time is divided into three intervals ([0,5min], [5min,12min], and [12min,15min]) in the above three control strategies. The cold air is sent to the room for reducing the indoor temperature. The indoor temperature is required to be lowered with respect to the initial temperature (θ0, 30°C) and regulated at (θset, 26°C) as soon as possible. In the last 3 minutes, the indoor cooling load increased significantly due to that the number of indoor participants increased (see Figure 10). The outdoor temperature is basically maintained at 30°C in the simulated 15 minutes.Table 1
Control strategies and types of pumps.
strategy pump working type control mode feedback coefficient switching state 1 A variable volume feedback/switching 0.0074 ON/OFF B variable volume feedback/switching NULL ON/OFF C variable volume feedback/switching 0.021 ON 2 A fixed volume switching NULL ON/OFF B fixed volume switching NULL ON/OFF C variable volume feedback 0.019 ON 3 A fixed volume switching NULL ON/OFF B fixed volume switching NULL ON/OFF C variable volume feedback 0.019 ONTable 2 pump rated volumeq0 (kg/s) range of variable volumeq (kg/s) A 0.01 0.005-0.01 B 0.02 0.015-0.02 C 0.033 0.0175-0.033 total 0.063 0.0315-0.063Table 3 pump rated volumeq0 (kg/s) range of variable volumeq (kg/s) A 0.01 0.01-0.01 B 0.02 0.02-0.02 C 0.033 0.0175-0.033 total 0.063 0.0315-0.063Table 4 Parameter Value Description m a (kg) 232.18 indoor air mass q w (kg/s) 0.149 rated volume of water system q a (kg/s) 0.022 rated volume of pumpA q b (kg/s) 0.044 rated volume of pumpB q c (kg/s) 0.083 rated volume of pumpC q s a (kg/s) 3.003 volume of sending air A 1 (m2) 56 area of walls A 2 (m2) 28 area of windows A 3 (m2) 0 area of roof C a (J/kg∗k) 1010 specific heat of air C w (J/kg∗k) 4180 specific heat of water K 1 (W/m2∗k) 0.049 heat transfer coefficient of walls K 2 (W/m2∗k) 0.051 heat transfer coefficient of windows K 3 (W/m2∗k) 0.05 heat transfer coefficient of roof Q q r (J) 20 latent heat load R r 0.11 return air rate Δ θ (°C) 5 temperature difference ε 0.89 transfer efficiency from water system to air system γ 0.095 coefficient of cooling capacity allocation θ i n i (°C) 30, 30, 30 initial temperature θ j (°C) 35, 35, 36 temperature of walls, windows, and roof, respectively θ s e t (°C) 26 setting temperatureFigure7 shows the switching dynamics of a continuous-time multiswitched system with parallel structure by the flow volume of the three pumps. In the first time interval, all the three pumps work as a subsystems-group in rated volume for reducing the temperature; in the middle time interval, only the pumpC works in variable flow mode under the feedback control; in the last time interval, both pumpAand pumpC work as a subsystems-group under feedback control for creasing cooling capacity because of the increase of indoor participants increased.Figure 7
Flow volume under strategy 1.Figure8 shows the switching of a mixed multiswitched system with parallel structure composed of one continuous-time subsystem and two discrete-time subsystems. In the first two time intervals, the dynamics of pumps is the same as Figure 7. In the last time interval, pumpA works at the rated volume as a discrete-time subsystem and the pumpC works in variable flow mode under feedback control as a continuous-time subsystem; in other words, the subsystems-group is composed with one continuous-time and one discrete-time subsystems. Figure 9 is similar to Figure 8, but the switching dynamics is different (in the middle time interval, all the three pumps work together; in the last time interval, only pumpC works).Figure 8
Flow volume under strategy 2.Figure 9
Flow volume under strategy 3.Figure 10
Change of cooling load.Figure11 shows the changes of indoor temperature under the three control strategies. The indoor temperature drops from the initial value (30°C) to the set point (26°C) in five minutes under the three strategies, because of rated volume by the whole water system. It is worth noting that, in the last 10 minutes, the indoor temperature is always stable near the set point under the ideal strategies (the strategies 1 and 2), even if the indoor cooling load changes significantly; but under the strategy 3, the indoor temperature expresses lower and higher instable rates in the last two time intervals, respectively.Figure 11
Indoor temperature under different strategies.


## 6. Conclusion

In this paper, a type of linear multiswitched system with parallel structure was proposed, and the framework and a switching unit were introduced. Based on this, various actual engineering applications were shown, which illustrated the properties of the system and differences with traditional switched systems. Next, the stability property for a type of linear multiswitched system with parallel structure is studied, whether in continuous-time, discrete-time, or a mixed situation. A subsystems-group as a basic switched unit instead of subsystem is proposed, the matrices of which are pairwise commutative based on some given conditions of subsystems. When all the subsystems are Hurwitz and Schur stable, there exists a common Lyapunov function for all the subsystems and subsystems-groups. Then, the switched system is exponentially stable for any arbitrary switching between the subsystems-groups. The results are extended to a parallel-like structure to obtain more ideal consequence of stability. A simulation example for refrigeration engineering application of the system is introduced as last, which shows the characteristics of the framework and stability.



---

*Source: 3834601-2018-08-08.xml*