Document not found! Please try again

Robust Finite-Time H-Infinity Control with Transients for Dynamic

0 downloads 0 Views 2MB Size Report
May 16, 2018 - have been proposed, including robust adaptive control [2], sliding mode .... we use max(β‹…)( min(β‹…) ) to denote the maximum (minimum).
Hindawi Mathematical Problems in Engineering Volume 2018, Article ID 2838749, 17 pages https://doi.org/10.1155/2018/2838749

Research Article Robust Finite-Time H-Infinity Control with Transients for Dynamic Positioning Ship Subject to Input Delay Xiaogong Lin, Kun Liang , Heng Li, Yuzhao Jiao , and Jun Nie College of Automation, Harbin Engineering University, Harbin, China Correspondence should be addressed to Kun Liang; [email protected] Received 1 January 2018; Revised 4 April 2018; Accepted 16 May 2018; Published 26 June 2018 Academic Editor: Sabri Arik Copyright Β© 2018 Xiaogong Lin et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. This paper presents the problem of robust finite-time 𝐻∞ control with transients for ocean surface vessels equipped with dynamic positioning (DP) system in presence of input delay. The main objective of this work is to design a finite-time 𝐻∞ state feedback controller, which ensures that all states of ship do not exceed a given threshold over a fixed time interval, with better robustness and transient performance subject to time-varying disturbance. Based on a novel augmented Lyapunov-Krasovskii-like function (LKLF) with triple integral terms and a method combining the Wirtinger inequality and reciprocally convex approach, a less conservative result is derived. In particular, an 𝐻∞ performance index with nonzero initial condition is introduced to attenuate the overconservatism caused by the assumption of zero initial condition and enhance the transient performance of ship subject to external disturbance. More precisely, the controller gain matrix for the DP system can be achieved by solving the linear matrix inequalities (LMIs), which can be easily facilitated by using some standard numerical packages. Finally, a numerical simulation for a ship is proposed to verify the effectiveness and less conservatism of the controller we designed.

1. Introductions With the increasing development in the ocean exploitation, dynamic positioning (DP) systems, regulating the horizontal position and heading of the vessel exclusively by means of active thrusters, have been developed for various marine and offshore applications such as drilling, salvage, pipelaying, and oil production [1]. To achieve expected trajectory tracking or positioning, various control strategies have been proposed, including robust adaptive control [2], sliding mode control [3], prescribed performance control [4], hybrid control [5], and neural network control [6]. However, in some applications, it is significant to maintain the vessel’s states under some bounds, particularly of which transient performance is emphasized, during a specific time interval, for example, when facing the matter of saturations or when the task of trajectory tracking should be fulfilled in a prescribed time interval. In [7], the proposed controller can maintain the bound of ship over an infinity time interval regardless of disturbance. However, it is inconvenient to analyze and to enhance the transient performance when the operation should be arrived in short time, such as rescuing

works. On the other hand, time delay is encountered in many dynamical systems and often leads to performance deterioration which engenders strongly growing interest in this topic in recent decades. In the DP system, the main kind of time delay is encountered in actuators [8], while another obvious kind of delay is the one produced between the sensors and the activation of the control mechanism [9]. The effect caused by time delay will be more significant in finite-time interval and it is the first motivation of this paper. The classical 𝐻∞ control theory [10] defines the control law, under the assumption of zero initial condition, providing the minimal value to the performance measure that is the worst-case norm of the regulated output over all exogenous signals and is applied to DP system successfully [11]. However, there exist some situations, when the initial state of ship is possibly nonzero, such as when saturation is activated or controller is switched to another one. It will cause an additional unknown disturbance [12] and the promising robustness is achieved at the expense of degraded nominal performance [13]. In [14], researchers introduced a performance measure that is the induced norm of the regulated output over all exogenous signals and initial states for finite

2 and infinite horizons. Unfortunately, to the best the authors’ knowledge, it has not been extended to the system with input delay. On another research frontier, various approaches, in the framework of finite-time boundedness for time-delay systems, have been developed to obtain the results with less conservatism [15, 16], usually indicated by the bound of state [17]. For DP system, engineers are willing to obtain control strategies which maintain the state of ship varying in a small region around the desired set-point or tracking path instead of the acceptable minimum state bound. Thus, the second motivate is to obtain a result which reduces the overconservatism caused by the assumption of zero initial condition and the loss information in the proof for time delay systems in a practical view. Based on the discussion aforementioned, the problem of robust finite-time 𝐻∞ control for DP system with input delay is studied. The main contribution of this paper lies in three aspects: Firstly, a finite-time 𝐻∞ controller with transients is designed for DP systems with input delay, by solving a couple of LMIs, which can guarantee the state of ship within a desired value over a fixed time interval in the presence of time-varying disturbances. Secondly, the concept of 𝐻∞ control under nonzero initial condition is introduced to time-delay system and its advantage on enhancing the transient performance will be shown in a numerical simulation compared to the 𝐻∞ with zero initial condition. In particular, the results are established in forms of LMI, which can be easily facilitated by using some standard numerical packages. Thirdly, a novel augmented LKLF with triple integral terms, which contains more information is constructed; meanwhile a method combining the Wirtinger inequality and reciprocally convex approach is applied to obtain a tight bound of the integral terms of quadratic functions which may lead to a less conservative result compared to previous works. Moreover, the practical significance of this method in engineering will be demonstrated later. Finally, a numerical simulation for a ship is proposed to verify the effectiveness and advantage of the controller we designed. The rest of this paper is organized as follows: In Section 2, the problem formulation of finite-time 𝐻∞ control is detailed for vessels while some definitions and lemmas are introduced as preparation. In Section 3, the method to design finite-time 𝐻∞ controller is proposed, and the proposed control schemes are simulated in Section 4. In Section 5, the conclusion is drawn. Notation. Throughout this paper, R𝑛 is the 𝑛-dimensional Euclidean vector space, and Rπ‘šΓ—π‘› denotes the set of all π‘š Γ— 𝑛 real matrices. For symmetric matrices 𝑋 and π‘Œ, 𝑋 > π‘Œ (respectively, 𝑋 β‰₯ π‘Œ ) means that 𝑋 βˆ’ π‘Œ is positive definite (respectively, positive semidefinite). The superscript β€œπ‘‡β€ represents the transpose. The symmetric terms in a symmetric matrix are denoted by β€œβˆ—β€. Moreover, we use πœ† max (β‹…)(πœ† min (β‹…) ) to denote the maximum (minimum) eigenvalue of a symmetric matrix.

2. Problem Formulation 2.1. Model of DP System. At first, DP system model with three-DOP under low speed can be described as [18]

Mathematical Problems in Engineering πœ‚Μ‡ (𝑑) = 𝐽 (πœ“) 𝜐 (𝑑)

(1)

π‘€πœΜ‡ (𝑑) + 𝐷𝜐 (𝑑) = 𝑒 (𝑑 βˆ’ β„Ž (𝑑)) + πœ” (𝑑)

(2)

where 𝜐 = (πœ‡ V π‘Ÿ) is a vector of velocities given in the body-fixed coordinate system and πœ‚ = (𝑋 π‘Œ πœ“) is the position and orientation of the vessel with respect to an inertial reference coordinate system. 𝑒(𝑑) is a control vector of forces and moments provided by the propulsion systems. β„Ž(𝑑) is a time-varying function that expresses the actuator delay Μ‡ < πœ‡ . πœ”(𝑑) is the disturbance input and satisfies β„Ž(𝑑) < β„Ž β„Ž(𝑑) ∞ of the system and satisfies the condition of ∫0 πœ”(𝑑)𝑇 πœ”(𝑑)𝑑𝑑 < 𝑑 . 𝐽(πœ“) is the transformation matrix between the inertial and body-fixed coordinate frames. The inertia matrix 𝑀 includes hydrodynamic added inertia, and 𝐷 is a strictly positive damping matrix due to linear wave drift damping and laminar friction. Then, the structures of the matrices 𝐽(πœ“), 𝑀, and 𝐷 can be explicitly given as follows: cos (πœ“) βˆ’sin (πœ“) 0 𝐽 (πœ“) = ( sin (πœ“) cos (πœ“) 0) 0

0

π‘š βˆ’ 𝑋𝑒̇

0

0

π‘š βˆ’ π‘ŒVΜ‡

0

π‘šπ‘₯𝑔 βˆ’ π‘Œπ‘ŸΜ‡

𝑀=(

𝑋𝑒 0 𝐷 = βˆ’( 0

1 0 π‘šπ‘₯𝑔 βˆ’ π‘Œπ‘ŸΜ‡)

(3)

𝐼𝑧 βˆ’ π‘π‘ŸΜ‡

0

π‘ŒV π‘Œπ‘Ÿ )

0 𝑁V π‘π‘Ÿ where 𝑋𝑒 ,π‘ŒV ,π‘Œπ‘Ÿ ,π‘π‘Ÿ ,𝑋𝑒̇ ,π‘ŒVΜ‡ ,π‘Œπ‘ŸΜ‡, and π‘π‘ŸΜ‡ are the hydrodynamic parameters of the vessel. π‘š is the mass of the ship, and 𝐼𝑧 is the moment of inertia about the yaw rotation. π‘₯𝑔 is the vertical distance from coordinate origin to center of gravity in body-fixed frame. To simplify the model, some assumptions are introduced first. Assumption 1. All the parameters of state are available. Assumption 2. The roll and pitch angles are small enough; for DP system, it is a reasonable assumption. Under assumption 2, we can obtain the simplified equation in form of state-space [19] π‘₯Μ‡ (𝑑) = 𝐴π‘₯ (𝑑) + 𝐡1 𝑒 (𝑑 βˆ’ β„Ž (𝑑)) + 𝐡2 πœ” (𝑑)

(4)

𝑧 (𝑑) = 𝐢π‘₯ (𝑑) + 𝐷𝑒 (𝑑 βˆ’ β„Ž (𝑑))

(5)

where π‘₯ (𝑑) = (πœ‚ (𝑑) 𝜐 (𝑑)) 𝐴=(

0

𝐼

)

0 βˆ’π‘€βˆ’1 𝐷

Mathematical Problems in Engineering 𝐡1 = 𝐡2 = (

0 βˆ’1

𝑀

3 fact that system state reaches the equilibrium point of system in a finite time [24].

). (6)

Under assumption 1, we design the full-state feedback as 𝑒 (𝑑 βˆ’ β„Ž (𝑑)) = 𝐾π‘₯ (𝑑 βˆ’ β„Ž (𝑑))

(7)

so the (1) and (2) can be rewritten as π‘₯Μ‡ (𝑑) = 𝐴π‘₯ (𝑑) + 𝐡1 𝐾π‘₯ (𝑑 βˆ’ β„Ž (𝑑)) + 𝐡2 πœ” (𝑑)

(8)

𝑧 (𝑑) = 𝐢π‘₯ (𝑑) + 𝐷𝐾π‘₯ (𝑑 βˆ’ β„Ž (𝑑)) .

(9)

Remark 3. In the situation of DP motion, the motion in heave, roll, and pitch will be ignored since we focus on the motion on the surface of sea. Meanwhile, velocity of vessel is low enough, so the Coriolis-Centripetal matrix and nonlinearities in damping matrix can be neglected [18]. Remark 4. In measurement subsystem of DP system, various sensors are installed to obtain the state of the vessel motion accurately, including global position system for the vessel’ s position, gyrocompass for the vessel’ s heading, and attitude sensor for the pitch and roll. Meanwhile, data fusion technology is applied to DP system to obtain more accurate state information of vessel motion [18]. 2.2. Preliminaries. In the sequel, some definitions and lemmas are introduced to obtain the results. Definition 5 ((FTB), see [20]). Given a positive definite matrix 𝑅 and three positive constants 𝑐1 , 𝑐2 , 𝑇 with 𝑐1 < 𝑐2 , the time-delay system (4) with 𝑒(𝑑) = 0 is said to be finite-time boundedness with respect to (𝑐1 𝑐2 𝑑 𝑇 𝜏 𝑅), Μ‡ if supβˆ’πœ 0 and a differentiable signal π‘₯ in [βˆ’π‘Ž, βˆ’π‘] β†’ R𝑛 , the following inequality holds: βˆ’π‘

βˆ’π‘Ž

βˆ’π‘ 1 β‹… Ξ 1 (𝑆1 , π‘₯ (βˆ’π‘) , π‘₯ (βˆ’π‘Ž) , ∫ π‘₯𝑇 (𝑒) 𝑑𝑒) π‘Ž βˆ’ 𝑏 βˆ’π‘Ž 𝑇

Ξ 1 (𝑆1 , 𝛼, 𝛽, 𝛾)

=

sup

β€–πœ”(𝑑)β€–2,[0,𝑇] 2 +π‘₯(0)𝑇 𝑆π‘₯(0)=0ΜΈ

[

‖𝑧 (𝑑)β€–2,[0,𝑇]

β€–πœ” (𝑑)β€–2,[0,𝑇] 2 + π‘₯ (0)𝑇 𝑆π‘₯ (0)

]

Lemma 10 (see [29]). For constant matrices 𝑆1 > 0, 𝑆2 > 0 and constant scalars π‘Ž > 𝑏, the following inequalities hold for all continuously function π‘₯ in [βˆ’π‘Ž, βˆ’π‘] β†’ R𝑛 : βˆ’π‘

π‘‘βˆ’π‘

βˆ’π‘Ž

𝑑+πœƒ

∫ ∫

π‘₯̇𝑇 (𝑒) 𝑆2 π‘₯Μ‡ (𝑒) π‘‘π‘’π‘‘πœƒ

β‰₯ Ξ 2 (𝑆2 , π‘₯𝑇 (𝑑 βˆ’ 𝑏) , βˆ’π‘

π‘‘βˆ’π‘

βˆ’π‘Ž

𝑑+πœƒ

β‹…βˆ« ∫ βˆ’π‘

𝑑+πœƒ

βˆ’π‘Ž

π‘‘βˆ’π‘Ž

∫ ∫

where 𝛾 is a prescribed scalar and 𝑆 is a weighting diagonal matrix that penalizes the effect caused by the initial state. Remark 7. It is necessary to distinguish between finite-time stability and finite-time attractiveness [21]. The first concept is to maintain system states within a given boundary in a specified time interval [22, 23], while the latter describes the

π‘‘βˆ’π‘ 1 1 ∫ π‘₯𝑇 (𝑒) 𝑑𝑒, π‘Ž βˆ’ 𝑏 π‘‘βˆ’π‘Ž (π‘Ž βˆ’ 𝑏)2

π‘₯𝑇 (𝑒) π‘‘π‘’π‘‘πœƒ)

(13)

𝑇

π‘₯Μ‡ (𝑒) 𝑆3 π‘₯Μ‡ (𝑒) π‘‘π‘’π‘‘πœƒ

β‰₯ Ξ 3 (𝑆3 , π‘₯𝑇 (𝑑 βˆ’ π‘Ž) ,

(10)

0,βˆ‘π‘– 𝛼𝑖 =1}

βˆ‘ 𝑖

1 𝑓 (𝑑) = βˆ‘π‘“π‘– (𝑑) + maxβˆ‘π‘”π‘–,𝑗 (𝑑) 𝑔𝑖,𝑗 (𝑑) 𝛼𝑖 𝑖 𝑖 𝑖=𝑗̸

subject to {𝑔𝑖,𝑗 : π‘…π‘š 󳨃󳨀→ 𝑅, 𝑔𝑗,𝑖 (𝑑) = 𝑔𝑖,𝑗 (𝑑) , [

Lemma 12 (see [31]). For any matrix 𝑆4 > 0 and a vector function π‘₯ in [βˆ’π‘Ž, βˆ’π‘] β†’ R𝑛 , if the integrals concerned are well defined, then the following inequality holds: βˆ’π‘

∫ π‘₯𝑇 (𝑠) 𝑆4 π‘₯ (𝑠) 𝑑𝑠 βˆ’π‘Ž

βˆ’π‘ βˆ’π‘ 1 β‰₯ (∫ π‘₯𝑇 (𝑠) 𝑑𝑠) 𝑆4 (∫ π‘₯ (𝑠) 𝑑𝑠) . π‘Ž βˆ’ 𝑏 βˆ’π‘Ž βˆ’π‘Ž

𝑓𝑖 (𝑑) 𝑔𝑖,𝑗 (𝑑) 𝑔𝑗,𝑖 (𝑑) 𝑓𝑗 (𝑑)

(15) ] β‰₯ 0} .

3π‘Š 𝑍2

(

βˆ—

)>0 3π‘Š

1 𝑐1 πœƒ2 + β„Ž2 𝑐1 πœƒ3 + πœŒβ„Žπ‘1 πœƒ4 + πœŒβ„Žπ‘1 πœƒ5 + πœŒβ„Ž3 𝑐1 πœƒ6 2 1 1 + πœŒβ„Ž3 𝑐1 πœƒ7 + πœŒβ„Ž3 𝑐1 πœƒ8 + π›½πœƒ9 𝑑 < π‘’βˆ’π›Όπ‘‡ 𝑐2 πœƒ1 6 3

(16)

where 2.3. Control Objective. Finally, the objective of this paper is to derive the control gain 𝐾 such that (1) the closed-loop system (8) with 𝑒(𝑑) = 0 is FTB; (2) under given nonzero initial condition, the closed 𝑇 system (8) and (9) guarantees that 𝐽 = ∫0 𝑧𝑇 (𝑠)𝑧(𝑠) βˆ’ 𝛾2 πœ”π‘‡ (𝑠)πœ”(𝑠) βˆ’ 𝛾2 (1/𝑇)π‘₯𝑇 (0)𝑆π‘₯(0)𝑑𝑠 < 0 for all nonzero πœ”(𝑑) and 𝑑 ∈ [0, 𝑇].

The designing of the robust finite-time 𝐻∞ controller for the DP system with input delay is divided into three steps. 3.1. FTB Analysis of DP System. Firstly, the result guaranteeing the FTB of DP system is established in this subsection. Theorem 13. Given five positive scalars 𝛼, 𝛽, 𝑐1 , 𝑐2 , 𝑇 and positive define matrix 𝑅, the finite-time boundedness problem of system (4) is solvable if there exist positive scalars πœƒπ‘™ (𝑙 = 1, 2, β‹… β‹… β‹… , 8) and matrices 𝑃11 > 0, 𝑃22 > 0, π‘Š > 0, 𝑄𝑖 > 0 (𝑖 = 1, 2), π‘ˆπ‘— > 0 (𝑗 = 1, 2), 𝑀 > 0, 𝑍1 , 𝑍2 with appropriate dimensions, satisfying the following conditions:

π‘Š 𝑍1

(

βˆ— π‘Š

)>0

1 + 𝐴𝑇 β„Ž2 (π‘ˆ1 + π‘ˆ2 + 2π‘Š) 𝐴 βˆ’ 4πœŒπ‘Š βˆ’ 3πœŒπ‘ˆ1 2 βˆ’ 9πœŒπ‘ˆ2 , Ξ©22 = βˆ’πœŒπ‘„2 βˆ’ 4πœŒπ‘Š βˆ’ 6πœŒπ‘ˆ2 , 1 Ξ©33 = 𝐾𝑇𝐡1 𝑇 β„Ž2 (2π‘Š + π‘ˆ1 + π‘ˆ2 ) 𝐡1 𝐾 βˆ’ 8πœŒπ‘Š 2

3. Main Results

Ξ© = (Ω𝑖𝑗 )8Γ—8 < 0

Ξ©11 = 𝐴𝑇 𝑃11 + 𝑃11 𝐴 βˆ’ 𝛼𝑃11 + 𝑄1 + 𝑄2

(17) (18)

+ πœŒπ‘1 𝑇 + πœŒπ‘1 βˆ’ πœŒπ‘2 𝑇 βˆ’ πœŒπ‘2 βˆ’ (1 βˆ’ πœ‡) 𝑄2 , 1 Ξ©44 = 𝐡2 𝑇 β„Ž2 (2π‘Š + π‘ˆ1 + π‘ˆ2 ) 𝐡2 βˆ’ 𝛽𝑀, 2 Ξ©66 = βˆ’12πœŒπ‘Š, Ξ©77 = βˆ’12πœŒπ‘Š, Ξ©55 = π›Όβ„Ž2 𝑃22 βˆ’ 18πœŒπ‘ˆ2 βˆ’ 6πœŒπ‘ˆ1 , Ξ©88 = βˆ’36πœŒπ‘ˆ1 βˆ’ 36πœŒπ‘ˆ2 , Ξ©12 = πœŒπ‘1 βˆ’ πœŒπ‘2 + 6πœŒπ‘ˆ2 , Ξ©15 = β„Žπ‘ƒ22 + 12πœŒπ‘ˆ2 , Ξ©16 = 6πœŒπ‘Š, Ξ©17 = 2πœŒπ‘2 ,

(19)

(20)

Mathematical Problems in Engineering

5

1 Ξ©13 = 𝑃11 𝐡1 𝐾 + 𝐴𝑇 β„Ž2 (2π‘Š + π‘ˆ1 + π‘ˆ2 ) 𝐡1 𝐾 βˆ’ 2πœŒπ‘Š 2 βˆ’ πœŒπ‘1 βˆ’ πœŒπ‘2 , 1 Ξ©14 = 𝑃11 𝐡2 + 𝐴𝑇 β„Ž2 (2π‘Š + π‘ˆ1 + π‘ˆ2 ) 𝐡2 , 2 Ξ©18 = βˆ’18πœŒπ‘ˆ2 + 6πœŒπ‘ˆ1 ,

̃𝑖 = π‘…βˆ’1/2 𝑄𝑖 π‘…βˆ’1/2 𝑄

(𝑖 = 1, 2, 3) ,

Μƒ = π‘…βˆ’1/2 π‘Šπ‘…βˆ’1/2 , π‘Š (𝑗 = 1, 2) ,

0 < πœƒ1 𝐼 < 𝑃11 < πœƒ2 𝐼,

Ξ©25 = βˆ’β„Žπ‘ƒ22 βˆ’ 6πœŒπ‘ˆ2 ,

0 < 𝑃22 < πœƒ3 𝐼,

Ξ©26 = 2πœŒπ‘2 𝑇 ,

0 < 𝑄1 < πœƒ4 𝐼,

Ξ©27 = 6πœŒπ‘Š,

0 < 𝑄2 < πœƒ5 𝐼,

Ξ©28 = 12πœŒπ‘ˆ2 , Ξ©34 = 𝐾 𝐡1

𝑃̃22 = π‘…βˆ’1/2 𝑃22 π‘…βˆ’1/2 ,

̃𝑗 = π‘…βˆ’1/2 π‘ˆπ‘— π‘…βˆ’1/2 π‘ˆ

Ξ©23 = βˆ’πœŒπ‘1 𝑇 βˆ’ πœŒπ‘2 𝑇 βˆ’ 2πœŒπ‘Š,

𝑇

𝑃̃11 = π‘…βˆ’1/2 𝑃11 π‘…βˆ’1/2 ,

0 < π‘Š < πœƒ6 𝐼,

𝑇1 2

2

β„Ž (2π‘Š + π‘ˆ1 + π‘ˆ2 ) 𝐡2 ,

0 < π‘ˆ1 < πœƒ7 𝐼,

Ξ©36 = 6πœŒπ‘Š + 2πœŒπ‘2 𝑇 ,

0 < π‘ˆ2 < πœƒ8 𝐼,

Ξ©37 = 6πœŒπ‘Š + 2πœŒπ‘2 ,

0 < 𝑀 < πœƒ9 𝐼. (21)

Ξ©58 = 12πœŒπ‘ˆ1 + 24πœŒπ‘ˆ2 , Ξ©67 = βˆ’4πœŒπ‘2 ,

Proof. Consider the candidate augment LKLF as follows:

Ξ©24 = Ξ©35 = Ξ©38 = Ξ©45 = Ξ©46 = Ξ©47 = Ξ©48 = Ξ©56 = Ξ©57 = Ξ©68 = Ξ©78 = 0, πœƒ1 = πœ† min (𝑃̃11 ) πœƒ2 = πœ† max (𝑃̃11 ) , πœƒ3 = πœ† max (𝑃̃22 ) , Μƒ1 ) , πœƒ4 = πœ† max (𝑄 Μƒ2 ) , πœƒ5 = πœ† max (𝑄 Μƒ , πœƒ6 = πœ† max (π‘Š) Μƒ1 ) , πœƒ7 = πœ† max (π‘ˆ Μƒ2 ) , πœƒ8 = πœ† max (π‘ˆ πœƒ9 = πœ† max (𝑀) , 𝜌 = π‘’π›Όβ„Ž ,

𝑉 (𝑑) = 𝑉1 (𝑑) + 𝑉2 (𝑑) + 𝑉3 (𝑑) + 𝑉4 (𝑑)

(22)

𝑉1 (𝑑) = πœ€ (𝑑)𝑇 π‘ƒπœ€ (𝑑)

(23)

𝑉2 (𝑑) = ∫

𝑑

π‘‘βˆ’β„Ž

+∫

𝑒𝛼(π‘‘βˆ’π‘ ) π‘₯ (𝑠)𝑇 𝑄1 π‘₯ (𝑠) 𝑑𝑠 (24)

𝑑

𝛼(π‘‘βˆ’π‘ )

𝑒

π‘‘βˆ’β„Ž(𝑑)

0

𝑉3 (𝑑) = β„Ž ∫ ∫

𝑑

βˆ’β„Ž 𝑑+πœƒ

0

0

0

π‘₯ (𝑠) 𝑄2 π‘₯ (𝑠) 𝑑𝑠

𝑒𝛼(π‘‘βˆ’π‘ ) π‘₯Μ‡ (𝑠)𝑇 π‘Šπ‘₯Μ‡ (𝑠) π‘‘π‘ π‘‘πœƒ

𝑑

𝑉4 (𝑑) = ∫ ∫ ∫ βˆ’β„Ž πœƒ

𝑇

𝑑+𝛽

πœƒ

(25)

𝑒𝛼(π‘‘βˆ’π‘ ) π‘₯Μ‡ (𝑠)𝑇 π‘ˆ1 π‘₯Μ‡ (𝑠) π‘‘π‘ π‘‘π›½π‘‘πœƒ

+∫ ∫ ∫

𝑑

βˆ’β„Ž βˆ’β„Ž 𝑑+𝛽

(26) 𝛼(π‘‘βˆ’π‘ )

𝑒

𝑇

π‘₯Μ‡ (𝑠) π‘ˆ2 π‘₯Μ‡ (𝑠) π‘‘π‘ π‘‘π›½π‘‘πœƒ

𝑑

where πœ€(𝑑) = (π‘₯(𝑑) βˆ«π‘‘βˆ’β„Ž π‘₯(𝑠)𝑑𝑠) 𝑃 = ( π‘ƒβˆ—11 𝑃022 ). To deal with the formulas conveniently, let us provide such definitions:

πœ‰ (𝑑) = (π‘₯ (𝑑) π‘₯ (𝑑 βˆ’ β„Ž) π‘₯ (𝑑 βˆ’ β„Ž (𝑑)) πœ” (𝑑) 0 β‹… β‹… β‹… 0) . 0 β‹… β‹… β‹… 0 1 ⏟⏟⏟⏟⏟⏟⏟⏟⏟ π‘™π‘˜ 𝑇 = (⏟⏟⏟⏟⏟⏟⏟⏟⏟ π‘˜βˆ’1

8βˆ’π‘˜

𝑑 π‘‘βˆ’β„Ž(𝑑) 1 1 1 0 𝑑 1 𝑑 π‘₯ (𝑠) 𝑑𝑠 π‘₯ (𝑠) 𝑑𝑠 2 ∫ ∫ π‘₯ (𝑠) π‘‘π‘ π‘‘πœƒ) (27) ∫ π‘₯ (𝑠) 𝑑𝑠 ∫ ∫ β„Ž π‘‘βˆ’β„Ž β„Ž (𝑑) π‘‘βˆ’β„Ž(𝑑) β„Ž βˆ’ β„Ž (𝑑) π‘‘βˆ’β„Ž β„Ž βˆ’β„Ž 𝑑+πœƒ

6

Mathematical Problems in Engineering

The time-derivative of 𝑉(𝑑) along system (4) can be bounded as

Ξ›1 = Ξ›2 = (

𝑉1Μ‡ (𝑑) = 𝛼𝑉1 (𝑑) + πœ€ Μ‡ (𝑑)𝑇 π‘ƒπœ€ (𝑑) + πœ€ (𝑑)𝑇 π‘ƒπœ€ Μ‡ (𝑑) βˆ’ 𝛼𝑉1 (𝑑) (28) 𝑉2Μ‡ (𝑑)

𝑇

𝑉3Μ‡ (𝑑) ≀ 𝛼𝑉3 (𝑑) + β„Ž2 π‘₯Μ‡ (𝑑)𝑇 π‘Šπ‘₯Μ‡ (𝑑) βˆ’ π‘’π›Όβ„Ž πœ›3 𝑇 Ξ› 3 πœ›3

𝑇

βˆ’ (1 βˆ’ β„ŽΜ‡ (𝑑)) π‘’π›Όβ„Ž(𝑑) π‘₯ (𝑑 βˆ’ β„Ž (𝑑))𝑇 𝑄2 π‘₯ (𝑑 βˆ’ β„Ž (𝑑))

where πœ›3 𝑇 = ( πœ›πœ›1 𝑇 ), Ξ› 3 = ( βˆ—βˆ— 2

1 𝑉4Μ‡ (𝑑) ≀ 𝛼𝑉4 (𝑑) + β„Ž2 π‘₯Μ‡ (𝑑)𝑇 π‘ˆ1 π‘₯Μ‡ (𝑑) 2

βˆ’ (1 βˆ’ πœ‡) π‘₯ (𝑑 βˆ’ β„Ž (𝑑))𝑇 𝑄2 π‘₯ (𝑑 βˆ’ β„Ž (𝑑))

1 + β„Ž2 π‘₯Μ‡ (𝑑)𝑇 π‘ˆ2 π‘₯Μ‡ (𝑑) βˆ’ π‘’π›Όβ„Ž πœ›4 𝑇 Ξ› 4 πœ›4 2

𝑉3Μ‡ (𝑑) = 𝛼𝑉3 (𝑑) + β„Ž2 π‘₯Μ‡ (𝑑)𝑇 π‘Š2 π‘₯Μ‡ (𝑑)

(30)

πœ›4 𝑇

1 = 𝛼𝑉4 (𝑑) + β„Ž2 π‘₯Μ‡ (𝑑)𝑇 π‘ˆ1 π‘₯Μ‡ (𝑑) 2 βˆ’π‘’

∫ ∫

𝑑

0

𝑑+𝛽

βˆ’β„Ž π‘‘βˆ’β„Ž

π‘₯Μ‡ (𝑠) π‘ˆ1 π‘₯Μ‡ (𝑠) 𝑑𝑠𝑑𝛽

6π‘ˆ1 βˆ’12π‘ˆ1

(31) Ξ›4 = (

π‘₯Μ‡ (𝑠)𝑇 π‘ˆ2 π‘₯Μ‡ (𝑠) 𝑑𝑠𝑑𝛽.

𝑉3Μ‡ (𝑑) = 𝛼𝑉3 (𝑑) + β„Ž2 π‘₯Μ‡ (𝑑)𝑇 π‘Šπ‘₯Μ‡ (𝑑) βˆ’ π‘’π›Όβ„Ž β„Ž ∫

π‘‘βˆ’β„Ž(𝑑) π‘‘βˆ’β„Ž(𝑑)

βˆ’ π‘’π›Όβ„Ž β„Ž ∫

π‘‘βˆ’β„Ž

π‘₯Μ‡ (𝑠)𝑇 π‘Šπ‘₯Μ‡ (𝑠) 𝑑𝑠

βˆ—

βˆ—

βˆ—

0

6π‘ˆ2 βˆ’12π‘ˆ2 βˆ—

).

36π‘ˆ2

𝑉 (𝑑) < 𝑒𝛼𝑑 𝑉 (0) + 𝛽 ∫ 𝑒𝛼(π‘‘βˆ’π‘ ) πœ”π‘‡ (𝑠) π‘€πœ” (𝑠) 𝑑𝑠 0

where 𝑇 𝑇

(𝑙1 βˆ’ 𝑙3 ) πœ›1 𝑇 = πœ‰π‘‡ (𝑑) ( ) 𝑇 𝑇 (𝑙1 + 𝑙3 𝑇 βˆ’ 2𝑙6 𝑇 ) 𝑇

(𝑙3 𝑇 βˆ’ 𝑙2 𝑇 )

(37)

𝑑

(32)

β„Ž β„Ž βˆ’ π‘’π›Όβ„Ž πœ›1 𝑇 Ξ› 1 πœ›1 βˆ’ π‘’π›Όβ„Ž πœ› 𝑇Λ πœ› β„Ž (𝑑) β„Ž βˆ’ β„Ž (𝑑) 2 2 2

𝑇

βˆ—

0

Assuming that Ξ© < 0 and integrating the left part of inequality (37),

≀ 𝛼𝑉3 (𝑑) + β„Ž2 π‘₯Μ‡ (𝑑)𝑇 π‘Šπ‘₯Μ‡ (𝑑)

𝑇

36π‘ˆ1

0

𝑉̇ βˆ’ 𝛼𝑉 βˆ’ π›½πœ”π‘‡ (𝑑) π‘€πœ” (𝑑) < πœ‰π‘‡ (𝑑) Ξ©πœ‰ (𝑑) .

π‘₯Μ‡ (𝑠)𝑇 π‘Šπ‘₯Μ‡ (𝑠) 𝑑𝑠

𝑇

βˆ—

0

(36)

Combining (28), (29), (34), and (35) with the definition of 𝜌 = π‘’π›Όβ„Ž , we can obtain

Invoking Lemma 9, we can obtain

𝑑

) ) ) )

𝑇 1 ( 𝑙2 𝑇 βˆ’ 𝑙5 𝑇 + 𝑙8 𝑇 ) ) ( 2

1 + β„Ž2 π‘₯Μ‡ (𝑑)𝑇 π‘ˆ2 π‘₯Μ‡ (𝑑) 2 βˆ’ π‘’π›Όβ„Ž ∫ ∫

𝑇 1 ( 𝑙1 𝑇 βˆ’ 𝑙8 𝑇 ) 2 𝑇 (𝑙2 𝑇 βˆ’ 𝑙5 𝑇 )

( ( = πœ‰π‘‡ (𝑑) ( (

𝑇

βˆ’β„Ž 𝑑+𝛽

𝑇

𝑇

(𝑙1 𝑇 βˆ’ 𝑙2 𝑇 )

π‘₯Μ‡ (𝑠)𝑇 π‘Š2 π‘₯Μ‡ (𝑠) 𝑑𝑠

0

(35)

where

𝑉4Μ‡ (𝑑)

π›Όβ„Ž

0 0 𝑍2 ). 3π‘Š

ity 𝑇

βˆ’ 𝑒 π‘₯ (𝑑 βˆ’ β„Ž) 𝑄1 π‘₯ (𝑑 βˆ’ β„Ž) + π‘₯ (𝑑) 𝑄2 π‘₯ (𝑑)

π‘‘βˆ’β„Ž

0 0 3π‘Š βˆ—

(34)

By applying Lemma 10 to (31), we can obtain the inequal-

𝑇

𝑑

π‘Š 𝑍1 π‘Š βˆ— βˆ— βˆ—

𝑇

(29)

≀ 𝛼𝑉2 (𝑑) + π‘₯ (𝑑)𝑇 𝑄1 π‘₯ (𝑑)

βˆ’ π‘’π›Όβ„Ž β„Ž ∫

). βˆ— 3π‘Š (33)

βˆ’ 𝑒 π‘₯ (𝑑 βˆ’ β„Ž) 𝑄1 π‘₯ (𝑑 βˆ’ β„Ž) + π‘₯ (𝑑) 𝑄2 π‘₯ (𝑑)

π›Όβ„Ž

0

In order to obtain a tighter bound of integral term, Lemma 11 is applied to (32) as follows:

= 𝛼𝑉2 (𝑑) + π‘₯ (𝑑)𝑇 𝑄1 π‘₯ (𝑑) π›Όβ„Ž

π‘Š

𝑇

πœ›2 = πœ‰ (𝑑) ( ) (𝑙3 𝑇 + 𝑙2 𝑇 βˆ’ 2𝑙7 𝑇 )

𝑇

π‘₯ (𝑑) 𝑃11 π‘₯ (𝑑) = π‘₯

𝑇

(𝑑) 𝑅1/2 𝑃̃11 𝑅1/2 π‘₯ (𝑑)

> πœ† min (𝑃̃11 ) π‘₯𝑇 (𝑑) 𝑅π‘₯ (𝑑) .

Μƒ 11 = 𝑃11 𝐴𝑇 + 𝐴𝑃11 βˆ’ 𝛼𝑃11 + 𝑄1 + 𝑄2 βˆ’ 4πœŒπ‘Š Ξ  βˆ’ 3πœŒπ‘ˆ1 βˆ’ 9πœŒπ‘ˆ2 ,

(40)

Μƒ 44 = βˆ’π›Ύ2 𝐼, Ξ  Μƒ 22 = βˆ’πœŒπ‘„2 βˆ’ 4πœŒπ‘Š βˆ’ 6πœŒπ‘ˆ2 , Ξ 

Therefore, conditions (17) to (20) can guarantee the FTB of system (4). This completes the proof. Remark 14. Among the existing approaches, there are two threads: one is to construct a novel LKLF that involves more information of delay; the other is to find a tighter estimation of upper bound for cross terms coming from the derivative of the LKLF. In this paper, these two techniques are applied to obtain the result. In addition, the less conservatism in practical engineering will be shown in the simulations which are always ignored in most literature sources.

Μƒ 33 = βˆ’8πœŒπ‘Š + πœŒπ‘1 𝑇 + πœŒπ‘1 βˆ’ πœŒπ‘2 𝑇 βˆ’ πœŒπ‘2 Ξ  βˆ’ (1 βˆ’ πœ‡) 𝑄2 , Μƒ 13 = 𝐡1 π‘Œ βˆ’ 2πœŒπ‘Š βˆ’ πœŒπ‘1 βˆ’ πœŒπ‘2 , Ξ  Μƒ 14 = 𝐡2 , Ξ  Μƒ 66 = βˆ’12πœŒπ‘Š, Ξ  Μƒ 77 = βˆ’12πœŒπ‘Š, Ξ 

3.2. Controller Design. In this subsection, we focus on the problem of finite-time 𝐻∞ state feedback designing based on Theorem 13, that is, designing a state feedback controller in the form of (7) such that the resulting DP system satisfies the control objective proposed in Section 1.

Μƒ 55 = π›Όβ„Ž2 𝑃22 βˆ’ 18πœŒπ‘ˆ2 βˆ’ 6πœŒπ‘ˆ1 , Ξ 

Theorem 15. For given positive 𝛼, 𝑐1 , 𝑐2 , 𝑇, 𝛾 and matrices 𝑅 > 0, 𝑆 β‰₯ 0, if there exist positive scalars πœƒπ‘™ (𝑙 = 1, 2, β‹… β‹… β‹… , 8) and matrices 𝑃11 > 0, 𝑃22 > 0, π‘Š > 0, 𝑄𝑖 > 0 (𝑖 = 1, 2), π‘ˆπ‘— > 0 (𝑗 = 1, 2), π‘Œ, 𝑍1 , 𝑍2 with appropriate dimensions, satisfying the following conditions:

Μƒ 15 = β„Žπ‘ƒ22 + 12πœŒπ‘ˆ2 , Ξ 

Μƒ = (Ξ  Μƒ 𝑖𝑗 ) 0 βˆ— π‘Š 3π‘Š 𝑍2

(

βˆ—

3π‘Š

)>0

(45)

(43)

Μƒ 88 = βˆ’36πœŒπ‘ˆ1 βˆ’ 36πœŒπ‘ˆ2 , Ξ  Μƒ 12 = πœŒπ‘1 βˆ’ πœŒπ‘2 + 6πœŒπ‘ˆ2 , Ξ 

Μƒ 16 = 6πœŒπ‘Š, Ξ  Μƒ 17 = 2πœŒπ‘2 , Ξ  Μƒ 25 = βˆ’β„Žπ‘ƒ22 βˆ’ 6πœŒπ‘ˆ2 , Ξ  Μƒ 18 = βˆ’18πœŒπ‘ˆ2 + 6πœŒπ‘ˆ1 , Ξ  Μƒ 26 = 2πœŒπ‘2 𝑇 , Ξ  Μƒ 23 = βˆ’πœŒπ‘1 𝑇 βˆ’ πœŒπ‘2 𝑇 βˆ’ 2πœŒπ‘Š, Ξ  Μƒ 27 = 6πœŒπ‘Š, Ξ 

(44)

Μƒ 28 = 12πœŒπ‘ˆ2 , Ξ 

(46)

8

Mathematical Problems in Engineering Μƒ 36 = 6πœŒπ‘Š + 2πœŒπ‘2 𝑇 , Ξ 

Ξ£99 = βˆ’π›Ύ2 𝑑,

Μƒ 37 = 6πœŒπ‘Š + 2πœŒπ‘2 , Ξ 

Ξ£12 = βˆšπ‘1 𝐻2 ,

Μƒ 58 = 12πœŒπ‘ˆ1 + 24πœŒπ‘ˆ2 , Ξ 

Ξ£13 = βˆšβ„Ž2 𝑐1 𝐻3 ,

Μƒ 67 = βˆ’4πœŒπ‘2 , Ξ 

Ξ£14 = βˆšπœŒβ„Žπ‘1 𝐻4 ,

Μƒ 19 = 1 β„Žπ‘ƒ11 𝐴𝑇 , Ξ  2

Ξ£15 = βˆšπœŒβ„Žπ‘1 𝐻5 ,

Μƒ 39 = 1 β„Žπ‘Œπ‘‡ 𝐡1 𝑇 , Ξ  2

1 Ξ£16 = √ πœŒβ„Ž3 𝑐1 𝐻6 , 2

Μƒ 49 = 1 β„Žπ΅2 𝑇 , Ξ  2

1 Ξ£17 = √ πœŒβ„Ž3 𝑐1 𝐻7 , 6

Μƒ 1,10 = 1 β„Žπ‘ƒ11 𝐴𝑇 , Ξ  2

1 Ξ£18 = √ πœŒβ„Ž3 𝑐1 𝐻8 , 3

Μƒ 3,10 = 1 β„Žπ‘Œπ‘‡ 𝐡1 𝑇 , Ξ  2

Ξ£19 = 𝑃11 ,

Μƒ 4,10 = 1 β„Žπ΅2 𝑇 , Ξ  2

(48) others

Μƒ 1,11 = β„Žπ‘ƒ11 𝐴𝑇 , Ξ 

Σ𝑖𝑗 = 0

Μƒ 3,11 = β„Žπ‘Œπ‘‡ 𝐡1 𝑇 , Ξ 

then a state feedback 𝐻∞ controller in form of (7) exists, such that (1) the closed-loop system (8) with 𝑒(𝑑) = 0 is FTB; (2) under given nonzero initial condition, the closed 𝑇 system (8) and (9) guarantees that 𝐽 = ∫0 𝑧𝑇 (𝑠)𝑧(𝑠) βˆ’ 𝛾2 πœ”π‘‡ (𝑠)πœ”(𝑠) βˆ’ 𝛾2 (1/𝑇)π‘₯𝑇 (0)𝑆π‘₯(0)𝑑𝑠 < 0 for all nonzero πœ”(𝑑) and 𝑑 ∈ [0, 𝑇]. βˆ’1 The controller gain can be calculated by 𝐾 = π‘Œπ‘ƒ11 .

Μƒ 4,11 = β„Žπ΅2 𝑇 , Ξ  Μƒ 1,12 = 𝑃11 𝐢𝑇 , Ξ  Μƒ 3,12 = π‘Œπ‘‡ 𝐷𝑇 , Ξ  Μƒ 99 = 1 π‘ˆ1 βˆ’ 𝑃11 , Ξ  2 Μƒ 10,10 = 1 π‘ˆ2 βˆ’ 𝑃11 , Ξ  2

Proof. In preparation for designing, we set 𝑀 = 𝛾2 𝐼, 𝛽 = 1 for (37); then (37) and (38) can be rewritten as follows:

Μƒ 11,11 = π‘Š βˆ’ 2𝑃11 , Ξ  Μƒ 12,12 = βˆ’π‘’βˆ’π›Όπ‘‡ 𝐼, Ξ 

Μ‚ (𝑑) 𝑉̇ βˆ’ 𝛼𝑉 βˆ’ 𝛾2 πœ”π‘‡ (𝑑) πœ” (𝑑) < πœ‰π‘‡ (𝑑) Ξ©πœ‰ (47)

𝑑

0

Μƒ 𝑖𝑗 = 0, Ξ 

Ξ£22 = βˆ’π»2 , Ξ£33 = βˆ’π»3 , Ξ£44 = βˆ’π»4 , Ξ£55 = βˆ’π»5 , Ξ£66 = βˆ’π»6 , Ξ£77 = βˆ’π»7 , Ξ£88 = βˆ’π»8 ,

(50)

𝑉 (𝑑) < 𝑒𝛼𝑑 𝑉 (0) + 𝛾2 ∫ 𝑒𝛼(π‘‘βˆ’π‘ ) πœ”π‘‡ (𝑠) πœ” (𝑠) 𝑑𝑠

others

Ξ£11 = βˆ’π‘’βˆ’π›Όπ‘‡ 𝑐2 𝐻1 ,

(49)

𝛼𝑇