Leader – Follower Formation Maneuvers for Multi-Robot Systems via Derivative and Integral Terminal Sliding Mode

This paper investigates the formation problem of multiple robots based on the leader–follower mechanism. At first, the dynamics of such a leader–follower framework are modeled. The input–output equations are depicted by calculating the relative degree of a leader–follower formation system. Furthermore, the derivative and integral terminal sliding mode controller is designed based on the relative degree. Since the formation system suffers from uncertainties, the nonlinear disturbance observer is adopted to deal with the uncertainties. The stability of the closed-loop control system is proven in the sense of Lyapunov. Finally, some numerical simulations are displayed to verify the feasibility and effectiveness by the designed controller and observer.


Introduction
In recent years, the coordination control scheme of multiple robots has drawn considerable attention in various fields [1].Multiple robots can be applied in many dangerous places to free the human being, including the earthquake rescue, the warehouse translations, and some tasks at nuclear power plants.A multi-robot system can be treated as a coupling network of some robots, where the robots communicate with each other to achieve some complex duties [2,3].Various investigations have been explored to achieve the coordination control of multiple robots.These investigations can be roughly classified into leader-follower formations [4][5][6][7][8], virtual structure mechanisms [9][10][11][12], graph-based approaches [13,14], and behavior-based methods [15,16].
The leader-follower formations are attractive in the coordination control of multi-robot systems.Partly, such formations benefit multiple robots because the formations can have guaranteed formation stability via control design [17].The basic control idea of the leader-follower mechanism is that multiple robots are divided into several leader-follower pairs.In the leader-follower mechanism, all follower robots share the same leader.In each pair, the leader robot moves along the predefined trajectory, while the follower robots track the leader with desired relative distance and angle.In the leader-follower system of multiple robots, only partial followers can obtain the state of the leader, and the interaction between follower robots and leader robot is local [18].Many control methods have been applied in the leader-follower multi-robot systems, such as sliding mode control (SMC) based on nonlinear disturbance observer [18], SMC [19], second-order SMC [20], adaptive control [21,22], predictive control [23], integral terminal SMC [24], and terminal SMC [25].
In actuality, it is unavoidable for any robots to be affected by uncertainties such as external disturbances, unmodeled dynamics, and parameter perturbations [26].The dynamics of multi-robot systems becomes uncertain, due to these uncertainties [27].These uncertainties can be categorized into unmatched uncertainties and matched uncertainties [28].However, the SMC method, as a strong robust tool, has invariant nature to the matched uncertainties when an SMC system enters into the sliding mode.Unfortunately, the effects of unmatched uncertainties cannot be suppressed by the SMC methods [29].The unmatched uncertainties can challenge the performance of the SMC system seriously.The characteristic of terminal SMC (T-SMC) has its nonlinear sliding surface.Compared with those traditional SMC approaches, the T-SMC method has faster convergence speed and higher accuracy.However, the T-SMC method has the singular problem due to its fractional function.Therefore, the derivative and integral T-SMC (DIT-SMC) method is proposed [30].The DIT-SMC method is of merit.Due to the existence of the integral term, the sliding mode of the DIT-SMC method starts on the derivative and integral terminal sliding mode surface.Moreover, the DIT-SMC method can guarantee the exact estimation of finite error convergence time, and resolve the singular problem of the T-SMC.On the other hand, the derivative term of the DIT-SMC method can reduce the nonlinear effects to the stability of a DIT-SMC system.
In the previous works [21][22][23][24][25][26][27][28][29][30], the assumption that the uncertainties have a known boundary is assumed.Concerning the formation maneuvers of multi-robot systems, the assumption is not mild.In fact, the boundary of uncertainties in multiple robots is hard to be known exactly in advance.In case of the lack of the important information, several serious problems may be raised in reality, for example, the decrease of the formation robustness, the deterioration of the formation performance, as far as the deficiency of the formation stability.In order to resolve the problem of the uncertainties, the nonlinear disturbance observer is adopted.The unknown unmatched uncertainties are estimated by the nonlinear disturbance observer.The technique of nonlinear disturbance observer (NDOB) can handle the unmatched uncertainties problem and improve the robustness of the formation control system.This paper deals with the formation problem of multiple robots with uncertainties.The control scheme combining derivative and integral terminal sliding mode and nonlinear disturbance observer is investigated.The derivative and integral terminal sliding mode method allows the system start on the sliding surface.The reaching time of sliding surface is eliminated.The matched uncertainties in formation system are suppressed by the DIT-SMC method.Under the mild assumption that the uncertainties have an unknown boundary, the NDOB is designed to estimate the unmatched uncertainties in the formation system.The estimate errors will converge to zero in the limited time by setting the parameter of NDOB.In the sense of Lyapunov, the system stability is guaranteed in spite of uncertainties.Finally, some numerical simulations are displayed to illustrate the feasibility and effectiveness.

Modeling of Single Robot
Shown by Figure 1, a unicycle-like robot is taken into account.The robot is round, with r in radius, and has two parallel wheels controlled independently by two DC motors.Because the robot is capable of simultaneous arbitrary rotation and translation in the horizontal plane, a three dimensional vector q = [x, y, θ] T is used to describe the robot.In Figure 1, (x, y) represents the translational coordinates of the robot, and is the center of the robot.The rotational coordinate is depicted by the variable θ.There are n robots in the formation system of multi-robot.Provided the pure rolling and noslipping condition, the ideal dynamic models of the nth robot are described by where n v , n ω are the linear velocities and angular velocities, respectively.
Differentiate (1) with respect to time t.Considering the parameter fluctuations, model uncertainties, and external disturbances, such as slipping and skidding effects in the formation system, the dynamic model of the nth robot is obtained by ( ) where un = [αn βn] T is the control input of nth robot.αn, βn are the acceleration and angular acceleration respectively, which are described by αn = Fn/mn, βn = τn/Jn.Here, Fn, mn, τn, and Jn donate the force, the nominal mass, the torque of the robot, and the nominal moment of inertia, respectively.∆n represents the parameter fluctuations, written by where n ε , n ε ′ represent the variant on the mass and the inertia.
( ) q q π  is described by [πnx πny πnθ] T , meaning the uncertainties and external disturbances in the lumped model.

Leader-Follower Formation Framework
In this section, the kinematics model of the leader-follower formation system is given.The leader-follower formation mechanism is displayed in Figure 2. In the leader-follower formation system, there is a leader robot, and others are selected as follower robots.The ith robot is set as the leader robot, and the kth robot is picked up as the representative of all follower robots.The relative distance lik and relative bearing angle ϕik between the leader robot and follower robot are defined in Figure 2. The relative distance lik means the distance between the center of the leader robot ith and the front castor of the follower robot kth, described by ( ) ( ) Here, (xi, yi) denotes the center of the leader robot i, and ( , ) x y represents the caster position of the follower robot k.The calculation of k x , k y has the following form There are n robots in the formation system of multi-robot.Provided the pure rolling and no-slipping condition, the ideal dynamic models of the nth robot are described by .
where v n , ω n are the linear velocities and angular velocities, respectively.Differentiate (1) with respect to time t.Considering the parameter fluctuations, model uncertainties, and external disturbances, such as slipping and skidding effects in the formation system, the dynamic model of the nth robot is obtained by where u n = [α n β n ] T is the control input of nth robot.α n , β n are the acceleration and angular acceleration respectively, which are described by α n = F n /m n , β n = τ n /J n .Here, F n , m n , τ n , and J n donate the force, the nominal mass, the torque of the robot, and the nominal moment of inertia, respectively.∆ n represents the parameter fluctuations, written by where ε n , ε n represent the variant on the mass and the inertia.π n q n , .
q n is described by [π nx π ny π nθ ] T , meaning the uncertainties and external disturbances in the lumped model.

Leader-Follower Formation Framework
In this section, the kinematics model of the leader-follower formation system is given.The leader-follower formation mechanism is displayed in Figure 2. In the leader-follower formation system, there is a leader robot, and others are selected as follower robots.The ith robot is set as the leader robot, and the kth robot is picked up as the representative of all follower robots.The relative distance l ik and relative bearing angle ϕ ik between the leader robot and follower robot are defined in Figure 2. The relative distance l ik means the distance between the center of the leader robot ith and the front castor of the follower robot kth, described by (3) Here, (x i , y i ) denotes the center of the leader robot i, and (x k , y k ) represents the caster position of the follower robot k.The calculation of x k , y k has the following form Here, r is the radius of the round robot, and (x k , y k ) denotes the center of the follower robot k.
Simultaneously, ψ ik is formulated by Here, θ i denotes the orientation angle of the leader robot i, Here, r is the radius of the round robot, and (xk, yk) denotes the center of the follower robot k.
Simultaneously, ik ψ is formulated by Here, θi denotes the orientation angle of the leader robot i, sin arctan cos In this paper, the derivative and integral terminal sliding mode controller is designed so that the follower robots can follow the leader robot with desired relative distance and angler.Therefore, the following conditions are satisfied: The collisions between the robots are avoided.There is no communication delay between the leader robot and the follower robot.Each follower robot knows its position, velocity, and corresponding information of the leader robot.
According to leader-follower formation mechanism, the robots move along a specified trajectory with desired relative distance and bearing angle.It is necessary to shape the dynamics of the leaderfollower formation system.Differentiate (3) and ( 5) twice with the respect to time t, and substitute (2) into the second derivative of (3) and (5).Define the state variable The dynamic model of the formation system has the form of ( , Here, ,1 ik B is a 2 × 2 matrix whose columns are smooth vector fields , Aik, Bik,1, Bik,2, h(xik) are described by . In this paper, the derivative and integral terminal sliding mode controller is designed so that the follower robots can follow the leader robot with desired relative distance and angler.Therefore, the following conditions are satisfied: The collisions between the robots are avoided.There is no communication delay between the leader robot and the follower robot.Each follower robot knows its position, velocity, and corresponding information of the leader robot.
According to leader-follower formation mechanism, the robots move along a specified trajectory with desired relative distance and bearing angle.It is necessary to shape the dynamics of the leader-follower formation system.Differentiate (3) and ( 5) twice with the respect to time t, and substitute (2) into the second derivative of (3) and (5).Define the state variable ψ ik .The dynamic model of the formation system has the form of .
Here, B ik,1 is a 2 × 2 matrix whose columns are smooth vector fields B ik,1K .h(x ik ) is the output equation of formation system.Here, , h(x ik ) are described by Here, ϕ ik = ψ ik + θ ik , d ik (t) denotes the uncertainties of the leader-follower formation system (6), written by F 1 , F 2 , P 1 , P 2 are depicted respectively by

Control Problem Formulation
Considering the dynamic mode ( 6) of the leader-follower formation system, a relative-degree is calculated by Here, r K (K = 1, 2) is the smallest integer so that the least one of the control inputs appears in y Here, r 1 = r 2 = 2, m = 2. L B ik,1K , L f are lie derivatives.Further, the input-output dynamic equation is depicted by .. Here, Here, I 3 , I 4 are 2 × 2 matrices, u i is the control input of ith robot, and other matrices are described by Hypothesis 1. H(x ik , d ik ) has the normal part G ik and nonlinear part G ik ∆ k , and meets the following in equation Here, δ ik > 1, H dik = H x d ik , I is the 2 × 2 square matrix,x d ik denotes the desired state vector in (6).
Remark 1. G ik ∆ k u k is the matched uncertainties in (13), meaning the parameter fluctuations of the follower robot k.The term d ik depicts the unmatched uncertainties in the leader-follower formation system, and consists of three parts.Due to the formation, framework ( 13) is applied to the follower robot, and the information of leader robot u i can't be matched.
The terms F ik , P ik denote the model uncertainties and external disturbances caused by slipping, friction, and obstacles etc., which are also hard to be matched.

Control Design
Due to the inherent characteristics of centralization, the scheme mainly depends on the leader robots and exists as the "single point of failure" problem.In order to develop derivative and integral terminal sliding mode approach to coordinate the leader robot i and follower robot k, a recursive structure of the terminal sliding function for high relative-degree MIMO systems (with r 1 , r 2 > 1) is designed as .γ ik,1 , γ ik,2 , λ ik,1 , λ ik,2 are all positive constants.p Kj > q Kj , here K, j = 1, 2. p Kj and q Kj are all odd positive constants.
Theorem 1. Considering the derivative and integral terminal sliding mode surface S ik (t) with the fractional function, the state error of formation system can reach the equilibrium point e = 0 at the limited time Here, t 1K (K = 1, 2) is the reaching time of terminal slide mode e D1K .
Proof.From ( 17), the sliding mode s ik starts on t = 0.Then, the equations e ik,D .
The converge time of sliding mode e ik,I1 and e ik,I2 can obtained by solving (19).
1−q 21 /p 21 λ ik,1 (1 − q 21 /p 21 ) , t 12 = e ik,D12 (0) In sliding mode s ik = 0, e ik,D 11 = −λ ik,1 e ik,I1 , and e ik,D 12 = −λ ik,2 e ik,I2 can always hold true.Therefore, the reaching time of e ik,D 11 and e ik,D 12 are the same as the convergence time of e ik,I1 and e ik,I2 , respectively.When e ik,D 11 = e ik,D 12 = 0, e ik,I1 and e ik,I2 will converge to zero successfully.At t = t 1K (K = 1, 2), the e ik,I1 and e ik,I2 are formulated by   e ik,D01 = 0 will spend the time 1−q 11 /p 11 1−q 12 /p 12 Since the sliding mode s ik = 0 consists of the derivative term and integral term, the time T K spending from s ik.K = 0 (K = 1, 2) to e ik,K = 0 (K = 1, 2) is the summation of the two terms.Since the fact that each sliding mode s ik,1 , s ik,2 is independent, the time spent for equilibrium point is the max of the T K .

Hypothesis 2. d ik
The derivative and integral terminal sliding mode control law is set as where sign(s ik ) = sign s ik,1 s ik,2 T , k eik , η ik are all the positive constant set by designer.d ik (t) is the upper bound of the unmatched uncertainties.
V ik < 0 can be guaranteed when k iek > d ik holds true.However, apart from Hypothesis 2, the unmatched uncertainties d ik in (13) are unknown, which means that the upper bound is also unknown.Therefore, it cannot select an appropriate parameter k iek to guarantee .
V ik < 0. Therefore, the stability of control system cannot be guaranteed.

DIT-SMC Design Based NDOB
In order to resolve the above problem, the nonlinear disturbance is proposed to estimate the uncertainties d ik in the leader-follower formation system (13).At first, the following assumption is taken into account.
Considering the formation dynamic model ( 6), the nonlinear disturbance observer is formulated by . Here are the state vector of nonlinear disturbance observer, the observer gain matrix set by designer, and the estimated value of unmatched uncertainties respectively.Define the estimate error vector as Differentiate e dik with respect to time t and take the Hypothesis 3 into account.Furthermore, the dynamics of e dik is presented as .
The solution of ( 30) is e dik = exp(−L ik B ik,2 t)e dik (0), which indicates the estimate error will exponentially converge to zero as t → ∞ if L ik B ik,2 is set as a positive constant.Here, e dik (0) is the initial state of e dik .
Considering input-output dynamics (13) and observer (29), the control law based on NDOB is determined by Theorem 2. Consider the dynamic model of leader-follower formation system (6), take the assumption 1, 2, 3, 4 into account, adopt the input-output model (13), design the derivative and integral terminal sliding mode surface (17) and nonlinear disturbance observer (29).If the derivative and integral terminal control law is set as (31), the leader-follower formation system with unmatched uncertainties is asymptotically stable when k iek > e * d .
Proof.Selecting the Lyapunov function as V ik = 1 2 s . s, differentiating V ik with the respect to time t and substituting the .s ik into the derivative of V ik yields .
Here, J ik = diag .e (P 11 /q 11 −1) D01 , .e (P 12 /q 12 −1) D02 . Since p 1K > 0, q 1K > 0, p 2K > 0, q 2K > 0, p 1K > q 1K , p 2k > q 2K exist in the controller,  = 0.According to the Hypothesis 1, ( 27), the second term and the third term can be deduced by In ( 33), the condition k eik > 1/(1 − χ ik ) can be picked up so that the φ ik s ik − k eik ( s ik 2 η ik > 0, k ik > e dik ∞ can be selected in the control design in order to ensure ≤ 0 is held true. . V ik < 0 can be picked up by deducing from (32)-(34).That illustrates the control law can asymptotically stabilize the leader-follower formation system by the derivative and integral terminal sliding mode.Therefore, the follower robots can trace the leader robot with the desired distance and angle steadily.The characteristics of DIT-SMC method are as follows: (1) the convergence time T k can be adjusted by the parameters of the control law; (2) the formation system starts on the derivative and integral terminal sliding surface; (3) the singular problem of T-SMC is avoided; (4) the derivative term can weaken the nonlinear effect.
In (31), the parameter must be assigned as a conservative value to guarantee the formation system stability.From (30), e dik can be exponentially convergent to 0 2×1 by selecting L ik , meaning that κ ik can be very small.Even if κ ik is assigned from a conservative perspective, its value may not be very large.That illustrates the DIT-SMC based NDOB control law protects the formation from the high switching frequency problem, and can substantially alleviate the chattering problem.

Numerical Simulations
Considering the dynamic model of the leader-follower formation system (6), the derivative and integral terminal sliding mode controller is proposed.There are three robots in the leader-follower framework, where the two follower robots track along with the leader robot.The radius of each robot is 0.05 m.The parameter fluctuations in formation system are determined by where i = 1 denotes the leader robot, and k = 2, 3 represent the two follower robots.The uncertainties and external disturbances in lumped model ( 6) are depicted by The parameters in control design are set as γ 12,1 = γ 12,2 = γ 13,1 = γ 13,2 = 4, λ 12,1 = λ 12,2 = λ 13,1 = λ 13,2 = 1, p 11 = p 12 = 9, q 11 = q 12 = 7, q 21 = 3, p 21 = 5, q 22 = 7, p  Figure 3 displays the moving curve of the formation system, which shows the robots are in a line while moving along the circular trajectory.In Figure 3, the solid point denotes the initial position of the formation robots.The arrows are the moving directions of the three robots.It is seen form the Figure 3 that the follower can track the leader robot with the desired distance and angle, while the leader robot tracks the circular trajectory.
In order to provide more insight into the system performance, some comparisons among the SMC method, the second-order SMC, the SMC based NDOB [18], and the DIT-SMC based NDOB are shown in Figures 4-6.The parameters of the sole SMC and SMC based NDOB method are presented in [18], and the parameters of second-order SMC are same as paper [20].The relative distance and relative angular between the leader robot 1 and the two follower robots 2, 3 are displayed in Figure 4. Comparing with the SMC method, the second-order SMC and the SMC based NDOB, the DIT-SMC based NDOB method has shorter convergence time and smoother than the other methods.SMC [18] second-order SMC [18] SMC with NDOB [18] DIT-SMC with NDOB   SMC [18] second-order SMC [18] SMC with NDOB [18] DIT-SMC with NDOB The control inputs of follower robot 3 are shown in the Figure 6, which denotes the acceleration and angular acceleration of follower robot 3. The accelerations of follower robot of follower robot using the three control methods are displayed in Figure 6a, while the angular acceleration using three methods are shown in Figure 6b.From Figures 5 and 6, the combination of the DIT-SMC and NDOB can benefit the decrease of the chattering phenomenon that is an inherent drawback of the SMC methodology.The control inputs of follower robot 3 are shown in the Figure 6, which denotes the acceleration and angular acceleration of follower robot 3. The accelerations of follower robot of follower robot using the three control methods are displayed in Figure 6a, while the angular acceleration using three methods are shown in Figure 6b.From Figures 5 and 6, the combination of the DIT-SMC and NDOB can benefit the decrease of the chattering phenomenon that is an inherent drawback of the SMC methodology.Figure 7 denotes the sliding mode vectors of two follower robots.As proven in the Theorem 1, the reaching time of sliding surface will be eliminated, and the error of formation system will reach to the equilibrium point in the finite time.From Figure 7, the formation system can enter the sliding mode in the beginning, which can guarantee the system stability.8b.From the Figure8, the estimate-error can converge to zero in the finite time.The value of estimate-error is max when t = 0, that is, the maximum is less than 0.5.However, the value of kik is selected 2. Therefore, the system stability can be guaranteed.The control inputs of follower robot 3 are shown in the Figure 6, which denotes the acceleration and angular acceleration of follower robot 3. The accelerations of follower robot of follower robot using the three control methods are displayed in Figure 6a, while the angular acceleration using three methods are shown in Figure 6b.From Figures 5 and 6, the combination of the DIT-SMC and NDOB can benefit the decrease of the chattering phenomenon that is an inherent drawback of the SMC methodology.
Figure 7 denotes the sliding mode vectors of two follower robots.As proven in the Theorem 1, the reaching time of sliding surface will be eliminated, and the error of formation system will reach to the equilibrium point in the finite time.From Figure 7, the formation system can enter the sliding mode in the beginning, which can guarantee the system stability.Figure 7 denotes the sliding mode vectors of two follower robots.As proven in the Theorem 1, the reaching time of sliding surface will be eliminated, and the error of formation system will reach to the equilibrium point in the finite time.From Figure 7, the formation system can enter the sliding mode in the beginning, which can guarantee the system stability.8b.From the Figure8, the estimate-error can converge to zero in the finite time.The value of estimate-error is max when t = 0, that is, the maximum is less than 0.5.However, the value of kik is selected 2. Therefore, the system stability can be guaranteed.8, the estimate-error can converge to zero in the finite time.The value of estimate-error is max when t = 0, that is, the maximum is less than 0.5.However, the value of k ik is selected 2. Therefore, the system stability can be guaranteed.

Conclusions
This paper investigates the formation control problem of multi-robot systems based on the leader-follower mechanism.The leader-follower formation system becomes uncertain because of some adverse effects, such as the parameter fluctuations, external disturbances, and so on.In order to estimate the uncertainties, a control scheme, combining the DIT-SMC and the NDOB, is proposed under the assumption that the uncertainties have an unknown boundary.The stability of the control scheme is proven in the light of Lyapunov theorem.Some simulation results are demonstrated to show the feasibility of the control scheme.

Conclusions
This paper investigates the formation control problem of multi-robot systems based on the leader-follower mechanism.The leader-follower formation system becomes uncertain because of some adverse effects, such as the parameter fluctuations, external disturbances, and so on.In order to estimate the uncertainties, a control scheme, combining the DIT-SMC and the NDOB, is proposed under the assumption that the uncertainties have an unknown boundary.The stability of the control scheme is proven in the light of Lyapunov theorem.Some simulation results are demonstrated to show the feasibility of the control scheme.

Figure 1 .
Figure 1.Sketches of the mobile robot.

Figure 1 .
Figure 1.Sketches of the mobile robot.

Figure 2 .
Figure 2. Sketches of the leader-follower coordinated framework.

Figure 2 .
Figure 2. Sketches of the leader-follower coordinated framework.

Figure 3
Figure3displays the moving curve of the formation system, which shows the robots are in a line while moving along the circular trajectory.In Figure3, the solid point denotes the initial position of the formation robots.The arrows are the moving directions of the three robots.It is seen form the Figure3that the follower can track the leader robot with the desired distance and angle, while the leader robot tracks the circular trajectory.

Figure 5
Figure5denotes the control input of follower robot 2 using different control method.In Figure5a, the acceleration of follower robot 2 are shown, while the angular accelerations of follower robot 2 are displayed in Figure5b.From Figure5a,b, the control input of DIT-SMC based NDOB is smoother than other control methods, which denotes the acceleration and angular acceleration are more stable.

Figure 5
Figure 5 denotes the control input of follower robot 2 using different control method.In Figure 5a, the acceleration of follower robot 2 are shown, while the angular accelerations of follower robot 2 are displayed in Figure 5b.From Figure 5a,b, the control input of DIT-SMC based NDOB is smoother than other control methods, which denotes the acceleration and angular acceleration are more stable.

Figure 5 Figure 5 .
Figure5denotes the control input of follower robot 2 using different control method.In Figure5a, the acceleration of follower robot 2 are shown, while the angular accelerations of follower robot 2 are displayed in Figure5b.From Figure5a,b, the control input of DIT-SMC based NDOB is smoother than other control methods, which denotes the acceleration and angular acceleration are more stable.

Figure 8
Figure 8 illustrates the elements of the estimate-error vectors, where the vectors ed12 are shown in the Figure 8a, and the vectors ed13 are shown in the Figure8b.From the Figure8, the estimate-error can converge to zero in the finite time.The value of estimate-error is max when t = 0, that is, the maximum is less than 0.5.However, the value of kik is selected 2. Therefore, the system stability can be guaranteed.

Figure 8
Figure 8  illustrates the elements of the estimate-error vectors, where the vectors ed12 are shown in the Figure8a, and the vectors ed13 are shown in the Figure8b.From the Figure8, the estimate-error can converge to zero in the finite time.The value of estimate-error is max when t = 0, that is, the maximum is less than 0.5.However, the value of kik is selected 2. Therefore, the system stability can be guaranteed.

Figure 8
Figure 8 illustrates the elements of the estimate-error vectors, where the vectors ed 12 are shown in the Figure 8a, and the vectors ed 13 are shown in the Figure 8b.From the Figure8, the estimate-error can converge to zero in the finite time.The value of estimate-error is max when t = 0, that is, the maximum is less than 0.5.However, the value of k ik is selected 2. Therefore, the system stability can be guaranteed.
∞ ≤ d ik .It means that the unmatched uncertainties of formation system (13) have a boundary.Differentiating the sliding function s ik with the respect to time t, and substituting (13) into the derivative of s ik can get .