<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">ICA</journal-id><journal-title-group><journal-title>Intelligent Control and Automation</journal-title></journal-title-group><issn pub-type="epub">2153-0653</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/ica.2016.71002</article-id><article-id pub-id-type="publisher-id">ICA-63383</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Computer Science&amp;Communications</subject></subj-group></article-categories><title-group><article-title>
 
 
  Flocking for Leader-Following Multi-Agent Systems with Time-Varying Delay
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>englan</surname><given-names>Sun</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Rui</surname><given-names>Wang</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Yongfu</surname><given-names>Li</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Feng</surname><given-names>Liu</given-names></name><xref ref-type="aff" rid="aff3"><sup>3</sup></xref></contrib></contrib-group><aff id="aff3"><addr-line>School of Automation, China University of Geosciences, Wuhan, China</addr-line></aff><aff id="aff2"><addr-line>College of Automation, Chongqing University of Posts and Telecommunications, Chongqing, China</addr-line></aff><aff id="aff1"><addr-line>Research Center of System Theory and Application, Chongqing University of Posts and Telecommunications, Chongqing, China</addr-line></aff><pub-date pub-type="epub"><day>04</day><month>01</month><year>2016</year></pub-date><volume>07</volume><issue>01</issue><fpage>9</fpage><lpage>15</lpage><history><date date-type="received"><day>28</day>	<month>November</month>	<year>2015</year></date><date date-type="rev-recd"><day>accepted</day>	<month>3</month>	<year>February</year>	</date><date date-type="accepted"><day>6</day>	<month>February</month>	<year>2016</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  This paper investigates the flocking problem in multi-agent system with time-varying delay and a virtual leader. Each agent here is subject to nonlinear dynamics. For the system, the corresponding algorithm with time-varying delay is proposed. Under the assumption that the initial network is connected, it is proved that the distance between agents is in the desired distance. The theoretical deduction shows that the stable flocking motion is achieved.
 
</p></abstract><kwd-group><kwd>Flocking</kwd><kwd> Time-Varying Delay</kwd><kwd> Nonlinear Dynamics</kwd><kwd> Multi-Agent System</kwd><kwd> Virtual Leader</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Flocking is a collective behavior of large number of interacting agents with a common group objective. Examples of these agents include birds, fish, penguins, ants, bees, and crowds. Many scientists from rather diverse disciplines, including physics, mathematics, control engineering and biology, have been interested in flocking problem [<xref ref-type="bibr" rid="scirp.63383-ref1">1</xref>] -[<xref ref-type="bibr" rid="scirp.63383-ref8">8</xref>] . The first well-known flocking model was proposed by Craig Reynolds [<xref ref-type="bibr" rid="scirp.63383-ref2">2</xref>] . Reynolds started with a boid model to build a simulated flock and introduced three rules (i.e., separation, cohesion and alignment rules) for flocking. Based on Reynolds’ three rules, flocking problems have been investigated from various perspectives [<xref ref-type="bibr" rid="scirp.63383-ref9">9</xref>] -[<xref ref-type="bibr" rid="scirp.63383-ref15">15</xref>] . In [<xref ref-type="bibr" rid="scirp.63383-ref16">16</xref>] , an artificial potential function is put forward and three algorithms are introduced. It provides a theoretical framework for the designing of flocking algorithms. Multi-agent flocking under topological interactions is considered, which define a notion of hierarchical structure in the interaction graph that establish conditions building upon previous work on multi-agent systems with switching communication networks in [<xref ref-type="bibr" rid="scirp.63383-ref17">17</xref>] . Using structurally balanced signed graph theory and a specified potential function, a stable bipartite flock formation is achieved for both virtual leader and leaderless situations in [<xref ref-type="bibr" rid="scirp.63383-ref18">18</xref>] . However, there is a common assumption that virtual leaders guide the flocking behaviors of the group. In this paper, leader follower flocking problem of multi-agent system is considered.</p><p>A flocking problem concerning multiple leaders in which followers use the position of flocking center to keep their connections is studied in [<xref ref-type="bibr" rid="scirp.63383-ref19">19</xref>] . In [<xref ref-type="bibr" rid="scirp.63383-ref20">20</xref>] , two leader-follower adaptive flocking algorithms are proposed with the combination of consensus and attraction/repulsion function respectively to solve the cohesive flocking problem and the formation flocking problem. Aiming at the group of autonomous agents consisting of multiple leader agents and multiple follower ones, a flocking behavior method with multiple leaders and a global trajectory was proposed in [<xref ref-type="bibr" rid="scirp.63383-ref21">21</xref>] . Yu et al., [<xref ref-type="bibr" rid="scirp.63383-ref22">22</xref>] give a distributed leader-follower algorithm considering the group consisting of one leader. In [<xref ref-type="bibr" rid="scirp.63383-ref23">23</xref>] , for the circumstance with a virtual leader, the agents would follow the virtual leader and achieve the same velocity asymptotically.</p><p>In practice, time delay is inevitable and would damage the stability of system. Jing et al. [<xref ref-type="bibr" rid="scirp.63383-ref23">23</xref>] investigate flocking problem of multi-agent systems with time delay and discuss systems with homogeneous and inhomogeneous time delay. Yang et al. [<xref ref-type="bibr" rid="scirp.63383-ref24">24</xref>] proposed an adaptive flocking algorithm for multi-agent system with time delay. It is proved that the distance between agents can be larger than a constant during the motion evolution by using the flocking algorithm. The authors investigate the flocking problem of multi-agent systems led by one active virtual leader with a directed topology containing time-varying coupling delays, which based only on the three classical assumptions for flocking systems in [<xref ref-type="bibr" rid="scirp.63383-ref25">25</xref>] . Because of these problems, this paper will study the flocking problems in the multi-agent system with a virtual leader and time-varying delay.</p><p>The rest of this paper is organized as follows. Some basic preliminaries and flocking algorithms are presented in Section 2. Section 3 gives the nonlinear leader-following multi-agent models. Algorithms and main results are presented in Section 4. Section 5 concludes the paper and offers suggestions for future work.</p></sec><sec id="s2"><title>2. Preliminaries</title><p>In this section, some related preliminary knowledge are introduced. For any vector <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x6.png" xlink:type="simple"/></inline-formula> denotes its transpose and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x7.png" xlink:type="simple"/></inline-formula> denotes the Euclidean norm. Let <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x8.png" xlink:type="simple"/></inline-formula> be a weighted undirected graph with the set of nodes <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x9.png" xlink:type="simple"/></inline-formula> and the set of agents <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x10.png" xlink:type="simple"/></inline-formula></p><p>Graphs with self-loops will not be considered in this paper. The weight adjacency matrix <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x11.png" xlink:type="simple"/></inline-formula> where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x12.png" xlink:type="simple"/></inline-formula> if <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x13.png" xlink:type="simple"/></inline-formula> otherwise, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x14.png" xlink:type="simple"/></inline-formula>An edge denoted by the pair <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x15.png" xlink:type="simple"/></inline-formula> represents a commu- nication link from node j to i. A path <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x16.png" xlink:type="simple"/></inline-formula> from node i to node j is a sequence of edges, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x17.png" xlink:type="simple"/></inline-formula> in which all nodes <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x18.png" xlink:type="simple"/></inline-formula> are distinct. An undirected graph is called connected if there is a path between each pair of distinct nodes. <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x19.png" xlink:type="simple"/></inline-formula>is the degree matrix whose diagonal elements are defined by <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x20.png" xlink:type="simple"/></inline-formula> The Laplacian matrix of graph <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x20.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x21.png" xlink:type="simple"/></inline-formula> is <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x20.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x21.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x22.png" xlink:type="simple"/></inline-formula> Then it has following properties [<xref ref-type="bibr" rid="scirp.63383-ref26">26</xref>] ,</p><p>1) The eigenvalues of L satisfy <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x23.png" xlink:type="simple"/></inline-formula> If the graph <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x23.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x24.png" xlink:type="simple"/></inline-formula> is connected, there is</p><disp-formula id="scirp.63383-formula176"><graphic  xlink:href="http://html.scirp.org/file/2-7900437x25.png"  xlink:type="simple"/></disp-formula><p>2) The Laplacian matrix L is a positive semi-definite matrix that satisfies the following sum-of-squares property:</p><disp-formula id="scirp.63383-formula177"><graphic  xlink:href="http://html.scirp.org/file/2-7900437x26.png"  xlink:type="simple"/></disp-formula><p>Lemma 1. [<xref ref-type="bibr" rid="scirp.63383-ref4">4</xref>] Suppose that <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x27.png" xlink:type="simple"/></inline-formula> is an undirected graph of order N, and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x28.png" xlink:type="simple"/></inline-formula> is a graph generated by adding some edges into the graph<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x29.png" xlink:type="simple"/></inline-formula>. Then <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x30.png" xlink:type="simple"/></inline-formula> where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x30.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x31.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x30.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x31.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x32.png" xlink:type="simple"/></inline-formula> are the symmetric Laplacian matrices of graphs <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x30.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x31.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x32.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x33.png" xlink:type="simple"/></inline-formula> and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x30.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x31.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x32.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x33.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x34.png" xlink:type="simple"/></inline-formula>, respectively.</p><p>Lemma 2. For any vectors <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x35.png" xlink:type="simple"/></inline-formula> the following matrix inequality holds: <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x36.png" xlink:type="simple"/></inline-formula></p></sec><sec id="s3"><title>3. Problems Formulation</title><p>Consider the multi-agent system described by</p><disp-formula id="scirp.63383-formula178"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-7900437x37.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x38.png" xlink:type="simple"/></inline-formula> are the position and velocity states of i<sup>th</sup> agent, respectively. <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x38.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x39.png" xlink:type="simple"/></inline-formula> is the nonlinear dynamic of agent i and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x38.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x39.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x40.png" xlink:type="simple"/></inline-formula> is the control input. Denote <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x38.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x39.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x40.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x41.png" xlink:type="simple"/></inline-formula> as the relative distance between agent i and agent j.</p><p>For the systems with virtual leader available, the dynamics of virtual leader is described as</p><disp-formula id="scirp.63383-formula179"><label>(2)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-7900437x42.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x43.png" xlink:type="simple"/></inline-formula> represent the position, velocity and control input vector of the virtual leader.</p><p>Assumption (A): There exists a positive constant <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x44.png" xlink:type="simple"/></inline-formula> satisfying</p><disp-formula id="scirp.63383-formula180"><label>(3)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-7900437x45.png"  xlink:type="simple"/></disp-formula><p>Supposed that all agents have the same sensing radius <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x46.png" xlink:type="simple"/></inline-formula> Then the neighboring set of agent i is denoted as <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x46.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x47.png" xlink:type="simple"/></inline-formula> Since the size of agent cannot be ignored usually that a minimum allowable distance r (collision distance) is considered in the model.</p><p>Definition 1: Given a constant<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x48.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x48.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x49.png" xlink:type="simple"/></inline-formula>is called a dynamic undirected graph with a time-varying set of links <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x48.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x49.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x50.png" xlink:type="simple"/></inline-formula> such that</p><p>1) Initial links are generated by <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x51.png" xlink:type="simple"/></inline-formula></p><p>2) If <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x52.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x52.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x53.png" xlink:type="simple"/></inline-formula> then <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x52.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x53.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x54.png" xlink:type="simple"/></inline-formula> is a new link to be added<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x52.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x53.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x54.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x55.png" xlink:type="simple"/></inline-formula>. It is called hysteresis effect and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x52.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x53.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x54.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x55.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x56.png" xlink:type="simple"/></inline-formula> is the hysteresis distance, which is crucial in preserving connectivity of the network;</p><p>3) If <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x57.png" xlink:type="simple"/></inline-formula> then<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x57.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x58.png" xlink:type="simple"/></inline-formula>.</p><p>The neighboring set of agent i is divided into four regions, named collision region, separation region, alignment region and attraction region, in which <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x59.png" xlink:type="simple"/></inline-formula> If <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x59.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x60.png" xlink:type="simple"/></inline-formula> agents i and j are in desired distance.</p><p>Definition 2. <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x61.png" xlink:type="simple"/></inline-formula>is a bounded function with respect to <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x61.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x62.png" xlink:type="simple"/></inline-formula> between agents i and j, which satisfies</p><p>1) <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x63.png" xlink:type="simple"/></inline-formula>when <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x63.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x64.png" xlink:type="simple"/></inline-formula></p><p>2) <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x65.png" xlink:type="simple"/></inline-formula>when <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x65.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x66.png" xlink:type="simple"/></inline-formula></p><p>3) <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x67.png" xlink:type="simple"/></inline-formula>otherwise.</p><p>Definition 3. The pairwise bounded potential function <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x68.png" xlink:type="simple"/></inline-formula> can be defined.</p><disp-formula id="scirp.63383-formula181"><label>(4)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-7900437x69.png"  xlink:type="simple"/></disp-formula><p>which satisfies</p><p>1) <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x70.png" xlink:type="simple"/></inline-formula>decreases with the increase of z when<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x70.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x71.png" xlink:type="simple"/></inline-formula>.</p><p>2) <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x72.png" xlink:type="simple"/></inline-formula>is increasing on<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x72.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x73.png" xlink:type="simple"/></inline-formula>. Obviously, the potential function <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x72.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x73.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x74.png" xlink:type="simple"/></inline-formula> reaches its minimum value 0 when <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x72.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x73.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x74.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x75.png" xlink:type="simple"/></inline-formula></p><p>According to conditions, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x76.png" xlink:type="simple"/></inline-formula>can be constructed as follows:</p><disp-formula id="scirp.63383-formula182"><graphic  xlink:href="http://html.scirp.org/file/2-7900437x77.png"  xlink:type="simple"/></disp-formula><p>where the parameters <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x78.png" xlink:type="simple"/></inline-formula> are positive constants. There are<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x78.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x79.png" xlink:type="simple"/></inline-formula>, and</p><disp-formula id="scirp.63383-formula183"><graphic  xlink:href="http://html.scirp.org/file/2-7900437x80.png"  xlink:type="simple"/></disp-formula></sec><sec id="s4"><title>4. Algorithms and Main Results</title><p>For system (1) with virtual leader (2), the flocking algorithm can be described by</p><disp-formula id="scirp.63383-formula184"><label>(5)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-7900437x81.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x82.png" xlink:type="simple"/></inline-formula> are positive constants. Denote <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x82.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x83.png" xlink:type="simple"/></inline-formula> respectively. Then the system can be described by</p><disp-formula id="scirp.63383-formula185"><label>(6)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-7900437x84.png"  xlink:type="simple"/></disp-formula><p>The control input (6) can be equivalently rewritten as</p><disp-formula id="scirp.63383-formula186"><label>(7)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-7900437x85.png"  xlink:type="simple"/></disp-formula><p>Denote <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x86.png" xlink:type="simple"/></inline-formula></p><p>Definition 4. Flocking motion with a virtual leader is said to be achieved asymptotically for systems (1) and (2), if for any initial state, there is<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x87.png" xlink:type="simple"/></inline-formula> <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x87.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x88.png" xlink:type="simple"/></inline-formula> <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x87.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x88.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x89.png" xlink:type="simple"/></inline-formula> <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x87.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x88.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x89.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x90.png" xlink:type="simple"/></inline-formula></p><p>To demonstrate the validity of control protocol (7), the following positive semi-definite function is constructed</p><disp-formula id="scirp.63383-formula187"><label>(8)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-7900437x91.png"  xlink:type="simple"/></disp-formula><p>where<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x92.png" xlink:type="simple"/></inline-formula>.</p><p>Theorem 1. Consider a multi-agent system modeled by dynamics (1) and (2), driven by control protocol (5). Suppose that the network is initially connected and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x93.png" xlink:type="simple"/></inline-formula> is bounded.When</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x94.png" xlink:type="simple"/></inline-formula>, then following statements hold,</p><p>1) <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x95.png" xlink:type="simple"/></inline-formula>is connected for all<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x95.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x96.png" xlink:type="simple"/></inline-formula>;</p><p>2) No collision occurs among agents for all<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x97.png" xlink:type="simple"/></inline-formula>;</p><p>3) Flocking motion with a virtual leader is achieved asymptotically.</p><p>Proof: Denote the topology switching time sequence as <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x98.png" xlink:type="simple"/></inline-formula> Without loss of generality, assume<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x98.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x99.png" xlink:type="simple"/></inline-formula>. Taking the time derivative of the Lyapunov function Q on <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x98.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x99.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x100.png" xlink:type="simple"/></inline-formula> gives</p><disp-formula id="scirp.63383-formula188"><graphic  xlink:href="http://html.scirp.org/file/2-7900437x101.png"  xlink:type="simple"/></disp-formula><p>From Lemma 2, there is</p><disp-formula id="scirp.63383-formula189"><graphic  xlink:href="http://html.scirp.org/file/2-7900437x102.png"  xlink:type="simple"/></disp-formula><p>For a positive constant k, one has<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x103.png" xlink:type="simple"/></inline-formula>, thus</p><disp-formula id="scirp.63383-formula190"><graphic  xlink:href="http://html.scirp.org/file/2-7900437x104.png"  xlink:type="simple"/></disp-formula><p>Assume that <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x105.png" xlink:type="simple"/></inline-formula> satisfy<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x105.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x106.png" xlink:type="simple"/></inline-formula>, there is</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x107.png" xlink:type="simple"/></inline-formula>,</p><p>which implies that</p><disp-formula id="scirp.63383-formula191"><graphic  xlink:href="http://html.scirp.org/file/2-7900437x108.png"  xlink:type="simple"/></disp-formula><p>By definition (2), one has<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x109.png" xlink:type="simple"/></inline-formula>. Therefore, no edge distance will be tend to R for<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x109.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x110.png" xlink:type="simple"/></inline-formula>, implying that no existing edges will be lost before time<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x109.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x110.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x111.png" xlink:type="simple"/></inline-formula>. Hence, new edges must be added into the network at<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x109.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x110.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x111.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x112.png" xlink:type="simple"/></inline-formula>. For a system consists of N agents, there are at most <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x109.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x110.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x111.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x112.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x113.png" xlink:type="simple"/></inline-formula> edges. At the initial instant<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x109.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x110.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x111.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x112.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x113.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x114.png" xlink:type="simple"/></inline-formula>, the system consists of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x109.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x110.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x111.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x112.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x113.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x114.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x115.png" xlink:type="simple"/></inline-formula> edges, then</p><disp-formula id="scirp.63383-formula192"><graphic  xlink:href="http://html.scirp.org/file/2-7900437x116.png"  xlink:type="simple"/></disp-formula><p>Hence there is no edge lost. In addition, from the definite of potential function, one has<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x117.png" xlink:type="simple"/></inline-formula>. Therefore, no collision occurs during<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x117.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x118.png" xlink:type="simple"/></inline-formula>.</p><p>Similar to the above analysis, taking the time derivative of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x119.png" xlink:type="simple"/></inline-formula> on every<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x119.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x120.png" xlink:type="simple"/></inline-formula>. By lemma 1, there is</p><disp-formula id="scirp.63383-formula193"><graphic  xlink:href="http://html.scirp.org/file/2-7900437x121.png"  xlink:type="simple"/></disp-formula><p>one has<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x122.png" xlink:type="simple"/></inline-formula>, which implies <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x122.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x123.png" xlink:type="simple"/></inline-formula> Thus no edge distance will tend to R for<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x122.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x123.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x124.png" xlink:type="simple"/></inline-formula>, implying that no edge will be lost before time <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x122.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x123.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x124.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x125.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x122.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x123.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x124.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x125.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x126.png" xlink:type="simple"/></inline-formula> is finite. Since <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x122.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x123.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x124.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x125.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x126.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x127.png" xlink:type="simple"/></inline-formula> is connected and no edge in <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x122.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x123.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x124.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x125.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x126.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x127.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x128.png" xlink:type="simple"/></inline-formula> is lost, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x122.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x123.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x124.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x125.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x126.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x127.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x128.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x129.png" xlink:type="simple"/></inline-formula>will be connected for all<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x122.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x123.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x124.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x125.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x126.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x127.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x128.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x129.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x130.png" xlink:type="simple"/></inline-formula>. This completes the proof of part (1).</p><p>Similarly, from<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x131.png" xlink:type="simple"/></inline-formula>, deducing that no edge distance will tend to r, for all<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x131.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x132.png" xlink:type="simple"/></inline-formula>. Thus collision is avoided during the whole process. This completes the proof of part (2). To proof part (3), assume that there are <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x131.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x132.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x133.png" xlink:type="simple"/></inline-formula> new edges being added to the evolving network at time<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x131.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x132.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x133.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x134.png" xlink:type="simple"/></inline-formula>. As no edges are lost for<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x131.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x132.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x133.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x134.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x135.png" xlink:type="simple"/></inline-formula>, and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x131.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x132.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x133.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x134.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x135.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x136.png" xlink:type="simple"/></inline-formula>. Therefore, the number of switching times of the system (1) is finite, which implies that the evolving network <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x131.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x132.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x133.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x134.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x135.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x136.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x137.png" xlink:type="simple"/></inline-formula> eventually becomes fixed. Denote the last topology switching as<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x131.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x132.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x133.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x134.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x135.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x136.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x137.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x138.png" xlink:type="simple"/></inline-formula>. Then Q is continuous and monotonously decreasing for<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x131.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x132.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x133.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x134.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x135.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x136.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x137.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x138.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x139.png" xlink:type="simple"/></inline-formula>. Hence the set</p><disp-formula id="scirp.63383-formula194"><graphic  xlink:href="http://html.scirp.org/file/2-7900437x140.png"  xlink:type="simple"/></disp-formula><p>is positively invariant, where</p><disp-formula id="scirp.63383-formula195"><graphic  xlink:href="http://html.scirp.org/file/2-7900437x141.png"  xlink:type="simple"/></disp-formula><p>and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x142.png" xlink:type="simple"/></inline-formula></p><p>Since <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x143.png" xlink:type="simple"/></inline-formula> is connected for all<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x143.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x144.png" xlink:type="simple"/></inline-formula>, one has<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x143.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x144.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x145.png" xlink:type="simple"/></inline-formula>, for all <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x143.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x144.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x145.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x146.png" xlink:type="simple"/></inline-formula> As<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x143.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x144.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x145.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x146.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x147.png" xlink:type="simple"/></inline-formula>, one</p><p>has<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x148.png" xlink:type="simple"/></inline-formula>. thus<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x148.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x149.png" xlink:type="simple"/></inline-formula>. Therefore, the set Ω is compact. It follows from LaSalle’s invariance</p><p>principle that if the initial condition lies in Ω, then the corresponding trajectories will converge to the largest invariant set inside the region</p><disp-formula id="scirp.63383-formula196"><graphic  xlink:href="http://html.scirp.org/file/2-7900437x150.png"  xlink:type="simple"/></disp-formula><p>From (8), <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x151.png" xlink:type="simple"/></inline-formula>if and only if<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x151.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x152.png" xlink:type="simple"/></inline-formula>, which implies that the velocities of all agents will converge to that of the virtual leader asymptotically.</p><p>Since<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x153.png" xlink:type="simple"/></inline-formula>, there is <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x153.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x154.png" xlink:type="simple"/></inline-formula> for all <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x153.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x154.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x155.png" xlink:type="simple"/></inline-formula> From (6), one has</p><disp-formula id="scirp.63383-formula197"><graphic  xlink:href="http://html.scirp.org/file/2-7900437x156.png"  xlink:type="simple"/></disp-formula><p>Thus, unless the inital configuration of the agents is close enough to the global minimum, almost every final configuration locally minimizes each agent’s global potential. which implies</p><disp-formula id="scirp.63383-formula198"><graphic  xlink:href="http://html.scirp.org/file/2-7900437x157.png"  xlink:type="simple"/></disp-formula><p>Then the flocking is achieved. This completes the proof of part (3), thus Theorem 1 hold.</p><p>Remark 1. If <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7900437x158.png" xlink:type="simple"/></inline-formula> is the constant delay, from the deduction above, Theorem 1 is also hold.</p></sec><sec id="s5"><title>5. Conclusion</title><p>This paper mainly discusses the flocking problem of multi-agent system with a virtual leader and time-varying delay. Unlike most existing flocking algorithms, each agent here is subject to nonlinear dynamics. The corresponding algorithms with the time-varying delay are proposed. Under the assumption that the initial network is connected, the theoretical deduction is made. The related topic over the directed network or the jointly connected network will be studied in future.</p></sec><sec id="s6"><title>Acknowledgements</title><p>We thank the Editor and the referee for their comments. This work was supported by the National Nature Sci- ence Foundation of China under Grants 61503053, 61472374 and 61304197, the Natural Science Foundation Pro- ject of CQ CSTC, China (Grant No. cstc2013jcyjA40018), the Youth Science Research Project of CQUPT, China (Grant Nos. A2012-78 and A2012-82), the Doctor Start-up Foundation of CQUPT, China (Grant Nos. A2012-23 and A2012-26), the Natural Science Fundation of CQJW, China (Grant Nos. KJ130506 and KJ1400435), and Training Programme Foundation for the Talents of Higher Education Commission. This support is greatly appreciated.</p></sec><sec id="s7"><title>Cite this paper</title><p>FenglanSun,RuiWang,YongfuLi,FengLiu, (2016) Flocking for Leader-Following Multi-Agent Systems with Time-Varying Delay. Intelligent Control and Automation,07,9-15. doi: 10.4236/ica.2016.71002</p></sec></body><back><ref-list><title>References</title><ref id="scirp.63383-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Toner, J. and Tu, Y. (1998) Flocks, Herds, and Schools: A Quantitative Theory of Flocking. Physical Review E, 58, 4828. http://dx.doi.org/10.1103/PhysRevE.58.4828</mixed-citation></ref><ref id="scirp.63383-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Reynolds, C.W. (1987) Flocks, Herds and Schools: A Distributed Behavioral Model. ACM Siggraph Computer Graphics, 21, 25-34. http://dx.doi.org/10.1145/37402.37406</mixed-citation></ref><ref id="scirp.63383-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Toner, J. and Tu, Y. (2005) Hydrodynamics and Phases of Flocks. Annals of Physics, 318, 170-244. http://dx.doi.org/10.1016/j.aop.2005.04.011</mixed-citation></ref><ref id="scirp.63383-ref4"><label>4</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Olfati-Saber</surname><given-names> R. </given-names></name>,<etal>et al</etal>. (<year>2006</year>)<article-title>Flocking for Multi-Agents with a Virtual Leader</article-title><source> IEEE Transactions on Automatic Control</source><volume> 51</volume>,<fpage> 410</fpage>-<lpage>420</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.63383-ref5"><label>5</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Crowther</surname><given-names> B. </given-names></name>,<etal>et al</etal>. (<year>2002</year>)<article-title>Flocking of Autonomous Unmanned Air Vehicles</article-title><source> Aeronautical Journal</source><volume> 107</volume>,<fpage> 99</fpage>-<lpage>109</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.63383-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Akyildiz, I.F., Su, W. and Cayirci, E. (2002) A Survey on Sensor Networks. IEEE Communications Magazine, 40, 102-114. http://dx.doi.org/10.1109/MCOM.2002.1024422</mixed-citation></ref><ref id="scirp.63383-ref7"><label>7</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Partridge</surname><given-names> B.L. </given-names></name>,<etal>et al</etal>. (<year>1984</year>)<article-title>The Chorus-Line Hypothesis of Maneuver in Avian Flocks</article-title><source> Nature</source><volume> 309</volume>,<fpage> 344</fpage>-<lpage>345</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.63383-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Gazi, V. and Passino, K.M. (2004) Stability Analysis of Social Foraging Swarms. IEEE Transactions on Systems, Man, and Cybernetics, Part B: Cybernetics, 34, 539-557. http://dx.doi.org/10.1109/TSMCB.2003.817077</mixed-citation></ref><ref id="scirp.63383-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Ren, W. and Atkins, E. (2007) Distributed Multi-Vehicle Coordinated Control via Local Information Exchange. International Journal of Robust and Nonlinear Control, 17, 1002-1033. http://dx.doi.org/10.1002/rnc.1147</mixed-citation></ref><ref id="scirp.63383-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Tanner, H.G., Jadbabaie, A. and Pappas, G.J. (2007) Flocking in Fixed and Switching Networks. IEEE Transactions on Automatic Control, 52, 863-868. http://dx.doi.org/10.1109/TAC.2007.895948</mixed-citation></ref><ref id="scirp.63383-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Su, H., Wang, X. and Lin, Z. (2009) Flocking of Multi-Agents with a Virtual Leader. IEEE Transactions on Automatic Control, 54, 293-307. http://dx.doi.org/10.1109/TAC.2008.2010897</mixed-citation></ref><ref id="scirp.63383-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Topaz, C.M. and Bertozzi, A.L. (2004) Swarming Patterns in a Two-Dimensional Kinematic Model for Biological Groups. SIAM Journal on Applied Mathematics, 65, 152-174. http://dx.doi.org/10.1137/S0036139903437424</mixed-citation></ref><ref id="scirp.63383-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Chen, Z. and Zhang, H.T. (2013) A Remark on Collective Circular Motion of Heterogeneous Multi-Agents. Automatica, 49, 1236-1241. http://dx.doi.org/10.1016/j.automatica.2013.01.017</mixed-citation></ref><ref id="scirp.63383-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">Zhang, H.T., Zhai, C. and Chen, Z. (2011) A General Alignment Repulsion Algorithm for Flocking of Multi-Agent Systems. IEEE Transactions on Automatic Control, 56, 430-445. http://dx.doi.org/10.1109/TAC.2010.2089652</mixed-citation></ref><ref id="scirp.63383-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">Lu, X.Q., Austin, F. and Chen, S.H. (2011) Flocking in Multi-Agent Systems with Active Virtual Leader and Time-Varying Delays Coupling. Communications in Nonlinear Science and Numerical Simulation, 16, 1014-1026.http://dx.doi.org/10.1016/j.cnsns.2010.05.004</mixed-citation></ref><ref id="scirp.63383-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">Olfati-Saber, R. (2006) Flocking for Multi-Agent Dynamic Systems: Algorithms and Theory. IEEE Transactions on Automatic Control, 51, 401-420. http://dx.doi.org/10.1109/TAC.2005.864190</mixed-citation></ref><ref id="scirp.63383-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">Martin, S. (2014) Multi-Agent Flocking under Topological Interactions. Systems and Control Letters, 69, 53-61.http://dx.doi.org/10.1016/j.sysconle.2014.04.004</mixed-citation></ref><ref id="scirp.63383-ref18"><label>18</label><mixed-citation publication-type="other" xlink:type="simple">Fan, M.C., Zhang, H.T. and Wang, M. (2014) Bipartite Flocking for Multi-Agent Systems. Communications in Nonlinear Science and Numerical Simulation, 19, 3313-3322. http://dx.doi.org/10.1016/j.cnsns.2013.10.009</mixed-citation></ref><ref id="scirp.63383-ref19"><label>19</label><mixed-citation publication-type="other" xlink:type="simple">Gu, D. and Wang, Z. (2009) Leader-Follower Flocking: Algorithms and Experiments. IEEE Transactions on Control Systems Technology, 17, 1211-1219.</mixed-citation></ref><ref id="scirp.63383-ref20"><label>20</label><mixed-citation publication-type="other" xlink:type="simple">Jia, Y. and Zhang, W. (2014) Distributed Adaptive Flocking of Robotic Fish System with a Leader of Bounded Unknown Input. International Journal of Control, Automation and Systems, 12, 1049-1058.http://dx.doi.org/10.1007/s12555-013-0518-6</mixed-citation></ref><ref id="scirp.63383-ref21"><label>21</label><mixed-citation publication-type="other" xlink:type="simple">Li, M., Liang, J. and Li, S. (2014) Flocking Behavior with Multiple Leaders and Global Trajectory. Journal of Central South University, 21, 2324-2333. http://dx.doi.org/10.1007/s11771-014-2184-0</mixed-citation></ref><ref id="scirp.63383-ref22"><label>22</label><mixed-citation publication-type="other" xlink:type="simple">Yu, W., Chen, G. and Cao, M. (2010) Distributed Leader-Follower Flocking Control for Multi-Agent Dynamical Systems with Time-Varying Velocities. Systems &amp; Control Letters, 59, 543-552. http://dx.doi.org/10.1016/j.sysconle.2010.06.014</mixed-citation></ref><ref id="scirp.63383-ref23"><label>23</label><mixed-citation publication-type="other" xlink:type="simple">Luo, J., Zhang, S., Kang, H., et al. (2013) Flocking Algorithms for Multi-Agent Systems with Time-Delay. Proceedings of the IEEE International Workshop on Microwave and Millimeter Wave Circuits and System Technology (MMWCST), Chengdu, 24-25 October 2013, 428-431.</mixed-citation></ref><ref id="scirp.63383-ref24"><label>24</label><mixed-citation publication-type="other" xlink:type="simple">Shi, H., Wang, L. and Chu, T. (2006) Virtual Leader Approach to Coordinated Control of Multiple Mobile Agents with Asymmetric Interactions. Physica D: Nonlinear Phenomena, 213, 51-65. http://dx.doi.org/10.1016/j.physd.2005.10.012</mixed-citation></ref><ref id="scirp.63383-ref25"><label>25</label><mixed-citation publication-type="other" xlink:type="simple">Lu, X., Austin, F. and Chen, S. (2011) Flocking in Multi-Agent Systems with Active Virtual Leader and Time-Varying Delays Coupling. Communications in Nonlinear Science and Numerical Simulation, 16, 1014-1026.http://dx.doi.org/10.1016/j.cnsns.2010.05.004</mixed-citation></ref><ref id="scirp.63383-ref26"><label>26</label><mixed-citation publication-type="other" xlink:type="simple">Yang, Z., Zhang, Q. and Chen, Z. (2013) A Novel Adaptive Flocking Algorithm for Multi-Agents System with Time Delay and Nonlinear Dynamics. Proceedings of the 32nd IEEE Chinese Control Conference (CCC), Xi’an, 26-28 July 2013, 998-1001.</mixed-citation></ref></ref-list></back></article>