Network Synchronization With Nonlinear Dynamics and Switching Interactions

Size: px
Start display at page:

Download "Network Synchronization With Nonlinear Dynamics and Switching Interactions"

Transcription

1 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 61, NO. 10, OCTOBER Network Synchronization With Nonlinear Dynamics an Switching Interactions Tao Yang, Ziyang Meng, Guoong Shi, Yiguang Hong, an Karl Henrik Johansson Abstract This technical note consiers the synchronization problem for networks of couple nonlinear ynamical systems uner switching communication topologies. Two types of nonlinear agent ynamics are consiere. The first one is non-expansive ynamics [stable ynamics with a convex Lyapunov function ϕ )] an the secon one is ynamics that satisfies a global Lipschitz conition. For the non-expansive case, we show that various forms of joint connectivity for communication graphs are sufficient for networks to achieve global asymptotic ϕ-synchronization. We also show that ϕ-synchronization leas to state synchronization provie that certain aitional conitions are satisfie. For the globally Lipschitz case, unlike the non-expansive case, joint connectivity alone is not sufficient for achieving synchronization. A sufficient conition for reaching global exponential synchronization is establishe in terms of the relationship between the global Lipschitz constant an the network parameters. Inex Terms Multi-agent systems, nonlinear agents, switching interactions, synchronization. I. INTRODUCTION We consier the synchronization problem for a network of couple nonlinear agents with agent set V = {1, 2,...,N}. Their interactions communications in the network) are escribe by a time-varying irecte graph G σt) =V, E σt) ), with σ :[0, ) Pas a piecewise constant signal, where P is a finite set of all possible graphs over V. The state of agent i V at time t is enote as x i t) R n an evolves as ẋ i = ft, x i ) j N i σt a ij t)x j x i ), i V 1) where ft, x i ):[0, ) R n R n is piecewise continuous in t an continuous in x i representing the uncouple inherent agent ynamics, N i σt = {j V:j, i) E σt) } is the set of agent i s neighbors Manuscript receive August 23, 2015; revise October 23, 2015; accepte November 2, Date of publication November 4, 2015; ate of current version September 23, This work was supporte in part by the Knut an Alice Wallenberg Founation an the Sweish Research Council. Recommene by Associate Eitor S. Zampieri. T. Yang is with the Pacific Northwest National Laboratory, Richlan, WA USA Tao.Yang@pnnl.gov). Z. Meng is with the Department of Precision Instrument an the State Key Laboratory of Precision Measurement Technology an Instruments, Tsinghua University, Beijing , China ziyangmeng@mail. tsinghua.eu.cn). G. Shi is with the College of Engineering an Computer Science, The Australian National University, Canberra, ACT 0200, Australia guoong.shi@anu.eu.au). Y. Hong is with the Key Laboratory of Systems an Control, Institute of Systems Science, Chinese Acaemy of Science, Beijing , China yghong@iss.ac.cn). K. H. Johansson is with the ACCESS Linnaeus Centre, School of Electrical Engineering, Royal Institute of Technology, Stockholm, Sween kallej@kth.se). Digital Object Ientifier /TAC at time t, ana ij t) is a piecewise continuous function marking the weight of ege j, i) at time t. Note that N i σt may be an empty set at certain time intervals. Systems of the form 1) have attracte consierable attention. Most works focus on the case where the communication graph G σt) is fixe, e.g., [1] [7]. It is shown that for the case where ft, x i ) satisfies a Lipschitz conition, synchronization is achieve for a connecte graph provie that the coupling strength is sufficiently large. However, for the case where the communication graph is time-varying, the synchronization problem becomes much more challenging an existing literature mainly focuses on a few special cases when ft, x i ) is linear, e.g., the single-integrator case [8] [10], the ouble-integrator case [11], an the neutrally stable case [12], [13]. Other stuies assume some particular structures for the communication graph [14] [16]. In particular, in [14], the authors focus on the case where the ajacency matrices associate with all communication graphs are simultaneously triangularizable. The authors of [15] consier switching communication graphs that are weakly connecte an balance at all times. A more general case where the switching communication graph frequently has a irecte spanning tree has been consiere in [16]. These special structures on the switching communication graph are rather restrictive compare to joint connectivity where the communication can be lost at any time. This technical note aims to investigate whether joint connectivity for switching communication graphs can yiel synchronization for the nonlinear ynamics 1). We istinguish two cases epening on whether the nonlinear agent ynamics ft, x i ) is expansive or not. For the non-expansive case, we focus on the case where the nonlinear agent ynamics is stable with a convex Lyapunov function ϕ ). We show that various forms of joint connectivity for communication graphs are sufficient for networks to achieve global asymptotic ϕ- synchronization, that is, the function ϕ of the agent state converges to a common value. We also show that ϕ-synchronization implies state synchronization provie that aitional conitions are satisfie. For the expansive case, we focus on when the nonlinear agent ynamics is globally Lipschitz, an establish a sufficient conition for networks to achieve global exponential synchronization in terms of a relationship between the Lipschitz constant an the network parameters. The remainer of this technical note is organize as follows. Section II presents the problem efinition an main results. Section III provies all technical proofs. Finally, concluing remarks are offere in Section IV. A. Problem Set-Up II. PROBLEM DEFINITION AND MAIN RESULTS In this technical note, we consier the network of couple nonlinear agents 1). The communications are moele by a time-varying irecte graph G σt) =V, E σt) ). The joint graph of G σt) in the time interval [t 1,t 2 ) with t 1 <t 2 is enote as G[t 1,t 2 = t [t1,t 2 )Gt) =V, t [t1,t 2 )E σt) ). For the communication graph, we introuce the following efinition IEEE. Personal use is permitte, but republication/reistribution requires IEEE permission. See for more information.

2 3104 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 61, NO. 10, OCTOBER 2016 Definition 1: i) G σt) is uniformly jointly strongly connecte if there exists a constant T>0 such that G[t, t T is strongly connecte for any t 0. ii) Assume that G σt) is unirecte for all t 0. G σt) is infinitely jointly connecte if G[t, is connecte for any t 0. We assume that there are constants 0 <a a such that a a ij t) a for all t 0. We also make a stanar well time assumption [17] on the switching signal σt): there is a lower boun τ D > 0 between two consecutive switching time instants of σt). We enote x =[x T 1,xT 2,...,xT N ]T R nn an assume that the initial time is t = t 0 0, an the initial state xt 0 )=x T 1 t 0 ),...,x T N t 0 T R nn. In this technical note, we are intereste in the following synchronization problems. Definition 2: The multi-agent system 1) achieves global asymptotic ϕ-synchronization, where ϕ : R n R is a continuously ifferentiable function, if for any initial state xt 0 ), there exists a constant xt 0, such that lim t ϕx i t = xt 0 for any i V an any t 0 0. Definition 3: i) The multi-agent system 1) achieves global asymptotic synchronization if lim t x i t) x j t = 0 for any i, j V, any t 0 0 an any xt 0 ) R nn. Multi-agent system 1) achieves global exponential synchronization if there exist γ 1 an λ>0 such that x i t) x j t) 2 γe λt t 0) x i t 0 ) x j t 0 ) 2, t t 0 2) for any t 0 0 an any xt 0 ) R nn. Remark 1: ϕ-synchronization is a type of output synchronization where the output of agent i Vischosen to be ϕx i ). It is relate to butifferent fromχ-synchronization [18], [19] since ϕ is afunction of an iniviual agent state while χ is a function of all agent states. B. Non-Expansive Dynamics In this section, we focus on when the nonlinear agent ynamics is non-expansive as inicate by the following assumption. Assumption 1: ϕ : R n R is a continuously ifferentiable positive efinite convex function satisfying 1). lim η ϕη) = ; 2). ϕη),ft, η) 0 for any η R n an any t 0. The following lemma shows how Assumption 1 enforces nonexpansive ynamics. Lemma 1: Let Assumption 1 hol. Along the multi-agent ynamics 1), i V ϕx i t is non-increasing for all t 0. We now state main results for the non-expansive case. Theorem 1: Let Assumption 1 hol. The multi-agent system 1) achieves global asymptotic ϕ-synchronization if G σt) is uniformly jointly strongly connecte. Theorem 2: Let Assumption 1 hol. Assume that G σt) is unirecte for all t t 0. The multi-agent system 1) achieves global asymptotic ϕ-synchronization if G σt) is infinitely jointly connecte. Remark 2: For the linear time-varying case ft, x) =At)x, if there exists a matrix P = P T > 0 such that PAt)A T t)p 0, t 0 3) then ϕx) =x T Pxfor x R n satisfies Assumption 1. For the linear time-invariant case ft, x) =Ax, the conition 3) is equivalent to that the matrix A is neutrally stable [13]. C. ϕ-synchronization vs. State Synchronization The following result establishes conitions uner which ϕ- synchronization implies state synchronization. Theorem 3: Let G σt) Gwith G being a fixe, strongly connecte igraph uner which the multi-agent system 1) achieves global asymptotic ϕ-synchronization for some positive efinite function ϕ : R n R. Let Assumption 1 hol. Moreover, assume that 1) ft, η) is boune for any t 0 an any η R n. 2) c 1 η 2 ϕη) c 2 η 2 for some 0 <c 1 c 2 ; 3) ϕ ) is strongly convex. Then the multi-agent system 1) achieves global asymptotic synchronization. D. Lipschitz Dynamics We consier also the case when the nonlinear agent ynamics is possibly expansive. We focus on when the ynamics satisfies the following global Lipschitz conition. Assumption 2: There exists a constant L>0 such that ft, η) ft, ζ) L η ζ, η, ζ R n, t 0. 4) Our main result for this case is given below. Theorem 4: Let Assumption 2 hol. Assume that G σt) is uniformly jointly strongly connecte. The multi-agent system 1) achieves global exponential synchronization if L<ρ /2, where ρ is a constant epening on the network parameters. Remark 3: Assumption 2 an its variants have been consiere in the literature for fixe communication graphs, e.g., [1] [7]. Compare with the existing literature, we here stuy a more challenging case, where the communication graphs are time-varying. Unlike the fixe case where the global Lipschitz conition is sufficient to guarantee synchronization, Theorem 4 establishe a sufficient synchronization conition relate to the Lipschitz constant an the network parameters. III. PROOFS OF MAIN RESULTS In this section, we provie proofs of the main results. A. Proof of Lemma 1 Recall that the upper Dini erivative of a function ht) :a, b) R at t is efine as D ht) = lim sup s 0 ht s) ht/s).the following lemma from [10], [20] is useful for the proof. Lemma 2: Let V i t, x): R R n R i =1,...,N) be continuously ifferentiable an V t, x) = i=1,...,n V i t, x). IfIt) = {i {1, 2,...,N} : V t, xt = V i t, xt} is the set of inices where the imum is reache at t, then D V t, xt = i It) V i t, xt. Denote It) ={i V: j V ϕx j t = ϕx i t}. Wefirst note that the convexity property of ϕ ) implies that [21, pp. 69] ϕη),ζ η ϕζ) ϕη), η, ζ R n. 5) It then follows from Lemma 2, Assumption 1ii) an 5) that: D ϕx it = ϕx i ),ft, x i ) a ij t)x j x i ) i V i It) i It) j N i σt j N i σt a ij t)ϕx j ) ϕx i 0

3 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 61, NO. 10, OCTOBER where the last inequality follows from ϕx j ) ϕx i ). This proves the lemma. B. Proof of Theorem 1 It follows from Lemma 1 that for any initial state xt 0 ) R nn, there exists a constant = xt 0 0, such that lim t i V ϕx i )=. We shall show that is the require constant in Definition 2 of ϕ-synchronization. We first note that, by Lemma 1, for all i V, there exist constants 0 α i β i, such that lim inf ϕ x it = α i, lim sup ϕ x i t = β i. t t Also note that it follows from lim t i V ϕx i t = that for any ε>0, there exists T 1 ε) > 0 such that: ϕ x i t [0, ε], i V, t T 1 ε). 6) The proof of Theorem 1 is base on a contraiction argument an relies on the following lemma. Lemma 3: Let Assumption 1 hol. Assume that G σt) is uniformly jointly strongly connecte. If there exists an agent k 0 V such that 0 α k0 <, then there exists 0 < ρ <1 an t such that for all i V, ϕx i t N 1)T 0 ρm 0 1 ρ) ε),where T 0 T 2τ D 7) with T given in Definition 1i) an τ D is the well time. Proof: LetusfirstefineM 0 α k0 β k0 )/2) <.Then there exists an infinite time sequence t 0 < t 1 < < t k < with lim k t k = such that ϕx t k = M 0 for all k =1, 2,...Let then t k be greater than or equal to T 1 ε), an enote it by t k0. We now prove the lemma by estimating an upper boun of the scalar function ϕx i ) agent by agent. The proof is base on a generalization of the metho propose in the proof of [22, Lemma 4.3] but with substantial ifferences on the ynamics an the Lyapunov function. Moreover, the convexity of ϕ ) plays an important role. Step 1. Focus on agent k 0. By using Assumption 1ii), 5), an 6), we obtain that for all t t k0 t ϕx k 0 t = ϕx k0 ),ft, x k0 ) a k0 jt)x j x k0 ) j N k0 σt j N k0 σt a k0 jt)ϕx j ) ϕ x k0 a N 1) ε ϕ x k0. It then follows that for all t t k0 : ϕ x k0 t e λ 1t t k0 ) ϕ xk0 t k0 1 e λ 1t t k0 ) ) ε) 8) where λ 1 = a N 1). Step 2. Consier agent k 1 k 0 such that k 0,k 1 ) E σt) for t [ t k0, t k0 T 0 ). The existence of such an agent can be shown as follows. Since G σt) is uniformly jointly strongly connecte, it is not har to see that there exists an agent k 1 k 0 Van t 1 t k0 such that k 0,k 1 ) E σt) for t [t 1,t 1 τ D ) [ t k0, t k0 T 0 ). From 8), we obtain that for all t [ t k0, t k0 N 1)T 0 ] ϕ x k0 t κ 0 ρm 0 1 ρ) ε) 9) where ρ = e λ 1N 1)T 0 = e a N 1) 2 T 0. We can similarly boun ϕx k1 t by consiering two ifferent cases: i) ϕx k1 t >ϕx k0 t for all t [t 1,t 1 τ D ) an ii) there exists a time instant t 1 [t 1,t 1 τ D ) such that ϕx k1 t 1 ϕx k0 t 1 κ 0. For both cases, we obtain that for all t [t 1 τ D, t k0 N 1)T 0 ] ϕ x k1 t ϕ 1 M 0 1 ϕ 1 ) ε). where ϕ 1 =1 μ)ρ 2. It then follows from 9) an 0 <ϕ 1 <ρ<1 that for all t [t 1 τ D, t k0 N 1)T 0 ]: ϕ x j tϕ 1 M 0 1 ϕ 1 ) ε), j {k 0,k 1 }. 10) Step 3. Consier agent k 2 {k 0,k 1 } such that there exists an ege from the set {k 0,k 1 } to the agent k 2 in E σt) for t [t 2,t 2 τ D ) [ t k0 T 0, t k0 2T 0 ). The existence of such an agent k 2 an t 2 follows similarly from the argument in Step 2. Similarly, we can boun ϕx k2 t by consiering two ifferent cases an obtain that for all t [t 2 τ D, t k0 N 1)T 0 ] ϕ x k2 t ϕ 2 M 0 1 ϕ 2 ) ε) 11) where ϕ 2 = 1 μ)ρ 2 ) 2. By combining 10) an 11), an using 0 <ϕ 2 <ϕ 1 < 1, we obtain that for all t [t 2 τ D, t k0 N 1)T 0 ] ϕ x j t ϕ 2 M 0 1 ϕ 2 ) ε), j {k 0,k 1,k 2 }. Step 4. By repeating the above process on time intervals [ t k0 2T 0, t k0 3T 0 ),..., [ t k0 N 2)T 0, t k0 N 1)T 0 ),weeventually obtain that for all i V ϕ x i t k0 N 1)T 0 ϕn 1 M 0 1 ϕ N 1 ) ε). where ϕ N 1 = 1 μ)ρ 2 ) N 1. The result of the lemma then follows by choosing ρ = ϕ N 1 an t = t k0. We are now reay to prove Theorem 1 by contraiction. Suppose that there exists an agent k 0 V such that 0 α k0 <.Itthen follows from Lemma 3 that ϕx i t N 1)T 0 < for all i V, provie that ε< ρ M 0 )/1 ρ. This contraicts the fact that lim t i V ϕx i )=. Thus, there oes not exist an agent k 0 V such that 0 α k0 <. Hence, lim t ϕx i t = for all i V. C. Proof of Theorem 2 The proof relies on the following lemma. Lemma 4: Let Assumption 1 hol. Assume that G σt) is infinitely jointly connecte. If there exists an agent k 0 Vsuch that 0 α k0 <, then there exist 0 < ρ <1 an t such that ϕ x i t τ D ) ) ρm 0 1 ρ) ε), i V. The proof of Lemma 4 is similar to that of Lemma 3 an base on estimating an upper boun for the scalar quantity ϕx i ) agent by agent. However, the metho by which we obtaine the sequence of agents k 0,k 1,...,k N 1 in the proof of Theorem 1 oes not apply uner the conition that G σt) is infinitely jointly connecte. We can however apply the strategy in [22]. We omit the etails an refer to [23] for a complete proof. The remaining proof of Theorem 2 follows from a contraiction argument an Lemma 4 in the same way as the proof of Theorem 1.

4 3106 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 61, NO. 10, OCTOBER 2016 D. Proof of Theorem 3 If the multi-agent system 1) reaches asymptotic ϕ-synchronization, i.e., lim t ϕx i t = for all i V, then for any ɛ>0, there exists a T ɛ > 0 such that Step 3. We first note that the strong convexity of ϕ ) implies that [21, pp. 459] there exists an m>0 such that ϕη),ζ η ϕζ) ϕη) m 2 η ζ 2, η, ζ R n. 14) ɛ ϕ x i t ɛ, i V, t T ɛ. 12) If =0 the esire conclusion hols trivially, i.e., lim t x i t) =0for all i Vue to the positive efiniteness of ϕ ). Inthe remainer of the proof we assume > 0. We shall prove the result by contraiction. Suppose that state synchronization is not achieve. Then there exist two agents i 0,j 0 V such that lim sup t x i0 t) x j0 t) > 0. In other wors, there exist an infinite time sequence t 1 < <t k < with lim k t k =, an a constant >0 such that x i0 t k ) x j0 t k ) =2 N 1 /c 1 for all k 1, wherec 1 is given in conition ii) of Theorem 3. We ivie the following analysis into three steps. Step 1. In this step, we prove the following crucial claim. Claim.Foranyt k, there are two agents i,j Vwith i,j ) E such that x i t k ) x j t k ) /c 1. We establish this claim via a recursive analysis. If either i 0,j 0 ) E or j 0,i 0 ) E then the result follows trivially. Otherwise we pick up another agent k 0 satisfying that there is an ege between k 0 an {i 0,j 0 }.Thisk 0 always exists since G is strongly connecte. Then either x k0 t k ) x i0 t k ) 2 N 2 /c 1 or x k0 t k ) x j0 t k ) 2 N 2 /c 1 must hol. Thus, we have again either establishe the claim, or we can continue to select another agent ifferent from i 0, j 0,ank 0 an repeat the argument. Since we have a finite number of agents, the esire claim hols. Furthermore, since there is a finite number of agent pairs, without loss of generality, we assume that the given agent pair i,j oes not vary for ifferent t k otherwise we can always select an infinite subsequence of t k for the following iscussions). Step 2. In this step, we establish a lower boun of x i t) x j t) 2 for a small time interval after a particular t k satisfying t k > T ɛ. From 12) an lim x ϕx) = giveninassumption1i), we see that x i t) an ϕx i t) x j t are boune for all i, j V an for all t T ɛ. It then follows from conition i) of Theorem 3 an 1) that: t ϕ x i x j ) ϕ x i x j ),ft, x i ) f t, x j ) ϕ x i x j ), ϕ x i x j ), k 1 N i σt k 2 N j σt a i k 1 t)x k1 x i ) a j k 2 t)x k2 x j ) L for all t t k an some L > 0. Without loss of generality we assume that ɛ 1. ThenL will be inepenent of ɛ. By plugging in the fact that x i t k ) x j t k ) /c1 an using the conition ii) of Theorem 3, we obtain that x i t) x j t) 2 2c 2, t [ t k,t k 2L ]. 13) By using Assumption 1ii) an 14), we obtain for t T ɛ t ϕx j ta j i t) ϕ x i ) ϕx j ) m 2 a j i t) x i t) x j t) 2 k N j σt\{i } a j kt) ϕx k ) ϕ x j ). 15) By using 12), 13), 15), an conition ii) of Theorem 3, we obtain that for t [t k,t k /2L )] which yiels ϕ x j t k t ϕ x j t 2N 1)a ɛ a m 4c 2 2L It is then straightforwar to see that if we take ϕ x j t k [ ɛ 2N 1)a ɛ a ] m. 4c 2 2L 2L { a m 2 ɛ<min, 32c 2 L < a m 2 a m 32c 2 L 16N 1)a c 2 However, this contraicts the efinition of ϕ-synchronization since t k is arbitrarily chosen. This completes the proof an the esire conclusion hols. E. Proof of Theorem 4 }. The proof is base on the convergence analysis of the scalar quantity where V t, xt = V ij t, xt= 1 2 e 2Lt t 0) x i t) x j t) 2, V ij t, xt 16) {i, j} V V. 17) Unlike the contraiction argument use in the proof of Theorem 1, where the convergence rate is unclear, here we explicitly characterize the convergence rate. The proof relies on the following lemmas. Lemma 5: Let Assumption 4 hol. Along the multi-agent ynamics 1), V t, xt is non-increasing for all t 0.

5 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 61, NO. 10, OCTOBER Proof: Let V 1 V 2 be the set containing all the noe pairs that reach the imum at time t, i.e., V 1 t) V 2 t) ={{i, j} V V: V ij t) =V t)}. It is not har to obtain that D V = 1 2 {i,j} V 1 V 2 e 2Lt t 0) x i x j ) T ft, x i ) ft, x j e 2Lt t 0) x i x j ) T e 2Lt t 0) x i x j ) T k 1 N i σt k 2 N j σt Le 2Lt t0) x i x j 2 {i,j} V 1 V 2 e 2Lt t 0) k 1 N i σt {i,j} V 1 V 2 a ik1 t) x j x k1 2 x i x j 2) e 2Lt t 0) k 2 N j σt a ik1 t)x k1 x i ) a jk2 t)x k2 x j ) a jk2 t) x i x k2 2 x i x j 2) k 1 N i σt k 2 N j σt a ik1 t)v jk1 V ij ) a jk2 t)v ik2 V ij ) 0 18) where the first equality follows from Lemma 2 an 1), the first inequality follows from 4) an ±ab a 2 b 2 )/2) for all a, b R, an the secon inequality follows from 17). Lemma 6: Let Assumption 4 hol. Assume that G σt) is uniformly jointly strongly connecte. Then there exists 0 < β <1 such that V ij NT0,xNT 0 ) ) βv, {i, j} V V where N = N 1, T 0 given by 7) an V = V t 0,xt 0. Proof: The proof is base on the convergence analysis of V ij t, xt) for all agent pairs {i, j} V Vin several steps, which is similar to the proof of Lemma 3. Without loss of generality, we assume that t 0 =0. We also sometimes enote V t, xt an V ij t, xt as V an V ij, respectively, for notational simplification. Step 1. We begin by consiering any agent i 1 V.SinceG σt) is uniformly jointly strongly connecte, we know that i 1 is the root an that there exists a time t 1 ananagenti 2 V\{i 1 } such that i 1,i 2 ) Euring t [t 1,t 1 τ D ) [0,T 0 ]. We first note that it follows from 16) an Lemma 5 that for all t [0, NT 0 ] V ij t, xt V t, xt V, {i, j} V V. 19) Taking the erivative of V ij along the trajectories of 1), we obtain that for all t [t 1,t 1 τ D ) V i1 i 2 = Le 2Lt t0) x i1 x i2 2 e 2Lt t0) x i1 x i2 ) T a i1 k 1 t)x k1 x i1 ) k 1 N i1 σt a i2 k 2 t)x k2 x i2 )f t, x i1 ) ft, x i2 k 2 N i2 σt a i1 kt)v i2 k 1 V i1 i 2 ) a i2 i 1 t)v i1 i 2 k 1 N i1 σt k 2 N i2 σt\{i 1 } a i2 k 2 t)v i1 k 2 V i1 i 2 ) N 1)a V V i1 i 2 ) a V i1 i 2 N 2)a V V i1 i 2 ) ) 2N 3)a = α V i1i2 V α where α =2N 3)a a. The first inequality follows from 4) an 17), while the secon inequality follows from 19). It thus follows that: V i1 i 2 t 1 τ D,xt 1 τ D ˆα 1 V 20) where ˆα 1 =1 a /α)1 e ατ D) 0, 1). Similarly we obtain that for all t [t 1 τ D, NT 0 ], Vi1 i 2 αv V i1 i 2 ),whereα =2N 1)a. It then follows from 20) that: V i1 i 2 t, xt α 1V, t [t 1 τ D, NT 0 ] 21) where α 1 =1 1 e ατ D )a /α)e αnt 0 0, 1). Step 2. Since G σt) is uniformly jointly strongly connecte, we know that that there exists a time instant t 2 an an arc from h V 1 {i 1,i 2 } to i 3 V\V 1 uring [t 2,t 2 τ D ) [T 0, 2T 0 ]. We then estimate an upper boun for V hi3 by consiering two ifferent cases: h = i 1 an h = i 2. We eventually obtain that for all t [t 2 τ D, NT 0 ] where V hi3 t, xt 1 β 2 ) V, h V 1. 22) β =1 e ατ D ) a α e αnt 0 0, 1). 23) It then follows from 21), 22), an 1 β 2 > 1 β = α 1 that for all t [2T 0, NT 0 ]: V i1 k t, xt 1 β 2 ) V, k V 2 \{i 1 } where V 2 {i 1,i 2,i 3 }. Step 3. By continuing the above process, we obtain that for all k V\{i 1 } V i1 k NT0,xNT 0 ) ) ) 1 β N V. 24) Step 4. Since G σt) is uniformly jointly strongly connecte, 24) hols for any i 1 V. By using the same analysis, we eventually obtain that for all i, j V, V ij NT 0,xNT 0 1 β N )V. Hence the result follows by choosing β =1 β N with β given by 23). We are now reay to prove Theorem 4. By using Lemma 6 an 16), we obtain that V t, xt β t/nt 0 V 1/ β)e ρ t V,where

6 3108 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 61, NO. 10, OCTOBER 2016 t/nt 0 enotes the largest integer that is not greater than t/nt 0 ) an ρ =1/NT 0 )ln1/ β). It then follows from 16) an 17) that: x i t) x j t) 2 1 β e ρ 2L)t x i 0) x j 0) 2. Hence, global exponential synchronization is achieve with γ = 1/ β) an λ = ρ 2L provie that ρ > 2L. This conclues the proof of the esire theorem. IV. CONCLUSION In this technical note, synchronization problems for networks with nonlinear agent ynamics an switching topologies have been investigate. Two types of nonlinear ynamics were consiere: nonexpansive an globally Lipschitz. For the non-expansive case, we foun that the convexity of the Lyapunov function plays a crucial rule in the analysis an showe that the uniformly joint strong connectivity is sufficient for achieving global asymptotic ϕ-synchronization. When communication graphs are unirecte, the infinitely joint connectivity is a sufficient synchronization conition. Moreover, we establishe conitions uner which ϕ-synchronization implies state synchronization. For the globally Lipschitz case, we foun that joint connectivity alone is not sufficient to achieve synchronization but establishe a sufficient synchronization conition. The propose conition reveals the relationship between the Lipschitz constant an the network parameters. The results of this technical note can be extene to leaerfollower networks, see [23] for etails. An interesting future irection is to stuy the synchronization problem for couple non-ientical nonlinear ynamics uner general switching topologies. REFERENCES [1] C. W. Wu an L. O. Chua, Application of graph theory to the synchronization in an array of couple nonlinear oscillators, IEEE Trans. Circuits Syst. I, Funam. Theory Appl., vol. 42, no. 8, pp , Aug [2] C.W.Wu, Synchronization in Complex Networks of Nonlinear Dynamical Systems. Singapore: Worl Scientific: Singapore, [3] V. Belykh, I. Belykh, an M. Hasler, Connection graph stability metho for synchronize couple chaotic systems, Physica D, vol. 195, no. 1/2, pp , Aug [4] P. DeLellis, M. D. Bernaro, an G. Russo, On QUAD, Lipschitz, contracting vector fiels for consensus an synchronization of networks, IEEE Trans. Circuits Syst. I, Reg. Papers, vol. 58, no. 3, pp , Mar [5] W. Yu, G. Chen, an M. Cao, Consensus in irecte networks of agents with nonlinear ynamics, IEEE Trans. Autom. Control, vol. 56, no. 6, pp , Jun [6] U. Münz, A. Papachristooulou, an F. Allgöwer, Consensus in multiagent systems with coupling elays an switching topology, IEEE Trans. Autom. Control, vol. 56, no. 12, pp , Dec [7] H. Liu, M. Cao, an C. W. Wu, Coupling strength allocation for synchronization in complex networks using spectral graph theory, IEEE Trans. Circuits Syst. I, Reg. Papers, vol. 61, no. 5, pp , May [8] A. Jababaie, J. Lin, an A. S. Morse, Coorination of groups of mobile autonomous agents using nearest neighbor rules, IEEE Trans. Autom. Control, vol. 48, no. 6, pp , Jun [9] L. Moreau, Stability of multiagent systems with time-epenent communication links, IEEE Trans. Autom. Control, vol. 50, no. 2, pp , Feb [10] Z. Lin, B. Francis, an M. Maggiore, State agreement for continuoustime couple nonlinear systems, SIAM J. Control Optim., vol. 46, no. 1, pp , [11] W. Ren an R. W. Bear, Consensus seeking in multiagent systems uner ynamically changing interaction topologies, IEEE Trans. Autom. Control, vol. 50, no. 5, pp , May [12] L. Scarovi an R. Sepulchre, Synchronization in networks of ientical linear systems, Automatica, vol. 45, no. 11, pp , Nov [13] Y. Su an J. Huang, Stability of a class of linear switching systems with applications to two consensus problem, IEEE Trans. Autom. Control, vol. 57, no. 6, pp , Jun [14] J. Zhao, D. J. Hill, an T. Liu, Synchronization of complex ynamical networks with switching topology: A switche system point of view, Automatica, vol. 45, no. 11, pp , Nov [15] J. Qin, H. Gao, an W. Zheng, Exponential synchronization of complex networks of linear systems an nonlinear oscillators: A unifie analysis, IEEE Trans. Neural Netw. Learn. Syst., vol. 61, no. 2, pp , Mar [16] G. Wen, Z. Duan, G. Chen, an W. Yu, Consensus tracking of multi-agent systems with Lipschitz-type noe ynamics an switching topologies, IEEE Trans. Circuits Syst. I, Reg. Papers, vol. 61, no. 2, pp , Feb [17] D. Liberzon an A. S. Morse, Basic problem in stability an esign of switche systems, IEEE Control Syst. Mag., vol. 19, no. 5, pp , Oct [18] J. Cortés, Distribute algorithms for reaching consensus on general functions, Automatica, vol. 44, no. 3, pp , Mar [19] X. Wang an Y. Hong, Distribute finite-time χ-consensus algorithms for multi-agent systems with variable coupling topology, J. Syst. Sci. Complexity, vol. 23, no. 2, pp , [20] J. Danskin, The theory of -min, with applications, SIAM J. Appl. Math., vol. 14, no. 6, pp , [21] S. Boy an L. Vanenberghe, Convex Optimization. New York, NY, USA: Cambrige Univ. Press, [22] G. Shi, K. H. Johansson, an Y. Hong, Reaching an optimal consensus: Dynamical systems that compute intersections of convex sets, IEEE Trans. Autom. Control, vol. 58, no. 3, pp , Mar [23] T. Yang, Z. Meng, G. Shi, Y. Hong, an K. H. Johansson, Network synchronization with nonlinear ynamics an switching interactions, Systems an Control, pp. 1 27, [Online]. Available: abs/

arxiv: v3 [cs.sy] 24 Aug 2015

arxiv: v3 [cs.sy] 24 Aug 2015 Network Synchronization with Nonlinear Dynamics and Switching Interactions Tao Yang, Ziyang Meng, Guodong Shi, Yiguang Hong and Karl Henrik Johansson arxiv:1401.6541v3 [cs.sy] 24 Aug 2015 Abstract This

More information

Multi-agent Systems Reaching Optimal Consensus with Time-varying Communication Graphs

Multi-agent Systems Reaching Optimal Consensus with Time-varying Communication Graphs Preprints of the 8th IFAC Worl Congress Multi-agent Systems Reaching Optimal Consensus with Time-varying Communication Graphs Guoong Shi ACCESS Linnaeus Centre, School of Electrical Engineering, Royal

More information

Laplacian Cooperative Attitude Control of Multiple Rigid Bodies

Laplacian Cooperative Attitude Control of Multiple Rigid Bodies Laplacian Cooperative Attitue Control of Multiple Rigi Boies Dimos V. Dimarogonas, Panagiotis Tsiotras an Kostas J. Kyriakopoulos Abstract Motivate by the fact that linear controllers can stabilize the

More information

IN the multiagent systems literature, the consensus problem,

IN the multiagent systems literature, the consensus problem, IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS II: EXPRESS BRIEFS, VOL. 63, NO. 7, JULY 206 663 Periodic Behaviors for Discrete-Time Second-Order Multiagent Systems With Input Saturation Constraints Tao Yang,

More information

Exponential Tracking Control of Nonlinear Systems with Actuator Nonlinearity

Exponential Tracking Control of Nonlinear Systems with Actuator Nonlinearity Preprints of the 9th Worl Congress The International Feeration of Automatic Control Cape Town, South Africa. August -9, Exponential Tracking Control of Nonlinear Systems with Actuator Nonlinearity Zhengqiang

More information

Lie symmetry and Mei conservation law of continuum system

Lie symmetry and Mei conservation law of continuum system Chin. Phys. B Vol. 20, No. 2 20 020 Lie symmetry an Mei conservation law of continuum system Shi Shen-Yang an Fu Jing-Li Department of Physics, Zhejiang Sci-Tech University, Hangzhou 3008, China Receive

More information

Continuous observer design for nonlinear systems with sampled and delayed output measurements

Continuous observer design for nonlinear systems with sampled and delayed output measurements Preprints of th9th Worl Congress The International Feeration of Automatic Control Continuous observer esign for nonlinear systems with sample an elaye output measurements Daoyuan Zhang Yanjun Shen Xiaohua

More information

A Note on Exact Solutions to Linear Differential Equations by the Matrix Exponential

A Note on Exact Solutions to Linear Differential Equations by the Matrix Exponential Avances in Applie Mathematics an Mechanics Av. Appl. Math. Mech. Vol. 1 No. 4 pp. 573-580 DOI: 10.4208/aamm.09-m0946 August 2009 A Note on Exact Solutions to Linear Differential Equations by the Matrix

More information

Periodic Behaviors in Multi-agent Systems with Input Saturation Constraints

Periodic Behaviors in Multi-agent Systems with Input Saturation Constraints 5nd IEEE Conference on Decision and Control December 10-13, 013. Florence, Italy Periodic Behaviors in Multi-agent Systems with Input Saturation Constraints Tao Yang, Ziyang Meng, Dimos V. Dimarogonas,

More information

Reachable Set Analysis for Dynamic Neural Networks with Polytopic Uncertainties

Reachable Set Analysis for Dynamic Neural Networks with Polytopic Uncertainties Commun. Theor. Phys. 57 (2012) 904 908 Vol. 57, No. 5, May 15, 2012 Reachable Set Analysis for Dynamic Neural Networks with Polytopic Uncertainties ZUO Zhi-Qiang ( ãö), CHEN Yin-Ping (í ), an WANG Yi-Jing

More information

Generalized-Type Synchronization of Hyperchaotic Oscillators Using a Vector Signal

Generalized-Type Synchronization of Hyperchaotic Oscillators Using a Vector Signal Commun. Theor. Phys. (Beijing, China) 44 (25) pp. 72 78 c International Acaemic Publishers Vol. 44, No. 1, July 15, 25 Generalize-Type Synchronization of Hyperchaotic Oscillators Using a Vector Signal

More information

Introduction to the Vlasov-Poisson system

Introduction to the Vlasov-Poisson system Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its

More information

Adaptive Gain-Scheduled H Control of Linear Parameter-Varying Systems with Time-Delayed Elements

Adaptive Gain-Scheduled H Control of Linear Parameter-Varying Systems with Time-Delayed Elements Aaptive Gain-Scheule H Control of Linear Parameter-Varying Systems with ime-delaye Elements Yoshihiko Miyasato he Institute of Statistical Mathematics 4-6-7 Minami-Azabu, Minato-ku, okyo 6-8569, Japan

More information

Lower bounds on Locality Sensitive Hashing

Lower bounds on Locality Sensitive Hashing Lower bouns on Locality Sensitive Hashing Rajeev Motwani Assaf Naor Rina Panigrahy Abstract Given a metric space (X, X ), c 1, r > 0, an p, q [0, 1], a istribution over mappings H : X N is calle a (r,

More information

Switching Time Optimization in Discretized Hybrid Dynamical Systems

Switching Time Optimization in Discretized Hybrid Dynamical Systems Switching Time Optimization in Discretize Hybri Dynamical Systems Kathrin Flaßkamp, To Murphey, an Sina Ober-Blöbaum Abstract Switching time optimization (STO) arises in systems that have a finite set

More information

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x) Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)

More information

Stability Analysis of Nonlinear Networks via M-matrix Theory: Beyond Linear Consensus

Stability Analysis of Nonlinear Networks via M-matrix Theory: Beyond Linear Consensus 2012 American Control Conference Fairmont Queen Elizabeth, Montréal, Canaa June 27-June 29, 2012 Stability Analysis of Nonlinear Networks via M-matrix Theory: Beyon Linear Consensus Airlie Chapman an Mehran

More information

6 General properties of an autonomous system of two first order ODE

6 General properties of an autonomous system of two first order ODE 6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x

More information

Math 342 Partial Differential Equations «Viktor Grigoryan

Math 342 Partial Differential Equations «Viktor Grigoryan Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite

More information

Lecture 2 Lagrangian formulation of classical mechanics Mechanics

Lecture 2 Lagrangian formulation of classical mechanics Mechanics Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,

More information

Lecture 10: October 30, 2017

Lecture 10: October 30, 2017 Information an Coing Theory Autumn 2017 Lecturer: Mahur Tulsiani Lecture 10: October 30, 2017 1 I-Projections an applications In this lecture, we will talk more about fining the istribution in a set Π

More information

Convergence of Random Walks

Convergence of Random Walks Chapter 16 Convergence of Ranom Walks This lecture examines the convergence of ranom walks to the Wiener process. This is very important both physically an statistically, an illustrates the utility of

More information

MULTI-AGENT TRACKING OF A HIGH-DIMENSIONAL ACTIVE LEADER WITH SWITCHING TOPOLOGY

MULTI-AGENT TRACKING OF A HIGH-DIMENSIONAL ACTIVE LEADER WITH SWITCHING TOPOLOGY Jrl Syst Sci & Complexity (2009) 22: 722 731 MULTI-AGENT TRACKING OF A HIGH-DIMENSIONAL ACTIVE LEADER WITH SWITCHING TOPOLOGY Yiguang HONG Xiaoli WANG Received: 11 May 2009 / Revised: 16 June 2009 c 2009

More information

Stability Analysis for Switched Systems with Sequence-based Average Dwell Time

Stability Analysis for Switched Systems with Sequence-based Average Dwell Time 1 Stability Analysis for Switched Systems with Sequence-based Average Dwell Time Dianhao Zheng, Hongbin Zhang, Senior Member, IEEE, J. Andrew Zhang, Senior Member, IEEE, Steven W. Su, Senior Member, IEEE

More information

TRAJECTORY TRACKING FOR FULLY ACTUATED MECHANICAL SYSTEMS

TRAJECTORY TRACKING FOR FULLY ACTUATED MECHANICAL SYSTEMS TRAJECTORY TRACKING FOR FULLY ACTUATED MECHANICAL SYSTEMS Francesco Bullo Richar M. Murray Control an Dynamical Systems California Institute of Technology Pasaena, CA 91125 Fax : + 1-818-796-8914 email

More information

ELEC3114 Control Systems 1

ELEC3114 Control Systems 1 ELEC34 Control Systems Linear Systems - Moelling - Some Issues Session 2, 2007 Introuction Linear systems may be represente in a number of ifferent ways. Figure shows the relationship between various representations.

More information

Lecture Introduction. 2 Examples of Measure Concentration. 3 The Johnson-Lindenstrauss Lemma. CS-621 Theory Gems November 28, 2012

Lecture Introduction. 2 Examples of Measure Concentration. 3 The Johnson-Lindenstrauss Lemma. CS-621 Theory Gems November 28, 2012 CS-6 Theory Gems November 8, 0 Lecture Lecturer: Alesaner Mąry Scribes: Alhussein Fawzi, Dorina Thanou Introuction Toay, we will briefly iscuss an important technique in probability theory measure concentration

More information

Tensors, Fields Pt. 1 and the Lie Bracket Pt. 1

Tensors, Fields Pt. 1 and the Lie Bracket Pt. 1 Tensors, Fiels Pt. 1 an the Lie Bracket Pt. 1 PHYS 500 - Southern Illinois University September 8, 2016 PHYS 500 - Southern Illinois University Tensors, Fiels Pt. 1 an the Lie Bracket Pt. 1 September 8,

More information

Optimal Control of Spatially Distributed Systems

Optimal Control of Spatially Distributed Systems Optimal Control of Spatially Distribute Systems Naer Motee an Ali Jababaie Abstract In this paper, we stuy the structural properties of optimal control of spatially istribute systems. Such systems consist

More information

. ISSN (print), (online) International Journal of Nonlinear Science Vol.6(2008) No.3,pp

. ISSN (print), (online) International Journal of Nonlinear Science Vol.6(2008) No.3,pp . ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.6(8) No.3,pp.195-1 A Bouneness Criterion for Fourth Orer Nonlinear Orinary Differential Equations with Delay

More information

Slide10 Haykin Chapter 14: Neurodynamics (3rd Ed. Chapter 13)

Slide10 Haykin Chapter 14: Neurodynamics (3rd Ed. Chapter 13) Slie10 Haykin Chapter 14: Neuroynamics (3r E. Chapter 13) CPSC 636-600 Instructor: Yoonsuck Choe Spring 2012 Neural Networks with Temporal Behavior Inclusion of feeback gives temporal characteristics to

More information

1 Heisenberg Representation

1 Heisenberg Representation 1 Heisenberg Representation What we have been ealing with so far is calle the Schröinger representation. In this representation, operators are constants an all the time epenence is carrie by the states.

More information

SYSTEMS OF DIFFERENTIAL EQUATIONS, EULER S FORMULA. where L is some constant, usually called the Lipschitz constant. An example is

SYSTEMS OF DIFFERENTIAL EQUATIONS, EULER S FORMULA. where L is some constant, usually called the Lipschitz constant. An example is SYSTEMS OF DIFFERENTIAL EQUATIONS, EULER S FORMULA. Uniqueness for solutions of ifferential equations. We consier the system of ifferential equations given by x = v( x), () t with a given initial conition

More information

Power Generation and Distribution via Distributed Coordination Control

Power Generation and Distribution via Distributed Coordination Control Power Generation an Distribution via Distribute Coorination Control Byeong-Yeon Kim, Kwang-Kyo Oh, an Hyo-Sung Ahn arxiv:407.4870v [math.oc] 8 Jul 204 Abstract This paper presents power coorination, power

More information

Criteria for Global Stability of Coupled Systems with Application to Robust Output Feedback Design for Active Surge Control

Criteria for Global Stability of Coupled Systems with Application to Robust Output Feedback Design for Active Surge Control Criteria for Global Stability of Couple Systems with Application to Robust Output Feeback Design for Active Surge Control Shiriaev, Anton; Johansson, Rolf; Robertsson, Aners; Freiovich, Leoni 9 Link to

More information

Dissipative numerical methods for the Hunter-Saxton equation

Dissipative numerical methods for the Hunter-Saxton equation Dissipative numerical methos for the Hunter-Saton equation Yan Xu an Chi-Wang Shu Abstract In this paper, we present further evelopment of the local iscontinuous Galerkin (LDG) metho esigne in [] an a

More information

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent

More information

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5

More information

Exponential asymptotic property of a parallel repairable system with warm standby under common-cause failure

Exponential asymptotic property of a parallel repairable system with warm standby under common-cause failure J. Math. Anal. Appl. 341 (28) 457 466 www.elsevier.com/locate/jmaa Exponential asymptotic property of a parallel repairable system with warm stanby uner common-cause failure Zifei Shen, Xiaoxiao Hu, Weifeng

More information

Proof of SPNs as Mixture of Trees

Proof of SPNs as Mixture of Trees A Proof of SPNs as Mixture of Trees Theorem 1. If T is an inuce SPN from a complete an ecomposable SPN S, then T is a tree that is complete an ecomposable. Proof. Argue by contraiction that T is not a

More information

Nonlinear Adaptive Ship Course Tracking Control Based on Backstepping and Nussbaum Gain

Nonlinear Adaptive Ship Course Tracking Control Based on Backstepping and Nussbaum Gain Nonlinear Aaptive Ship Course Tracking Control Base on Backstepping an Nussbaum Gain Jialu Du, Chen Guo Abstract A nonlinear aaptive controller combining aaptive Backstepping algorithm with Nussbaum gain

More information

Permanent vs. Determinant

Permanent vs. Determinant Permanent vs. Determinant Frank Ban Introuction A major problem in theoretical computer science is the Permanent vs. Determinant problem. It asks: given an n by n matrix of ineterminates A = (a i,j ) an

More information

An Optimal Algorithm for Bandit and Zero-Order Convex Optimization with Two-Point Feedback

An Optimal Algorithm for Bandit and Zero-Order Convex Optimization with Two-Point Feedback Journal of Machine Learning Research 8 07) - Submitte /6; Publishe 5/7 An Optimal Algorithm for Banit an Zero-Orer Convex Optimization with wo-point Feeback Oha Shamir Department of Computer Science an

More information

TOEPLITZ AND POSITIVE SEMIDEFINITE COMPLETION PROBLEM FOR CYCLE GRAPH

TOEPLITZ AND POSITIVE SEMIDEFINITE COMPLETION PROBLEM FOR CYCLE GRAPH English NUMERICAL MATHEMATICS Vol14, No1 Series A Journal of Chinese Universities Feb 2005 TOEPLITZ AND POSITIVE SEMIDEFINITE COMPLETION PROBLEM FOR CYCLE GRAPH He Ming( Λ) Michael K Ng(Ξ ) Abstract We

More information

E-Companion to The Evolution of Beliefs over Signed Social Networks

E-Companion to The Evolution of Beliefs over Signed Social Networks OPERATIONS RESEARCH INFORMS E-Companion to The Evolution of Beliefs over Signed Social Networks Guodong Shi Research School of Engineering, CECS, The Australian National University, Canberra ACT 000, Australia

More information

A new four-dimensional chaotic system

A new four-dimensional chaotic system Chin. Phys. B Vol. 19 No. 12 2010) 120510 A new four-imensional chaotic system Chen Yong ) a)b) an Yang Yun-Qing ) a) a) Shanghai Key Laboratory of Trustworthy Computing East China Normal University Shanghai

More information

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013 Survey Sampling Kosuke Imai Department of Politics, Princeton University February 19, 2013 Survey sampling is one of the most commonly use ata collection methos for social scientists. We begin by escribing

More information

A new proof of the sharpness of the phase transition for Bernoulli percolation on Z d

A new proof of the sharpness of the phase transition for Bernoulli percolation on Z d A new proof of the sharpness of the phase transition for Bernoulli percolation on Z Hugo Duminil-Copin an Vincent Tassion October 8, 205 Abstract We provie a new proof of the sharpness of the phase transition

More information

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments 2 Conference on Information Sciences an Systems, The Johns Hopkins University, March 2, 2 Time-of-Arrival Estimation in Non-Line-Of-Sight Environments Sinan Gezici, Hisashi Kobayashi an H. Vincent Poor

More information

GLOBAL SOLUTIONS FOR 2D COUPLED BURGERS-COMPLEX-GINZBURG-LANDAU EQUATIONS

GLOBAL SOLUTIONS FOR 2D COUPLED BURGERS-COMPLEX-GINZBURG-LANDAU EQUATIONS Electronic Journal of Differential Equations, Vol. 015 015), No. 99, pp. 1 14. ISSN: 107-6691. URL: http://eje.math.txstate.eu or http://eje.math.unt.eu ftp eje.math.txstate.eu GLOBAL SOLUTIONS FOR D COUPLED

More information

TRACKING CONTROL OF MULTIPLE MOBILE ROBOTS: A CASE STUDY OF INTER-ROBOT COLLISION-FREE PROBLEM

TRACKING CONTROL OF MULTIPLE MOBILE ROBOTS: A CASE STUDY OF INTER-ROBOT COLLISION-FREE PROBLEM 265 Asian Journal of Control, Vol. 4, No. 3, pp. 265-273, September 22 TRACKING CONTROL OF MULTIPLE MOBILE ROBOTS: A CASE STUDY OF INTER-ROBOT COLLISION-FREE PROBLEM Jurachart Jongusuk an Tsutomu Mita

More information

Consensus Tracking for Multi-Agent Systems with Nonlinear Dynamics under Fixed Communication Topologies

Consensus Tracking for Multi-Agent Systems with Nonlinear Dynamics under Fixed Communication Topologies Proceedings of the World Congress on Engineering and Computer Science Vol I WCECS, October 9-,, San Francisco, USA Consensus Tracking for Multi-Agent Systems with Nonlinear Dynamics under Fixed Communication

More information

Complex Laplacians and Applications in Multi-Agent Systems

Complex Laplacians and Applications in Multi-Agent Systems 1 Complex Laplacians and Applications in Multi-Agent Systems Jiu-Gang Dong, and Li Qiu, Fellow, IEEE arxiv:1406.186v [math.oc] 14 Apr 015 Abstract Complex-valued Laplacians have been shown to be powerful

More information

On the Cauchy Problem for Von Neumann-Landau Wave Equation

On the Cauchy Problem for Von Neumann-Landau Wave Equation Journal of Applie Mathematics an Physics 4 4-3 Publishe Online December 4 in SciRes http://wwwscirporg/journal/jamp http://xoiorg/436/jamp4343 On the Cauchy Problem for Von Neumann-anau Wave Equation Chuangye

More information

A Weak First Digit Law for a Class of Sequences

A Weak First Digit Law for a Class of Sequences International Mathematical Forum, Vol. 11, 2016, no. 15, 67-702 HIKARI Lt, www.m-hikari.com http://x.oi.org/10.1288/imf.2016.6562 A Weak First Digit Law for a Class of Sequences M. A. Nyblom School of

More information

Agmon Kolmogorov Inequalities on l 2 (Z d )

Agmon Kolmogorov Inequalities on l 2 (Z d ) Journal of Mathematics Research; Vol. 6, No. ; 04 ISSN 96-9795 E-ISSN 96-9809 Publishe by Canaian Center of Science an Eucation Agmon Kolmogorov Inequalities on l (Z ) Arman Sahovic Mathematics Department,

More information

WELL-POSEDNESS OF A POROUS MEDIUM FLOW WITH FRACTIONAL PRESSURE IN SOBOLEV SPACES

WELL-POSEDNESS OF A POROUS MEDIUM FLOW WITH FRACTIONAL PRESSURE IN SOBOLEV SPACES Electronic Journal of Differential Equations, Vol. 017 (017), No. 38, pp. 1 7. ISSN: 107-6691. URL: http://eje.math.txstate.eu or http://eje.math.unt.eu WELL-POSEDNESS OF A POROUS MEDIUM FLOW WITH FRACTIONAL

More information

A Sketch of Menshikov s Theorem

A Sketch of Menshikov s Theorem A Sketch of Menshikov s Theorem Thomas Bao March 14, 2010 Abstract Let Λ be an infinite, locally finite oriente multi-graph with C Λ finite an strongly connecte, an let p

More information

Robust Forward Algorithms via PAC-Bayes and Laplace Distributions. ω Q. Pr (y(ω x) < 0) = Pr A k

Robust Forward Algorithms via PAC-Bayes and Laplace Distributions. ω Q. Pr (y(ω x) < 0) = Pr A k A Proof of Lemma 2 B Proof of Lemma 3 Proof: Since the support of LL istributions is R, two such istributions are equivalent absolutely continuous with respect to each other an the ivergence is well-efine

More information

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control 19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior

More information

MAINTAINING LIMITED-RANGE CONNECTIVITY AMONG SECOND-ORDER AGENTS

MAINTAINING LIMITED-RANGE CONNECTIVITY AMONG SECOND-ORDER AGENTS MAINTAINING LIMITED-RANGE CONNECTIVITY AMONG SECOND-ORDER AGENTS KETAN SAVLA, GIUSEPPE NOTARSTEFANO AND FRANCESCO BULLO Abstract. In this paper we consier a-hoc networks of robotic agents with ouble integrator

More information

The total derivative. Chapter Lagrangian and Eulerian approaches

The total derivative. Chapter Lagrangian and Eulerian approaches Chapter 5 The total erivative 51 Lagrangian an Eulerian approaches The representation of a flui through scalar or vector fiels means that each physical quantity uner consieration is escribe as a function

More information

Distributed Force/Position Consensus Tracking of Networked Robotic Manipulators

Distributed Force/Position Consensus Tracking of Networked Robotic Manipulators 180 IEEE/CAA JOURNAL OF AUTOMATICA SINICA, VOL. 1, NO. 2, APRIL 2014 Distribute Force/Position Consensus Tracking of Networke Robotic Manipulators Lijiao Wang Bin Meng Abstract In this paper, we aress

More information

From Local to Global Control

From Local to Global Control Proceeings of the 47th IEEE Conference on Decision an Control Cancun, Mexico, Dec. 9-, 8 ThB. From Local to Global Control Stephen P. Banks, M. Tomás-Roríguez. Automatic Control Engineering Department,

More information

On colour-blind distinguishing colour pallets in regular graphs

On colour-blind distinguishing colour pallets in regular graphs J Comb Optim (2014 28:348 357 DOI 10.1007/s10878-012-9556-x On colour-blin istinguishing colour pallets in regular graphs Jakub Przybyło Publishe online: 25 October 2012 The Author(s 2012. This article

More information

Optimal Control of Spatially Distributed Systems

Optimal Control of Spatially Distributed Systems Optimal Control of Spatially Distribute Systems Naer Motee an Ali Jababaie Abstract In this paper, we stuy the structural properties of optimal control of spatially istribute systems. Such systems consist

More information

Time Headway Requirements for String Stability of Homogeneous Linear Unidirectionally Connected Systems

Time Headway Requirements for String Stability of Homogeneous Linear Unidirectionally Connected Systems Joint 48th IEEE Conference on Decision an Control an 8th Chinese Control Conference Shanghai, PR China, December 6-8, 009 WeBIn53 Time Heaway Requirements for String Stability of Homogeneous Linear Uniirectionally

More information

The Wiener Index of Trees with Prescribed Diameter

The Wiener Index of Trees with Prescribed Diameter 011 1 15 4 ± Dec., 011 Operations Research Transactions Vol.15 No.4 The Wiener Inex of Trees with Prescribe Diameter XING Baohua 1 CAI Gaixiang 1 Abstract The Wiener inex W(G) of a graph G is efine as

More information

Table of Common Derivatives By David Abraham

Table of Common Derivatives By David Abraham Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec

More information

A nonlinear inverse problem of the Korteweg-de Vries equation

A nonlinear inverse problem of the Korteweg-de Vries equation Bull. Math. Sci. https://oi.org/0.007/s3373-08-025- A nonlinear inverse problem of the Korteweg-e Vries equation Shengqi Lu Miaochao Chen 2 Qilin Liu 3 Receive: 0 March 207 / Revise: 30 April 208 / Accepte:

More information

Distributed Robust Consensus of Heterogeneous Uncertain Multi-agent Systems

Distributed Robust Consensus of Heterogeneous Uncertain Multi-agent Systems Distributed Robust Consensus of Heterogeneous Uncertain Multi-agent Systems Zhongkui Li, Zhisheng Duan, Frank L. Lewis. State Key Laboratory for Turbulence and Complex Systems, Department of Mechanics

More information

Minimum-time constrained velocity planning

Minimum-time constrained velocity planning 7th Meiterranean Conference on Control & Automation Makeonia Palace, Thessaloniki, Greece June 4-6, 9 Minimum-time constraine velocity planning Gabriele Lini, Luca Consolini, Aurelio Piazzi Università

More information

arxiv: v2 [math.dg] 16 Dec 2014

arxiv: v2 [math.dg] 16 Dec 2014 A ONOTONICITY FORULA AND TYPE-II SINGULARITIES FOR THE EAN CURVATURE FLOW arxiv:1312.4775v2 [math.dg] 16 Dec 2014 YONGBING ZHANG Abstract. In this paper, we introuce a monotonicity formula for the mean

More information

SINGULAR PERTURBATION AND STATIONARY SOLUTIONS OF PARABOLIC EQUATIONS IN GAUSS-SOBOLEV SPACES

SINGULAR PERTURBATION AND STATIONARY SOLUTIONS OF PARABOLIC EQUATIONS IN GAUSS-SOBOLEV SPACES Communications on Stochastic Analysis Vol. 2, No. 2 (28) 289-36 Serials Publications www.serialspublications.com SINGULAR PERTURBATION AND STATIONARY SOLUTIONS OF PARABOLIC EQUATIONS IN GAUSS-SOBOLEV SPACES

More information

BEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS. Mauro Boccadoro Magnus Egerstedt Paolo Valigi Yorai Wardi

BEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS. Mauro Boccadoro Magnus Egerstedt Paolo Valigi Yorai Wardi BEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS Mauro Boccaoro Magnus Egerstet Paolo Valigi Yorai Wari {boccaoro,valigi}@iei.unipg.it Dipartimento i Ingegneria Elettronica

More information

On Dwell Time Minimization for Switched Delay Systems: Free-Weighting Matrices Method

On Dwell Time Minimization for Switched Delay Systems: Free-Weighting Matrices Method On Dwell Time Minimization for Switched Delay Systems: Free-Weighting Matrices Method Ahmet Taha Koru Akın Delibaşı and Hitay Özbay Abstract In this paper we present a quasi-convex minimization method

More information

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+

More information

Ramsey numbers of some bipartite graphs versus complete graphs

Ramsey numbers of some bipartite graphs versus complete graphs Ramsey numbers of some bipartite graphs versus complete graphs Tao Jiang, Michael Salerno Miami University, Oxfor, OH 45056, USA Abstract. The Ramsey number r(h, K n ) is the smallest positive integer

More information

On the number of isolated eigenvalues of a pair of particles in a quantum wire

On the number of isolated eigenvalues of a pair of particles in a quantum wire On the number of isolate eigenvalues of a pair of particles in a quantum wire arxiv:1812.11804v1 [math-ph] 31 Dec 2018 Joachim Kerner 1 Department of Mathematics an Computer Science FernUniversität in

More information

Monotonicity for excited random walk in high dimensions

Monotonicity for excited random walk in high dimensions Monotonicity for excite ranom walk in high imensions Remco van er Hofsta Mark Holmes March, 2009 Abstract We prove that the rift θ, β) for excite ranom walk in imension is monotone in the excitement parameter

More information

Exponential Energy Decay of Solutions for a Transmission Problem With Viscoelastic Term and Delay

Exponential Energy Decay of Solutions for a Transmission Problem With Viscoelastic Term and Delay mathematics Article Exponential Energy Decay of Solutions for a Transmission Problem With Viscoelastic Term an Delay Danhua Wang *, Gang Li an Biqing Zhu College of Mathematics an Statistics, Nanjing University

More information

Iterated Point-Line Configurations Grow Doubly-Exponentially

Iterated Point-Line Configurations Grow Doubly-Exponentially Iterate Point-Line Configurations Grow Doubly-Exponentially Joshua Cooper an Mark Walters July 9, 008 Abstract Begin with a set of four points in the real plane in general position. A to this collection

More information

NOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy,

NOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy, NOTES ON EULER-BOOLE SUMMATION JONATHAN M BORWEIN, NEIL J CALKIN, AND DANTE MANNA Abstract We stuy a connection between Euler-MacLaurin Summation an Boole Summation suggeste in an AMM note from 196, which

More information

Tracking control for multi-agent consensus with an active leader and variable topology

Tracking control for multi-agent consensus with an active leader and variable topology Automatica 42 (2006) 1177 1182 wwwelseviercom/locate/automatica Brief paper Tracking control for multi-agent consensus with an active leader and variable topology Yiguang Hong a,, Jiangping Hu a, Linxin

More information

Quantum mechanical approaches to the virial

Quantum mechanical approaches to the virial Quantum mechanical approaches to the virial S.LeBohec Department of Physics an Astronomy, University of Utah, Salt Lae City, UT 84112, USA Date: June 30 th 2015 In this note, we approach the virial from

More information

A note on asymptotic formulae for one-dimensional network flow problems Carlos F. Daganzo and Karen R. Smilowitz

A note on asymptotic formulae for one-dimensional network flow problems Carlos F. Daganzo and Karen R. Smilowitz A note on asymptotic formulae for one-imensional network flow problems Carlos F. Daganzo an Karen R. Smilowitz (to appear in Annals of Operations Research) Abstract This note evelops asymptotic formulae

More information

Tractability results for weighted Banach spaces of smooth functions

Tractability results for weighted Banach spaces of smooth functions Tractability results for weighte Banach spaces of smooth functions Markus Weimar Mathematisches Institut, Universität Jena Ernst-Abbe-Platz 2, 07740 Jena, Germany email: markus.weimar@uni-jena.e March

More information

Multiple-mode switched observer-based unknown input estimation for a class of switched systems

Multiple-mode switched observer-based unknown input estimation for a class of switched systems Multiple-mode switched observer-based unknown input estimation for a class of switched systems Yantao Chen 1, Junqi Yang 1 *, Donglei Xie 1, Wei Zhang 2 1. College of Electrical Engineering and Automation,

More information

Nested Saturation with Guaranteed Real Poles 1

Nested Saturation with Guaranteed Real Poles 1 Neste Saturation with Guarantee Real Poles Eric N Johnson 2 an Suresh K Kannan 3 School of Aerospace Engineering Georgia Institute of Technology, Atlanta, GA 3332 Abstract The global stabilization of asymptotically

More information

Lecture 5. Symmetric Shearer s Lemma

Lecture 5. Symmetric Shearer s Lemma Stanfor University Spring 208 Math 233: Non-constructive methos in combinatorics Instructor: Jan Vonrák Lecture ate: January 23, 208 Original scribe: Erik Bates Lecture 5 Symmetric Shearer s Lemma Here

More information

Systems & Control Letters

Systems & Control Letters Systems & ontrol Letters ( ) ontents lists available at ScienceDirect Systems & ontrol Letters journal homepage: www.elsevier.com/locate/sysconle A converse to the eterministic separation principle Jochen

More information

The chromatic number of graph powers

The chromatic number of graph powers Combinatorics, Probability an Computing (19XX) 00, 000 000. c 19XX Cambrige University Press Printe in the Unite Kingom The chromatic number of graph powers N O G A A L O N 1 an B O J A N M O H A R 1 Department

More information

Formation Control and Network Localization via Distributed Global Orientation Estimation in 3-D

Formation Control and Network Localization via Distributed Global Orientation Estimation in 3-D Formation Control and Network Localization via Distributed Global Orientation Estimation in 3-D Byung-Hun Lee and Hyo-Sung Ahn arxiv:1783591v1 [cssy] 1 Aug 17 Abstract In this paper, we propose a novel

More information

Perturbation Analysis and Optimization of Stochastic Flow Networks

Perturbation Analysis and Optimization of Stochastic Flow Networks IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. XX, NO. Y, MMM 2004 1 Perturbation Analysis an Optimization of Stochastic Flow Networks Gang Sun, Christos G. Cassanras, Yorai Wari, Christos G. Panayiotou,

More information

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21 Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting

More information

INVERSE PROBLEM OF A HYPERBOLIC EQUATION WITH AN INTEGRAL OVERDETERMINATION CONDITION

INVERSE PROBLEM OF A HYPERBOLIC EQUATION WITH AN INTEGRAL OVERDETERMINATION CONDITION Electronic Journal of Differential Equations, Vol. 216 (216), No. 138, pp. 1 7. ISSN: 172-6691. URL: http://eje.math.txstate.eu or http://eje.math.unt.eu INVERSE PROBLEM OF A HYPERBOLIC EQUATION WITH AN

More information

Distributed Adaptive Synchronization of Complex Dynamical Network with Unknown Time-varying Weights

Distributed Adaptive Synchronization of Complex Dynamical Network with Unknown Time-varying Weights International Journal of Automation and Computing 3, June 05, 33-39 DOI: 0.007/s633-05-0889-7 Distributed Adaptive Synchronization of Complex Dynamical Network with Unknown Time-varying Weights Hui-Na

More information

Indirect Adaptive Fuzzy and Impulsive Control of Nonlinear Systems

Indirect Adaptive Fuzzy and Impulsive Control of Nonlinear Systems International Journal of Automation an Computing 7(4), November 200, 484-49 DOI: 0007/s633-00-053-7 Inirect Aaptive Fuzzy an Impulsive Control of Nonlinear Systems Hai-Bo Jiang School of Mathematics, Yancheng

More information

Passivity-Based Position Consensus of Multiple Mechanical Integrators with Communication Delay

Passivity-Based Position Consensus of Multiple Mechanical Integrators with Communication Delay To appear in ACC 010 Passivity-Base Position Consensus of Multiple Mechanical Integrators with Communication Delay Ke Huang an Dongjun Lee Abstract We present a consensus framework for multiple variable-rate

More information

Global Solutions to the Coupled Chemotaxis-Fluid Equations

Global Solutions to the Coupled Chemotaxis-Fluid Equations Global Solutions to the Couple Chemotaxis-Flui Equations Renjun Duan Johann Raon Institute for Computational an Applie Mathematics Austrian Acaemy of Sciences Altenbergerstrasse 69, A-44 Linz, Austria

More information