Open access peer-reviewed chapter

Robust Adaptive Cooperative Control for Formation-Tracking Problem in a Network of Non-Affine Nonlinear Agents

By Muhammad Nasiruddin bin Mahyuddin and Ali Safaei

Submitted: December 17th 2016Reviewed: April 20th 2017Published: September 13th 2017

DOI: 10.5772/intechopen.69352

Downloaded: 792

Abstract

In this chapter, a decentralized cooperative control protocol is proposed with application to any network of agents with non-affine nonlinear multi-input-multi-output (MIMO) dynamics. Here, the main purpose of cooperative control protocol is to track a time-variant reference trajectory while maintaining a desired formation. The reference trajectory is defined to a leader, which has at least one information connection with one of the agents in the network. The design procedure includes a robust adaptive law for estimating the unknown nonlinear terms of each agent’s dynamics in a model-free format, that is, without the use of any regressors. Moreover, an observer is designed to have an approximation on the values of control parameters for the leader at the agents without connection to the leader. The entire design procedure is analysed successfully for the stability using Lyapunov stability theorem. Finally, the simulation results for the application of the proposed method on a network of nonholonomic wheeled mobile robots (WMR) are presented. Desirable leader-following tracking and geometric formation control performance have been successfully demonstrated through simulated group of wheeled mobile robots.

Keywords

  • cooperative protocol
  • formation control
  • decentralized control
  • robust adaptive law
  • distributed observer
  • mobile robot
  • non-affine nonlinear system

1. Introduction

Great attention has been paid to the problems of the multi-agent network ranging from consensus, collective behaviours of flocks and swarms, formation control of multi-robot systems, leader-following, algebraic connectivity of complex network, rendezvous, containment and so on [16]. The formation control problem is an interesting issue in biology, automatic control, robotics, artificial intelligence and so on, which requires each agent to move according to the prescribed trajectory. Various control strategies have been formulated to achieve the group control objectives.

The systems are usually in nonlinear form due to unpredictable environmental disturbances, unmodelled dynamics or other uncertainties. A class of nonlinear first-order multi-agent systems with external disturbances consensus problem was discussed in Ref. [7], whereas other works that involve second-order and higher order nonlinear multi-agent systems are reported in Refs. [8] and [9], respectively. Wang et al. [10] reported the design of distributed state/output feedback cooperative control approaches for uncertain multi-agents in undirected communication graphs. This is later extended to a condition of directed graphs containing a spanning tree [11]. To remedy the problem of a non-affine system for a general class, several reported works such as Ref. [12] employ a direct adaptive approach using an artificial neural network (ANN) to approximate an ideal controller. By employing a system transformation, a non-affine system can be transformed into an affine system as demonstrated in Ref. [11]. However, the transformation technique to convert a multi-agent non-affine system to a multi-agent affine system is still new and open to further studies which are to be discussed in this chapter.

Hou et al. [13] illustrate the method of dealing with non-affine multi-agent system by incorporating dynamic surface control or DSC but it is limited to a single-input-single-output (SISO) type of system, that is, with one control input. A similar approach is reported in Ref. [14] where the distributed dynamic surface design approach is used to design local consensus controllers using the transformation to convert the system to an affine strict-feedback multi-agent system. The work is also limited to a single control input per agent.

In this chapter, several novel contributions can be highlighted, that is, the introduction of transformation techniques from a non-affine multi-agent system to an affine multi-agent system for a network of generic nonlinear multi-input-multi-output (MIMO) systems, that is, a single agent may have more than one control input and more than one output. The second contribution to be highlighted in the chapter is the estimation of nonlinear terms in the dynamics without requiring the linear-in-parameter condition (LIP), that is, the dependence on any model regressor is elevated. The lumped nonlinear function existing in the model agent can be estimated online despite time-varying characteristics. This implies that the estimation is model free. By virtue of a sigma-modified adaptive law with projection algorithm that drives the estimation using the cooperative consensus error, the unknown nonlinear function can be reconstructed. The proposed cooperative control scheme requires a robust adaptive observer which can reconstruct the control signal from all agents to be used in the consensus formation control. Owing to the robustification term in the observer, the control signals can be estimated in finite time. The proposed robust adaptive formation control is to be exemplified in a form of simulation of multi nonholonomic mobile robots with differential drive configurations. They are commissioned to follow the leader trajectory while at the same time required to maintain predefined geometric formation guaranteeing safe inter-agent separation.

The chapter is organized into preliminaries, problem definition, design procedure of the proposed robust adaptive formation control algorithm, simulated results and lastly the conclusion of the chapter.

2. Preliminaries

2.1. Mean value theorem

Suppose that the function Fis continuous on the closed interval [a,b]and differentiable on the open interval (a,b)(i.e. Fis Lipschitz). Then, there is a point X0in the open interval (a,b)at which [15]

F˙(X0)= F(b)F(a)baE1

In physical terms, the mean value theorem says that the average velocity of a moving object during an interval of time is equal to the instantaneous velocity at some moment in the interval [15].

2.2. Kronecker product

The Kronecker product of matrices ARm×nand BRp×qis defined as [16]

AB=[a11Ba1nBam1BamnB] E2

which satisfies the following properties [16]

(AB)(CD)=(AC)(BD)
(AB)T=ATBTE3
A(B+C)=AB+AC

2.3. Schur complement lemma

For any constant symmetric matrix S=[S11S12S12TS22], the following statements are equivalent [17]

   S>0   S11>0 . S22S12TS111S12>0   S22>0 . S11S12S221S12T>0E4

2.4. Graph theory preliminaries

Consider a network consisting of N agents. Let G(V,E,A)be a graph with the set of N nodes V={ν1,ν2,, νN}, a set of edges E={eij}RN×Nand associated adjacency matrix A=(aij)RN×N. An edge eij in Gis a link between a pair of nodes (νj,νi), representing the flow of information from νj(as parent) to νi(as child). The eij is in existence if and only if aij>0. The graph is undirected, that is, the eij and eij in Gare considered to be the same. We name νiand νjas neighbors if eijE . A path is defined as a sequence of connected edges in a graph. A graph is connected if there is a path between every pair of the nodes. The degree matrix DL=diag{d1,d2,… ,dN}RN×N, where each di is the input degree to each node, which is equal to the number of all edges through it (i.e. di=j=1:Naij). Hence, we can define Laplacian Matrix (L) as below [16, 18, 19]

L=DLAE5

Furthermore, we can define an adjacency matrix for the leader as follows

B=diag{b1,b2,… , bN}RN×NE6

where each bi indicates the existence of a communication link between the leader and each agent [16, 18, 19]. Besides, we would have,

H=L+BE7

3. Problem definition

Consider a network of N agents with general non-affine nonlinear dynamics for each of them. The problem is to design a set of decentralized control protocols for all agents to enhance a desired formation in the state space and also track a reference trajectory on state variables. Here, a virtual node is considered as the leader, which knows the desired trajectory and has at least one communication link with the agents in the network. It means that some agents are unaware about the leader states and also their control inputs. The whole problem in a general format can be considered as a platform for any possible state space in diverse applications.

For a MIMO system, one can define the following general nonlinear formulation

x˙i1=h1(xi)+R1(xi)+f1(xi,ui)x˙i2=h2(xi)+R2(xi)+f2(xi,ui)x˙in=hn(xi)+Rt(xi)+fn(xi,ui)E8

where n is the number of states for the system, t is the total number of nonlinear terms in the system (which tn), xiRnis the states vector, uiRmis the input (or control parameters) vector, m is the number of control parameters, hj for j=[1,n]is any linear combination on xi, Rj for j=[1,n]is any Lipschitz continuous nonlinear function on xi and fj for j=[1,n]is any Lipschitz continuous nonlinear function on both xi and ui. The last term defines the non-affine property of the system which represents the completely coupled inter-relation between states and control parameters. Each agent dynamic can be represented in matrix form as follows

X˙i=CXi+Ri+FiXi=[xi1,xi2,… ,xin]T C :constant matrixRi=[R1(xi),R2(xi),… ,Rt(xi)]T , tnFi=[F1(xi,ui),F2(xi,ui),… ,Fn(xi,ui)]TE9

where CRn×nis a constant matrix including the multipliers for each state. The elements of C define the dependence of each state’s derivative to the other states.

For a network of N of similar agents (or systems), dynamics for each agent i can be represented by Eq. (9). Also, the dynamic of the leader node can be proposed by this format. The difference is that the control parameters for the leader are defined with respect to a time-varying reference trajectory, that is

x˙01=h1(x0)+h1(u0)x˙02=h2(x0)+h2(u0)x˙0n=hn(x0)+hn(u0)E10

where hjfor j=[1,n]is any linear combination on the leader control parameters (i.e. reference trajectory u0). Actually, the reference trajectory is a set of inputs which provide certain dynamics in state space for the leader agent. The leader dynamics can be represented in the matrix form as the following:

X0˙=CX0+Du0E11
X0=[x01,x02,,x0n]T , u0=[u01,u02,,u0m]TC & D :constant matrices

where DRn×mis a constant matrix including the multipliers for each control parameters.

Moreover, the desired formation among the agents in a network can be presented by a set of constant values F(RN×Rn), which determines the relative distance between agents in the state space.

The problem is to enhance Famong the network agents and track the reference trajectory defined by (x0, u0) at the leader node with inter-agent communication topology defined by the communication graph.

4. Design procedure for robust adaptive cooperative control protocol

This section is dedicated to presenting the design process for cooperative control protocol, an observer to estimate the control parameters of the leader at each agent and a robust adaptive law to estimate the nonlinear terms at each agent. The design process is initiated by dealing with the non-affinity property of the agents.

4.1. Dealing with non-affinity property

Using the mean-value theorem presented in Section 1, for the nonlinear functions fj, which has a coupled terms of xi and ui, we have [19]

fj(xi,ui)u|u=u*=μ=fj(xi,ui)fj(xi,ui¯)uiu¯i , u¯i<u*<uE12

and without any loss of generality we can consider μ = 1 and u¯iis any constant value.

fj(xi,ui)=ui+qj(xi)qj(xi)=fj(xi,u¯i)μu¯i E13

where qj(xi)is an unknown nonlinear function depending only on xi. As can be seen, the non-affine nonlinear function fj(xi,ui)is converted to an affine form. Now, the dynamics of each agent can be modified as

x˙i1=h1(xi)+R1(xi)+h1(ui)+q1(xi)x˙i2=h2(xi)+R2(xi)+h2(ui)+q2(xi)x˙in=hn(xi)+Rt(xi)+hN(ui)+qt(xi)E14

Considering

gj(xi)= Rj(xi)+qj(xi) , j[1,t] , tnE15

where gj(xi)is an unknown nonlinear function depending on xi, the matrix format for each agent dynamics can be presented as

X˙i=CXi+Dui+D1GiD & D1 :constant matricesGi=[g1(xi),g2(xi),,gt(xi)]TE16

where DRn×mis a constant matrix including the multipliers for each control parameter. Actually, the elements of D define the dependence of each state’s derivative to each control parameters. Moreover, D1Rn×tis a diagonal matrix defining the existence of nonlinear functions in the equation for derivative of each state. Elements of D1 can only be one or zero. It should be noted that since tn, we may have some states’ derivatives which do not include any nonlinear terms.

In the following subsections, the elements of Gi, which define the unknown nonlinear functions on each state’s derivative, would be estimated (adapted) online using consensus error of the network.

4.2. Cooperative protocol for formation and tracking problem

For a network of N agents with the dynamics described by Eq. (16), we can have a lumped formulation for the dynamics of all agents using the Kronecker product,

X˙=(INC)X+(IND)U+(IND1)GX=XNn×1=[X1,X2,, XN]T , U=UNm×1=[u1,u2,, uN]TG=GNt×1=[G1,G2,, GN]T , IN=diag{1,1,,1}RN×NE17

For this network, we can define the combined formation and tracking errors in a single formulation in relation to the neighbouring information available to each agent i via the communication graph [16]

ei=j=1Naij((XiXj)(ΔiΔj))+bi((XiX0)(ΔiΔ0))E18

where ΔRn×1is the vector of desired values for states of agents and also the leader. We can consider ei as the consensus error for agent i. Hence

ei=j=1Naij((XiΔi)(XjΔj))+bi((XiΔi)(X0Δ0))E19

By changing the variables, we have

ei=j=1Naij(ZiZj)+bi(ZiZ0)Zi=XiΔiZj=XjΔjZ0=X0Δ0E20

Trying to lump the consensus errors of all agents in an N-array format, we have

E=(HIn)Z(BZ0)1¯Z=ZNn×1=[Z1,Z2,, ZN]TIn=diag{1,1,,1}Rn×n ,1¯=[1,1, ,1]T  RN×1E21

Besides, considering Eq. (17), we can have an N-array form for dynamics of agents in the changed variables space

Z˙=(INC)Z+(IND)U+(IND1)GE22

If the consensus errors of all agents converge to zero, then both formation and tracking objectives are reached, that is

limtE= 0¯E23

Here, the cooperative protocol U is designed using the Lyapunov stability theorem to ensure Eq. (23) is reached. Consider the following Lyapunov function

V=12ETEE24

Then,

V˙=ET((HIn)Z˙(BZ˙0)1¯)V˙=ET((HIn)(INC)Z+(HIn)(IND)U+(HIn)(IND1)G(BZ˙0)1¯)E25

Considering Eq. (3), we have

(HIn)(IND)=(HD)(HIn)(IND1)=(HD1)E26

Besides, using Eqs. (3) and (21), we have

(HIn)(INC)Z=(INC)E+ (BCZ0)1¯E27

Then, Eq. (25) leads to,

V˙=ET((INC)E+(BCZ0)1¯+(HD)U+(HD1)G(BZ˙0)1¯)E28

Forcing V˙<0and referring to Eq. (11), we have

(INC)E+(BDu0)1¯+(HD)U+(HD1)G= PEP=PT>0 , PRNn×NnE29

Hence,

(HD)U= (P+(INC))E(BDu0)1¯(HD1)G E30

Based on Lyapunov stability theorem, using URNm×1in Eq. (30) as the cooperative control protocol will ensure that V˙<0and that E reaches zero asymptotically. Hence, the objectives in formation problem and tracking problem have been accomplished. Expressing the control signal at agent level for agent i

j=1NHijDuj=(Pi+C)eibiDu0j=1NHijD1GjE31
Pi=P(k*,r*) , k*,r*=[((i1)×n+1):(i×n)] , Hij=H(i,j)

and then

HiiDui=(Pi+C)eibiDu0j=1NHijD1Gj j=1iNHijDuj E32

Finally, the control parameter for agent i can be presented as the following

ui=1Hii(DTD)1DT((Pi+C)eibiDu0j=1NHijD1Gj j=1iNHijDuj )E33

Here, a pseudo-inverse method is employed on D.

There are two required conditions on achieving this goal, which are explained in the following assumptions.

Assumption 1. The communication graph should be undirected and connected. It means sufficient information can be available on agents.

Assumption 2. The dynamics of each agent should be completely controllable, that is D matrix should be full rank. It leads us to a state transformation in some applications.

Looking at the proposed cooperative control protocol in Eq. (33), there are two terms, which are not totally available to all agents:

  1. uj (fourth term in the prentices in Eq. (33)), which is the control parameter for the neighbouring agent at the current moment.

  2. Gj (third term in the prentices in Eq. (33)), which includes the unknown nonlinear terms for dynamics of neighbouring agents.

By reaching consensus on the states of agents, we can conclude that the control parameters of each agent has converged to the values of leader control parameters [20]

limt(uju0)= 0¯ , j[1,N]E34

Hence, the control parameters for the neighbouring agent (uj) are approximated by the control parameter of the leader, which in turn will be observed locally at each agent. It means that each agent has its own estimation on u0 and sends it to the neighbouring agents as its control parameter. The observed data will be transmitted to the neighbouring agents via communication graph to compute the control protocols.

The unknown nonlinear terms (Gj) also will be estimated using the consensus error of each agent. Similarly, the adapted data are shared with neighbouring agents through the communication graph.

4.3. Observer design for leader control parameters

Here, the objective is to have consensus on the value of u0 among the all agents in the network. For this objective, we can define the following consensus error for each agent

Δci= j=1Naij(T^iT^j)+bi(T^iu0) E35

where T^iRm×1is the observed vector at agent i for the leader control parameter, and again the aij and bi are the elements of adjacency matrix for the communication graph in the network. Eq. (35) can be represented in a lumped format as the following

Δc=(HIm)T^(Bu0)1¯Δc=ΔcNm×1=[Δc1,Δc2,, ΔcN]TT^=T^Nm×1=[T^1,T^2,, T^N]TE36

If the equation

limtΔc= 0¯E37

is satisfied, we can say that the observation objective is achieved. Considering the following Lyapunov function, we have

V1=12ΔcTΔcE38

Then,

V1˙=ΔcT((HIm)T^˙(Bu˙0)1¯)E39

Since the summation of all elements in each row of the Laplacian matrix is zero, we can say that

(Lu˙0)1¯=0E40

and recalling Eq. (7), Eq. (39) can be written as following,

V1˙=ΔcT(HIm)T^˙ΔcT(Hu˙0)1¯E41

Considering T^˙=Δc+T^', we have

V1˙=ΔcT(HIm)Δc+ΔcT(HIm)T^'ΔcT(Hu˙0)1¯E42

where since (HIm)is the positive definite recalling the Schur Complement Lemma, the first term is surely negative. To achieve V1˙<0, we should show that

V˙11= ΔcT(HIm)T^'ΔcT(Hu˙0)1¯0.E43

Recalling Eq. (3), we have

(Hu0˙)=(HIm)(INu˙0)E44

Hence, the Eq. (43) is,

V˙11=ΔcT(HIm)T^'ΔcT(HIm)(INu˙0)1¯V˙11 ΔcT(HIm)T^'+||ΔcT(HIm)|| (INU˙0M)1¯E45

where U˙0Mis the upper band or maximum absolute value for u˙0. This value should be available beforehand. Now, we should only show that

ΔcT(HIm)T^'+||ΔcT(HIm)|| (INU˙0M)1¯=0E46

Hence,

ΔcT(HIm)T^'=||ΔcT(HIm)|| (INU˙0M)1¯ΔcT(HIm)T^'=ΔcT(HIm) sign(ΔcT(HIm))(INU˙0M)1¯E47

where sign(ΔcT(HIm)) RNm×Nmis a diagonal matrix whose diagonal elements are the signs of each element in ΔcT(HIm)R1×Nm. Finally, since we have

(ΔcΔcT(HIm))1ΔcΔcT(HIm)=INImE48

the second term in T^˙=Δc+T^', is

T^'= sign(ΔcT(HIm))(INU˙0M)1¯E49

and recalling Eq. (36), the rate for the observed parameter is

T^˙=(HIm)T^+(Bu0)1¯  sign(ΔcT(HIm))(INU˙0M)1¯.E50

By using T^˙from Eq. (50), we can have V1˙0, which in turn shows that the consensus error on observation (i.e. Δc) is stable in accordance to the Lyapunov stability theorem. It is obvious that the observed values for u0(i.e. T^) at each agent are computed iteratively using the rate value proposed in Eq. (50).

The lumped format for rate of observer parameter in Eq. (50) can be presented for each agent as the following

T^˙i=Δci(r=1msign(yir)×u˙0Mr)yi=j=1NHijΔcj=[yi1,yi2,,yim] , U˙0M=[u˙0M1,u˙0M2,,u˙0Mm]E51

where Δciis defined as in Eq. (35).

4.4. Adaptive law design for unknown nonlinear terms in each agent dynamics

In this subsection, the objective is to estimate the values of unknown nonlinear terms in each agent dynamics (i.e. G in Eq. (30)). Since, there is not any data available on exact values of G, the estimation error for adaptation process is not available. Hence, the adaptation should be handled using the output error which in this problem is the consensus error (i.e. E in Eq. (21)).

Considering the consensus error in Eq. (21) and the agent dynamics according to Eq. (22), the derivative for consensus error is

E˙=(INC)E+(BDu0)1¯+(HD)U+(HD1)G E52

where G here is the exact value for nonlinear terms. If we put the designed cooperative control protocol (from Eq. (30))

(HD)U= (P+(INC))E(BDu0)1¯(HD1)G^ E53

with G^is the adapted value for the unknown nonlinear terms, into Eq. (52), we have

E˙=PE+(HD1)G˜ , G˜=GG^ E54

Using the following positive definite Lyapunov function

V2=12ETE+12G˜TΓ1G˜E55

where ΓRNn×Nnis a positive definite matrix, we have

V2˙=ETE˙+G˜TΓ1G˜˙V2˙=ETPE+ET(HD1)G˜+G˜TΓ1G˜˙E56

where the first term in the last equation is the negative definite. To show V2˙<0, we have

ET(HD1)G˜+G˜TΓ1G˜˙=0E57

Then,

G˜TΓ1G˜˙=ET(HD1)G˜E58

which in turn leads to this adaptive law

G^˙=G˜˙=+Γ(HTD1T)EΓ=diag{γ1,γ2,,γN} , γi=diag{γi1,γi2,,γit} , tnE59

Considering the Lyapunov stability theorem for the function in Eq. (55), if G^is updated using the rate value proposed in Eq. (59) iteratively, G˜converges to zeros asymptotically. It means that the adapted parameter G^will converge to the actual value of the nonlinear terms in agent dynamics. One of the important issues of the proposed adaptive law in Eq. (59) is that it is not required to include any set of nonlinear basis functions as regressors in the adaptive law. It is only based on the consensus error of the network, which may have sufficient information to tune the adaptive parameter.

Since the adapted signals are always vulnerable for being distracted and diverged by unknown terms, two robusting methods are provided to make the designed adaptive law robust against the divergence [21].

  1. Parameter projection method

    G^˙={Γ(HTD1T)E ,    if G^TG^< M0TM0¯(IΓGGTGTΓG)Γ(HTD1T)E ,  otherwiseE60
    M0¯=[M01¯,M02¯,,M0N¯]T , M0i¯=[M01,M02,,M0t]T , tn

where M0iis chosen so that M0i|gi|. The value for M0 should be defined beforehand. The algorithm is named as parameter projection in the literature [21].

  • σ-modification or leakage method;

    G^˙=+Γ((HTD1T)EρG^) , ρ>0RE61

  • Hence, the complete robust adaptive control for estimating the nonlinear terms in each agent’s dynamics is presented as the following

    G^˙={Γ(HTD1T)EρΓG^ , if G^TG^< M0TM0¯(IΓGGTGTΓG)(Γ(HTD1T)EρΓG^) ,   otherwiseE62
    M0¯=[M01¯,M02¯,,M0N¯] , M0i¯=[M01,M02,,M0t] , tn 

    The lumped format for the rate of adaptive parameter in Eq. (60) can be presented for agent i as the following

    G^˙i={γi(j=1NQijejρG^i) , if G^iTG^i<M0TM0¯(InγiGiGiTGiTγiGi)γi(j=1NQijejρG^i), otherwiseE63
    Q=(HTD1T) , MRNt×NnQij=Q(k*,r*) , k*=[((i1)×t+1):(i×t)] , r*=[((j1)×n+1):(j×n)]

    5. Application: wheeled mobile robot

    In this section, application of the proposed cooperative control protocol on a team including three nonholonomic wheeled mobile robots (WMRs) is presented. The robots are moving on a smooth planar surface with a constraint on the speed (Figure 1). They can only move in the direction of their attitudes and speed in the perpendicular direction is zero. This is a nonholonomic constraint. Few number of researches can be found in literatures, which deal with the cooperative control of the multi-agent of WMRs taking account of each agent’s WMR dynamics [22, 23].

    Figure 1.

    A diagram for kinematics of a nonholonomic planar wheeled robot.

    5.1. Problem definition

    Here, the kinematics and dynamics for motion of ith WMR are considered as the following

    x˙i=υicosθi , y˙i=υisinθi , θ˙i=ωiυ˙i=1mFi , ω˙i=1JTiE64

    where xi and yi represent the position of a single WMR in the inertial coordinate system, θi is the orientation of the WMR, υiis the translational speed in the WMR’s pose direction and ωi is the angular speed of WMR about the Z axis. Also, m and J are the mass and moment of inertia for WMR. Moreover, Fi and Ti are the force and torque generated by the electric motors disclosed in each wheel of WMR. The last parameters are the control parameters for motion of each WMR. By transforming the kinematics of WMR to a local coordinate system fixed to the WMR, [24]

    [xi1xi2xi3]=[cosθisinθi0sinθicosθi0001][xiyiθi]E65

    Then by considering xi4=υiand xi5=ωi,we have

    x˙i1=xi4+xi5xi2 , x˙i2=xi5xi1x˙i3=xi5 , x˙i4=ui1 , x˙i5=ui2E66

    where ui1=1mFiand ui2=1JTi. The state-space system can be represented in matrix form similar to Eq. (16), as the following

    X˙i=CXi+Dui+D1GiXi=[xi1,xi2,xi3,xi4,xi5]T, ui=[ui1,ui2]T, Gi=[xi5xi2,xi5xi1]TC=[0001000000000000001000000] , D=[0000010001] , D1=[1001000000] E67

    As can be seen, D is not full rank. According to assumption 2, we need a change of variables to have D in the full-rank form. Recalling the idea of the back-stepping method [25] we have

    δi1=υisi1 , δi2=ωisi2E68

    Applying the back-stepping method

    si3=ui1s˙i1 , si4=ui2s˙i2E69

    we have

    x˙i1=δi1+δi2xi2+si1+xi2si2x˙i2=δi2xi1xi1si2x˙i3=δi2+si2 , δ˙i1=si3 , δ˙i2=si4E70

    Then, the state-space representation of a single WMR can be represented in following format

    X¯˙i=C¯ X¯i+D¯ u¯i+D1¯ Gi¯X¯i=[xi1,xi2,xi3,δi1,δi2]T , u¯i=[si1,si2,si3,si4]TG¯i=[(δi2xi2+qi1(xi2)) ,(δi2xi1+qi2(xi1))]TC¯=C , D¯=[11000100000100010001] , D¯1=D1 E71

    which has a full rank D¯matrix. Hence, assumption 2 is satisfied and the proposed cooperative controller can be implemented. Hence, we have five state variables, four control parameters and two nonlinear terms for each WMR. At each agent within the network, the nonlinear terms will be adapted using Eq. (63) and the control parameters of the leader will be observed using Eq. (51).

    Here, the desired formation is a rectangle with four agents and four equal edges. The length of each edge is equal and is r. The virtual leader is positioned at the centroid of the geometry (Figure 2). Moreover, the communication graph for this network is shown in Figure 2. The leader information is only available to agent 1. Hence, the adjacency matrices are defined as the following

    Figure 2.

    (Left) A diagram for the desired positions of four agents in a network; (right) the communication graph for a network of four agents and a leader.

    A=[0100101000100110] , DL=[1000020000002001] , B=[1000000000000000]E72

    There is a well-known reference trajectory for this problem in the literature [20], which is presented as the following,

    x0=υrωrsinθ0 , y0=υrωrcosθ0 , θ0=ωrtE73

    where υrand ωr can be any known time-varying functions. Usually, these functions are considered as constant values. In Eq. (73), t is time.

    5.2. Simulation results

    The simulation for the problem defined in Section 5.1 is performed by MATLAB/Simulink. The constant values for running the simulation are presented in Table 1.

    ParameterValue
    Mass of each agent (M)1 kg
    Inertia of each agent (J)1 kg/m2
    Relative position of agents in the network (r)4 m
    Reference velocity (υr)5 m/s2
    Reference angular velocity (ωr)0.25 rad/s
    The adaptation rates (γ1, γ2)0.01 & 0.1
    The leakage factor (ρ)100
    The maximum value for rate of u0 (U˙0M)ones (4,1)
    The maximum value for adapted signal (M0¯)10 × ones (2, 1)

    Table 1.

    The constant parameters for simulation of a network of WMRs.

    Moreover, the values of Pi as the gain values for cooperative control protocol at each agent (see Eq. (33)) are as follows

    P1=diag{10,10,100,10,10}, P2=diag{10,10,12,10,10}P3=diag{10,10,30,10,10}, P4=diag{10,10,55,10,10}E74

    The values in Pi are determined in a way to ensure that the whole matrix P is positive definite and the sufficient transient performance of the whole network is achieved.

    The simulation results for this problem are presented in the following figures. The position of all agents in the X-Y plane is shown in Figure 3. The consensus on both reference trajectory and the desired formation can be seen. Actually, the desired formation is achieved gradually. In addition, the position of the centroid of all agents is compared with the reference trajectory in Figure 4. Moreover, the signals for translational and angular speeds of agent 4 are presented in Figure 5. Finally, the observed data for control parameters of the leader and also the adapted nonlinear terms at agent 4 are shown in Figures 6 and 7. Appropriate performance of proposed algorithms can be inferred by these figures.

    Figure 3.

    The reference trajectory (red) and position of agents in the desired formation (agent #1: blue, agent #2: green, agent #3: black and agent #4: yellow).

    Figure 4.

    The reference trajectory and position of the centroid of the agents in the desired formation.

    Figure 5.

    Translational and angular speed of agent #4.

    Figure 6.

    Observed data for control parameters of the leader at agent #4.

    Figure 7.

    Adapted nonlinear terms at agent #4.

    6. Conclusion

    This chapter is dedicated to the design procedure of a cooperative control protocol for any network consisting of agents with non-affine nonlinear dynamics and multi-input multi-output structure. The main goal is to satisfy a tracking problem for the whole network while maintaining a predefined formation topology in the state space of the agents’ dynamics. The proposed design procedure is including an adaptive law incorporated with a robustification method to estimate the unknown nonlinear terms in the agents’ dynamics. In addition, an observer is designed using the consensus-type error for estimating the leader’s control parameters at each agent. Since there are no complete information links between the leader and all agents, the observed control parameters of the leader are required at each agent to construct the cooperative control protocol. The entire design procedure is analysed successfully for the stability using Lyapunov stability theorem. The presented simulation results for a team of wheeled mobile robots show the appropriate performance of the proposed method.

    Acknowledgments

    This chapter is supported by the Fundamental Research Grant Scheme (FRGS) under Grant No. (FRGS-203/PELECT/6071290), awarded by Ministry of Education Malaysia. Besides, this chapter is under a USM-TWAS Postgraduate Fellowship.

    How to cite and reference

    Link to this chapter Copy to clipboard

    Cite this chapter Copy to clipboard

    Muhammad Nasiruddin bin Mahyuddin and Ali Safaei (September 13th 2017). Robust Adaptive Cooperative Control for Formation-Tracking Problem in a Network of Non-Affine Nonlinear Agents, Multi-agent Systems, Jorge Rocha, IntechOpen, DOI: 10.5772/intechopen.69352. Available from:

    chapter statistics

    792total chapter downloads

    More statistics for editors and authors

    Login to your personal dashboard for more detailed statistics on your publications.

    Access personal reporting

    Related Content

    This Book

    Next chapter

    Management of a Single-User Multi-Robot Teleoperated System for Maintenance in Offshore Plants

    By Seungyeol Lee and Jeon-Il Moon

    Related Book

    First chapter

    Introductory Chapter: Geographic Information Systems and Science

    By Cláudia M. Viana, Patrícia Abrantes and Jorge Rocha

    We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

    More About Us