The flocking control of multiagent system is a new type of decentralized control method, which has aroused great attention. The paper includes a detailed research in terms of distance constrained based adaptive flocking control for multiagent system with time delay. Firstly, the program on the adaptive flocking with time delay of multiagent is proposed. Secondly, a kind of adaptive controllers and updating laws are presented. According to the Lyapunov stability theory, it is proved that the distance between agents can be larger than a constant during the motion evolution. What is more, velocities of each agent come to the same asymptotically. Finally, the analytical results can be verified by a numerical example.
1. Introduction
In recent years, the research on the flocking behavior of multiagent system has attracted great attention. For a series of agents which can apply some simple rules and limited information of neighbors to organize into a coordinated state is called flocking phenomenon. There exist many forms of flocking behavior in nature, for example, flocking of birds, swarming of bacteria, and so on [1, 2]. With the development of technology and the importance of real application, the study of the flocking behaviors of multiagent system has caused attention from a lot of different backgrounds, for example, biology, computer science, physics, and so on [3–8].
There are a plenty of existing works contributing to the flocking problems. Three heuristic rules leading to emergence of the first computer animation of flocking were first reported by Reynolds in 1987 [5]. The essential flocking rules depict how a personal agent maneuvers based on the local flock mates’ positions as well as velocities. Vicsek et al. designed a simple flocking model of multiagents which can all move with the same speed but with different directions in the plane. The Vicsek model is a special version of a pattern introduced previously by Reynolds [6]. Jadbabaie et al. put forward the rigorous proof of the convergence for Vicsek model [7]. There are many generalized species of this model, such as a leader follower strategy which means one agent acts as a group leader and the other agents would just follow the leader and keep the cohesion/separation/alignment rules. Recently, various desired state flocking motions are deduced in most cases [9–13]. Zavlanos et al., especially, proposed connectivity preserving controllers, by designing novel interagent potentials, to realize the flocking of multiagent system under the initial connectivity assumption [12, 13].
The time delays of systems are a very common phenomenon in real life. Many factors, for example, finite signal transmission speeds and memory effects, can cause time delay in spreading and communication. Therefore, it should be considered to design the control scheme for multiagent system with time delay. The effect of exchange delays for consensus problems and formation problems has been discussed [14–19]. According to the matrix theory and the frequency analysis, Su et al. obtained the desired moving model with a delay-dependent formation control algorithm [20]. Yang et al. studied the virtual potential approach for stabilising a group of agents at a desired formation [21]. Adaptive control is a kind of very important method in the control of complex nonlinear systems [22–24]. A good adaptive control can adapt to the changes in a large range of parameters of controlled system. It can not only maintain stable operation of the system but also keep the optimal in degree. The literatures focus on flocking with collision avoidance. However, the distance between multiagents is required to be larger than a constant r in reality. For example, bird to incite wings must have their own space. Robot teams and UAV (unmanned air vehicles) in the formation movement in order to avoid a collision must consider their size. So, in the flocking control, it is not enough only to require the distance between them to be greater than zero. Motivated by this fact and on the basis of the abovementioned works, distance constrained based adaptive flocking control for multiagent system with time delay is presented in this paper. The innovation of this paper is mainly in the following aspects: the adaptive controller being designed to achieve the adaptive flocking of multiagents with time delays and keeping the distance between multiagents to be larger than a constant r. The stability of the adaptive flocking of the multiagent system with time delays is analyzed theoretically. A sufficient condition is given for the stability of the adaptive control system.
The rest of the paper was structured from the following aspects. In Section 2, the multiagent flocking problem and some preliminaries used throughout this paper are introduced. In Section 3, a controller is designed based on adaptive flocking control laws. In Section 4, the main theory results that velocities of each agent come to the same asymptotically and the distance between agents required to be larger than a constant r are proved. In Section 5, a simulation case is also presented to verify the effectiveness of our theoretical results. In Section 6, the full text content is summarized and the further research in the aspects of adaptive flocking of multiagent networks is investigated.
2. Preliminaries
A set of N(N≥2) agents moving in an n-dimensional Euclidean space are considered. The dynamics of each agent is characterized by the following dynamic system (see [25]):(1)x˙it=vit,v˙it=fvit+uit,i=1,2,…,N,where xi(t)∈R2n and vi(t)∈R2n are the position vector and the velocity of agent i, independently, f(vi(t)) is a nonlinear vector-valued continuous function which describes the intrinsic dynamics of agent i, and ui(t)∈R2n is the control input acting on agent i. Especially, the virtual leader for multiagent system (1) is a special agent described by(2)x˙lt=vlt,v˙lt=fvlt,where xl(t)∈R2n and vl(t)∈R2n are the position and velocity vector of the virtual leader, respectively. In this paper, an assumption has been made that all agents can get the information of the virtual leader. The information switching between multiagents with the leader exhibits time delays.
Define error vector(3)ei=vi-vl,i=1,2,…,N.
Based on the definition of ei, we have the following equations:(4)vjt-τ-vlt-τ=ejt-τi,j=1,2,…,N,(5)ei-ej=vi-vl-vj-vl=vi-vji,j=1,2,…,N.
Assumption 1.
Assume that there exists a nonnegative constant α satisfying (6)fvi-fvl2≤αei2,i=1,2,…,N.
Definition 2.
Communication radius R2 is defined as the biggest distance from which multiagents can get information of other agents. Hysteresis radius R1 is the distance in which a new edge will not be added to the graph until the distance between any two agents which are not connected decreases to R1. Safe radius r is the distance needed by the agent own activities.
Definition 3 (dynamic graphs [10]).
We call G(t)=(V,E(t)) a dynamic graph which consists of a set of vertices V={1,2,…,N} indexed by the set of agents and a time varying set of links Et={(i,j)∣i,j∈V}, such that, for any 0<r<R1<R2,
if r<‖xi(t)-xj(t)‖≤R1, then (i,j)∈E(t),
if R2≤‖xit-xj(t)‖, then (i,j)∉E(t).
Dynamic graphs G(t) meet the conditions that (i,j)∈E(t) if and only if (j,i)∈E(t) are called undirected, which constitute the key point of the paper. If any vertices i and j in an undirected graph G(t) are joined by a link (i,j)∈E(t), we can call them adjacent or neighbors at time t and denote them by i~j. Let L(t)=A(t)-Δ(t) be equivalently represented by a dynamic negative Laplacian matrix, where A(t)=(aij(t)) is a weighted adjacency matrix of the graph G(t), which satisfies that aij(t)>0 if i~j, aii=-∑j=1,j≠iNaij, and aij(t)=0, otherwise. A topological invariance of graphs, that is, graph connectivity, is of great interest for this paper. The switching process of dynamic graphs can be shown in Figure 1.
The switching process of dynamic graphs according to Definition 3.
Definition 4 (graph connectivity).
We say that a dynamic graph G(t) is connected at time t if there exists a path, that is, a sequence of distinct vertices such as the fact that consecutive vertices are adjacent, between any two vertices in G(t).
Hence, the problem which is mentioned can be formally stated as follows.
Definition 5 (flocking).
Consider the set of connected graphs CN on N vertices determining the distribution of control laws ui(t) for all agents i. If the initial network G(t0)∈CN and the initial distance of each agent is bigger than r, G(t)∈CN for all time t≥t0, and the distance of multiagents satisfies ‖xij‖>r. Moreover, the velocity of each agent asymptotically approaches the same.
3. The Design for Virtual Leader Based Adaptive Flocking Control Laws
In view of any dynamic graph G(t)=(V,E(t)), we can define the set of control laws,(7)ui=-2∑j=1Naij∇xiVij-2∑j=1,j≠iNaijvit-τ-vjt-τ-diei,iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii=1,2,…,N,d˙i=kieiTei,where di is the updating parameters, {ki,i=1,2,…,N} are positive constants, τ is the coupling time delay, and ∇xiVij is a direction vector of the negative gradient of an artificial potential function defined by the following equation:(8)Vijxij=1R22-xij2+1xij2-r2,xij∈r,R2with xij=xi-xj, where Vij grows unbound when ‖xij‖→R2-, or ‖xij‖→r+.
Based on the definition of Vij, we have the following equation:(9)∇xijVijxij=∇xiVijxij=-∇xjVijxijiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii,j=1,2,…,N.
From equations of (4), (7) and aii=-∑j=1,j≠iNaij, it was concluded that (10)ui=-2∑j=1Naij∇xiVij-2∑j=1,j≠iNaijvit-τ-vjt-τ-diei=-2∑j=1Naij∇xiVij-2∑j=1,j≠iNaijvit-τ+2∑j=1,j≠iNaijvjt-τ-diei=-2∑j=1Naij∇xiVij+2aiivit-τ+2∑j=1,j≠iNaijvjt-τ-diei=-2∑j=1Naij∇xiVij+2∑j=1Naijvjt-τ-diei=-2∑j=1Naij∇xiVij+2∑j=1Naijvjt-τ-2∑j=1Naijvlt-τ-diei=-2∑j=1Naij∇xiVij+2∑j=1Naijejt-τ-dieiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii=1,2,…,N.
From (1), (2), and (10), we can get the error equation(11)e˙i=fvi-fvl-2∑j=1Naij∇xiVij+2∑j=1Naijejt-τ-dieiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii=1,2,…,N.
We will give some main results for our proposed scheme in next chapter.
4. The Main Theory Results
For the considered multiagent system, we can assume that G(t) switches at time tp for p=1,2,…. That is to say, G(t) is a fixed graph in each time interval [tp-1,tp). At the same time, we can define a switching signal g(t):[t0,∞)→CN associated with connected graphs, and we can have the following conclusion.
Theorem 6.
For the multiagent systems (1) and (2), the switching signal g(t) satisfies E(tp)⊆E(tq) for any pair of switching times tp≤tq under the control laws (10).
Proof.
Consider the following semipositive definite function:(12)VG=12∑i=1NeiTei+12∑i=1Ndi-d^i2ki+∑i=1N∑j=1NaijVij+∑i=1N∫-τ0eiTt+ηeit+ηdη,where d^i is a sufficiently large positive constant.
The generalized time derivative of VG is(13)V˙G=∑i=1NeiTe˙i+∑i=1Ndi-d^ikid˙i+∑i=1N∑j=1NaijV˙ij+∑i=1NeiTei-eiTt-τeit-τ.Due to (9), we have (14)∑i=1N∑j=1NaijV˙ij=∑i=1N∑j=1Naijx˙ijT∇xijVij=∑i=1N∑j=1Naijx˙iT∇xijVij-x˙jT∇xijVij=∑i=1N∑j=1Naijx˙iT∇xiVij+x˙jT∇xjVij=∑i=1N∑j=1NaijviT∇xiVij+vjT∇xjVij=∑i=1N∑j=1NaijviT-vjT∇xiVij.Then, from (5), we have (15)∑i=1N∑j=1NaijviT-vjT∇xiVij=∑i=1N∑j=1NaijeiT-ejT∇xiVij=∑i=1N∑j=1NaijeiT∇xiVij+aijejT∇xjVij=2∑i=1N∑j=1NaijeiT∇xiVij.The first two parts of (13) can be rewritten as the following forms: (16)∑i=1NeiTe˙i+∑i=1Ndi-d^ikid˙i=∑i=1NeiTfvi-fvl-2∑j=1Naij∇xiVijmmmmimm+2∑j=1Naijejt-τ-diei+∑i=1Ndi-d^ieiTei=∑i=1NeiTfvi-fvl-2∑i=1N∑j=1NaijeiT∇xiVij+2∑i=1N∑j=1NaijeiTejt-τ-∑i=1NdieiTei+∑i=1Ndi-d^ieiTei=∑i=1NeiTfvi-fvl-2∑i=1N∑j=1NaijeiT∇xiVij+2∑i=1N∑j=1NaijeiTejt-τ-∑i=1NeiTd^iei=-2∑i=1N∑j=1NaijeiT∇xiVij+2∑i=1N∑j=1NaijeiTejt-τ+∑i=1NeiTfvi-fvl-d^iei≤-2∑i=1N∑j=1NaijeiT∇xiVij+2∑i=1N∑j=1NaijeiTejt-τ+∑i=1Nα-d^ieiTei.Thus, we can obtain (17)V˙G≤-2∑i=1N∑j=1NaijeiT∇xiVij+2∑i=1N∑j=1NaijeiTejt-τ+∑i=1Nα-d^ieiTei+2∑i=1N∑j=1NaijeiT∇xiVij+∑i=1NeiTei-eiTt-τeit-τ=2∑i=1N∑j=1NaijeiTejt-τ+∑i=1NeiTei-eiTt-τeit-τ+∑i=1Nα-d^ieiTei≤2∑i=1N∑j=1NaijeiTejt-τ-∑i=1NeiTt-τeit-τ+∑i=1NeiTei+∑i=1Nα-d^ieiTei=2eTΓet-τ-et-τTet-τ+eTΛe≤eTΓTΓe+eTΛe=eTΓTΓ+Λe,where (18)e=e12,e22,…,eN2T,et-τ=e1t-τ2,e2t-τ2,…,eNt-τ2T,Γ=aijN×N,Λ=α-d^1+10⋱0α-d^N+1.
As α and |aij| are nonnegative, one can select suitable positive constants d^i(1≤i≤N) to make(19)Λ+ΓTΓ≤0.So, it follows that V˙G≤0. Hence, V˙G is negative semidefinite for any signal g(t).
Because V˙G is negative semidefinite, VG cannot increase for all time. In addition, the VG is bounded. In other words, it means that, for any (i,j)∈E(t), Vij is bounded. On the one hand, if for some (i,j)∈E(t), ‖xij‖→R2, Vij(xij)→∞. Thus, by continuity of VG, it follows that ‖xij‖<R2, for all (i,j)∈E(t) and t∈[tp,tp+1]. In other words, all links in VG are preserved between switching times, which implies that E(tp)⊆E(tq). Apply this recursive argument to complete the proof. A similar argument for the case where xij→r+ can be used to establish that the distance of multiagents satisfies ‖xij‖>r for all time t≥t0.
Clearly, Theorem 6 explains that the switching signal will satisfy G(t)∈CN for all time t>t0, if G(t0)∈CN. We, specially, have the following corollary.
Corollary 7.
Under control law (7), the total number of switching times of system (1) is finite.
Proof.
The size of the set of links E(t0) forms an increasing sequence and (20)supt≥t0Et-Et0=nn-12-n-1,where n-1 is the number of links in G(t0) if it is minimally connected, that is, if it is a tree, and n(n-1)/2 corresponds to the number of links in a complete graph.
Theorem 8.
Take a group of N agents with dynamics (1), (2), each steered by protocol (10) into account. Suppose that initial network G(t) is connected and the initial distance of each agent is bigger than r. Then, G(t) is connected for all t≥0, the velocity of each agent asymptotically approaches the same, and the distance of each agent is bigger than r for all t≥0.
Proof.
The number of switching times of the closed loop system is finite by Corollary 7. And by Theorem 6 if G(t0) is connected, we can draw the conclusion that G(t) is connected for all time t>t0 and so eventually G(t)→G∈CN.
Referring to the analysis method of [22–24], we will give the proof of the conclusions that the velocity of each agent asymptotically approaches the same and the distance of each agent is bigger than r for all t≥0.
From (17) and (19), we know that VG is monotonically decreasing and having a lower bound. So, it concludes that VG asymptotically converges to a finite nonnegative value. And, from (17), we have V˙G≤eT(Λ+ΓTΓ)e. When d^i,i=1,2,…,N are sufficiently large positive constants, we can see that there exists a nonnegative constant λ satisfying(21)V˙G≤-λeTe.Thus, we have(22)e2≤-V˙Gλ.
Because limt→+∞∫0t-V˙G/λdη exists, we can derive from (22) that limt→+∞∫0te2dη exists by the generalized integral principle.
We can define (23)VG1=12∑i=1Ndi-d^i2ki,VG2=∑i=1N∑j=1NaijVij,VG3=∑i=1N∫-τ0eiTt+ηeit+ηdη.Then, we can get the following equation:(24)VG=12∑i=1NeiTei+VG1+VG2+VG3.
The generalized time derivative of VG1 is V˙G1=∑i=1N(di-d^i)eiTei. It is negative when {d^i,i=1,2,…,N} are sufficiently large positive constants. We can conclude that VG1 asymptotically converges to a finite nonnegative value. On the other hand, VG2 asymptotically converges to a finite nonnegative value based on the definition of Vij (refer to (8)). At the same time, the generalized time derivative of VG3 is V˙G3=∑i=1N[eiTei-eiT(t-τ)ei(t-τ)] which asymptotically converges to zero. Therefore, we can derive that VG3 asymptotically converges nonnegative value. By the boundedness of VG and (24), therefore, we can derive that limt→+∞∑i=1NeiTei exists and is a nonnegative real number. In what follows, we will prove that(25)limt→+∞∑i=1NeiTei=0.
If this is not true, we have (26)limt→+∞∑i=1NeiTei=μ1>0.
Then, there obviously exists a real number M>0 and 0<μ<μ1, such that ∑i=1NeiTei>μ>0 for t≥M. From (21), we can get(27)V˙G≤-λ∑i=1NeiTei<-λμ,t≥M.By integrating (25) with respect to t over the time period M to +∞, we can obtain (28)-VGM≤VG+∞-VGM=∫M+∞V˙Gtdt<-∫M+∞λμdt=-∞.This yields a contradiction, and so (29)limt→+∞∑i=1NeiTei=0which implies that limt→+∞ei=0. Moreover, it is natural to see limt→+∞vi=vl. Thus, in the steady state, all agent velocities asymptotically become the same. Furthermore, in the steady state, since ei=0 for i=1,2,…,N, we can get(30)ddtxi-xj2=2xi-xjTvi-vj=0,and so the distances between agents are invariant.
Finally, we prove that the distance of each gent is bigger than r for all t≥0. From above proof, we have that VG is bounded for all t≥t0. Also from the definite of potential function, if xij→r, we can get that VG is unbounded. Therefore, the distance of each agent is bigger than r for all t≥0.
5. Simulation Results
In this section, an example was presented to show the effectiveness of our proposed algorithm. The potential function Vij is defined in (8) with R2=5, r=1, and with hysteresis R1=4.5. The simulation is performed with 20 agents in R22. Each agent’s initial velocities, labeled with dots, are randomly selected in the unit square. The initial positions distributed on the perimeter of a circle of radius d=1.5 on average (Figure 2). In particular, we assume the velocity is v0(t)=[sin(t),cot(t)]T and send a signal to one agent which is selected randomly in this simulation. Solid curves connecting every agent indicate the vicinal traveled paths, while arrows indicate the agents velocities. The networks connectivity is guaranteed for all time and the groups asymptotic flocking is achieved, where τ=0.5. Figure 2 is the initial state. Figure 3 indicates the final steady state configuration. Figure 4 is the velocity error plot which is used for describing the errors between agent actual velocities, where the velocity error curves were indicated by the solid lines, and it demonstrates that all agent velocities asymptotically tend to the same. Figure 5 is changing curve of control laws. Figure 6 is the distance ‖xij‖ between agent i and agent j and agent i,j∈N(t). From which one can see that the distance between all agents is larger than r=1 for all t≥0. This simulation results have well verified the theory results.
Initial configuration.
End configuration.
Velocity errors.
Changing curve of control laws.
Distance of multiagents.
6. Conclusions
In this paper, the adaptive flocking of multiagents with time delay is studied. A novel adaptive flocking control method for multiagents is proposed, and the control law is designed depending on functions of the state information and the external signal. By the control law, all agents can follow the virtual leader and can ensure freedom from collisions between neighboring agents. Some theoretical results are attained, and a numerical example is given to show the practicability of the proposed method. The distance constrained based adaptive flocking control for multiagent system with time delay is proposed in this paper, which has received the expected results. But there are still some problems to be resolved in future research. Firstly, the research of the distance constrained based adaptive flocking control for multiagent system with time delay is achieved on the basics of the strongly connected network. We can study the adaptive flocking control for multiagent system based on the networks which contains a directed spanning tree. Secondly, we can apply the algorithm proposed in this paper to the concrete platform of multiple mobile robots or wireless sensor network, which can combine theory with practice better.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work is supported by National Natural Science Foundation of China (61174094, 61039001), National Natural Science Foundation of Tianjin (14JCYBJC18700), Technology and Innovation Fund Project of Civil Aviation University Of China (SY-1415), and Basic Research Projects of High Education (ZXH2010D011, ZXH2012 K002).
ShawE.Fish in schools19758484045TangM. S.MaoX. J.GuessoumZ.Research on an infectious disease transmission by flocking birds20132013719682310.1155/2013/1968232-s2.0-84880167204DoK. D.Flocking for multiple elliptical agents with limited communication ranges201127593194210.1109/TRO.2011.21594102-s2.0-80053648192FerraraA.VecchioC.Collision avoidance strategies and coordinated control of passenger vehicles200749447549210.1007/s11071-006-9110-42-s2.0-34547398091ReynoldsC. W.Flocks, herds, and schools: a distributed behavioral model21Proceedings of the 14th Annual Conference on Computer Graphics and Interactive Techniques, (SIGGRAPH '87)1987ACM2534VicsekT.CzirókA.Ben-JacobE.CohenI.ShochetO.Novel type of phase transition in a system of self-driven particles19957561226122910.1103/PhysRevLett.75.12262-s2.0-0001730998JadbabaieA.LinJ.MorseA. S.Coordination of groups of mobile autonomous agents using nearest neighbor rules2003486988100110.1109/TAC.2003.812781MR19862662-s2.0-0038548185TonerJ.TuY.RamaswamyS.Hydrodynamics and phases of flocks2005318117024410.1016/j.aop.2005.04.011MR2148645ZBL1126.823472-s2.0-20444507542Olfati-SaberR.Flocking for multi-agent dynamic systems: algorithms and theory200651340142010.1109/TAC.2005.864190MR22056792-s2.0-33644997478TannerH. G.JadbabaieA.PappasG. J.Stable flocking of mobile agents, part I: fixed topologyProceedings of the 42nd IEEE Conference on Decision and ControlDecember 2003Maui, Hawaii, USA201020152-s2.0-1542648889TannerH. G.JadbabaieA.PappasG. J.Stable flocking of mobile agents, part II: dynamic topologyProceedings of the IEEE Conference on Decision and Control2003Maui, Hawaii, USA20162021ZavlanosM. M.JadbabaieA.PappasG. J.Flocking while preserving network connectivityProceedings of the 46th IEEE Conference on Decision and Control (CDC '07)December 2007New Orleans, La, USA2919292410.1109/CDC.2007.44345302-s2.0-62749130828TannerH. G.JadbabaieA.PappasG. J.Flocking in fixed and switching networks200752586386810.1109/TAC.2007.895948MR23242462-s2.0-34247397566MoreauL.Stability of continuous-time distributed consensus algorithms2004School of Electronics and Computer Science, University of Southamptonhttp://arxiv.org/pdf/math/0306426.pdfOlfati-SaberR.MurrayR. M.Consensus problems in networks of agents with switching topology and time-delays20044991520153310.1109/TAC.2004.834113MR20869162-s2.0-4644244041HuJ. P.HongY. G.Leader-following coordination of multi-agent systems with coupling time delays2007374285386310.1016/j.physa.2006.08.0152-s2.0-33751083607XiaoF.WangL.State consensus for multi-agent systems with switching topologies and time-varying delays200679101277128410.1080/00207170600825097MR2252185ZBL050626272-s2.0-33746219707LiuC.-L.TianY.-P.Formation control of multi-agent systems with heterogeneous communication delays200940662763610.1080/00207720902755762MR25410002-s2.0-67651211614RistianS.CesareF.Formation control over delayed communication networksProceedings of the IEEE International Conference on Robotics and AutomationMay 2008Pasadena, Calif, USA536538SuH.ZhangN.ChenM. Z.WangH.WangX.Adaptive flocking with a virtual leader of multiple agents governed by locally Lipschitz nonlinearity201314179880610.1016/j.nonrwa.2012.06.010MR29698742-s2.0-84866398950YangZ.ZhangQ.JiangZ.ChenZ.Flocking of multi-agents with time delay201243112125213410.1080/00207721.2011.564675MR2982561WangJ.-L.WuH.-N.GuoL.Novel adaptive strategies for synchronization of linearly coupled neural networks with reaction-diffusion terms201425242944010.1109/TNNLS.2013.22760862-s2.0-84893695040WangJ. L.WuH. N.Synchronization and adaptive control of an array of linearly coupled reaction-diffusion neural networks with hybrid coupling201444813501361WangJ.-L.WuH.-N.Adaptive output synchronization of complex delayed dynamical networks with output coupling201414217418110.1016/j.neucom.2014.04.0502-s2.0-84901863271WangJ.-L.WuH.-N.Leader-following formation control of multi-agent systems under fixed and switching topologies201285669570510.1080/00207179.2012.662720MR29223602-s2.0-84862070311