This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and

Size: px
Start display at page:

Download "This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and"

Transcription

1 Ths artcle appeared n a ournal publshed by Elsever The attached copy s furnshed to the author for nternal non-commercal research and educaton use, ncludng for nstructon at the authors nsttuton and sharng wth colleagues Other uses, ncludng reproducton and dstrbuton, or sellng or lcensng copes, or postng to personal, nsttutonal or thrd party webstes are prohbted n most cases authors are permtted to post ther verson of the artcle (eg n Word or Tex form to ther personal webste or nsttutonal repostory Authors requrng further nformaton regardng Elsever s archvng and manuscrpt polces are encouraged to vst:

2 Systems & Control Letters 59 ( Contents lsts avalable at ScenceDrect Systems & Control Letters ournal homepage: wwwelsevercom/locate/sysconle Partal state consensus for networs of second-order dynamc agents Feng Xao a,,, Long Wang b, Je Chen a a School of Automaton, Beng nsttute of Technology, Beng 0008, Chna b ntellgent Control Laboratory, Center for Systems and Control, College of Engneerng, and Key Laboratory of Machne Percepton (Mnstry of Educaton, Peng Unversty, Beng 0087, Chna a r t c l e n f o a b s t r a c t Artcle hstory: Receved 4 March 00 Receved n revsed form 3 August 00 Accepted 4 September 00 Avalable onlne 3 October 00 Keywords: Mult-agent systems Partal state consensus Tme-delay Swtchng topology Ths paper addresses the partal state consensus problem of mult-agent systems wth second-order agent dynamcs and proposes an asynchronous dstrbuted consensus protocol for the case wth swtchng nteracton topology, tme-varyng delays and ntermttent nformaton transmsson Partal state consensus means reachng an agreement asymptotcally wth each other on part, but not all, of each ndvdual s states, where the concerned states usually cannot be decoupled from the other ones Partal state consensus has ts broad applcatons n the coordnaton of mult-robot systems, dstrbuted tas management, and dstrbuted estmaton for sensor networs, etc Ths paper assumes that poston-le states are the only detectable nformaton transmtted over the networ and velocty-le states are the ey quanttes of nterest, whch are requred to be equalzed We frst gve the asynchronous dstrbuted protocol based on the delayed poston-le state nformaton and then provde ts convergence result wth respect to velocty-le states t s shown that f the unon of the nteracton topology across the tme nterval wth a gven length always contans a spannng tree, then the proposed protocol wll solve the partal state (velocty-le state consensus problem asymptotcally 00 Elsever BV All rghts reserved ntroducton Consensus theory, as a basc and fundamental research topc n dstrbuted coordnaton of mult-agent systems, has receved consderable attenton from researchers n recent years because of ts potental applcatons n cooperatve control of unmanned ar vehcles, formaton control of moble robots, desgn of sensor networs, swarm-based computng, etc [ 5] t requres that all agents reach an agreement on certan quanttes of nterest The shared common value may be the antcpated atttude n multple spacecraft algnment, poston and velocty n flocng control, or processng rate n dstrbuted tas management Due to the complexty of networ dynamcs, a large amount of wors assume that each ndvdual follows a frst-order dfferental equaton [6,7] And there s also a fracton of the lterature concentratng on the networs wth second-order agent dynamcs For nstance, n [8], Xe and Wang consdered the state consensus problem n networs of multple double-ntegrators Ths wor was supported by NSFC ( , and and the Beng Educaton Commttee Cooperaton Buldng Foundaton Proect Correspondng author E-mal addresses: fxao@eceualbertaca, fengxao@pueducn (F Xao He s now vstng the Department of Electrcal and Computer Engneerng, Unversty of Alberta, Edmonton, Alberta T6G V4, Canada under fxed and swtchng topologes and desgned the protocols ensurng that poston-le states converge to a common statc value asymptotcally n [9], Hong et al used a set of frstorder agents to trac an actve second-order leader, where the consensus state may be tme-varyng Ren proposed and analyzed consensus algorthms for networs of double-ntegrators n [0,] and second-order lnear harmonc oscllators n [] under the assumpton that each agent can fully or partally access ts neghbors relatve states Along ths lne, fnte-tme agreement for multple double-ntegrators was dscussed n [3], where the proposed protocols are non-smooth and ther consensus property was proved by the theory of fnte-tme homogeneous systems Under these algorthms, all agents states, ncludng postonle and velocty-le states, wll reach consensus asymptotcally However, n many practcal stuatons, t often occurs that only a small part of state varables of each agent are the ey quanttes of nterest that are requred to be coordnated and usually cannot be decoupled from the other ones n these cases, the concept of consensus n the tradtonal sense wll not be sutable Partal state consensus means reachng an agreement asymptotcally wth each other on part, but not all, of each ndvdual s states, and n ts studes, we may face many new challenges caused by the constrant of coupled agreement and nonagreement varables, such as characterzng the partal state consensus problem n the framewor of coordnaton of full state varables, dstrbuted estmatons for nterestng states of neghborng agents va the /$ see front matter 00 Elsever BV All rghts reserved do:006/sysconle

3 776 F Xao et al / Systems & Control Letters 59 ( nonagreement coupled nformaton, and so on Those challenges greatly ncrease the dffculty of desgn and analyss of partal state consensus protocols As the frst step towards the general study on partal state consensus, ths paper consders the smplest case and nvestgates the networs of second-order dynamc agents t s assumed that poston-le states are the only detectable nformaton transmtted over the networ and velocty-le states are the quanttes of nterest, whch are requred to be equalzed Next, we show several scenaros where the velocty-le state s the only quantty of nterest The frst example ncludes the varous versons of the Vcse model [7,4,5], where the velocty consensus s a prerequste for further study Also n the theoretcal study, n [6], Barbarossa and Scutar proposed a decentralzed sensor networ scheme capable of reachng a globally optmum maxmum-lelhood estmate through self-synchronzaton of nonlnearly coupled dynamcal systems Although each node n the networ they studed s a frst-order dynamcal system, the fnal agreement estmate s related to the state dervatve of each sensor, namely, the authors studed the velocty-le state consensus problem n essence n applcatons, one example s the congeston control of the nternet t s desrable that each router coordnates ts data-processng rate to be consstent wth ts neghbors accordng to the amount of date processed n the latest tme, snce velocty coordnaton can greatly reduce the queung delay and pacage loss rate, and enhance the processng effcency Another example s the decomposton of complex tass Each agent also should coordnate ts processng rate to mprove wor effcency The man contrbuton of ths paper s to provde an effectve partal state consensus control strategy, vald for the case wth swtchng nteracton topology, tme-varyng delays and ntermttent nformaton transmsson We frst gve the desgn result of the dstrbuted coordnaton protocol based on delayed poston-le state nformaton Then by usng the tools from graph theory and nonnegatve matrx theory, we show that f the unon of the nteracton topology across the tme nterval wth some gven length always contans a spannng tree, the partal state consensus problem wll be solvable Moreover, the studed system s an asynchronous one, whch means that each agent does not necessarly approxmate ts neghbors velocty-le states for ts local feedbac at the same tme-steps by a global cloc Ths paper s organzed as follows Prelmnary notons are assembled n Secton The consdered problem s formulated n Secton 3 The man result s presented n Secton 4 and ts techncal proof s postponed to Secton 5 Fnally, concludng remars are summarzed n Secton 6 Notatons: Throughout ths paper, let = A = A A A, denotng the left product of matrces, let be the dentty matrx and let = [,,, ] T wth compatble dmensons We wrte A B f A B s nonnegatve Prelmnary notons n graph theory A drected graph G conssts of vertex set V(G = {v, v,, v n } and edge set E(G V(G V(G The edges such as (v, v are called self-loops A path n drected graph G from v to v s a sequence v, v,, v of fnte vertces such that (v l, v l+ E(G for l =,,, Drected graph G s sad to have a spannng tree f there exsts a vertex, called the root, such that t can be connected to any other vertces through paths The unon of a group of drected graphs G,, wth a common vertex set V, s also a drected graph wth the vertex set V and wth the edge set gven by E(G, where s the ndex set of the group A weghted drected graph G(C s a drected graph G together wth a nonnegatve weght matrx C = [c ] R n n such that (v, v E(G c > 0 And n ths case, c s called the weght of edge (v, v 3 Problem formulaton The system studed n ths paper conssts of n autonomous agents, labeled through n All the agents share a common state space R Let x and v denote the poston-le and veloctyle states of agent respectvely and suppose that agent, =,,, n, s wth the followng second-order dynamcs ẋ (t = v (t v (t = u (t where u (t s a local state feedbac, called protocol, to be desgned based on the nformaton receved by agent from ts neghbors Poston-le varables x may represent postons, worloads, etc, and they are the only nformaton transmtted over the networ Suppose that each agent can communcate wth some other agents whch are defned as ts neghbors We use a drected graph G wth vertex set V(G = {v, v,, v n } to represent the communcaton topology Vertex v represents agent and edge (v, v E(G f and only f there exsts an avalable nformaton channel from agent to agent Because of lmted detecton range of agents, exstence of obstacles, or external nterference n sgnals, the communcaton topology s usually dynamcally changng We denote the changng topology by G(t Defnton (Partal State Consensus n the Second-Order Case Gven protocol u, =,,, n, u or ths mult-agent system s sad to solve the velocty-le state consensus problem, that s, a partal state consensus problem, f for any ntal states, there exsts a common asymptotcally stable equlbrum pont v R for all agents wth respect to velocty-le states, such that lm t v (t = v for all Remar n the case of the networs of second-order dynamc agents, governed by Eq (, the poston-le state consensus mples the velocty-le state consensus, and thus the veloctyle state consensus s the only partal state consensus n the strct sense So ths paper focuses on the latter case only 3 Velocty-le state estmaton The obectve of ths paper s to propose an effectve protocol ensurng the solvablty of the partal state consensus problem under relaxable condtons To acheve ths end, we next gve a smple velocty-le state estmaton strategy based on the receved poston-le state nformaton Assume that agent detects ts neghbors poston-le states ntermttently at ts update tmes t, 0 t,, t,, and assume that the update tme sequence satsfes the followng assumpton: (A there exst common lower and upper bounds Ť u, ˆTu for the length of tme nterval between any two consecutve update tmes, such that 0 < Ť u t + t ˆTu for any N and any {,,, n} The above assumpton mples that the update tme sequence s strctly ncreasng, unbounded and wth no fnte accumulaton ponts, namely, lm t = By the propertes of the update tme sequence, the studed system can be classfed nto two categores: the one wth the property that t = t for all,, s called the synchronous system, and the other one wthout the precedng property s called the asynchronous system Ths paper wll focus on the asynchronous case At update tme t, agent may get only some of ts neghbors states because of the exstence of communcaton tme-delays Assume that each agent s equpped wth on-board memory to store nformaton receved from neghbors The nformaton may (

4 F Xao et al / Systems & Control Letters 59 ( be wth tme-delays To estmate the velocty-le states, t s further assumed that the tme-delays are detectable, e, the nformaton s tme-stamped Denote the avalable data of agent at tme t by D (t, whch s composed of the nformaton receved at prevous fnte update tmes To be consstent wth the protocol proposed n the next subsecton, we gve another defnton of neghbors: Defnton (Neghbor Set N E (t and Estmaton of Velocty-Le States Let T E > 0 be gven For any and any update tme t, f there exsts at least one tme-par (t, t such that (A x (t, x (t D (t ; (A3 t T E t < t t, then agent s called a neghbor of agent on the tme nterval [t, t + Denote the neghbor set of agent n ths sense at tme t by N E (t At update tme t, f N E (t, then agent selects one such tme-par (t, t, satsfyng Assumptons (A and (A3, and updates the estmaton about the velocty-le state of agent by v E (t = x (t x (t t t, t [t, t + ( n the above defnton, parameter T E s a pre-gven parameter, nown by all agents and ndependent of parameters Ť u and ˆTu n Assumpton (A, and t represents the maxmum allowable tmedelay, guaranteeng that the data used n the desgned protocol are suffcently new Moreover, the selecton of T E affects the number of elements n the neghbor set N E (t and the structure of nteracton topology G E (t, defned n the next subsecton And t also affects the fnal value of the consensus state n the process of data selecton n approxmatng velocty-le states, t may happen that there exst more than one tme-pars, le (t, t, satsfyng Assumptons (A and (A3 n ths case, agent can choose one tme-par randomly or by the most-recentdata strategy n the latter case, x (t s the most recent postonle state nformaton about agent n D (t and x (t s the most recent poston-le state nformaton satsfyng Assumptons (A and (A3 Obvously, n ths case, f N E (t N E (t +, then t,+ t and t,+ t t can be shown by smulatons that the most-recent-data strategy s more lely to result n a hgher convergence rate Notce that the choce of tme-par (t, t n fact affects the fnal value of consensus state and dfferent agents may choose dfferent polces n a dstrbuted manner for the choce among multple tme-pars, satsfyng Assumptons (A and (A3 Moreover, a lower bound requrement of t t can be added to reduce the effect of nose Remar Note that estmaton v E (t and thus protocol (3 are only dependent on the tme dfference t t and state dsplacement across the tme nterval [t ] Therefore, synchronzaton of all, t agents clocs s not a necessary condton, and t and t can be replaced by the ones decded by the local cloc of agent or accordng to practcal stuatons However, t s requred that all agents should evolve n the same tme scale From ths vewpont, the update tme sequence t, 0 t,, =,,, n, whch s decded by the global cloc, can be seen as the one decded by local clocs of agents, wthout losng the consensus property of protocol (3 Remar f all agents can get ther neghborng agents postonle states wthout tme-delays, then the assumpton of tmestamped nformaton can be removed One example s the velocty consensus control of multple robots, where x (t represents the poston of agent n ths case, we can suppose that at detectng tmes t and t, agent can detect the relatve postons of ts neghborng agent, namely, x (t x (t and x (t x (t, respectvely f the dsplacement of agent over the tme nterval [t, t] s obtanable by agent, then x (t x (t can be gotten by (x (t x (t (x (t x (t + (x (t x (t 3 Partal state consensus protocol Wth the above preparatons, we now propose the followng dstrbuted partal state consensus protocol : u (t = N E (t W (t N E (t W (t (ve (t v (t, t [t, t +, (3 where W (t are called weghtng factors [7], taen from a gven compact set W consstng of postve real numbers Obvously, u (t and v (t are not smooth but they are pecewse dfferentable wth respect to tme t By ths fact and by the asynchrony of update tmes, the dfferental Md-Value Theorem does not hold, n other words, for the estmaton v E (t gven by Eq (, there may not exst t [t, t ] such that v (t = v E (t Ths shows that system ( s not equvalent to the frst-order case wth tme-delays studed n [4] We end ths secton wth a further dscusson on the communcaton topology We now that the communcaton topology G(t only represents the avalable nformaton flow among agents, whereas t does not ndcate whether the neghbors nformaton s used n the feedbac The followng defnton of nteracton topology reflects the relatonshp determned by the nter-usage of nformaton Defnton 3 (nteracton Topology G E (t The vertex set of nteracton topology G E (t s {v, v,, v n }, representng the n agents respectvely, and (v, v E(G E (t f and only f N E (t 4 Convergence result Ths secton presents the convergence result about the system under protocol (3 and ts techncal proof s postponed to the next secton Theorem 4 f there exsts some T > 0, such that for any t 0, the unon of nteracton topology G E (t over tme nterval [t, t+t] always contans a spannng tree, and asynchronous system ( satsfes Assumptons (A (A3, then protocol (3 solves the partal state (velocty-le state consensus problem asymptotcally Remar The suffcent condton provded n the above theorem s a mld and less conservatve one Here, contanng a spannng tree s n fact to guarantee that the nformaton of at least one agent can flow to the entre networs drectly or ndrectly n the lterature, the typcal nteracton graph condtons ensurng solvablty of consensus problems under the assocated proposed protocols can be generally classfed nto three categores The frst s that the nteracton topology s always connected (n the bdrectonal case or always has a spannng tree (n the undrectonal case Ths s the most conservatve one The second s that the perodcal unon of nteracton graph s always Ths paper assumes that f N E (t =, then u (t = 0

5 778 F Xao et al / Systems & Control Letters 59 ( connected or always contans a spannng tree The suffcent condton stated n Theorem 4 belongs to ths category The last one s the unon of all forthcomng nteracton graphs contans a spannng tree and ths s the mldest one for the case under tme-dependent nteracton topology On the other hand, there s a tradeoff between the nteracton graph condton and the basc setup of the studed system Generally speang, the stronger basc assumpton s expected to have a relatvely mlder suffcent nteracton graph condton Tll now, to the best of our nowledge, the only few results that can get the last mldest condton usually concern the specal system n the bdrectonal case, see the wor of Moreau [7] To determne whether the partal state consensus problem s solvable by the topology of nformaton channel G(t, we have the followng corollary: Corollary 5 Assume that the communcaton topology G(t s tmenvarant and contans a spannng tree, and assume that each agent can obtan ts neghbors (determned by G(t poston-le states wth bounded tme-delays for any tme, that s, there exsts a maxmum transmsson tme-delay T max and f there exsts an nformaton channel from agent to agent, then, at any tme t, agent can obtan agent s poston-le state, denoted by x (t, wth the property that 0 t t T max Then there exsts an avalable dstrbuted control rule n the form of protocol (3, satsfyng Assumptons (A (A3, and t solves the partal state (velocty-le state consensus problem asymptotcally Proof To prove the result, t suffces to fnd a possble dstrbuted control rule, satsfyng the assumptons assumed by Theorem 4 Suppose that the maxmum transmsson tme-delay s T max and let the update tme sequence of agent, =,,, n, be t 0, t = t 0 + T max,, t + = t + T max, Then t satsfes Assumpton (A wth Ť u = ˆTu = T max We further assume that agent measures all ts neghbors poston-le state nformaton n the tme ntervals (t, t, =,, For wth (v, v E(G, denote the obtaned nformaton related to agent n tme nterval (t, t by x (t Then t < t < t and thus t > t For t [t, t +, 4, let v E (t = x (t x (t t t < t < t t, we have, satsfy Assumptons (A Then t s well defned Snce t 4 that t and t, n the place of t and t and (A3 wth T E = 4T max Furthermore, the above assumpton also mples that for t max t 4, E(GE (t = E(G(t Therefore, under the above proposed control rule, the system solves the partal state consensus problem asymptotcally 5 Techncal proof Ths secton performs the convergence analyss on asynchronous system ( based on the nonnegatve matrx approach, whch s an effectve way to show the consensus property of multagent systems wth swtchng topology and tme-varyng delays The frst subsecton summarzes some ey lemmas, establshed n [4,5,8] They descrbe the convergence property of the product of a compact set of SA (Stochastc, ndecomposable and Aperodc matrces and gve relaxable condtons ensurng a stochastc matrx to be an SA matrx, respectvely The second subsecton collects all event tmes and merges them nto a sngle ordered tme sequence, denoted by t 0, t,, t, The evoluton of velocty-le states and ther estmaton values s studed wth respect to the above newly defned tme sequence Ths approach s called the Analytc Synchronzaton method n [9] Then we defne a set of new varables v A (t + by (x (t + x (t /(t + t, =,,, n t s shown that the estmated velocty-le states v E ( can be represented by a convex combnaton of varables v A ( By nvestgatng the relatonshp among these varables, we ntroduce an augmented (m + n-dmensonal state varable y(, and then transform the contnuous-tme system ( nto ts equvalent dscrete-tme system (7, where m s a nonnegatve nteger, determned by Assumptons (A and (A3 Thus the partal state consensus problem can be treated as a full state consensus problem equvalently However, the state matrx Ξ( of the dscretetme system has ts specal structure, whch s dfferent from that of the exstng ones nvestgated n the lterature [4,5,7,5,8] Fortunately, by constructng and characterzng two compact sets M and H n the thrd subsecton, whch nclude all possble state matrces of system (7 and all possble products of a fnte number of state matrces at consecutve tme-steps of system (7, respectvely, we can apply the ey lemmas, presented n the frst subsecton, and get the convergence result Here, we emphasze that although the proof steps are smlar to that taen n [4], n other words, they are all based on the presented ey lemmas, the dfferences between the contrbutons of the two papers are also obvous Frst, they study two dstnct nds of problems The feasblty of the proof wth the help of Lemmas 6 8 owes much to the sllful choce of state varable y( Second, the proof detals are also dfferent Ths paper gves the only arguments that are needed to be clarfed and omts the obvous ones that can be learnt from other lterature 5 Key lemmas n ths subsecton, we frst gve the defnton of SA matrx and then lst three mportant lemmas, whch are useful n provng the man result A stochastc matrx A s called ndecomposable and aperodc (SA f there exsts a column vector ν such that lm A = ν T Lemma 6 ([5, Lemma 5] Let A be a compact set, consstng of n n SA matrces f for any and any A, A,, A A (repettons permtted, = A s SA, then for any gven nfnte matrx sequence A, A,, A, (repettons permtted, there exsts a column vector ν such that lm A = ν T = Lemma 7 ([8, Lemma ] Let A be a stochastc matrx f G(A contans a spannng tree wth the property that the assocated root vertex has a self-loop n G(A, then A s SA Lemma 8 ([4, Lemma 8] Let A 0, A,, A m be n n nonnegatve matrces, let A 0 A A m D = (m+n (m+n let 0 Q 0 = 0 0 (m+n (m+n

6 and let Q = D + Q 0 for any {,,, m} f G( m = A contans a spannng tree, then G(Q also contans a spannng tree wth the property that the assocated root vertex has a self-loop n G(Q 5 Equvalent representaton n ths subsecton, we employ the Analytc Synchronzaton method and get the dscrete-tme Eq (7, that s, an equvalent representaton of the orgnal contnuous-tme system The basc dea of Analytc Synchronzaton s to study the asynchronous system by usng a sutably defned dscrete-tme synchronous system, evolvng on the collecton of event tmes of all orgnal subsystems [9] Frst, for symbolc smplcty, we generalze the defnton of v E (t and ntroduce a weght matrx A(t = [a (t] R n n, where v E (t s generalzed by f N E (t, then v E (t s defned by Eq (; f =, then v E (t = v (t; 3 n other cases, v E (t = 0, and A(t s defned by f N E (t, W (t a (t = W s (t, N E (t s N E (t 0, otherwse, f N E (t =,, = a (t = 0, otherwse, t [t, t + F Xao et al / Systems & Control Letters 59 ( t [t, t + By the above defnton, v E (t,, s a pece-wse constant functon of tme t, and f gnore the weght of each edge and selfloops n G(A(t, then G(A(t and G E (t represent the same nteracton topology Notcng that W (t W, we get that all possble A(t consttute a compact set A and ther nonnegatve entres are not less than mn{w : w W}/((n max{w : w W} Now collect all tme {t, t, t : =,,, n, N E (t, N} and relabel the nonnegatve elements of them by t 0, t,, t,, n ncreasng order such that t 0 = 0, t < t + Here, we assume that t 0 = 0 for all {,,, n} ndeed, wthout ths assumpton, we can consder the dynamcs of agents after tme max t 0 and get the same convergence result For smplcty, denote t + t by τ and denote (x (t + x (t /(t + t by v A (t + (superscrpt A means average velocty, for N, =,,, n Next we study the evoluton of varables v (t and v A (t wth respect to We wll prove that ther reachng an agreement mples the solvablty of the partal state consensus problem Solvng Eq (3 gves that n v (t = e (t t v (t + ( e (t t a (t ve (t x (t x (t t t = e (t t t t + v (t e (t t t t = n = a (t ve (t where t < t t + t can be observed that 0 < ( e (t t /(t t < for any t < t t + The above equaton mples that (4 n v (t + = e τ v (t + ( e τ = v A (t + = e τ v (t τ + e τ n τ = a (t s v E (t s a (t s v E (t s where s N such that t s t < t + t s + Ths part gves an expresson v E (t s n Eq (5 n terms of v A ( Suppose that N E (t s, t s = t l and t s = t p By Assumptons (A and (A3, there exsts an m N, ndependent of,, (cf Lemma n [4], such that t m t l < t p t s for m Then v E (t s = x (t p x (t l t p t l = τ p v A (t p + τ p v A (t p + + τ l v A (t l + τ p + τ p + + τ l (6 whch means that v E (t s s a convex combnaton of v A (t l +, v A (t l +,, v A (t p From the fact that p l m, t follows that some of ts coeffcents are not less than /m To represent the evoluton Eq (5 n matrx form, we ntroduce the augmented state varable y( = [v(t T, v A (t T, v A (t T,, v A (t m+ T ] T for m, where v(t = [v (t, v (t,, v n (t ] T and v A (t = [v A (t, v A (t,, v A n (t ] T Combnng Eqs (5 and (6 yelds that y( + = Ξ(y(, (7 where Ξ( R (m+n (m+n s defned n Box and A γ ( = ], γ {0,,, m}, are defned by [a γ f N E (t s, τ γ a (t s, a γ τ = p + τ p + + τ l N E (t s, p + γ l 0, otherwse f N E (t s =, a γ, =, γ = 0 = 0, otherwse n the followng, the matrx gven n Box, wth the above structure, s denoted by M(τ, A 0 (, A (,, A m ( to ndcate ts dependence on parameters τ, A 0 (, A (,, A m ( By the above defnton, t can be easly obtaned that m Lemma 9 γ =0 A γ ( = A(t and thus Ξ( s stochastc; f N E (t, there exsts γ m, such that aγ mn{w : w W}/(m(n max{w : w W} Remar Eq (7 can be seen as the dynamcal equaton of a dscrete-tme mult-agent system consstng of (m + n agents under the tme-varyng nteracton topology G(Ξ(, and varable y( s the column vector, staced wth the states of the (m + n agents These states are the veloctyle states v (t, =,,, n, and the average velocty v A (t, v A (t,, v A (t m+, =,,, n n Lemma, we wll show that ther convergence to a consensus state as wll lead to that v (t, =,,, n, reach consensus asymptotcally as t Consensus problems of dscretetme mult-agent systems were wdely studed by researchers (5

7 780 F Xao et al / Systems & Control Letters 59 ( e τ + ( e τ A0 ( ( e τ A ( ( e τ Am ( ( e τ Am ( e τ τ + ( e τ A τ 0 ( ( e τ A τ ( ( e τ A τ m ( ( e τ A τ m ( Ξ( = 0 0 Box However, ths dscrete-tme model cannot be covered by the exstng ones, such as those studed n [7,7,0], because the dagonal entres of state matrx Ξ( are not all larger than zero and parameters τ may tae any value n (0, ˆTu ], see [4] for detaled dscussons Defnton 0 (Full State Consensus The dscrete-tme system (7 s sad to solve a (full state consensus problem f for any ntal state, there exsts an asymptotcally stable equlbrum pont y, y R, such that lm y( = y, n other words, lm v (t = lm v A (t = y, =,,, n The followng lemma states the relatonshp between contnuous-tme system ( and dscrete-tme system (7 Lemma System ( solves a partal state (velocty-le state consensus problem f and only f system (7 solves a full state consensus problem Proof The necessty s obvous and we only prove the suffcency Suppose that lm v (t = lm v A (t = v for all =,,, n Thus by Eqs (4 and (6 and Assumpton (A, lm t v (t = v The next subsecton wll prove Theorem 4 by showng system (7 solves a full state consensus problem 53 Convergence analyss of system (7 Ths subsecton conssts of three parts The frst part gves an equvalent representaton of the condton assumed n Theorem 4 The second part characterzes the propertes of state matrx Ξ( of system (7 by two compact matrx sets M and H The last part gves the proof of Theorem 4 Frst, the followng lemma restates the condton provded n Theorem 4 n an equvalent way Lemma f there exsts T > 0, such that for any t 0, the unon of nteracton topology G E (t across the tme nterval [t, t + T] always contans a spannng tree, then there exst a postve nteger h and a postve real number T h wth the followng property: for any N, there exsts a subset of {t, t +,, t +h }, denoted by T, such that the unon of G E (t on T contans a spannng tree, and for any t l T, T h τ l ˆTu Proof t s a straghtforward consequence of Assumptons (A and (A3 and the detals are omtted, see Lemma 9 n [4] for smlar dscussons n order to mae use of Lemma 6 to perform the convergence analyss, we next construct two compact sets M and H The frst compact set M ncludes all possble state matrces of system (7, whch s defned by M = M(ς, N 0, N,, N m : 0 ς ˆTu, N 0, N,, N m are nonnegatve, and there exsts some A A, such that N 0 + N + + N m = A where we use the conventon ( e ς /ς ς=0 = lm ς 0 ( e ς /ς = The second compact set H ncludes all possble products of h state matrces at consecutve tme-steps of system (7, whch s defned by h H = M(ς, N, 0 N,, N : m = M(ς, M, and there exsts a subset of {,,, h}, denoted by T, such that for any m T, T h ς ˆTu and G contans a spannng tree T where h and T h are gven n Lemma Lemma 3 M and H are compact sets, and for any m, Ξ( M +h and l= Ξ(l H; for any H H, H s SA, and for any N, f H, H,, H H (repettons permtted, then = H s SA Proof (a The compactness of M follows from the fact that set A and nterval [0, Ť u ] are compact; the compactness of H follows from that M s compact; all possble choces of T and the spannng tree are fnte; 3 the product of fnte matrces s a contnuous functon; 4 all the nonnegatve entres of matrces n A are lower bounded by mn{w : w W}/((n max{w : w W} (b The concluson that Ξ( M follows drectly from the defntons of Ξ( and M, and Lemma 9 By Lemmas 9 and, G( m l T γ =0 (A γ (l contans a spannng tree, and thus +h l= Ξ(l H (c Let H = h M(ς =, N, 0 N,, N m and T be the assocated subset of {,,, h} such that for any T, T h ς ˆTu and G m T =0 N contans a spannng tree Let Q0 be the same as the Q 0 n Lemma 8, let N 0 N N m D = and let ε = nf{e ς, ( e ς /ς : ς (0, ˆTu ]} Then ε > 0 and h M(ς, = =0 h εq0 + ( e ς D = ε h Q 0 h + ε h N h ( e ς Q h 0 D Q 0 = mn{ε h, ε h ( e T h } Q 0 h + T D Q 0 h

8 F Xao et al / Systems & Control Letters 59 ( where the last nequalty follows from Q 0 D D and from the fact that T h ς ˆTu for T Let the frst n rows of D Q h 0 be B, 0 B,, m B, where B R n n, = 0,,, m m Then = B = m =0 N and thus m T =0 B = m T =0 N Snce G m T =0 N contans a spannng tree, G m T =0 B also contans a spannng tree Let D R (m+n (m+n, wth the same frst n rows as T D Q 0 h and wth the other rows beng zeros Then h M(ς, mn{ε h, ε h ( e T h }(Q h 0 + D = By Lemma 8 and the fact that Q 0 = Q 0 m f m, G(Q 0 h + D contans a spannng tree wth the property that the assocated root vertex has a self-loop, and so s G(H Snce H s stochastc, by Lemma 7, H s SA By the same arguments, we have that for any H, H,, H H, = H s SA Proof of Theorem 4 Ths part only proves that system (7 solves a full state consensus problem t follows from Lemmas 6 and 3 that there exsts ν R (m+n such that ph lm p Ξ(m + l = ν T 3 (8 For any N, there exsts p such that p h < (p + h And snce matrx Ξ(m + l s stochastc, lm Ξ(m + l ν T = lm Ξ(m + l l=p h p h Ξ(m + l ν T Snce M s compact, we get that matrx l=p h Ξ(m + l n the above equaton belongs to a bounded set Furthermore, t follows from Eq (8 that lm p h p h Ξ(m + l ν T = lm p Thus, lm Ξ(m + l ν T = 0, whch yelds that lm y( = lm Ξ(m + ly(m = ν T y(m, Ξ(m + l ν T = 0 that s, system (7 solves a full state consensus problem 6 Concluson Ths paper consdered the partal state consensus problem n networs of second-order agents wth dynamcally changng nteracton topologes and tme-varyng communcaton tmedelays and proposed an asynchronous partal state consensus protocol, whose valdty can be guaranteed under the relaxable condton that poston-le states as the only nformaton are transmtted ntermttently over the networ Nevertheless, there stll exst some other nterestng problems that need to be addressed, such as the desgn and analyss of a full state consensus protocol n the framewor of ths paper References [] R Olfat-Saber, JA Fax, RM Murray, Consensus and cooperaton n networed mult-agent systems, Proceedngs of the EEE 95 ( [] W Ren, Mult-vehcle consensus wth a tme-varyng reference state, Systems & Control Letters 56 ( [3] Z J, Z Wang, H Ln, Z Wan, nterconnecton topologes for multagent coordnaton under leader follower framewor, Automatca 45 ( [4] F Xao, L Wang, Asynchronous consensus n contnuous-tme mult-agent systems wth swtchng topology and tme-varyng delays, EEE Transactons on Automatc Control 53 ( [5] F Xao, L Wang, Consensus protocols for dscrete-tme mult-agent systems wth tme-varyng delays, Automatca 44 ( [6] R Olfat-Saber, RM Murray, Consensus problems n networs of agents wth swtchng topology and tme-delays, EEE Transactons on Automatc Control 49 ( [7] W Ren, RW Beard, Consensus seeng n multagent systems under dynamcally changng nteracton topologes, EEE Transactons on Automatc Control 50 ( [8] G Xe, L Wang, Consensus control for a class of networs of dynamc agents, nternatonal Journal of Robust and Nonlnear Control 7 ( [9] Y Hong, J Hu, L Gao, Tracng control for mult-agent consensus wth an actve leader and varable topology, Automatca 4 ( [0] W Ren, On consensus algorthms for double-ntegrator dynamcs, EEE Transactons on Automatc Control 53 ( [] W Ren, Collectve moton from consensus wth Cartesan coordnate couplng, EEE Transactons on Automatc Control 54 ( [] W Ren, Synchronzaton of coupled harmonc oscllators wth local nteracton, Automatca 44 ( [3] X Wang, Y Hong, Fnte-tme consensus for mult-agent networs wth second-order agent dynamcs, n: Proceedngs of the 7th World Congress, The nternatonal Federaton of Automatc Control, 008, pp [4] T Vcse, A Czro, E Ben-Jacob, Cohen, O Schochet, Novel type of phasetranston n a system of self-drven partcles, Physcal Revew Letters 75 ( [5] A Jadbabae, J Ln, AS Morse, Coordnaton of groups of moble autonomous agents usng nearest neghbor rules, EEE Transactons on Automatc Control 48 ( [6] S Barbarossa, G Scutar, Decentralzed maxmum-lelhood estmaton for sensor networs composed of nonlnearly coupled dynamcal systems, EEE Transactons on Sgnal Processng 55 ( [7] L Moreau, Stablty of multagent systems wth tme-dependent communcaton lns, EEE Transactons on Automatc Control 50 ( [8] F Xao, L Wang, State consensus for mult-agent systems wth swtchng topologes and tme-varyng delays, nternatonal Journal of Control 79 ( [9] J Ln, AS Morse, BDO Anderson, The mult-agent rendezvous problem Part : the asynchronous case, SAM Journal on Control and Optmzaton 46 ( [0] D Angel, PA Blman, Stablty of leaderless dscrete-tme mult-agent systems, Mathematcs of Control, Sgnals, and Systems 8 ( f <, then l= Ξ(l =

Changing Topology and Communication Delays

Changing Topology and Communication Delays Prepared by F.L. Lews Updated: Saturday, February 3, 00 Changng Topology and Communcaton Delays Changng Topology The graph connectvty or topology may change over tme. Let G { G, G,, G M } wth M fnte be

More information

Dynamic Systems on Graphs

Dynamic Systems on Graphs Prepared by F.L. Lews Updated: Saturday, February 06, 200 Dynamc Systems on Graphs Control Graphs and Consensus A network s a set of nodes that collaborates to acheve what each cannot acheve alone. A network,

More information

More metrics on cartesian products

More metrics on cartesian products More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of

More information

Adaptive Consensus Control of Multi-Agent Systems with Large Uncertainty and Time Delays *

Adaptive Consensus Control of Multi-Agent Systems with Large Uncertainty and Time Delays * Journal of Robotcs, etworkng and Artfcal Lfe, Vol., o. (September 04), 5-9 Adaptve Consensus Control of Mult-Agent Systems wth Large Uncertanty and me Delays * L Lu School of Mechancal Engneerng Unversty

More information

Affine transformations and convexity

Affine transformations and convexity Affne transformatons and convexty The purpose of ths document s to prove some basc propertes of affne transformatons nvolvng convex sets. Here are a few onlne references for background nformaton: http://math.ucr.edu/

More information

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could

More information

Outline. Communication. Bellman Ford Algorithm. Bellman Ford Example. Bellman Ford Shortest Path [1]

Outline. Communication. Bellman Ford Algorithm. Bellman Ford Example. Bellman Ford Shortest Path [1] DYNAMIC SHORTEST PATH SEARCH AND SYNCHRONIZED TASK SWITCHING Jay Wagenpfel, Adran Trachte 2 Outlne Shortest Communcaton Path Searchng Bellmann Ford algorthm Algorthm for dynamc case Modfcatons to our algorthm

More information

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number

More information

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons

More information

Lecture 12: Discrete Laplacian

Lecture 12: Discrete Laplacian Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly

More information

The Order Relation and Trace Inequalities for. Hermitian Operators

The Order Relation and Trace Inequalities for. Hermitian Operators Internatonal Mathematcal Forum, Vol 3, 08, no, 507-57 HIKARI Ltd, wwwm-hkarcom https://doorg/0988/mf088055 The Order Relaton and Trace Inequaltes for Hermtan Operators Y Huang School of Informaton Scence

More information

College of Computer & Information Science Fall 2009 Northeastern University 20 October 2009

College of Computer & Information Science Fall 2009 Northeastern University 20 October 2009 College of Computer & Informaton Scence Fall 2009 Northeastern Unversty 20 October 2009 CS7880: Algorthmc Power Tools Scrbe: Jan Wen and Laura Poplawsk Lecture Outlne: Prmal-dual schema Network Desgn:

More information

APPENDIX A Some Linear Algebra

APPENDIX A Some Linear Algebra APPENDIX A Some Lnear Algebra The collecton of m, n matrces A.1 Matrces a 1,1,..., a 1,n A = a m,1,..., a m,n wth real elements a,j s denoted by R m,n. If n = 1 then A s called a column vector. Smlarly,

More information

Errors for Linear Systems

Errors for Linear Systems Errors for Lnear Systems When we solve a lnear system Ax b we often do not know A and b exactly, but have only approxmatons  and ˆb avalable. Then the best thng we can do s to solve ˆx ˆb exactly whch

More information

A new construction of 3-separable matrices via an improved decoding of Macula s construction

A new construction of 3-separable matrices via an improved decoding of Macula s construction Dscrete Optmzaton 5 008 700 704 Contents lsts avalable at ScenceDrect Dscrete Optmzaton journal homepage: wwwelsevercom/locate/dsopt A new constructon of 3-separable matrces va an mproved decodng of Macula

More information

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle

More information

Distributed Exponential Formation Control of Multiple Wheeled Mobile Robots

Distributed Exponential Formation Control of Multiple Wheeled Mobile Robots Proceedngs of the Internatonal Conference of Control, Dynamc Systems, and Robotcs Ottawa, Ontaro, Canada, May 15-16 214 Paper No. 46 Dstrbuted Exponental Formaton Control of Multple Wheeled Moble Robots

More information

Difference Equations

Difference Equations Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1

More information

Foundations of Arithmetic

Foundations of Arithmetic Foundatons of Arthmetc Notaton We shall denote the sum and product of numbers n the usual notaton as a 2 + a 2 + a 3 + + a = a, a 1 a 2 a 3 a = a The notaton a b means a dvdes b,.e. ac = b where c s an

More information

DO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED.

DO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED. EE 539 Homeworks Sprng 08 Updated: Tuesday, Aprl 7, 08 DO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED. For full credt, show all work. Some problems requre hand calculatons.

More information

Numerical Heat and Mass Transfer

Numerical Heat and Mass Transfer Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and

More information

Problem Set 9 Solutions

Problem Set 9 Solutions Desgn and Analyss of Algorthms May 4, 2015 Massachusetts Insttute of Technology 6.046J/18.410J Profs. Erk Demane, Srn Devadas, and Nancy Lynch Problem Set 9 Solutons Problem Set 9 Solutons Ths problem

More information

Perfect Competition and the Nash Bargaining Solution

Perfect Competition and the Nash Bargaining Solution Perfect Competton and the Nash Barganng Soluton Renhard John Department of Economcs Unversty of Bonn Adenauerallee 24-42 53113 Bonn, Germany emal: rohn@un-bonn.de May 2005 Abstract For a lnear exchange

More information

MMA and GCMMA two methods for nonlinear optimization

MMA and GCMMA two methods for nonlinear optimization MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons

More information

Linear Approximation with Regularization and Moving Least Squares

Linear Approximation with Regularization and Moving Least Squares Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...

More information

The lower and upper bounds on Perron root of nonnegative irreducible matrices

The lower and upper bounds on Perron root of nonnegative irreducible matrices Journal of Computatonal Appled Mathematcs 217 (2008) 259 267 wwwelsevercom/locate/cam The lower upper bounds on Perron root of nonnegatve rreducble matrces Guang-Xn Huang a,, Feng Yn b,keguo a a College

More information

princeton univ. F 17 cos 521: Advanced Algorithm Design Lecture 7: LP Duality Lecturer: Matt Weinberg

princeton univ. F 17 cos 521: Advanced Algorithm Design Lecture 7: LP Duality Lecturer: Matt Weinberg prnceton unv. F 17 cos 521: Advanced Algorthm Desgn Lecture 7: LP Dualty Lecturer: Matt Wenberg Scrbe: LP Dualty s an extremely useful tool for analyzng structural propertes of lnear programs. Whle there

More information

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng

More information

Cooperative Output Regulation of Linear Multi-agent Systems with Communication Constraints

Cooperative Output Regulation of Linear Multi-agent Systems with Communication Constraints 2016 IEEE 55th Conference on Decson and Control (CDC) ARIA Resort & Casno December 12-14, 2016, Las Vegas, USA Cooperatve Output Regulaton of Lnear Mult-agent Systems wth Communcaton Constrants Abdelkader

More information

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:

More information

Maximizing the number of nonnegative subsets

Maximizing the number of nonnegative subsets Maxmzng the number of nonnegatve subsets Noga Alon Hao Huang December 1, 213 Abstract Gven a set of n real numbers, f the sum of elements of every subset of sze larger than k s negatve, what s the maxmum

More information

Improved delay-dependent stability criteria for discrete-time stochastic neural networks with time-varying delays

Improved delay-dependent stability criteria for discrete-time stochastic neural networks with time-varying delays Avalable onlne at www.scencedrect.com Proceda Engneerng 5 ( 4456 446 Improved delay-dependent stablty crtera for dscrete-tme stochastc neural networs wth tme-varyng delays Meng-zhuo Luo a Shou-mng Zhong

More information

EEE 241: Linear Systems

EEE 241: Linear Systems EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they

More information

Convexity preserving interpolation by splines of arbitrary degree

Convexity preserving interpolation by splines of arbitrary degree Computer Scence Journal of Moldova, vol.18, no.1(52), 2010 Convexty preservng nterpolaton by splnes of arbtrary degree Igor Verlan Abstract In the present paper an algorthm of C 2 nterpolaton of dscrete

More information

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4) I. Classcal Assumptons Econ7 Appled Econometrcs Topc 3: Classcal Model (Studenmund, Chapter 4) We have defned OLS and studed some algebrac propertes of OLS. In ths topc we wll study statstcal propertes

More information

Asymptotics of the Solution of a Boundary Value. Problem for One-Characteristic Differential. Equation Degenerating into a Parabolic Equation

Asymptotics of the Solution of a Boundary Value. Problem for One-Characteristic Differential. Equation Degenerating into a Parabolic Equation Nonl. Analyss and Dfferental Equatons, ol., 4, no., 5 - HIKARI Ltd, www.m-har.com http://dx.do.org/.988/nade.4.456 Asymptotcs of the Soluton of a Boundary alue Problem for One-Characterstc Dfferental Equaton

More information

MATH 241B FUNCTIONAL ANALYSIS - NOTES EXAMPLES OF C ALGEBRAS

MATH 241B FUNCTIONAL ANALYSIS - NOTES EXAMPLES OF C ALGEBRAS MATH 241B FUNCTIONAL ANALYSIS - NOTES EXAMPLES OF C ALGEBRAS These are nformal notes whch cover some of the materal whch s not n the course book. The man purpose s to gve a number of nontrval examples

More information

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction ECONOMICS 5* -- NOTE (Summary) ECON 5* -- NOTE The Multple Classcal Lnear Regresson Model (CLRM): Specfcaton and Assumptons. Introducton CLRM stands for the Classcal Lnear Regresson Model. The CLRM s also

More information

Physics 5153 Classical Mechanics. Principle of Virtual Work-1

Physics 5153 Classical Mechanics. Principle of Virtual Work-1 P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal

More information

Yong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 )

Yong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 ) Kangweon-Kyungk Math. Jour. 4 1996), No. 1, pp. 7 16 AN ITERATIVE ROW-ACTION METHOD FOR MULTICOMMODITY TRANSPORTATION PROBLEMS Yong Joon Ryang Abstract. The optmzaton problems wth quadratc constrants often

More information

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION Advanced Mathematcal Models & Applcatons Vol.3, No.3, 2018, pp.215-222 ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EUATION

More information

Valuated Binary Tree: A New Approach in Study of Integers

Valuated Binary Tree: A New Approach in Study of Integers Internatonal Journal of Scentfc Innovatve Mathematcal Research (IJSIMR) Volume 4, Issue 3, March 6, PP 63-67 ISS 347-37X (Prnt) & ISS 347-34 (Onlne) wwwarcournalsorg Valuated Bnary Tree: A ew Approach

More information

U.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017

U.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017 U.C. Berkeley CS94: Beyond Worst-Case Analyss Handout 4s Luca Trevsan September 5, 07 Summary of Lecture 4 In whch we ntroduce semdefnte programmng and apply t to Max Cut. Semdefnte Programmng Recall that

More information

Power law and dimension of the maximum value for belief distribution with the max Deng entropy

Power law and dimension of the maximum value for belief distribution with the max Deng entropy Power law and dmenson of the maxmum value for belef dstrbuton wth the max Deng entropy Bngy Kang a, a College of Informaton Engneerng, Northwest A&F Unversty, Yanglng, Shaanx, 712100, Chna. Abstract Deng

More information

Winter 2008 CS567 Stochastic Linear/Integer Programming Guest Lecturer: Xu, Huan

Winter 2008 CS567 Stochastic Linear/Integer Programming Guest Lecturer: Xu, Huan Wnter 2008 CS567 Stochastc Lnear/Integer Programmng Guest Lecturer: Xu, Huan Class 2: More Modelng Examples 1 Capacty Expanson Capacty expanson models optmal choces of the tmng and levels of nvestments

More information

Time-Varying Systems and Computations Lecture 6

Time-Varying Systems and Computations Lecture 6 Tme-Varyng Systems and Computatons Lecture 6 Klaus Depold 14. Januar 2014 The Kalman Flter The Kalman estmaton flter attempts to estmate the actual state of an unknown dscrete dynamcal system, gven nosy

More information

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0 MODULE 2 Topcs: Lnear ndependence, bass and dmenson We have seen that f n a set of vectors one vector s a lnear combnaton of the remanng vectors n the set then the span of the set s unchanged f that vector

More information

NP-Completeness : Proofs

NP-Completeness : Proofs NP-Completeness : Proofs Proof Methods A method to show a decson problem Π NP-complete s as follows. (1) Show Π NP. (2) Choose an NP-complete problem Π. (3) Show Π Π. A method to show an optmzaton problem

More information

Société de Calcul Mathématique SA

Société de Calcul Mathématique SA Socété de Calcul Mathématque SA Outls d'ade à la décson Tools for decson help Probablstc Studes: Normalzng the Hstograms Bernard Beauzamy December, 202 I. General constructon of the hstogram Any probablstc

More information

Composite Hypotheses testing

Composite Hypotheses testing Composte ypotheses testng In many hypothess testng problems there are many possble dstrbutons that can occur under each of the hypotheses. The output of the source s a set of parameters (ponts n a parameter

More information

NUMERICAL DIFFERENTIATION

NUMERICAL DIFFERENTIATION NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the

More information

The Two-scale Finite Element Errors Analysis for One Class of Thermoelastic Problem in Periodic Composites

The Two-scale Finite Element Errors Analysis for One Class of Thermoelastic Problem in Periodic Composites 7 Asa-Pacfc Engneerng Technology Conference (APETC 7) ISBN: 978--6595-443- The Two-scale Fnte Element Errors Analyss for One Class of Thermoelastc Problem n Perodc Compostes Xaoun Deng Mngxang Deng ABSTRACT

More information

A Local Variational Problem of Second Order for a Class of Optimal Control Problems with Nonsmooth Objective Function

A Local Variational Problem of Second Order for a Class of Optimal Control Problems with Nonsmooth Objective Function A Local Varatonal Problem of Second Order for a Class of Optmal Control Problems wth Nonsmooth Objectve Functon Alexander P. Afanasev Insttute for Informaton Transmsson Problems, Russan Academy of Scences,

More information

For now, let us focus on a specific model of neurons. These are simplified from reality but can achieve remarkable results.

For now, let us focus on a specific model of neurons. These are simplified from reality but can achieve remarkable results. Neural Networks : Dervaton compled by Alvn Wan from Professor Jtendra Malk s lecture Ths type of computaton s called deep learnng and s the most popular method for many problems, such as computer vson

More information

Discrete Mathematics. Laplacian spectral characterization of some graphs obtained by product operation

Discrete Mathematics. Laplacian spectral characterization of some graphs obtained by product operation Dscrete Mathematcs 31 (01) 1591 1595 Contents lsts avalable at ScVerse ScenceDrect Dscrete Mathematcs journal homepage: www.elsever.com/locate/dsc Laplacan spectral characterzaton of some graphs obtaned

More information

Stability and Stabilization for Discrete Systems with Time-varying Delays Based on the Average Dwell-time Method

Stability and Stabilization for Discrete Systems with Time-varying Delays Based on the Average Dwell-time Method Proceedngs of the 29 IEEE Internatonal Conference on Systems, an, and Cybernetcs San Antono, TX, USA - October 29 Stablty and Stablzaton for Dscrete Systems wth Tme-varyng Delays Based on the Average Dwell-tme

More information

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2 Salmon: Lectures on partal dfferental equatons 5. Classfcaton of second-order equatons There are general methods for classfyng hgher-order partal dfferental equatons. One s very general (applyng even to

More information

Notes on Frequency Estimation in Data Streams

Notes on Frequency Estimation in Data Streams Notes on Frequency Estmaton n Data Streams In (one of) the data streamng model(s), the data s a sequence of arrvals a 1, a 2,..., a m of the form a j = (, v) where s the dentty of the tem and belongs to

More information

Computing Correlated Equilibria in Multi-Player Games

Computing Correlated Equilibria in Multi-Player Games Computng Correlated Equlbra n Mult-Player Games Chrstos H. Papadmtrou Presented by Zhanxang Huang December 7th, 2005 1 The Author Dr. Chrstos H. Papadmtrou CS professor at UC Berkley (taught at Harvard,

More information

General viscosity iterative method for a sequence of quasi-nonexpansive mappings

General viscosity iterative method for a sequence of quasi-nonexpansive mappings Avalable onlne at www.tjnsa.com J. Nonlnear Sc. Appl. 9 (2016), 5672 5682 Research Artcle General vscosty teratve method for a sequence of quas-nonexpansve mappngs Cuje Zhang, Ynan Wang College of Scence,

More information

Kernel Methods and SVMs Extension

Kernel Methods and SVMs Extension Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general

More information

Convergence of random processes

Convergence of random processes DS-GA 12 Lecture notes 6 Fall 216 Convergence of random processes 1 Introducton In these notes we study convergence of dscrete random processes. Ths allows to characterze phenomena such as the law of large

More information

Norm Bounds for a Transformed Activity Level. Vector in Sraffian Systems: A Dual Exercise

Norm Bounds for a Transformed Activity Level. Vector in Sraffian Systems: A Dual Exercise ppled Mathematcal Scences, Vol. 4, 200, no. 60, 2955-296 Norm Bounds for a ransformed ctvty Level Vector n Sraffan Systems: Dual Exercse Nkolaos Rodousaks Department of Publc dmnstraton, Panteon Unversty

More information

Economics 101. Lecture 4 - Equilibrium and Efficiency

Economics 101. Lecture 4 - Equilibrium and Efficiency Economcs 0 Lecture 4 - Equlbrum and Effcency Intro As dscussed n the prevous lecture, we wll now move from an envronment where we looed at consumers mang decsons n solaton to analyzng economes full of

More information

Introductory Cardinality Theory Alan Kaylor Cline

Introductory Cardinality Theory Alan Kaylor Cline Introductory Cardnalty Theory lan Kaylor Clne lthough by name the theory of set cardnalty may seem to be an offshoot of combnatorcs, the central nterest s actually nfnte sets. Combnatorcs deals wth fnte

More information

Inductance Calculation for Conductors of Arbitrary Shape

Inductance Calculation for Conductors of Arbitrary Shape CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors

More information

Graph Reconstruction by Permutations

Graph Reconstruction by Permutations Graph Reconstructon by Permutatons Perre Ille and Wllam Kocay* Insttut de Mathémathques de Lumny CNRS UMR 6206 163 avenue de Lumny, Case 907 13288 Marselle Cedex 9, France e-mal: lle@ml.unv-mrs.fr Computer

More information

Lecture 4. Instructor: Haipeng Luo

Lecture 4. Instructor: Haipeng Luo Lecture 4 Instructor: Hapeng Luo In the followng lectures, we focus on the expert problem and study more adaptve algorthms. Although Hedge s proven to be worst-case optmal, one may wonder how well t would

More information

Consensus of Multi-Agent Systems by Distributed Event-Triggered Control

Consensus of Multi-Agent Systems by Distributed Event-Triggered Control Preprnts of the 19th World Congress The Internatonal Federaton of Automatc Control Consensus of Mult-Agent Systems by Dstrbuted Event-Trggered Control Wenfeng Hu, Lu Lu, Gang Feng Department of Mechancal

More information

Random Walks on Digraphs

Random Walks on Digraphs Random Walks on Dgraphs J. J. P. Veerman October 23, 27 Introducton Let V = {, n} be a vertex set and S a non-negatve row-stochastc matrx (.e. rows sum to ). V and S defne a dgraph G = G(V, S) and a drected

More information

Design and Optimization of Fuzzy Controller for Inverse Pendulum System Using Genetic Algorithm

Design and Optimization of Fuzzy Controller for Inverse Pendulum System Using Genetic Algorithm Desgn and Optmzaton of Fuzzy Controller for Inverse Pendulum System Usng Genetc Algorthm H. Mehraban A. Ashoor Unversty of Tehran Unversty of Tehran h.mehraban@ece.ut.ac.r a.ashoor@ece.ut.ac.r Abstract:

More information

U.C. Berkeley CS294: Spectral Methods and Expanders Handout 8 Luca Trevisan February 17, 2016

U.C. Berkeley CS294: Spectral Methods and Expanders Handout 8 Luca Trevisan February 17, 2016 U.C. Berkeley CS94: Spectral Methods and Expanders Handout 8 Luca Trevsan February 7, 06 Lecture 8: Spectral Algorthms Wrap-up In whch we talk about even more generalzatons of Cheeger s nequaltes, and

More information

Digital Signal Processing

Digital Signal Processing Dgtal Sgnal Processng Dscrete-tme System Analyss Manar Mohasen Offce: F8 Emal: manar.subh@ut.ac.r School of IT Engneerng Revew of Precedent Class Contnuous Sgnal The value of the sgnal s avalable over

More information

Appendix B. The Finite Difference Scheme

Appendix B. The Finite Difference Scheme 140 APPENDIXES Appendx B. The Fnte Dfference Scheme In ths appendx we present numercal technques whch are used to approxmate solutons of system 3.1 3.3. A comprehensve treatment of theoretcal and mplementaton

More information

Linear Regression Analysis: Terminology and Notation

Linear Regression Analysis: Terminology and Notation ECON 35* -- Secton : Basc Concepts of Regresson Analyss (Page ) Lnear Regresson Analyss: Termnology and Notaton Consder the generc verson of the smple (two-varable) lnear regresson model. It s represented

More information

The Quadratic Trigonometric Bézier Curve with Single Shape Parameter

The Quadratic Trigonometric Bézier Curve with Single Shape Parameter J. Basc. Appl. Sc. Res., (3541-546, 01 01, TextRoad Publcaton ISSN 090-4304 Journal of Basc and Appled Scentfc Research www.textroad.com The Quadratc Trgonometrc Bézer Curve wth Sngle Shape Parameter Uzma

More information

Distributed Multi-Agent Coordination: A Comparison Lemma Based Approach

Distributed Multi-Agent Coordination: A Comparison Lemma Based Approach 2011 Amercan Control Conference on O'Farrell Street, San Francsco, CA, USA June 29 - July 01, 2011 Dstrbuted Mult-Agent Coordnaton: A Comparson Lemma Based Approach Yongcan Cao and We Ren Abstract In ths

More information

Randić Energy and Randić Estrada Index of a Graph

Randić Energy and Randić Estrada Index of a Graph EUROPEAN JOURNAL OF PURE AND APPLIED MATHEMATICS Vol. 5, No., 202, 88-96 ISSN 307-5543 www.ejpam.com SPECIAL ISSUE FOR THE INTERNATIONAL CONFERENCE ON APPLIED ANALYSIS AND ALGEBRA 29 JUNE -02JULY 20, ISTANBUL

More information

Structure and Drive Paul A. Jensen Copyright July 20, 2003

Structure and Drive Paul A. Jensen Copyright July 20, 2003 Structure and Drve Paul A. Jensen Copyrght July 20, 2003 A system s made up of several operatons wth flow passng between them. The structure of the system descrbes the flow paths from nputs to outputs.

More information

Assortment Optimization under MNL

Assortment Optimization under MNL Assortment Optmzaton under MNL Haotan Song Aprl 30, 2017 1 Introducton The assortment optmzaton problem ams to fnd the revenue-maxmzng assortment of products to offer when the prces of products are fxed.

More information

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family IOSR Journal of Mathematcs IOSR-JM) ISSN: 2278-5728. Volume 3, Issue 3 Sep-Oct. 202), PP 44-48 www.osrjournals.org Usng T.O.M to Estmate Parameter of dstrbutons that have not Sngle Exponental Famly Jubran

More information

Canonical transformations

Canonical transformations Canoncal transformatons November 23, 2014 Recall that we have defned a symplectc transformaton to be any lnear transformaton M A B leavng the symplectc form nvarant, Ω AB M A CM B DΩ CD Coordnate transformatons,

More information

PARTICIPATION FACTOR IN MODAL ANALYSIS OF POWER SYSTEMS STABILITY

PARTICIPATION FACTOR IN MODAL ANALYSIS OF POWER SYSTEMS STABILITY POZNAN UNIVE RSITY OF TE CHNOLOGY ACADE MIC JOURNALS No 86 Electrcal Engneerng 6 Volodymyr KONOVAL* Roman PRYTULA** PARTICIPATION FACTOR IN MODAL ANALYSIS OF POWER SYSTEMS STABILITY Ths paper provdes a

More information

Research Article Relative Smooth Topological Spaces

Research Article Relative Smooth Topological Spaces Advances n Fuzzy Systems Volume 2009, Artcle ID 172917, 5 pages do:10.1155/2009/172917 Research Artcle Relatve Smooth Topologcal Spaces B. Ghazanfar Department of Mathematcs, Faculty of Scence, Lorestan

More information

Lecture 17 : Stochastic Processes II

Lecture 17 : Stochastic Processes II : Stochastc Processes II 1 Contnuous-tme stochastc process So far we have studed dscrete-tme stochastc processes. We studed the concept of Makov chans and martngales, tme seres analyss, and regresson analyss

More information

A note on almost sure behavior of randomly weighted sums of φ-mixing random variables with φ-mixing weights

A note on almost sure behavior of randomly weighted sums of φ-mixing random variables with φ-mixing weights ACTA ET COMMENTATIONES UNIVERSITATIS TARTUENSIS DE MATHEMATICA Volume 7, Number 2, December 203 Avalable onlne at http://acutm.math.ut.ee A note on almost sure behavor of randomly weghted sums of φ-mxng

More information

Lecture 10 Support Vector Machines II

Lecture 10 Support Vector Machines II Lecture 10 Support Vector Machnes II 22 February 2016 Taylor B. Arnold Yale Statstcs STAT 365/665 1/28 Notes: Problem 3 s posted and due ths upcomng Frday There was an early bug n the fake-test data; fxed

More information

Global Sensitivity. Tuesday 20 th February, 2018

Global Sensitivity. Tuesday 20 th February, 2018 Global Senstvty Tuesday 2 th February, 28 ) Local Senstvty Most senstvty analyses [] are based on local estmates of senstvty, typcally by expandng the response n a Taylor seres about some specfc values

More information

Additional Codes using Finite Difference Method. 1 HJB Equation for Consumption-Saving Problem Without Uncertainty

Additional Codes using Finite Difference Method. 1 HJB Equation for Consumption-Saving Problem Without Uncertainty Addtonal Codes usng Fnte Dfference Method Benamn Moll 1 HJB Equaton for Consumpton-Savng Problem Wthout Uncertanty Before consderng the case wth stochastc ncome n http://www.prnceton.edu/~moll/ HACTproect/HACT_Numercal_Appendx.pdf,

More information

EXPANSIVE MAPPINGS. by W. R. Utz

EXPANSIVE MAPPINGS. by W. R. Utz Volume 3, 978 Pages 6 http://topology.auburn.edu/tp/ EXPANSIVE MAPPINGS by W. R. Utz Topology Proceedngs Web: http://topology.auburn.edu/tp/ Mal: Topology Proceedngs Department of Mathematcs & Statstcs

More information

Applied Nuclear Physics (Fall 2004) Lecture 23 (12/3/04) Nuclear Reactions: Energetics and Compound Nucleus

Applied Nuclear Physics (Fall 2004) Lecture 23 (12/3/04) Nuclear Reactions: Energetics and Compound Nucleus .101 Appled Nuclear Physcs (Fall 004) Lecture 3 (1/3/04) Nuclear Reactons: Energetcs and Compound Nucleus References: W. E. Meyerhof, Elements of Nuclear Physcs (McGraw-Hll, New York, 1967), Chap 5. Among

More information

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body Secton.. Moton.. The Materal Body and Moton hyscal materals n the real world are modeled usng an abstract mathematcal entty called a body. Ths body conssts of an nfnte number of materal partcles. Shown

More information

The Study of Teaching-learning-based Optimization Algorithm

The Study of Teaching-learning-based Optimization Algorithm Advanced Scence and Technology Letters Vol. (AST 06), pp.05- http://dx.do.org/0.57/astl.06. The Study of Teachng-learnng-based Optmzaton Algorthm u Sun, Yan fu, Lele Kong, Haolang Q,, Helongang Insttute

More information

Linear, affine, and convex sets and hulls In the sequel, unless otherwise specified, X will denote a real vector space.

Linear, affine, and convex sets and hulls In the sequel, unless otherwise specified, X will denote a real vector space. Lnear, affne, and convex sets and hulls In the sequel, unless otherwse specfed, X wll denote a real vector space. Lnes and segments. Gven two ponts x, y X, we defne xy = {x + t(y x) : t R} = {(1 t)x +

More information

Welfare Properties of General Equilibrium. What can be said about optimality properties of resource allocation implied by general equilibrium?

Welfare Properties of General Equilibrium. What can be said about optimality properties of resource allocation implied by general equilibrium? APPLIED WELFARE ECONOMICS AND POLICY ANALYSIS Welfare Propertes of General Equlbrum What can be sad about optmalty propertes of resource allocaton mpled by general equlbrum? Any crteron used to compare

More information

Refined Coding Bounds for Network Error Correction

Refined Coding Bounds for Network Error Correction Refned Codng Bounds for Network Error Correcton Shenghao Yang Department of Informaton Engneerng The Chnese Unversty of Hong Kong Shatn, N.T., Hong Kong shyang5@e.cuhk.edu.hk Raymond W. Yeung Department

More information

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U)

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U) Econ 413 Exam 13 H ANSWERS Settet er nndelt 9 deloppgaver, A,B,C, som alle anbefales å telle lkt for å gøre det ltt lettere å stå. Svar er gtt . Unfortunately, there s a prntng error n the hnt of

More information

Asynchronous Periodic Event-Triggered Coordination of Multi-Agent Systems

Asynchronous Periodic Event-Triggered Coordination of Multi-Agent Systems 017 IEEE 56th Annual Conference on Decson and Control (CDC) December 1-15, 017, Melbourne, Australa Asynchronous Perodc Event-Trggered Coordnaton of Mult-Agent Systems Yaohua Lu Cameron Nowzar Zh Tan Qng

More information

Lecture 4: Constant Time SVD Approximation

Lecture 4: Constant Time SVD Approximation Spectral Algorthms and Representatons eb. 17, Mar. 3 and 8, 005 Lecture 4: Constant Tme SVD Approxmaton Lecturer: Santosh Vempala Scrbe: Jangzhuo Chen Ths topc conssts of three lectures 0/17, 03/03, 03/08),

More information

Edge Isoperimetric Inequalities

Edge Isoperimetric Inequalities November 7, 2005 Ross M. Rchardson Edge Isopermetrc Inequaltes 1 Four Questons Recall that n the last lecture we looked at the problem of sopermetrc nequaltes n the hypercube, Q n. Our noton of boundary

More information