Finite‐time dissipative control for time‐delay Markov jump systems with conic‐type non‐linearities under guaranteed cost controller and quantiser
2020; Institution of Engineering and Technology; Volume: 15; Issue: 4 Linguagem: Inglês
10.1049/cth2.12031
ISSN1751-8652
AutoresXiang Zhang, Yanyan Yin, Hai Wang, Shuping He,
Tópico(s)Neural Networks Stability and Synchronization
ResumoIET Control Theory & ApplicationsVolume 15, Issue 4 p. 489-498 ORIGINAL RESEARCH PAPEROpen Access Finite-time dissipative control for time-delay Markov jump systems with conic-type non-linearities under guaranteed cost controller and quantiser Xiang Zhang, Xiang Zhang Key Laboratory of Intelligent Computing and Signal Processing (Ministry of Education), School of Electrical Engineering and Automation, Anhui University, Hefei, ChinaSearch for more papers by this authorYanyan Yin, Yanyan Yin School of Electrical Engineering, Computing and Mathematical Sciences, Curtin University, Perth, AustraliaSearch for more papers by this authorHai Wang, Hai Wang Discipline of Engineering and Energy, Murdoch University, Murdoch, AustraliaSearch for more papers by this authorShuping He, Corresponding Author Shuping He shuping.he@ahu.edu.cn Key Laboratory of Intelligent Computing and Signal Processing (Ministry of Education), School of Electrical Engineering and Automation, Anhui University, Hefei, China Correspondence Shuping He, Key Laboratory of Intelligent Computing and Signal Processing (Ministry of Education), School of Electrical Engineering and Automation, Anhui University, Hefei, China. Email: shuping.he@ahu.edu.cnSearch for more papers by this author Xiang Zhang, Xiang Zhang Key Laboratory of Intelligent Computing and Signal Processing (Ministry of Education), School of Electrical Engineering and Automation, Anhui University, Hefei, ChinaSearch for more papers by this authorYanyan Yin, Yanyan Yin School of Electrical Engineering, Computing and Mathematical Sciences, Curtin University, Perth, AustraliaSearch for more papers by this authorHai Wang, Hai Wang Discipline of Engineering and Energy, Murdoch University, Murdoch, AustraliaSearch for more papers by this authorShuping He, Corresponding Author Shuping He shuping.he@ahu.edu.cn Key Laboratory of Intelligent Computing and Signal Processing (Ministry of Education), School of Electrical Engineering and Automation, Anhui University, Hefei, China Correspondence Shuping He, Key Laboratory of Intelligent Computing and Signal Processing (Ministry of Education), School of Electrical Engineering and Automation, Anhui University, Hefei, China. Email: shuping.he@ahu.edu.cnSearch for more papers by this author First published: 29 December 2020 https://doi.org/10.1049/cth2.12031Citations: 2AboutSectionsPDF ToolsRequest permissionExport citationAdd to favoritesTrack citation ShareShare Give accessShare full text accessShare full-text accessPlease review our Terms and Conditions of Use and check box below to share full-text version of article.I have read and accept the Wiley Online Library Terms and Conditions of UseShareable LinkUse the link below to share a full-text version of this article with your friends and colleagues. Learn more.Copy URL Share a linkShare onFacebookTwitterLinked InRedditWechat Abstract For a class of conic-type non-linear time-delay Markov jump systems, the asynchronous dissipative output feedback controller based on the guaranteed cost control and quantiser is designed in this study. In real applications, the system and the controller modes are always non-synchronous, so we introduce the hidden Markov model to solve this problem. Furthermore, we define three novel auxiliary variables and use quantisers to accomplish the output feedback controller design. Then, the finite-time boundedness and strict dissipativity of the closed-loop systems are guaranteed by sufficient conditions, and the controller also meets the guaranteed cost-control performance. By solving a set of linear matrix inequalities, we get the controller gains, the guaranteed cost control performance index J*, and the dissipative performance index α. Finally, the correctness and feasibility of this designed approach are demonstrated by a given example. 1 INTRODUCTION In the 1960s, the conception of Markov jump systems (MJSs) [1] has been proposed. With the random jumping structure, MJSs can be seen as a class of hybrid systems and have always received considerable attention. Due to the randomness in data and structure, MJSs have been widely applied, such as multiagent systems [2], electrohydraulic servo systems [3], medical prognosis [4] and networked systems [5]. But in engineering applications, it is difficult to synchronise the system and the controller modes due to some unavoidable errors and delays. For this asynchronous phenomenon different from synchronisation [6], the hidden Markov model (HMM) might be helpful which introduces a random process to estimate the Markov process [7]. In [8], the finite-time boundedness (FTB) and H∞ performance of time-varying MJSs were investigated by designing an HMM-based controller. For time-delay MJSs (TDMJSs), the asynchronous controller was designed to guarantee the stochastic stability and H∞ performance [9]. For fuzzy MJSs, the asynchronous filtering was designed to guarantee the FTB and H ∞ performance [10]. In [11], the authors designed an HMM-based sliding mode controller to ensure the stochastic stability and dissipative performance for TDMJSs. On the other hand, many achievements have been made in the design of robust stabilising controllers [12]. Although it can make the system stable, the upper bound of the controller performance cannot be guaranteed. For this problem, the guaranteed cost control (GCC) strategy has been put forward [13]. For discrete-time switched singular systems, the robust GCC problem was investigated [14]. For fuzzy MJSs, the GCC performance was guaranteed by designing a quantised asynchronous controller [15]. The event-triggered GCC strategy was investigated in [16, 17]. In practical applications, because of disturbances and errors, the non-linear characteristics are indispensable. As a special kind of non-linear dynamics, conic-type non-linearities have great representativeness and they are widely used in engineering, such as locally sinusoidal non-linearities, dead-zone non-linearities and so forth. In fact, we can consider Lipschitz non-linearity as a special conic-type non-linearity. For discrete-time conic non-linear MJSs [18], the authors investigated the FTB and H ∞ performance by designing an asynchronous controller with a DC-motor model to show the feasibility. The stability analysis of the conic non-linear systems was studied [19, 20]. In [21], the fault detection of conic non-linear systems was investigated. The dissipative theory [22], since its proposal by Willems, has received much attention and has been widely used, for example, circuit analysis [23-26], neural networks [27-30], filtre design [31-34]. In [35], the authors achieved the stochastic stability and extended dissipativity of fuzzy switched systems by designing an asynchronous controller. For discrete-time fuzzy MJSs, the stochastic stability and strict dissipativity were studied by state-feedback controlling [36]. For linear MJSs, the authors designed an asynchronous time-delay controller to investigate the dissipativity [37]. In [38], the HMM-based dissipative control scheme for discrete-time TDMJSs was studied. In practice, the packet loss and time delay will reduce the stability and performance of the system. Due to the limited transmission rate of network, and the ability to map continuous signals to discrete sets, quantiser is needed and it is cheaper, more reliable and convenient. For the results of quantiser, the readers can refer [9, 10, 37]. To our best knowledge, the asynchronous dissipative control problem for TDMJSs with conic-type non-linearities under guaranteed cost controller and quantiser has not been fully studied. In this study, we introduced the HMM and quantiser to design the output feedback controller for a class of TDMJSs with conic-type non-linearities. The following points reflect the main contributions of this study: In order to make the matrix inequalities solvable and to reduce the computational complexity, we defined the three novel auxiliary variables to accomplish the controller design. By solving a set of linear matrix inequalities (LMIs), sufficient conditions are given to guarantee the FTB and strict dissipativity of the closed-loop systems, and the upper bound of the controller performance is also guaranteed. By the given liquid monopropellant rocket motor model with a pressure feeding system, the correctness and feasibility of the designed strategy are guaranteed. In the following, Table 1 introduces the presented notations in this study. TABLE 1. The notations Notation Denotes E{⋅} The mathematical expectation operator εmax(A) The maximum eigenvalue of A εmin(A) The minimum eigenvalue of A ℜn n-dimensional Euclidean space ℜn×m n×m real matrix diag{A B} Block-diagonal matrix of A and B I Unit matrix A − 1 Matrix inverse A T Matrix transpose * Symmetric matrix Her(A) The sum of A and transposition of A 2 SYSTEM DESCRIPTION AND PRELIMINARIES Consider the following TDMJSs with conic-type non-linearities: x ̇ ( t ) = f ( x ( t ) , x ( t − τ ) , ω ( t ) ) + D r ( t ) u ( t ) , z ( t ) = E r ( t ) x ( t ) + F r ( t ) x ( t − τ ) + G r ( t ) ω ( t ) , (1)where x ( t ) ∈ ℜ n is the system state, u ( t ) ∈ ℜ m is the controlled input, ω ( t ) ∈ ℜ q is the external disturbance with ω T ( t ) ω ( t ) ≤ ϖ ( t ) , z ( t ) ∈ ℜ p is the controlled output. f ( x ( t ) , x ( t − τ ) , ω ( t ) ) is an unknown non-linear function by the following dynamics conic sector: f ( x ( t ) , x ( t − τ ) , ω ( t ) ) − [ A r ( t ) x ( t ) + B r ( t ) x ( t − τ ) + C r ( t ) ω ( t ) ] ≤ A a r ( t ) x ( t ) + A b r ( t ) x ( t − τ ) + A c r ( t ) ω ( t ) . (2)The values of the Markov stochastic process { r ( t ) , t ≥ 0 } are in a finite set L = { 1 , 2 , … , L } with the transition rate matrix Π = [ λ s l ] given by P { r ( t + Δ t ) = l | r ( t ) = s } = λ s l Δ t + o ( Δ t ) , s ≠ l 1 + λ s s Δ t + o ( Δ ) , s = l (3)where Δ t satisfies lim Δ t → 0 o ( Δ t ) Δ t = 0 , λ s l ≥ 0 represents the jump rate from mode s at time t to mode l at time t + Δ t and λ s s = − Σ s = 1 , s ≠ l L λ s l . Letting r ( t ) = s and combining inequality in Equation (2), we get the following TDMJSs: x ̇ ( t ) = A s x ( t ) + B s x ( t − τ ) + C s ω ( t ) + D s u ( t ) + g s ( x ( t ) , x ( t − τ ) , ω ( t ) ) , z ( t ) = E s x ( t ) + F s x ( t − τ ) + G s ω ( t ) , (4)where g s ( x ( t ) , x ( t − τ ) , ω ( t ) ) = f ( x ( t ) , x ( t − τ ) , ω ( t ) ) − [ A s x ( t ) + B s x ( t − τ ) + C s ω ( t ) ] . Then, the following inequality holds: | | g s ( x ( t ) , x ( t − τ ) , ω ( t ) ) | | 2 ≤ | | A a s x ( t ) + A b s x ( t − τ ) + A c s ω ( t ) | | 2 (5)In this study, the HMM-based controller is designed by q ( t ) = K δ ( t ) z ( t ) , (6)where K δ ( t ) ∈ ℜ m × n is the controller gain to be designed and the stochastic jump process δ ( t ) is under the range of O = { 1 , 2 , … , O } . The conditional probability matrix Φ = [ ϕ s v ] is shown as P r = { δ ( t ) = v | r ( t ) = s } = ϕ s v , (7)where Σ v = 1 O ϕ s v = 1 . Then the logarithmic quantiser is defined as N i ( e ) = { N 1 ( e 1 ) ⋯ N i ( e i ) } , (8) N i ( e ) = b i j , 1 1 + χ i b i j < e ≤ 1 1 − χ i b i j 0 , e = 0 − N i ( − e ) , e < 0 (9)where i means the all quantisers, q i means the ith component, e and b i j = ι i b 0 with 0 < ι i < 1 , b 0 > 0 are the input and output of quantisers, respectively. The relationships between ι i and χ i are given by χ i = 1 − ι i 1 + ι i . In addition, − χ i e ≤ N i ( e ) − e ≤ χ i e shows the boundedness of the quantisation error, which can be expressed as N ( e ) − e = Ω i e , Ω i e ∈ [ − χ i , χ i ] . (10)It should be noted that Equation (10) always holds for t. Combining Equations (8) and (10), we get N v ( q v ( t ) ) = ( I + Δ v ( t ) ) q v ( t ) , (11)where Δ v ( t ) = diag { Ω 1 v ( t ) , … , Ω iv ( t ) } with Ω m v ( t ) ∈ [ − χ m v , χ m v ] , m = { 1 , 2 , … , i } . By controller in Equation (6) and quantiser in Equation (11), we can get the following controller: u ( t ) = ( I + Δ v ( t ) ) q v ( t ) . (12)Substituting controller in Equation (12) into TDMJSs in Equation (4), we obtain the following closed-loop TDMJSs: x ̇ ( t ) = ( A s + D s ( I + Δ v ( t ) ) K v E s ) x ( t ) + ( B s + D s ( I + Δ v ( t ) ) K v F s ) x ( t − τ ) + ( C s + D s ( 1 + Δ v ( t ) ) K v G s ) ω ( t ) + g s ( x ( t ) , x ( t − τ ) , ω ( t ) ) , z ( t ) = E s x ( t ) + F s x ( t − τ ) + G s ω ( t ) . (13) Remark 1.The mode manipulation of Markov chain is significant in engineering application of MJSs, but the controller cannot acquire the mode r ( t ) and some inaccuracy will be caused. In this study, we bring δ ( t ) as the controller mode to solve the non-synchronous phenomenon, and Equation (7) indicates the relationships between δ ( t ) and r ( t ) . In controller design, K δ ( t ) and δ ( t ) are only related to Δ v , which can reflect the hide information. Then, the closed-loop TDMJSs in Equation (13) can be regarded as a double random process. The following GCC performance index is introduced to design the controller in Equation (12): J = ∫ 0 T x T ( t ) R 1 x ( t ) + u T ( t ) R 2 u ( t ) d t , (14)where R1 and R2 are the given positive-definite matrices and J < J ∗ holds. J ∗ represents the minimal upper bound of the GCC performance index. The energy supply function of the closed-loop TDMJSs (13) is described by J ( z ( t ) , ω ( t ) , T ) = ∫ 0 T E { S ( z ( t ) , ω ( t ) ) } d t . (15)The supply rate is represented by S ( z ( t ) , ω ( t ) ) = z T ( t ) U z ( t ) + 2 z T ( t ) G ω ( t ) + ω T ( t ) V ω ( t ) . The real matrices U, G, V are known and V = V T , U = U T < 0 with − U = U ̲ T U ̲ holds. Definition 1.Given a time interval [ 0 , T ] [39], positive sacalars a 1 , a 2 with a 2 > a 1 and a weighting matrix S > 0 , the FTB will be guaranteed for the closed-loop TDMJSs in Equation (13) with ( a 1 , a 2 , T , S , d ) , where d ≥ 0 with ∫ 0 T ω T ( t ) ω ( t ) d t ≤ d , if x T ( 0 ) S x ( 0 ) ≤ a 1 ⇒ E { x T ( t ) S x ( t ) < a 2 } , ∀ t ∈ { 0 , T } . (16) Remark 2.From Definition 1, we find that the definitions between the FTB and Lyapunov asymptotic stability are different. Lyapunov asymptotic stability considers the infinite-time interval behaviour of dynamic systems. However, the main concern of FTB is the transient performance in a specified time interval. In a word, FTB is easier to satisfy than Lyapunov asymptotic stability.Definition 2. Under zero initial condition, if the given scalars α > 0 and T > 0 satisfy the following inequality: J ( z ( t ) , ω ( t ) , T ) > α ∫ 0 T ω T ( t ) ω ( t ) d t , (17) the closed-loop TDMJSs in Equation (13) is strictly ( U, G, V)–α-dissipative. Lemma 1.Given two real matrices with suitable dimensions X and Y, there exists a constant ε > 0 and vectors x , y ∈ ℜ n , such that 2 x T X Y y ≤ ε − 1 x T X T X x + ε y T Y T Y y . (18) 3 MAIN RESULTS Here, we set the controller gain as [8, 40] K v = W v H v − 1 , (19)where W v and H v are unknown matrices to be designed. Then, we define three novel auxiliary variables as θ ( t ) = Y s − 1 x ( t ) , (20) θ ( t − τ ) = Y s − 1 x ( t − τ ) , (21) ξ ( t ) = E s θ ( t ) + F s θ ( t − τ ) − H v − 1 z ( t ) , (22)where Y s ∈ ℜ n × n are a set of positive-definite symmetric matrices. Considering Equations (19) to (22), we can get x ̇ ( t ) = ( A s Y s + D s ( I + Δ v ( t ) ) W v E s ) θ ( t ) + ( B s Y s + D s ( I + Δ v ( t ) ) W v F s ) θ ( t − τ ) + C s ω ( t ) − D s ( I + Δ v ( t ) ) W v ξ ( t ) + g s ( x ( t ) , x ( t − τ ) , ω ( t ) ) . (23)Then, we will propose sufficient conditions to ensure the FTB of the closed-loop TDMJSs in Equation (13) and investigate the GCC performance. Theorem 1.Under the given scalars γ s > 0 , the FTB of the closed-loop TDMJSs in Equation (13) with ( a 1 , a 2 , T , S , d ) is guaranteed, and the GCC performance index holds J ∗ = x T ( 0 ) Y s − 1 x ( 0 ) + d , if for any s ∈ L and v ∈ O , there exists a set of mode-dependent scalars ρ s v , σ s > 0 and Y s > 0 satisfying the following LMIs: Ψ < 0 , (24) Σ < 0 , (25) S < Y s − 1 < σ s S , (26) e γ s T σ s a 1 + d γ s ( 1 − e γ s T ) < a 2 , (27) where Ψ = N 1 N 2 N 3 , Σ = Z 1 Z 2 Z 3 , N 1 = Z 1 + Z 2 − γ s Y s Z 3 Z 4 − Y s Z 5 ∗ Z 6 , N 2 = C s Y s A a s T Q 1 0 Y s A b s T 0 − ρ s v G s 0 0 , N 3 = − I A c s T 0 − ε − 1 I 0 ∗ Q 2 , Z 1 = Z 1 + Z 2 Z 3 Z 4 C s − Y s Z 5 0 ∗ Z 6 − ρ s v G s 0 0 − I , Z 2 = Q 3 Y s Y s A a s T Q 1 Q 4 0 Y s A b s T 0 Q 5 0 0 0 0 0 A c s T 0 , Z 3 = − R 2 − 1 0 0 0 0 − R 1 − 1 0 0 0 0 − ε − 1 0 0 0 0 Q 2 , Z 1 = ( λ s s + 1 ) Y s + ε − 1 , Z 2 = H e r ( Σ v = 1 O ϕ s v ( A s Y s + D s ( I + Δ v ( t ) ) W v E s ) ) , Z 3 = B s Y s + Σ v = 1 O ϕ s v ( A s Y s + D s ( I + Δ v ( t ) ) W v F s ) , Z 4 = ρ s v ( E s T H v T − Y s E s T ) − Σ v = 1 O ϕ s v ( D s ( I + Δ v ( t ) ) W v ) , Z 5 = ρ s v ( F s T H v T − Y s F s T ) , Z 6 = − ρ s v ( H e r ( H v ) ) , Q 1 = Y s [ λ s 1 , … λ s s − 1 , λ s s + 1 , … λ s L ] , Q 2 = − diag { Y 1 ⋯ Y s − 1 , Y s + 1 ⋯ Y L } , Q 3 = E s T W v T ( I + Δ v ( t ) ) T , Q 4 = F s T W v T ( I + Δ v ( t ) ) T , Q 5 = − W v T ( I + Δ v ( t ) ) T . Proof.We construct the stochastic Lyapunov functional candidate as V ( x ( t ) ) = x T ( t ) Y s − 1 x ( t ) + ∫ − τ 0 x T ( t + T ) Y s − 1 x ( t + T ) d T . (28) We define Λ V as the weak infinitesimal generator, and get Λ V ( x ( t ) ) = x T ( t ) ( Σ l = 1 L λ s l Y l − 1 ) x ( t ) + 2 Σ v = 1 O ϕ s v x ̇ T ( t ) Y s − 1 x ( t ) + x T ( t ) Y s − 1 x ( t ) − x T ( t − τ ) Y s − 1 x ( t − τ ) . (29)Considering Equations (19) to (22), Equation (29) can be written as Λ V ( x ( t ) ) = θ T ( t ) Y s ( Σ l = 1 O λ s l Y l − 1 ) Y s θ ( t ) + 2 Σ v = 1 O ϕ s v [ ( A s Y s + D s ( I + Δ v ( t ) ) W v E s ) θ ( t ) + ( B s Y s + D s ( I + Δ v ( t ) ) W v F s ) θ ( t − τ ) + C s ω ( t ) − D s ( I + Δ v ( t ) ) W v ξ ( t ) + g s ( x ( t ) , x ( t − τ ) , ω ( t ) ) ] T θ ( t ) + x ( t ) T Y s − 1 x ( t ) − x T ( t − τ ) Y s − 1 x ( t − τ ) , (30)and the following relationship holds: 0 = H v E s θ ( t ) + H v F s θ ( t − τ ) − H v ξ ( t ) − z ( t ) = ( H v E s − E s Y s ) θ ( t ) + ( H v F s − F s Y s ) θ ( t − τ ) − H v ξ ( t ) − G s ω ( t ) . (31)By Lemma 1 and inequality in Equation (5), we get 2 g s T ( x ( t ) , x ( t − τ ) , ω ( t ) ) θ ( t ) ≤ ε − 1 θ T ( t ) θ ( t ) + ε g s T ( x ( t ) , x ( t − τ ) , ω ( t ) ) g s ( x ( t ) , x ( t − τ ) , ω ( t ) ) ≤ ε − 1 θ T ( t ) θ ( t ) + ε [ A a s Y s θ ( t ) + A b s Y s θ ( t − τ ) + A c s ω ( t ) ] T [ A a s Y s θ ( t ) + A b s Y s θ ( t − τ ) + A c s ω ( t ) ] . (32)For the given scalars γ s > 0 , we define J 1 ( t ) = E { Λ V ( x ( t ) ) − γ s V ( x ( t ) ) − ω T ( t ) ω ( t ) } . (33)By Equations (30) to (33), the following relationship holds: J 1 ( t ) = E { Λ V ( x ( t ) ) − γ s V ( x ( t ) ) − ω T ( t ) ω ( t ) } + 2 ρ s v ξ T ( t ) [ ( H v E s − E s Y s ) θ ( t ) + ( H v F s − F s Y s ) θ ( t − τ ) − H v ξ ( t ) − G s ω ( t ) ] ≤ η T ( t ) Ψ 1 η ( t ) + ε [ A a s Y s θ ( t ) + A b s Y s θ ( t − τ ) + A c s ω ( t ) ] T [ A a s Y s θ ( t ) + A b s Y s θ ( t − τ ) + A c s ω ( t ) ] , (34)where η T ( t ) = [ θ T ( t ) θ T ( t − τ ) ξ T ( t ) ω T ( t ) ] , Ψ 1 = Z 7 − Z 2 − γ s Y s Z 3 Z 4 C s ∗ − Y s Z 5 0 ∗ ∗ Z 6 − ρ s v G s ∗ ∗ ∗ − I , Z 7 = Y s ( Σ l = 1 L λ s l Y l − 1 ) Y s + Y s + ε − 1 . By using Schur complement for inequality in Equation (34) and combining inequality in Equation (24), we can ensure J 1 ( t ) < 0 . Then, we have E { Λ V ( x ( t ) ) } < γ s V ( x ( t ) ) + ω T ( t ) ω ( t ) . (35)Multiplying inequality in Equation (35) by e − γ s t and taking integration from 0 to t, we obtain e − γ s t E { V ( x ( t ) ) } − E { V ( 0 ) } < ∫ 0 t e − γ s t ω T ( t ) ω ( t ) d t . (36)By γ s > 0 and t ∈ [ 0 , T ] , we get E { V ( x ( t ) ) } < e γ s t E { V ( 0 ) } + e γ s t d ∫ 0 t e − γ s t d t < e γ s t x T ( 0 ) Y s − 1 x ( 0 ) + d γ s ( 1 − e − γ s t ) ≤ e γ s T x T ( 0 ) Y s − 1 x ( 0 ) + d γ s ( 1 − e − γ s T ) . (37)Then, the following relationship holds: E { x T ( t ) S x ( t ) } < e γ s T ε m a x ( S − 1 2 Y s − 1 S − 1 2 x T ( 0 ) S x ( 0 ) ) + d γ s ( 1 − e − γ s T ) ε m i n ( S − 1 2 Y s − 1 S − 1 2 ) . (38)From inequalities in Equations (26) and (27), we can ensure ε m a x ( S − 1 2 Y s − 1 S − 1 2 ) < σ s and ε m i n ( S − 1 2 Y s − 1 S − 1 2 ) > 1 , which gives E { x T ( t ) S x ( t ) < a 2 } . From Definition 1, the FTB of the closed-loop TDMJSs in Equation (13) is guaranteed. Then, we will investigate the GCC performance and define J 2 ( t ) = Λ V ( x ( t ) ) + x T ( t ) R 1 x ( t ) + u T ( t ) R 2 u ( t ) − ω T ( t ) ω ( t ) . (39)Substituting Equation (12) into Equation (39) and combining Equations (30) to (32), we obtain J 2 ( t ) < 0 by inequality in Equation (25). Then, taking integration it from 0 to T and recalling to Definition 1, we get J < J ∗ = x T ( 0 ) Y s − 1 x ( 0 ) + d . (40)The proof is completed. In the next theorem, we will ensure the strict dissipativity of the closed-loop TDMJSs in Equation (13). Theorem 2.The FTB and strict dissipativity of the closed-loop TDMJSs in Equation (13) is guaranteed, and the controller meets the GCC performance, if for any s ∈ L and v ∈ O , there exists a set of mode-dependent scalars ρ s v , Y s > 0 satisfying Equations (24) to (27) and the following matrix inequality: Ξ < 0 , (41) where Ξ = Γ 1 Γ 2 ∗ Γ 3 , Γ 1 = Z 1 + Z 2 Z 3 Z 4 ∗ − Y s Z 5 ∗ ∗ Z 6 , Γ 2 = Z 8 Y s A a s T Q 6 Q 1 Z 9 Y s A b s T Q 7 0 − ρ s v G s 0 0 0 , Γ 3 = Z 10 A c s T Q 8 0 ∗ − ε − 1 I 0 0 ∗ ∗ − I 0 ∗ ∗ ∗ Q 2 , Z 8 = C s − Y s E s T G , Z 9 = − Y s F s T G , Z 10 = α I − V − H e r ( G s T G ) , Q 6 = { ϕ s 1 Y s E s T U , ϕ s 2 Y s E s T U , … , ϕ s o Y s E s T U } , Q 7 = { ϕ s 1 Y s F s T U , ϕ s 2 Y s F s T U , … , ϕ s o Y s F s T U } , Q 8 = { ϕ s 1 G s T U , ϕ s 2 G s T U , … , ϕ s o G s T U } . Proof.We define J 3 ( t ) = E { Λ V ( x ( t ) ) } − S ( z ( t ) , ω ( t ) ) + α ω T ( t ) ω ( t ) . (42) Combining Equations (20) to (22), the supply rate can be written as S ( z ( t ) , ω ( t ) ) = [ E s Y s θ ( t ) + F s Y s θ ( t − τ ) + G s ω ( t ) ] T × U [ E s Y s θ ( t ) + F s Y s θ ( t − τ ) + G s ω ( t ) ] + 2 [ E s Y s θ ( t ) + F s Y s θ ( t − τ ) + G s ω ( t ) ] T G ω ( t ) + ω T ( t ) V ω ( t ) . (43)From inequalities in Equations (30) to (32) and Equations (42) and (43), we have J 3 ( t ) = E { Λ V ( x ( t ) ) } − S ( z ( t ) , ω ( t ) ) + α ω T ( t ) ω ( t ) + 2 ρ s v ξ T ( t ) [ ( H v E s − E s Y s ) θ ( t ) + ( H v F s − F s Y s ) θ ( t − τ ) − H v ξ ( t ) − G s ω ( t ) ] ≤ η T ( t ) Ξ 1 η ( t ) + ε [ A a s Y s θ ( t ) + A b s Y s θ ( t − τ ) + A c s ω ( t ) ] T [ A a s Y s θ ( t ) + A b s Y s θ ( t − τ ) + A c s ω ( t ) ] , (44)where Ξ 1 = Z 7 + Z 2 − Z 11 Z 3 − Z 12 Z 4 Z 8 − Z 13 ∗ − Y s − Z 14 Z 5 Z 9 − Z 15 ∗ ∗ Z 6 − ρ s v G s ∗ ∗ ∗ Z 10 − Z 16 , Z 11 = Y s E s T U E s Y s , Z 12 = Y s E s T U F s Y s , Z 13 = Y s E s T U G s , Z 14 = Y s F s T U F s Y s , Z 15 = Y s F s T U G s , Z 16 = G s T U G s . By using Schur complement for inequality in Equation (44) and combining inequality in Equation (41), we can ensure J 3 ( t ) < 0 . Then, integrating it with zero initial conditions, we get E { V ( x ( t ) ) − ∫ 0 T S ( z ( t ) , ω ( t ) ) d t + ∫ 0 T α ω T ( t ) ω ( t ) d t } < 0 . (45)Considering V ( x ( t ) ) > 0 , we obtain ∫ 0 T E { S ( z ( t ) , ω ( t ) ) } d t > α ∫ 0 T ω T ( t ) ω ( t ) d t . (46)By inequality in Equation (17), the FTB and the strict dissipativiy of the closed-loop TDMJSs in Equation (13) are guaranteed, and the controller meets the GCC performance. The proof is completed. Remark 3. In order to avoid the existence of the Y s and Y s − 1 in the matrix inequalities at the same time, we define two novel auxiliary variables in Equations (20) and (21). Then, in order to guarantee that LMIs in Equations (24), (25) and (41) are established, we need to make the principal diagonal be negative. We define the novel auxiliary variable in Equation (22). For the unsolvable form Y s ( Σ l = 1 L λ s l Y l − 1 ) Y s , we use Schur complement to transform it into Q 1 and Q 2 shown in LMIs in Equations (24), (25) and (41). 4 SIMULATION EXPERIMENTS In this section, we will introduce the liquid monopropellant rocket motor model with a pressure-feeding system [41, 42] to demonstrate the correctness and feasibility of the designed approach shown as x 1 ( t ) = ( β s − 1 ) x 1 − β s x 1 ( t − τ ) + x 3 ( t − τ ) , x 2 ( t ) = 1 π ϝ 1 s [ − x 4 ( t ) + u ( t ) + ω 2 ( t ) ] , x 3 ( t ) = 1 ( 1 − π ) ϝ 1 s [ − x 3 ( t ) + x 4 ( t ) − ϝ 2 s x 1 ( t ) ] , x 4 ( t ) = 1 ϝ
Referência(s)