Integral sliding mode-based formation control of multiple unertain robots via nonlinear disturbane observer

Abstract

This article proposes a control scheme for formation of maneuvers of a team of mobile robots. The control scheme integrates the integral sliding mode control method with the nonlinear disturbance observer technique. The leader–follower formation dynamics suffer from uncertainties originated from the individual robots. The uncertainties challenge the formation control of such robots. Assuming that the uncertainties are unknown but bounded, an nonlinear disturbance observer-based observer is utilized to approximate them. The observer outputs feed on an integral sliding mode control-based controller. The controller and observer are integrated into the control scheme to realize formation maneuvers despite uncertainties. The formation stability is analyzed by means of the Lyapunov’s theorem. In the sense of Lyapunov, not only the convergence of the approximation errors is guaranteed but also such a control scheme can asymptotically stabilize the formation system. Compared to the results by the sole integral sliding mode control, some simulations are presented to demonstrate the feasibility and performance of the control scheme.

Introduction

Today, robots are changing the way we live and work. ¹ They free human being at dangerous places. Manipulation at disaster sites, nuclear plants and warehouses, is often accomplished with robots. ² To achieve some complex jobs in rescue, manufacturing, and space exploration, several robots have to cooperate with each other such that a multi-robot system consists of them. ³ The multi-robot system can deal with some tasks that are beyond the individual capacities. ⁴ Some potential advantages are often referred, such as efficiency, reliability, and flexibility. In the multi-robot system, each robot plays a unique role. To schedule the robots, coordination of such a system rises up. ⁵

Formation is one of the typical cooperative tasks. The task can be treated as a special case of consensus problems. ^{6 –8} Concerning the multi-robot system, the robots have to form up into some predefined shapes in order to achieve a given formation task. The background of formation originates from the real world. For instant, the robots need to maneuver in some formations when they move in very complex environments.

To achieve multi-robot formations, coordinated control mechanisms should be predefined. Up to now, several kinds of coordinated control mechanisms have been reported. These mechanisms can roughly be classified as behavior-based mechanism, virtual structure mechanism, and leader–follower mechanism. ³ Among the three categories, the category of leader–follower mechanism has become popular because its dynamics are experimentally modelled, but the internal formation stability can be theoretically guaranteed. ⁵ Based on the leader–follower mechanism, various control methods have been developed for multi-robot systems, that is, robust control, ⁹ model predictive control, ¹⁰ feedback linearization control, ¹¹ backstepping control, ¹² decentralized control, ¹³ to name but a few. A systematic review on this topic is presented by Kala. ¹⁴

As a nonlinear feedback design tool, the sliding mode control (SMC) technique ¹⁵ is attractive because of its invariance, meaning that matched uncertainties in an SMC system can be suppressed by the invariance. Some SMC-based methods have been addressed to solve the formation control problem of multi-robot systems, that is, integral SMC, ^5,16 first-order SMC, ¹⁷ terminal SMC, ¹⁸ second-order SMC, ¹⁹ intelligent SMC, ^20,21 and so on. These contributions have verified the feasibility of the SMC methodology for multi-robot formations.

Compared to other SMC methods, the integral SMC design can guarantee an integral SMC system against matched uncertainties from the initial time, indicating that an entire response of the system is of invariance. ⁵ Some successful applications of the integral SMC method have been reported in industries, for example, electric vehicle system, ²² aerospace engineering, ²³ and so on.

Uncertainties exist everywhere. Recall the multi-robot system, where each individual robot may suffer from uncertainties, including but not limited to external disturbances, unmodelled dynamics and parameter perturbations. Affected by the uncertainties, formation dynamics of the multi-robot system become uncertain. Inherently, the formation uncertainties are unmatched such that the invariance of SMC cannot suppress them. In previous works about the SMC-based formations, ^{5,16 –20} the formation uncertainties are considered because they adversely affect the formation stability. Two solutions can be summarized from the works. One solution is to discuss the formation stability by means of graph theory. ^16,21 The other is to analyze the formation stability in light of the Lyapunov’s theorem. ^{5,17 –20} To guarantee the formation stability, both the solutions need to assume that the formation uncertainties are bounded by a known boundary. Unfortunately, the boundary is rather hard to exactly measure or know in reality. The lack of such a boundary may result in severe problems, that is, decrease of the formation robustness, deterioration of the formation performance as far as deficiency of the formation stability. To obtain the important information, it is desired to adaptively approximate the formation uncertainties.

The technique of nonlinear disturbance observer (NDO) has been proven to be effective in handling uncertainties and improving robustness. ²⁴ The applications of NDO have been investigated by some actual cases. ²⁵ Such a technique can be considered as an alternative to attack the issue of formation uncertainties. So far, the academic problem of how to eliminate the adverse effects of the formation uncertainties by NDO still remains unsolved and problematic.

Concerning the formation problem of an uncertain multi-robot system, this article proposes a control scheme that integrates one controller with one observer, where the controller is designed by the integral SMC method and the observer is structured by the NDO technique. One of crucial issues induced by the integration is how to ensure the integral SMC-based formation maneuvers with guaranteed stability in spite of the formation uncertainties. To settle the issue, an integral SMC law of the controller is developed by the Lyapunov’s theorem as well as the estimation errors of the observer. Some results are presented to demonstrate that the formation maneuvers are of asymptotic stability.

Modelling

One single robot

Consider a unicycle-like robot in Figure 1. The round robot is with a radius of r. It is equipped with two driven wheels that are symmetrically mounted on both sides of an axis crossing the center. Each wheel of the robot is independently controlled by one DC motor. Define a vector $\tau ={\left[x\hspace{0.17em}y\hspace{0.17em}\theta \right]}^T$ to describe the robot’s posture, where (x y) are the center in the inertial frame and θ denotes the orientation angle of a mobile frame.

Figure 1.

One single robot.

In a multi-robot system, there are N identical robots displayed in Figure 1. Considering the pure rolling and nonslip conditions, the kinematic model describing the n th robot has the form of

${\dot{\tau }}_n=\left[\begin{array}{c}{\dot{x}}_n\\ {\dot{y}}_n\\ {\dot{\theta }}_n\end{array}\right]=\left[\begin{array}{cc}\cos {\theta }_n& 0\\ \text{sin}{\theta }_n& 0\\ 0& 1\end{array}\right]\cdot \left[\begin{array}{c}{v}_n\\ {\omega }_n\end{array}\right]$

1

The conditions indicate (1) subjects to $\left[\begin{array}{ccc}-\text{sin}{\theta }_n& \cos {\theta }_n& 0\end{array}\right]\cdot {\dot{\tau }}_n=0$ . v_n in equation (1) is the linear velocity of this robot. Accordingly, ω_n means the angular velocity. Taking uncertainties into considerations, the dynamic model of the robot can be depicted by

$\left[\begin{array}{c}{\ddot{x}}_n\\ {\ddot{y}}_n\\ {\ddot{\theta }}_n\end{array}\right]=\left[\begin{array}{c}-{\dot{y}}_n{\dot{\theta }}_n\\ {\dot{x}}_n{\dot{\theta }}_n\\ 0\end{array}\right]+\left[\begin{array}{cc}\cos {\theta }_n& 0\\ \text{sin}{\theta }_n& 0\\ 0& 1\end{array}\right]{\mathbf{u}}_n+\left[\begin{array}{cc}\cos {\theta }_n& 0\\ \text{sin}{\theta }_n& 0\\ 0& 1\end{array}\right]{\Delta }_n{\mathbf{u}}_n+{\pi }_n({\tau }_n,{\dot{\tau }}_n)$

2

In (2), ${\mathbf{u}}_n={\left[{\alpha }_n\hspace{0.17em}{\beta }_n\right]}^T$ is the control inputs vector. Here, ${\alpha }_n=F_n/m_n$ and ${\beta }_n={\sigma }_n/J_n$ are the acceleration and the angular acceleration, respectively. m_n and J_n indicate the mass and moment of inertia of this robot and their values are nominal. F_n and σ_n mean the force and torque. Δ_n indicates the internal uncertainties from the robot, that is, parameter perturbations, described by $\text{diag}\{{\epsilon }_n,\hspace{0.17em}{{\epsilon }^{\prime }}_n\}$ , here ${{\epsilon }^{\prime }}_n$ and ε_n represent the variations in the moment of inertia and the mass, respectively. ${\pi }_n\left({\tau }_n,{\dot{\tau }}_n\right)$ includes all the external disturbances, that is, slipping and skidding, defined by ${\left[{\pi }_{nx}\hspace{0.17em}{\pi }_{ny}{\pi }_{n\theta }\right]}^T$ .

Leader–follower mechanism

Having modelled the n th robot, a formation mechanism of the group should be defined to coordinate the robots and achieve their formation maneuvers. According to the leader–follower mechanism, appoint arbitrary one and name it leader. The remaining N − 1 robots are titled followers. Such a hierarchical mechanism creates N − 1 pairs, where the sole leader traces a predefined trajectory and all the followers takes charge of tracking the leader. To form up into a given formation, each follower has to keep track of its desired distance and relative angle from the leader. Figure 2 illustrates the mechanism by a pair, where indexes i and k represent the leader and the follower, respectively.

Figure 2.

A leader–follower pair.

Given a receiver located at the leader’s center and a sensor equipped at the follower’s front castor, the distance l_ik is defined by $l_{ik}=\sqrt{{\left(x_i-x_k-r\cos {\theta }_k\right)}^2+{\left(y_i-y_k-r\text{sin}{\theta }_k\right)}^2}$ and the relative angle ψ_ik is determined by ${\psi }_{ik}=\pi -{\theta }_i+{\zeta }_{ik}$ , where ζ_ik has a form of ${\zeta }_{ik}=\arctan \frac{y_i-y_k-r\text{sin}{\theta }_k}{x_i-x_k-r\cos {\theta }_k}$ .

The individual uncertainties in equation (2) have unfavorable effects on the leader–follower mechanism, leading to the formation uncertainties. Such derivative uncertainties adversely affect the formation stability. Owning to their randomness and complexity, the formation uncertainties are rather hard to be modelled so that they challenge the formation control design very much. Motivated by the challenge, the article develops a control scheme under a mild assumption that the formation uncertainties are bounded by an unknown boundary. In order to concentrate upon the development, some ideal conditions in the pair are taken into considerations: (a) no collision, (b) no communication delay, and (c) the follower knows the leader’s posture as well as its own.

The leader–follower mechanism in Figure 2 can be modelled by differentiating l_ik and ψ_ik twice with respect to time. Finally, the formation dynamics can have the form of

${\dot{\mathbf{x}}}_{ik}={\mathbf{A}}_{ik}{\mathbf{x}}_{ik}+{\mathbf{B}}_{ik\mathrm{,1}}{\mathbf{u}}_k+{\mathbf{B}}_{ik\mathrm{,2}}{\mathbf{d}}_{ik}\left({\mathbf{x}}_{ik},{\tau }_i,\hspace{0.17em}{\tau }_k,\hspace{0.17em}{\dot{\tau }}_i,\hspace{0.17em}{\dot{\tau }}_k\right)$

3

In equation (3), ${\mathbf{\text{x}}}_{ik}={\left[x_1\hspace{0.17em}{x}_2\hspace{0.17em}{x}_3\hspace{0.17em}{x}_4\right]}^T$ where $x_1=l_{ik}$ , $x_2={\dot{l}}_{ik}$ , $x_3={\psi }_{ik}$ , and $x_4={\dot{\psi }}_{ik}$ . Moreover, ${\mathbf{A}}_{ik}=\left[\begin{array}{cccc}0& 1& 0& 0\\ 0& 0& 0& 0\\ 0& 0& 0& 1\\ 0& 0& 0& 0\end{array}\right]$ , ${\mathbf{B}}_{ik\mathrm{,1}}=\left[\begin{array}{cc}0& 0\\ \cos {\varphi }_{ik}& r\text{sin}{\varphi }_{ik}\\ 0& 0\\ -\frac{\text{sin}{\varphi }_{ik}}{l_{ik}}& \frac{r\cos {\varphi }_{ik}}{l_{ik}}\end{array}\right]$ , ${\mathbf{B}}_{ik\mathrm{,2}}={\left[\begin{array}{cccc}0& 1& 0& 0\\ 0& 0& 0& 1\end{array}\right]}^T$ , and ${\mathbf{d}}_{ik}\left({\mathbf{\text{x}}}_{ik},\hspace{0.17em}{\tau }_i,\hspace{0.17em}{\tau }_k,\hspace{0.17em}{\dot{\tau }}_i,\hspace{0.17em}{\dot{\tau }}_k\right)$ abbreviated by d _ik is given by ${\mathbf{G}}_{ik}{\Delta }_k{\mathbf{u}}_k+{\mathbf{L}}_{ik}\left({\mathbf{I}}_2+{\Delta }_i\right){\mathbf{u}}_i+{\mathbf{F}}_{ik}+{\mathbf{P}}_{ik}$ , where I₂ denotes a 2 × 2 identity matrix, ${\mathbf{G}}_{ik}=\left[\begin{array}{cc}\cos {\varphi }_{ik}& r\text{sin}{\varphi }_{ik}\\ \frac{-\text{sin}{\varphi }_{ik}}{l_{ik}}& \frac{r\cos {\varphi }_{ik}}{l_{ik}}\end{array}\right]$ , ${\mathbf{L}}_{ik}=\left[\begin{array}{cc}-\cos {\psi }_{ik}& 0\\ \frac{\text{sin}{\psi }_{ik}}{l_{ik}}& -1\end{array}\right]$ , ${\varphi }_{ik}={\psi }_{ik}+{\theta }_{ik}$ , ${\theta }_{ik}={\theta }_i-{\theta }_k$ , ${\mathbf{F}}_{ik}={\left[F_1\hspace{0.17em}{F}_2\right]}^T$ , ${\mathbf{P}}_{ik}={\left[P_1\hspace{0.17em}{P}_2\right]}^T$ , F₁, F₂, P₁, and P₂ are formulated by

$F_1={\left({\dot{\psi }}_{ik}\right)}^2{l}_{ik}+2{\dot{\psi }}_{ik}{\dot{\theta }}_i{l}_{ik}+{\left({\dot{\theta }}_i\right)}^2{l}_{ik}-r\cos {\varphi }_{ik}{\left({\dot{\theta }}_k\right)}^2-\left({\dot{y}}_k{\dot{\theta }}_k-{\dot{y}}_i{\dot{\theta }}_i\right)\cos \left({\psi }_{ik}+{\theta }_i\right)-\left({\dot{x}}_i{\dot{\theta }}_i-{\dot{x}}_k{\dot{\theta }}_k\right)\text{sin}\left({\psi }_{ik}+{\theta }_i\right)$

$\begin{array}{ll}{F}_2=\hfill & \frac{-\left({\dot{y}}_k{\dot{\varphi }}_{ik}-{\dot{\psi }}_{ik}{\dot{y}}_i\right)\text{sin}\left({\psi }_{ik}+{\theta }_i\right)-\left({\dot{x}}_k{\dot{\varphi }}_{ik}-{\dot{\psi }}_{ik}{\dot{x}}_i\right)\cos \left({\psi }_{ik}+{\theta }_i\right)-r{\dot{\theta }}_k{\dot{\varphi }}_{ik}\text{sin}{\varphi }_{ik}}{l_{ik}}\hfill \\ \hspace{0.17em}\hfill & +\frac{{\dot{l}}_{ik}\left[\left({\dot{y}}_i-{\dot{y}}_k\right)\cos \left({\psi }_{ik}+{\theta }_i\right)-\left({\dot{x}}_i-{\dot{x}}_k\right)\text{sin}\left({\psi }_{ik}+{\theta }_i\right)-r{\dot{\theta }}_k\cos {\varphi }_{ik}\right]}{l_{ik}{\hspace{0.17em}}^2}\hfill \end{array}$

$P_1=-\left({\pi }_{ix}-{\pi }_{kx}\right)\cos \left({\psi }_{ik}+{\theta }_i\right)-\left({\pi }_{iy}-{\pi }_{ky}\right)\text{sin}\left({\psi }_{ik}+{\theta }_i\right)+r{\pi }_{k\theta }\text{sin}{\varphi }_{ik}$

$P_2=\frac{\left({\pi }_{ix}-{\pi }_{kx}\right)\text{sin}\left({\psi }_{ik}+{\theta }_i\right)-\left({\pi }_{iy}-{\pi }_{ky}\right)\cos \left({\psi }_{ik}+{\theta }_i\right)+r{\pi }_{k\theta }\cos {\varphi }_{ik}-l_{ik}{\pi }_{i\theta }}{l_{ik}}$

From equation (3), $\text{rank}\left[{\mathbf{B}}_{ik\mathrm{,1}}{\mathbf{A}}_{ik}{\mathbf{B}}_{ik\mathrm{,1}}\right]=4$ is of full rank, indicating that equation (3) is controllable so that the poles of the closed-loop formation system can be arbitrarily assigned. d _ik aggregates all the formation uncertainties and is inherently unmatched. It is mild that d _ik is hypothetically bounded by a unknown boundary. Unfortunately, the assumption contradicts the formation control design because the formation stability depends on the value of boundary. Usually, a conservative boundary can be empirically predefined to guarantee the formation stability. However, such a conservative boundary may adversely affect the formation performance. Especially, a drawback is that the chattering phenomenon may be induced. To fill the gap between the unknown boundary and the formation performance, the NDO technique is employed so that such an observer serves as a compensator to approximate d _ik .

Formation design

Nonlinear disturbance observer

Assumption 1

• ${\Vert {\mathbf{d}}_{ik}\Vert }_{\infty }⩽d_{ik}^{*}$ , where $d_{ik}^{*}>0$ is constant but unknown and ${\Vert \cdot \Vert }_{\infty }$ denotes ∞-norm.

Assumption 2

• d _ik is slowly time-varying, indicating ${\dot{\mathbf{d}}}_{ik}\simeq {\mathbf{0}}_2{=[00]}^T$ .

Consider the formation dynamics (3) and design the following NDO-based observer (4). ²⁴

$\{\begin{array}{l}{\dot{\mathbf{p}}}_{ik}=-{\mathbf{L}}_{ik}{\mathbf{B}}_{ik\mathrm{,2}}{\mathbf{p}}_{ik}-{\mathbf{L}}_{ik}\left({\mathbf{B}}_{ik\mathrm{,2}}{\mathbf{L}}_{ik}{\mathbf{x}}_{ik}+{\mathbf{A}}_{ik}{\mathbf{x}}_{ik}+{\mathbf{B}}_{ik\mathrm{,1}}{\mathbf{u}}_k\right)\hfill \\ {\widehat{\mathbf{d}}}_{ik}={\mathbf{p}}_{ik}+{\mathbf{L}}_{ik}{\mathbf{x}}_{ik}\hfill \end{array}$

4

In (4), ${\mathbf{p}}_{ik}\in {\Re }^{2\times 1}$ is the internal states of the observer, the gain matrix ${\mathbf{L}}_{ik}\in {\Re }^{2\times 4}$ is designed such that all the eigenvalues of H _ik are positive (Here, ${\mathbf{H}}_{ik}={\mathbf{L}}_{ik}{\mathbf{B}}_{ik\mathrm{,2}}$ ) and ${\widehat{\mathbf{d}}}_{ik}$ is the approximation of d _ik .

Define an error vector ${\mathbf{e}}_{dik}={\mathbf{d}}_{ik}-{\widehat{\mathbf{d}}}_{ik}$ , differentiate the error vector with respect to time, substitute equation (4) into the derivative of e _dik , and take assumptions 1 and 2 into account. We can obtain equation (5).

$\begin{array}{c}{\dot{\mathbf{e}}}_{dik}={\dot{\mathbf{d}}}_{ik}-{\dot{\widehat{\mathbf{d}}}}_{ik}\simeq -{\dot{\mathbf{p}}}_{ik}-{\mathbf{L}}_{ik}{\dot{\mathbf{x}}}_{ik}\\ ={\mathbf{H}}_{ik}{\mathbf{\text{p}}}_{ik}+{\mathbf{L}}_{ik}\left({\mathbf{B}}_{ik\mathrm{,2}}{\mathbf{L}}_{ik}{\mathbf{\text{x}}}_{ik}+{\mathbf{A}}_{ik}{\mathbf{\text{x}}}_{ik}+{\mathbf{B}}_{ik\mathrm{,1}}{\mathbf{\text{u}}}_k\right)-{\mathbf{L}}_{ik}\left({\mathbf{A}}_{ik}{\mathbf{\text{x}}}_{ik}+{\mathbf{B}}_{ik\mathrm{,1}}{\mathbf{\text{u}}}_k+{\mathbf{B}}_{ik\mathrm{,2}}{\mathbf{d}}_{ik}\right)\\ ={\mathbf{H}}_{ik}\left({\widehat{\mathbf{d}}}_{ik}-{\mathbf{L}}_{ik}{\mathbf{\text{x}}}_{ik}\right)+{\mathbf{L}}_{ik}\left({\mathbf{B}}_{ik\mathrm{,2}}{\mathbf{L}}_{ik}{\mathbf{\text{x}}}_{ik}+{\mathbf{A}}_{ik}{\mathbf{\text{x}}}_{ik}+{\mathbf{B}}_{ik\mathrm{,1}}{\mathbf{\text{u}}}_k\right)-{\mathbf{L}}_{ik}\left({\mathbf{A}}_{ik}{\mathbf{\text{x}}}_{ik}+{\mathbf{B}}_{ik\mathrm{,1}}{\mathbf{\text{u}}}_k+{\mathbf{B}}_{ik\mathrm{,2}}{\mathbf{d}}_{ik}\right)\\ ={\mathbf{H}}_{ik}\left({\widehat{\mathbf{d}}}_{ik}-{\mathbf{d}}_{ik}\right)=-{\mathbf{H}}_{ik}{\mathbf{e}}_{dik}\end{array}$

5

The solution of equation (5) is ${\mathbf{e}}_{dik}=\exp \left(-{\mathbf{H}}_{ik}t\right){\mathbf{e}}_{dik}\left(0\right)$ , where ${\mathbf{e}}_{dik}\left(0\right)$ is the initial conditions at t = 0. Owing to ${\lambda }_{\text{min}}\left({\mathbf{H}}_{ik}\right)>0$ (here ${\lambda }_{\text{min}}\left({\mathbf{H}}_{ik}\right)$ is the minimum eigenvalue of $\left({\mathbf{H}}_{ik}\right)$ , this fact indicates that the error vector e _dik is exponentially convergent to 0₂ as t → ∞ (better performance than asymptotic stability).

Integral sliding mode controller

To coordinate the pair in Figure 2 by the integral SMC method, a nominal integral sliding surface vector (6) has to be developed at first.

${\delta }_{ik}={\mathbf{C}}_{ik}\left[\left({\mathbf{x}}_{ik}-{\mathbf{x}}_{ik}^0\right)-{\displaystyle {\int }_0^t}\left({\mathbf{A}}_{ik}-{\mathbf{B}}_{ik\mathrm{,1}}{\mathbf{K}}_{ik}\right)\cdot {\mathbf{x}}_{ik}\left(\varsigma \right)\text{d}\varsigma \right]$

6

where ${\delta }_{ik}\in {\Re }^{2\times 1}$ defined by ${\left[{\delta }_{ik\mathrm{,1}\times 1}\hspace{0.17em}{\delta }_{ik\mathrm{,2}\times 1}\right]}^T$ denotes the nominal integral sliding surfaces. ${\mathbf{x}}_{ik}^0$ means the initial states at t = 0; ${\mathbf{C}}_{ik}\in {\Re }^{2\times 4}$ is a positive constant matrix and K _ik is calculated by pole assignment to make all the eigenvalues of matrix $\left({\mathbf{A}}_{ik}-{\mathbf{B}}_{ik\mathrm{,1}}{\mathbf{K}}_{ik}\right)$ rigidly located at the left-half complex plane.

Differentiate the vector δ_ik in equation (6) with respect to time, substitute equation (3) into the derivative of δ_ik and ${\dot{\delta }}_{ik}$ becomes ${\dot{\delta }}_{ik}={\mathbf{C}}_{ik}{\mathbf{B}}_{ik\mathrm{,1}}{\mathbf{K}}_{ik}{\mathbf{x}}_{ik}+{\mathbf{C}}_{ik}{\mathbf{B}}_{ik\mathrm{,1}}{\mathbf{u}}_k+{\mathbf{C}}_{ik}{\mathbf{B}}_{ik\mathrm{,2}}{\mathbf{d}}_{ik}$ . Suggested by Utkin, ¹⁵ the equivalent control law u_keq independently governs the integral SMC system on the sliding mode stage. The expression of u_keq can be derived from ${\dot{\delta }}_{ik}=0_2$ , formulated by $-{\left({\mathbf{C}}_{ik}{\mathbf{B}}_{ik\mathrm{,1}}\right)}^{-1}\left[{\mathbf{C}}_{ik}{\mathbf{B}}_{ik\mathrm{,1}}{\mathbf{K}}_{ik}{\mathbf{x}}_{ik}+{\mathbf{C}}_{ik}{\mathbf{B}}_{ik\mathrm{,2}}{\mathbf{d}}_{ik}\right]$ . Since d _ik is completely unknown except its boundedness, the unknown term d _ik has to be replaced by its approximation ${\widehat{\mathbf{d}}}_{ik}$ . Finally, u _keq is rearranged by

${\mathbf{u}}_{k\text{eq}}=-{\left({\mathbf{C}}_{ik}{\mathbf{B}}_{ik\mathrm{,1}}\right)}^{-1}\left[{\mathbf{C}}_{ik}{\mathbf{B}}_{ik\mathrm{,1}}{\mathbf{K}}_{ik}{\mathbf{x}}_{ik}+{\mathbf{C}}_{ik}{\mathbf{B}}_{ik\mathrm{,2}}{\widehat{\mathbf{d}}}_{ik}\right]$

7

Concerning the uncertain term, a real sliding mode surfaces vector is defined as

${\mathbf{s}}_{ik}={\delta }_{ik}+{\widehat{\mathbf{d}}}_{ik}={\mathbf{C}}_{ik}\left[\left({\mathbf{x}}_{ik}-{\mathbf{x}}_{ik}^0\right)-{\displaystyle {\int }_0^t}\left({\mathbf{A}}_{ik}-{\mathbf{B}}_{ik\mathrm{,1}}{\mathbf{K}}_{ik}\right)\cdot {\mathbf{x}}_{ik}\left(\varsigma \right)\text{d}\varsigma \right]+{\widehat{\mathbf{d}}}_{ik}$

8

where ${\mathbf{s}}_{ik}\in {\Re }^{2\times 1}$ defined by ${\left[s_{ik\mathrm{,1}\times 1}\hspace{0.17em}{s}_{ik\mathrm{,2}\times 1}\right]}^T$ denotes the real integral sliding surfaces.

According to the designed control scheme, the control law is defined by ${\mathbf{u}}_k={\mathbf{u}}_{k\text{eq}}+{\mathbf{u}}_{k\text{sw}}$ , where ${\mathbf{u}}_{k\text{sw}}$ is the switching control law. Design ${\mathbf{u}}_{k\text{sw}}=-{\left({\mathbf{C}}_{ik}{\mathbf{B}}_{ik\mathrm{,1}}\right)}^{-1}\left[{\overline{\kappa }}_{ik}\text{sign}\left({\mathbf{s}}_{ik}\right)+{\overline{\eta }}_{ik}{\mathbf{s}}_{ik}\right]$ , where $\overline{\kappa }$ and $\overline{\eta }$ are positive constants and $\text{sign}\left({\mathbf{s}}_{ik}\right)\in {\Re }^{2\times 1}$ defined by ${\left[\text{sign}\left(s_{ik\mathrm{,1}\times 1}\right)\text{sign}\left(s_{ik\mathrm{,2}\times 1}\right)\right]}^T$ is the signum functions vector. Finally, the control law u _k can be formulated by

${\mathbf{u}}_k=-{\left({\mathbf{C}}_{ik}{\mathbf{B}}_{ik\mathrm{,1}}\right)}^{-1}\left[{\mathbf{C}}_{ik}{\mathbf{B}}_{ik\mathrm{,1}}{\mathbf{K}}_{ik}{\mathbf{x}}_{ik}+{\mathbf{J}}_{ik}{\widehat{\mathbf{d}}}_{ik}+{\overline{\kappa }}_{ik}\text{sign}\left({\mathbf{s}}_{ik}\right)+{\overline{\eta }}_{ik}{\mathbf{s}}_{ik}\right]$

9

where ${\mathbf{J}}_{ik}={\mathbf{C}}_{ik}{\mathbf{B}}_{ik\mathrm{,2}}$ .

System structure

The system structure is presented in Figure 3. From Figure 3, the state variables, composed of l_ik, ψ_ik, and their first derivatives, are located at the feedback channel. They feed on the NDO-based observer to estimate the approximation ${\widehat{\mathbf{d}}}_{ik}$ . Furthermore, the state variables construct the vector s _ik . The designed control scheme employs s _ik , x _ik , and ${\widehat{\mathbf{d}}}_{ik}$ to generate the final control inputs vector u _k . Finally, u _k is applied to the follower k in order to maintain the desired distance and relative angle with respect to the leader i during formation maneuvers.

Figure 3.

Structure of the control scheme.

System stability

Assumption 3

• ${\Vert e_{dik}\Vert }_{\infty }⩽e_{dik}^{*}$ , where $e_{dik}^{*}>0$ is constant but unknown.

Theorem 1

Consider the formation dynamics (3), take assumptions 2 and 3 into account, design the NDO-based observer (4), define the real sliding surfaces (8), and utilize the integral SMC-based control law (9). The closed-loop formation control system structured by Figure 3 is asymptotically stable if ${\overline{\kappa }}_{ik}>\left||{\mathbf{H}}_{ik}+{\mathbf{J}}_{ik}\right|{|}_{\infty }{e}_{dik}^{*}$ .

Proof

Differentiate s _ik in equation (8) with respect to time. Substituting equation (3) into the derivative of s _ik yields

${\dot{\mathbf{s}}}_{ik}={\mathbf{C}}_{ik}{\mathbf{B}}_{ik\mathrm{,1}}{\mathbf{K}}_{ik}{\mathbf{x}}_{ik}+{\mathbf{C}}_{ik}{\mathbf{B}}_{ik\mathrm{,1}}{\mathbf{u}}_k+{\mathbf{J}}_{ik}{\mathbf{d}}_{ik}+{\dot{\widehat{\mathbf{d}}}}_{ik}$

10

Considering the control law (9), equation (10) becomes

${\dot{\mathbf{s}}}_{ik}=-{\overline{\kappa }}_{ik}\text{sign}\left({\mathbf{s}}_{ik}\right)-{\overline{\eta }}_{ik}{\mathbf{s}}_{ik}+{\mathbf{J}}_{ik}\left({\mathbf{d}}_{ik}-{\widehat{\mathbf{d}}}_{ik}\right)+{\dot{\widehat{\mathbf{d}}}}_{ik}$

11

From equation (5), ${\dot{\mathbf{e}}}_{dik}(t)+{\mathbf{H}}_{ik}{\mathbf{e}}_{dik}(t)=0$ holds true which means ${\dot{\widehat{\mathbf{d}}}}_{ik}=-{\mathbf{H}}_{ik}\left({\widehat{\mathbf{d}}}_{ik}-{\mathbf{d}}_{ik}\right)$ . Finally, ${\dot{\mathbf{s}}}_{ik}$ has a form of

${\dot{\mathbf{s}}}_{ik}=-{\overline{\kappa }}_{ik}\text{sign}\left({\mathbf{s}}_{ik}\right)-{\overline{\eta }}_{ik}{\mathbf{s}}_{ik}+\left({\mathbf{J}}_{ik}+{\mathbf{H}}_{ik}\right)\left({\mathbf{d}}_{ik}-{\widehat{\mathbf{d}}}_{ik}\right)$

12

Select a Lyapunov candidate function $V(t)={\Vert {\mathbf{s}}_{ik}\Vert }_2$ where ${\Vert \cdot \Vert }_2$ means 2-norm. Differentiate $V(t)$ with respect to time. The derivative of V can be written by $\dot{V}(t)=\frac{{\mathbf{s}}_{ik}^T{\dot{\mathbf{s}}}_{ik}}{{\Vert {\mathbf{s}}_{ik}\Vert }_2}$ . Replace $\dot{\mathbf{s}}$ in the derivative of V by equation (12). Equation (13) can be obtained

$\begin{array}{c}\dot{V}\left(t\right)=\frac{{\mathbf{s}}_{ik}^T}{{\Vert {\mathbf{s}}_{ik}\Vert }_2}\left[-{\overline{\kappa }}_{ik}\text{sign}\left({\mathbf{s}}_{ik}\right)-{\overline{\eta }}_{ik}{\mathbf{s}}_{ik}+\left({\mathbf{J}}_{ik}+{\mathbf{H}}_{ik}\right)\left({\mathbf{d}}_{ik}-{\widehat{\mathbf{d}}}_{ik}\right)\right]=-{\overline{\kappa }}_{ik}\frac{{\mathbf{s}}_{ik}^T\text{sign}\left({\mathbf{s}}_{ik}\right)}{{\Vert {\mathbf{s}}_{ik}\Vert }_2}-{\overline{\eta }}_{ik}\frac{{\mathbf{s}}_{ik}^T{\mathbf{s}}_{ik}}{{\Vert {\mathbf{s}}_{ik}\Vert }_2}+\frac{{\mathbf{s}}_{ik}^T}{{\Vert {\mathbf{s}}_{ik}\Vert }_2}\left({\mathbf{J}}_{ik}+{\mathbf{H}}_{ik}\right)\left({\mathbf{d}}_{ik}-{\widehat{\mathbf{d}}}_{ik}\right)\\ =-{\overline{\kappa }}_{ik}\frac{{\Vert {\mathbf{s}}_{ik}\Vert }_1}{{\Vert {\mathbf{s}}_{ik}\Vert }_2}-{\overline{\eta }}_{ik}\frac{{\Vert {\mathbf{s}}_{ik}\Vert }_2^2}{{\Vert {\mathbf{s}}_{ik}\Vert }_2}+\frac{{\mathbf{s}}_{ik}^T}{{\Vert {\mathbf{s}}_{ik}\Vert }_2}\left({\mathbf{J}}_{ik}+{\mathbf{H}}_{ik}\right){\mathbf{e}}_{dik}⩽-{\overline{\kappa }}_{ik}\frac{{\Vert {\mathbf{s}}_{ik}\Vert }_1}{{\Vert {\mathbf{s}}_{ik}\Vert }_2}-{\overline{\eta }}_{ik}\frac{{\Vert {\mathbf{s}}_{ik}\Vert }_2^2}{{\Vert {\mathbf{s}}_{ik}\Vert }_2}+\frac{{\mathbf{s}}_{ik}^T}{{\Vert {\mathbf{s}}_{ik}\Vert }_2}\left[\begin{array}{c}{\Vert \left({\mathbf{J}}_{ik}+{\mathbf{H}}_{ik}\right){\mathbf{e}}_{dik}\Vert }_{\infty }\\ {\Vert \left({\mathbf{J}}_{ik}+{\mathbf{H}}_{ik}\right){\mathbf{e}}_{dik}\Vert }_{\infty }\end{array}\right]\\ ⩽-{\overline{\kappa }}_{ik}\frac{{\Vert {\mathbf{s}}_{ik}\Vert }_1}{{\Vert {\mathbf{s}}_{ik}\Vert }_2}-{\overline{\eta }}_{ik}\frac{{\Vert {\mathbf{s}}_{ik}\Vert }_2^2}{{\Vert {\mathbf{s}}_{ik}\Vert }_2}+\frac{{\Vert {\mathbf{s}}_{ik}\Vert }_1}{{\Vert {\mathbf{s}}_{ik}\Vert }_2}{\Vert \left({\mathbf{J}}_{ik}+{\mathbf{H}}_{ik}\right){\mathbf{e}}_{dik}\Vert }_{\infty }⩽-{\overline{\kappa }}_{ik}\frac{{\Vert {\mathbf{s}}_{ik}\Vert }_1}{{\Vert {\mathbf{s}}_{ik}\Vert }_2}-{\overline{\eta }}_{ik}\frac{{\Vert {\mathbf{s}}_{ik}\Vert }_2^2}{{\Vert {\mathbf{s}}_{ik}\Vert }_2}+\frac{{\Vert {\mathbf{s}}_{ik}\Vert }_1}{{\Vert {\mathbf{s}}_{ik}\Vert }_2}{\Vert \left({\mathbf{J}}_{ik}+{\mathbf{H}}_{ik}\right)\Vert }_{\infty }{\mathbf{e}}_{dik}^{*}\\ =-\left[{\overline{\kappa }}_{ik}-{\Vert \left({\mathbf{J}}_{ik}+{\mathbf{H}}_{ik}\right)\Vert }_{\infty }{\mathbf{e}}_{dik}^{*}\right]\frac{{\Vert {\mathbf{s}}_{ik}\Vert }_1}{{\Vert {\mathbf{s}}_{ik}\Vert }_2}-{\overline{\eta }}_{ik}\frac{{\Vert {\mathbf{s}}_{ik}\Vert }_2^2}{{\Vert {\mathbf{s}}_{ik}\Vert }_2}\end{array}$

13

where ${\Vert \cdot \Vert }_1$ means 1-norm.

Pick up ${\overline{\kappa }}_{ik}>\left||{\mathbf{H}}_{ik}+{\mathbf{J}}_{ik}\right|{|}_{\infty }{e}_{dik}^{*}$ such that $\dot{V}\left(t\right)<0$ exists. Concerning $V\left(t\right)\ge 0$ , the closed-loop formation system is asymptotically stable in the sense of Lyapunov.

Considering assumption 3, $e_{dik}^{*}$ is unknown, meaning that it is hard to determine $\overline{\kappa }$ . To guarantee the system stability, a conservative value of $e_{dik}^{*}$ should be assigned. On this aspect, there seem no benefits earned from such an NDO-based integral SMC scheme. However, ${\mathbf{e}}_{ik}$ originated from the scheme is exponentially convergent, indicating that a small value of $e_{dik}^{*}$ can be chosen. According to theorem 2, ${\overline{\kappa }}_{ik}$ can be decreased to improve the chattering phenomenon. Thus, the presented control scheme can improve the formation performance.

Simulation results

To demonstrate the performance of the integral SMC-based control scheme, this section implements some numerical simulations on a multi-robot platform and discusses the results. The platform is composed of three mobile robots. The robots are identical (see Figure 1) with radius 0.05 m, where one leader is numbered by 1 and two followers are numbered by 2 and 3, respectively.

The leader in the formation mechanism is to take charge of tracking its predefined trajectories; its control can be treated as tracking control of a single robot and it is equipped with a kinematic-based sliding mode controller. ²⁶ The designed control scheme is applied to the other two followers in order that the two robots can maintain the desired distances and the desired relative angles.

Δ_i and Δ_k in equation (3) are considered as $\text{diag}\{0.3\text{rand}()-0.15\hspace{0.17em}0.3\text{rand}()-0.15\}$ , where i = 1 delegates the leader, k = 2, 3 means the followers, and rand() is a MATLAB R2010a function to generate a random number. Such consideration indicates that the largest parameter variations are up to 15% of the nominal ones. P _ik contains ${\pi }_{ix}-{\pi }_{kx}$ and ${\pi }_{iy}-{\pi }_{ky}$ in equation (3) which shows that it is not representative to define ${\pi }_i={\pi }_k$ . Here, we consider ${\pi }_{ix}={\pi }_{iy}={\pi }_{i\theta }=0.5\text{sin}\left(2\text{π}t\right)$ and ${\pi }_{kx}={\pi }_{ky}={\pi }_{k\theta }=0.8\text{sin}\left(\text{π}t\right)$ .

The parameter matrices K _ik and C _ik of the real sliding surface vector s _ik in equation (8) have to be predefined. K _ik needs to be designed such that the eigenvalues of matrix ${\mathbf{A}}_{ik}-{\mathbf{B}}_{ik}{\mathbf{K}}_{ik}$ in equation (8) are those in the specified vector ${\left[\hspace{0.17em}-1.7500\pm 0.6614i\hspace{0.17em}-1.4500\pm 0.8930i\hspace{0.17em}\right]}^T$ . Finally, K _ik is determined by ${\mathbf{K}}_{ik}=\left[\begin{array}{cccc}\text{3.5cos}\left({\varphi }_{ik}\right)& \text{3.5cos}\left({\varphi }_{ik}\right)& -\text{2.9}{l}_{ik}\text{sin}\left({\varphi }_{ik}\right)& -\text{2.9}{l}_{ik}\text{sin}\left({\varphi }_{ik}\right)\\ \frac{\text{3.5sin}\left({\varphi }_{ik}\right)}{d}& \frac{\text{3.5sin}\left({\varphi }_{ik}\right)}{d}& \frac{\text{2.9}{l}_{ik}\text{cos}\left({\varphi }_{ik}\right)}{d}& \frac{\text{2.9}{l}_{ik}\text{cos}\left({\varphi }_{ik}\right)}{d}\end{array}\right]$ . The other matrix C _ik is determined by $\left[\begin{array}{cccc}120& 88& 0& 0\\ 0& 0& 90& 30\end{array}\right]$ . The two parameters ${\overline{\kappa }}_{ik}$ and ${\overline{\eta }}_{ik}$ of the switching control law are determined by 1 and 15, respectively. Concerning the observer, its gain matrix is chosen as $\left[\begin{array}{cccc}0& 15& 0& 18\\ 0& 15& 0& 20\end{array}\right]$ by trial and error.

Consider the following formation task. The robots move along straight lines and form up into a string formation. Their formation maneuvers take place at t = 8 s. Thereafter, the alternative formation becomes triangular. The initial postures of the three mobile robots are located at ${\tau }_1={\left[0\hspace{0.17em}\text{m}\hspace{0.17em}0\hspace{0.17em}\text{m}\hspace{0.17em}0\hspace{0.17em}\text{rad}\right]}^T$ , ${\tau }_2={\left[0\hspace{0.17em}\text{m}\hspace{0.17em}0.4\hspace{0.17em}\text{m}\hspace{0.17em}0\hspace{0.17em}\text{rad}\right]}^T$ , and ${\tau }_3={\left[0\hspace{0.17em}\text{m}\hspace{0.17em}-0.3\hspace{0.17em}\text{m}\hspace{0.17em}0\hspace{0.17em}\text{rad}\right]}^T$ , respectively. Accordingly, the initial states of the followers can be calculated by ${\mathbf{x}}_{12}^0={\left[0.4\hspace{0.17em}\text{m}\hspace{0.17em}0\hspace{0.17em}\text{m}/\text{s}\hspace{0.17em}\frac{\pi }{2}\text{rad}\hspace{0.17em}0\hspace{0.17em}\text{rad}/\text{s}\right]}^T$ and ${\mathbf{x}}_{13}^0={\left[0.3\hspace{0.17em}\text{m}\hspace{0.17em}0\hspace{0.17em}\text{m}/\text{s}\hspace{0.17em}-\frac{\pi }{2}\text{rad}\hspace{0.17em}0\hspace{0.17em}\text{rad}/\text{s}\right]}^T$ . To achieve the given formation maneuvers, the desired states during the first 8 s are ${\mathbf{x}}_{12}^d={\left[0.2\hspace{0.17em}\text{m}\hspace{0.17em}0\hspace{0.17em}\text{m}/\text{s}\hspace{0.17em}\frac{\pi }{2}\text{rad}\hspace{0.17em}0\hspace{0.17em}\text{rad}/\text{s}\right]}^T$ and ${\mathbf{x}}_{13}^d={\left[0.15\hspace{0.17em}\text{m}\hspace{0.17em}0\hspace{0.17em}\text{m}/\text{s}\hspace{0.17em}-\frac{\pi }{2}\text{rad}\hspace{0.17em}0\hspace{0.17em}\text{rad}/\text{s}\right]}^T$ . The desired states during the next 8 s are ${\mathbf{x}}_{12}^d={\left[0.5\hspace{0.17em}\text{m}\hspace{0.17em}0\hspace{0.17em}\text{m}/\text{s}\hspace{0.17em}\frac{2\pi }{3}\text{rad}\hspace{0.17em}0\hspace{0.17em}\text{rad}/\text{s}\right]}^T$ and ${\mathbf{x}}_{13}^d={\left[0.5\hspace{0.17em}\text{m}\hspace{0.17em}0\hspace{0.17em}\text{m}/\text{s}\hspace{0.17em}-\frac{2\pi }{3}\text{rad}\hspace{0.17em}0\hspace{0.17em}\text{rad}/\text{s}\right]}^T$ .

Figure 4 displays the group of mobile robots forms up into a triangular formation from a string one while moving in straight lines, where those filled circles denote the initial positions of the robots and those unfilled circles indicate the robots’ positions in the dynamic process. In each circle, the arrow means the orientation angle of that robot. To demonstrate their postures and formation, the solid lines in Figure 4 bond the robots together at the same moment. In the simulation, the desired linear velocity $v_1^d$ and the desired angular velocity ${\omega }_1^d$ of the leader are set by $0.5\hspace{0.17em}\text{m}\cdot {\text{s}}^{-1}$ and $0\hspace{0.17em}\text{rad}\cdot {\text{s}}^{-1}$ , respectively. From Figure 4, the two followers can follow the leader and the given formation task has been achieved. Concerning the motor load, the acceleration and the angular acceleration of the two followers are confined to $\parallel {\alpha }_k\parallel \le \frac{\text{π}}{2}\text{m}\cdot {\text{s}}^{-2}$ and $\parallel {\beta }_k\parallel \le 30\hspace{0.17em}\text{rad}\cdot {\text{s}}^{-2}$ , respectively.

Figure 4.

Maneuvers of the multi-robot system from a string formation to a triangular one while moving along straight lines

In the process of achieving the formation task, some comparisons between the adaptive integral SMC and the sole integral SMC are shown in Figures 5, 6, and 7, where the uncertainties in the sole integral SMC-based formation control system are given a conservative boundary.

Figure 5.

Performance comparisons between the integral SMC-based formation systems. (a) Distance l₁₂, (b) relative angle ψ₁₂, (c) distance l₁₃, and (d) relative angle ψ₁₃. SMC: sliding mode control.

Figure 6.

Performance comparisons of u₂. (a) Acceleration of follower 2 α₂ by NDO-based integral SMC, (b) angular acceleration of follower 2 β₂ by NDO-based integral SMC, (c) acceleration of follower 2 α₂ by integral SMC, and (d) angular acceleration of follower 2 β₂ by integral SMC. SMC: sliding mode control; NDO: nonlinear disturbance observer.

Figure 7.

Performance comparisons of u₃. (a) Acceleration of follower 3 α₃ by NDO-based integral SMC, (b) angular acceleration of follower 3 β₃ by NDO-based integral SMC, (c) acceleration of follower 3 α₃ by integral SMC, and (d) angular acceleration of follower 3 β₃ by integral SMC. SMC: sliding mode control; NDO: nonlinear disturbance observer.

Figure 5 depicts the system states of the two followers, that is, the distance l_ik and the relative angle ψ_ik, where the dot curves are desired, the solid curves denote the performances by the presented control scheme (short for NDO-based integral SMC), and the dash curves delegate the performances by the sole integral SMC methd (short for integral SMC). Both the SMC systems share the same integral sliding surfaces. Although the control system can achieve the formation maneuvers by the sole integral SMC, the control system is more sensitive to uncertainties because it has no ability to adaptively resist unknown uncertainties. In contrast, the NDO-based integral SMC-based formation system is more robust and its curves are smoother in the dynamic process.

Figure 6 shows the control inputs of the follower 2 by the two integral SMC-based methods. In Figure 6, the curves by the sole integral SMC dramatically jump back and forth. The phenomenon originates from chattering, a drawback of SMC. The jumping points in Figure 6 are different in each of simulations because Δ_i and Δ_k are random. Compared to the results of the two methods, the curves by the NDO-based integral SMC are much smoother because the NDO-based observer can suppress uncertainties so that the control performances can be improved. Some similar results can be draw from Figure 7.

Figure 8 displays each elements of the real sliding surface vector s _ik . As proven, the real sliding surface vectors are of asymptotic stability. The designed control scheme combines the integral SMC method and the NDO technique. The integral SMC-based controller and the NDO-based observer work together to realize the robust formation control for the group of uncertain multiple robots. Figure 9 illustrates the approximation errors ${\widehat{\mathbf{e}}}_{ik}$ . From Figure 9, the NDO-based observer has noticeable errors at the outset. However, the errors are exponentially convergent to guarantee the formation stability and to improve the formation performance in the presence of uncertainties.

Figure 8.

Sliding-surface vectors and their elements. (a) s_12,1, (b) s_{12, 2}, (c) s_{13, 1}, and (d) s_{13, 2}.

Figure 9.

Approximation errors e _ik . (a) e_12,1, (b) e_{12, 2}, (c) e_{13, 1}, and (d) e_{13, 2}.

Conclusions

This article has investigated the formation control problem for swarms of uncertain robots. By the leader–follower mechanism, individual uncertainties from the robots chain to the formation dynamics. The formation uncertainties bounded by an unknown boundary trouble the formation control problem. To resist the formation uncertainties when forming up the robots, the integral SMC-based control scheme is addressed. The control scheme consists of an integral sliding mode controller and a NDO. The theoretical analysis proves that the formation control system in the presence of uncertainties is of asymptotic stability. The control scheme has achieved the formation maneuvers by a multi-robot platform. The simulation results have demonstrated the effectiveness of the control scheme.

The control scheme benefits us a lot because it can approximate the boundary of uncertainties so that the formation control system can have the guaranteed stability. The effort to consider the communication topology of the multi-robot systems via the control scheme is our continuous interest and this field is still in progress.