1 / 15

I.4 Polyhedral Theory

I.4 Polyhedral Theory.

avian
Download Presentation

I.4 Polyhedral Theory

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. I.4 Polyhedral Theory

  2. Objective of Study: want to know how to describe the convex hull of the solution set to the IP problem S = {xZ+n: Ax  b} using linear inequalities (at least approximately). Which inequalities are necessary? How can we identify the inequalities necessary (or strong) to describe the convex hull? • Def 1.1: Given S  Rn, x  Rn is a convex combination of points of S if there exists {xi}i=1t  S such that x = i=1t ixi, i=1t i = 1,   R+t. Convex hull of S: set of all convex combination of points in S (inside description) Intersection of all convex sets containing S (outside description) • linearly dependent: {x1, … , xk} some xi can be described as a linear combination of the remaining xj‘s. linearly independent: x1, … , xkRn linearly independent  i=1k ixi = 0 implies i = 0 i. (i.e. not linearly dependent)

  3. Subspace; The set closed under addition and scalar multiplication (or arbitrary linear combination) of elements in the set. Given a set A  Rn, Inside description: (subspace generated by A  Rn) Outside description: linear hull of A  Rn, (intersection of all subspaces containing A) • Prop: x1, … , xk  Rn linearly independent and x0 =  ixi, Then (1) all I’s unique and (2) {x1, … , xk}{x0}\{xj} linearly independent  j  0. • Def: Rank of a matrix A: mn : maximum number of linearly independent rows of A (= maximum number of linearly independent columns of A) Basis of A  Rn : linearly independent subset of A which generates all of A. (minimal generating set in A, maximal independent set in A) Rank of A  Rn : basis size

  4. Basis equicardinality property • Prop 1.2: The following statements are equivalent: a) {xRn: Ax = b}   b) rank(A) = rank(A, b) •  ixi is an affine combination if  i = 1. Affine space: the set closed under affine combination, i.e. x1, … , xk L, i=1k i = 1   ixi  L. • Note: affine space L = S + {a}, for some a  L, S: subspace. • Constrained form: subspace S = {x: Ax = 0}, affine space L = {x: Ax b} Affine span, affine hull of A  Rn. • Def: x1, … , xk  Rnaffinely dependent if some xi can be expressed as an affine combination of remaining xj, j  i. Otherwise affinely independent. • Def 1.4: x1, … , xk  Rnaffinelyindependent if the unique solution of i=1k ixi = 0, i=1k i= 0 is i = 0 for i = 1, … , k.

  5. Prop 1.3: The following statements are equivalent: • x1, … , xkRn are affinely independent. • x2 – x1, … , xk – x1 are linearly independent. • (x1, -1), … , (xk, -1)  Rn+1 are linearly independent. • Def: • Affine rank of A  Rn is the maximum number of affinely independent points in A. • dim(L) = dim(S), L is affine space and L = S + {a} • Maximum number of affinely independent points in Rn is n+1. (n linearly independent points + 0 vector)

  6. Prop 1.4: If {xRn: Ax = b}  , maximum number of affinely independent solutions of Ax = b is n+1-rank(A). (Compare with rank(A) + nullity(A) = n) • Solution of Ax = b is translation of solution of Ax = 0. Solution set of Ax = 0 is null space (orthogonal subspace) of rows of A whose dimension is n – rank(A)  affine rank is (n+1) – rank(A). • Def: p  Rn, H subspace, then the projection of p on H is q  H such that p-q  H. S  Rn, the projection of S on H is denoted by projH(S) = {q: q is projection of p on H for some p  S}.

  7. 2.Definitions of Polyhedra and Dimension • Def: Polyhedron is the set of points that satisfy a finite number of linear inequalities, i.e. P = {xRn : Ax  b}. (outside description) Bounded polyhedron is polytope (convex hull of finitely many points) T  Rn is a convex set if x1, x2  T implies that x1 + (1- )x2  T for all 0    1. Cone C  Rn : x  C  x  C,   R+1 (we only consider convex, polyhedral cones) • Prop: Polyhedron is a convex set. • Prop: P = {xRn : Ax  0} is a cone. • Def: A polyhedron P is of dimension k, denoted dim(P) = k, if the maximum number of affinely independent points in P is k+1. (dimension of the smallest affine space containing P.) • Def: A polyhedron P  Rn is full-dimensional if dim(P) = n.

  8. Notation: M = {1, … , m}, m: number of constraints M= = {iM: aix = bi,  xP} M = {iM: aix< bi, for some xP} = M\M= (A=, b=), (A, b) are corresponding rows of (A, b) P = {xRn: A=x = b=, Ax  b}       • Note that if i M, then (ai, bi) cannot be written as a linear combination of the rows of (A=, b=). • Def: inner point x: aix< bi, for all i M • Def: interior point x: aix < bi, for all i M • Prop 2.3: Every nonempty polyhedron P has an inner point. Pf) If M = , every point of P is inner. For each i M, there exists xiP such that aixi< bi. Let x* = (1/|M|)iM xi  P, then akx*< bk, for all k  M, hence inner point. 

  9. Prop 2.4: P  Rn (P  ), then dim(P) + rank(A=, b=) = n. (see Prop 1.4) Pf) Suppose rank(A=) = rank(A=, b=) = n-k, 0  k  n.  dim{x: A=x = 0} = k   k linearly independent points, say y1, … , yk. Let x* be an inner point (existence guaranteed) • x* + yi  P for small  > 0 and x*, x* + y1, … , x* + yk, affinely independent (since (x* + yi) – x* linearly independent) • dim(P)  k  dim(P) + rank(A=, b=)  n Now suppose dim(P) = k and x0, x1, … , xk are affinelyindep. points of P.  xi – x0 are linearly independent and A=(xi – x0) = 0  nullity(A=)  k  rank(A=) = rank(A=, b=)  n-k  dim(P) + rank(A=, b=)  n  • Cor: P is full-dimensional if and only if P has an interior point.

  10. 3. Describing Polyhedra by Facets • Def: x  0 [or (, 0)] is called a valid inequality for P if it is satisfied by all x  P. • (, 0) valid if and only if max{x: x  P}  0 • Def: (, 0) valid. F = {x  P: x = 0} is called a face of P and (, 0) represents F. F is said to be proper if F   and F  P. • F    max{x: x  P} = 0 If F  , we say that (, 0) supports P. • Prop: F  P nonempty face of P   c  Rn such that cx is maximized over P precisely on F.

  11. Prop 3.1: P = {xRn : Ax  b} with equality set M=  M, F is a nonempty face of P. Then, F is a polyhedron and F = {xRn : aix= bi, iMF=, aix  bi, iMF}, where MF=  M=, MF = M \ M=. The number of distinct faces of P is finite. Pf) Let F be the set of optimal solutions to 0 = max{ x: Ax  b}. Let u* be an optimal solution to min{ ub: uA = , u  0}, I* = { i: ui* > 0}. Let F* = {xRn : aix= bi, iI*, aix  bi, iM \ I* } Show F = F* 1) Suppose x  F*, then x = u*Ax = iI* ui*aix = iI* ui*bi = 0 (u*A =  from u* dual feasible)  x  F, hence F*  F. 2) Suppose x  P \ F*, then akx bk for some k  I*. • x = iI* ui*aix< iI* ui*bi = 0  x  F, hence F  F*. (showed x F*  x  F, i.e. x  F  x  F* ) From 1), 2), F = F*, F polyhedron. Since F  P, the equality set (AF, bF) of F must have the required property. Finally, since M is finite, possible equality set MF= is finite, so the number of distinct faces is finite. 

  12. Def: F is a facet if dim(F) = dim(P) – 1. • Prop 3.2: If F is a facet of P,  akx  bk, k  M representing F. Pf) dim(F) = dim(P) - 1  rank(AF=, bF=) = rank(A=, b=) + 1.  • Prop 3.3: For each facet F of P, one of the inequalities representing F is necessary in the description of P. Pf) Let PF be the polyhedron obtained by dropping inequalities representing F. Show PF \ P  . Let x* be an inner point of F, and arx  br be an inequality representing F. ar linearly independent of rows of A=  does not exist x such that xA= = ar  By thm of alternatives,  y such that A=y = 0, ary > 0 (thm of alternatives for subspaces) x* inner point of F  aix* < bi,  i  M \{inequalities representing F} Now ai(x* + y) = aix* + aiy = bi,  i  M= ai(x* + y) = aix* + aiy< bi,  i  M \{inequalities representing F} ar(x* + y) = arx* + ary > br  x* + y  PF \ P for small  > 0.  (For pf of thm of alt, may consider (P) min 0x, xA = ar, (D) max ary, Ay = 0)

  13. Prop 3.4: Every inequality representing a face of P of dimension less than dim(P) – 1 is irrelevant to the description of P. • When two inequalities (1, 01) and (2, 02) are equivalent in the description of P? {x: A=x = b=, x  0} = {x: A=x = b=, ( + A=)x  0+ b=} >0,  R|M=| Hence equivalent if (2, 02) = (1, 01) + (A=, b=) for some >0, R|M=| • Thm 3.5: • P full-dimensional  P has a unique representation (to within positive scalar multiplication) by a finite set of inequalities. • If dim(P) = n-k, k>0, then P = {xRn : aix= bi, i = 1, … , k, aix bi, i = k+1, … , k+t}. For i= 1, … , k, (ai,bi) are a maximal set of linearly independent rows of (A=, b=) For i = k+1, … , k+t, (ai,bi) is any inequality from the equivalence class of inequalities representing Fi.

  14. Thm 3.6: F = {x  P: x = 0} proper face of P. Then the following two statements are equivalent: • F is a facet of P. • If x = 0  x  F, then (, 0) = (A=, 0 + b=) for some   R1 and   R|M=| Pf) 2)  1) : Let L = {(, 0)Rn+1: (, 0) satisfies (, 0)=(  A=, 0 + b=)} (generated set) L’ = {(, 0)Rn+1: x = 0,  x  F} (constrained set) L  L’ since x  A=x = 0 + b=  x  F. By hypothesis of 2), L’  L  L = L’ (L, L’ subspaces) Suppose dim(P) = n –k  rank(A=, b=) = k  dim(L) = k+1 (F proper face  (, 0) linearly independent of (A=, b=)) Suppose x1, … , xr maximal affinely independent points in F. (continued)

  15. D = rank r Consider (, 0)DT = 0 Maximum number of affinely independent solution is (n+1) + 1 – rank(D) = n+2-r  dim(L’) = n+1-r = k+1  r = n-k  F is a facet. 1)  2): As above, L  L’. Show L = L’ Suppose dim(P) = (n-k), F facet   n-k affinely independent points in F  Similarly dim(L’) = k+1, dim(L) = k+1 and L  L’  L = L’ 

More Related