arxiv: v2 [math.oc] 2 Mar 2015

Size: px
Start display at page:

Download "arxiv: v2 [math.oc] 2 Mar 2015"

Transcription

1 Disribued even-riggered communicaion for dynamic average consensus in neworked sysems Solmaz S. Kia a Jorge Corés b Sonia Marínez b a Deparmen of Mechanical and Aerospace Engineering, Universiy of California a Irvine, Irvine, CA 9697, USA arxiv:44.747v [mah.oc] Mar 5 b Deparmen of Mechanical and Aerospace Engineering, Universiy of California a San Diego, La Jolla, CA 993, USA Absrac This paper presens disribued algorihmic soluions ha employ opporunisic iner-agen communicaion o achieve dynamic average consensus. In our soluions each agen is endowed wih a local crierion ha enables i o deermine wheher o broadcas is sae o is neighbors. Our saring poin is a coninuous-ime disribued coordinaion sraegy ha, under coninuous-ime communicaion, achieves pracical asympoic racking of he dynamic average of he ime-varying agens reference inpus. Then, for his algorihm, depending on he direced or undireced naure of he ime-varying ineracions and under suiable conneciviy condiions, we propose wo differen disribued even-riggered communicaion laws ha prescribe agen communicaions a discree ime insans in an opporunisic fashion. In boh cases, we esablish posiive lower bounds on he iner-even imes of each agen and characerize heir dependence on he algorihm design parameers. This analysis allows us o rule ou he presence of Zeno behavior and characerize he asympoic correcness of he resuling implemenaions. Several simulaions illusrae he resuls. Key words: cooperaive conrol, dynamic average consensus, even-riggered communicaion, weigh-balanced direced graphs. Inroducion The dynamic average consensus problem consiss of designing a disribued algorihm ha allows a group of agens o rack he average of individual ime-varying reference inpus, one per agen. This problem has applicaions in numerous areas ha involve disribued sensing and filering, including disribued racking [Yang e al., 7], muli-robo coordinaion [Yang e al., 8], sensor fusion [Olfai-Saber, 7, Olfai-Saber and Shamma, 5], and disribued esimaion [Carron e al., 3]. Our goal here is o develop algorihmic soluions o he dynamic average consensus problem which rely on agens locally deciding when o share informaion wih heir neighbors in an opporunisic fashion for greaer efficiency and energy savings. By opporunisic, we mean ha he ransmission of informaion o he neighbors should happen a imes when i is needed o preserve he sabiliy and convergence of he coordinaion algorihm. A preliminary version appears a he IEEE Conference on Decision and Conrol as [Kia e al., 4a]. addresses: solmaz@uci.edu (Solmaz S. Kia), cores@ucsd.edu (Jorge Corés), soniamd@ucsd.edu (Sonia Marínez). Lieraure review: The lieraure of cooperaive conrol has proposed dynamic average consensus algorihms ha are execued eiher in coninuous-ime [Bai e al.,, Freeman e al., 6, Kia e al., 4b, Olfai- Saber and Shamma, 5, Spanos e al., 5] or in fixed sepsize discree-ime [Kia e al., 4b, Zhu and Marínez, ]. Coninuous-ime algorihms operae under he assumpion of coninuous agen-o-agen informaion sharing. Alhough discree-ime algorihms are more amenable o pracical implemenaion, use of a fixed communicaion sep-size, which should be designed o address also rarely occurring wors-case siuaions, can be a waseful use of he nework resources. In addiion, in hese discree-ime algorihms communicaion and compuaion sepsizes are ied ogeher, resuling in poenially a conservaive sepsize for communicaion imes. This can resul in cosly implemenaions, as performing communicaion usually requires more energy han compuaion. Moreover, he assumpion of periodic, synchronous communicaion is unrealisic in many scenarios involving cyber-physical sysems, as processors are subjec o naural delays and errors which may deviae hem from he perfec operaional condiions he sraegies are designed for. Even-riggered communicaion offers a way o address hese shorcomings by prescribing in an opporunisic Preprin submied o Auomaica 8 Ocober 8

2 way he imes for informaion sharing and allowing individual agens o deermine hese auonomously. In recen years, an increasing body of work ha seeks o rade compuaion and decision-making for less communicaion, sensing or acuaion effor while guaraneeing a desired level of performance has emerged, see e.g., [Heemels e al.,, Mazo and Tabuada,, Wang and Lemmon, ]. Closes o he problem considered here are he works ha sudy even-riggered communicaion laws for saic average consensus, see e.g., [Dimarogonas e al.,, Fan e al., 3, Garcia e al., 3, Nowzari and Corés, 4, Seyboh e al., 3] and references herein. Saemen of conribuions: We propose novel algorihmic soluions o he dynamic average consensus problem ha employ opporunisic sraegies o deermine he communicaion imes among neighboring agens. The basic idea is ha agens share heir informaion wih neighbors when he uncerainy in he oudaed informaion is such ha he monoonic convergen behavior of he overall nework can no longer be guaraneed. To realize his concep, depending on he conneciviy properies of he ineracion opology, we propose and characerize he asympoic correcness of wo differen disribued even-riggered communicaion laws. Our leas sringen conneciviy condiions are modeled by a ime-varying, weigh-balanced piecewise consan digraph which is joinly srongly conneced over an infinie sequence of coniguous and uniformly bounded ime inervals. In he second scenario, we consider ineracion opologies modeled by a ime-varying, piecewise coninuous undireced conneced graph, which allows us o furher refine our analyical guaranees. By esablishing posiive lower bounds on he iner-even imes of each agen for boh cases, we also show ha he proposed disribued even-riggered communicaion laws are free from Zeno behavior (he undesirable siuaion where an infinie number of communicaion rounds are riggered in a finie amoun of ime). Finally, we analyze he dependence of he iner-even imes on he algorihm design parameers. This characerizaion provides guidelines on he rade-offs beween he minimum iner-even imes for communicaion and he performance and energy efficiency of he proposed algorihms. We demonsrae hrough several comparaive simulaion sudies he advanages of our proposed even-riggered communicaion sraegies over schemes ha rely on coninuousime communicaion as well as discree-ime communicaion wih fixed sepsize. Noaion and erminology We le R, R, R >, and Z denoe he se of real, nonnegaive real, posiive real, and nonnegaive ineger, respecively. The ranspose of a marix A is A. We le n (resp. n ) denoe he vecor of n ones (resp. n zeros). We le Π n = I n n n n, where I n is he n n ideniy marix. When clear from he conex, we do no specify he marix dimensions. For u R d, u = u u is he sandard Euclidean norm. For u R, u is is absolue value. For a ime-varying measurable locally essenially bounded signal u : R R m, we denoe by u ess he essenial supremum norm. For a scalar signal u, we use u ess insead. For vecors u,..., u m, we le u = (u,..., u m ) represen he aggregaed vecor. In a neworked sysem, we disinguish he local variables a each agen by a superscrip, e.g., x i is he local sae of agen i. If p i R d is a variable of agen i, he aggregae of a nework wih N agens is p = (p,..., p N ) (R d ) N. Graph heory: Here, we briefly review some basic conceps from graph heory and linear algebra following [Bullo e al., 9]. A direced graph, or simply a digraph, is a pair G = (V, E), where V = {,..., N} is he node se and E V V is he edge se. For an edge (i, j) E, i is called an in-neighbor of j and j is called an ou-neighbor of i. We le N i denoe he se of ou-neighbors of i V. A graph is undireced if (i, j) E when (j, i) E. A direced pah is a sequence of nodes conneced by edges. A digraph is called srongly conneced if for every pair of verices here is a direced pah connecing hem. Given digraphs G i = (V, E i ), i {,..., m}, heir union is he graph n i= G i = (V, E E E m ). A weighed digraph is a riple (V, E, A), where G = (V, E) is a digraph and A R N N is a weighed adjacency marix wih he propery ha a ij > if (i, j) E and a ij =, oherwise. A weighed digraph is undireced if a ij = a ji for all i, j V. We refer o a srongly conneced and undireced graph as conneced. The weighed ou- and in-degrees of a node i are, respecively, d i ou = j= a ij and d i in = j= a ji. A digraph is weigh-balanced if a each node i V, he weighed ou- and in-degrees coincide (alhough hey migh be differen across differen nodes). The (ou-) Laplacian marix is L = D ou A, where D ou = Diag(d ou,, d N ou) R N N. Noe ha L N =. A digraph is weigh-balanced if and only if T NL = if and only if Sym(L) = (L LT ) is posiive semi-definie. Based on he srucure of L, a leas one of he eigenvalues of L, denoed by λ,..., λ N, is zero and he res of hem have nonnegaive real pars. We le λ = and R(λ i ) R(λ j ), for i < j, where R( ) denoes he real par of a complex number. We denoe he eigenvalues of Sym(L) by ˆλ,..., ˆλ N. For a srongly conneced and weigh-balanced digraph, zero is a simple eigenvalue of boh L and Sym(L). In his case, we order he eigenvalues of Sym(L) as = ˆλ < ˆλ ˆλ 3... ˆλ N. Throughou he paper, we deal wih ime-varying digraphs wih fixed verex se. A ime-varying digraph R G() = (V, E(), A()) is piecewise coninuous (respecively, piecewise consan) if he map A() is piecewise coninuous (respecively, piecewise consan) from he righ. In such case, we denoe by {s k } k Z

3 he ime insans a which his map is disconinuous and refer o hem as swiching imes. By convenion, s =. A ime-varying digraph R G() has uniformly bounded weighs if, for all R, < a a ij () [a, ā], wih < a ā, if (j, i) E(), and a ij = oherwise. A ime-varying digraph R G() is srongly conneced if each G() is srongly conneced, and is joinly srongly conneced over [, ) if [, )G() is srongly conneced. A piecewise consan ime-varying digraph R G() is recurrenly joinly srongly conneced if he sequence of inerswiching imes {s k s k } k Z is uniformly lower bounded and here exiss an infinie sequence of coniguous uniformly bounded inervals {[s kj, s kj )} j Z, wih s k = s, such ha G() is joinly srongly conneced over [s kj, s kj ), for all j Z. Finally, a imevarying digraph R G() is weigh-balanced if each G() is weigh-balanced. For a piecewise consan recurrenly joinly srongly conneced and weighbalanced digraph wih uniformly bounded weighs, we le {s i k } k Z {s k } k Z be he imes when an agen i V acquires a new in-neighbor and we define L =sup{ L R }, ˆλ =inf{ˆλ (L ) R }, d i ou =sup{d i ou() R }, i V, where L is he Laplacian of G(). Noe ha if G has uniformly bounded weighs and is weigh-balanced and srongly conneced, hen ˆλ >. If G is a ime-varying conneced graph, we use he noaion λ insead of ˆλ. The following resul, aken from [Kia e al., 4b, Lemma 4.5], is useful when dealing wih recurrenly joinly srongly conneced digraphs. When such digraphs are weigh-balanced wih uniformly bounded weighs, here exis ˆλ σ > and ρ > such ha e βr L R( ) ρ e βˆλ σ( ),, () for any β >. If he digraph is addiionally srongly conneced, hen () is saisfied wih ρ = and ˆλ σ = ˆλ. 3 Nework model and problem saemen Consider a nework of N agens wih single-inegraor dynamics, ẋ i = g i, i V, where x i R is he agreemen sae and g i R is he driving command of agen i. Our consideraion of simple dynamics is moivaed by he fac ha he sae of he agens does no necessarily correspond o some physical quaniy, bu insead o some logical variable on which agens perform compuaion and processing. Each agen i V has access o a imevarying reference signal r i : R R. Agens ransmi informaion o oher agens hrough wireless communicaion and heir ineracion opology is modeled by a ime-varying weighed digraph G. An edge (i, j) from i o j a ime means ha agen j can send informaion o agen i a. For convenience, we le ˆx i denoe he las sampled sae of agen i V. We le { i k} k Z = { i k} k Z {s i k} k Z R denoe he sequence of updae imes of agen i V. Here, { i k } k Z is he se of imes a which he sae of he agen is sampled, ha is ˆx i () = x i ( i k ) for [ i k, i k ), while {si k } k Z is he se of swiching imes of he underlying communicaion digraph a which he agen acquires an in-neighbor. According o his model, agen i broadcass ˆx i o is in-neighbors a each ime i k if di in (i k ) > and also a he imes si k. The variable x i () = ˆx i () x i () denoes he mismach beween he las sampled sae and he sae of agen i a ime. When he graph is no fixed, he model above assumes ha individual agens are made aware of he ideniy of newly acquired in-neighbors (so ha hey can communicae o hem heir las ransmied sae) and depared ou-neighbors (so ha hey can remove heir sae from heir compuaions). Under he nework model described above, our goal is o design a disribued algorihm ha allows each agen o asympoically rack he average of he reference inpus N j= rj () across he group. The algorihm design amouns o specifying, for each agen i V, a suiable disribued driving command g i : R N i R ogeher wih a mechanism for riggering communicaion wih is in-neighbors in an opporunisic fashion. By disribued, we mean ha each agen only needs o receive informaion from is ou-neighbors o evaluae g i and he communicaion riggering law. A key requiremen on he communicaion riggering mechanism is ha he resuling nework evoluion is free from Zeno behavior, i.e., does no exhibi an infinie amoun of communicaion rounds in any finie amoun of ime. 4 Coninuous-ime compuaion wih disribued even-riggered communicaion Here, we presen our soluion o he problem saed in Secion 3. Our saring poin is he coninuous-ime algorihm for dynamic average consensus proposed in our previous work [Kia e al., 4b], v i = αβ j= a ij(x i x j ), (a) ẋ i =ṙ i α(x i r i ) β a ij(x i x j ) v i, (b) j= for each i V. In (), α, β R > are design parameers. Noe ha he execuion of his algorihm requires coninuous agen-o-agen sharing of he variable x. Remark 4. (Knowledge of derivaive of reference signals): Ineresingly, () can be execued wihou knowledge of he ime derivaives of he reference signals. In fac, wih he change of variables x = x r, () reads v i = αβ j= a ij( x i r i x j r j ), x i = α x i β a ij( x i r i x j r j ) v i, j= x i = x i r i, i V. 3

4 Agens operae on heir corresponding componens of v and x, and inerchange wih heir neighbors he corresponding componens of x. We only use he represenaion () for convenience in our echnical analysis laer. The following resul summarizes, for reference, he asympoic correcness guaranees of (). Theorem 4. (Convergence of () over weigh-balanced and recurrenly joinly srongly conneced digraphs [Kia e al., 4b]): Assume he agen inpus saisfy Π N ṙ ess = γ <. Le he communicaion opology be a weigh-balanced and recurrenly joinly srongly conneced ime-varying digraph G wih uniformly bounded weighs. Then, for any α, β R >, he evoluion of he algorihm () over G iniialized a z i (), v i () R wih i= vi () = is bounded and saisfies, for i V, lim sup x i () N j= rj () ρ γ, (3) βˆλ σ wih ˆλ σ and ρ saisfying (). Given he nework model of Secion 3, where he ransmission of informaion is limied o discree insans of ime, we propose here he following implemenaion of () wih discree-ime communicaion, v i = αβ j= a ij(ˆx i ˆx j ), (4a) ẋ i = ṙ i α(x i r i ) β a ij(ˆx i ˆx j ) v i, (4b) j= for each i V. Our ask is o provide individual agens wih riggers ha allow hem o deermine in an opporunisic fashion when o ransmi informaion o heir in-neighbors. The design of such riggers is challenging because riggers need o be disribued, so ha agens can check hem wih he informaion available o hem from heir ou-neighbors, hey mus guaranee he absence of Zeno behavior, and hey have o ensure he nework achieves dynamic average consensus even hough agens operae wih oudaed informaion while inpus are changing wih ime. 4. Compac-form algorihm represenaions Here, we presen wo equivalen compac-form represenaions of he algorihm (4) for analysis purposes. Firs, consider he change of variable w =v απ N r, y =x r, r = N j= rj N. (5a) which ransforms he algorihm () ino (compac form) ẇ = αβl y αβl x απ N ṙ, ẏ = αy βl y βl x Π N ṙ w, (6a) (6b) Noe ha y is he aggregaed racking error vecor. Here, we have used L ˆx = L (x x) = L yl x wih x = ˆx x. We use he orhonormal ransfer marix T R N N T=[r R], r = N N, r R=, R R = I N. (7) o obain our second represenaion below ha separaes ou he consan dynamics of he algorihm, q = r w, q :N = αr y R w, z = T y. (8) We pariion he new variable z as (z, z :N ), where z R. Then, if he nework ineracion opology is weighbalanced, he algorihm (6) can be wrien as, q =, (9a) q :N = αq :N, (9b) ż = αz q, (9c) ż :N = βr L Rz :N βr L xr ṙ q :N. (9d) We close his secion by describing he relaionship beween he iniial condiions of he variables for each represenaion. We invoke hese relaions in our analysis below. Noe ha q :N = R (αyw) = R (α(x r)v). Then, given x(), v() R N wih i= vi () =, and using r Π N = and RR = Π N = Π N, q () = r w() = r v() =, q :N () = απ N (x() r()) v(), z () = r y() = r (x() r()), z :N () = Π N (x() r()). (a) (b) (c) (d) 4. Communicaion riggering law for weigh-balanced and recurrenly joinly srongly conneced digraphs In his secion, for neworks wih ime-varying digraph ineracions, we inroduce a disribued even-riggered mechanism ha agens can employ o deermine heir sequence of communicaion imes. For each agen, he execuion of his mechanism relies on is local variables. This naurally resuls in asynchronous schedules of communicaion, which poses addiional analysis challenges. Neverheless, we are able o overcome hem in he following resul which saes ha he closed-loop nework execuion is free from Zeno behavior and guaraneed o achieve pracical dynamic average consensus. Theorem 4.3 (Convergence of (4) over recurrenly joinly srongly conneced and weigh-balanced digraph wih asynchronous disribued even-riggered communicaion): Assume ha he inpu of each agen i V saisfies ṙ i ess = κ i <, while Π N ṙ ess = γ <. Le he communicaion opology be a weigh-balanced and recurrenly joinly srongly conneced ime-varying digraph {G} {sk } wih uniformly bounded weighs. For ɛ R N >, consider an implemenaion of he algorihm (4) over {G} {sk }, where agens communicae according o he model described in Secion 3 wih he sampling imes { i k } k Z of agen i V, saring a i =, deermined by i k =argmax{ [ i k, ) x i ( i k) x i () ɛ i }. () Then, for any α, β >, he evoluion saring from x i () R and v i () R wih i= vi ()= saisfies lim sup x i () N j= rj () γβ L ɛ βˆλ σ ρ, () 4

5 for i V wih an exponenial rae of convergence of min{α, βˆλ σ }. Here, ˆλ σ and ρ saisfy (). Furhermore, he iner-execuion imes of () for each agen i V are lower bounded by τ i = ( α ln αɛ ) i c i, (3) where c i = κ i (α β d i ou) η r (x() r()) and Π N (α(x() r()) v()) αη, (4) η = ρ (γβ L ɛ ) Π βˆλ N (x() r()) ρ q σ :N () (( βˆλ σ α βˆλ σ α ) β ˆλσ α β ˆλσ ( βˆλ σ α ) α α β ˆλσ ), if βˆλ σ α,, if βˆλ βˆλ σ = α. σe Hence, { k } k Z = N i= k Z i k has no accumulaion poin and he execuion of (4) over {G} {sk } is Zeno-free. PROOF. Given an iniial condiion, le [, T ) be he maximal inerval on which here is no accumulaion poin in he se of updae imes { k } k Z. Noe ha T >, since he number of agens is finie and, for each i V, ɛ i > and x i () = ˆx i () x i () =. The dynamics (4), under he even-riggered communicaion scheme (), has a unique soluion in he ime inerval [, T ). Our firs sep is o show ha he rajecory says bounded in [, T ). Consider he compac-form represenaion (9) of he algorihm. Given i= vi () = and (), for R, q ()=, q :N ()=q :N ()e αț z ()=z ()e α. (5) Hence hese variables are bounded. To bound z :N (), we look ino he soluion of (9d) by subsiuing q :N () = q :N ()e α and considering ( x, Rṙ) as exogenous inpus. In he ime inerval [, T ) ha his soluion exiss, one has z :N () = Φ(, )z :N () Φ(, τ) e ατ q :N ()dτ Φ(, τ)(βr L τ x(τ) R ṙ(τ))dτ where Φ(, τ) = e βr L R( τ). Given he evenriggered communicaion law (), we have x ɛ. Then, for [, T ) we have z :N () ρe βˆλ σ z :N () (6) γ β L ɛ ρ ( e βˆλσ ) ρ q :N () βˆλ { σ (e βˆλ σ e α ), if βˆλ α βˆλ σ α, σ e βˆλ σ, if βˆλ σ = α. Here, we also used R ṙ = Πṙ γ and (). Taking he maximum of each erm in he righhand side of (6) and using (), we have z :N () η, [, T ), (7) were he consan η is given in he saemen. Given ha he number of he agens is finie, { k } k Z is free of Zeno if boh {s k } k Z and { i k } k Z, i {,..., n}, are uniformly lower bounded. Since he former fac holds by assumpion, we nex esablish he laer. Tha is, for each agen, we esablish a lower bound on he iner-execuion imes of () by deermining a lower bound on he amoun of ime i akes for agen i V o have ˆx i x i evolve from o ɛ i. Noe ha y() = z() η r (x() r()), (8) for all [, T ), where we have used (), (5) and (7). On he oher hand, since r w() = for all R by (5), we have Π N w() = w(). Muliplying he second equaion in (8) by R and using RR = Π N, we obain w() = Rq :N () αrz :N (). Using now R R = I N, we deduce w() q :N () α z :N () απ N (x() r()) v() αη, for all [, T ), where we have again used (), (5) and (7). Nex noice ha, for each i V, using (4) and (5), we have d d ˆxi x i = (ˆxi x i ) ẋ i ˆx i x i ẋ i = ṙ i α(x i r i ) β a ij(ˆx i ˆx j ) v i j= = ṙ i α(x i N j= rj ) β a ij(ŷ i ŷ j ) w i j= κ i α ˆx i x i α ŷ i β a ij( ŷ i ŷ j ) w i. j= Then, by noing ha y i y and w i w, and using he bounds esablished above on y and w, we obain d d ˆxi x i α ˆx i x i c i, (9) where c i is given in (4). Nex, using he Comparison Lemma, cf. Khalil [], and ˆx i = x i ( i k ), we deduce ˆx i x i () ci i α (eα( k ) ), i k. () Therefore, he ime i akes ˆx i x i o reach ɛ i is lower bounded by τ i > as given in (3). This fac also implies ha T =. To see his, we reason by conradicion, i.e., suppose T <. Then, he sequence of imes { k } k Z has an accumulaion poin a T. Because here is a finie number of agens, his implies ha here is an agen i V for which { i k } k Z has an accumulaion poin a T. This implies ha i ransmis infiniely ofen in he ime inerval [T, T ) for any (, T ]. Given ha he swiching imes of he digraph are uniformly lower bounded, his conradics he fac ha he iner-even 5

6 imes of () are lower bounded by τ i > on [, T ). Finally, noe ha x i N j= rj x r = y = z, () for i V, so he bound (8) holds for all R. From (5), lim z () = wih an exponenial convergence rae α. From (6), lim sup z :N () (γ β L ɛ )/(βˆλ ) wih an exponenial rae of convergence min{α, βˆλ }. Combining hese facs wih () yields (), and his concludes he proof. No surprisingly, he ulimae convergence error bound () is worse han he bound (3) obained when agens communicae coninuously. The rigger () does no use he full agen sae and hence can be seen as an oupu feedback even-riggered conroller, see e.g., [Donkers and Heemels, ], for which guaraneeing he exisence of lower-bounded iner-execuion imes is in general difficul. Remark 4.4 (Iner-even imes as a funcion of he design parameers): The lower bound τ i in (3) along wih knowledge of he lower bound on he swiching imes of he opology allows o compue bounds on he maximum number of communicaion rounds (and associaed energy spen) by each agen i V (and hence he nework) during a given ime inerval. This lower bound depends on he various ingrediens as follows: τ i is an increasing funcion of ɛ i and a decreasing funcion of α and c i. Through he laer variable, he bound also depends on he graph opology and he design parameer β. Given he definiion of c i, one can deduce ha he faser an inpu of an agen is changing (larger κ i ) or he farher he agen iniially sars from he average of he inpus, he more ofen ha agen would need o rigger communicaion. The connecion beween nework performance and communicaion overhead can also be noed here. Increasing β or decreasing ɛ i o improve he bound () resuls in smaller iner-even imes. Given ha he convergence rae of (4) under () is min{α, βˆλ σ }, decreasing α o increase he iner-even imes slows down convergence. 4.3 Communicaion riggering law for ime-varying conneced undireced graphs Here, we design a disribued even-riggered communicaion law for (4) over neworks wih ime-varying conneced undireced graph ineracion opologies. While he resuls of he previous secion are valid for hese opologies, here we show ha he srucural properies of he Laplacian marix in he undireced case allows he alernaive even-riggered law o have longer inereven imes wih similar racking performance. Proposiion 4.5 (Convergence of (4) over ime-varying conneced undireced graphs wih asynchronous disribued even-riggered communicaion): Assume ha he inpu of each agen i V saisfies ṙ i ess = κ i <, while Π N ṙ ess = γ <. Le he communicaion opology be a conneced, piecewise coninuous ime-varying undireced graph {G} sk wih uniformly bounded weighs. For ɛ R N >, consider an implemenaion of he algorihm (4) over {G} {sk }, where agens communicae according o he model described in Secion 3 wih he sampling imes { i k } k Z of agen i V, saring a i =, deermined by i k = argmax{ [ i k, ) ˆx i () x i () () 4d i ou() j= a ij() ˆx i () ˆx j () ɛ i 4d i ou() }. Then, for any α, β R >, he evoluion saring from x i () R and v i () R wih i= vi () = saisfies lim sup x i () N j= rj () γ (3) ( γ ) ɛ, λ for i V. Furhermore, he iner-execuion imes of agen i V are lower bounded by τ i = ( α ln αɛ ) i, (4) c i di ou where c i is given in (4), wih η subsiued by { ζ = max Π N (x() r()), α Π N (x() r())v() γ ( απ N (x() r())v() γ ) ɛ λ }. Hence, { k } k Z = N i= k Z i k has no accumulaion poin and he execuion of (4) over {G} {sk } is Zeno-free. PROOF. Given an iniial condiion, le [, T ) be he maximal inerval on which here is no accumulaion poin in he se of even imes { k } k Z = N i= k Z i k. The expressions in (5) are equally valid in his case. To bound z :N (), consider V (z :N ) = z :N z :N. The derivaive of V (z :N ) along he rajecories of (9d) can be upper bounded, for [, T ), as V z :N q :N () e α z :Nz :N β(z :NR L Rz :N z :N R L x) γ z :N. For convenience, le s= z :N R L Rz :N z :N R L x. Using RR = Π N, L Π N = Π N L = L, and (8), we obain s = x L x x L x. Given ha he communicaion graph is conneced for all, we can show (see Kia e al. [4b] for deails) s (4d i ou() ˆx i x i a ij() ˆx i ˆx j ), i= j= 6

7 which, ogeher wih () implies ha s ɛ for [, T ) under he even-riggered communicaion law (). Therefore, for θ (, ), we have V ( θ) z :Nz :N θz :Nz :N λ ɛ γ z :N ( z :N q :N () e α ) ( θ)z :Nz :N r, (5) where r = z :N q :N () θz :N z :N λ ɛ γ z :N. Noice ha r < for z :N q :N() γ θ θ ( q:n () θ γ θ ) ɛ λ θ = ˆζ. Hence, for [, T ), as long as z :N () ˆζ, one has V ( θ)z :N z :N. Combining his inequaliy wih he definiion of V, and considering he limiing case θ, we deduce ha, for any z :N () R N and [, T ), z :N () max{ z :N (), ˆζ} = ζ. (6) Nex, following he same argumens as in he proof of Theorem 4.3, one can esablish a lower bound on he iner-execuion imes of any agen i V. To do his, we deermine a lower bound on he ime i akes i o have ˆx i x i evolve from o ɛ i /( di ou ) (noe he conservaiveness in his sep as we disregard he firs erm on he righhand side of ()). The resul of his analysis yields (), bu wih he value of c i given in he saemen (i.e., he expression in (4) wih η subsiued by ζ as defined in (6)). Therefore, he ime i akes ˆx i x i o reach ɛ i /( di ou ) is lower bounded by τ i > as given in (4). This fac also implies T =. Finally, () ogeher wih (5) and (6), imply ha for i V, x i () N j= rj () ζ r (x() r()), for R. Moreover, since T =, from (5), we have V ( θ)z :Nz :N r(), for R, where r() = z :N q :N () e α θz :N z :N λ ɛ γ z :N. Noe ha r() < for z :N q :N() e α γ θ θ ( q:n () e α θ γ θ ) ɛ λ θ = ζ(). Therefore, for R, as long as z :N () ζ(), V ( θ)z :N z :N. x i N j= uj Agens Fig.. Execuions of (4) wih he even-riggered communicaion law described in Secion 3 whose sampling rule is () and of () wih coninuous-ime communicaion. The nework is a weigh-balanced ime-varying digraph of 5 agens wih uni weighs, where for [, 5) i is a fixed ring digraph, for [5, 5) every seconds i is a single conneced pair of nodes and hen for [5, ) i is a ring digraph again. The inpus are u () =.5 sin(.8), u () =.5 sin(.7).5 cos(.6), u 3 () = sin(.), u 4 ()=aan(.5), u 5 ()=. cos(). The op plo shows he racking error and he boom one shows he communicaion imes of each agen. Black dashed (resp. gray solid) lines correspond o he even-riggered sraegy wih he sampling law () wih ɛ i =. (resp. coninuous-ime communicaion ()). In boh cases, α = and β = 4. For each agen, indicaes a sampling (and broadcas) ime and indicaes a ime when he agen acquires a new in-neighbor. As a resul, lim sup z :N () γ ( γ ) ɛ. λ Here, we used θ. On he oher hand, lim z () =. Combining hese facs wih () yields (3), concluding he proof. We should poin ou ha, he guaraneed lower bound (4) on he iner-even-imes is more conservaive han sricly necessary. This is because in our developmen, in order o decouple he analysis of he lower bound on he iner-even imes of each agen from is ou-neighbors, we have negleced he effec of he erm 4d i ou () j= a ij() ˆx i () ˆx j () in (). The simulaions of Secion 5 show he implemenaion of () resuling in iner-even imes longer han he ones of he even-riggered law (). 7

8 5 Simulaions In his secion, we illusrae he performance of he coordinaion algorihm (4) under he even-riggered communicaion laws () and () over a recurrenly joinly srongly conneced digraph (cf. Fig. ), a ring graph (cf. Fig. ) and a ime-varying conneced graph (cf. Fig. 3). Figure shows a small degradaion beween he racking performance of he algorihm (4) wih he even-riggered communicaion law () and he algorihm () wih coninuous-ime communicaion. In he even-riggered implemenaion, he number of imes ha agens {,, 3, 4, 5} communicae in he ime inerval [, ] is (4, 49, 44, 3, 4), respecively. The large error observed in he ime inerval [5, 5) is expeced as in his ime period every wo seconds only wo agens are communicaing wih each oher. Naurally, hese wo agens end o converge o he average of heir inpus and he res of he agens, being oblivious o he inpus of he oher agens, follow heir own inpu. Figure compares he algorihm (4) wih evenriggered communicaion () and he Euler discreizaions of he algorihm () and he proporional-inegral (PI) dynamic average consensus algorihm proposed in [Freeman e al., 6]. We se he parameers of he PI algorihm so ha is ulimae racking error is similar o ha of (). For he discreizaions, we use he larges possible fixed sepsize δ =.39 for he PI algorihm (beyond his value he algorihm diverges) and we use he sepsize δ =. for he algorihm () (from [Kia e al., 4b], convergence is guaraneed if δ (, min{α, β (d ou max) }), which for his example resuls in δ (,.5)). The number of imes ha agens {,, 3, 4, 5} communicae in he ime inerval [, ] is (39, 4, 4, 4, 39), respecively, when implemening even-riggered communicaion (). This is significanly less han he communicaion used by each agen in he Euler discreizaions of () (/. 66 rounds) and he PI algorihm (/.39 5 rounds). Figure 3 shows he execuion of (4) wih he evenriggered communicaion laws () and () over a ime-varying conneced graph. For each agen i {,, 3, 4, 5}, we choose ɛ i for each law so ha he summand in he righ-hand side of he rigger (ɛ i for (), ɛ i /( di ou ) for ()) amouns o he same quaniy. The plos show similar racking performance for boh algorihms, wih he law () inducing less han half communicaion han (). In fac, he number of imes ha agens {,, 3, 4, 5} communicae in he ime inerval [, ] is (43, 56, 7, 57, 55) under () and (, 3, 33, 6, 3) under (). 6 Conclusions We have sudied he muli-agen dynamic average consensus problem over neworks where iner-agen communicaion akes place a discree ime insans in an oppor- x i N j= uj x i N j= uj Agens Fig.. Comparison beween he algorihm (4) employing he even-riggered communicaion law described in Secion 3 wih sampling rule () and he Euler discreizaions of he algorihm () and he proporional-inegral (PI) dynamic average consensus algorihm proposed in [Freeman e al., 6]. For he firs wo, we se α = and β = 4. For he laer, we se γ = 5, L P = L and L I = 4L. The nework is a conneced ring graph of 5 agens wih uni weighs and he inpus are he same of Figure. In he op plo, black (resp. gray) lines correspond o he even-riggered law () wih ɛ i/( d i ou) =. (resp. he Euler discreizaion of he algorihm () wih fixed sepsize δ =.). The middle plo shows he response of he Euler discreizaion of he PI algorihm wih fixed sepsize δ =.39. The horizonal lines in boh he op and middle plos show he ±.5 error bound for reference. The boom plo shows he communicaion imes of each agen using he even-riggered sraegy. unisic fashion. Our saring poin has been our previously developed coninuous-ime dynamic average consensus algorihm which is known o converge exponenially o a small neighborhood of he nework s inpus average. We have proposed wo differen disribued evenriggered laws ha agens can employ o rigger communicaion wih neighbors, depending on wheher he ineracion opology is described by a weigh-balanced and recurrenly joinly srongly conneced digraph or a ime-varying conneced undireced graph. In boh cases, we have esablished he correcness of he algorihm and showed ha a posiive lower bound on he iner-even imes of each agen exiss, ruling ou he presence of Zeno behavior. Fuure work will be devoed o furher relaxing he conneciviy requiremens on he ineracion opol- 8

9 x i N j= uj Agens Fig. 3. Execuions of (4) wih he sampling rule () and () wih he even-riggered communicaion law described in Secion 3. The nework is a ime-varying graph of 5 agens corresponding o a conneced ring graph wih uni weighs where one edge breaks every 3 seconds. The inpus are u () =.5 sin() /( ), u () =.5 sin() /( ) 4, u 3 ()=.5 sin() /( ) 3 5, u 4 ()=.5 sin() e 4, u 5 () =.5 sin() aan().5. The op plo shows he racking error wih he gray solid (resp. black dashed) lines correspond o he law () wih ɛ i =. (resp. () wih ɛ i/( d i ou) =.). In boh cases, α = β =. The boom plo shows he communicaion imes of each agen wih o (resp. ) markers corresponding o he law () (resp. ()). For boh cases, shows he broadcass associaed o he acquisiion of new in-neighbors. ogy (ying hem in o he evoluion of he dynamic inpus available o he agens), he improvemen of he pracical convergence guaranees using ime-varying hresholds in he rigger design, he use of agen absracions in he developmen of self-riggered communicaion laws, and he synhesis of oher disribued riggers ha individual agens can evaluae auonomously and lead o a more efficien use of he limied nework resources. References H. Bai, R. A. Freeman, and K. M. Lynch. Robus dynamic average consensus of ime-varying inpus. In IEEE In. Conf. on Decision and Conrol, pages 34 39, Alana, GA, USA, December. F. Bullo, J. Corés, and S. Marínez. Disribued Conrol of Roboic Neworks. Applied Mahemaics Series. Princeon Universiy Press, 9. ISBN Available a hp:// A. Carron, M. Todescao, R. Carli, and L. Schenao. Adapive consensus-based algorihms for fas esimaion from relaive measuremens. In IFAC Workshop on Disribued Esimaion and Conrol in Neworked Sysems, Koblenz, Germany, March 3. D. V. Dimarogonas, E. Frazzoli, and K. H. Johansson. Disribued even-riggered conrol for muli-agen sysems. IEEE Transacions on Auomaic Conrol, 57(5):9 97,. M. C. F. Donkers and W. P. M. H. Heemels. Oupu-based even-riggered conrol wih guaraneed L -gain and improved and decenralized even-riggering. IEEE Transacions on Auomaic Conrol, 57(6):36 376,. Y. Fan, G. Feng, Y Wang, and C. Song. Disribued evenriggered conrol of muli-agen sysems wih combinaional measuremens. Auomaica, 49():67 675, 3. R. A. Freeman, P. Yang, and K. M. Lynch. Sabiliy and convergence properies of dynamic average consensus esimaors. In IEEE In. Conf. on Decision and Conrol, pages , 6. E. Garcia, Y. Cao, H. Yuc, P. Ansaklis, and D. Casbeer. Decenralised even-riggered cooperaive conrol wih limied communicaion. Inernaional Journal of Conrol, 86 (9): , 3. W. P. M. H. Heemels, K.H Johansson, and P. Tabuada. An inroducion o even-riggered and self-riggered conrol. In IEEE In. Conf. on Decision and Conrol, pages , Maui, HI,. H. Khalil. Nonlinear Sysems. Prenice Hall,. S. S. Kia, J. Corés, and S. Marínez. Dynamic average consensus wih disribued even-riggered communicaion. In IEEE In. Conf. on Decision and Conrol, Los Angeles, CA, 4a. Submied. S. S. Kia, J. Corés, and S. Marínez. Dynamic average consensus under limied conrol auhoriy and privacy requiremens. Inernaional Journal on Robus and Nonlinear Conrol, 4b. To appear. M. Mazo and P. Tabuada. Decenralized even-riggered conrol over wireless sensor/acuaor neworks. IEEE Transacions on Auomaic Conrol, 56():456 46,. C. Nowzari and J. Corés. Zeno-free, disribued evenriggered communicaion and conrol for muli-agen average consensus. In American Conrol Conference, Porland, OR, 4. To appear. R. Olfai-Saber. Disribued Kalman filering for sensor neworks. In IEEE In. Conf. on Decision and Conrol, pages , New Orleans, LA, December 7. R. Olfai-Saber and J. S. Shamma. Consensus filers for sensor neworks and disribued sensor fusion. In IEEE In. Conf. on Decision and Conrol and European Conrol Conference, pages , Seville, Spain, December 5. G. S. Seyboh, D. V. Dimarogonas, and K. H. Johansson. Even-based broadcasing for muli-agen average consensus. Auomaica, 49():45 5, 3. D. P. Spanos, R. Olfai-Saber, and R. M. Murray. Dynamic consensus on mobile neworks. In IFAC World Congress, Prague, Czech Republic, July 5. X. Wang and M. D. Lemmon. Even-riggering in disribued neworked conrol sysems. IEEE Transacions on Auomaic Conrol, 56(3):586 6,. P. Yang, R. A. Freeman, and K. M. Lynch. Disribued cooperaive acive sensing using consensus filers. In IEEE In. Conf. on Roboics and Auomaion, pages 45 4, Roma, Ialy, April 7. P. Yang, R. A. Freeman, and K. M. Lynch. Muli-agen coordinaion by decenralized esimaion and conrol. IEEE Transacions on Auomaic Conrol, 53():48? 496, 8. M. Zhu and S. Marínez. Discree-ime dynamic average consensus. Auomaica, 46():3 39,. 9

Distributed event-triggered communication for dynamic average consensus in networked systems

Distributed event-triggered communication for dynamic average consensus in networked systems Disribued even-riggered communicaion for dynamic average consensus in neworked sysems Solmaz S. Kia a Jorge Corés b Sonia Marínez b a Deparmen of Mechanical and Aerospace Engineering, Universiy of California

More information

Distributed event-triggered communication for dynamic averageconsensusinnetworkedsystems

Distributed event-triggered communication for dynamic averageconsensusinnetworkedsystems Disribued even-riggered communicaion for dynamic averageconsensusinneworkedsysems SolmazS.Kia a JorgeCorés a SoniaMarínez a a Deparmen of Mechanical and Aerospace Engineering, Universiy of California a

More information

Saturation-tolerant average consensus with controllable rates of convergence

Saturation-tolerant average consensus with controllable rates of convergence Sauraion-oleran average consensus wih conrollable raes of convergence Solmaz S. Kia Jorge Cores Sonia Marinez Absrac This paper considers he saic average consensus problem for a muli-agen sysem and proposes

More information

Singularly perturbed algorithms for dynamic average consensus

Singularly perturbed algorithms for dynamic average consensus Singularly perurbed algorihms for dynamic average consensus Solmaz S. Kia Jorge Corés Sonia Marínez Absrac This paper proposes wo coninuous-ime dynamic average consensus algorihms for neworks wih srongly

More information

arxiv: v1 [cs.dc] 8 Mar 2012

arxiv: v1 [cs.dc] 8 Mar 2012 Consensus on Moving Neighborhood Model of Peerson Graph arxiv:.9v [cs.dc] 8 Mar Hannah Arend and Jorgensen Jos Absrac In his paper, we sudy he consensus problem of muliple agens on a kind of famous graph,

More information

Adaptation and Synchronization over a Network: stabilization without a reference model

Adaptation and Synchronization over a Network: stabilization without a reference model Adapaion and Synchronizaion over a Nework: sabilizaion wihou a reference model Travis E. Gibson (gibson@mi.edu) Harvard Medical School Deparmen of Pahology, Brigham and Women s Hospial 55 h Conference

More information

Vehicle Arrival Models : Headway

Vehicle Arrival Models : Headway Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where

More information

Distributed convex optimization via continuous-time coordinationalgorithmswithdiscrete-timecommunication

Distributed convex optimization via continuous-time coordinationalgorithmswithdiscrete-timecommunication Disribued convex opimizaion via coninuous-ime coordinaionalgorihmswihdiscree-imecommunicaion SolmazS.Kia a JorgeCorés a SoniaMarínez a a Deparmen of Mechanical and Aerospace Engineering, Universiy of California

More information

Lecture 20: Riccati Equations and Least Squares Feedback Control

Lecture 20: Riccati Equations and Least Squares Feedback Control 34-5 LINEAR SYSTEMS Lecure : Riccai Equaions and Leas Squares Feedback Conrol 5.6.4 Sae Feedback via Riccai Equaions A recursive approach in generaing he marix-valued funcion W ( ) equaion for i for he

More information

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,

More information

Chapter 2. First Order Scalar Equations

Chapter 2. First Order Scalar Equations Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.

More information

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3 and d = c b - b c c d = c b - b c c This process is coninued unil he nh row has been compleed. The complee array of coefficiens is riangular. Noe ha in developing he array an enire row may be divided or

More information

Stability and Bifurcation in a Neural Network Model with Two Delays

Stability and Bifurcation in a Neural Network Model with Two Delays Inernaional Mahemaical Forum, Vol. 6, 11, no. 35, 175-1731 Sabiliy and Bifurcaion in a Neural Nework Model wih Two Delays GuangPing Hu and XiaoLing Li School of Mahemaics and Physics, Nanjing Universiy

More information

Expert Advice for Amateurs

Expert Advice for Amateurs Exper Advice for Amaeurs Ernes K. Lai Online Appendix - Exisence of Equilibria The analysis in his secion is performed under more general payoff funcions. Wihou aking an explici form, he payoffs of he

More information

10. State Space Methods

10. State Space Methods . Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he

More information

2.160 System Identification, Estimation, and Learning. Lecture Notes No. 8. March 6, 2006

2.160 System Identification, Estimation, and Learning. Lecture Notes No. 8. March 6, 2006 2.160 Sysem Idenificaion, Esimaion, and Learning Lecure Noes No. 8 March 6, 2006 4.9 Eended Kalman Filer In many pracical problems, he process dynamics are nonlinear. w Process Dynamics v y u Model (Linearized)

More information

The Asymptotic Behavior of Nonoscillatory Solutions of Some Nonlinear Dynamic Equations on Time Scales

The Asymptotic Behavior of Nonoscillatory Solutions of Some Nonlinear Dynamic Equations on Time Scales Advances in Dynamical Sysems and Applicaions. ISSN 0973-5321 Volume 1 Number 1 (2006, pp. 103 112 c Research India Publicaions hp://www.ripublicaion.com/adsa.hm The Asympoic Behavior of Nonoscillaory Soluions

More information

STATE-SPACE MODELLING. A mass balance across the tank gives:

STATE-SPACE MODELLING. A mass balance across the tank gives: B. Lennox and N.F. Thornhill, 9, Sae Space Modelling, IChemE Process Managemen and Conrol Subjec Group Newsleer STE-SPACE MODELLING Inroducion: Over he pas decade or so here has been an ever increasing

More information

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon 3..3 INRODUCION O DYNAMIC OPIMIZAION: DISCREE IME PROBLEMS A. he Hamilonian and Firs-Order Condiions in a Finie ime Horizon Define a new funcion, he Hamilonian funcion, H. H he change in he oal value of

More information

di Bernardo, M. (1995). A purely adaptive controller to synchronize and control chaotic systems.

di Bernardo, M. (1995). A purely adaptive controller to synchronize and control chaotic systems. di ernardo, M. (995). A purely adapive conroller o synchronize and conrol chaoic sysems. hps://doi.org/.6/375-96(96)8-x Early version, also known as pre-prin Link o published version (if available):.6/375-96(96)8-x

More information

Notes on Kalman Filtering

Notes on Kalman Filtering Noes on Kalman Filering Brian Borchers and Rick Aser November 7, Inroducion Daa Assimilaion is he problem of merging model predicions wih acual measuremens of a sysem o produce an opimal esimae of he curren

More information

Lecture 4 Notes (Little s Theorem)

Lecture 4 Notes (Little s Theorem) Lecure 4 Noes (Lile s Theorem) This lecure concerns one of he mos imporan (and simples) heorems in Queuing Theory, Lile s Theorem. More informaion can be found in he course book, Bersekas & Gallagher,

More information

EXERCISES FOR SECTION 1.5

EXERCISES FOR SECTION 1.5 1.5 Exisence and Uniqueness of Soluions 43 20. 1 v c 21. 1 v c 1 2 4 6 8 10 1 2 2 4 6 8 10 Graph of approximae soluion obained using Euler s mehod wih = 0.1. Graph of approximae soluion obained using Euler

More information

Guest Lectures for Dr. MacFarlane s EE3350 Part Deux

Guest Lectures for Dr. MacFarlane s EE3350 Part Deux Gues Lecures for Dr. MacFarlane s EE3350 Par Deux Michael Plane Mon., 08-30-2010 Wrie name in corner. Poin ou his is a review, so I will go faser. Remind hem o go lisen o online lecure abou geing an A

More information

Course Notes for EE227C (Spring 2018): Convex Optimization and Approximation

Course Notes for EE227C (Spring 2018): Convex Optimization and Approximation Course Noes for EE7C Spring 018: Convex Opimizaion and Approximaion Insrucor: Moriz Hard Email: hard+ee7c@berkeley.edu Graduae Insrucor: Max Simchowiz Email: msimchow+ee7c@berkeley.edu Ocober 15, 018 3

More information

SOLUTIONS TO ECE 3084

SOLUTIONS TO ECE 3084 SOLUTIONS TO ECE 384 PROBLEM 2.. For each sysem below, specify wheher or no i is: (i) memoryless; (ii) causal; (iii) inverible; (iv) linear; (v) ime invarian; Explain your reasoning. If he propery is no

More information

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004 ODEs II, Lecure : Homogeneous Linear Sysems - I Mike Raugh March 8, 4 Inroducion. In he firs lecure we discussed a sysem of linear ODEs for modeling he excreion of lead from he human body, saw how o ransform

More information

On-line Adaptive Optimal Timing Control of Switched Systems

On-line Adaptive Optimal Timing Control of Switched Systems On-line Adapive Opimal Timing Conrol of Swiched Sysems X.C. Ding, Y. Wardi and M. Egersed Absrac In his paper we consider he problem of opimizing over he swiching imes for a muli-modal dynamic sysem when

More information

A DELAY-DEPENDENT STABILITY CRITERIA FOR T-S FUZZY SYSTEM WITH TIME-DELAYS

A DELAY-DEPENDENT STABILITY CRITERIA FOR T-S FUZZY SYSTEM WITH TIME-DELAYS A DELAY-DEPENDENT STABILITY CRITERIA FOR T-S FUZZY SYSTEM WITH TIME-DELAYS Xinping Guan ;1 Fenglei Li Cailian Chen Insiue of Elecrical Engineering, Yanshan Universiy, Qinhuangdao, 066004, China. Deparmen

More information

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes Represening Periodic Funcions by Fourier Series 3. Inroducion In his Secion we show how a periodic funcion can be expressed as a series of sines and cosines. We begin by obaining some sandard inegrals

More information

Class Meeting # 10: Introduction to the Wave Equation

Class Meeting # 10: Introduction to the Wave Equation MATH 8.5 COURSE NOTES - CLASS MEETING # 0 8.5 Inroducion o PDEs, Fall 0 Professor: Jared Speck Class Meeing # 0: Inroducion o he Wave Equaion. Wha is he wave equaion? The sandard wave equaion for a funcion

More information

12: AUTOREGRESSIVE AND MOVING AVERAGE PROCESSES IN DISCRETE TIME. Σ j =

12: AUTOREGRESSIVE AND MOVING AVERAGE PROCESSES IN DISCRETE TIME. Σ j = 1: AUTOREGRESSIVE AND MOVING AVERAGE PROCESSES IN DISCRETE TIME Moving Averages Recall ha a whie noise process is a series { } = having variance σ. The whie noise process has specral densiy f (λ) = of

More information

Recursive Least-Squares Fixed-Interval Smoother Using Covariance Information based on Innovation Approach in Linear Continuous Stochastic Systems

Recursive Least-Squares Fixed-Interval Smoother Using Covariance Information based on Innovation Approach in Linear Continuous Stochastic Systems 8 Froniers in Signal Processing, Vol. 1, No. 1, July 217 hps://dx.doi.org/1.2266/fsp.217.112 Recursive Leas-Squares Fixed-Inerval Smooher Using Covariance Informaion based on Innovaion Approach in Linear

More information

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Simulaion-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Week Descripion Reading Maerial 2 Compuer Simulaion of Dynamic Models Finie Difference, coninuous saes, discree ime Simple Mehods Euler Trapezoid

More information

ON QUANTIZATION AND COMMUNICATION TOPOLOGIES IN MULTI-VEHICLE RENDEZVOUS 1. Karl Henrik Johansson Alberto Speranzon,2 Sandro Zampieri

ON QUANTIZATION AND COMMUNICATION TOPOLOGIES IN MULTI-VEHICLE RENDEZVOUS 1. Karl Henrik Johansson Alberto Speranzon,2 Sandro Zampieri ON QUANTIZATION AND COMMUNICATION TOPOLOGIES IN MULTI-VEHICLE RENDEZVOUS 1 Karl Henrik Johansson Albero Speranzon, Sandro Zampieri Deparmen of Signals, Sensors and Sysems Royal Insiue of Technology Osquldas

More information

Distributed Linear Supervisory Control

Distributed Linear Supervisory Control 3rd IEEE Conference on Decision and Conrol December -7, Los Angeles, California, USA Disribued Linear Supervisory Conrol Ali Khanafer, Tamer Başar, and Daniel Liberzon Absrac In his work, we propose a

More information

Matrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality

Matrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality Marix Versions of Some Refinemens of he Arihmeic-Geomeric Mean Inequaliy Bao Qi Feng and Andrew Tonge Absrac. We esablish marix versions of refinemens due o Alzer ], Carwrigh and Field 4], and Mercer 5]

More information

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is UNIT IMPULSE RESPONSE, UNIT STEP RESPONSE, STABILITY. Uni impulse funcion (Dirac dela funcion, dela funcion) rigorously defined is no sricly a funcion, bu disribuion (or measure), precise reamen requires

More information

Sliding Mode Extremum Seeking Control for Linear Quadratic Dynamic Game

Sliding Mode Extremum Seeking Control for Linear Quadratic Dynamic Game Sliding Mode Exremum Seeking Conrol for Linear Quadraic Dynamic Game Yaodong Pan and Ümi Özgüner ITS Research Group, AIST Tsukuba Eas Namiki --, Tsukuba-shi,Ibaraki-ken 5-856, Japan e-mail: pan.yaodong@ais.go.jp

More information

The expectation value of the field operator.

The expectation value of the field operator. The expecaion value of he field operaor. Dan Solomon Universiy of Illinois Chicago, IL dsolom@uic.edu June, 04 Absrac. Much of he mahemaical developmen of quanum field heory has been in suppor of deermining

More information

Matlab and Python programming: how to get started

Matlab and Python programming: how to get started Malab and Pyhon programming: how o ge sared Equipping readers he skills o wrie programs o explore complex sysems and discover ineresing paerns from big daa is one of he main goals of his book. In his chaper,

More information

Mean-square Stability Control for Networked Systems with Stochastic Time Delay

Mean-square Stability Control for Networked Systems with Stochastic Time Delay JOURNAL OF SIMULAION VOL. 5 NO. May 7 Mean-square Sabiliy Conrol for Newored Sysems wih Sochasic ime Delay YAO Hejun YUAN Fushun School of Mahemaics and Saisics Anyang Normal Universiy Anyang Henan. 455

More information

Continuous Time. Time-Domain System Analysis. Impulse Response. Impulse Response. Impulse Response. Impulse Response. ( t) + b 0.

Continuous Time. Time-Domain System Analysis. Impulse Response. Impulse Response. Impulse Response. Impulse Response. ( t) + b 0. Time-Domain Sysem Analysis Coninuous Time. J. Robers - All Righs Reserved. Edied by Dr. Rober Akl 1. J. Robers - All Righs Reserved. Edied by Dr. Rober Akl 2 Le a sysem be described by a 2 y ( ) + a 1

More information

Zürich. ETH Master Course: L Autonomous Mobile Robots Localization II

Zürich. ETH Master Course: L Autonomous Mobile Robots Localization II Roland Siegwar Margaria Chli Paul Furgale Marco Huer Marin Rufli Davide Scaramuzza ETH Maser Course: 151-0854-00L Auonomous Mobile Robos Localizaion II ACT and SEE For all do, (predicion updae / ACT),

More information

International Journal of Scientific & Engineering Research, Volume 4, Issue 10, October ISSN

International Journal of Scientific & Engineering Research, Volume 4, Issue 10, October ISSN Inernaional Journal of Scienific & Engineering Research, Volume 4, Issue 10, Ocober-2013 900 FUZZY MEAN RESIDUAL LIFE ORDERING OF FUZZY RANDOM VARIABLES J. EARNEST LAZARUS PIRIYAKUMAR 1, A. YAMUNA 2 1.

More information

6.2 Transforms of Derivatives and Integrals.

6.2 Transforms of Derivatives and Integrals. SEC. 6.2 Transforms of Derivaives and Inegrals. ODEs 2 3 33 39 23. Change of scale. If l( f ()) F(s) and c is any 33 45 APPLICATION OF s-shifting posiive consan, show ha l( f (c)) F(s>c)>c (Hin: In Probs.

More information

2. Nonlinear Conservation Law Equations

2. Nonlinear Conservation Law Equations . Nonlinear Conservaion Law Equaions One of he clear lessons learned over recen years in sudying nonlinear parial differenial equaions is ha i is generally no wise o ry o aack a general class of nonlinear

More information

Linear Response Theory: The connection between QFT and experiments

Linear Response Theory: The connection between QFT and experiments Phys540.nb 39 3 Linear Response Theory: The connecion beween QFT and experimens 3.1. Basic conceps and ideas Q: How do we measure he conduciviy of a meal? A: we firs inroduce a weak elecric field E, and

More information

Optimality Conditions for Unconstrained Problems

Optimality Conditions for Unconstrained Problems 62 CHAPTER 6 Opimaliy Condiions for Unconsrained Problems 1 Unconsrained Opimizaion 11 Exisence Consider he problem of minimizing he funcion f : R n R where f is coninuous on all of R n : P min f(x) x

More information

Predator - Prey Model Trajectories and the nonlinear conservation law

Predator - Prey Model Trajectories and the nonlinear conservation law Predaor - Prey Model Trajecories and he nonlinear conservaion law James K. Peerson Deparmen of Biological Sciences and Deparmen of Mahemaical Sciences Clemson Universiy Ocober 28, 213 Ouline Drawing Trajecories

More information

Chapter 3 Boundary Value Problem

Chapter 3 Boundary Value Problem Chaper 3 Boundary Value Problem A boundary value problem (BVP) is a problem, ypically an ODE or a PDE, which has values assigned on he physical boundary of he domain in which he problem is specified. Le

More information

Research Article Event-Based Control for Average Consensus of Wireless Sensor Networks with Stochastic Communication Noises

Research Article Event-Based Control for Average Consensus of Wireless Sensor Networks with Stochastic Communication Noises Discree Dynamics in Naure and Sociey Volume 2013, Aricle ID 610169, 8 pages hp://dx.doi.org/10.1155/2013/610169 Research Aricle Even-Based Conrol for Average Consensus of Wireless Sensor Neworks wih Sochasic

More information

2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes

2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes Some common engineering funcions 2.7 Inroducion This secion provides a caalogue of some common funcions ofen used in Science and Engineering. These include polynomials, raional funcions, he modulus funcion

More information

Optima and Equilibria for Traffic Flow on a Network

Optima and Equilibria for Traffic Flow on a Network Opima and Equilibria for Traffic Flow on a Nework Albero Bressan Deparmen of Mahemaics, Penn Sae Universiy bressan@mah.psu.edu Albero Bressan (Penn Sae) Opima and equilibria for raffic flow 1 / 1 A Traffic

More information

Finish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details!

Finish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details! MAT 257, Handou 6: Ocober 7-2, 20. I. Assignmen. Finish reading Chaper 2 of Spiva, rereading earlier secions as necessary. handou and fill in some missing deails! II. Higher derivaives. Also, read his

More information

Application of a Stochastic-Fuzzy Approach to Modeling Optimal Discrete Time Dynamical Systems by Using Large Scale Data Processing

Application of a Stochastic-Fuzzy Approach to Modeling Optimal Discrete Time Dynamical Systems by Using Large Scale Data Processing Applicaion of a Sochasic-Fuzzy Approach o Modeling Opimal Discree Time Dynamical Sysems by Using Large Scale Daa Processing AA WALASZE-BABISZEWSA Deparmen of Compuer Engineering Opole Universiy of Technology

More information

SPECTRAL EVOLUTION OF A ONE PARAMETER EXTENSION OF A REAL SYMMETRIC TOEPLITZ MATRIX* William F. Trench. SIAM J. Matrix Anal. Appl. 11 (1990),

SPECTRAL EVOLUTION OF A ONE PARAMETER EXTENSION OF A REAL SYMMETRIC TOEPLITZ MATRIX* William F. Trench. SIAM J. Matrix Anal. Appl. 11 (1990), SPECTRAL EVOLUTION OF A ONE PARAMETER EXTENSION OF A REAL SYMMETRIC TOEPLITZ MATRIX* William F Trench SIAM J Marix Anal Appl 11 (1990), 601-611 Absrac Le T n = ( i j ) n i,j=1 (n 3) be a real symmeric

More information

Multi-Agent Trajectory Tracking with Self-Triggered Cloud Access

Multi-Agent Trajectory Tracking with Self-Triggered Cloud Access 2016 IEEE 55h Conference on Decision and Conrol CDC ARIA Resor & Casino December 12-14, 2016, Las Vegas, USA Muli-Agen Traecory Tracking wih Self-Triggered Cloud Access Anonio Adaldo, Davide Liuzza, Dimos

More information

A Decentralized Second-Order Method with Exact Linear Convergence Rate for Consensus Optimization

A Decentralized Second-Order Method with Exact Linear Convergence Rate for Consensus Optimization 1 A Decenralized Second-Order Mehod wih Exac Linear Convergence Rae for Consensus Opimizaion Aryan Mokhari, Wei Shi, Qing Ling, and Alejandro Ribeiro Absrac This paper considers decenralized consensus

More information

Math 334 Fall 2011 Homework 11 Solutions

Math 334 Fall 2011 Homework 11 Solutions Dec. 2, 2 Mah 334 Fall 2 Homework Soluions Basic Problem. Transform he following iniial value problem ino an iniial value problem for a sysem: u + p()u + q() u g(), u() u, u () v. () Soluion. Le v u. Then

More information

A Primal-Dual Type Algorithm with the O(1/t) Convergence Rate for Large Scale Constrained Convex Programs

A Primal-Dual Type Algorithm with the O(1/t) Convergence Rate for Large Scale Constrained Convex Programs PROC. IEEE CONFERENCE ON DECISION AND CONTROL, 06 A Primal-Dual Type Algorihm wih he O(/) Convergence Rae for Large Scale Consrained Convex Programs Hao Yu and Michael J. Neely Absrac This paper considers

More information

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature On Measuring Pro-Poor Growh 1. On Various Ways of Measuring Pro-Poor Growh: A Shor eview of he Lieraure During he pas en years or so here have been various suggesions concerning he way one should check

More information

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t... Mah 228- Fri Mar 24 5.6 Marix exponenials and linear sysems: The analogy beween firs order sysems of linear differenial equaions (Chaper 5) and scalar linear differenial equaions (Chaper ) is much sronger

More information

Average consensus by gossip algorithms with quantized communication

Average consensus by gossip algorithms with quantized communication Average consensus by gossip algorihms wih quanized communicaion Paolo Frasca, Ruggero Carli, Fabio Fagnani, and Sandro Zampieri Absrac This work sudies how he randomized gossip algorihm can solve he average

More information

An Introduction to Malliavin calculus and its applications

An Introduction to Malliavin calculus and its applications An Inroducion o Malliavin calculus and is applicaions Lecure 5: Smoohness of he densiy and Hörmander s heorem David Nualar Deparmen of Mahemaics Kansas Universiy Universiy of Wyoming Summer School 214

More information

Hybrid Control and Switched Systems. Lecture #3 What can go wrong? Trajectories of hybrid systems

Hybrid Control and Switched Systems. Lecture #3 What can go wrong? Trajectories of hybrid systems Hybrid Conrol and Swiched Sysems Lecure #3 Wha can go wrong? Trajecories of hybrid sysems João P. Hespanha Universiy of California a Sana Barbara Summary 1. Trajecories of hybrid sysems: Soluion o a hybrid

More information

On Boundedness of Q-Learning Iterates for Stochastic Shortest Path Problems

On Boundedness of Q-Learning Iterates for Stochastic Shortest Path Problems MATHEMATICS OF OPERATIONS RESEARCH Vol. 38, No. 2, May 2013, pp. 209 227 ISSN 0364-765X (prin) ISSN 1526-5471 (online) hp://dx.doi.org/10.1287/moor.1120.0562 2013 INFORMS On Boundedness of Q-Learning Ieraes

More information

Designing Information Devices and Systems I Spring 2019 Lecture Notes Note 17

Designing Information Devices and Systems I Spring 2019 Lecture Notes Note 17 EES 16A Designing Informaion Devices and Sysems I Spring 019 Lecure Noes Noe 17 17.1 apaciive ouchscreen In he las noe, we saw ha a capacior consiss of wo pieces on conducive maerial separaed by a nonconducive

More information

Synchronization in Networks of Identical Linear Systems

Synchronization in Networks of Identical Linear Systems Proceedings of he 47h IEEE Conference on Decision and Conrol Cancun, Meico, Dec. 9-11, 28 Synchronizaion in Neworks of Idenical Linear Sysems Luca Scardovi and Rodolphe Sepulchre Absrac The paper invesigaes

More information

control properties under both Gaussian and burst noise conditions. In the ~isappointing in comparison with convolutional code systems designed

control properties under both Gaussian and burst noise conditions. In the ~isappointing in comparison with convolutional code systems designed 535 SOFT-DECSON THRESHOLD DECODNG OF CONVOLUTONAL CODES R.M.F. Goodman*, B.Sc., Ph.D. W.H. Ng*, M.S.E.E. Sunnnary Exising majoriy-decision hreshold decoders have so far been limied o his paper a new mehod

More information

CHAPTER 2 Signals And Spectra

CHAPTER 2 Signals And Spectra CHAPER Signals And Specra Properies of Signals and Noise In communicaion sysems he received waveform is usually caegorized ino he desired par conaining he informaion, and he undesired par. he desired par

More information

LAPLACE TRANSFORM AND TRANSFER FUNCTION

LAPLACE TRANSFORM AND TRANSFER FUNCTION CHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION Professor Dae Ryook Yang Spring 2018 Dep. of Chemical and Biological Engineering 5-1 Road Map of he Lecure V Laplace Transform and Transfer funcions

More information

On Distributed Estimation for Sensor Networks

On Distributed Estimation for Sensor Networks On Disribued Esimaion for Sensor Neworks Albero Speranzon, Carlo Fischione, and Karl Henrik Johansson Absrac Disribued esimaors for sensor neworks are discussed. The considered problem is on how o rack

More information

Two Popular Bayesian Estimators: Particle and Kalman Filters. McGill COMP 765 Sept 14 th, 2017

Two Popular Bayesian Estimators: Particle and Kalman Filters. McGill COMP 765 Sept 14 th, 2017 Two Popular Bayesian Esimaors: Paricle and Kalman Filers McGill COMP 765 Sep 14 h, 2017 1 1 1, dx x Bel x u x P x z P Recall: Bayes Filers,,,,,,, 1 1 1 1 u z u x P u z u x z P Bayes z = observaion u =

More information

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems.

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems. Mah 2250-004 Week 4 April 6-20 secions 7.-7.3 firs order sysems of linear differenial equaions; 7.4 mass-spring sysems. Mon Apr 6 7.-7.2 Sysems of differenial equaions (7.), and he vecor Calculus we need

More information

PENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD

PENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD PENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD HAN XIAO 1. Penalized Leas Squares Lasso solves he following opimizaion problem, ˆβ lasso = arg max β R p+1 1 N y i β 0 N x ij β j β j (1.1) for some 0.

More information

T L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB

T L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB Elecronic Companion EC.1. Proofs of Technical Lemmas and Theorems LEMMA 1. Le C(RB) be he oal cos incurred by he RB policy. Then we have, T L E[C(RB)] 3 E[Z RB ]. (EC.1) Proof of Lemma 1. Using he marginal

More information

Robotics I. April 11, The kinematics of a 3R spatial robot is specified by the Denavit-Hartenberg parameters in Tab. 1.

Robotics I. April 11, The kinematics of a 3R spatial robot is specified by the Denavit-Hartenberg parameters in Tab. 1. Roboics I April 11, 017 Exercise 1 he kinemaics of a 3R spaial robo is specified by he Denavi-Harenberg parameers in ab 1 i α i d i a i θ i 1 π/ L 1 0 1 0 0 L 3 0 0 L 3 3 able 1: able of DH parameers of

More information

EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE

EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE Version April 30, 2004.Submied o CTU Repors. EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE Per Krysl Universiy of California, San Diego La Jolla, California 92093-0085,

More information

Average Number of Lattice Points in a Disk

Average Number of Lattice Points in a Disk Average Number of Laice Poins in a Disk Sujay Jayakar Rober S. Sricharz Absrac The difference beween he number of laice poins in a disk of radius /π and he area of he disk /4π is equal o he error in he

More information

5. Stochastic processes (1)

5. Stochastic processes (1) Lec05.pp S-38.45 - Inroducion o Teleraffic Theory Spring 2005 Conens Basic conceps Poisson process 2 Sochasic processes () Consider some quaniy in a eleraffic (or any) sysem I ypically evolves in ime randomly

More information

Hamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation:

Hamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation: M ah 5 7 Fall 9 L ecure O c. 4, 9 ) Hamilon- J acobi Equaion: Weak S oluion We coninue he sudy of he Hamilon-Jacobi equaion: We have shown ha u + H D u) = R n, ) ; u = g R n { = }. ). In general we canno

More information

Some Basic Information about M-S-D Systems

Some Basic Information about M-S-D Systems Some Basic Informaion abou M-S-D Sysems 1 Inroducion We wan o give some summary of he facs concerning unforced (homogeneous) and forced (non-homogeneous) models for linear oscillaors governed by second-order,

More information

Final Spring 2007

Final Spring 2007 .615 Final Spring 7 Overview The purpose of he final exam is o calculae he MHD β limi in a high-bea oroidal okamak agains he dangerous n = 1 exernal ballooning-kink mode. Effecively, his corresponds o

More information

Reading from Young & Freedman: For this topic, read sections 25.4 & 25.5, the introduction to chapter 26 and sections 26.1 to 26.2 & 26.4.

Reading from Young & Freedman: For this topic, read sections 25.4 & 25.5, the introduction to chapter 26 and sections 26.1 to 26.2 & 26.4. PHY1 Elecriciy Topic 7 (Lecures 1 & 11) Elecric Circuis n his opic, we will cover: 1) Elecromoive Force (EMF) ) Series and parallel resisor combinaions 3) Kirchhoff s rules for circuis 4) Time dependence

More information

Stabilization of NCSs: Asynchronous Partial Transfer Approach

Stabilization of NCSs: Asynchronous Partial Transfer Approach 25 American Conrol Conference June 8-, 25 Porland, OR, USA WeC3 Sabilizaion of NCSs: Asynchronous Parial Transfer Approach Guangming Xie, Long Wang Absrac In his paper, a framework for sabilizaion of neworked

More information

4. Advanced Stability Theory

4. Advanced Stability Theory Applied Nonlinear Conrol Nguyen an ien - 4 4 Advanced Sabiliy heory he objecive of his chaper is o presen sabiliy analysis for non-auonomous sysems 41 Conceps of Sabiliy for Non-Auonomous Sysems Equilibrium

More information

Event-triggered PI control subject to actuator saturation

Event-triggered PI control subject to actuator saturation IFAC Conference on Advances in PID Conrol Brescia (Ialy), March 28-30, 2012 Even-riggered PI conrol subjec o acuaor sauraion D. Lehmann K. H. Johansson Auomaic Conrol Laboraory, KTH Royal Insiue of Technology,

More information

Cash Flow Valuation Mode Lin Discrete Time

Cash Flow Valuation Mode Lin Discrete Time IOSR Journal of Mahemaics (IOSR-JM) e-issn: 2278-5728,p-ISSN: 2319-765X, 6, Issue 6 (May. - Jun. 2013), PP 35-41 Cash Flow Valuaion Mode Lin Discree Time Olayiwola. M. A. and Oni, N. O. Deparmen of Mahemaics

More information

GENERALIZATION OF THE FORMULA OF FAA DI BRUNO FOR A COMPOSITE FUNCTION WITH A VECTOR ARGUMENT

GENERALIZATION OF THE FORMULA OF FAA DI BRUNO FOR A COMPOSITE FUNCTION WITH A VECTOR ARGUMENT Inerna J Mah & Mah Sci Vol 4, No 7 000) 48 49 S0670000970 Hindawi Publishing Corp GENERALIZATION OF THE FORMULA OF FAA DI BRUNO FOR A COMPOSITE FUNCTION WITH A VECTOR ARGUMENT RUMEN L MISHKOV Received

More information

A Shooting Method for A Node Generation Algorithm

A Shooting Method for A Node Generation Algorithm A Shooing Mehod for A Node Generaion Algorihm Hiroaki Nishikawa W.M.Keck Foundaion Laboraory for Compuaional Fluid Dynamics Deparmen of Aerospace Engineering, Universiy of Michigan, Ann Arbor, Michigan

More information

The Arcsine Distribution

The Arcsine Distribution The Arcsine Disribuion Chris H. Rycrof Ocober 6, 006 A common heme of he class has been ha he saisics of single walker are ofen very differen from hose of an ensemble of walkers. On he firs homework, we

More information

Undetermined coefficients for local fractional differential equations

Undetermined coefficients for local fractional differential equations Available online a www.isr-publicaions.com/jmcs J. Mah. Compuer Sci. 16 (2016), 140 146 Research Aricle Undeermined coefficiens for local fracional differenial equaions Roshdi Khalil a,, Mohammed Al Horani

More information

Variational Iteration Method for Solving System of Fractional Order Ordinary Differential Equations

Variational Iteration Method for Solving System of Fractional Order Ordinary Differential Equations IOSR Journal of Mahemaics (IOSR-JM) e-issn: 2278-5728, p-issn: 2319-765X. Volume 1, Issue 6 Ver. II (Nov - Dec. 214), PP 48-54 Variaional Ieraion Mehod for Solving Sysem of Fracional Order Ordinary Differenial

More information

Timed Circuits. Asynchronous Circuit Design. Timing Relationships. A Simple Example. Timed States. Timing Sequences. ({r 6 },t6 = 1.

Timed Circuits. Asynchronous Circuit Design. Timing Relationships. A Simple Example. Timed States. Timing Sequences. ({r 6 },t6 = 1. Timed Circuis Asynchronous Circui Design Chris J. Myers Lecure 7: Timed Circuis Chaper 7 Previous mehods only use limied knowledge of delays. Very robus sysems, bu exremely conservaive. Large funcional

More information

In this chapter the model of free motion under gravity is extended to objects projected at an angle. When you have completed it, you should

In this chapter the model of free motion under gravity is extended to objects projected at an angle. When you have completed it, you should Cambridge Universiy Press 978--36-60033-7 Cambridge Inernaional AS and A Level Mahemaics: Mechanics Coursebook Excerp More Informaion Chaper The moion of projeciles In his chaper he model of free moion

More information

CHAPTER 10 VALIDATION OF TEST WITH ARTIFICAL NEURAL NETWORK

CHAPTER 10 VALIDATION OF TEST WITH ARTIFICAL NEURAL NETWORK 175 CHAPTER 10 VALIDATION OF TEST WITH ARTIFICAL NEURAL NETWORK 10.1 INTRODUCTION Amongs he research work performed, he bes resuls of experimenal work are validaed wih Arificial Neural Nework. From he

More information

Testing for a Single Factor Model in the Multivariate State Space Framework

Testing for a Single Factor Model in the Multivariate State Space Framework esing for a Single Facor Model in he Mulivariae Sae Space Framework Chen C.-Y. M. Chiba and M. Kobayashi Inernaional Graduae School of Social Sciences Yokohama Naional Universiy Japan Faculy of Economics

More information

A Dynamic Model of Economic Fluctuations

A Dynamic Model of Economic Fluctuations CHAPTER 15 A Dynamic Model of Economic Flucuaions Modified for ECON 2204 by Bob Murphy 2016 Worh Publishers, all righs reserved IN THIS CHAPTER, OU WILL LEARN: how o incorporae dynamics ino he AD-AS model

More information

Adaptive Leader-following Consensus of Multiple Uncertain Rigid Spacecraft Systems

Adaptive Leader-following Consensus of Multiple Uncertain Rigid Spacecraft Systems Preprins of he 19h World Congress he Inernaional Federaion of Auomaic Conrol Cape own, Souh Africa. Augus 24-29, 214 Adapive Leader-following Consensus of Muliple Uncerain Rigid Spacecraf Sysems He Cai

More information