u i log v i + 1 T (u v).

Size: px
Start display at page:

Download "u i log v i + 1 T (u v)."

Transcription

1 where f(v) = n v log v s the negatve entropy of v. Soluton. The negatve entropy s strctly convex and dfferentable on R n ++, hence f(u) > f(v) + f(v) T (u v) for all u, v R n ++ wth u v. Evaluatng both sdes of the nequalty, we obtan u log u > v log v + (log v + 1)(u v ) = u log v + 1 T (u v). Re-arrangng ths nequalty gves the desred result Convex-concave functons and saddle-ponts. We say the functon f : R n R m R s convex-concave f f(x, z) s a concave functon of z, for each fxed x, and a convex functon of x, for each fxed z. We also requre ts doman to have the product form dom f = A B, where A R n and B R m are convex. (a) Gve a second-order condton for a twce dfferentable functon f : R n R m R to be convex-concave, n terms of ts Hessan 2 f(x, z). (b) Suppose that f : R n R m R s convex-concave and dfferentable, wth f( x, z) = 0. Show that the saddle-pont property holds: for all x, z, we have f( x, z) f( x, z) f(x, z). Show that ths mples that f satsfes the strong max-mn property: sup z nf x f(x, z) = nf x sup z f(x, z) (and ther common value s f( x, z)). (c) Now suppose that f : R n R m R s dfferentable, but not necessarly convexconcave, and the saddle-pont property holds at x, z: Soluton. for all x, z. Show that f( x, z) = 0. f( x, z) f( x, z) f(x, z) (a) The condton follows drectly from the second-order condtons for convexty and concavty: t s 2 xxf(x, z) 0, 2 zzf(x, z) 0, for all x, z. In terms of 2 f, ths means that ts 1, 1 block s postve semdefnte, and ts 2, 2 block s negatve semdefnte. (b) Let us fx z. Snce xf( x, z) = 0 and f(x, z) s convex n x, we conclude that x mnmzes f(x, z) over x,.e., for all z, we have f( x, z) f(x, z). Ths s one of the nequaltes n the saddle-pont condton. We can argue n the same way about z. Fx x, and note that zf( x, z) = 0, together wth concavty of ths functon n z, means that z maxmzes the functon,.e., for any x we have f( x, z) f( x, z). (c) To establsh ths we argue the same way. If the saddle-pont condton holds, then x mnmzes f(x, z) over all x. Therefore we have f x( x, z) = 0. Smlarly, snce z maxmzes f( x, z) over all z, we have f z( x, z) = 0.

2 Exercses Examples 3.15 A famly of concave utlty functons. For 0 < α 1 let u α(x) = xα 1 α, wth dom u α = R +. We also defne u 0(x) = log x (wth dom u 0 = R ++). (a) Show that for x > 0, u 0(x) = lm α 0 u α(x). (b) Show that u α are concave, monotone ncreasng, and all satsfy u α(1) = 0. These functons are often used n economcs to model the beneft or utlty of some quantty of goods or money. Concavty of u α means that the margnal utlty (.e., the ncrease n utlty obtaned for a fxed ncrease n the goods) decreases as the amount of goods ncreases. In other words, concavty models the effect of sataton. Soluton. (a) In ths lmt, both the numerator and denomnator go to zero, so we use l Hoptal s rule: (d/dα)(x α 1) x α log x lm uα(x) = lm = lm = log x. α 0 α 0 (d/dα)α α 0 1 (b) By nspecton we have The dervatve s gven by u α(1) = 1α 1 α u α(x) = x α 1, = 0. whch s postve for all x (snce 0 < α < 1), so these functons are ncreasng. To show concavty, we examne the second dervatve: u α(x) = (α 1)x α 2. Snce ths s negatve for all x, we conclude that u α s strctly concave For each of the followng functons determne whether t s convex, concave, quasconvex, or quasconcave. (a) f(x) = e x 1 on R. Soluton. Strctly convex, and therefore quasconvex. Also quasconcave but not concave. (b) f(x 1, x 2) = x 1x 2 on R Soluton. The Hessan of f s 2 f(x) = whch s nether postve semdefnte nor negatve semdefnte. Therefore, f s nether convex nor concave. It s quasconcave, snce ts superlevel sets are convex. It s not quasconvex. (c) f(x 1, x 2) = 1/(x 1x 2) on R Soluton. The Hessan of f s {(x 1, x 2) R 2 ++ x 1x 2 α} 2 f(x) = 1 x 1x 2 2/(x 2 1 ) 1/(x 1x 2) 1/(x 1x 2) 2/x 2 2 Therefore, f s convex and quasconvex. It s not quasconcave or concave., 0

3 (d) f(x 1, x 2) = x 1/x 2 on R Soluton. The Hessan of f s 2 f(x) = 0 1/x 2 2 1/x 2 2 2x 1/x 3 2 whch s not postve or negatve semdefnte. Therefore, f s not convex or concave. It s quasconvex and quasconcave (.e., quaslnear), snce the sublevel and superlevel sets are halfspaces. (e) f(x 1, x 2) = x 2 1/x 2 on R R ++. Soluton. f s convex, as mentoned on page 72. (See also fgure 3.3). Ths s easly verfed by workng out the Hessan: 2 f(x) = 2/x 2 2x 1/x 2 2 2x 1/x 2 2 2x 2 1/x 3 2 = (2/x 2) 1 2x 1/x 2 1 2x1/x 2 0. Therefore, f s convex and quasconvex. It s not concave or quasconcave (see the fgure). (f) f(x 1, x 2) = x α 1 x 1 α 2, where 0 α 1, on R Soluton. Concave and quasconcave. The Hessan s 2 α(α 1)x α 2 1 x 1 α 2 α(1 α)x α 1 1 x α 2 f(x) = α(1 α)x α 1 1 x α 2 (1 α)( α)x α 1 x α 1 2 = α(1 α)x α 1 x 1 α 1/x 2 1 1/x 1x 2 2 1/x 1x 2 1/x 2 2 = α(1 α)x α 1 x 1 α 1/x1 1/x1 2 1/x 2 1/x 2 0. f s not convex or quasconvex Suppose p < 1, p 0. Show that the functon ( f(x) = x p ) 1/p wth dom f = R n ++ s concave. Ths ncludes as specal cases f(x) = ( n x1/2 ) 2 and the harmonc mean f(x) = ( n 1/x) 1. Hnt. Adapt the proofs for the log-sum-exp functon and the geometrc mean n Soluton. The frst dervatves of f are gven by f(x) ( ) 1 p f(x) = ( x p x )(1 p)/p x p 1 =. The second dervatves are 2 f(x) x x j = 1 p x ( ) p ( f(x) f(x) x x j ) 1 p x T ( ) = 1 p f(x) 2 1 p f(x) x x j for j, and 2 f(x) x 2 = 1 p f(x) ( ) f(x) 2 1 p ( 1 p f(x) x 2 x x ) 1 p.

4 whch s a nonnegatve sum of the convex functons x 1, x 1 + x 2, x 1 + x 2 + x 3,..., x 1 + x x r. (b) Let T (x, ω) denote the trgonometrc polynomal Show that the functon T (x, ω) = x 1 + x 2 cos ω + x 3 cos 2ω + + x n cos(n 1)ω. f(x) = 2π 0 log T (x, ω) dω s convex on {x R n T (x, ω) > 0, 0 ω 2π}. Soluton. The functon g(x, ω) = log(x 1 + x 2 cos ω + x 3 cos 2ω + + +x n cos(n 1)ω) s convex n x for fxed ω. Therefore s convex n x. f(x) = 2π 0 g(x, ω)dω 3.20 Composton wth an affne functon. Show that the followng functons f : R n R are convex. (a) f(x) = Ax b, where A R m n, b R m, and s a norm on R m. Soluton. f s the composton of a norm, whch s convex, and an affne functon. (b) f(x) = (det(a 0 + x 1A x na n)) 1/m, on {x A 0 + x 1A x na n 0}, where A S m. Soluton. f s the composton of the convex functon h(x) = (det X) 1/m and an affne transformaton. To see that h s convex on S m ++, we restrct h to a lne and prove that g(t) = det(z + tv ) 1/m s convex: g(t) = (det(z + tv )) 1/m = (det Z) 1/m (det(i + tz 1/2 V Z 1/2 )) 1/m m = (det Z) 1/m ( (1 + tλ )) 1/m where λ 1,..., λ m denote the egenvalues of Z 1/2 V Z 1/2. We have expressed g as the product of a negatve constant and the geometrc mean of 1 + tλ, = 1,..., m. Therefore g s convex. (See also exercse 3.18.) (c) f(x) = tr (A 0 + x 1A x na n) 1, on {x A 0 +x 1A 1 + +x na n 0}, where A S m. (Use the fact that tr(x 1 ) s convex on S m ++; see exercse 3.18.) Soluton. f s the composton of tr X 1 and an affne transformaton x A 0 + x 1A x na n Pontwse maxmum and supremum. Show that the followng functons f : R n R are convex.

5 Exercses (a) f(x) = max,...,k A () x b (), where A () R m n, b () R m and s a norm on R m. Soluton. f s the pontwse maxmum of k functons A () x b (). Each of those functons s convex because t s the composton of an affne transformaton and a norm. (b) f(x) = r x on R n, where x denotes the vector wth x = x (.e., x s the absolute value of x, componentwse), and x s the th largest component of x. In other words, x 1, x 2,..., x n are the absolute values of the components of x, sorted n nonncreasng order. Soluton. Wrte f as f(x) = r x = max 1 1 < 2 < < x x r r n whch s the pontwse maxmum of n!/(r!(n r)!) convex functons Composton rules. Show that the followng functons are convex. (a) f(x) = log( log( m eat x+b )) on dom f = {x m eat x+b < 1}. You can use the fact that log( n ey ) s convex. Soluton. g(x) = log( m eat x+b ) s convex (composton of the log-sum-exp functon and an affne mappng), so g s concave. The functon h(y) = log y s convex and decreasng. Therefore f(x) = h( g(x)) s convex. (b) f(x, u, v) = uv x T x on dom f = {(x, u, v) uv > x T x, u, v > 0}. Use the fact that x T x/u s convex n (x, u) for u > 0, and that x 1x 2 s convex on R Soluton. We can express f as f(x, u, v) = u(v x T x/u). The functon h(x 1, x 2) = x 1x 2 s convex on R 2 ++, and decreasng n each argument. The functons g 1(u, v, x) = u and g 2(u, v, x) = v x T x/u are concave. Therefore f(u, v, x) = h(g(u, v, x)) s convex. (c) f(x, u, v) = log(uv x T x) on dom f = {(x, u, v) uv > x T x, u, v > 0}. Soluton. We can express f as f(x, u, v) = log u log(v x T x/u). The frst term s convex. The functon v x T x/u s concave because v s lnear and x T x/u s convex on {(x, u) u > 0}. Therefore the second term n f s convex: t s the composton of a convex decreasng functon log t and a concave functon. (d) f(x, t) = (t p x p p) 1/p where p > 1 and dom f = {(x, t) t x p}. You can use the fact that x p p/u p 1 s convex n (x, u) for u > 0 (see exercse 3.23), and that x 1/p y 1 1/p s convex on R 2 + (see exercse 3.16). Soluton. We can express f as ( )) 1/p ( ) 1/p f(x, t) = (t p 1 t x p p = t 1 1/p t x p p. t p 1 t p 1 Ths s the composton of h(y 1, y 2) = y 1/p 1 y 1 1/p 2 (convex and decreasng n each argument) and two concave functons g 1(x, t) = t 1 1/p, g 2(x, t) = t x p p t p 1.

6 (e) f(x, t) = log(t p x p p) where p > 1 and dom f = {(x, t) t > x p}. You can use the fact that x p p/u p 1 s convex n (x, u) for u > 0 (see exercse 3.23). Soluton. Express f as f(x, t) = log t p 1 log(t x p p/t p 1 ) = (p 1) log t log(t x p p/t p 1 ). The frst term s convex. The second term s the composton of a decreasng convex functon and a concave functon, and s also convex Perspectve of a functon. (a) Show that for p > 1, f(x, t) = x1 p + + x n p t p 1 = x p p t p 1 s convex on {(x, t) t > 0}. Soluton. Ths s the perspectve functon of x p p = x 1 p + + x n p. (b) Show that f(x) = Ax + b 2 2 c T x + d s convex on {x c T x + d > 0}, where A R m n, b R m, c R n and d R. Soluton. Ths functon s the composton of the functon g(y, t) = y T y/t wth an affne transformaton (y, t) = (Ax + b, c T x + d). Therefore convexty of f follows from the fact that g s convex on {(y, t) t > 0}. For convexty of g one can note that t s the perspectve of x T x, or drectly verfy that the Hessan 2 I/t y/t 2 g(y, t) = y T /t y T y/t 3 s postve semdefnte, snce T v I/t y/t 2 v w y T /t y T y/t 3 w for all v and w. = tv yw 2 2/t Some functons on the probablty smplex. Let x be a real-valued random varable whch takes values n {a 1,..., a n} where a 1 < a 2 < < a n, wth prob(x = a ) = p, = 1,..., n. For each of the followng functons of p (on the probablty smplex {p R n + 1 T p = 1}), determne f the functon s convex, concave, quasconvex, or quasconcave. (a) E x. Soluton. E x = p 1a p na n s lnear, hence convex, concave, quasconvex, and quasconcave (b) prob(x α). Soluton. Let j = mn{ a α}. Then prob(x α) = n p, Ths s a lnear =j functon of p, hence convex, concave, quasconvex, and quasconcave. (c) prob(α x β). Soluton. Let j = mn{ a α} and k = max{ a β}. Then prob(α x β) = k p. Ths s a lnear functon of p, hence convex, concave, quasconvex, =j and quasconcave.

7 Exercses we have so p >q (p q ) = ( p q (p q ) d mp(p, q) = (1/2) (p q ) (1/2) (p q ) p >q p q = (1/2) p q = (1/2) p q 1. Ths makes t very clear that d mp s convex. The best way to nterpret ths result s as an nterpretaton of the l 1-norm for probablty dstrbutons. It states that the l 1-dstance between two probablty dstrbutons s twce the maxmum dfference n probablty, over all events, of the dstrbutons More functons of egenvalues. Let λ 1(X) λ 2(X) λ n(x) denote the egenvalues of a matrx X S n. We have already seen several functons of the egenvalues that are convex or concave functons of X. The maxmum egenvalue λ 1(X) s convex (example 3.10). The mnmum egenvalue λ n(x) s concave. The sum of the egenvalues (or trace), tr X = λ 1(X) + + λ n(x), s lnear. The sum of the nverses of the egenvalues (or trace of the nverse), tr(x 1 ) = n 1/λ(X), s convex on Sn ++ (exercse 3.18). The geometrc mean of the egenvalues, (det X) 1/n = ( n λ(x))1/n, and the logarthm of the product of the egenvalues, log det X = n log λ(x), are concave on X S n ++ (exercse 3.18 and page 74). In ths problem we explore some more functons of egenvalues, by explotng varatonal characterzatons. (a) Sum of k largest egenvalues. Show that k λ(x) s convex on Sn. Hnt. HJ85, page 191 Use the varatonal characterzaton k λ (X) = sup{tr(v T XV ) V R n k, V T V = I}. Soluton. The varatonal characterzaton shows that f s the pontwse supremum of a famly of lnear functons tr(v T XV ). (b) Geometrc mean of k smallest egenvalues. Show that ( n =n k+1 λ(x))1/k s concave on S n ++. Hnt. MO79, page 513 For X 0, we have ( n ) 1/k λ (X) = 1 k nf{tr(v T XV ) V R n k, det V T V = 1}. =n k+1 Soluton. f s the pontwse nfmum of a famly of lnear functons tr(v T XV ). (c) Log of product of k smallest egenvalues. Show that n log λ(x) s concave =n k+1 on S n ++. Hnt. MO79, page 513 For X 0, { n k } λ (X) = nf (V T XV ) V R n k, V T V = I. =n k+1 ),

8 Soluton. f s the pontwse nfmum of a famly of concave functons log (V T XV ) = log(v T XV ) Dagonal elements of Cholesky factor. Each X S n ++ has a unque Cholesky factorzaton X = LL T, where L s lower trangular, wth L > 0. Show that L s a concave functon of X (wth doman S n ++). Hnt. L can be expressed as L = (w z T Y 1 z) 1/2, where Y z z T s the leadng submatrx of X. Soluton. The functon f(z, Y ) = z T Y 1 z wth dom f = {(z, Y ) Y 0} s convex jontly n z and Y. To see ths note that Y z (z, Y, t) ep f Y 0, z T 0, t so ep f s a convex set. Therefore, w z T Y 1 z s a concave functon of X. Snce the squareroot s an ncreasng concave functon, t follows from the composton rules that l kk = (w z T Y 1 z) 1/2 s a concave functon of X. Operatons that preserve convexty 3.28 Expressng a convex functon as the pontwse supremum of a famly of affne functons. In ths problem we extend the result proved on page 83 to the case where dom f R n. Let f : R n R be a convex functon. Defne f : R n R as the pontwse supremum of all affne functons that are global underestmators of f: f(x) = sup{g(x) g affne, g(z) f(z) for all z}. (a) Show that f(x) = f(x) for x nt dom f. (b) Show that f = f f f s closed (.e., ep f s a closed set; see A.3.3). Soluton. (a) The pont (x, f(x)) s n the boundary of ep f. (If t were n nt ep f, then for small, postve ɛ we would have (x, f(x) ɛ) ep f, whch s mpossble.) From the results of 2.5.2, we know there s a supportng hyperplane to ep f at (x, f(x)),.e., a R n, b R such that w a T z + bt a T x + bf(x) for all (z, t) ep f. Snce t can be arbtrarly large f (z, t) ep f, we conclude that b 0. Suppose b = 0. Then a T z a T x for all z dom f whch contradcts x nt dom f. Therefore b > 0. Dvdng the above nequalty by b yelds t f(x) + (a/b) T (x z) for all (z, t) ep f. Therefore the affne functon g(z) = f(x) + (a/b) T (x z) s an affne global underestmator of f, and hence by defnton of f, f(x) f(x) g(x). However g(x) = f(x), so we must have f(x) = f(x).

9 (a) Let A = conv B. Snce B A, we obvously have S B(y) S A(y). Suppose we have strct nequalty for some y,.e., y T u < y T v for all u B and some v A. Ths leads to a contradcton, because by defnton v s the convex combnaton of a set of ponts u B,.e., v = θu, wth θ 0, θ = 1. Snce y T u < y T v for all, ths would mply y T v = θ y T u < θ y T v = y T v. We conclude that we must have equalty S B(y) = S A(y). (b) Follows from (c) Follows from S A+B(y) = sup{y T (u + v) u A, v B} = sup{y T u u A} + sup{y T v u B} = S A(y) + S B(y). S A B(y) = sup{y T u u A B} = max{sup{y T u u A}, sup{y T v u B} = max{s A(y), S B(y)}. (d) Obvously, f A B, then S A(y) S B(y) for all y. We need to show that f A B, then S A(y) > S B(y) for some y. Suppose A B. Consder a pont x A, x B. Snce B s closed and convex, x can be strctly separated from B by a hyperplane,.e., there s a y 0 such that y T x > y T x for all x B. It follows that S B(y) < y T x S A(y). Conjugate functons 3.36 Derve the conjugates of the followng functons. (a) Max functon. f(x) = max,...,n x on R n. Soluton. We wll show that { f 0 f y 0, 1 T y = 1 (y) = otherwse. We frst verfy the doman of f. Frst suppose y has a negatve component, say y k < 0. If we choose a vector x wth x k = t, x = 0 for k, and let t go to nfnty, we see that x T y max x = ty k, so y s not n dom f. Next, assume y 0 but 1 T y > 1. We choose x = t1 and let t go to nfnty, to show that x T y max x = t1 T y t

10 Exercses s unbounded above. Smlarly, when y 0 and 1 T y < 1, we choose x = t1 and let t go to nfnty. The remanng case for y s y 0 and 1 T y = 1. In ths case we have x T y max x for all x, and therefore x T y max x 0 for all x, wth equalty for x = 0. Therefore f (y) = 0. (b) Sum of largest elements. f(x) = r x on R n. Soluton. The conjugate s { f 0 0 y 1, 1 T y = r (y) = otherwse, We frst verfy the doman of f. Suppose y has a negatve component, say y k < 0. If we choose a vector x wth x k = t, x = 0 for k, and let t go to nfnty, we see that x T y f(x) = ty k, so y s not n dom f. Next, suppose y has a component greater than 1, say y k > 1. If we choose a vector x wth x k = t, x = 0 for k, and let t go to nfnty, we see that x T y f(x) = ty k t, so y s not n dom f. Fnally, assume that 1 T x r. We choose x = t1 and fnd that x T y f(x) = t1 T y tr s unbounded above, as t or t. If y satsfes all the condtons we have x T y f(x) for all x, wth equalty for x = 0. Therefore f (y) = 0. (c) Pecewse-lnear functon on R. f(x) = max,...,m(a x + b ) on R. You can assume that the a are sorted n ncreasng order,.e., a 1 a m, and that none of the functons a x + b s redundant,.e., for each k there s at least one x wth f(x) = a k x + b k. Soluton. Under the assumpton, the graph of f s a pecewse-lnear, wth breakponts (b b +1)/(a +1 a ), = 1,..., m 1. We can wrte f as ( f (y) = sup x ) xy max (ax + b),...,m We see that dom f = a 1, a m, snce for y outsde that range, the expresson nsde the supremum s unbounded above. For a y a +1, the supremum n the defnton of f s reached at the breakpont between the segments and + 1,.e., at the pont (b +1 b )/(a +1 a ), so we obtan f y a (y) = b (b +1 b ) a +1 a where s defned by a y a +1. Hence the graph of f s also a pecewse-lnear curve connectng the ponts (a, b ) for = 1,..., m. Geometrcally, the epgraph of f s the epgraphcal hull of the ponts (a, b ).

11 (d) Power functon. f(x) = x p on R ++, where p > 1. Repeat for p < 0. Soluton. We ll use standard notaton: we defne q by the equaton 1/p + 1/q = 1,.e., q = p/(p 1). We start wth the case p > 1. Then x p s strctly convex on R +. For y < 0 the functon yx x p acheves ts maxmum for x > 0 at x = 0, so f (y) = 0. For y > 0 the functon acheves ts maxmum at x = (y/p) 1/(p 1), where t has value y(y/p) 1/(p 1) (y/p) p/(p 1) = (p 1)(y/p) q. Therefore we have f (y) = { 0 y 0 (p 1)(y/p) q y > 0. For p < 0 smlar arguments show that dom f = R ++ and f (y) = p q ( y/p)q. (e) Geometrc mean. f(x) = ( x ) 1/n on R n ++. Soluton. The conjugate functon s { ( f (y) = 0 f y 0, ( y)) 1/n 1/n otherwse. We frst verfy the doman of f. Assume y has a postve component, say y k > 0. Then we can choose x k = t and x = 1, k, to show that x T y f(x) = ty k + y t 1/n s unbounded above as a functon of t > 0. Hence the condton y 0 s ndeed requred. Next assume that y 0, but ( ( y))1/n < 1/n. We choose x = t/y, and obtan ( 1/n x T y f(x) = tn t ( )) 1 y as t. Ths demonstrates that the second condton for the doman of f s also needed. Now assume that y 0 and ( ( y)) 1/n 1/n, and x 0. The arthmetcgeometrc mean nequalty states that k x T y n ( ) 1/n ( ) 1/n ( y x ) 1 x, n.e., x T y f(x) wth equalty for x = 1/y. Hence, f (y) = 0. (f) Negatve generalzed logarthm for second-order cone. f(x, t) = log(t 2 x T x) on {(x, t) R n R x 2 < t}. Soluton. f (y, u) = 2 + log 4 log(u 2 y T y), dom f = {(y, u) y 2 < u}. We frst verfy the doman. Suppose y 2 u. Choose x = sy, t = s( x 2 + 1) > s y 2 su, wth s 0. Then y T x + tu > sy T y su 2 = s(u 2 y T y) 0,

12 Exercses so y x + tu goes to nfnty, at a lnear rate, whle the functon log(t 2 x T x) goes to as log s. Therefore y T x + tu + log(t 2 x T x) s unbounded above. Next, assume that y 2 < u. Settng the dervatve of y T x + ut + log(t 2 x T x) wth respect to x and t equal to zero, and solvng for t and x we see that the maxmzer s 2y x = u 2 y T y, t = 2u u 2 y T y. Ths gves f (y, u) = ut + y T x + log(t 2 x T x) = 2 + log 4 log(y 2 u t u) Show that the conjugate of f(x) = tr(x 1 ) wth dom f = S n ++ s gven by f (Y ) = 2 tr( Y ) 1/2, dom f = S n +. Hnt. The gradent of f s f(x) = X 2. Soluton. We frst verfy the doman of f. Suppose Y has egenvalue decomposton Y = QΛQ T = λ q q T wth λ 1 > 0. Let X = Q dag(t, 1,..., 1)Q T = tq 1q1 T + n =2 qqt. We have tr XY tr X 1 = tλ 1 + λ 1/t (n 1), whch grows unboundedly as t. Therefore Y dom f. Next, assume Y 0. If Y 0, we can fnd the maxmum of =2 tr XY tr X 1 by settng the gradent equal to zero. We obtan Y = X 2,.e., X = ( Y ) 1/2, and f (Y ) = 2 tr( Y ) 1/2. Fnally we verfy that ths expresson remans vald when Y 0, but Y s sngular. Ths follows from the fact that conjugate functons are always closed,.e., have closed epgraphs Young s nequalty. Let f : R R be an ncreasng functon, wth f(0) = 0, and let g be ts nverse. Defne F and G as F (x) = x 0 f(a) da, G(y) = y 0 g(a) da. Show that F and G are conjugates. Gve a smple graphcal nterpretaton of Young s nequalty, xy F (x) + G(y). Soluton. The nequalty xy F (x) + G(y) has a smple geometrc meanng, llustrated below.

Solutions HW #2. minimize. Ax = b. Give the dual problem, and make the implicit equality constraints explicit. Solution.

Solutions HW #2. minimize. Ax = b. Give the dual problem, and make the implicit equality constraints explicit. Solution. Solutons HW #2 Dual of general LP. Fnd the dual functon of the LP mnmze subject to c T x Gx h Ax = b. Gve the dual problem, and make the mplct equalty constrants explct. Soluton. 1. The Lagrangan s L(x,

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

Math1110 (Spring 2009) Prelim 3 - Solutions

Math1110 (Spring 2009) Prelim 3 - Solutions Math 1110 (Sprng 2009) Solutons to Prelm 3 (04/21/2009) 1 Queston 1. (16 ponts) Short answer. Math1110 (Sprng 2009) Prelm 3 - Solutons x a 1 (a) (4 ponts) Please evaluate lm, where a and b are postve numbers.

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

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

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

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

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

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

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,

More information

First day August 1, Problems and Solutions

First day August 1, Problems and Solutions FOURTH INTERNATIONAL COMPETITION FOR UNIVERSITY STUDENTS IN MATHEMATICS July 30 August 4, 997, Plovdv, BULGARIA Frst day August, 997 Problems and Solutons Problem. Let {ε n } n= be a sequence of postve

More information

6) Derivatives, gradients and Hessian matrices

6) Derivatives, gradients and Hessian matrices 30C00300 Mathematcal Methods for Economsts (6 cr) 6) Dervatves, gradents and Hessan matrces Smon & Blume chapters: 14, 15 Sldes by: Tmo Kuosmanen 1 Outlne Defnton of dervatve functon Dervatve notatons

More information

6.854J / J Advanced Algorithms Fall 2008

6.854J / J Advanced Algorithms Fall 2008 MIT OpenCourseWare http://ocw.mt.edu 6.854J / 18.415J Advanced Algorthms Fall 2008 For nformaton about ctng these materals or our Terms of Use, vst: http://ocw.mt.edu/terms. 18.415/6.854 Advanced Algorthms

More information

Some basic inequalities. Definition. Let V be a vector space over the complex numbers. An inner product is given by a function, V V C

Some basic inequalities. Definition. Let V be a vector space over the complex numbers. An inner product is given by a function, V V C Some basc nequaltes Defnton. Let V be a vector space over the complex numbers. An nner product s gven by a functon, V V C (x, y) x, y satsfyng the followng propertes (for all x V, y V and c C) (1) x +

More information

Math 217 Fall 2013 Homework 2 Solutions

Math 217 Fall 2013 Homework 2 Solutions Math 17 Fall 013 Homework Solutons Due Thursday Sept. 6, 013 5pm Ths homework conssts of 6 problems of 5 ponts each. The total s 30. You need to fully justfy your answer prove that your functon ndeed has

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

Exercise Solutions to Real Analysis

Exercise Solutions to Real Analysis xercse Solutons to Real Analyss Note: References refer to H. L. Royden, Real Analyss xersze 1. Gven any set A any ɛ > 0, there s an open set O such that A O m O m A + ɛ. Soluton 1. If m A =, then there

More information

Lecture 20: Lift and Project, SDP Duality. Today we will study the Lift and Project method. Then we will prove the SDP duality theorem.

Lecture 20: Lift and Project, SDP Duality. Today we will study the Lift and Project method. Then we will prove the SDP duality theorem. prnceton u. sp 02 cos 598B: algorthms and complexty Lecture 20: Lft and Project, SDP Dualty Lecturer: Sanjeev Arora Scrbe:Yury Makarychev Today we wll study the Lft and Project method. Then we wll prove

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

Appendix B. Criterion of Riemann-Stieltjes Integrability

Appendix B. Criterion of Riemann-Stieltjes Integrability Appendx B. Crteron of Remann-Steltes Integrablty Ths note s complementary to [R, Ch. 6] and [T, Sec. 3.5]. The man result of ths note s Theorem B.3, whch provdes the necessary and suffcent condtons for

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

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

Spectral Graph Theory and its Applications September 16, Lecture 5

Spectral Graph Theory and its Applications September 16, Lecture 5 Spectral Graph Theory and ts Applcatons September 16, 2004 Lecturer: Danel A. Spelman Lecture 5 5.1 Introducton In ths lecture, we wll prove the followng theorem: Theorem 5.1.1. Let G be a planar graph

More information

( ) 2 ( ) ( ) Problem Set 4 Suggested Solutions. Problem 1

( ) 2 ( ) ( ) Problem Set 4 Suggested Solutions. Problem 1 Problem Set 4 Suggested Solutons Problem (A) The market demand functon s the soluton to the followng utlty-maxmzaton roblem (UMP): The Lagrangean: ( x, x, x ) = + max U x, x, x x x x st.. x + x + x y x,

More information

Bezier curves. Michael S. Floater. August 25, These notes provide an introduction to Bezier curves. i=0

Bezier curves. Michael S. Floater. August 25, These notes provide an introduction to Bezier curves. i=0 Bezer curves Mchael S. Floater August 25, 211 These notes provde an ntroducton to Bezer curves. 1 Bernsten polynomals Recall that a real polynomal of a real varable x R, wth degree n, s a functon of the

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

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

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

8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS

8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS SECTION 8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS 493 8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS All the vector spaces you have studed thus far n the text are real vector spaces because the scalars

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

Bézier curves. Michael S. Floater. September 10, These notes provide an introduction to Bézier curves. i=0

Bézier curves. Michael S. Floater. September 10, These notes provide an introduction to Bézier curves. i=0 Bézer curves Mchael S. Floater September 1, 215 These notes provde an ntroducton to Bézer curves. 1 Bernsten polynomals Recall that a real polynomal of a real varable x R, wth degree n, s a functon of

More information

REAL ANALYSIS I HOMEWORK 1

REAL ANALYSIS I HOMEWORK 1 REAL ANALYSIS I HOMEWORK CİHAN BAHRAN The questons are from Tao s text. Exercse 0.0.. If (x α ) α A s a collecton of numbers x α [0, + ] such that x α

More information

arxiv: v2 [math.ca] 24 Sep 2010

arxiv: v2 [math.ca] 24 Sep 2010 A Note on the Weghted Harmonc-Geometrc-Arthmetc Means Inequaltes arxv:0900948v2 [mathca] 24 Sep 200 Gérard Maze, Urs Wagner e-mal: {gmaze,uwagner}@mathuzhch Mathematcs Insttute Unversty of Zürch Wnterthurerstr

More information

Math 426: Probability MWF 1pm, Gasson 310 Homework 4 Selected Solutions

Math 426: Probability MWF 1pm, Gasson 310 Homework 4 Selected Solutions Exercses from Ross, 3, : Math 26: Probablty MWF pm, Gasson 30 Homework Selected Solutons 3, p. 05 Problems 76, 86 3, p. 06 Theoretcal exercses 3, 6, p. 63 Problems 5, 0, 20, p. 69 Theoretcal exercses 2,

More information

= = = (a) Use the MATLAB command rref to solve the system. (b) Let A be the coefficient matrix and B be the right-hand side of the system.

= = = (a) Use the MATLAB command rref to solve the system. (b) Let A be the coefficient matrix and B be the right-hand side of the system. Chapter Matlab Exercses Chapter Matlab Exercses. Consder the lnear system of Example n Secton.. x x x y z y y z (a) Use the MATLAB command rref to solve the system. (b) Let A be the coeffcent matrx and

More information

CSCE 790S Background Results

CSCE 790S Background Results CSCE 790S Background Results Stephen A. Fenner September 8, 011 Abstract These results are background to the course CSCE 790S/CSCE 790B, Quantum Computaton and Informaton (Sprng 007 and Fall 011). Each

More information

APPROXIMATE PRICES OF BASKET AND ASIAN OPTIONS DUPONT OLIVIER. Premia 14

APPROXIMATE PRICES OF BASKET AND ASIAN OPTIONS DUPONT OLIVIER. Premia 14 APPROXIMAE PRICES OF BASKE AND ASIAN OPIONS DUPON OLIVIER Prema 14 Contents Introducton 1 1. Framewor 1 1.1. Baset optons 1.. Asan optons. Computng the prce 3. Lower bound 3.1. Closed formula for the prce

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

Topic 5: Non-Linear Regression

Topic 5: Non-Linear Regression Topc 5: Non-Lnear Regresson The models we ve worked wth so far have been lnear n the parameters. They ve been of the form: y = Xβ + ε Many models based on economc theory are actually non-lnear n the parameters.

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

MA 323 Geometric Modelling Course Notes: Day 13 Bezier Curves & Bernstein Polynomials

MA 323 Geometric Modelling Course Notes: Day 13 Bezier Curves & Bernstein Polynomials MA 323 Geometrc Modellng Course Notes: Day 13 Bezer Curves & Bernsten Polynomals Davd L. Fnn Over the past few days, we have looked at de Casteljau s algorthm for generatng a polynomal curve, and we have

More information

Georgia Tech PHYS 6124 Mathematical Methods of Physics I

Georgia Tech PHYS 6124 Mathematical Methods of Physics I Georga Tech PHYS 624 Mathematcal Methods of Physcs I Instructor: Predrag Cvtanovć Fall semester 202 Homework Set #7 due October 30 202 == show all your work for maxmum credt == put labels ttle legends

More information

Math 702 Midterm Exam Solutions

Math 702 Midterm Exam Solutions Math 702 Mdterm xam Solutons The terms measurable, measure, ntegrable, and almost everywhere (a.e.) n a ucldean space always refer to Lebesgue measure m. Problem. [6 pts] In each case, prove the statement

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

The Second Anti-Mathima on Game Theory

The Second Anti-Mathima on Game Theory The Second Ant-Mathma on Game Theory Ath. Kehagas December 1 2006 1 Introducton In ths note we wll examne the noton of game equlbrum for three types of games 1. 2-player 2-acton zero-sum games 2. 2-player

More information

Problem Do any of the following determine homomorphisms from GL n (C) to GL n (C)?

Problem Do any of the following determine homomorphisms from GL n (C) to GL n (C)? Homework 8 solutons. Problem 16.1. Whch of the followng defne homomomorphsms from C\{0} to C\{0}? Answer. a) f 1 : z z Yes, f 1 s a homomorphsm. We have that z s the complex conjugate of z. If z 1,z 2

More information

Another converse of Jensen s inequality

Another converse of Jensen s inequality Another converse of Jensen s nequalty Slavko Smc Abstract. We gve the best possble global bounds for a form of dscrete Jensen s nequalty. By some examples ts frutfulness s shown. 1. Introducton Throughout

More information

Math 261 Exercise sheet 2

Math 261 Exercise sheet 2 Math 261 Exercse sheet 2 http://staff.aub.edu.lb/~nm116/teachng/2017/math261/ndex.html Verson: September 25, 2017 Answers are due for Monday 25 September, 11AM. The use of calculators s allowed. Exercse

More information

Chapter 7 Generalized and Weighted Least Squares Estimation. In this method, the deviation between the observed and expected values of

Chapter 7 Generalized and Weighted Least Squares Estimation. In this method, the deviation between the observed and expected values of Chapter 7 Generalzed and Weghted Least Squares Estmaton The usual lnear regresson model assumes that all the random error components are dentcally and ndependently dstrbuted wth constant varance. When

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

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

SELECTED PROOFS. DeMorgan s formulas: The first one is clear from Venn diagram, or the following truth table:

SELECTED PROOFS. DeMorgan s formulas: The first one is clear from Venn diagram, or the following truth table: SELECTED PROOFS DeMorgan s formulas: The frst one s clear from Venn dagram, or the followng truth table: A B A B A B Ā B Ā B T T T F F F F T F T F F T F F T T F T F F F F F T T T T The second one can be

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

Calculus of Variations Basics

Calculus of Variations Basics Chapter 1 Calculus of Varatons Bascs 1.1 Varaton of a General Functonal In ths chapter, we derve the general formula for the varaton of a functonal of the form J [y 1,y 2,,y n ] F x,y 1,y 2,,y n,y 1,y

More information

Lecture Notes on Linear Regression

Lecture Notes on Linear Regression Lecture Notes on Lnear Regresson Feng L fl@sdueducn Shandong Unversty, Chna Lnear Regresson Problem In regresson problem, we am at predct a contnuous target value gven an nput feature vector We assume

More information

1 Matrix representations of canonical matrices

1 Matrix representations of canonical matrices 1 Matrx representatons of canoncal matrces 2-d rotaton around the orgn: ( ) cos θ sn θ R 0 = sn θ cos θ 3-d rotaton around the x-axs: R x = 1 0 0 0 cos θ sn θ 0 sn θ cos θ 3-d rotaton around the y-axs:

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

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

Mathematical Economics MEMF e ME. Filomena Garcia. Topic 2 Calculus

Mathematical Economics MEMF e ME. Filomena Garcia. Topic 2 Calculus Mathematcal Economcs MEMF e ME Flomena Garca Topc 2 Calculus Mathematcal Economcs - www.seg.utl.pt/~garca/economa_matematca . Unvarate Calculus Calculus Functons : X Y y ( gves or each element X one element

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

SELECTED SOLUTIONS, SECTION (Weak duality) Prove that the primal and dual values p and d defined by equations (4.3.2) and (4.3.3) satisfy p d.

SELECTED SOLUTIONS, SECTION (Weak duality) Prove that the primal and dual values p and d defined by equations (4.3.2) and (4.3.3) satisfy p d. SELECTED SOLUTIONS, SECTION 4.3 1. Weak dualty Prove that the prmal and dual values p and d defned by equatons 4.3. and 4.3.3 satsfy p d. We consder an optmzaton problem of the form The Lagrangan for ths

More information

MATH 281A: Homework #6

MATH 281A: Homework #6 MATH 28A: Homework #6 Jongha Ryu Due date: November 8, 206 Problem. (Problem 2..2. Soluton. If X,..., X n Bern(p, then T = X s a complete suffcent statstc. Our target s g(p = p, and the nave guess suggested

More information

k t+1 + c t A t k t, t=0

k t+1 + c t A t k t, t=0 Macro II (UC3M, MA/PhD Econ) Professor: Matthas Kredler Fnal Exam 6 May 208 You have 50 mnutes to complete the exam There are 80 ponts n total The exam has 4 pages If somethng n the queston s unclear,

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

MATH Sensitivity of Eigenvalue Problems

MATH Sensitivity of Eigenvalue Problems MATH 537- Senstvty of Egenvalue Problems Prelmnares Let A be an n n matrx, and let λ be an egenvalue of A, correspondngly there are vectors x and y such that Ax = λx and y H A = λy H Then x s called A

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

The Geometry of Logit and Probit

The Geometry of Logit and Probit The Geometry of Logt and Probt Ths short note s meant as a supplement to Chapters and 3 of Spatal Models of Parlamentary Votng and the notaton and reference to fgures n the text below s to those two chapters.

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

How Strong Are Weak Patents? Joseph Farrell and Carl Shapiro. Supplementary Material Licensing Probabilistic Patents to Cournot Oligopolists *

How Strong Are Weak Patents? Joseph Farrell and Carl Shapiro. Supplementary Material Licensing Probabilistic Patents to Cournot Oligopolists * How Strong Are Weak Patents? Joseph Farrell and Carl Shapro Supplementary Materal Lcensng Probablstc Patents to Cournot Olgopolsts * September 007 We study here the specal case n whch downstream competton

More information

On the Multicriteria Integer Network Flow Problem

On the Multicriteria Integer Network Flow Problem BULGARIAN ACADEMY OF SCIENCES CYBERNETICS AND INFORMATION TECHNOLOGIES Volume 5, No 2 Sofa 2005 On the Multcrtera Integer Network Flow Problem Vassl Vasslev, Marana Nkolova, Maryana Vassleva Insttute of

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

A 2D Bounded Linear Program (H,c) 2D Linear Programming

A 2D Bounded Linear Program (H,c) 2D Linear Programming A 2D Bounded Lnear Program (H,c) h 3 v h 8 h 5 c h 4 h h 6 h 7 h 2 2D Lnear Programmng C s a polygonal regon, the ntersecton of n halfplanes. (H, c) s nfeasble, as C s empty. Feasble regon C s unbounded

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

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

PHYS 705: Classical Mechanics. Calculus of Variations II

PHYS 705: Classical Mechanics. Calculus of Variations II 1 PHYS 705: Classcal Mechancs Calculus of Varatons II 2 Calculus of Varatons: Generalzaton (no constrant yet) Suppose now that F depends on several dependent varables : We need to fnd such that has a statonary

More information

Supplement: Proofs and Technical Details for The Solution Path of the Generalized Lasso

Supplement: Proofs and Technical Details for The Solution Path of the Generalized Lasso Supplement: Proofs and Techncal Detals for The Soluton Path of the Generalzed Lasso Ryan J. Tbshran Jonathan Taylor In ths document we gve supplementary detals to the paper The Soluton Path of the Generalzed

More information

f(x,y) = (4(x 2 4)x,2y) = 0 H(x,y) =

f(x,y) = (4(x 2 4)x,2y) = 0 H(x,y) = Problem Set 3: Unconstraned mzaton n R N. () Fnd all crtcal ponts of f(x,y) (x 4) +y and show whch are ma and whch are mnma. () Fnd all crtcal ponts of f(x,y) (y x ) x and show whch are ma and whch are

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

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

Module 2. Random Processes. Version 2 ECE IIT, Kharagpur

Module 2. Random Processes. Version 2 ECE IIT, Kharagpur Module Random Processes Lesson 6 Functons of Random Varables After readng ths lesson, ou wll learn about cdf of functon of a random varable. Formula for determnng the pdf of a random varable. Let, X be

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

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

Causal Diamonds. M. Aghili, L. Bombelli, B. Pilgrim

Causal Diamonds. M. Aghili, L. Bombelli, B. Pilgrim Causal Damonds M. Aghl, L. Bombell, B. Plgrm Introducton The correcton to volume of a causal nterval due to curvature of spacetme has been done by Myrhem [] and recently by Gbbons & Solodukhn [] and later

More information

332600_08_1.qxp 4/17/08 11:29 AM Page 481

332600_08_1.qxp 4/17/08 11:29 AM Page 481 336_8_.qxp 4/7/8 :9 AM Page 48 8 Complex Vector Spaces 8. Complex Numbers 8. Conjugates and Dvson of Complex Numbers 8.3 Polar Form and DeMovre s Theorem 8.4 Complex Vector Spaces and Inner Products 8.5

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

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

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

763622S ADVANCED QUANTUM MECHANICS Solution Set 1 Spring c n a n. c n 2 = 1.

763622S ADVANCED QUANTUM MECHANICS Solution Set 1 Spring c n a n. c n 2 = 1. 7636S ADVANCED QUANTUM MECHANICS Soluton Set 1 Sprng 013 1 Warm-up Show that the egenvalues of a Hermtan operator  are real and that the egenkets correspondng to dfferent egenvalues are orthogonal (b)

More information

14 Lagrange Multipliers

14 Lagrange Multipliers Lagrange Multplers 14 Lagrange Multplers The Method of Lagrange Multplers s a powerful technque for constraned optmzaton. Whle t has applcatons far beyond machne learnng t was orgnally developed to solve

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

THE CHINESE REMAINDER THEOREM. We should thank the Chinese for their wonderful remainder theorem. Glenn Stevens

THE CHINESE REMAINDER THEOREM. We should thank the Chinese for their wonderful remainder theorem. Glenn Stevens THE CHINESE REMAINDER THEOREM KEITH CONRAD We should thank the Chnese for ther wonderful remander theorem. Glenn Stevens 1. Introducton The Chnese remander theorem says we can unquely solve any par of

More information

Complex Numbers. x = B B 2 4AC 2A. or x = x = 2 ± 4 4 (1) (5) 2 (1)

Complex Numbers. x = B B 2 4AC 2A. or x = x = 2 ± 4 4 (1) (5) 2 (1) Complex Numbers If you have not yet encountered complex numbers, you wll soon do so n the process of solvng quadratc equatons. The general quadratc equaton Ax + Bx + C 0 has solutons x B + B 4AC A For

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

e - c o m p a n i o n

e - c o m p a n i o n OPERATIONS RESEARCH http://dxdoorg/0287/opre007ec e - c o m p a n o n ONLY AVAILABLE IN ELECTRONIC FORM 202 INFORMS Electronc Companon Generalzed Quantty Competton for Multple Products and Loss of Effcency

More information

Chat eld, C. and A.J.Collins, Introduction to multivariate analysis. Chapman & Hall, 1980

Chat eld, C. and A.J.Collins, Introduction to multivariate analysis. Chapman & Hall, 1980 MT07: Multvarate Statstcal Methods Mke Tso: emal mke.tso@manchester.ac.uk Webpage for notes: http://www.maths.manchester.ac.uk/~mkt/new_teachng.htm. Introducton to multvarate data. Books Chat eld, C. and

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

DIFFERENTIAL FORMS BRIAN OSSERMAN

DIFFERENTIAL FORMS BRIAN OSSERMAN DIFFERENTIAL FORMS BRIAN OSSERMAN Dfferentals are an mportant topc n algebrac geometry, allowng the use of some classcal geometrc arguments n the context of varetes over any feld. We wll use them to defne

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