Multi-Agent Trajectory Tracking with Self-Triggered Cloud Access

Size: px
Start display at page:

Download "Multi-Agent Trajectory Tracking with Self-Triggered Cloud Access"

Transcription

1 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 V. Dimarogonas, Member, IEEE and Karl H. Johansson, Fellow, IEEE Absrac This paper presens a cloud-suppored conrol algorihm for coordinaed raecory racking of neworked auonomous agens. The moivaing applicaion is he coordinaed conrol of Auonomous Underwaer Vehicles. The conrol obecive is o have he vehicles rack a reference raecory while keeping an assigned formaion. Raher han relying on iner-agen communicaion, which is inerdiced underwaer, coordinaion is achieved by leing he agens inermienly access a shared informaion reposiory hosed on a cloud. An even-based law is proposed o schedule he accesses of each agen o he cloud. We show ha, wih he proposed scheduling of he cloud accesses, he agens achieve he required coordinaion obecive. Numerical simulaions corroborae he heoreical resuls. I. INTRODUCTION Neworked vehicle sysems have araced a noable amoun of research in he pas few decades [1] [3]. In mos applicaions, employing a eam of vehicle agens insead of a single plaform has numerous advanages. For example, a group of agens usually provides robusness wih respec o he failure of a single agen in he group. When sampling a propery in a region of he space, a eam of mobile agens will provide a larger number of samples and increased daa redundancy. Also, cerain asks may be srucurally impossible o perform wih a single agen. However, he use of a flee ineviably brings abou he problem of coordinaing he vehicles. Muli-agen coordinaion is paricularly challenging in he case of Auonomous Underwaer Vehicles AUVs because of heir limied communicaion, sensing and localizaion capabiliies [4], [5]. AUVs have numerous applicaions, including, o name us a few, oceanographic surveys, mine search, inspecion of underwaer srucures, and measuremen of chemical properies in a waer body [6]. Underwaer communicaion may be implemened by means of acousic modems, bu such modems are nooriously expensive, power-hungry and limied in boh radius and bandwidh. Underwaer posiioning is also difficul, since good inerial sensors are very expensive, and acousic posiioning sysems A. Adaldo, D. V. Dimarogonas and K. H. Johansson are wih he deparmen of Auomaic Conrol and ACCESS Linnaeus Cener, School of Elecrical Engineering, KTH Royal Insiue of Technology, Osquldas väg 10, 10044, Sockholm, Sweden; s: {adaldo,dimos,kalle}@kh.se. D. Liuzza is wih Universiy of Sannio in Beneveno, Deparmen of Engineering, Piazza Roma 21, 82100, Beneveno, Ialy; davide.liuzza@unisannio.i. This work has received funding from he European Union Horizon 2020 Research and Innovaion Programme under he Gran Agreemen No , AEOROWORKS, from he Swedish Foundaion for Sraegic Research, from he Swedish Research Council, and from he Knu och Alice Wallenberg foundaion. have a limied range. Moreover, GPS is no available underwaer, and a vehicle has o surface whenever i needs o ge a posiion fix [7]. To deal wih communicaion consrains, even- and self-riggered conrol designs [8] can be applied o neworked muli-agen sysems [9]. In his paper, selfriggered muli-agen conrol is considered in combinaion wih he suppor of a shared informaion reposiory hosed on a cloud. Namely, he cloud is inermienly accessed by he agens according o a self-riggered proocol. Since direc communicaion among he agens is inerdiced when hey are underwaer, he agens only exchange daa hrough he cloud reposiory, which is accessed asynchronously. The moivaing applicaion is a leader-following raecory racking ask for a formaion of AUVs subec o disurbances. In he radiional even- and self-riggered neworked conrol, oneo-one communicaion needs o be esablished a leas when an agen needs o updae is conrol signal. Conversely, wih cloud-based approaches [10] [12], he vehicles exchange informaion wihou opening a communicaion channel beween each oher. Each vehicle simply leaves is informaion on he cloud for he ohers o download laer. A cloudsuppored conrol archiecure for muli-agen coordinaion was proposed by he auhors in [12], where he problem of driving a eam of vehicles o a saic formaion is addressed. In his paper, he approach is furher developed o address muli-agen leader-follower racking problems, under more general nework opologies. Two differen coordinaion obecives, namely pracical and asympoic convergence o a given formaion, are formulaed mahemaically, and graphheoreical resuls are used o show ha he proposed cloudsuppored sraegy achieves he desired obecives, despie only using oudaed informaion. Edge-space analysis [13], [14] is used o address direced nework opologies, hus allowing for leader-follower coordinaion. II. PRELIMINARIES In his paper, denoes he euclidean norm of a vecor or he corresponding induced norm of a marix. Moreover, {A} i,k denoes he enry of A in he i-h row and k-h column, while {A} i denoes he row vecor corresponding o he i-h row of A. The null vecor in R n is denoed as 0 n. The se of he posiive inegers is denoed as N, while N 0 = N {0}. A digraph is a uple, wih = {1,, N} and {, i i,, i }. The elemens of and are called respecively verexes and edges of he graph. A pah from verex o verex i is a sequence of verexes saring wih and ending wih i such ha any wo consecuive verexes in he sequence consiue an edge. A spanning ree is a digraph /16/$ IEEE 2207

2 , wih N 1 edges such ha here exiss a node r wih a pah o any oher node. The node r is called he roo of he spanning ree. A digraph, is said o conain a spanning ree if, is a spanning ree for some subse of. Consider a digraph,, and le he edges be denoed as = {e 1,, e M }. The incidence marix of he digraph is defined as B R N M such ha 1 if e k =, i for some, {B} i,k = 1 if e k = i, for some, 0 oherwise. The in-incidence marix is defined as B R N M such ha {B } i,k = {B} i,k if {B} i,k {0, 1} and {B } i,k = 0 if {B} i,k = 1. For a spanning ree, he edge Laplacian [14] is defined as L e = B B. For a digraph, ha conains a spanning ree, bu is no iself a spanning ree, le =, wih, being a spanning ree. Wihou loss of generaliy, le = {e 1,, e N 1 } and = {e N,, e M }. Pariion he incidence and in-incidence marices accordingly as B = [B, B ] and B = [B, B ] respecively. Then, B has a lef-pseudoinverse B [15], and he reduced edge Laplacian is defined as L r = B B + B B B. 1 For any spanning ree, he edge Laplacian is posiive definie, while for any graph ha conains a spanning ree, bu is no iself a spanning ree, he reduced edge Laplacian is posiive definie [14]. III. PROBLEM SETTING Consider a muli-agen sysem composed of N agens indexed as i {1,, N} =, wih kinemaics described by { x i = u i + d i 0, i, 2 x i 0 = x i,0 i. where x i R n represens he sae of agen i, x i,0 is he iniial sae, u i represens a conrol inpu, and d i represens a disurbance signal. The conrol obecive is o have all he agens follow a desired reference raecory r wihin a given olerance. Such conrol obecive is referred o as pracical consensus o he raecory r, and can be formalized as follows. Definiion III.1. The muli-agen sysem 2 is said o achieve pracical consensus o he reference raecory r wih olerance ε > 0 if lim sup x i r ε for all i. Remark III.1. In erms of our moivaing applicaion, we have n = 2, each agen represens an AUV, and x i + b i R 2 represens a waypoin for vehicle i, where b i is a consan bias vecor. In his way, pracical consensus of x 1,, x N corresponds o convergence of he vehicles o a formaion abou r defined by he bias vecors b 1,, b N. However, he analysis remains valid for any n N. The reference raecory can be measured only by a subse of he agens, which are referred o as he leaders in he muli-agen sysem. In his work, we assume ha he agens canno exchange any direc informaion wih each oher. This models he scenarios where, as in our AUV seup, communicaion among he agens is physically inerdiced. In order o exchange informaion, he agens can only upload and download daa on a shared reposiory hosed on a cloud. Namely, when i is conneced o he cloud, an agen can deposi some informaion, and, a he same ime, download some informaion ha was previously uploaded by some oher agens. On he oher hand, when i is no conneced o he cloud, an agen canno exchange informaion a all. For he purposes of his work, an agen s access o he cloud can be considered an insananeous even, while communicaion proocol problems, such as delays and packe losses, are lef ou of he scope of his work. The cloud is modelled as a shared resource wih limied hroughpu and sorage capaciy, and hanks o he conrol algorihm ha we are going o define, i is accessed only inermienly and asynchronously, and he amoun of daa sored herein does no grow over ime. For each agen i, we define he sequence { i,k } k N0 of he agen s accesses o he cloud. Namely, i,k wih k N denoes he ime when agen i accesses he cloud for he kh ime, while convenionally i,0 = 0 for all i. When an agen i accesses he cloud a ime i,k, i also riggers a measuremen of is curren sae, which we denoe as x i,k : x i,k = x i,k. 3 If agen i is a leader, i also produces a measuremen of he curren value of he reference raecory, which we denoe as r i,k : r i,k = r i,k. 4 The disurbance signals and he reference raecory saisfy he following assumpion. Assumpion III.1. The disurbance signals d i in 2 and he reference raecory r saisfy d i δ i and r δ 0, where δ i = δ i,0 δ i, e λ δ +δ i,, i {0, 1,, N}, and δ i,0, δ i,, λ δ, are known non-negaive consans for all i {0, 1,, N}. Our goal is o propose a conrol sraegy such ha each agen uses he informaion acquired from he cloud o aain pracical consensus as by Definiion III.1. To compleely specify he conrol sraegy, we need o define: he conrol signals for each agen, he informaion ha is uploaded and downloaded by each agen when accessing he cloud, and a law for scheduling he cloud accesses. Firs, we define he conrol signals u i. In he proposed conrol sraegy, each signal u i is piecewise consan, and i is updaed upon agen s i cloud accesses, i.e., u i = u i,k [ i,k, i,k

3 AGENT TIME POSITION CONTROL NEXT 1 1,l1 x 1,l1 u 1,l1 1,l ,l2 x 2,l2 u 2,l2 2,l2 +1 N N,lN x N,lN u N,lN N,lN +1 Table III. Schemaic represenaion of he daa sored in he cloud a a generic ime insan 0. Namely, he conrol signals are compued as follows: u i,k = c p i r i,k x i,k + x i,k x i,k, 6 i where c > 0 is a conrol gain, p i = 1 if i and p i = 0 oherwise, i {i}, and x i,k is an esimae of he sae of agen done by agen i a ime i,k. Such esimae is defined laer in his secion. Nex, we define he informaion uploaded and downloaded by each agen when accessing he cloud. When agen i accesses he cloud a ime i,k, i uploads: he curren ime i,k, he measuremen of is curren sae x i,k, he conrol signal u i,k ha is going o be applied unil he following access, and he scheduled ime i,k+1 of he following access. When hese values are uploaded, hey overwrie hose ha were uploaded by he same agen upon he previous access. In his way, he amoun of daa conained in he cloud remains consan. Namely, for each agen, he cloud conains he informaion ha was uploaded upon ha agen s mos recen access. Denoing as l i he index of he mos recen access of agen i before ime, i.e. l i = max{k N 0 i,k }, a abular represenaion of he daa conained in he cloud a a generic ime is given in Table III. Before uploading is own informaion, agen i downloads and sores he informaion corresponding o he agens i. Such informaion is used by agen i o consruc he esimaes x i,k for i ha are used for compuing he conrol signal 6, and also o schedule is following access o he cloud. Namely, he esimae x i,k is compued as follows: x i,k = x,l i,k + u,l i,k i,k,l i,k. 7 Noe ha such esimae coincides wih he sae ha agen would have a ime i,k if no disurbances were acing on i in he ime inerval [,l i,k, i,k. Finally, le us define he rule for scheduling he agens accesses o he cloud. Each agen schedules is own accesses recursively, according o he following rule: i,k+1 = inf{ > i,k Δ i,k ζ i,k σ i,k ς i }, 8 Δ i,k = δ i,k i τdτ, 9 { } ςq ζ i,k = min, 10 q i q 2c q σ i,k =c i + p i u i,k i,k u,h min{,,h +1} i,k i + i + p i Δ i,k + p i δ 0 τdτ i,k + Δ,h + μ τdτ i,h +1 i >,h +1, 11 ς i =ς i,0 ς i, e λ ς +ς i,, 12 where ς i,0, ς i, and λ ς are given non-negaive consans for all i, h = l i,k, δ i and ρ are defined in Assumpion III.1, and μ τ is a bounded scalar signal o be given laer in he paper. The expression of C i,k emerges from he analysis conduced in he following secion, and herefore, will be clarified laer. Noe however ha funcions 9 12 can be compued locally by agen i a ime i,k by using he informaion acquired from he cloud a ha ime, and knowing q for q i q. The funcions ς i wih i are referred o as hreshold funcions, since ς i a hreshold ha σ i,k mus overcome o rigger he cloud access i,k+1 of agen i. Remark III.2. In erms of our moivaing applicaion, an agen s accesses he cloud correspond o he imes when an AUV comes o he waer surface. A posiion measuremen corresponds o GPS fix ha a vehicle can obain while on he waer surface. On he oher hand, when a vehicle is underwaer, i canno communicae wih oher vehicles or access GPS. Neverheless, i has o find a conrol value and he nex surfacing insan coping wih he fac ha in he fuure oher vehicles will surface and updae heir conrol inpu o a ye unknown value. Remark III.3. The fundamenal difference beween he proposed conrol sraegy and he maoriy of he exising selfriggered coordinaion sraegies for muli-agen sysems is ha, in he proposed sraegy, an agen does no require oher agens o exchange informaion when i needs o updae is conrol signal. Conversely, when an agen needs o updae is conrol signal, i uses he informaion ha is already available in he cloud, i.e., ha was previously uploaded by he oher agens upon heir own access imes. IV. MAIN RESULT In his secion, we show how he muli-agen sysem 2, under he conrol algorihm defined by 5 12, can achieve pracical consensus o he reference raecory as by Definiion III.1. Firs, we need o inroduce a digraph induced by he ses and i wih i ha capures he opology of he informaion exchanges processed hrough he cloud. Definiion IV.1. Consider he muli-agen sysem 2 under he conrol law 6. Le = {0} and = {, i i, i } {0, i i }. We say ha he digraph, is he digraph associaed wih he muli-agen sysem. Moreover, we denoe he edges of he digraph as = {e 1,, e M }. Noe ha, i for some i, if and only if agen i downloads he informaion uploaded by agen 2209

4 , while 0, i if and only if agen i is a leader, i.e., if i receives informaion abou he reference raecory. Therefore, he digraph, represens he opology of he informaion exchanges ha are processed hrough he cloud. The following assumpion ensures ha he informaion abou he reference raecory can reach all he agens in he sysem. Assumpion IV.1. The digraph, associaed wih he muli-agen sysem 2 conains a spanning ree wih roo in he verex 0. Namely, we le, wihou loss of generaliy, =, where = {e 1,, e N 1 }, = {e N,, e M }, and, is a spanning ree wih roo in he verex 0. For each edge e l =, i, wih l {1,, M}, we le y l = x x i if, and y l = r x i if = 0. In oher words, y l is he mismach beween he saes of he wo agens whose indexes appear in he edge e l. Le y = [y 1,, y N 1 ], y = [y N,, y M ], and y = [y, y ]. 13 Le B and B be respecively he incidence and inincidence marices of,, and le hem be pariioned as B = [B, B ] and B = [B, B ] according o how is pariioned ino and. Noe ha, leing x = [r, x 1,, x N ], we have y = B I nx 14 y = B I n x. 15 Under Assumpion IV.1, B has a lef-pseudoinverse B see [15] for furher deails, herefore, from 14 and 15, i follows ha y = B B I n y. 16 Finally, le he reduced edge Laplacian L r of, be defined as in Secion II. If, is iself a spanning ree, le y = y, B = B and L r = L e, where L e is also defined in Secion II. Nex, le us inroduce some signals which shall be used in he convergence analysis. Consider he signals v i = c p i r x i + x x i, 17 i for all i. Noe ha v i can be obained from 6 by subsiuing he measuremens r i,k, x i,k and he esimaes x i,k respecively wih r i, x i and x. Le v = [0 n, v 1,, v N ] so ha we can rewrie 17 compacly as v = cb I n y. 18 Le u i be he mismach beween he acual conrol inpu u i and v i for each i, i.e., u i = u i v i, 19 and le u = [0 n, u 1,, u N ]. We are now in he posiion o sae our firs convergence resul. Theorem IV.1. Consider he muli-agen sysem 2, and le Assumpions III.1 and IV.1 hold. If u i ς i for all [0, and all i, hen here exis α, λ > 0 such ha y η for all [0,, where η =α η 0 e cλ + B e cλ τ δτ + ςτ dτ, 0 20 η 0 = y 0, δ = [δ 0, δ 1,, δ N ] and ς = [0, ς 1,, ς N ]. Proof. Subsiuing 19 ino 2, we have x i = v i + u i + d i. 21 Leing d = [ r, d 1,, d N ], 21 can be rewrien compacly as x = v + u + d. 22 Lef-muliplying boh sides of 22 by B I n, and using 14 and 18 and he properies of he Kronecker produc, we have y = cb B I n y B I n u + d, 23 Subsiuing 13 ino 23, observing ha 16 holds hanks o Assumpion IV.1, and using he properies of he Kronecker produc, we have y = cl r I n y B I n u + d, 24 where L r is he reduced edge Laplacian of,, as defined in 1. The Laplace soluion of 24 reads y = e cl y 0 e cl τ B u + ddτ, 0 25 where we have denoed L = L r I n and B = B I n for breviy. Taking norms of boh sides in 25, and using he riangular inequaliy, he properies of he Kronecker produc, Assumpion III.1, and he hypohesis u i ς i for all [0, and all i, we have y e cl y 0 + B e cl τ δτ + ςτ dτ, 0 26 for all [0,, where δ and ς are defined in he heorem saemen. Since L r is posiive definie, L = L r I n is Hurwiz, and herefore here exis α, λ > 0 such ha e cl α e cλ The proof is concluded by subsiuing 27 ino 26. Remark IV.1. The posiive scalar λ mus be smaller han min{eigl r }, bu can be chosen as close o ha as desired. If L r is diagonalizable, one can choose λ = min{eigl r } and α = V V 1, where L r = V ΛV 1 and Λ is diagonal [16]. Corollary IV.1. Under he same hypoheses as Theorem IV.1, we have u i μ i, for all [0, and all i, where μ i = β i η + ς i,

5 β i = c {B + B B B } i, and η is defined in 20. Proof. Subsiuing 13 ino 18, and using 16, we have v i = c{b + B B B } i I n y. Taking norms of boh sides, and using he Cauchy-Swarz inequaliy, we have v i β i y. 29 From 19, using he riangular inequaliy, we have u i v i + u i. 30 Using 29 and he hypohesis u i ς i ino 30 concludes he proof. The nex sep in our analysis is o show ha he condiion u i ς i holds for all 0 and all i if he conrol algorihm defined by 5 12 is applied wih η given by 20. This is formalized in he following heorem. Theorem IV.2. Consider he muli-agen sysem 2, le Assumpions III.1 and IV.1 hold, and le he sysem be conrolled by he algorihm defined by 5 12, 20 and 28. Then he closed-loop sysem does no exhibi Zeno behavior and u i ς i holds for all 0 and all i. Proof. Subsiuing 6 and 17 ino 19, we have u i =c p i r i,k x i,k + x i,k x i,k i p i r x i x x i i 31 for [ i,k, i,k+1, where x i,k, r i,k and x i,k are defined in 3, 4 and 7 respecively. Reordering he erms in 31, we have u i =c i + p i x i x i,k cp i r r i,k c x x i,k 32 i for [ i,k, i,k+1. Firs consider he erm x i x i,k in 32. Inegraing 2 in [ i,k,, and using 3, 5 and 6, we have x i x i,k = u i,k i,k + i,k d i τdτ. 33 Now consider he erm r r i,k in 32. Using 4, we can wrie r r i,k = i,k rτdτ. 34 Finally, consider he erms x x i,k in 32. For hese erms we need o disinguish wo cases, namely,h +1 and >,h +1. Noice ha he laer case corresponds o he fac ha agen updaes is conrol inpu o a value unknown o agen i. In he firs case, inegraing 2 for agen in [,h,, using 3 and 7, and noing ha u = u,h for [,h,,h +1, we have x x i,k =u,h i,k +,h d τdτ, [,h,,h In he second case, similar observaions lead o x x i,k =u,h,h +1 i,k +,h +1 u τdτ + Subsiuing ino 32 yields u i =c i + p i u i,k i,k u,h min{,,h +1 } i,k i + i + p i d i τ p i rdτ i,k d τdτ,h i i,h +1,h d τdτ, >,h +1.,h u τdτ. 37 Taking norms of boh sides in 37, and using he riangular inequaliy and Assumpion III.1, we have u i c i + p i u i,k i,k u,h min{,,h +1 } i,k i + i + p i Δ i,k + p i δ 0 τdτ i,k + Δ,h + u τ dτ i,h +1 i,h +1, 38 for [ i,k, i,k+1. Nex, suppose by conradicion ha some agen i a some ime [ i,k, i,k+1 aains u i > ς i, while u q ς q for all [0, and all q. Then, using Corollary IV.1, we have u τ β ητ + ς τ [0, i. 39 Evaluaing 38 for =, and using 39, we have u i σ i,k, 40 where σ i,k is defined in 11. By 40, u i > ς i implies σ i,k > ς i. Bu his is a conradicion, since he scheduling rule 8 12 and 20 guaranees ha σ i,k ς i for all [ i,k, i,k+1 and all k N 0. To exclude ha he sysem exhibis Zeno behavior, consider he condiions 8 ha rigger he cloud accesses. From 9, we see ha he riggering condiion Δ i,k ζ i,k requires i,k ς q, 2cδ i,0 q for some q, and 2211

6 herefore, i canno generae Zeno behavior. Nex, consider he condiion σ i,k ς i. Evaluaing 11 for = i,k we have σ i,k i,k = c i Δ,h i,k. 41 Recalling ha i,k [,h,,h +1, and noing ha 8 guaranees Δ,h < ς i 2c i for all [,h,,h +1, from 41 we have σ i,k i,k ς i i,k Differeniaing boh sides of 11, and using 9, he coninuiy of σ i,k and he riangular inequaliy, we have σ i,k σ i,k i,k + i,k s i,k τdτ, 43 where s i,k =c i + p i u i,k + δ i + p i δ 0 + <,h +1 u,h + i δ +,h +1 μ 44 From Corollary IV.1, we have u i,k = u i μ i and u,h = u μ for all i such ha <,h +1. Since μ = β η+ς is upper-bounded, δ δ,0 and ρ ρ 0 for all 0, from 44, we have s i,k c i + p i μ i + δ i,0 + p i ρ 0 + μ + δ,0, 45 where μ denoes he maximum value aained by he funcion μ. Denoing he righ-hand side of 45 as s i,k, and subsiuing 42 and 45 ino 43, we have σ i,k ς i i,k 2 + s i,k i,k 46 From 46, a necessary condiion for having σ i,k ς i is ς i i,k 2 + s i,k i,k ς i. 47 Observing ha ς i = ς i i,k ς i, e λ ς i,k +ς i,, we can rewrie 47 as ς i i,k 2 + s i,k i,k ς i i,k ς i, e λ ς i,k +ς i,. For any ς i i,k > ς i, 0 and any λ ς > 0, he posiive soluions in he unknown τ of he equaion ς i i,k 2 + s i,k τ ς i i,k ς i, e λςτ +ς i, is lower-bounded. Therefore, condiion 47 canno be saisfied for i,k arbirarily small. Consequenly, he riggering condiion ς i,k ς i canno generae Zeno behavior eiher. We can conclude ha he closed-loop sysem defined by 2, 5 12 and 20 does no exhibi Zeno behavior. Remark IV.2. Agen i can compue μ for by 20 and 28, and herefore by only using some neighborhood informaion on he nework opology β for and he iniial condiions y 0. Theorems IV.1 and IV.2 amoun o our main resul, which is formalized as follows. Theorem IV.3. Consider he muli-agen sysem 2, le Assumpions III.1 and IV.1 hold, and le he sysem be conrolled by he algorihm defined by 5 12 and 20. Then he closed-loop sysem does no exhibi Zeno behavior and achieves pracical consensus as by Definiion III.1, wih olerance ε = max i { m i }η, where m i is he number of edges in he shores pah from verex 0 o verex i in he graph,, and η = lim η = α B δ +ς, where δ cλ = [δ 0,, δ 1,,, δ N, ] and ς = [0, ς 1,,, ς N, ]. Proof. From Theorems IV.1 and IV.2, we have y η for all 0, where η is defined by 20. Leing, we have herefore lim sup y η. Finally, observing ha r x i m i y yields he desired resul. V. ASYMPTOTIC CONVERGENCE If he disurbances vanish quickly enough, and he reference raecory converges quickly enough o a fixed poin, hen he proposed algorihm, wih only small adusmens, is capable o drive all he agens o he reference poin asympoically. In his case, he following assumpion is needed. Assumpion V.1. Assumpion III.1 holds wih ρ, δ 1,,, δ N, = 0 and λ δ < c min{eigl r }. wih In his scenario, he hreshold funcions are chosen as ς i = ς i,0 e λ ς, 48 0 < λ ς < λ δ < c min{eigl r }. 49 Noe ha Theorem IV.1 and Corollary IV.1 sill hold. Moreover, solving he inegral in 20, and using Assumpion V.1 and 49, we have η η e λ ς, 50 where η = α y 0 + δ0 + ς0, 51 cλ λ δ cλ λ ς and δ, ς are defined in he saemen of Theorem IV.1. In he following heorem we show ha his version of he proposed algorihm is sill Zeno-free. Theorem V.1. Consider he muli-agen sysem 2, le Assumpions IV.1 and V.1 hold, and le he sysem be conrolled by he algorihm defined by 5 11, 20 and 48, wih λ ς saisfying 49. Then, he closed-loop sysem does no exhibi Zeno behavior and u i ς i holds for all 0 and all i. Proof. Reasoning as in Theorem IV.2, we can show ha he scheduling rule 8 11, 20 and 48 guaranees u i ς i for any 0. To show ha he closed-loop sysem is Zeno free, consider again he condiion 8 ha riggers he cloud accesses. From 9, we can see ha he riggering condiion Δ i,k ζ i,k requires δ q,0 e λ δ i,k 1 e λ δ i,k λ δ 2212

7 ς q,0 2c q e λ ς i,k e λ ς i,k for some q. By 49, he previous inequaliy implies ς q,0 2c q + δ q,0 e λ ς i,k δ q,0 λ δ Therefore, he condiion Δ i,k ζ i,k canno generae Zeno behavior. Nex, consider he condiion σ i,k ς i. Wih similar reasoning as in Theorem IV.2, we can show ha 42 and 43 sill hold. Using Corollary IV.1 and 48 and 50, and recalling ha [ i,k, i,k+1, we have u i,k = u i μ i β i η + ς i,0 e λ ς and u,h = λ δ. u μ β η + ς,0 e λ ς for all i such ha <,h +1. Replacing hese wo inequaliies ino 44, and using 48 50, we have s i,k c i + p i β i η + ς i,0 + δ i,0 + p i ρ β η + ς,0 e λ ς [ i,k, i,k+1, where η is defined in 51. Subsiuing 42 and 52 in 43, we have σ i,k ς i i,k 2 + ξ i,k λ ς e λ ς i,k 1 e λ ς i,k, 53 where ξ i,k denoes he coefficien ha muliplies e λ ς in 52. From 53, a necessary condiion for having σ i,k ς i is ς i i,k 2 + ξ i,k λ ς e λ ς i,k 1 e λ ς i,k ς i. 54 Observing ha ς i = ς i i,k e λ ς i,k and ς i i,k = ς i,0 e λ ς i,k, 54 can be rewrien as ς i,0 2 + ξ i,k ς λ i,0 + ξ i,k e λ ς i,k, ς λ ς which has lower-bounded soluions in he unknown i,k. Hence, he riggering condiion σ i,k ς i canno generae Zeno behavior eiher. We can conclude ha he closed-loop sysem does no exhibi Zeno behavior. Theorems IV.1 and V.1 amoun o our asympoic convergence resul, which is formalized as follows. Theorem V.2. Consider he muli-agen sysem 2, le Assumpions IV.1 and V.1 hold, and le he sysem be conrolled by he algorihm defined by 5 11, 20 and 48, wih λ ς saisfying 49. Then he closed-loop sysem does no exhibi Zeno behavior and lim r x i = 0 for all i. Proof. Reasoning as in Theorem IV.3, we have r x i m i η, where m i is defined in Theorem IV.3. From 50, leing in he previous inequaliy yields he desired convergence resul. VI. NUMERICAL SIMULATIONS In his secion, we presen wo numerical simulaions of he proposed algorihm, which demonsrae respecively pracical consensus and asympoic convergence. We consider a sysem made up of N = 4 agens plus he reference raecory, wih graph opology as in Figure 1. Noe ha for Fig. 1. Topology of he muli-agen sysem used in he simulaions. The node r represens he reference raecory r 1 x 1 1 x1 2 x1 3 x Fig. 2. Resuls of he firs simulaion scenario. Top: firs posiion variable x 1 i for each agen and r 1 for he reference raecory. Boom: error norm x i r for each agen. his opology Assumpion IV.1 is saisfied, and we have λ = min{eigl r } = 0.53, B = 1.90 and α = For he firs simulaion, we le c = 1.8, ς = [0.0, , 0.48, 0.57, 0.67] 5.0 e , δ = [99.9, 1.0, 2.0, 3.0, 4.0] 10 3 for all 0. We le he derivaive of he reference raecory and he disurbances be sinusoidal, namely r = δ 0 [cos2πf 0, sin2πf 0 ], d i = δ i [cos2πf i, sin2πf i ]. Noe ha wih his choice Assumpion III.1 is saisfied. The frequencies are chosen as f i = i T for i {0, 1,, N}, where T is he simulaion ime. The simulaion runs for [0.0, 2.5], wih a fixed sep of 10 4 he physical ime scale can be chosen according o he paricular applicaion. Figures 2 and 4 show he resuls of his simulaion. We have a oal of 305 updaes corresponding o an average iner-updae ime of 0.033, wo order of magniudes larger han he simulaion sep. Hence, he simulaion corroboraes he absence of Zeno behavior. For he second simulaion, we le c = 1.0, ς = [0.0, , 0.48, 0.57, 0.67] 10.0 e 1.0, δ = [99.9, 1.0, 2.0, 3.0, 4.0] 10 3 e 1.1. Noe ha, wih hese choices, Assumpion V.1 and 48 are saisfied, and he derivaive of he reference raecory and he disurbances vanish asympoically. The simulaion runs for [0.0, 10.0], wih he same sep as before. Figures 3 and 4 show he resuls of his simulaion. In his case, we have a oal of 1418 updaes, corresponding o an average inereven ime of , sill wo order of magniudes larger han he simulaion sep. Hence, he simulaion corroboraes 2213

8 i i Fig. 4. Top: updaes in he firs simulaion scenario for each agen. Boom: updaes in he second simulaion scenario for [9.0, 10.0] for each agen r 1 x 1 1 x1 2 x1 3 x Fig. 3. Resuls of he second simulaion scenario. Top: firs posiion variable x 1 i for each agen and r 1 for he reference raecory. Boom: error norm x i r for each agen. he absence of Zeno behavior. VII. CONCLUSIONS AND FUTURE DEVELOPMENTS A cloud-suppored conrol algorihm for leader-follower raecory racking in a nework of mobile agens under disurbances has been proposed. The considered seup allows muli-agen coordinaion in case of inerdiced communicaion among he agens. Specifically, he scenario of conrolling a formaion of AUVs has been considered as a moivaing example. The conrol algorihm overcomes he limiaion of having a pre-assigned raecory for he whole flee as well as synchronizing he surfacing of he agens [7]. Sufficien condiions for boh bounded and asympoic convergence have been idenified, in erms of he opology of he informaion exchanges wih he cloud, and of he scheduling of he conrol updaes. The proposed conrol algorihm is effecive in guaraneeing he overall sabiliy despie each agen receiving oudaed informaion and no knowing oher agens fuure conrol inpus. Fuure work will furher develop he approach of he paper considering more complex agen dynamics and conrol obecives. Also, non-idealiies in he cloud access, such as delays and packe losses, will be aken ino accoun. REFERENCES [1] W. Ren, Disribued Consensus in Muli-vehicle Cooperaive Conrol, IEEE Conrol Sysems Magazine, April [2] R. Olfai-Saber, Flocking for muli-agen dynamic sysems: Algorihms and heory, IEEE Transacions on Auomaic Conrol, vol. 51, no. 3, [3] J. A. Fax and R. M. Murray, Graph laplacians and sabilizaion of vehicle formaions, IEEE Transacions on Auomaic Conrol, vol. 15, no. 1, [4] L. Suers, H. Liu, C. Tilman, and D. J. Brown, Navigaion echnologies for auonomous underwaer vehicles, IEEE Transacions on Sysems, Man and Cyberneics Par C: Applicaions and Reviews, vol. 38, no. 4, [5] L. Paull, S. Saeedi, M. Seo, and H. Li, AUV navigaion and localizaion: A review, IEEE Journal of Oceanic Engineering, vol. 39, no. 1, [6] J. W. Nicholson and A. J. Healey, The Presen Sae of Auonomous Underwaer Vehicle AUV Applicaions and Technologies, Marine Technology Sociey Journal, vol. 42, no. 1, [7] E. Fiorelli, N. E. Leonard, P. Bhaa, D. A. Paley, R. Bachmayer, and D. M. Fraanoni, Muli-AUV conrol and adapive sampling in Monerey Bay, IEEE Journal of Oceanic Engineering, vol. 31, no. 4, pp , [8] W. P. M. H. Heemels, K. H. Johansson, and P. Tabuada, An Inroducion o Even-riggered and Self-riggered Conrol, in IEEE Conference on Decision and Conrol, [9] D. V. Dimarogonas, E. Frazzoli, and K. H. Johansson, Disribued even-riggered conrol for muli-agen sysems, IEEE Transacions on Auomaic Conrol, vol. 57, no. 5, pp , [10] M. T. Hale and M. Egersed, Differenially privae cloud-based muliagen opimizaion wih consrains, in American Conrol Conference, [11] P. Pandey, D. Pompili, and J. Yi, Dynamic Collaboraion beween Neworked Robos and Clouds in Resource-Consrained Environmens, IEEE Transacions on Auomaion Science and Engineering, vol. 12, no. 2, pp , [12] A. Adaldo, D. Liuzza, D. V. Dimarogonas, and K. H. Johansson, Conrol of Muli-Agen Sysems wih Even-Triggered Cloud Access, in European Conrol Conference, [13] D. Zelazo, A. Rahmani, and M. Mesbahi, Agreemen via he edge Laplacian, in IEEE Conference on Decision and Conrol, [14] Z. Zeng, X. Wang, and Z. Zheng, Nonlinear Consensus under Direced Graph via he Edge Laplacian, in Chinese Conrol and Decision Conference, [15] M. Mesbahi and M. Egersed, Graph Theoreic Mehods in Muliagen Neworks. Princeon Univerisiy Press, [16] R. A. Horn and C. R. Johnson, Marix Analysis, 2nd ed. Cambridge Universiy Press,

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

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

15. Vector Valued Functions

15. Vector Valued Functions 1. Vecor Valued Funcions Up o his poin, we have presened vecors wih consan componens, for example, 1, and,,4. However, we can allow he componens of a vecor o be funcions of a common variable. For example,

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

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

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

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

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

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

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

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

Some Ramsey results for the n-cube

Some Ramsey results for the n-cube Some Ramsey resuls for he n-cube Ron Graham Universiy of California, San Diego Jozsef Solymosi Universiy of Briish Columbia, Vancouver, Canada Absrac In his noe we esablish a Ramsey-ype resul for cerain

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

POSITIVE SOLUTIONS OF NEUTRAL DELAY DIFFERENTIAL EQUATION

POSITIVE SOLUTIONS OF NEUTRAL DELAY DIFFERENTIAL EQUATION Novi Sad J. Mah. Vol. 32, No. 2, 2002, 95-108 95 POSITIVE SOLUTIONS OF NEUTRAL DELAY DIFFERENTIAL EQUATION Hajnalka Péics 1, János Karsai 2 Absrac. We consider he scalar nonauonomous neural delay differenial

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

Global Synchronization of Directed Networks with Fast Switching Topologies

Global Synchronization of Directed Networks with Fast Switching Topologies Commun. Theor. Phys. (Beijing, China) 52 (2009) pp. 1019 1924 c Chinese Physical Sociey and IOP Publishing Ld Vol. 52, No. 6, December 15, 2009 Global Synchronizaion of Direced Neworks wih Fas Swiching

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

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

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still. Lecure - Kinemaics in One Dimension Displacemen, Velociy and Acceleraion Everyhing in he world is moving. Nohing says sill. Moion occurs a all scales of he universe, saring from he moion of elecrons in

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

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

CERTAIN CLASSES OF SOLUTIONS OF LAGERSTROM EQUATIONS

CERTAIN CLASSES OF SOLUTIONS OF LAGERSTROM EQUATIONS SARAJEVO JOURNAL OF MATHEMATICS Vol.10 (22 (2014, 67 76 DOI: 10.5644/SJM.10.1.09 CERTAIN CLASSES OF SOLUTIONS OF LAGERSTROM EQUATIONS ALMA OMERSPAHIĆ AND VAHIDIN HADŽIABDIĆ Absrac. This paper presens sufficien

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

Notes for Lecture 17-18

Notes for Lecture 17-18 U.C. Berkeley CS278: Compuaional Complexiy Handou N7-8 Professor Luca Trevisan April 3-8, 2008 Noes for Lecure 7-8 In hese wo lecures we prove he firs half of he PCP Theorem, he Amplificaion Lemma, up

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

11!Hí MATHEMATICS : ERDŐS AND ULAM PROC. N. A. S. of decomposiion, properly speaking) conradics he possibiliy of defining a counably addiive real-valu

11!Hí MATHEMATICS : ERDŐS AND ULAM PROC. N. A. S. of decomposiion, properly speaking) conradics he possibiliy of defining a counably addiive real-valu ON EQUATIONS WITH SETS AS UNKNOWNS BY PAUL ERDŐS AND S. ULAM DEPARTMENT OF MATHEMATICS, UNIVERSITY OF COLORADO, BOULDER Communicaed May 27, 1968 We shall presen here a number of resuls in se heory concerning

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

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients Secion 3.5 Nonhomogeneous Equaions; Mehod of Undeermined Coefficiens Key Terms/Ideas: Linear Differenial operaor Nonlinear operaor Second order homogeneous DE Second order nonhomogeneous DE Soluion o homogeneous

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

An introduction to the theory of SDDP algorithm

An introduction to the theory of SDDP algorithm An inroducion o he heory of SDDP algorihm V. Leclère (ENPC) Augus 1, 2014 V. Leclère Inroducion o SDDP Augus 1, 2014 1 / 21 Inroducion Large scale sochasic problem are hard o solve. Two ways of aacking

More information

Echocardiography Project and Finite Fourier Series

Echocardiography Project and Finite Fourier Series Echocardiography Projec and Finie Fourier Series 1 U M An echocardiagram is a plo of how a porion of he hear moves as he funcion of ime over he one or more hearbea cycles If he hearbea repeas iself every

More information

Comments on Window-Constrained Scheduling

Comments on Window-Constrained Scheduling Commens on Window-Consrained Scheduling Richard Wes Member, IEEE and Yuing Zhang Absrac This shor repor clarifies he behavior of DWCS wih respec o Theorem 3 in our previously published paper [1], and describes

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

Random Walk with Anti-Correlated Steps

Random Walk with Anti-Correlated Steps Random Walk wih Ani-Correlaed Seps John Noga Dirk Wagner 2 Absrac We conjecure he expeced value of random walks wih ani-correlaed seps o be exacly. We suppor his conjecure wih 2 plausibiliy argumens and

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

( ) ( ) 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

Challenge Problems. DIS 203 and 210. March 6, (e 2) k. k(k + 2). k=1. f(x) = k(k + 2) = 1 x k

Challenge Problems. DIS 203 and 210. March 6, (e 2) k. k(k + 2). k=1. f(x) = k(k + 2) = 1 x k Challenge Problems DIS 03 and 0 March 6, 05 Choose one of he following problems, and work on i in your group. Your goal is o convince me ha your answer is correc. Even if your answer isn compleely correc,

More information

arxiv: v1 [math.gm] 4 Nov 2018

arxiv: v1 [math.gm] 4 Nov 2018 Unpredicable Soluions of Linear Differenial Equaions Mara Akhme 1,, Mehme Onur Fen 2, Madina Tleubergenova 3,4, Akylbek Zhamanshin 3,4 1 Deparmen of Mahemaics, Middle Eas Technical Universiy, 06800, Ankara,

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

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

KEY. Math 334 Midterm III Winter 2008 section 002 Instructor: Scott Glasgow

KEY. Math 334 Midterm III Winter 2008 section 002 Instructor: Scott Glasgow KEY Mah 334 Miderm III Winer 008 secion 00 Insrucor: Sco Glasgow Please do NOT wrie on his exam. No credi will be given for such work. Raher wrie in a blue book, or on your own paper, preferably engineering

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

4.6 One Dimensional Kinematics and Integration

4.6 One Dimensional Kinematics and Integration 4.6 One Dimensional Kinemaics and Inegraion When he acceleraion a( of an objec is a non-consan funcion of ime, we would like o deermine he ime dependence of he posiion funcion x( and he x -componen of

More information

We just finished the Erdős-Stone Theorem, and ex(n, F ) (1 1/(χ(F ) 1)) ( n

We just finished the Erdős-Stone Theorem, and ex(n, F ) (1 1/(χ(F ) 1)) ( n Lecure 3 - Kövari-Sós-Turán Theorem Jacques Versraëe jacques@ucsd.edu We jus finished he Erdős-Sone Theorem, and ex(n, F ) ( /(χ(f ) )) ( n 2). So we have asympoics when χ(f ) 3 bu no when χ(f ) = 2 i.e.

More information

Ordinary Differential Equations

Ordinary Differential Equations Ordinary Differenial Equaions 5. Examples of linear differenial equaions and heir applicaions We consider some examples of sysems of linear differenial equaions wih consan coefficiens y = a y +... + a

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

MATH 128A, SUMMER 2009, FINAL EXAM SOLUTION

MATH 128A, SUMMER 2009, FINAL EXAM SOLUTION MATH 28A, SUMME 2009, FINAL EXAM SOLUTION BENJAMIN JOHNSON () (8 poins) [Lagrange Inerpolaion] (a) (4 poins) Le f be a funcion defined a some real numbers x 0,..., x n. Give a defining equaion for he Lagrange

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

A Necessary and Sufficient Condition for the Solutions of a Functional Differential Equation to Be Oscillatory or Tend to Zero

A Necessary and Sufficient Condition for the Solutions of a Functional Differential Equation to Be Oscillatory or Tend to Zero JOURNAL OF MAEMAICAL ANALYSIS AND APPLICAIONS 24, 7887 1997 ARICLE NO. AY965143 A Necessary and Sufficien Condiion for he Soluions of a Funcional Differenial Equaion o Be Oscillaory or end o Zero Piambar

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

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

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

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

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

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

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

Solutions from Chapter 9.1 and 9.2

Solutions from Chapter 9.1 and 9.2 Soluions from Chaper 9 and 92 Secion 9 Problem # This basically boils down o an exercise in he chain rule from calculus We are looking for soluions of he form: u( x) = f( k x c) where k x R 3 and k is

More information

On Oscillation of a Generalized Logistic Equation with Several Delays

On Oscillation of a Generalized Logistic Equation with Several Delays Journal of Mahemaical Analysis and Applicaions 253, 389 45 (21) doi:1.16/jmaa.2.714, available online a hp://www.idealibrary.com on On Oscillaion of a Generalized Logisic Equaion wih Several Delays Leonid

More information

Chapter 8 The Complete Response of RL and RC Circuits

Chapter 8 The Complete Response of RL and RC Circuits Chaper 8 The Complee Response of RL and RC Circuis Seoul Naional Universiy Deparmen of Elecrical and Compuer Engineering Wha is Firs Order Circuis? Circuis ha conain only one inducor or only one capacior

More information

Solutions to Assignment 1

Solutions to Assignment 1 MA 2326 Differenial Equaions Insrucor: Peronela Radu Friday, February 8, 203 Soluions o Assignmen. Find he general soluions of he following ODEs: (a) 2 x = an x Soluion: I is a separable equaion as we

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

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

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

Lecture Notes 2. The Hilbert Space Approach to Time Series

Lecture Notes 2. The Hilbert Space Approach to Time Series Time Series Seven N. Durlauf Universiy of Wisconsin. Basic ideas Lecure Noes. The Hilber Space Approach o Time Series The Hilber space framework provides a very powerful language for discussing he relaionship

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

RANDOM LAGRANGE MULTIPLIERS AND TRANSVERSALITY

RANDOM LAGRANGE MULTIPLIERS AND TRANSVERSALITY ECO 504 Spring 2006 Chris Sims RANDOM LAGRANGE MULTIPLIERS AND TRANSVERSALITY 1. INTRODUCTION Lagrange muliplier mehods are sandard fare in elemenary calculus courses, and hey play a cenral role in economic

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

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

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

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

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

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

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

SUFFICIENT CONDITIONS FOR EXISTENCE SOLUTION OF LINEAR TWO-POINT BOUNDARY PROBLEM IN MINIMIZATION OF QUADRATIC FUNCTIONAL

SUFFICIENT CONDITIONS FOR EXISTENCE SOLUTION OF LINEAR TWO-POINT BOUNDARY PROBLEM IN MINIMIZATION OF QUADRATIC FUNCTIONAL HE PUBLISHING HOUSE PROCEEDINGS OF HE ROMANIAN ACADEMY, Series A, OF HE ROMANIAN ACADEMY Volume, Number 4/200, pp 287 293 SUFFICIEN CONDIIONS FOR EXISENCE SOLUION OF LINEAR WO-POIN BOUNDARY PROBLEM IN

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

Properties Of Solutions To A Generalized Liénard Equation With Forcing Term

Properties Of Solutions To A Generalized Liénard Equation With Forcing Term Applied Mahemaics E-Noes, 8(28), 4-44 c ISSN 67-25 Available free a mirror sies of hp://www.mah.nhu.edu.w/ amen/ Properies Of Soluions To A Generalized Liénard Equaion Wih Forcing Term Allan Kroopnick

More information

Differential Equations

Differential Equations Mah 21 (Fall 29) Differenial Equaions Soluion #3 1. Find he paricular soluion of he following differenial equaion by variaion of parameer (a) y + y = csc (b) 2 y + y y = ln, > Soluion: (a) The corresponding

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

Electrical and current self-induction

Electrical and current self-induction Elecrical and curren self-inducion F. F. Mende hp://fmnauka.narod.ru/works.hml mende_fedor@mail.ru Absrac The aricle considers he self-inducance of reacive elemens. Elecrical self-inducion To he laws of

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

Signal and System (Chapter 3. Continuous-Time Systems)

Signal and System (Chapter 3. Continuous-Time Systems) Signal and Sysem (Chaper 3. Coninuous-Time Sysems) Prof. Kwang-Chun Ho kwangho@hansung.ac.kr Tel: 0-760-453 Fax:0-760-4435 1 Dep. Elecronics and Informaion Eng. 1 Nodes, Branches, Loops A nework wih b

More information

Approximation Algorithms for Unique Games via Orthogonal Separators

Approximation Algorithms for Unique Games via Orthogonal Separators Approximaion Algorihms for Unique Games via Orhogonal Separaors Lecure noes by Konsanin Makarychev. Lecure noes are based on he papers [CMM06a, CMM06b, LM4]. Unique Games In hese lecure noes, we define

More information

Research Article Existence and Uniqueness of Periodic Solution for Nonlinear Second-Order Ordinary Differential Equations

Research Article Existence and Uniqueness of Periodic Solution for Nonlinear Second-Order Ordinary Differential Equations Hindawi Publishing Corporaion Boundary Value Problems Volume 11, Aricle ID 19156, 11 pages doi:1.1155/11/19156 Research Aricle Exisence and Uniqueness of Periodic Soluion for Nonlinear Second-Order Ordinary

More information

Optimal Path Planning for Flexible Redundant Robot Manipulators

Optimal Path Planning for Flexible Redundant Robot Manipulators 25 WSEAS In. Conf. on DYNAMICAL SYSEMS and CONROL, Venice, Ialy, November 2-4, 25 (pp363-368) Opimal Pah Planning for Flexible Redundan Robo Manipulaors H. HOMAEI, M. KESHMIRI Deparmen of Mechanical Engineering

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

8. Basic RL and RC Circuits

8. Basic RL and RC Circuits 8. Basic L and C Circuis This chaper deals wih he soluions of he responses of L and C circuis The analysis of C and L circuis leads o a linear differenial equaion This chaper covers he following opics

More information

= ( ) ) or a system of differential equations with continuous parametrization (T = R

= ( ) ) or a system of differential equations with continuous parametrization (T = R XIII. DIFFERENCE AND DIFFERENTIAL EQUATIONS Ofen funcions, or a sysem of funcion, are paramerized in erms of some variable, usually denoed as and inerpreed as ime. The variable is wrien as a funcion of

More information

R t. C t P t. + u t. C t = αp t + βr t + v t. + β + w t

R t. C t P t. + u t. C t = αp t + βr t + v t. + β + w t Exercise 7 C P = α + β R P + u C = αp + βr + v (a) (b) C R = α P R + β + w (c) Assumpions abou he disurbances u, v, w : Classical assumions on he disurbance of one of he equaions, eg. on (b): E(v v s P,

More information

System of Linear Differential Equations

System of Linear Differential Equations Sysem of Linear Differenial Equaions In "Ordinary Differenial Equaions" we've learned how o solve a differenial equaion for a variable, such as: y'k5$e K2$x =0 solve DE yx = K 5 2 ek2 x C_C1 2$y''C7$y

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

The L p -Version of the Generalized Bohl Perron Principle for Vector Equations with Infinite Delay

The L p -Version of the Generalized Bohl Perron Principle for Vector Equations with Infinite Delay Advances in Dynamical Sysems and Applicaions ISSN 973-5321, Volume 6, Number 2, pp. 177 184 (211) hp://campus.ms.edu/adsa The L p -Version of he Generalized Bohl Perron Principle for Vecor Equaions wih

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

Longest Common Prefixes

Longest Common Prefixes Longes Common Prefixes The sandard ordering for srings is he lexicographical order. I is induced by an order over he alphabe. We will use he same symbols (,

More information

EECE251. Circuit Analysis I. Set 4: Capacitors, Inductors, and First-Order Linear Circuits

EECE251. Circuit Analysis I. Set 4: Capacitors, Inductors, and First-Order Linear Circuits EEE25 ircui Analysis I Se 4: apaciors, Inducors, and Firs-Order inear ircuis Shahriar Mirabbasi Deparmen of Elecrical and ompuer Engineering Universiy of Briish olumbia shahriar@ece.ubc.ca Overview Passive

More information

Ann. Funct. Anal. 2 (2011), no. 2, A nnals of F unctional A nalysis ISSN: (electronic) URL:

Ann. Funct. Anal. 2 (2011), no. 2, A nnals of F unctional A nalysis ISSN: (electronic) URL: Ann. Func. Anal. 2 2011, no. 2, 34 41 A nnals of F uncional A nalysis ISSN: 2008-8752 elecronic URL: www.emis.de/journals/afa/ CLASSIFICAION OF POSIIVE SOLUIONS OF NONLINEAR SYSEMS OF VOLERRA INEGRAL EQUAIONS

More information

A New Perturbative Approach in Nonlinear Singularity Analysis

A New Perturbative Approach in Nonlinear Singularity Analysis Journal of Mahemaics and Saisics 7 (: 49-54, ISSN 549-644 Science Publicaions A New Perurbaive Approach in Nonlinear Singulariy Analysis Ta-Leung Yee Deparmen of Mahemaics and Informaion Technology, The

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

MATH 4330/5330, Fourier Analysis Section 6, Proof of Fourier s Theorem for Pointwise Convergence

MATH 4330/5330, Fourier Analysis Section 6, Proof of Fourier s Theorem for Pointwise Convergence MATH 433/533, Fourier Analysis Secion 6, Proof of Fourier s Theorem for Poinwise Convergence Firs, some commens abou inegraing periodic funcions. If g is a periodic funcion, g(x + ) g(x) for all real x,

More information

THE 2-BODY PROBLEM. FIGURE 1. A pair of ellipses sharing a common focus. (c,b) c+a ROBERT J. VANDERBEI

THE 2-BODY PROBLEM. FIGURE 1. A pair of ellipses sharing a common focus. (c,b) c+a ROBERT J. VANDERBEI THE 2-BODY PROBLEM ROBERT J. VANDERBEI ABSTRACT. In his shor noe, we show ha a pair of ellipses wih a common focus is a soluion o he 2-body problem. INTRODUCTION. Solving he 2-body problem from scrach

More information

INDEPENDENT SETS IN GRAPHS WITH GIVEN MINIMUM DEGREE

INDEPENDENT SETS IN GRAPHS WITH GIVEN MINIMUM DEGREE INDEPENDENT SETS IN GRAPHS WITH GIVEN MINIMUM DEGREE JAMES ALEXANDER, JONATHAN CUTLER, AND TIM MINK Absrac The enumeraion of independen ses in graphs wih various resricions has been a opic of much ineres

More information