arxiv: v2 [cs.sy] 19 Sep 2018

Transcription

1 Synchronzaton of Kuramoto oscators n a bdrectona frequency-dependent tree network Matn Jafaran, Xne Y, Mohammad Pran, Henrk Sandberg, Kar Henrk Johansson arxv: v [cs.sy] 19 Sep 18 Abstract Ths paper studes the synchronzaton of a fnte number of Kuramoto oscators n a frequency-dependent bdrectona tree network. We assume that the coupng strength of each nk n each drecton s equa to the product of a common coeffcent and the exogenous frequency of ts correspondng source oscator. We derve a suffcent condton for the common coupng strength n order to guarantee frequency synchronzaton n tree networks. Moreover, we dscuss the dependency of the obtaned bound on both the graph structure and the way that exogenous frequences are dstrbuted. Further, we present an appcaton of the obtaned resut by means of an event-trggered agorthm for achevng frequency synchronzaton n a star network assumng that the common coupng coeffcent s gven. I. INTRODUCTION Oscaton s the fundamenta functon behnd the operaton of many compex networks, ncudng boogca and neura networks [17], []. The we-ceebrated Kuramoto oscator [13] has been a paradgm for studyng nterconnected oscators. Kuramoto oscator has been orgnay desgned to study the synchronzaton of couped oscators n chemca networks, and t has been wdey used n other dscpnes ncudng synchronzaton n bran networks [3]. From a technca pont of vew, synchronzaton of Kuramoto oscators have been wdey studed n the terature. The man resuts have focused on the orgna mode of Kuramoto n a compete graph [1] where bounds on the crtca coupng are derved [], [5] to provde suffcent and necessary condtons for frequency synchronzaton. Other reevant probems have aso been studed, for exampe synchronzaton of oscators over genera connected graphs [9], synchronzaton wth tme-varyng exogenous frequences [7], and custer synchronzaton of Kuramoto oscators n a connected, weghted and undrected graphs [6]. Man contrbutons: Ths paper consders frequency synchronzaton of Kuramoto oscators n a bdrectona frequencydependent tree network. We are motvated by the nterest behnd studyng synchronzaton between dfferent areas of a compex bran-ke network. Our choce of studyng tree networks s encouraged by observatons that arge-scae nterarea connectvty n the bran can be approxmated as a tree network [19]. The dea of frequency-dependent coupng stems from the evdence of frequency-dependent synaptc Ths work was supported by the Knut and Ace Waenberg Foundaton, the Swedsh Strategc Research Foundaton and the Swedsh Research Counc. The authors are wth the Automatc Contro Department, Schoo of Eectrca Engneerng and Computer Scence, KTH Roya Insttute of Technoogy, Stockhom, Sweden. Ema: matnj@kth.se, xney@kth.se, pran@kth.se, hsan@kth.se, kaej@kth.se. coupng between neurons e.g. [14]. In the presentaton of the current paper, we are confned wthn the mathematca framework and graph theory to represent our mode wthout a drect usage of termnooges from the doman of neuroscence. We consder a tree graph and assume that each nk (edge) of the graph s bdrectona and ts weght n each drecton depends on a common coupng term as we as the exogenous frequency of the oscator at ts source node []. Whe the oscators exogenous frequences are dfferent from each other, there s a common coupng coeffcent whch affects a nks equay. In other words, each oscator s connected to ts neghbors wth κω, where κ denotes the common stffness and ω vares for each oscator. We derve a suffcent condton on the bound of κ such that the network acheves frequency synchronzaton. We show and dscuss the dependency of the obtaned bound on the exogenous frequences as we as graph structure. Compared wth [] where star graphs wth dentca eaf frequences are consdered, we consder a genera cass of tree graphs where noda exogenous frequences are dfferent from each other and aso use dfferent anaytca toos from contro theory. Compared to the prevous work, e.g. [], [9], [5], [6], [7], we are consderng a frequency-dependent dynamcs for Kuramoto oscators whch has not been consdered before. Moreover, we study synchronzaton n a tree graph (a non-compete graph) and derve a condton on the coupng bound whch depends on the exogenous frequences and the graph structure. In addton to studyng frequency synchronzaton, we present an event-based agorthm for synchronzaton n star networks, whch s a speca case of a tree network, assumng a specfed κ whch may not necesary satsfy the suffcent condton for synchronzaton. Compared wth [16], we consder a dfferent underyng dynamcs for the oscatory network and desgn a dfferent agorthm. Ths paper s organzed as foows. Secton II presents premnares and probem formuaton. Secton III gves a suffcent condton for the common coupng strength n order to acheve frequency synchronzaton. An eventtrggered agorthm for synchronzaton n a star network s presented and anayzed n Secton IV. Secton V presents smuaton resuts and Secton VI concudes the paper.

2 II. PRELIMINARIES AND PROBLEM FORMULATION For a connected undrected graph G(V, E), the node-set V corresponds to n nodes and the edge-set E V V corresponds to m edges. The ncdence matrx B n m assocated to G(V, E) descrbes whch nodes are couped by an edge. Each eement of B s defned as foows +1 f node s at the postve end of edge b = 1 f node s at the negatve end of edge otherwse, where the abeng of the nodes can be done n an arbtrary fashon. The matrx L = BB T s caed the graph Lapacan and L g = B T B s the edge-lapacan. If the underyng graph s connected, the egenvaues of the Lapacan matrx are ordered as = λ 1 (L) < λ (L)... λ n (L), where λ (L) s caed the agebrac connectvty of the network. If the undrected graph s a tree, a egenvaues of L g (whch are n 1 egenvaues) are equa to the nonzero egenvaues of L [15]. In ths paper, we consder trees whch are a subcass of connected graphs wthout cyces,.e., any two nodes are connected by exacty one unque path. The edge Lapacan of a tree graph s nvertbe [15]. Defnton 1 A subset S R n s sad to be forward nvarant wth respect to the dfferenta equaton ẋ = f(x) provded that each souton x( ) wth x() S has the property that x(t) S for a postve t n the doman of defnton of x( ) [18]. Notaton Symbo 1 n s a n-dmensona vector and 1 m m represents a m m matrx whose eements are a equa to 1. The notaton x,j s equvaenty used for x x j. The notaton θ,k ndcates that the node of graph G s connected to the edge k. The mnmum egenvaue of the postve defnte matrx M s denoted by λ 1 (M). A. Probem formuaton Consder n oscators communcatng over a connected and bdrectona graph. The orgna Kuramoto mode foows θ = ω κ sn(θ θ j ), (1) j N where θ R, ω > are the phases and exogenous frequences of oscator, and N denotes the set of neghborng nodes of node. The parameter κ >, κ R s the constant coeffcent of the coupng strength of a nks of the graph. We now contnue wth a dfferent mode where the dynamcs of each node foows [] θ = ω κω sn(θ θ j ). () j N Ths mode can be nterpreted as a bdrectona communcaton where the weghts of coupng of each edge at each drecton depends on the frequency of the source node. Ths makes the nteracton topoogy a drected and weghted graph. The mode of the network n compact form s θ = ω(1 n κb sn(b T θ)), (3) where ω dag(ω 1, ω,..., ω n ) s a dagona matrx such that ω > represents the exogenous frequency of node, θ [θ 1, θ,..., θ n ] T, and sn functon acts eement-wse. As presented n [5], the notons of synchronzaton ncude phase cohesveness,.e., θ θ j η, phase-synchronzaton,.e., θ = θ j,, j, frequency synchronzaton,.e., θ = θ j,, j. In ths paper we frst characterze the suffcent condton for the coupng strength such that for κ >, the phase cohesveness and frequency synchronzaton occur provded that the nta condtons of the oscators are wthn a prescrbed bound. Second, we present an event-trggered mechansm to acheve frequency synchronzaton n a star network. III. THE COUPLING STRENGTH Ths secton studes a suffcent condton for the coupng strength κ to guarantee phase-cohesveness and frequency synchronzaton for the network presented by (3). Assumpton 1 The communcaton topoogy for the network wth node dynamcs n () s a tree. Assumpton The nta reatve phase θ () θ j () [ η, η], j N, where η >, η = π ε for some ε >. Assumpton 3 A exogenous frequences ω are strcty postve,.e., ω ζ >. Snce we are nterested n frequency synchronzaton,.e. θ = θ j, and θ θ j η, et us frst wrte the compact reatve phase dynamcs, B T θ, by mutpyng (3) wth B T as foows B T θ = B T ω1 n κ B T ωb sn(b T θ). (4) Lemma 1 Under Assumptons 1 and 3, the matrx B T ωb s postve defnte and ts smaest egenvaue s postve,.e., λ 1 (B T ωb) >. Proof: Consder x T B T ωbx, x R n and defne y = Bx. We have, x T B T ωbx = y T ωy. Snce ω >, y T ωy ω mn y T y. Thus, x T B T ωbx ω mn y T y. Now y T y = f and ony f y =. If we show that y = mpes x =, then B T ωb >. We argue as foows. Assume that y = Bx =. Then, B T Bx =. Snce for a tree graph, B T B s nvertbe, we concude that x = whch ends the proof. Reca Assumpton and for a gven ε > defne S = {θ R n : θ θ j η, η = π ε, (, j) E}. (5)

3 In addton, take V (θ) = sn T ( BT θ ) sn( BT θ ) and defne S = {θ R n : V (θ) c(ε)} (6) such that S S s the argest eve set of V (θ) that fts n S. Notce that c(ε) > R s defned for the gven ε. π ε γ V k < V k = Proposton 1 For the reatve ange dynamcs n (4) under Assumptons 1-3, consder the set S n (6) and the set S n (5). Then S S s forward nvarant for system (4) provded that κ ω max λ 1(B T ωb) cos(ε), where ω max = max (,j) E ω ω j. Fg. 1. V k < The sgn of V k for θ,j [ π + ε, π ε]. Proof: Consder V = sn T ( BT θ ) sn( BT θ ) as the Lyapunov functon and defne S as n (6). Wrte V n the eement-wse form, V = j sn ( θ θj ), and cacuate V. We obtan V = sn( θ θ j ) cos( θ θ j )( θ θ j ) j = sn(θ θ j )( θ θ j ) (7) j = sn T (B T θ) B T θ. Repacng (4) n the above gves V = sn T (B T θ)b T (ω1 n κωb sn(b T θ)) = sn T (B T θ)b T ω1 n κ sn T (B T θ)b T ωb sn(b T θ). (8) Notce that B T ω1 n ω max1 m where ω max denotes the maxmum of ω ω j and m = n 1 s the number of edges of a tree graph. Defne y = sn(b T θ). Hence, ω maxy T 1 m y T B T ω1 n ω maxy T 1 m, where ω maxy T = ω max y T and y T s the eementwse absoute vaue. Aso, from Lemma 1, we have that λ 1 (B T ωb) >. Appyng the above smpfcatons n (8), we obtan V ω maxy T 1 m κλ 1 (B T ωb)y T y (α y β1 m ) T (α y β1 m ) + mβ m ((α y k β) β ), k=1 where α = κλ 1 (B T ωb), β = ω max, and y k = κλ 1(B T ωb) sn(θ,k θ j,k ) wth k denotng the edge connectng two nodes and j. Defne V k = ((α y k β) β ). Based on (9), we have V k hods f y k = or (9) y k = sn(θ,k θ j,k ) > β, β >, α >, (1) α ω max λ 1(B T ωb) sn(θ,k θ j,k ). that s κ Snce sn(θ,k θ j,k ) can be very sma then there s a regon where V k > (see Fgure 1). Hence, f κ s arger than κ ω max λ 1(B T ωb) sn γ, Vk w be negatve on the set γ < θ,j < π. Now, f we take γ cose enough to π ε, the suffcent condton for S to be forward nvarant s that ω κ max λ 1 (B T ωb) cos(ε), (11) whch ends the proof. For a tree graph, the symmetrc matrx A = B T ωb, whch s n the form of a weghted edge Lapacan, has the foowng structure ω,1 + ω j,1 ω 1,... ω 1,m A m m =...., ω m,1 ω m,... ω,m + ω j,m (1) where ω,k, ω k,p denote the frequency of node connected to the nk k and the frequency of the shared node of two nks k, and p, respectvey. The foowng proposton provdes bounds on λ 1 (B T ωb) based on the network topoogy and exogenous frequences whch can be used n (11). Proposton The mnmum egenvaue of B T ωb, λ 1 (B T ωb), for a tree structure s ower bounded by { } max ω mn λ (BB T ), mn }{{} {( d )ω,k + ( d j )ω j,k } k E }{{} () () λ 1 (B T ωb) mn k E {ω,k + ω j,k }, where d s the degree of node of the underyng graph. (13) Proof: We motvate the reasonng behnd each of the three eements n (13). Reca that based on Lemma 1, x T B T ωbx ax T x where a >. ) We have ω mn I n ω ω max I n. Thus, B T ωb > ω mn B T B > ω mn λ 1 (B T B)I n. Snce for a tree graph B T B s nvertbe [15], we have B T B > λ (L)I n whch gves the frst bound that s rather conservatve. ) The second bound s obtaned as a ower bound on the smaest snguar vaue of a genera symmetrc matrx A m m [1] m σ mn (A) mn{ a, a,j }. j=1,j For the case of A = B T ωb as n (1), the above bound w be a ower bound for ts smaest egenvaue. We have a kk = ω,k +ω j,k where k s the edge whch connects node

4 to node j. Aso, m j=1,j a,j = (d 1)ω,k + (d j 1)ω j,k, and t gves the resut. The upper bound comes from the Rayegh quotent nequaty [8] s obtaned by choosng x = e k where e k s the k-th vector of the canonca bass and k = arg mn k E {ω,k + ω j,k }. A comparson between the tghtness of the bounds proposed n (13) s dscussed n Secton III-A. ω 3 ω 8 ω 3 ω ω 4 ω 1 ω 8 ω 1 ω (a) (c) ω 7 ω 6 ω 5 ω 7 (b) ω 6 ω 5 ω 1 ω ω 3 ω 4 ω 5 ω 6 ω 7 ω 4 Proposton 3 Under Assumptons 1-3, the network n (4) acheves frequency synchronzaton,.e. θ = θ j,, j V f the condton on κ n (11) hods. Proof: Consder the reatve phase dynamcs as n (4). Defne z = B T θ. Cacuatng ż, we obtan ż = κaw cos θ z, (14) where A = B T ωb and W cos θ s a m m dagona matrx wth dagona eements equa to cos(θ,k θ j,k ) where k {1,..., m} denotes the ndex of edge k of the graph. Consder the Lyapunov functon V = z T A 1 z. Assume that the condton on κ n (11) hods. Thus, from proposton 1, we have cos(θ,k θ j,k ) sn(ε) >. Hence, sn(ε)i m W cos θ I m. We obtan V κ sn(ε)z T z. Thus, the system n (14) s unformy asymptotcay stabe whch ends the proof. Remark 1 Notce that f ω = ωi m, ω (, + ), then both the phase and frequency synchronzaton w be acheved κ > [9]. Ths case s dentca wth the orgna Kuramoto mode (as n (1)) wth dentca frequences. A. Exampes: Effects of the graph structure, sze and dstrubuton of exogenous frequences on the suffcent κ Ths secton frst presents some exampes to show that not ony the magntudes of exogenous frequences affect the suffcent bound of κ but aso the way that these frequences are dstrbuted. Moreover, we provde some exampes to compare the tghtness of the bounds proposed n (13) wth respect to the graph structure and sze. Exampe 1 (Effects of the dstrbuton of exogenous frequences): We provde some exampes to show that not ony the magntude of ω max affects κ but aso the way that the exogenous frequences are dstrbuted. We consder the effect of changng a snge exogenous frequency ω on λ 1 (B T ωb) and consequenty on (11). Suppose that the frequences of a nodes are 1, except one whch s 1. We want to assgn ths frequency ω = 1 to one of the nodes n the network such that the resutng egenvaue λ 1 (B T ωb) s maxmzed. For the case of a star graph, f we pace t n the center (hub), λ 1 (B T ωb) = 1 s obtaned and f we put t on one of the eaves, t gves λ 1 (B T ωb) =.16. For the graph shown n Fg. -b, whch conssts of two stars wth the same sze connected va a snge edge, due to the symmetry there are two possbtes; ether ω 1 = 1 or ω 8 = 1. For the former case, λ 1 (B T ωb) =.67 s obtaned, whe the atter Fg.. Dfferent ocatons of the vunerabe node n three specfc graphs. case gves λ 1 (B T ωb) = For the ne graph shown n Fg. -c, the argest vaue for λ 1 (B T ωb) takes pace when we assgn ω 4 = 1, whch gves λ 1 (B T ωb) = 1.98 and t decreases by gong from the center of the ne to one the ends whch gves λ 1 (B T ωb) =.86. We observe that n the above exampes the optma pace of the vunerabe node shoud be one of the we-known centrates. Notce that the vaues of λ 1 n these exampes are reported based on the exacty cacuated λ 1 (B T ωb) for the gven graphs. Usng the estmated vaue of λ 1 n (13) eads to a smar concuson. For exampe, for the star graph n Fg. -a, the bound n (13) gves 5 λ 1 11 f ω 1 = 1 and 1 λ 1 11 f the exogenous frequency of one of eaves s equa to one. Exampe (Comparng bounds n (13) and the effect of network sze and structure): For the graph shown n Fg. 3- a the ower bound () gves 1, whe bound () s 99. Hence, bound () s tghter. If we keep ncreasng the number of eaves up to 1, the ower bound () s st tghter than () (for Fg. 3-b, () gves 1 and () gves ). For graph Fg. 3-c, bound () gves 1 and bound () gves 49. Hence, dependng on the network structure, ether of the two ower bounds become tghter. We shoud note that the argest agebrac connectvty among a trees beongs to star graphs, whch s 1. For most of the tree structures, λ (L) scaes wth the sze of the network, e.g., ne graphs. Fg. 3-d shows the roe of network sze on the scang of bound (). In ths exampe bound () gves zero. The graph topoogy s a ne graph and we know that for these ( graphs λ (L) = 1 cos( π n ) [1]. Thus, bound () gves ω mn 1 cos( π n )) whch bgger than zero, athough t goes to zero as the network sze grows. IV. APPLICATION: AN EVENT-BASED ALGORITHM FOR FREQUENCY SYNCHRONIZATION IN A STAR NETWORK In ths secton, we present an appcaton of usng the suffcent coupng strength obtaned n Proposton 1 by means of a centrazed event-based agorthm (e.g. [4]) for synchronzaton n a star network. Foowng the prevous secton, we are nterested n frequency synchronzaton and phase-cohesveness,.e. θ j θ < η and θ = θ j. To ths purpose, we assume that the coupng stffness κ s gven. Hence, to derve the network towards frequency synchronzaton, we manpuate ω usng the resuts of the

5 prevous secton such that for each edge θ θ j < η (η = π ε, ε > ) s enforced. Consder a star graph (see Fg. -a) wth the node dynamcs as foows θ h = ω h (1 κ sn(θ h θ )), θ = (ω + α )(1 + κ sn(θ h θ )), (15) where the centra node of the star s caed hub, denoted by θ h, and other nodes are caed eaf, denoted by θ. The probem s how to desgn α n order to acheve our goa. One approach coud be based on desgnng an α controer to contnuouesy reguate ω, e.g. [11]. In ths paper, however, we are nterested n an event-based approach whch does not requre updatng the contro acton for a tmes. Assumpton 4 We assume that 1) the exogenous frequences ω j wth j {h, } are sow enough to be estmated as constant for a arge enough perod of tme, ) reatve phases and ther dervatves,.e., θ θ h and θ θ h, are known to the hub. Under the above assumpton, we desgn a centrazed event-trggered agorthm such that the hub updates α = (α,..., α n ) T for achevng frequency synchronzaton. For system (15), we defne a set of trggerng tmes t h, t1 h, t h,... at whch the vector α gets updated, such that α(t) = α(t k h ), t [tk h, tk+1 h ). Before presentng the agorthm, we frst defne the events (E 1, E ) and the requred cacuatons for updatng α. Event E1 s actvated f the reatve phase of nodes and h s arger than a prescrbed mt. In ths case, the eaf exofrequency w be updated by adjustng α whe the hub exo-frequency s kept unchanged. Event E 1 s actvated f there exsts one edge whch meets the trggerng condton. Hence, E 1 = n E1, where E1 denotes event E 1 for node =. Besdes event E 1, event E s desgned to update the trggerng condton of each edge n order to avod chatterng (repettve swtchngs) of E1 [1] (see Remark ). We wrte E = n E. = Event E 1 (Update of α ): {,..., n} s.t. θ h θ > ω 4 = 1 ω 1 = 1 ω = 1 ω 3 = 1 ω 1 = 1 ω 4 = 1 ω 4 = 1 ω 1 = 5 ω 3 = 1 ω = 1 ω 1 = 1 ω 99 = 1 ω 3 = 1 ω = 1 (a) (b) (c) (d) ω 1 ω ω 3 ω 4 ω 1 Fg. 3. Exampes whch show the tghtness of bounds proposed n (13) and the roe of the network scang and structure on each ower bound. η, wth η = π ε >. Let α denote α (t k h ) and α+ denote α (t k+1 h ) for node. The update of α (.e., α + ) shoud guarantee that the nk h, s contrbutng to the decrease of the overa Lyapunov functon of the system (see the proof of Proposton 1). Hence, the suffcent condton on κ shoud ocay hod. Consderng (11), α + shoud ocay guarantee that κ > ω h (ω + α + ) λ 1 (B T ωb) cos(ε). In the above, consderng cos(ε) guarantees that V w be ocay negatve for θ h, > π ε, and thus t w negatve for θ h, > η, η = π ε (see Fg 1). To ocay estmate λ 1 (B T ωb), whch s requred for updatng α, we use the estmaton based on ω mn λ (BB T ) (see n (13)). Defne ωmn L = mn{ω h, ω }, where ωmn L denotes the oca estmaton of ω mn and ω s the desred vaue for ω + α +. Snce ω s a desgn choce, we opt for the case where ω > ω h (motvated by the exampes n Secton III-A). Notce that for a star graph λ (BB T ) = 1. Hence, from (11), κ > ω ω h ω h cos(ε) shoud hod. The hub then cacuates and updates ω = ω h (1 + cos(ε)), (16) α + = ω (ω + α ). (17) Notce that n cacuaton of ω, t s assumed that ω h s known to the hub. In fact, the structure of the star graph together wth Assumpton 4- aow the hub to cacuate both ω h and ω by usng the frst and second dervatves of θ h θ. Now, we present our event-trggered agorthm as foows. Consder Assumpton 4. Let t denote the nta tme and {t k h } k=1 denote the trggerng tmes of the hub determned by E1 or E (see Remark ). Agorthm 1: 1: Choose η > and ε >. : Intaze α =, t h = t and k =. 3: From tme s = t k h, the hub contnuousy senses θ h θ to detect E = E 1 E, where E 1 = n E1 and E = = n E, such that τ = nf{r s : s.t. E1 or E}, = where E1 mpes {, θ h θ > η} and E mpes {, θ h θ < η ε}. 4: If E, the hub determnes t k+1 h = τ, and for every whch meets E1, the hub updates α based on (17) and repaces the defnton of E1 wth E. Aso, for every j whch meets E j, the hub repaces the defnton of Ej wth E j 1. The hub goes back to Step 3. Remark [Avodng chatterng] After trggerng E 1, the same event can be mmedatey trggered snce there w be a t t the condton of event E 1 s voated. To prevent ths behavor, event E s ntroduced. In fact, E w temporary repace E 1 and thus w avod chatterng of ths event. Notce that athough t s possbe that two dfferent edges trgger

6 d3/dt d3/dt an event at the same tme (for exampe edge trggers E1 and edge j trggers E j smutaneousy), t s mpossbe that E1, E occur for the same edge smutaneousy. System (15) wth the above event-trggered agorthm can be represented as the foowng hybrd system wth state (θ, α) such that the contnuous evouton of the system obeys θ = ω α (1 n κb sn(b T θ)), α =, (18) where ω α s a dagona matrx whose dagona eements are ω α 1 = ω h and ω α = ω + α. Aso, f there s a nk whch meets the jump condton, the foowng dscrete transton occurs θ + = θ, α + = ω (ω + α ). (19) Proposton 4 There s a ower bound on the nter-trggerng tmes of the soutons to hybrd system (18)-(19) consderng a star graph wth Assumptons -4. Moreover, the set S S, wth S n (6) and S n (5), s forward nvarant for system (18)-(19). Proof: Frst we prove that S s forward nvarant. Snce at the swtches the state θ stays unchanged, we take a smar Lyapunov functon as n the proof of Proposton 1. Snce event E 1 prevents the phase dfferences to grow arger than η, and aso the correctng terms α are desgned such that the suffcent condton on κ s respected, then based on the argument n the proof of Proposton 1, the set S s forward nvarant for (18)-(19). To prove the exstence of a non-zero dwe tme, we argue that an edge cannot trgger E 1 and E at the same tme (see Remark ). Aso, after trggerng ether E 1, at east a dstance (ε) shoud be paved wth a bounded veocty (snce θ s bounded) whch mpes E does not nstantaneousy happen after E 1. The atter ndcates that there exsts a non-zero dwe tme between each two trggerng tmes. Notce that at each trggerng tme, t s possbe to have more than one edge that trgger an event but number of nodes f fnte. Hence, there s no zeno behavor and the soutons can evove n tme. Proposton 5 Under Assumptons -4, a star network wth the dynamcs as n (18)-(19) wth the event-trggered Agorthm 1 acheves frequency synchronzaton. Proof: Defne z = B T θ. Foowng (14), consder ż = κa α W cos θ z, () where A α = B T ω σ B and W cos θ s defned as before. We now argue that there exsts a tme T at whch a noda exogenous frequences are ether updated or w stay unchanged. We reason as foows. If θ h θ exceeds the desgned upper-bound, an event w be generated and the exogenous frequency of node w be repaced by ω. If no event s generated by node, we concude that θ h θ s sma enough and wthn the desred bound. Snce, the number of nodes are fnte, there exts a fnte tme T at whch the exogenous frequences w not get updated anymore. We abe the exogenous noda frequency of node after tme T by ω. Take V = z T M 1 z wth M = ω h 1 m m +(ω )I m, where M > has the structure of A α (see (1)), as the Lyapunov canddate. Notce that for t T, M = A α snce ω σ w be constant. Thus, at jumps, θ (hence z) w stay unchanged. Aso, { t T, whch ends the proof. V κ sn εz T z}, V. SIMULATION RESULTS Ths secton presents smuaton resuts for a network of four oscators over a star and a ne graph topoogy. The nta condton for the oscators s set to θ() = [ π 4, π 1, π, π 5 ]. We smuate both star and ne networks for two sets of exogenous frequences ω = [, 3,, 1] and ω = [1, 1, 5, 6]. For the star graph node 1 concdes wth the hub and for the ne graph node 1 and 4 are termna nodes. Tabe V shows the exact vaue of the mnmum egenvaue for each of the cases (the ncdence matrx for the star graph s denoted by B s and for the ne graph wth B ), together wth the bound obtaned based on (13) and a suffcent bound for κ. The atter s cacuated usng the exact vaue of λ mn reported n the second coumn of the tabe. As shown, the ower bound of κ wth ω, where the hub frequency s mnmum, s smaer (for both star and ne graph) than ω (see Secton III-A). Choce λ mn Estmaton κ n (13) Bs T ωb s Bs T ωb s B T ωb B T ωb Tabe 1 Fgure 4 (Fg. 5) shows the reatve phases and noda frequences for two sets of exogenous frequences ω and ω for a star (ne) graph. We used κ = 5 for a four cases. As shown n Fgure 5, the ne graph aso acheves frequency synchronzaton wth a κ smaer than the suffcent bound obtaned n Tabe V. In a cases, reatve phases converge to a non-zero vaue and a noda frequences reach a consensus. To smuate the resuts of the event-trggered Fg. 4. A star graph wth 4 nodes and two sets of exogenous frequences ω (eft) and ω (rght). agorthm, we frst take a star graph wth 4 nodes wth

7 d 3 /dt d3/dt d3/dt Fg. 5. A ne graph wth 4 nodes and two sets of exogenous frequences ω (eft) and ω (rght). the same nta condtons as the above exampes. We set ω = [, 18, 16, 6] and κ = 1.1. The resuts s shown n Fgure 6-(1). As shown the frequences de-synchronze and the reatve phases are unbounded. For the same network, we use the event-trggered contro wth Agorthm 1, and set = 1.1, η = π π 1, and ε = π 1. The resuts are shown n Fgure 6-(). As shown the oscators synchronze. 3 h Fg. 6. A star graph wth (1) frequency de-synchronzaton (eft) and () event-trggered controer (rght). VI. CONCLUSIONS Ths paper has studed the synchronzaton of a fnte number of Kuramoto oscators n a tree network where the coupng strength of each nk between every two oscators n each drecton s weghted by a common coeffcent and the exogenous frequency of ts correspondng source oscator. We have drven a suffcent condton for the common coupng strength and showed ts dependency on both exogenous frequences and graph structure. We have aso provded an exampe of the appcaton of the obtaned suffcent bound to acheve frequency synchronzaton n a star network for whch the vaue of κ was gven. Future avenues ncude aowng zero and tme-varyng exogenous frequences, and extendng the event-trggered agorthm for a more genera cass of graphs. [4] D.V. Dmarogonas, E. Frazzo, and K.H. Johansson. Dstrbuted event-trggered contro for mut-agent systems. IEEE Transactons on Automatc Contro, 57(5): , 1. [5] F. Dörfer and F. Buo. On the crtca coupng for Kuramoto oscators. SIAM Journa on Apped Dynamca Systems, 1(3):17 199, 11. [6] C. Favaretto, A. Cenedese, and F. Pasquaett. Custer synchronzaton n networks of Kuramoto oscators. IFAC-PapersOnLne, 5(1): , 17. [7] A. Franc, A. Chaet, and W. Pasas-Lépne. Phase-ockng between Kuramoto oscators: robustness to tme-varyng natura frequences. In Decson and Contro (CDC), 1 49th IEEE Conference on, pages IEEE, 1. [8] R. A. Horn and C. R. Johnson. Matrx Anayss. Cambrdge Unversty Press, 199. [9] A. Jadbabae, N. Motee, and M. Barahona. On the stabty of the Kuramoto mode of couped nonnear oscators. In Proceedngs of the Amercan Contro Conference, voume 5, pages IEEE, 4. [1] M. Jafaran. Robust consensus of uncyces usng ternary and hybrd controers. Internatona Journa of Robust and Nonnear Contro, 7(17): , 17. [11] M. Jafaran, E. Vos, C. De Perss, J. Scherpen, and A. J. van der Schaft. Dsturbance rejecton n formaton keepng contro of nonhoonomc wheeed robots. Internatona Journa of Robust and Nonnear Contro, 6(15): , 16. [1] C.R. Johnson. A gersgorn-type ower bound for the smaest snguar vaue. Lnear Agebra and ts Appcatons, 11:1 7, [13] Y. Kuramoto. Chemca oscatons, waves, and turbuence, voume 19. Sprnger Scence & Busness Meda, 1. [14] H. Markram, A. Gupta, A. Uze, Y. Wang, and M. Tsodyks. Informaton processng wth frequency-dependent synaptc connectons. Neuroboogy of earnng and memory, 7(1-):11 11, [15] M. Mesbah and M. Egerstedt. Graph theoretc methods n mutagent networks. Prnceton Unversty Press, 1. [16] A.V. Proskurnkov and M. Cao. Synchronzaton of puse-couped oscators and cocks under mnma connectvty assumptons. IEEE Transactons on Automatc Contro, 6(11): , 17. [17] P. Sacre and R. Sepuchre. Senstvty anayss of oscator modes n the space of phase-response curves: Oscators as open systems. IEEE Contro Systems, 34():5 74, 14. [18] E.D. Sontag. Structure and stabty of certan chemca networks and appcatons to the knetc proofreadng mode of t-ce receptor sgna transducton. IEEE Transactons on Automatc Contro, 46(7):18 147, 1. [19] C.J. Stam, P. Teware, E. Van Deen, E.C.W. Van Straaten, A. Hebrand, and P. Van Meghem. The trees and the forest: characterzaton of compex bran networks wth mnmum spannng trees. Internatona Journa of Psychophysoogy, 9(3):19 138, 14. [] E. Steur, I. Tyukn, and H. Njmejer. Sem-passvty and synchronzaton of dffusvey couped neurona oscators. Physca D: Nonnear Phenomena, 38(1):119 18, 9. [1] S.H. Strogatz. From Kuramoto to Crawford: exporng the onset of synchronzaton n popuatons of couped oscators. Physca D: Nonnear Phenomena, 143(1-4):1,. [] C. Xu, Y. Sun, J. Gao, T. Qu, Z. Zheng, and S. Guan. Synchronzaton of phase oscators wth frequency-weghted coupng. Scentfc reports, 6, 16. REFERENCES [1] N. M. M. Abreu. Od and new resuts on agebrac connectvty of graphs. SIAM Journ. Matrx Ana. App, 43:53 73, 6. [] N. Chopra and M.W. Spong. On exponenta synchronzaton of Kuramoto oscators. IEEE transactons on Automatc Contro, 54(): , 9. [3] D. Cumn and C.P. Unsworth. Generasng the Kuramoto mode for the study of neurona synchronsaton n the bran. Physca D: Nonnear Phenomena, 6(): , 7.