Degree Fluctuations and the Convergence Time of Consensus Algorithms

Size: px
Start display at page:

Download "Degree Fluctuations and the Convergence Time of Consensus Algorithms"

Transcription

1 Degree Fluctuatons and the Convergence Tme of Consensus Algorthms Alex Olshevsy John N. Tstsls Abstract We consder a consensus algorthm n whch every node n a tme-varyng undrected connected graph assgns equal weght to each of ts neghbors. Under the assumpton that the degree of any gven node s constant n tme, we show that the algorthm acheves consensus wthn a gven accuracy ɛ on n nodes n tme O(n 3 ln(n/ɛ)). Because there s a drect relaton between consensus algorthms n tme-varyng envronments and nhomogeneous random wals, our result also translates nto a general statement on such random wals. Moreover, we gve smple proofs that the convergence tme becomes exponentally large n the number of nodes n under slght relaxatons of the above assumptons. We prove that exponental convergence tme s possble for consensus algorthms on fxed drected graphs, and we use an example of Cao, Spelman, and Morse to gve a smple argument that the same s possble f the constant degrees assumpton s even slghtly relaxed. I. INTRODUCTION Consensus algorthms are a class of teratve update schemes that are commonly used as buldng blocs for the desgn of dstrbuted control laws. Ther man advantage s robustness n the presence of tme varyng envronments and and unexpected communcaton ln falures. Consensus algorthms have attracted sgnfcant nterest n a varety of contexts such as dstrbuted optmzaton [8],[7] or coverage control [2], and many other contexts nvolvng networs n whch central control s absent and communcaton capabltes are tmevaryng. Whle the convergence propertes of consensus algorthms n tme-varyng envronments are well understood, much less s nown about the correspondng convergence tmes. An nspecton of the classcal convergence proofs ([3], [3]) leads to convergence tme upper bounds that grow exponentally wth the number of nodes. It s then natural to loo for condtons under whch the convergence tme only grows polynomally, and ths s the subject of ths paper. In our man result, we show that a consensus algorthm n whch every node assgns equal weght to each of ts neghbors n an undrected, connected graph (where the graph can be tme-varyng) has polynomal convergence tme f the degree of any gven node s constant n tme. Because there s a drect relaton between consensus algorthms n tme-varyng envronments and nonhomogeneous random wals, our result also translates nto a general statement on such random wals. Research partally supported by the NSF under grant CMMI Alex Olshevsy s wth the Department of Mechancal and Aerospace Engneerng, Prnceton Unversty, Emal: aolshevs@prnceton.edu John N. Tstsls s wth the Laboratory for Informaton and Decson Systems, Department of Electrcal Engneerng and Computer Scence, Massachusetts Insttute of Technology. Emal: jnt@mt.edu. A. Model, notaton, and bacground In ths subsecton, we defne our notaton, the model of nterest, and some bacground on consensus algorthms. Gven a drected graph G, we wll use N (G) to denote the set {j (, j) s an edge} of drect successors of node n G, and d (G) to denote the cardnalty of N (G). Gven a sequence of drected graphs G(0), G(),..., G( ), we wll use the smpler notaton N (t), d (t) n place of N (G(t)), d (G(t)), and we wll mae a smlar smplfcaton for other varables of nterest. We are nterested n analyzng a consensus algorthm n whch a node assgns equal weght to each one of ts neghbors. We consder n nodes and assume that at each dscrete tme t, node stores a real number x (t). We let x(t) = (x (t),..., x n (t)). For any gven sequence of drected graphs G(0), G(), G(2),..., and any ntal vector x(0), the algorthm s descrbed by the update equaton x (t + ) = d (t) j N (t) whch can also be wrtten n the form x(t + ) = A(t)x(t), x j (t), =,..., n, () for a sutably defned sequence of matrces A(0), A(),..., A(t ). The graphs G(t), whch appear n the above update rule through d (t) and N (t), correspond to nformaton flow among the agents; the edge (, j) s present n G(t) f and only f agent uses the value x j (t) of agent j n ts update at tme t. To reflect the fact that every agent always has access to ts own nformaton, we assume that every graph G(t) contans all the self-loops (, ). Note that we have [A(t)] j > 0 f and only f (, j) s an edge n G(t). It s well nown ([8], [3]) that, subject to some natural condtons on the graph sequence, every component of x(t) converges to a common value. In ths paper, we focus on the convergence rate of ths process n some natural settngs. To quantfy the progress of the algorthm towards consensus, we wll use the functon S(x) = max x mn x. For any ɛ > 0, we wll say that a sequence of graphs G(0), G(),..., G( ) (alternatvely, a sequence of matrces A(0), A(),..., A( )) results n ɛ-consensus f S(x()) ɛs(x(0)) for all ntal vectors x(0). We wll focus on graph sequences n whch every graph G(t) s bdrectonal, meanng that f (, j) s an edge n

2 G(t), then so s (j, ). In practce, graphs that capture nformaton flows are often bdrectonal. For example, G(t) s bdrectonal f: () G(t) contans all the edges between agents that are physcaly wthn some dstance of each other; () G(t) contans all the edges between agents that have lneof-sght vews of each other; () G(t) contans the edges correspondng to pars of agents that can send messages to each other usng a protocol that reles on acnowledgements. It s an mmedate consequence of exstng convergence proofs ([3], [3]) that any sequence of Cn n ln(/ɛ) connected bdrectonal graphs, wth self-loops at every node, results n ɛ-consensus. Here, C s a constant that does not depend on the problem parameters n and ɛ. We are nterested n smple condtons under whch the undesrable O(n n ) scalng becomes polynomal n n. B. Our results Our contrbutons are twofold. Frst, n Secton II, we prove our man result. Theorem. Consder a sequence G(0), G(),..., G( ) of connected bdrectonal graphs, wth self-loops at each node. Suppose, furthermore, that the degree of each node stays constant n tme,.e., d (t) = d (t ),, t, t. If the length of the graph sequence s at least n 3 ln 2n2 ɛ, then ɛ-consensus s acheved. To put Theorem n perspectve, we note that polynomal convergence tmes were only nown for the cases where: (a) the graphs G(t) are the same at each tme t, bdrectonal, connected, wth all self-loops present [4]; or (b) n addton to some natural connectvty assumptons, the underlyng teraton matrces are doubly stochastc, whch, for the consensus algorthms consdered here amounts to an assumpton that each graph G(t) s regular [6]. Theorem can be vewed as a generalzaton of the above two results. In Secton III, we gve an nterpretaton of our results n terms of Marov chans. Theorem can be nterpreted as provdng a suffcent condton for a random wal on a tmevaryng graph to forget ts ntal dstrbuton n polynomal tme. In Secton IV, we show through examples that relaxng the assumptons of Theorem even slghtly can lead to a convergence tme whch s exponental n n. Specfcally, we show the followng. () If we do not requre the graphs G(t) to be bdrectonal, exponental convergence tme s possble, even f the graphs G(t) do not change wth tme; () If we replace the assumpton that each d (t) s ndependent of t wth the weaer assumpton that the sorted degree sequence (say, n non-ncreasng order) s ndependent of t (thus allowng nodes to swap degrees), exponental convergence tme s possble. Whle ths fact was nown (although unpublshed) [5], our contrbuton s to provde a smple proof. The presentaton of our counterexamples captalzes on the correspondence wth Marov chans dscussed n Secton III. In summary: for connected bdrectonal graphs wth selfloops, unchangng degrees s a suffcent condton for polynomal tme convergence, but relaxng t even slghtly by ether allowng the nodes to swap degrees or by losng ln symmetry leads to the possblty of exponental convergence tme. C. Prevous wor There s consderable and growng lterature on the convergence tme of consensus algorthms. We only menton papers that are closest to our own wor, omttng references to the lterature on varous aspects of consensus convergence tmes that we do not address here, such as topology desgn, performance n geometrc random graphs, etc. Worst-case upper bounds on the convergence tmes of consensus algorthms have been establshed n [8], [6], [7], [], [2], [], [0], [9]. The papers [8], [6], [7] consdered a settng slghtly more general than ours, and establshed exponental upper bounds. The papers [], [2] addressed the convergence tmes of consensus algorthms n terms of spannng trees that capture the nformaton flow between the nodes. It was observed that n several cases ths approach produces tght estmates of the convergence tmes. The wors [], [0], [9] tae a geometrc approach, and consder the convergence tme n a somewhat dfferent model, nvolvng nteractons between geographc nearest neghbors. These papers fnd that the convergence tme s qute hgh (ether sngly exponental or terated exponental, dependng on the model). The orgnal papers [8], [3] also consdered the effect of delays on convergence; some more recent wor on ths subject may be found n [7] and [4]. Our wor dffers from these papers n that our convergence tme bounds are polynomal n n. To the best of our nowledge, polynomal bounds on the partcular consensus algorthm consdered n ths paper had been derved earler only n [4] and [6]. Our wor encompasses a much wder class of stuatons than [4], whch requred the graphs G(t) to be constant n tme. Moreover, our results may be vewed as complementary to those n [6], whch proved a polynomal convergence tme bound for averagng algorthms, nvolvng doubly stochastc matrces. II. PROOF OF THEOREM As n the statement of Theorem, we assume that we are gven a sequence of bdrectonal connected graphs G(0), G(),..., wth self-loops at each node, and such that d (t) s the same n each G(t). We wll thus drop the parameter t and refer to the degree of node smply as d. Observe that d > 0 due to the connectvty of the graphs G(t).

3 We wll use G to refer to the class of bdrectonal connected graphs wth self-loops at every node such that the degree of node s d. We let D be the n n dagonal matrx whose th dagonal entry s d. We wll use E (G) to denote the edge set of a graph G. We wll sometmes fnd t convenent to use the notaton E(G) to refer to the set of unordered pars (, j) such that the ordered pars (, j) and (j, ) belong to E (G). Defnton. We defne the nner product, d by n x, y d = d x y. = Note that because d > 0 for all, x, y d s a vald nner product. Defnton 2. Gven a drected graph G, we defne the update matrx A(G) by { /d (G), f j N (G), [A(G)] j = 0, otherwse. We use A(t) as a shorthand for A(G(t)), so that Eq. () can be wrtten as x(t + ) = A(t)x(t). (2) Conversely, gven an update matrx A of the above form, we wll use G(A) to denote the graph G whose update matrx s A. We use N (A) as shorthand for N (G(A)); the quanttes d (A), E(A), and E (A) are defned smlarly. Fnally, we use A to denote the set of update matrces A(G) assocated wth graphs G G. Note that DA s the adjacency matrx assocated wth the graph correspondng to A. Gven that we restrct to bdrectonal graphs, DA s symmetrc for every A A. Lemma. For any A A, we have x, Ay d = (,j) E (A) x y j. Proof: We have x, Ay d = x T DAy = n = j N (A) x y j = (,j) E (A) x y j Lemma 2. Each A A s self-adjont wth respect to the nner product, d. Proof: Usng the fact that DA s symmetrc and D s dagonal, we have x, Ay d = x T DAy = x T (DA) T y = x T A T Dy = Ax, y d. Lemma 2 s the reason for ntroducng the nner product, d. The fact that matrces n the set A are self-adjont plays a central role n the analyss of the algorthm (2). One of ts consequences s that matrces n A have real egenvalues. We use the notaton λ (A) to denote the th largest egenvalue of a matrx A A. Note that every A A s a stochastc matrx, and therefore λ (A) =. Next, we dentfy a weghted average that s preserved by the teraton x(t + ) = A(t)x(t). For any x, we let x = x, d, d. We observe that for any A A, x, d = x T D = x T DA = x T A T D = Ax, d, where the second equalty used the fact that A s a stochastc matrx. Therefore, Ax = x, A A. Wth these prelmnares n place, we now proceed to the man part of our analyss, whch s based on the Lyapunov functon n V (x) = x x, x x d = d (x x) 2. = The next lemma quantfes the decrease of V ( ) when a vector x s multpled by some matrx A A. Lemma 3. For any A A and any vector x, we have V (Ax) λ 2 (A 2 )V (x). Proof: Fx some A A and some vector x. Snce Ax = x, the lemma asserts that Ax x, Ax x d λ 2 (A 2 ) x x, x x d. Let y = x x, and note that y, d = 0. It therefore suffces to show that Ay, Ay d λ 2 (A 2 ) y, y d, for all y wth y, d = 0; or, equvalently, that y, y d Ay, Ay d y, y d λ 2 2(A), f y, d = 0 and y 0. To establsh Eq. (3), we note that y, y d Ay, Ay d y, y d y, A 2 y d mn = mn y, d =0 y, y d y, d =0 y, y d y, (I A 2 )y d = mn, y, d =0 y, y d where n the mnma above we only consder nonzero vectors y; note that these mnma are attaned t suffces to consder vectors y on the unt ball, a compact set. Now, observe that I A 2 s also self-adjont under the nner product, d. Moreover, snce A (and, therefore, A 2 as well) s stochastc, the smallest egenvalue of I A 2 s 0, wth an assocated egenvector of. Consequently, by the Courant-Fscher varatonal characterzaton of egenvalues, the expresson on the rght s the second smallest egenvalue of I A 2 : y, y d Ay, Ay d mn = λ n (I A 2 ) = λ 2 (A 2 ), y, d =0 y, y d (3)

4 whch concludes the proof. Thus, to bound how much V (x) decreases at each step, t suffces to obtan an upper bound on λ 2 (A 2 ), for matrces A A. Ths can be done usng the next lemma, whch s the the man result of [4]. We nclude a short proof for completeness. Lemma 4. Let A A, and let l be the dameter of the graph G(A). Then, λ2 (A 2 )= max{ λ n (A), λ 2 (A)} where d max s the largest of the degrees d. nd max l, Proof: The frst nequalty follows because the egenvalues of A 2 are the squares of the egenvalues of A. For the second nequalty, usng agan the Courant-Fsher characterzaton, and some easy algebra, we have λ 2 (A) = max x, d =0 x,x d = = mn x, d =0 x,x d = x, Ax d = (,j) E(A) max x, d =0 x,x d = (,j) E (A) (x x j ) 2. Thus, t suffces to show that (x x j ) 2 mn x, d =0 x,x d = (,j) E(A) nd max l x x j Towards ths purpose, we carry out a varaton of an argument frst used n [4]. Fx some x that satsfes x, d = 0 and x, x d =. Wthout loss of generalty, we assume that () node has the largest value of d x 2, () node has the smallest value of x, and () the shortest path from node to s (, 2), (2, 3),..., (, ). The condton x, x d = mples that d x 2 /n, and consequently that x / nd max ; the requrement x, d = 0 mples that x < 0. Thus, x x / nd max, whch we wrte as (x x 2 ) + (x 2 x 3 ) + + (x x ) Applyng the Cauchy-Schwarz nequalty, we get = (x x + ) 2 nd max. ndmax. We then use the fact that l, to obtan the clamed bound on λ 2 (A). As for λ n, we observe that the dagonal entres of A are at least /n and the row sums are. The Gershgorn crcle theorem mmedately gves λ n + /n, whch s stronger than the bound we have clamed. We can now complete the proof of Theorem. Lemma 3 descrbes the decrease n the varance V (x(t)) n terms of λ 2 (A 2 (t)), and Lemma 4 gves us a way to upper bound the latter quantty. Proof of Theorem : Usng Lemmas 3 and 4, and the bounds d max n, l n, we see that for every A A, and every x, we have V (Ax) λ 2 (A 2 )V (x) = max{ λ n (A) 2, λ 2 2(A)}V (x) ( ) 2V (x). n 3 Because the defnton of ɛ-consensus s n terms of S(x) rather than V (x), we need to relate these two quanttes. On the one hand, for every x, we have n n V (x) = d (x x) 2 n (x x) 2 n 2 S 2 (x). = = On the other hand, for every x, we have V (x) max (x x) 2 4 (max x mn x ) 2 = 4 S2 (x). Suppose that t n 3 ln(2n 2 /ɛ). Then, S(x(t)) ( 4V (x(t)) 2 ) n 3 ln(2n 2 /ɛ) V (x(0)) n 3 2n 2( n 3 ) n 3 ln(2n 2 /ɛ) S(x(0)) 2n 2 e ln(2n2 /ɛ) S(x(0)) = ɛs(x(0)). (We have used here the nequalty ( c/n) n e c, for c > 0.) III. MARKOV CHAIN INTERPRETATION In ths secton, we gve an alternatve nterpretaton of the convergence tme of a consensus algorthm n terms of nhomogeneous Marov chans. In the next secton, we wll use ths nterpretaton to gve some examples of graph sequences that do not satsfy Theorem and whch have exponentally large convergence tmes. We consder an nhomogeneous Marov chan whose transton probablty matrx at tme s A(). We fx some tme t and defne P = A(0)A() A(t ). Ths s the assocated t-step transton probablty matrx: the j-th entry of P, denoted by p j, s the probablty that the state at tme t s j, gven that the ntal state s. Let p be the vector whose th component s p ; thus p T s the th row of P. We address a queston whch s generc n the study of Marov chans, namely, whether the chan eventually forgets ts ntal state,.e., whether for all, j, p p j converges to zero as t ncreases, and f so, at what rate. We wll say that the sequence of matrces A(0), A(),..., A(t ) s ɛ-forgetful f for all, j, we have 2 p p j ɛ. The above quantty, 2 max,j p p j s nown as the coeffcent of ergodcty of the matrx P, and appears often n the study of consensus algorthms (see, for example, [8]). The matrx P can also be nterpreted n terms of consensus updates: an ntal vector x s multpled by the matrces A(t

5 ), A(t 2),..., A(0), to produce the vector P x; note that the dfferent matrces are now appled n the reverse order. The result that follows relates the tmes to acheve ɛ-consensus or ɛ-forgetfulness, and s essentally the same as Proposton 4.5 of [5]. IV. SOME COUNTEREXAMPLES A. The bdrectonalty requrement n Theorem cannot be relaxed Proposton 2. The sequence of matrces A(0), A(),..., A(t ) s ɛ-forgetful f and only f the sequence of matrces A(t ), A(t 2),..., A(0) results n ɛ-consensus (.e., S(P x) ɛs(x), for every vector x.) Proof: Suppose that the matrx sequence A(0), A(),..., A(t ) s ɛ-forgetful,.e., that 2 p p j ɛ, for all and j. Gven a vector x, let c = (max x + mn x )/2. Note that x c = (max x mn x )/2 = S(x)/2. We then have [P x] [P x] j = (p p j )x = (p p j )(x c) p p j x c 2ɛS(x)/2 = ɛs(x). Snce ths s true for every and j, we obtan S(P x) ɛs(x), and the sequence A(t ), A(t 2),..., A(0) results n ɛ- consensus. Conversely, suppose that the sequence of matrces A(t ), A(t 2),..., A(0) results n ɛ-consensus. Fx some and j. Let x be a vector whose th component s /2 f p p j and /2 otherwse. Note that S(x) =. We have 2 p p j = (p T p T j )x = [P x] [P x] j ɛs(x) = ɛ, where the last nequalty made use of the ɛ-consensus assumpton. Thus, the sequence of matrces A(0), A(),..., A(t ) s ɛ-forgetful. We wll use Proposton 2 for the specal case of Marov chans that are random wals. Gven a drected graph G(t), we let, as before, N (t) be the set of nodes j for whch the edge (, j) s present. In the random wal assocated wth ths graph, f the state at tme t s, the state at tme t+ s chosen to be one of the elements of N (t), wth equal probablty. We let A(t) be the assocated transton probablty matrx. We wll say that a sequence of graphs s ɛ-forgetful whenever the correspondng sequence of transton probablty matrces s ɛ-forgetful. Proposton 2 allows us to renterpret Theorem as follows: random wals on tme-varyng bdrectonal connected graphs wth self-loops forget ther ntal dstrbuton n a polynomal number of steps. Proposton 2 also has a corollary whch we wll use later. It s based on the followng observaton: concatenatng two sequence of graphs, each of whch acheves ɛ-consensus, results n a sequence whch acheves ɛ 2 -consensus. Corollary 3. Suppose that a sequence of graphs s ɛ-forgetful. Then, concatenatng ths sequence wth tself tmes results n a sequence whch s ɛ -forgetful. Fg.. The graph used n Proposton 4. In ths subsecton, we show that t s mpossble to drop the assumpton that the graph s bdrectonal, even when the graph does not change wth tme. Proposton 4. Let G be the graph shown n Fgure. Consder the graph sequence consstng of G, repeated t tmes. For ths graph sequence to result n (/8)-consensus, we must have t 2 n/2 /6. Proof: Suppose that ths graph sequence of length t results n ɛ-consensus. By Proposton 2, t s (/8)-forgetful. In partcular, p p n /4 (recall that p j stands for the t-step transton probablty from to j). For n/2, let T be the frst tme that a random wal that starts at state vsts the bottom part of the graph, and let δ be the probablty that T s less than or equal to t. Note that condtonal on startng at and never transtonng to the bottom half, the probablty of beng at state at tme t s /2; thus, p 2 ( δ ). Furthermore, p n δ by symmetry. Therefore, 4 p p n 2 ( δ ) δ = δ, whch yelds δ /6. Usng a straghtforward couplng argument, we have δ δ /6, for = 2,..., n/2. By vewng perods of length t as a sngle attempt to get to the bottom half of the graph, wth each attempt havng probablty at least /6 to succeed, we conclude that E[T ] 6t. On the other hand, T s just the frst tme untl the wal moves to the rght n/2 consecutve tmes. Each tme that the random wal returns to state, we have a new tral wth probablty of success equal to 2 n/2. Thus, 2 n/2 E[T ] 6t, whch yelds the desred result. Remar: It s not hard to see that f we add a self-loop to each node n Fgure, the result holds wth mnor modfcatons.

6 E[T ] 4t, whch yelds the desred result. Fg. 2. The top-left fgure shows graph G(0); top-rght shows G(); bottom-left shows G((n/2) 2); bottom-rght shows G((n/2) ). As these fgures llustrate, G(t + ) s obtaned by applyng a crcular shft to each half of G(t). All edges are bdrectonal, and every node has a self-loop whch s not shown. For aesthetc reasons, nstead of labelng the nodes as,..., n, we label them wth,..., n/2 and,..., (n/2). B. The unchangng degrees condton n Theorem cannot be relaxed In ths subsecton, we show that t s mpossble to relax the condton of unchangng degrees n Theorem. In partcular, f we only mpose the slghtly weaer condton that the sorted degree sequence (the non-ncreasng lst of node degrees) does not change wth tme, the tme to acheve ɛ-consensus can grow exponentally wth n. Ths s an unpublshed result of Cao, Spelman, and Morse [5]; we provde here a smple proof. Proposton 5. Let n be even and let t be an nteger multple of n/2. Consder the graph sequence of length t = n/2, consstng of perodc repettons of the reversal of the lengthn/2 sequence descrbed n Fgure 2. For ths graph sequence to result n (/4)-consensus, we must have t 2 (n/2) /8. Proof: Suppose that ths graph sequence of length t results n (/4)-consensus. Then Proposton 2 mples that the sequence of length n/2 consstng of perodc repettons of the length n/2 sequence descrbed n Fgure 2 s (/4)- forgetful. Let p j be the assocated t-step transton probabltes. Let T be the tme that t taes for a random wal that starts at state n/2 to cross nto the rght-hand sde part of the graph. Let δ be the probablty that T s less than or equal to t. Let R = {... (n/2) }. We have j R p (n/2),j δ and, usng symmetry, j R p (n/2),j δ. Usng the fact that the graph sequence s (/4)-forgetful n the frst nequalty below, we have 2 j R p (n/2),j p ((n/2),j p (n/2),j j R j R ( δ) δ = 2δ, p (n/2),j whch yelds δ /4. Argung as n the proof of Proposton 4 (node n/2 s the least favorable non-central startng state n the left-hade sde of the graph), we obtan E[T ] 4t. Note that for a wal that starts at state n/2 to cross nto the rght-hand sde part of the graph, t must frst tae a selfloop (n/2) consecutve tmes. Consequently, 2 (n/2) V. CONCLUSIONS The man contrbuton of ths paper s Theorem, whch shows that consensus algorthms converge n polynomal tme for a large class of graph sequences. We also gave smple proofs showng that ths fndng s fragle, and even slght relaxatons of the hypotheses cause the concluson to fal. A smlar result s avalable for the case of doubly stochastc update matrces A(t), and our result can be vewed as complementary [6]. Interestngly, both results rely on a sutable quadratc Lyapunov functon as well as on the fact that all update matrces share a common left egenvector. REFERENCES [] D. Angel, P.-A. Blman, Tght estmates for convergence of some nonstatonary consensus algorthms, Systems and Control Letters, vol. 57, no. 2, pp , [2] D. Angel, P.-A. Blman, Convergence speed of unsteady dstrbuted consensus: decay estmate along the settlng spannng-trees, SIAM Journal on Control and Optmzaton, vol. 48, no., pp. -32, [3] D. P. Bertseas and J. N. Tstsls, Parallel and Dstrbuted Computaton: Numercal Methods, Prentce Hall, 989. [4] P.-A. Blman, A. Nedc, A. Ozdaglar, Rate of convergence for consensus wth delays, Proceedngs of the 47th IEEE Conference on Decson and Control, Cancun, Mexco, Dec [5] Mng Cao, personal communcaton. [6] M. Cao, A. S. Morse, B. D. O. Anderson, Reachng a consensus n a dynamcally changng envronment: a graphcal approach, SIAM Journal on Control and Optmzaton, vol. 47, no. 2, pp , [7] M. Cao, A. S. Morse, B. D. O. Anderson, Reachng a consensus n a dynamcally changng envronment: convergence rates, measurement delays, and asynchronous events, SIAM Journal on Control and Optmzaton, vol. 47, no. 2, pp , [8] M. Cao, D. Spelman, A. S. Morse, A lower bound on convergence of a dstrbuted networ consensus slgorthm, Proceedngs of the 44th IEEE Conference on Decson and Control and European Control Conference, Madrd, Span, Dec [9] B. Chazelle, Analytcal tools for natural algorthms, Proceedngs of the Frst Symposum on Innovatons n Computer Scence, Bejng, Chna, Jan [0] B. Chazelle, The convergence of brd flocng, Proceedngs of the 26th Annual Symposum on Computatonal Geometry, Snowbrd, USA, Jun [] B. Chazelle, Natural algorthms, Proceedngs of the ACM-SIAM Symposum on Dscrete Algorthms, New Yor, USA, Jan [2] C. Gao, J. Cortes, F. Bullo, Notes on averagng over acyclc graphs and dscrete coverage control, Automatca, vol. 44, no. 8, pp , [3] A. Jadbabae, J. Ln, and A. S. Morse, Coordnaton of groups of moble autonomous agents usng nearest neghbor rules, IEEE Transactons on Automatc Control, vol. 48, no. 3, pp , [4] H. J. Landau, A. M. Odlyzo, Bounds for egenvalues of certan stochastc matrces, Lnear Algebra and Its Applcatons, vol. 38, pp. 5-5, 98. [5] D. A. Levn, Y. Peres, and E. L. Wlmer, Marov Chans and Mxng Tmes, Amercan Mathematcal Socety, [6] A. Nedc, A. Olshevsy, A. Ozdaglar, and J. N. Tstsls, On dstrbuted averagng algorthms and quantzaton effects, IEEE Transactons on Automatc Control, vol. 54, no., pp , [7] A. Nedc, A. Ozdaglar, Dstrbuted subgradent methods for mult-agent optmzaton, IEEE Transactons on Automatc Control, vol. 54, no., pp. 48-6, [8] J. N. Tstsls, D. P. Bertseas, and M. Athans, Dstrbuted asynchronous determnstc and stochastc gradent optmzaton algorthms, IEEE Transactons on Automatc Control, vol. 3, no. 9, 986, pp

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

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

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

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

Lecture Space-Bounded Derandomization

Lecture Space-Bounded Derandomization Notes on Complexty Theory Last updated: October, 2008 Jonathan Katz Lecture Space-Bounded Derandomzaton 1 Space-Bounded Derandomzaton We now dscuss derandomzaton of space-bounded algorthms. Here non-trval

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

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

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

2.3 Nilpotent endomorphisms

2.3 Nilpotent endomorphisms s a block dagonal matrx, wth A Mat dm U (C) In fact, we can assume that B = B 1 B k, wth B an ordered bass of U, and that A = [f U ] B, where f U : U U s the restrcton of f to U 40 23 Nlpotent endomorphsms

More information

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 12 10/21/2013. Martingale Concentration Inequalities and Applications

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 12 10/21/2013. Martingale Concentration Inequalities and Applications MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.65/15.070J Fall 013 Lecture 1 10/1/013 Martngale Concentraton Inequaltes and Applcatons Content. 1. Exponental concentraton for martngales wth bounded ncrements.

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

Lecture 10: May 6, 2013

Lecture 10: May 6, 2013 TTIC/CMSC 31150 Mathematcal Toolkt Sprng 013 Madhur Tulsan Lecture 10: May 6, 013 Scrbe: Wenje Luo In today s lecture, we manly talked about random walk on graphs and ntroduce the concept of graph expander,

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

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

Lecture 3. Ax x i a i. i i

Lecture 3. Ax x i a i. i i 18.409 The Behavor of Algorthms n Practce 2/14/2 Lecturer: Dan Spelman Lecture 3 Scrbe: Arvnd Sankar 1 Largest sngular value In order to bound the condton number, we need an upper bound on the largest

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

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

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

BOUNDEDNESS OF THE RIESZ TRANSFORM WITH MATRIX A 2 WEIGHTS

BOUNDEDNESS OF THE RIESZ TRANSFORM WITH MATRIX A 2 WEIGHTS BOUNDEDNESS OF THE IESZ TANSFOM WITH MATIX A WEIGHTS Introducton Let L = L ( n, be the functon space wth norm (ˆ f L = f(x C dx d < For a d d matrx valued functon W : wth W (x postve sem-defnte for all

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

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

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

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

Stanford University CS359G: Graph Partitioning and Expanders Handout 4 Luca Trevisan January 13, 2011

Stanford University CS359G: Graph Partitioning and Expanders Handout 4 Luca Trevisan January 13, 2011 Stanford Unversty CS359G: Graph Parttonng and Expanders Handout 4 Luca Trevsan January 3, 0 Lecture 4 In whch we prove the dffcult drecton of Cheeger s nequalty. As n the past lectures, consder an undrected

More information

1 GSW Iterative Techniques for y = Ax

1 GSW Iterative Techniques for y = Ax 1 for y = A I m gong to cheat here. here are a lot of teratve technques that can be used to solve the general case of a set of smultaneous equatons (wrtten n the matr form as y = A), but ths chapter sn

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

C/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1

C/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1 C/CS/Phy9 Problem Set 3 Solutons Out: Oct, 8 Suppose you have two qubts n some arbtrary entangled state ψ You apply the teleportaton protocol to each of the qubts separately What s the resultng state obtaned

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

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

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

Eigenvalues of Random Graphs

Eigenvalues of Random Graphs Spectral Graph Theory Lecture 2 Egenvalues of Random Graphs Danel A. Spelman November 4, 202 2. Introducton In ths lecture, we consder a random graph on n vertces n whch each edge s chosen to be n the

More information

Open Systems: Chemical Potential and Partial Molar Quantities Chemical Potential

Open Systems: Chemical Potential and Partial Molar Quantities Chemical Potential Open Systems: Chemcal Potental and Partal Molar Quanttes Chemcal Potental For closed systems, we have derved the followng relatonshps: du = TdS pdv dh = TdS + Vdp da = SdT pdv dg = VdP SdT For open systems,

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

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

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

Polynomials. 1 More properties of polynomials

Polynomials. 1 More properties of polynomials Polynomals 1 More propertes of polynomals Recall that, for R a commutatve rng wth unty (as wth all rngs n ths course unless otherwse noted), we defne R[x] to be the set of expressons n =0 a x, where a

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

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

U.C. Berkeley CS278: Computational Complexity Professor Luca Trevisan 2/21/2008. Notes for Lecture 8

U.C. Berkeley CS278: Computational Complexity Professor Luca Trevisan 2/21/2008. Notes for Lecture 8 U.C. Berkeley CS278: Computatonal Complexty Handout N8 Professor Luca Trevsan 2/21/2008 Notes for Lecture 8 1 Undrected Connectvty In the undrected s t connectvty problem (abbrevated ST-UCONN) we are gven

More information

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity LINEAR REGRESSION ANALYSIS MODULE IX Lecture - 30 Multcollnearty Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur 2 Remedes for multcollnearty Varous technques have

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

Matrix Approximation via Sampling, Subspace Embedding. 1 Solving Linear Systems Using SVD

Matrix Approximation via Sampling, Subspace Embedding. 1 Solving Linear Systems Using SVD Matrx Approxmaton va Samplng, Subspace Embeddng Lecturer: Anup Rao Scrbe: Rashth Sharma, Peng Zhang 0/01/016 1 Solvng Lnear Systems Usng SVD Two applcatons of SVD have been covered so far. Today we loo

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

Solutions Homework 4 March 5, 2018

Solutions Homework 4 March 5, 2018 1 Solutons Homework 4 March 5, 018 Soluton to Exercse 5.1.8: Let a IR be a translaton and c > 0 be a re-scalng. ˆb1 (cx + a) cx n + a (cx 1 + a) c x n x 1 cˆb 1 (x), whch shows ˆb 1 s locaton nvarant and

More information

Prof. Dr. I. Nasser Phys 630, T Aug-15 One_dimensional_Ising_Model

Prof. Dr. I. Nasser Phys 630, T Aug-15 One_dimensional_Ising_Model EXACT OE-DIMESIOAL ISIG MODEL The one-dmensonal Isng model conssts of a chan of spns, each spn nteractng only wth ts two nearest neghbors. The smple Isng problem n one dmenson can be solved drectly n several

More information

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve

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

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

Min Cut, Fast Cut, Polynomial Identities

Min Cut, Fast Cut, Polynomial Identities Randomzed Algorthms, Summer 016 Mn Cut, Fast Cut, Polynomal Identtes Instructor: Thomas Kesselhem and Kurt Mehlhorn 1 Mn Cuts n Graphs Lecture (5 pages) Throughout ths secton, G = (V, E) s a mult-graph.

More information

Finding Dense Subgraphs in G(n, 1/2)

Finding Dense Subgraphs in G(n, 1/2) Fndng Dense Subgraphs n Gn, 1/ Atsh Das Sarma 1, Amt Deshpande, and Rav Kannan 1 Georga Insttute of Technology,atsh@cc.gatech.edu Mcrosoft Research-Bangalore,amtdesh,annan@mcrosoft.com Abstract. Fndng

More information

Section 8.3 Polar Form of Complex Numbers

Section 8.3 Polar Form of Complex Numbers 80 Chapter 8 Secton 8 Polar Form of Complex Numbers From prevous classes, you may have encountered magnary numbers the square roots of negatve numbers and, more generally, complex numbers whch are the

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

TAIL BOUNDS FOR SUMS OF GEOMETRIC AND EXPONENTIAL VARIABLES

TAIL BOUNDS FOR SUMS OF GEOMETRIC AND EXPONENTIAL VARIABLES TAIL BOUNDS FOR SUMS OF GEOMETRIC AND EXPONENTIAL VARIABLES SVANTE JANSON Abstract. We gve explct bounds for the tal probabltes for sums of ndependent geometrc or exponental varables, possbly wth dfferent

More information

Feature Selection: Part 1

Feature Selection: Part 1 CSE 546: Machne Learnng Lecture 5 Feature Selecton: Part 1 Instructor: Sham Kakade 1 Regresson n the hgh dmensonal settng How do we learn when the number of features d s greater than the sample sze n?

More information

Chapter 13: Multiple Regression

Chapter 13: Multiple Regression Chapter 13: Multple Regresson 13.1 Developng the multple-regresson Model The general model can be descrbed as: It smplfes for two ndependent varables: The sample ft parameter b 0, b 1, and b are used to

More information

COS 521: Advanced Algorithms Game Theory and Linear Programming

COS 521: Advanced Algorithms Game Theory and Linear Programming COS 521: Advanced Algorthms Game Theory and Lnear Programmng Moses Charkar February 27, 2013 In these notes, we ntroduce some basc concepts n game theory and lnear programmng (LP). We show a connecton

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

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

U.C. Berkeley CS294: Beyond Worst-Case Analysis Handout 6 Luca Trevisan September 12, 2017 U.C. Berkeley CS94: Beyond Worst-Case Analyss Handout 6 Luca Trevsan September, 07 Scrbed by Theo McKenze Lecture 6 In whch we study the spectrum of random graphs. Overvew When attemptng to fnd n polynomal

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

Finding Primitive Roots Pseudo-Deterministically

Finding Primitive Roots Pseudo-Deterministically Electronc Colloquum on Computatonal Complexty, Report No 207 (205) Fndng Prmtve Roots Pseudo-Determnstcally Ofer Grossman December 22, 205 Abstract Pseudo-determnstc algorthms are randomzed search algorthms

More information

Games of Threats. Elon Kohlberg Abraham Neyman. Working Paper

Games of Threats. Elon Kohlberg Abraham Neyman. Working Paper Games of Threats Elon Kohlberg Abraham Neyman Workng Paper 18-023 Games of Threats Elon Kohlberg Harvard Busness School Abraham Neyman The Hebrew Unversty of Jerusalem Workng Paper 18-023 Copyrght 2017

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

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

Perron Vectors of an Irreducible Nonnegative Interval Matrix

Perron Vectors of an Irreducible Nonnegative Interval Matrix Perron Vectors of an Irreducble Nonnegatve Interval Matrx Jr Rohn August 4 2005 Abstract As s well known an rreducble nonnegatve matrx possesses a unquely determned Perron vector. As the man result of

More information

The equation of motion of a dynamical system is given by a set of differential equations. That is (1)

The equation of motion of a dynamical system is given by a set of differential equations. That is (1) Dynamcal Systems Many engneerng and natural systems are dynamcal systems. For example a pendulum s a dynamcal system. State l The state of the dynamcal system specfes t condtons. For a pendulum n the absence

More information

The Minimum Universal Cost Flow in an Infeasible Flow Network

The Minimum Universal Cost Flow in an Infeasible Flow Network Journal of Scences, Islamc Republc of Iran 17(2): 175-180 (2006) Unversty of Tehran, ISSN 1016-1104 http://jscencesutacr The Mnmum Unversal Cost Flow n an Infeasble Flow Network H Saleh Fathabad * M Bagheran

More information

Hidden Markov Models & The Multivariate Gaussian (10/26/04)

Hidden Markov Models & The Multivariate Gaussian (10/26/04) CS281A/Stat241A: Statstcal Learnng Theory Hdden Markov Models & The Multvarate Gaussan (10/26/04) Lecturer: Mchael I. Jordan Scrbes: Jonathan W. Hu 1 Hdden Markov Models As a bref revew, hdden Markov models

More information

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS Avalable onlne at http://sck.org J. Math. Comput. Sc. 3 (3), No., 6-3 ISSN: 97-537 COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

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

Solution Thermodynamics

Solution Thermodynamics Soluton hermodynamcs usng Wagner Notaton by Stanley. Howard Department of aterals and etallurgcal Engneerng South Dakota School of nes and echnology Rapd Cty, SD 57701 January 7, 001 Soluton hermodynamcs

More information

Singular Value Decomposition: Theory and Applications

Singular Value Decomposition: Theory and Applications Sngular Value Decomposton: Theory and Applcatons Danel Khashab Sprng 2015 Last Update: March 2, 2015 1 Introducton A = UDV where columns of U and V are orthonormal and matrx D s dagonal wth postve real

More information

x = , so that calculated

x = , so that calculated Stat 4, secton Sngle Factor ANOVA notes by Tm Plachowsk n chapter 8 we conducted hypothess tests n whch we compared a sngle sample s mean or proporton to some hypotheszed value Chapter 9 expanded ths to

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

10-801: Advanced Optimization and Randomized Methods Lecture 2: Convex functions (Jan 15, 2014)

10-801: Advanced Optimization and Randomized Methods Lecture 2: Convex functions (Jan 15, 2014) 0-80: Advanced Optmzaton and Randomzed Methods Lecture : Convex functons (Jan 5, 04) Lecturer: Suvrt Sra Addr: Carnege Mellon Unversty, Sprng 04 Scrbes: Avnava Dubey, Ahmed Hefny Dsclamer: These notes

More information

2E Pattern Recognition Solutions to Introduction to Pattern Recognition, Chapter 2: Bayesian pattern classification

2E Pattern Recognition Solutions to Introduction to Pattern Recognition, Chapter 2: Bayesian pattern classification E395 - Pattern Recognton Solutons to Introducton to Pattern Recognton, Chapter : Bayesan pattern classfcaton Preface Ths document s a soluton manual for selected exercses from Introducton to Pattern Recognton

More information

A PROBABILITY-DRIVEN SEARCH ALGORITHM FOR SOLVING MULTI-OBJECTIVE OPTIMIZATION PROBLEMS

A PROBABILITY-DRIVEN SEARCH ALGORITHM FOR SOLVING MULTI-OBJECTIVE OPTIMIZATION PROBLEMS HCMC Unversty of Pedagogy Thong Nguyen Huu et al. A PROBABILITY-DRIVEN SEARCH ALGORITHM FOR SOLVING MULTI-OBJECTIVE OPTIMIZATION PROBLEMS Thong Nguyen Huu and Hao Tran Van Department of mathematcs-nformaton,

More information

Complete subgraphs in multipartite graphs

Complete subgraphs in multipartite graphs Complete subgraphs n multpartte graphs FLORIAN PFENDER Unverstät Rostock, Insttut für Mathematk D-18057 Rostock, Germany Floran.Pfender@un-rostock.de Abstract Turán s Theorem states that every graph G

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

CS : Algorithms and Uncertainty Lecture 17 Date: October 26, 2016

CS : Algorithms and Uncertainty Lecture 17 Date: October 26, 2016 CS 29-128: Algorthms and Uncertanty Lecture 17 Date: October 26, 2016 Instructor: Nkhl Bansal Scrbe: Mchael Denns 1 Introducton In ths lecture we wll be lookng nto the secretary problem, and an nterestng

More information

LECTURE 9 CANONICAL CORRELATION ANALYSIS

LECTURE 9 CANONICAL CORRELATION ANALYSIS LECURE 9 CANONICAL CORRELAION ANALYSIS Introducton he concept of canoncal correlaton arses when we want to quantfy the assocatons between two sets of varables. For example, suppose that the frst set of

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

Integrals and Invariants of Euler-Lagrange Equations

Integrals and Invariants of Euler-Lagrange Equations Lecture 16 Integrals and Invarants of Euler-Lagrange Equatons ME 256 at the Indan Insttute of Scence, Bengaluru Varatonal Methods and Structural Optmzaton G. K. Ananthasuresh Professor, Mechancal Engneerng,

More information

Maximum Likelihood Estimation of Binary Dependent Variables Models: Probit and Logit. 1. General Formulation of Binary Dependent Variables Models

Maximum Likelihood Estimation of Binary Dependent Variables Models: Probit and Logit. 1. General Formulation of Binary Dependent Variables Models ECO 452 -- OE 4: Probt and Logt Models ECO 452 -- OE 4 Maxmum Lkelhood Estmaton of Bnary Dependent Varables Models: Probt and Logt hs note demonstrates how to formulate bnary dependent varables models

More information

Formulas for the Determinant

Formulas for the Determinant page 224 224 CHAPTER 3 Determnants e t te t e 2t 38 A = e t 2te t e 2t e t te t 2e 2t 39 If 123 A = 345, 456 compute the matrx product A adj(a) What can you conclude about det(a)? For Problems 40 43, use

More information

The Expectation-Maximization Algorithm

The Expectation-Maximization Algorithm The Expectaton-Maxmaton Algorthm Charles Elan elan@cs.ucsd.edu November 16, 2007 Ths chapter explans the EM algorthm at multple levels of generalty. Secton 1 gves the standard hgh-level verson of the algorthm.

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

Chapter 8 Indicator Variables

Chapter 8 Indicator Variables Chapter 8 Indcator Varables In general, e explanatory varables n any regresson analyss are assumed to be quanttatve n nature. For example, e varables lke temperature, dstance, age etc. are quanttatve n

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

Inexact Newton Methods for Inverse Eigenvalue Problems

Inexact Newton Methods for Inverse Eigenvalue Problems Inexact Newton Methods for Inverse Egenvalue Problems Zheng-jan Ba Abstract In ths paper, we survey some of the latest development n usng nexact Newton-lke methods for solvng nverse egenvalue problems.

More information

ECE559VV Project Report

ECE559VV Project Report ECE559VV Project Report (Supplementary Notes Loc Xuan Bu I. MAX SUM-RATE SCHEDULING: THE UPLINK CASE We have seen (n the presentaton that, for downlnk (broadcast channels, the strategy maxmzng the sum-rate

More information

Estimation: Part 2. Chapter GREG estimation

Estimation: Part 2. Chapter GREG estimation Chapter 9 Estmaton: Part 2 9. GREG estmaton In Chapter 8, we have seen that the regresson estmator s an effcent estmator when there s a lnear relatonshp between y and x. In ths chapter, we generalzed the

More information

Lecture 4: Universal Hash Functions/Streaming Cont d

Lecture 4: Universal Hash Functions/Streaming Cont d CSE 5: Desgn and Analyss of Algorthms I Sprng 06 Lecture 4: Unversal Hash Functons/Streamng Cont d Lecturer: Shayan Oves Gharan Aprl 6th Scrbe: Jacob Schreber Dsclamer: These notes have not been subjected

More information

THE ARIMOTO-BLAHUT ALGORITHM FOR COMPUTATION OF CHANNEL CAPACITY. William A. Pearlman. References: S. Arimoto - IEEE Trans. Inform. Thy., Jan.

THE ARIMOTO-BLAHUT ALGORITHM FOR COMPUTATION OF CHANNEL CAPACITY. William A. Pearlman. References: S. Arimoto - IEEE Trans. Inform. Thy., Jan. THE ARIMOTO-BLAHUT ALGORITHM FOR COMPUTATION OF CHANNEL CAPACITY Wllam A. Pearlman 2002 References: S. Armoto - IEEE Trans. Inform. Thy., Jan. 1972 R. Blahut - IEEE Trans. Inform. Thy., July 1972 Recall

More information

Notes prepared by Prof Mrs) M.J. Gholba Class M.Sc Part(I) Information Technology

Notes prepared by Prof Mrs) M.J. Gholba Class M.Sc Part(I) Information Technology Inverse transformatons Generaton of random observatons from gven dstrbutons Assume that random numbers,,, are readly avalable, where each tself s a random varable whch s unformly dstrbuted over the range(,).

More information

Approximate Smallest Enclosing Balls

Approximate Smallest Enclosing Balls Chapter 5 Approxmate Smallest Enclosng Balls 5. Boundng Volumes A boundng volume for a set S R d s a superset of S wth a smple shape, for example a box, a ball, or an ellpsod. Fgure 5.: Boundng boxes Q(P

More information

Solutions to exam in SF1811 Optimization, Jan 14, 2015

Solutions to exam in SF1811 Optimization, Jan 14, 2015 Solutons to exam n SF8 Optmzaton, Jan 4, 25 3 3 O------O -4 \ / \ / The network: \/ where all lnks go from left to rght. /\ / \ / \ 6 O------O -5 2 4.(a) Let x = ( x 3, x 4, x 23, x 24 ) T, where the varable

More information

Report on Image warping

Report on Image warping Report on Image warpng Xuan Ne, Dec. 20, 2004 Ths document summarzed the algorthms of our mage warpng soluton for further study, and there s a detaled descrpton about the mplementaton of these algorthms.

More information

A new Approach for Solving Linear Ordinary Differential Equations

A new Approach for Solving Linear Ordinary Differential Equations , ISSN 974-57X (Onlne), ISSN 974-5718 (Prnt), Vol. ; Issue No. 1; Year 14, Copyrght 13-14 by CESER PUBLICATIONS A new Approach for Solvng Lnear Ordnary Dfferental Equatons Fawz Abdelwahd Department of

More information

Linear Algebra and its Applications

Linear Algebra and its Applications Lnear Algebra and ts Applcatons 4 (00) 5 56 Contents lsts avalable at ScenceDrect Lnear Algebra and ts Applcatons journal homepage: wwwelsevercom/locate/laa Notes on Hlbert and Cauchy matrces Mroslav Fedler

More information

Grover s Algorithm + Quantum Zeno Effect + Vaidman

Grover s Algorithm + Quantum Zeno Effect + Vaidman Grover s Algorthm + Quantum Zeno Effect + Vadman CS 294-2 Bomb 10/12/04 Fall 2004 Lecture 11 Grover s algorthm Recall that Grover s algorthm for searchng over a space of sze wors as follows: consder the

More information