Decentralised Formation Control and Stability Analysis for Multi-vehicle Cooperative Manoeuvre

Aolei Yang Wasif NaeemMinrui Fei

I.INTRODUCTION

TSTUDY of multi-agent systems(MAS)is concerned with the development and analysis of control architectures and sophisticated artificial intelligence problem-solving for both single-and multiple-agent systems.A typical topic of research in MAS is the multi-vehicle cooperative formation control problem,which is to control the displacement and/or attitude of each vehicle within some workspace to accomplish a common task.Since a vehicle can be considered as an actual implementation of an agent,the word“vehicle”in this paper will be equivalent to the word“agent”.

It can be argued that a group of vehicles can perform tasks faster and more efficiently compared to a single vehicle.For example,a group of robots can be employed to move a large object which may be difficult or impossible for a single robot;multiple unmanned aerial vehicle(UAV)and autonomous underwater vehicle(AUV)formations achieve better area coverage for reconnaissance,thereby improving mission success[1].In recent years,research in cooperative flight control of multi-aircraft systems has attracted growing interest.The implications of the formation flight are very significant,not just for fuel and energy saving but for future applications across civilian and military domains[2].For instance,formation flight can be used to handle the increase in air traffic around airports by means of civil aeroplane formations taking-off and landing to improve the efficiency of runway use.Applications in the military field include aerial refueling,include aerial refueling,aircraft logistics,air formation systems.

More specifically,there are four cooperative formation control strategies of importance:behaviour-based strategy,leader following approach,virtual structure method and artificial potential field approach.

The behaviour-based approach is decentralised and the control action of each element in the formation is a weighted average of the control for each behaviour[3].Its main advantage is that the collision avoidance problem can be easily dealt with due to the existing reactions between vehicles.However,the whole system is more complex and difficult to be analysed mathematically.

In the leader-following approach[4],some vehicles are considered as leaders,while others act as followers.The leader vehicle maintains the given trajectory while the followers track a fixed relative distance from the designated neighbouring vehicles[5].This technique is decentralised,conceptually simple and widely employed in multi-vehicle formation control.Its major drawback is that the leader becomes a single point of failure,as well as producing error propagation that incurs poor responses for the followers.

The virtual structure approach is centralised and describes the whole formation as a single rigid body[6].Each vehicle has its own relative position in the body and tracks the desired trajectory which is translated from the desired motion of the rigid body[7].Its strength is that all elements in the formation hold the same transient without error propagation,and a tight formation can then be implemented and maintained.However,it contains a single point of failure which is undesirable,since collision avoidance is not possible due to the absence of interaction between the vehicles,as well as a lack of explicit feedback to the formation.

The standard artificial potential field approach is based on applying the negative gradient of an artificial potential function(APF)as control inputs to drive the overall formation to convergence[8].This technique is decentralised,reactive,and can be beneficially combined with the behaviour-based method to make use of both benefits.However,it suffers from a problem with local minima and there is computational difficulty in designing a feasible artificial potential function.

1)Stability analysis approach

Any closed-loop formation control methodology must naturally be guaranteed to be stable,especially where there is a potential of loss of human life and property.Most research papers in this area make use of either algebraic graph theory or Lyapunov stability theory.

The use of graph theory in the analysis of interconnected systems is not new,and the current interest is in applying graph-theoretic concepts to the multi-vehicle formation problem.An application of graph theory was discussed in[9],where a directed graph was used to represent the communication network and to relate its topology with forma tion stability.Desai et al.in[10]presented a framework for describing the behaviours of robots in a formation,representing possible control graphs and the coordination of transitions with formation changes from one geometry to another.Additionally,in[11],elementary graph theory and linear matrix inequalities(LMI)were combined with logic-based switching to shed light on the various control strategies which are feasible in the leader follo wing framework for formation flying of multiple spacecraft.Furthermore,graphs have also been applied to describe the interconnection topology in a formation[12],to re flect the control structure,information flow,and error propagation.

Lyapunov stability theory is an alternative way to analyse the overall formation stability.In[13],a formation Lyapunov stability function was defined as a weighted sum of the control Lyapunov function for each vehicle to support the formation stability analysis.In[8],the idea of relative-position-based formation stability was proposed and the Lyapunov method was also used to design the decentralised controllers,along with an extended LMI to analyse the conditions required for formation stability.Moreover,Hong et al.in[14]proposed a Lyapunov-based approach to give a sufficient condition to make all the agents converge to a common value,and a common Lyapunov function was explicitly constructed in the case of switching jointly connected topologies.In another paper[15],a combination of the above two strategies was employed,where a few concepts from graph theory were borrowed to evaluate the controllability and observability of the individual system and a vector Lyapunov method was then applied to prove the stability of the multi-vehicle formation.

2)Contribution of this paper

The purpose of this paper is to propose a formation stability conjecture which is used to conceive a new methodology for decentralised multi-vehicle formation control and to analyse its associated formation stability.Although differing from the previous strategies,it necessarily builds on them and has the desired characteristics.

Based on the formation stability conjecture,the“extensiondecomposition-aggregation(EDA)”process[16]is adopted to support the decentralised formation controller design and formation stability analysis.A virtual additional system(VAS)is first employed to build the relationship between isolated vehicles and to produce an overall formation control system involving formation information.Decomposition to the overall formation system is then carried out to translate this complex formation system into a group of feasible subsystems interacting via boundary conditions.These decomposed subsystems are referred to as individual augmented subsystems(IASs).A scalar Lyapunov function is next selected as an index for representing the stability of each IAS.These indices are f i nally aggregated to analyse the subsystem interactions and the stability of the overall formation using Lyapunov theory.It is concluded that the overall formation system is not only asymptotically but also exponentially stable in the sense of Lyapunov within a neighbourhood of the desired formation,if all the IASs can be stabilised by using any approach.It is essential to note that the proposed methodology on multivehicle formation control and stability analysis in this paper can be easily extended into research of other topic in multiagent system.

The remainder of the paper is organised as follows.Section II outlines the problem formulation.Section III presents the formation stability conjecture and its specific EDA scheme,whilst Section IV describes the relevant stability analysis of the overall formation system.Extensive simulation results are presented in Section V to illustrate the feasibility and verify the formation stability.The paper finally concludes with some discussion and suggestions for future work.

II.PROBLEM FORMULATION

It is essential to note that only using the formation geometry between the vehicles is not suf ficient to define a formation,since the motion of each vehicle in 3-D space has six DOF defining the position and attitude,and unstable attitude will definitely affect the above relative distance.On the other hand,it is believed that when a formation change is needed,the related attitude should generally be regulated by a formation control algorithm until all the vehicles achieve their new desired formation geometry.In other words,the primary function of formation controller is to maintain the stability of the relative positions by controlling the states of each vehicle(such as aircraft attitude angles).

Motivated by this,a complete formation definition(FD)in Euclidean space should consist of the relative distances,and states of all the vehicles.Here,is to denote the state of theith vehicle,whilst the overall state of the formation system withNvehicles is defined as a vector=by combining the states of all the individual vehicles.Taking these into account,the FD is expressed as(1)in terms of the relative distances and each vehicle′s states.

where the subscript“real”represents the real-time formation,and“rrij”means relative distance between theith vehicle and thejth vehicle.

Similar to the state definitionof the overall formation system,the input vector of the overall formation system is defined aswhich consists of input vectors of all the individual vehicles.Based on the above definitions,a common and intuitive definition to formation stability is now presented below.

Formation stability definition 1.The formation withNvehicles is asymptotically stable in the sense of Lyapunov if the following is satisfied[1],

wheretis the elapsed time,defined by(3)is called cumulative formation error,and denotes the summation over the individual formation errors.Lerepresents the equilibrium value of,and from the viewpoint of an optimisation problem,Leshould be the minimum of,i.e.,.

In(3),is a cost function and represents the individual formation error between theith vehicle andjth vehicle,and the final termis a penalty on the inputs.

It is obvious that the above definition illustrates the formation stability from a global perspective,but may not be convenient for analysing the stability of the decentralised formation control strategy.Here,based on the formation definition in(1),another definition is given as follows.

Formation stability definition 2.The formation withNvehicles is asymptotically stable in the sense of Lyapunov if the followings are satisfied,

wheretis the elapsed time,Feindicates the equilibrium shape of the overall formation,andconsists of the equilibrium states of,i=1,2,···,N.This formation stability definition means that both the states of all the vehicles and the relative distances between them should converge to steady states.

It is obvious thatFrealin(4)shows a global expression of the overall formation.Based on the adopted decentralised control architecture,it is assumed thatFrealcan be decomposed into a set of local-formation representations,i.e.,

whereFrealicorresponds to the local-formation representation of theith vehicle.Suppose that it can be expressed by a vectorto reflect the variation of the associated local-formation.The state vectorof theith vehicle can then be augmented as

whereis considered as the augmented state vector derived from the original one,of theith vehicle.Based on this,formation stability definition 2 is then rede fined as follows.

Formation stability definition 3.The formation withNvehicles is asymptotically stable in the sense of Lyapunov if the following is satisfied,

whereconsists of all the state vectors of the augmented systems derived from the original vehicle systems,anddenotes the equilibrium states including all the individual equilibrium states of the augmented systems.Although the above three definitions show different forms,they are equivalent to each other in essence.

III.NOVEL FORMATION CONTROL STRATEGY

Most formation control approaches mentioned previously are direct methods and are based on direct distance or position errors.Such methods are easy to understand,but are not necessarily suitable for analysing and implementing a practical formation control system.Here,a formation stability conjecture is proposed to conceive a new methodology for dealing with multi-vehicle decentralised formation control problem.

A.Formation Stability Conjecture

As to the problem of formation control of multiple vehicles,it is obvious that ensuring the stability of all the original vehicles in a formation is not suf ficient to prove the stability of the overall formation because equation(4)cannot be guaranteed.However,are there corresponding augmented systems from the original vehicles to support the analysis of formation stability?If the answer is yes and the stability of all the augmented systems is able to guarantee the stability of the overall formation system,the process of solving the formation control problem can be greatly simplified.Hence,motivated by this idea,we propose a conjecture as follows:

Formation stability conjecture.The overall formation system is guaranteed to be stable if the following three requirements are satisfied.

1)There exists a group of dynamic systems(DSs)which are augmented from corresponding original vehicles.

2)All the DSs interacting with each other completely constitute the overall formation control system.

3)Each individual DS can be stabilised by the respective decentralised controller.

The relationship between the conjecture and the analysis of the formation stability is intuitively given in Fig.1,where stage 1 will be discussed in Section IV,the stage 2 is supported by the formation stability definition 3.It is essential to point out that the DS represents a system that needs to be determined or designed.Recalling the definition of the augmented states in(8),a DS can be considered as the augmented system of each original vehicle,and the combination of states of all the DSs can also be denoted by,wherexximeans the augmented state vector of theith DS.If all the DSs can be stabilised by their decentralised controller,it means that the stability requirement in(8)is guaranteed,and consequently the whole formation system is stable.

Based on the proposed conjecture,the analysis process of the overall formation stability can be brie fly listed as the following procedure:

Step 1.Attempt to seek or determine feasible DSs according to a desired speci fication.

Step 2.Design individual controllers to stabilise all the achieved DSs.

Step 3.Carry out simulation analyses to verify the performance of the overall formation system.If the achieved performance is not satisfied with the desired target,return to Step 1.

Step 4.If the performance in Step 3 is satisfactory,analyse the stability of the whole formation mathematically.

In the above descriptions,the determination of the DSs is the strategic point,and the stability of the overall formation is completely dependent on them.In other words,if the suitable DSs are found,mathematical stability analysis and implementation of the decentralised formation controller become feasible and convenient.

Fig.1.Analysis process of the formation stability based on the conjecture.

Based on the above conjecture and the subsequent procedure,the“extension-decomposition-aggregation(EDA)”proposed in[16]is considered as an implementation of the formation stability conjecture.In contrast to[16],this paper focuses on mathematically analyzing the stability of multivehicle formation based on this EDA.This will be discussed in the next subsection.

B.Extension-decomposition-aggregation Strategy for Formation Control

The design framework of the EDA scheme is illustrated in Fig.2.The introduced virtual additional system(VAS)is:1)to build a relationship between the isolated vehicles,leading to a new overall formation control system;2)to involve the desired formation variables or parameters,which could be combined with the related individual vehicle model;3)to support the subsequent decomposition,and simplify the stability analysis to the overall formation.

Note that the VAS is merely an algorithm to act as an“interaction bridge”providing each vehicle with the capability of sensing its local-formation states,which can then be used for the formation controller design.

Since the overall formation system involving many variables is difficult to handle as a whole,it is natural to decompose it into several local-formation subsystems.However,it should be noted that there is no general systematic procedure for decomposing such a complex dynamical control system.Here,using physical insight,the overall formation system is decomposed intoNindividual subsystems,each being called individual augmented subsystem(IAS)since it combines the individual vehicle model and the local-formation variables.

In light of the above definitions and descriptions,the procedure for analysing the whole formation stability is specified as follows:

1)All the decomposed subsystems or IASs are used to design the decentralised formation controllers,the interactions between them being considered as the bounded exogenous disturbances.

2)The formation stability of each vehicle,involving several state variables,is represented by a single variable—the Lyapunov function(say,Viforith vehicle),which is considered as an index for each subsystem representing local-formation stability.This approach simpli fies the stability problem,although it sacrifices the detailed information about the size of variation of each separate state variable.

3)These indices(Vi,i=1,2,···,N)are aggregated and the interactions between them are studied to analyse mathematically the stability of the overall formation,neglecting the interconnections between decomposed subsystems.

Broadly speaking,theNIASs are associated with theNscalar Lyapunov functions(Vi,i=1,2,···,N),and each determines the stability property of an individual IAS.These scalar functions are considered as the components of the overall Lyapunov functionV.A differential inequality is formed in terms of this function,using the original inequalities of the decomposed IASs.The stability of the overall formation control system can then be determined by considering only the differential inequalities involving theNLyapunov functions.This brings about a considerable reduction in the dimensionality of the formation stability problem.

IV .FORMATION STABILITY ANALYSIS

Recalling the EDA scheme shown in Fig.2,the coupled multiple pendulum(CMP)[17]system is suitable to be considered as a candidate of the VAS,where the variation ofεiis a reflection of the formation error,and there also exists an important equation expressed as(9)[18]

whereTijdenotes the torque on theith pendulum exerted by the interactingjth pendulum.The achieved equations above crucially show that if all the individual pendulums are stable,the resultant torque of the whole CMP is equal to zero whatever the interconnection formation topology is.This useful result can be used to support the overall formation stability analysis of the multi-vehicle formation.

It is assumed that the overall formation control system,Sconsisting of the CMP system as the VAS,can be linearised and decomposed intoNcoupled subsystems or so-called IAS,as

whereanddescribe the interactions of theith IAS with its neighbouring(i−1)th and(i+1)th IASs,respectively.Now consider the stability of the corresponding decoupled linear time-invariant subsystemsby ignoring the above interactions.

For the stability ofin(11),a scalar function

is defined as a candidate Lyapunov function.Thus,for any choice of the positive definite matricesQi,there exist positive definite matricesPias the solution to the Lyapunov matrix equation(13).

The time derivative of the functionalong the solution of subsystemis therefore given by

whereFrom the above equations,the range estimates forandare given by[19]

where

and the properties of a positive semi-definite symmetric matrix(H)shown in(17)and(18)are used in achieving(16),

whereλMandλmare the maximum and minimum eigenvalues of the indicated matrix,respectively.Note that the coef ficientsci1,ci2,ci3,andci4are positive because of the positive definite property of matricesPiandQi.

Now the Lyapunov functionVi(xxi)is used as the stability index for each linear subsystem.Furthermore,the derivativealong the solutions of the interacting IAS(SSSi)in(10)is given by

Note that the first term,(grad Vi)TAi xxxi,is associated with the decoupled system(11),and the expression,,follows from the first inequality in(15).So,there exists an inequality given by

The stability of the overall formation system,is expressed by considering the interactions between theseVi(xxi),i=1,2,···,N.The function given in(21)defines the Lyapunov functionforSS,which is the weighted sum of all the individual Lyapunov functions,Vi(xxi),i=1,2,···,N,

where the choice of coefficientski,i=1,2,···,Nreflects that vehicles in a formation may have completely different dynamics.Its time derivativeis given as

By taking the norm on the right-hand side of(22)and using the estimates in(15)and(20),theis then given

wherekmax=max{k1,k2,···,kN},▽max=max??????and.

Moreover,depending on the characteristics of the CMP,the interaction terms,hij,i/=j,i,j=1,2,···,Nin the model of(11)correspond to the torques of the CMP system.The result achieved in(9)can be used to generate the following equation:

If all the decomposed IASs can be stabilised,these terms satisfy the condition in(24).By incorporating(24)in(23),the following differential inequality is obtained.

The above inequality is valid for allt≥t0,and its solution is given by

whereis the initial value offor the overall formation system,andV10,V20,···,VN0denote the initial values corresponding to the subsystems(,i=1,2,···,N).Then,using the first inequality in(15)for the functionsas well as the definition(21),the following inequalities hold,

wherecM=max＞0,andi=1,2,···,N,denote the initial states of all the IASs.As for the whole formation systemSin(10),a solution,···,,is defined to combine the solutions of,i=1,2,···,N.Now the relationship between their norms is given by

Combining(27)and(28)generates two inequalities,

and applying the above inequality to(26)produces the final result for showing the stability of the overall formation.Thus,

whereandXXXis the solution of the whole system,SS,for the initial conditions(t0,XXX0).Based on the formation stability definition 3,the inequality in(30)indicates that the solution of the overall formation system exponentially converges to the equilibrium pointXXXe=C‖XXX0‖,which is associated with the neighbourhood of the desired formation.

Based on this unique and significant result,the EDA strategy will be applied to multi-UAV formation flight and the corresponding decentralised formation controllers are implemented to guarantee the formation stability.Simulation studies were performed to illustrate the feasibility and stability of the overall formation.

V.APPLICATION:MULTI-UAVFORMATION MANOEUVRE

The performance of the proposed methodology was evaluated by using the 6-DOF Aerosonde UAV nonlinear simulink model[20].The autonomous Aerosonde has a length of 1.74m,a wingspan of 2.87m,maximum payload capacity of 5kg and endurance of up to 30hours.It was designed for applications including long-range weather data acquisition and reconnaissance over oceanic and remote areas.

A.Implementation of The Decentralised Formation Controller

As presented in[17],the longitudinal and lateral autopilots framework can be employed to design the corresponding longitudinal and lateral formation controller.Using the LMI-basedH∞control methodology[21],the calculated state-feedback gain matrices of the longitudinal and lateral formation dynamics are given by

Furthermore,the longitudinal and lateral compensators were designed heuristically to eliminate the corresponding formation errors,and their parameters are listed in Table I.

TABLE I PI PARAMETERS OF ALL THE COMPENSATORS

To perform the computer simulations and evaluate the performance of the proposed formation control strategy,the nonlinear Aerosonde model was utilised with all the controller gains and parameters as shown previously.

B.Maintenance and Change of Multi-UAV Formation

To illustrate the performance of the maintenance and change during multi-UAV motion,several general assumptions and definitions are necessary to be given first:

1)The reference airspeed for all the vehicles with respect to the inertial reference frame is 26m/s,and the constraint on the airspeed is 15～50m/s.Also,the defined constraints on the pitch and roll angles are given as±20deg and±45deg,respectively.

2)All the UAVs were initially flying with the desired formation geometry as shown in Fig.3(a),where the formation communication topology is assumed to be bidirectional partialmesh.

3)During the formation changes,the formation topology is maintained.

This topological structure of the multi-UAV formation is intuitively displayed in Fig.3,where the diamond symbol labelled with number indicates the corresponding UAV,whereas UAV 2 acts as the reference vehicle,and,i,j=1,2,3,4,i/=jis to describe the reaction relationship between the UAVs.The double-headed arrow means that there exists the bidirectional reaction or interaction between the related UAVs.The altitude formation or the relative coordinates in theZLaxis are indicated by the neighboured annotations(pLiz,i=1,2,3,4).Alternatively,the desired formation change process can be mathematically expressed by where subscript“d”means the desired formation representation and superscript“p”represents“position coordination”in the local-formation reference frame.

The designed scenario of the formation manoeuvre is illustrated as:1)The group of UAVs is tasked to navigate in 3-D space with the constant reference heading angleψref=0rad;2)Fig.3(a)shows the initial formation arrangement with the defined bidirectional partial-mesh formation topology;3)All UAVs were flying steadily at an altitudeH=800m in the first stage.Whent=20s,the formation started to climb until their altitudes were equal to 900m;4)Figs.3(b)and 3(c)illustrate the 3-D formation changes att=80s and att=120s respectively.

Fig.4.Motional trajectories of the UAVs formation in 3-D.

For the whole simulation,the resulting 3-D motional trajectories of the UAV group are shown in Fig.4.In order to facilitate the observation,the above 3-D flight trajectories are separately projected onto the horizontal and vertical planes.The resulting horizontal projection is shown in Fig.5,and the vertical one is shown in Fig.6 which represents the altitude dynamics.Furthermore,the achieved states and individual formation errors are displayed in Figs.7 and 8,respectively.

Based on these figures,the following observations could be made:

1)As shown in Fig.5,the planar formation geometry satisfies the desired formation shapes:the desired formations defined in Figs.3(b)and 3(c)were obtained at aboutt=120s andt=210s,respectively.As shown in Fig.6,the climbing flight and the altitude formation change were achieved with a satisfying performance.

2)The forward,lateral and altitude formations of the multi-UAVs were stable and followed the reference formation′s commands closely.This is an argument for supporting the result of the overall formation stability achieved in Section IV.

3)Variations in the velocities of the UAVs could be observed during the formation changes.This was necessary to rapidly achieve the new reference formation,for instance,UAV 1 needed to speed up to align itself with UAV 3 whent＞80s.Moreover,all the vehicles in the formation attained the reference speed of 26m/s once the new formation was achieved.

Fig.5.Planar projections of 3-D trajectories of UAVs formation corresponding to Fig.3.

4)As depicted in Fig.7,it is clear that the airspeed and attitude(pitch and roll)of each UAV remained within the specified constraints during the cooperative flight.

5)The individual formation errors are shown in Fig.8.During the formation changes,the transient formation errors asymptotically converged to zero after the new formation was achieved.This means that all steady-state formation errors were effectively eliminated by the supplemented compensators.

VI.CONCLUDING REMARKS

Fig.6.Altitude dynamics of the four UAVs.

Fig.8.Individual formation errors of the four UAVs.

In this paper,based on the given formation definition and the three types of formation stability definitions,the formation stability conjecture is first proposed,and the EDA strategy is then specified to support the decentralised formation controller design.Throughout the subsequent stability analysis,it is proven that with the support of the CMP system,the stabilisation of all the IASs is sufficient to guarantee the stability of the overall formation system.The designed decentralised controllers are used to reject the impact of the formation changes and maintain the stability of each IAS,and subsequently guarantee the overall formation stability.Simulation studies on Aerosonde UAVs have been carried out to verify the achieved formation stability result.Future work includes formation stability analysis under dynamic formation topologies,and application of the methodology presented in this paper to other multi-agent problem-solving.

ACKNOWLEDGEMENT

The authors would like to acknowledge the EPSRC under UK-China Science Bridge grant(EP/G042594/1)and Flight Test Center of Commercial Aircraft Corporation of China(COMAC).

REFERENCES

[1]Murray R M.Recent research in cooperative control of multi-vehicle systems.Journal of Dynamic Systems,Measurement,and Control,2007,129(5):571−583

[2]Turpin M,Michael N,Kumar V.Trajectory design and control for aggressive formation flight with quad-rotors.Autonomous Robots,2012,33(1−2):143−156

[3]Monteiro S,Bicho E.Attractor dynamics approach to formation control:theory and application.Autonomous Robots,2010,29(3−4):331−355

[4]Wang J L,Wu H N.Leader-following formation control of multi-agent systems under fixed and switching topologies.International Journal of Control,2012,85(6):695−705

[5]Chen G,Lewis F L.Leader-following control for multiple inertial agents.International Journal of Robust and Nonlinear Control,2011,21(8):925−942

[6]Mehrjerdi H,Ghommam J,Saad M.Nonlinear coordination control for a group of mobile robots using a virtual structure.Mechatronics,2011,21(7):1147−1155

[7]Askari A,Mortazavi M,Talebi H.UAV formation control via the virtual structure approach.Journal of Aerospace Engineering,doi:10.1061/(ASCE)AS.1943-5525.0000351

[8]Xue D,Yao J,Chen G,Yu Y L.Formation control of networked multi-agent system.IET Control Theory and Application,2010,4(10):2168−2176

[9]Fax J A,Murray R M.Information flow and cooperative control of vehicle formations.IEEE Transactions on Automatic Control,2004,49(9):1465−1476

[10]Desai J P,Ostrowski J P,Kumar V.Modeling and control of formations of nonholonomic mobile robots.IEEE Transactions on Robotics and Automation,2001,17(6):905−908

[11]Mesbahi M,Hadaegh F Y.Formation flying control of multiple spacecraft via graphs,matrix inequalities,and switching.Journal of Guidance,Control,and Dynamics,2001,24(2):369−377

[12]Olfati-Saber R,Murray R M.Distributed structural stabilization and tracking for formations of dynamic multiagents.In:Proceedings of the 41st IEEE Conference on Decision and Control.Las Vegas,USA:IEEE,2002.209−215

[13]Ogren P,Egerstedt M,Hu X M.A control Lyapunov function approach to multiagent coordination.IEEE Transactions on Robotics and Automation,2002,18(5):847−851

[14]Hong Y G,Gao L X,Cheng D Z,Hu J P.Lyapunov-based approach to multiagent systems with switching jointly connected interconnection.IEEE Transactions on Automatic Control,2007,52(5):943−948

[15]Feddema J T,Lewis C,Schoenwald D A.Decentralized control of cooperative robotic vehicles:theory and application.IEEE Transactions on Robotics and Automation,2002,18(5):852−864

[16]Yang A L,Naeem W,Irwin G W,Li K.Novel decentralised formation control for unmanned vehicles.In:Proceedings of the 2012 IEEE Intelligent Vehicles Symposium(IV 2012).Alcala´ deHenares,Spain:IEEE,2012.13−18

[17]Yang A L,Naeem W,Irwin G W,Li K.A decentralised control strategy for formation flight of unmanned aerial vehicles.In:Proceedings of the 2012 UKACC International Conference on Control.Cardiff,UK:IEEE,2012.345−350

[18]Yang A L,Naeem W,Irwin G W,Li K.Stability analysis and implementation of a decentralised formation control strategy for unmanned vehicles.IEEE Transactions on Control Systems Technology,doi:10.1109/TCST.2013.2259168

[19]Khalil H K.Nonlinear Systems(3rd edition).New Jersey:Prentice Hall,2001

[20]Marius Niculescu.AeroSim Aeronautical Simulation Block Set V1.2,Users Guide.Oregon,US:Unmanned Dynamics,LLC,2003

[21]iljak D D,StipanoviD M.Robust stabilization of nonlinear systems:the LMI approach.Mathematical Problems in Engineering,2000,6(5):461−493