Disturbance estimation based tracking control for periodic piecewise time‐varying delay systems
2020; Institution of Engineering and Technology; Volume: 15; Issue: 3 Linguagem: Inglês
10.1049/cth2.12055
ISSN1751-8652
AutoresR. Sakthivel, S. Harshavarthini, Nasser‐eddine Tatar,
Tópico(s)Iterative Learning Control Systems
ResumoIET Control Theory & ApplicationsVolume 15, Issue 3 p. 459-471 ORIGINAL RESEARCH PAPEROpen Access Disturbance estimation based tracking control for periodic piecewise time-varying delay systems Rathinasamy Sakthivel, Corresponding Author krsakthivel@yahoo.com Department of Applied Mathematics, Bharathiar University, Coimbatore, India Correspondence Rathinasamy Sakthivel, Department of Applied Mathematics, Bharathiar University, Coimbatore–641046, India. Email: krsakthivel@yahoo.comSearch for more papers by this authorShanmugam Harshavarthini, Department of Applied Mathematics, Bharathiar University, Coimbatore, IndiaSearch for more papers by this authorNasser-Eddine Tatar, Department of Mathematics and Statistics, King Fahd University of Petroleum and Minerals, Dhahran, Saudi ArabiaSearch for more papers by this author Rathinasamy Sakthivel, Corresponding Author krsakthivel@yahoo.com Department of Applied Mathematics, Bharathiar University, Coimbatore, India Correspondence Rathinasamy Sakthivel, Department of Applied Mathematics, Bharathiar University, Coimbatore–641046, India. Email: krsakthivel@yahoo.comSearch for more papers by this authorShanmugam Harshavarthini, Department of Applied Mathematics, Bharathiar University, Coimbatore, IndiaSearch for more papers by this authorNasser-Eddine Tatar, Department of Mathematics and Statistics, King Fahd University of Petroleum and Minerals, Dhahran, Saudi ArabiaSearch for more papers by this author First published: 16 December 2020 https://doi.org/10.1049/cth2.12055AboutSectionsPDF 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 onEmailFacebookTwitterLinked InRedditWechat Abstract This work is mainly focused on the design of uncertainty and disturbance estimator-based tracking control scheme for a class of periodic piecewise time-varying delay systems with uncertainties and disturbances. Precisely, the considered system consists of several subsystems, where the fixed dwell time and switching sequence are utilised to achieve the coordination between each subsystem. By constructing Lyapunov–Krasovskii functional with time-varying periodic piecewise matrices and employing the matrix polynomial approach, an uncertainty and disturbance estimator-based periodic tracking control design is developed such that the trajectories of the considered system are driven to track the reference model. Moreover, the lumped disturbances so called the summation of uncertainties and exogenous disturbances are estimated with the aid of strictly proper low-pass filter. Finally, a numerical example is presented to verify the effectiveness and advantages of the developed theoretical results. 1 INTRODUCTION Periodicity is a common factor encountered in practice, wherein the certain characteristics are repeated in regular period of time. More specifically, the periodic behaviour occurs in many engineering problems such as helicopter rotors, DC–DC converters, transmission lines with distinct loads, vehicle suspension, rotor bearing systems, and power electronic equipment's with multi-function converters. Due to its significance and widespread applications, the stability analysis of periodic systems has received much attention among many research communities [1-3]. Periodic systems with finite number of subsystems in a unit period of time are referred to as periodic piecewise systems (PPSs) [4, 5]. Generally, the PPSs are composed of linear time-invariant subsystems and hence it may lose dynamic characteristics of the original periodic systems. Thus, the study of PPSs with time-varying subsystems attains considerable attention among researchers [6-10]. Notably in [6], the authors developed a state feedback control law with time-varying gain matrices to ensure the exponential stability of PPSs with time-varying subsystems. Furthermore in [7], a novel matrix polynomial lemma is proposed to linearize the non-linearities raised in between the time-varying matrices, wherein the system and controller is constituted with distinct width of time-varying coefficients which bring more flexibility to convex optimisation tools. In the work of Zhu et al. [9], the exponential stability of the periodic piecewise time-varying systems (PPTVSs) is ensured by Lyapunov stability theory. Li et al. [10] proposed an H ∞ control with time-varying polynomial controller gain matrix that ensures the exponential stability of PPTVSs with polynomial subsystems. In practice, the occurrence of delay is inexorable and it is the source of instability or performance degradation of various kinds of dynamical control systems [11-13]. Therefore, it is necessary and important to deal with the complications resulting from the combination of time-varying delay and periodicity in various practical systems. Considerable efforts devoted to the study of stability and control problems of delayed PPSs are reported in the literature [14-16]. The authors in [14] proposed an H ∞ control scheme for a class of PPSs with time-invariant subsystems, disturbances and time-varying delay, which ensures the exponential stability of the system. Xie and Lam [15] developed a guaranteed cost periodic control scheme for delayed PPSs with time-invariant subsystems. In [16], the non-fragile controller with additive and multiplicative gain perturbations is designed to stabilise the closed-loop delayed PPSs with time-varying subsystems, but it is limited to constant delays. However, to the best of the present author's knowledge, no work on the combination of PPSs with time-varying subsystems and time-varying delays has been reported in the literature. It should be noted that the exact mathematical model of the system is primarily required to achieve the desired system performance, but in practice, the exact modelling of physical systems is not always guaranteed because of the differences or errors between the mathematical model and the real physical systems, which are commonly referred to as model uncertainties. Moreover, the existence of uncertainties is categorised by several factors such as modelling errors, variation in real parameters, aging of system components, temperature fluctuations, and operating point shifting. Furthermore, unexpected errors in modelling the real-world systems happen frequently. Due to the limited knowledge of physical phenomena, uncertainties are regarded as unknown dynamics which are frequently a source of instability and very difficult to design a robust control law. Therefore, it is necessary and important to formulate the system with model uncertainties. A great number of fruitful works on the stability and stabilisation analysis of uncertain dynamical systems have been discussed in [17-20]. On the other hand, the existence of abrupt disturbances are unpredictable and may severely affect the dynamics of the system [21, 22]. Because of the combination of uncertainties and disturbances, much effort on the control design is needed for the stabilisation of system trajectories. In this regard, it is more significant to estimate the uncertainties, non-linearities and exogenous disturbances in control systems. At the early stage, several control approaches were proposed for various dynamical systems to compensate the disturbance effects. Most of the control approaches require some assumptions about the unpredictable aspects like bound values of unknown uncertainties and disturbances, but in practice they are not always measurable. In order to conquer the above complexity, an effective technique called UDE-based strategy was proposed in [23] and it is widely used for many dynamical systems [23-29]. The UDE-based control approach proposed by Zhong and Rees [23] ensures better tracking performance of the linear time-invariant systems and also estimates the lumped disturbances with respect to the appropriate bandwidth selection of the low-pass filter. In unpredictable environments, the desired behaviour of the system trajectories is one of the key objectives to maintain better performance of the control systems. The tracking control problem is one of the most eminent and essential research topics which has enormous applications in science and engineering [30-35]. Owing to these perspectives, the authors in [31] developed an H ∞ tracking control scheme with periodic control gain parameters for a class of PPSs with time-invariant subsystems, wherein the desired tracking performance of the system states to the reference state is attained within a finite-period of time. Unfortunately, only one published paper addressed the problem in designing the tracking control for PPSs with time-invariant subsystems. To the best of the present author's knowledge, the tracking control problem of PPSs with time-varying subsystems has not yet been investigated. In order to fulfill such research gap and also inspired by the significance of the works in [6-10, 23 31, 36, 37], a disturbance and uncertainty estimation based tracking control is designed for a class of uncertain time-varying delayed PPTVSs with disturbances. To be precise, the PPTV control gain matrix is obtained by solving a set of delay-dependent stability conditions, which is derived by considering periodically time-varying Lyapunov–Krasovskii functional. Especially, the matrix polynomial approach proposed in [7] is utilised to linearize the non-linearities raised in between the time-varying matrices. Further, the main contributions of this work are listed as follows: 1. Tracking control design with periodic behaviour is developed for a class of uncertain PPTVSs with time-varying delay and disturbances for the first time. 2. Inspired by the work in [31], a new periodic piecewise time-varying reference model with time-varying delay and bounded input is considered for the addressed system. 3. To tackle the impact of lumped disturbances, that is the total sum of uncertainties and unexpected disturbances, the low-pass filter is implemented with appropriate bandwidth value and to facilitate the complexities in derivations. Moreover, here the general assumptions on lumped disturbances are relaxed together with less tuning parameters. 4. The proposed periodic control design is more effective for state tracking, which is endorsed by presenting simulation results. 2 PROBLEM FORMULATION AND PRELIMINARIES Consider a class of uncertain linear PPTVSs with time-varying state delay and disturbances in the following compact form: ψ ̇ ( t ) = ( A ( t ) + δ A ( t ) ) ψ ( t ) + ( A γ ( t ) + δ A γ ( t ) ) ψ ( t − γ ( t ) ) + B ( t ) u ( t ) + B ω ( t ) ω ( t ) , (1)where ψ ( t ) ∈ R n ψ , u ( t ) ∈ R m and ω ( t ) ∈ L 2 [ 0 , ∞ ) , respectively, represent the state, control and disturbance vectors; A ( t ) , A γ ( t ) , B ( t ) and B ω ( t ) are known time-varying appropriate dimensioned periodic matrices with the fundamental period T p , that is, A ( t ) = A ( t + ρ T p ) , A γ ( t ) = A γ ( t + ρ T p ) , B ( t ) = B ( t + ρ T p ) and B ω ( t ) = B ω ( t + ρ T p ) , ∀ t > 0 , ρ = 0 , 1 , 2 , … ; δ A ( t ) and δ A γ ( t ) represent the unknown time-varying uncertainties of the considered system (1); γ ( t ) denotes the delay in system dynamics which satisfies the conditions 0 ≤ γ ( t ) ≤ γ ¯ and γ ̇ ( t ) ≤ γ ̂ , wherein the parameters γ ¯ > 0 and γ ̂ < 1 are known scalars. Now, the each periodic time interval [ ρ T p , ( ρ + 1 ) T p ) is partitioned into S sub-intervals [ ρ T p + t i − 1 , ρ T p + t i ) , i ∈ S = { 1 , 2 , … , S } . In this connection, for t ∈ T i = Δ [ ρ T p + t i − 1 , ρ T p + t i ) , the time-varying system parameters of the ith subsystem with dwell time T i = t i − t i − 1 are denoted by A i ( t ) = A i + ξ i ( t ) ( A i + 1 − A i ) ; δ A i ( t ) = δ A ¯ i ( t ) + ξ i ( t ) ( δ A ¯ i + 1 ( t ) − δ A ¯ i ( t ) ; A γ i ( t ) = A γ i + ξ i ( t ) ( A γ ( i + 1 ) − A γ i ) ; δ A γ i ( t ) = δ A ¯ γ i ( t ) + ξ i ( t ) ( δ A ¯ γ ( i + 1 ) ( t ) − δ A ¯ γ i ( t ) ) ; B i ( t ) = B i + ξ i ( t ) ( B i + 1 − B i ) ; B ω i ( t ) = B ω i + ξ i ( t ) ( B ω ( i + 1 ) − B ω i ) , where ξ i ( t ) = t − ( ρ T p + t i − 1 ) T i , A i , A γ i , B i , B ω i , i ∈ S are appropriate dimensioned known constant matrices and δ A ¯ i ( t ) , δ A ¯ γ i ( t ) represent the uncertainties of system dynamics over each sub-interval with the following structure: δ A ¯ i ( t ) δ A ¯ γ i ( t ) = H i Γ i ( t ) E A i E B i , (2)where H i , E A i and E B i are appropriate dimensioned matrices, and Γ i T ( t ) Γ i ( t ) ≤ I . The system (1) is then reformulated as follows: ψ ̇ ( t ) = ( A i ( t ) + δ A i ( t ) ) ψ ( t ) + ( A γ i ( t ) + δ A γ i ( t ) ) ψ ( t − γ ( t ) ) + B i ( t ) u ( t ) + B ω i ( t ) ω ( t ) , ∀ t ∈ [ ρ T p + t i − 1 , ρ T p + t i ) . (3) Assumption 1.The periodic time-varying coefficient matrix B i ( t ) of control input vector is assumed to be full column rank matrix. Remark 1.Inspired by the work in [4], the system uncertainties are described in periodic piecewise form δ A i ( t ) = δ A ¯ i ( t ) + ξ i ( t ) ( δ A ¯ i + 1 ( t ) − δ A ¯ i ( t ) and δ A γ i ( t ) = δ A ¯ γ i ( t ) + ξ i ( t ) × ( δ A ¯ γ ( i + 1 ) ( t ) − δ A ¯ γ i ( t ) ) . Further, the structure defined in (2) consists of partially known matrices H i ( t ) , E A i and E B i . Unlike the conventional description of uncertainties, in this work, the coupling of time-varying term ξ i ( t ) and unknown matrix Γ i ( t ) exist in the uncertainties. It ends up as a non-convex variable that is difficult to decouple in stability analysis, which brings more technical challenges. Moreover, this non-convex variable affects the fixed dwell time switching. Remark 2.The addressed PPTVS is formulated by the inherent dynamic characteristics of practical time-varying systems with periodically varying elements. In particular, the system (1) is partitioned into finite number of sub-intervals with fixed dwell time T i . The selection of T i , i = 1 , 2 , … , S is based on the variation of the system dynamics. Subsequently, the corresponding subsystem over each sub-interval is approximated by an averaged model [2]. Based on this concept, the continuous-time periodic system (1) is remodelled as PPTVS in (3) with PPTV matrix coefficients. In order to study the state tracking control problem, for all t ∈ [ ρ T p + t i − 1 , ρ T p + t i ) , let us consider the following PPTV reference model: ψ ̇ r ( t ) = A r i ( t ) ψ r ( t ) + A r γ i ( t ) ψ r ( t − γ ( t ) ) + B r i ( t ) ω r ( t ) , (4)where ψ r ( t ) and ω r ( t ) are the reference state and uniformly bounded input vector, respectively; A r i ( t ) = A r i + ξ i ( t ) ( A r ( i + 1 ) − A r i ) , A r γ i ( t ) = A r γ i + ξ i ( t ) ( A r γ ( i + 1 ) − A r γ i ) and B r i ( t ) = B r i + ξ i ( t ) ( B r ( i + 1 ) − B r i ) are the given appropriate dimensioned PPTV matrices with the fundamental period T p . The main purpose of this work is to develop a robust tracking control u ( t ) , ensuring the asymptotic convergence of the error vector ε ( t ) = ψ r ( t ) − ψ ( t ) , that is, the states of the system (3) asymptotically track the states of the reference model (4). In other words, the error dynamics model ε ̇ ( t ) = ( A r i ( t ) + K i ( t ) ) ε ( t ) + A r γ i ( t ) ε ( t − γ ( t ) ) (5)should be asymptotically stable, where K i ( t ) is the appropriate dimensioned error feedback PPTV gain matrix to be determined. Combining the equations (3)-(5), we obtain u ( t ) = B i + ( t ) [ A r i ( t ) ψ ( t ) + A r γ i ( t ) ψ ( t − γ ( t ) ) + B r i ( t ) ω r ( t ) − K i ( t ) ε ( t ) − A i ( t ) ψ ( t ) − A γ i ( t ) ψ ( t − γ ( t ) ) − u ω ( t ) ] , (6)where B i + ( t ) = ( B i T ( t ) B i ( t ) ) − 1 B i T ( t ) is the pseudo inverse of B i ( t ) and u ω ( t ) = δ A i ( t ) ψ ( t ) + δ A γ i ( t ) ψ ( t − γ ( t ) ) + B ω i ( t ) × ω ( t ) is the collection of both undesirable uncertainties and unknown disturbances in the addressed PPTVS (3). It should be noted that the control law defined in the equation (6) is not readily applicable because of the presence of unknown terms of the system (3). Thus, the unknown lumped disturbance u ω ( t ) is observed with the aid of the dynamics of system (3): u ω ( t ) = ψ ̇ ( t ) − [ A i ( t ) ψ ( t ) + A γ i ( t ) ψ ( t − γ ( t ) ) + B i ( t ) u ( t ) ] . (7) However, the state derivative ψ ̇ ( t ) in the relation (7) will cause algebraic loops. To overcome this issue, the UDE-based control strategy in [23] is adopted to estimate the lumped disturbance (7). Precisely, u ω ( t ) is estimated by using the low-pass filter with appropriate bandwidth. In this connection, a low-pass filter Q g ( s ) = 1 1 + f s is utilised to estimate the unknown lumped disturbance with respect to the appropriate bandwidth f of the filter. Now, the lumped disturbance u ω ( t ) can be accurately estimated by passing it into the developed low-pass filter Q g ( s ) . Then, the estimated lumped disturbance is given by u ω e ( t ) = u ω ( t ) ∗ q g ( t ) , (8)where u ω ( t ) = ψ ̇ ( t ) − ( A i ( t ) ψ ( t ) + A γ i ( t ) ψ ( t − γ ( t ) ) + B i ( t ) × u ( t ) ) , ` ` ∗ ′ ′ denotes the convolution operator and q g ( t ) represents the impulse response of the filter Q g ( s ) . Then, by substituting the estimated lumped disturbance u ω e ( t ) in (6), we get the UDE-based control law as follows: u ( t ) = B i + ( t ) [ A r i ( t ) ψ ( t ) + A r γ i ( t ) ψ ( t − γ ( t ) ) + B r i ( t ) ω r ( t ) − K i ( t ) ε ( t ) − A i ( t ) ψ ( t ) − A γ i ( t ) ψ ( t − γ ( t ) ) − [ ψ ̇ ( t ) − ( A i ( t ) ψ ( t ) + A γ i ( t ) ψ ( t − γ ( t ) ) + B i ( t ) u ( t ) ] ∗ q g ( t ) ) ] ] . (9)By employing the linear property of the convolution operator and inverse Laplace transform, the UDE-based control law in the above equation is reformulated as u ( t ) = B i + ( t ) [ L − 1 { 1 1 − Q g ( s ) } ∗ [ A r i ( t ) ψ ( t ) + A r γ i ( t ) × ψ ( t − γ ( t ) ) + B r i ( t ) ω r ( t ) − K i ( t ) ε ( t ) ] − A i ( t ) ψ ( t ) − A γ i ( t ) ψ ( t − γ ( t ) ) − L − 1 { s Q g ( s ) 1 − Q g ( s ) } ∗ ψ ( t ) ] , (10)where L − 1 { . } denotes the inverse Laplace operator. Now, by substituting the low-pass filter Q g ( s ) in the above equation and by simple calculation, we obtain u ( t ) = B i + ( t ) { − A i ( t ) ψ ( t ) − A γ i ( t ) ψ ( t − γ ( t ) ) + A r i ( t ) ψ r ( t ) + A r γ i ( t ) ψ r ( t − γ ( t ) ) + B r i ( t ) ω r ( t ) − A r i ( t ) ε ( t ) + A r γ i ( t ) ε ( t − γ ( t ) ) − K i ( t ) ε ( t ) + 1 f [ ε ( t ) − ∫ 0 t [ ( A r i ( s ) + K i ( s ) ) ε ( s ) + A r γ i ( s ) ε ( s − γ ( s ) ) ] d s ] } . (11) Remark 3.It should be noted that the lumped disturbance may degrade the system performance. The proposed UDE-based controller effectively estimates the lumped disturbance and ensures better tracking performance of the considered PPTVS. More specifically, the significant application of the PPTVSs in many real-world problems and the necessity of disturbance rejection motivates this work to design the UDE-based tracking control scheme for PPTVSs. The proposed UDE-based control scheme consists of three distinct components, namely, reference model, filter design and error feedback gain. In particular, the reference model is used to determine the set-point response. The filter approach is utilised to estimate the unpredictable factors of the addressed system and the disturbance response is determined by the error feedback gain. Owing to this perspective, the developed control law (11) paves a way to obtain a better tracking performance between the states of the actual system and those of the reference model. Besides, the following lemmas are required to derive the desired stability constraints. Lemma 1. ([[7]])The matrix polynomial function ϖ : [ 0 , 1 ] 2 → R is defined by ϖ ( ξ ¯ 1 , ξ ¯ 2 ) = Ξ 0 + ξ ¯ 1 Ξ 1 + ξ ¯ 2 Ξ 2 + ξ ¯ 1 ξ ¯ 2 Ξ 3 , where ξ ¯ 1 , ξ ¯ 2 ∈ [ 0 , 1 ] and Ξ 0 , Ξ 1 , Ξ 2 , Ξ 3 are real-valued symmetric matrices. If Ξ 0 < 0 , Ξ 0 + Ξ 1 < 0 and Ξ 0 + Ξ 2 < 0 , Ξ 0 + Ξ 1 + Ξ 2 + Ξ 3 < 0 , then the matrix polynomial ϖ ( ξ ¯ 1 , ξ ¯ 2 ) < 0 . Lemma 2. ([[13]])The inequality ∫ u v θ T ( s ) Υ θ ( s ) d s ≥ 1 v − u ς T R ς holds for all continuously differentiable function θ : [ u , v ] → R n and positive definite matrix Υ, where ς = ∫ u v θ ( s ) d s ∫ u v ∫ s v θ ( r ) d r d s and R = 4 Υ − 6 v − u Υ ∗ 12 ( v − u ) 2 Υ . 3 MAIN RESULT This section will focus on deriving sufficient conditions in the form of linear matrix inequalities ensuring the tracking performance between the states of the addressed system (1) and the desired reference model (4). In other words, the asymptotic stabilisation of the PPTV error system (5) with the aid of the developed UDE-based periodic tracking control (11) is ensured. More specifically, Lemma 1 (matrix polynomial approach) is incorporated for linearising the non-linearities raised in between the periodic time-varying matrices. Before establishing the stability conditions, we consider the continuous time-varying positive-definite matrix P ( t ) with the periodic nature of period T p , that is, P ( t ) = P ( t + ρ T p ) , ∀ t > 0 , ρ = 0 , 1 , 2 , ⋯ . Further, the each subsystem belonging to the sub-interval [ ρ T p + t i − 1 , ρ T p + t i ) is partitioned into M i ∈ N segments of length m i = T i M i . Based on the above consideration, the positive-definite matrix P i ( t ) is defined as follows: P i ( t ) = P i , j + θ i , j ( t ) ( P i , j + 1 − P i , j ) , (12)where θ i , j ( t ) = t − ( ρ T p + t i − 1 + j m i ) m i ∈ [ 0 , 1 ) and P i , M i = P i + 1 , 0 , P S , M i = P i , 0 , j = 0 , 1 , … , M i − 1 . Moreover, the schematic representation of PPTV positive-definite matrix P ( t ) is given in Figure 1. Theorem 1.The PPTV error system (5) with Assumption 1 is considered along with the known scalars γ ¯ > 0 and γ ̂ . If there exist a symmetric positive-definite matrices P i , j , R k and appropriate dimensioned matrices Y i , j , i = 1 , 2 , … , S , j = 1 , 2 , … , M i − 1 , k = 1 , 2 , 3 , such that the following LMI constraints hold: Π i , j , 0 < 0 , (13) Π i , j , 0 + Π i , j , 1 < 0 , (14) Π i , j , 0 + Π i , j , 2 < 0 , (15) Π i , j , 0 + Π i , j , 1 + Π i , j , 2 + Π i , j , 3 < 0 , (16)where Π 11 i , j , 0 = s y m ( P i , j A r i + Y i , j ) + P i , j + 1 − P i , j m i + R 1 + γ ¯ R 2 + γ ¯ 2 R 3 , Π 12 i , j , 0 = P i , j A r γ i , Π 22 i , j , 0 = − ( 1 − γ ̂ ) R 1 , Π 33 i , j , 0 = − γ ¯ R 2 , Π 44 i , j , 0 = − 4 R 3 , Π 45 i , j , 0 = 6 γ ¯ R 3 , Π 55 i , j , 0 = − 12 γ ¯ 2 R 3 , Π 11 i , j , 1 = s y m ( P i , j [ A r ( i + 1 ) − A r i ] ) , Π 12 i , j , 1 = P i , j [ A r γ ( i + 1 ) − A r γ i ] , Π 11 i , j , 2 = s y m ( [ P i , j + 1 − P i , j ] A r i + Y i , j + 1 − Y i , j ) , Π 12 i , j , 2 = [ P i , j + 1 − P i , j ] A r γ i , Π 11 i , j , 3 = s y m ( [ P i , j + 1 − P i , j ] × [ A r ( i + 1 ) − A r i ] and Π 12 i , j , 3 = [ P i , j + 1 − P i , j ] [ A r γ ( i + 1 ) − A r γ i ] , then the system (5) is asymptotically stable. Moreover, to obtain the error feedback PPTV gain matrix, the following relation together with the time-varying matrix functions P i ( t ) in (12) is utilised: K i ( t ) = P i − 1 ( t ) Y i ( t ) , (17) Y i ( t ) = Y i , j + θ i , j ( t ) ( Y i , j + 1 − Y i , j ) . (18) FIGURE 1Open in figure viewerPowerPoint Evolution of PPTV positive-definite matrix P ( t ) Proof.In order to establish the desired delay-dependent stability constraint that ensures the required result, the time-varying Lyapunov–Krasovskii functional candidate for the error system (5) with respect to the PPTV matrix P ( t ) defined in (12) is formulated in the following form: υ ( t ) = ε T ( t ) P ( t ) ε ( t ) + ∫ t − γ ( t ) t ε T ( s ) R 1 ε ( s ) d s + γ ¯ ∫ t − γ ¯ t ε T ( s ) R 2 ε ( s ) d s + γ ¯ ∫ − γ ¯ 0 ∫ t + s t ε T ( u ) R 3 ε ( u ) d u d s , (19)where P i , j , R 1 , R 2 and R3 are the positive-definite matrices.The time derivative of the above considered time-varying Lyapunov–Krasovskii functional along the solution of (5) is computed as υ ̇ ( t ) ≤ 2 ε T ( t ) P i ( t ) ε ̇ ( t ) + ε T ( t ) [ P ̇ i ( t ) + R 1 + γ ¯ R 2 + γ ¯ 2 R 3 ] ε ( t ) − ( 1 − γ ̂ ) ε T ( t − γ ( t ) ) R 1 ε ( t − γ ( t ) ) − γ ¯ ε T ( t − γ ¯ ) R 2 ε ( t − γ ¯ ) − γ ¯ ∫ t − γ ¯ t ε T ( s ) R 3 ε ( s ) d s = ε T ( t ) [ P i ( t ) A r i ( t ) + P i ( t ) K i ( t ) + A r i T ( t ) P i ( t ) + K i T ( t ) P i ( t ) + P ̇ i ( t ) + R 1 + γ ¯ R 2 + γ ¯ 2 R 3 ] ε ( t ) + ε T ( t ) P i ( t ) A r γ i ( t ) ε ( t − γ ( t ) ) + ε T ( t − γ ( t ) ) A r γ i T ( t ) × P i ( t ) ε ( t ) − ( 1 − γ ̂ ) ε T ( t − γ ( t ) ) R 1 ε ( t − γ ( t ) ) − ε T ( t − γ ¯ ) γ ¯ R 2 ε ( t − γ ¯ ) − γ ¯ ∫ t − γ ¯ t ε T ( s ) R 3 ε ( s ) d s . (20)Further, by employing Lemma 2 to the single integral term in (20), we get − γ ¯ ∫ t − γ ¯ t ε T ( s ) R 3 ε ( s ) d s ≤ ψ ∼ T ( t ) − 4 R 3 6 γ ¯ R 3 ∗ − 12 γ ¯ 2 R 3 ψ ∼ ( t ) , (21)where ψ ∼ T ( t ) = ∫ t − γ ¯ t ε T ( s ) d s ∫ − γ ¯ 0 ∫ t + s t ε T ( u ) d u d s .By calculating the time derivative of the matrix function P ( t ) in (12), we obtain P ̇ i ( t ) = P i , j + 1 − P i , j m i , ∀ t > 0 . (22)Now by combining the inequalities (20)–(22), we get υ ̇ ( t ) ≤ ψ ¯ T ( t ) [ Π ] 5 × 5 ψ ¯ ( t ) , (23)where ψ ¯ T ( t ) = ε T ( t ) ε T ( t − γ ( t ) ) ε T ( t − γ ¯ ) ψ ∼ T ( t ) , Π 11 = P i ( t ) A r i ( t ) + Y i ( t ) + A r i T ( t ) P i ( t ) + Y i T ( t ) + P i , j + 1 − P i , j m i + R 1 + γ ¯ R 2 + γ ¯ 2 R 3 , Π 12 = P i ( t ) A r γ i ( t ) , Π 22 = − ( 1 − γ ̂ ) R 1 , Π 33 = − γ ¯ R 2 , Π 44 = − 4 R 3 , Π 45 = 6 γ ¯ R 3 , Π 55 = − 12 γ ¯ 2 R 3 and Y i ( t ) = P i ( t ) K i ( t ) .By simple manipulation, the time-varying piecewise periodic matrices in Π with the fundamental period T p can be expanded in the form of the following matrix polynomial: Π = Π i , j , 0 + ξ i ( t ) Π i , j , 1 + θ i , j ( t ) Π i , j , 2 + ξ i ( t ) θ i , j ( t ) Π i , j , 3 , (24)where ξ i ( t ) , θ i , j ( t ) ∈ [ 0 , 1 ) and the elements of the matrices Π i , j , 0 , Π i , j , 1 , Π i , j , 2 and Π i , j , 3 are the same as in the statement of the theorem. By using Lemma 1 and the conditions (13)–(16) in the equation (24), we obtain the matrix Π < 0 .
Referência(s)