Numerical Optimization

Similar documents
EC 521 MATHEMATICAL METHODS FOR ECONOMICS. Lecture 2: Convex Sets

Math 5593 Linear Programming Lecture Notes

Lecture 2 - Introduction to Polytopes

Lecture 2. Topology of Sets in R n. August 27, 2008

Lecture 5: Duality Theory

Division of the Humanities and Social Sciences. Convex Analysis and Economic Theory Winter Separation theorems

Convex Sets. CSCI5254: Convex Optimization & Its Applications. subspaces, affine sets, and convex sets. operations that preserve convexity

However, this is not always true! For example, this fails if both A and B are closed and unbounded (find an example).

Convex Sets. Pontus Giselsson

Convex Geometry arising in Optimization

Conic Duality. yyye

Convexity: an introduction

REVIEW OF FUZZY SETS

Locally convex topological vector spaces

Chapter 4 Concepts from Geometry

FACES OF CONVEX SETS

Advanced Operations Research Techniques IE316. Quiz 1 Review. Dr. Ted Ralphs

2. Convex sets. x 1. x 2. affine set: contains the line through any two distinct points in the set

Math 414 Lecture 2 Everyone have a laptop?

Lecture 3: Convex sets

2. Convex sets. affine and convex sets. some important examples. operations that preserve convexity. generalized inequalities

maximize c, x subject to Ax b,

POLYHEDRAL GEOMETRY. Convex functions and sets. Mathematical Programming Niels Lauritzen Recall that a subset C R n is convex if

60 2 Convex sets. {x a T x b} {x ã T x b}

Lecture 2 September 3

Convex Optimization. 2. Convex Sets. Prof. Ying Cui. Department of Electrical Engineering Shanghai Jiao Tong University. SJTU Ying Cui 1 / 33

Convex Optimization. Convex Sets. ENSAE: Optimisation 1/24

Lecture 19 Subgradient Methods. November 5, 2008

CS522: Advanced Algorithms

Combinatorial Geometry & Topology arising in Game Theory and Optimization

Shiqian Ma, MAT-258A: Numerical Optimization 1. Chapter 2. Convex Optimization

MTAEA Convexity and Quasiconvexity

Convex Optimization - Chapter 1-2. Xiangru Lian August 28, 2015

Convex Sets (cont.) Convex Functions

Lecture 6: Faces, Facets

Topological properties of convex sets

MA4254: Discrete Optimization. Defeng Sun. Department of Mathematics National University of Singapore Office: S Telephone:

CS675: Convex and Combinatorial Optimization Spring 2018 Convex Sets. Instructor: Shaddin Dughmi

Lecture 4: Convexity

Lecture 2: August 29, 2018

Convex Optimization Lecture 2

2 Geometry Solutions

COM Optimization for Communications Summary: Convex Sets and Convex Functions

Linear Programming in Small Dimensions

EE/ACM Applications of Convex Optimization in Signal Processing and Communications Lecture 6

Lecture 19: Convex Non-Smooth Optimization. April 2, 2007

Open and Closed Sets

THREE LECTURES ON BASIC TOPOLOGY. 1. Basic notions.

Lecture 2: August 29, 2018

DM545 Linear and Integer Programming. Lecture 2. The Simplex Method. Marco Chiarandini

Intersection Cuts with Infinite Split Rank

Convex sets and convex functions

Introductory Operations Research

Convex sets and convex functions

Polar Duality and Farkas Lemma

Introduction to Modern Control Systems

Lecture 2 Convex Sets

Simplex Algorithm in 1 Slide

Point Cloud Simplification and Processing for Path-Planning. Master s Thesis in Engineering Mathematics and Computational Science DAVID ERIKSSON

CS 473: Algorithms. Ruta Mehta. Spring University of Illinois, Urbana-Champaign. Ruta (UIUC) CS473 1 Spring / 29

CS599: Convex and Combinatorial Optimization Fall 2013 Lecture 4: Convex Sets. Instructor: Shaddin Dughmi

Mathematical and Algorithmic Foundations Linear Programming and Matchings

CMU-Q Lecture 9: Optimization II: Constrained,Unconstrained Optimization Convex optimization. Teacher: Gianni A. Di Caro

EXTREME POINTS AND AFFINE EQUIVALENCE

MATH 890 HOMEWORK 2 DAVID MEREDITH

A mini-introduction to convexity

CS 372: Computational Geometry Lecture 10 Linear Programming in Fixed Dimension

11 Linear Programming

Lecture 15: The subspace topology, Closed sets

Lecture 3. Corner Polyhedron, Intersection Cuts, Maximal Lattice-Free Convex Sets. Tepper School of Business Carnegie Mellon University, Pittsburgh

Aspects of Convex, Nonconvex, and Geometric Optimization (Lecture 1) Suvrit Sra Massachusetts Institute of Technology

Lecture 4: Rational IPs, Polyhedron, Decomposition Theorem

Dubna 2018: lines on cubic surfaces

Section 5 Convex Optimisation 1. W. Dai (IC) EE4.66 Data Proc. Convex Optimisation page 5-1

Mathematical Programming and Research Methods (Part II)

LP Geometry: outline. A general LP. minimize x c T x s.t. a T i. x b i, i 2 M 1 a T i x = b i, i 2 M 3 x j 0, j 2 N 1. where

1 Linear Programming. 1.1 Optimizion problems and convex polytopes 1 LINEAR PROGRAMMING

Convexity. 1 X i is convex. = b is a hyperplane in R n, and is denoted H(p, b) i.e.,

Lecture 5: Properties of convex sets

AMS : Combinatorial Optimization Homework Problems - Week V

Lecture: Convex Sets

Walheer Barnabé. Topics in Mathematics Practical Session 2 - Topology & Convex

In this chapter we introduce some of the basic concepts that will be useful for the study of integer programming problems.

Linear programming and duality theory

Polytopes Course Notes

Graph Connectivity in Sparse Subspace Clustering

Topology 550A Homework 3, Week 3 (Corrections: February 22, 2012)

College of Computer & Information Science Fall 2007 Northeastern University 14 September 2007

Interval-Vector Polytopes

Three Dimensional Geometry. Linear Programming

Chapter 8. Voronoi Diagrams. 8.1 Post Oce Problem

Convexity Theory and Gradient Methods

Lecture 2: August 31

Applied Lagrange Duality for Constrained Optimization

L-CONVEX-CONCAVE SETS IN REAL PROJECTIVE SPACE AND L-DUALITY

IE 5531: Engineering Optimization I

LECTURE 7 LECTURE OUTLINE. Review of hyperplane separation Nonvertical hyperplanes Convex conjugate functions Conjugacy theorem Examples

Search and Intersection. O Rourke, Chapter 7 de Berg et al., Chapter 11

AM 221: Advanced Optimization Spring 2016

ORIE 6300 Mathematical Programming I September 2, Lecture 3

Transcription:

Convex Sets Computer Science and Automation Indian Institute of Science Bangalore 560 012, India. NPTEL Course on

Let x 1, x 2 R n, x 1 x 2. Line and line segment Line passing through x 1 and x 2 : {y R n : y = λx 1 + (1 λ)x 2, λ R} Line Segment, LS[x 1, x 2 ]: {y R n : y = λx 1 + (1 λ)x 2, λ [0, 1]}

Affine sets Definition A set X R n is affine if for any x 1, x 2 X, and λ R, λx 1 + (1 λ)x 2 X. If X R n is an affine set, x 1, x 2,..., x k X and i λ i = 1, then the point λ 1 x 1 + λ 2 x 2 +... + λ k x k X. Examples Solution set of linear equations: {x R n : Ax = b} where A R m n, b R m. A subspace or a translated subspace

Result If X is an affine set and x 0 X, then {x x 0 : x X} forms a subspace. Proof. Let Y = {x x 0 : x X}. To show that Y is a subspace, we need to show, αy 1 + βy 2 Y for any y 1, y 2 Y and α, β R. Let y 1, y 2 Y and α, β R. Therefore, y 1 + x 0 X and y 2 + x 0 X. αy 1 + βy 2 = α(y 1 + x 0 ) + β(y 2 + x 0 ) (α + β)x 0. αy 1 + βy 2 + x 0 = α(y 1 + x 0 ) + β(y 2 + x 0 ) + (1 α β)x 0. Since X is an affine set, αy 1 + βy 2 + x 0 X. Thus, αy 1 + βy 2 Y Y is a subspace.

Definition Let X = {x 1, x 2,..., x k }. A point x is said to be an affine combination of points in X if x = k λ i x i, i=1 k λ i = 1. i=1 Definition Let X R n. The set of all affine combinations of points in X is called the affine hull of X. aff(x) = { k λ i x i : x 1,..., x k X, i i=1 λ i = 1}

Convex Sets Definition A set C R n is convex if for any x 1, x 2 C and any scalar λ with 0 λ 1, we have Examples Some simple convex sets λx 1 + (1 λ)x 2 C.

Some examples of nonconvex sets

Convex sets Examples: The empty set φ, any singleton set {x 0 } and R n are convex subsets of R n. A closed ball in R n, B[x 0, r] and an open ball in R n, B(x 0, r), are convex sets. Any affine set is a convex set. A line segment is a convex set, but not an affine set.

Possible to construct new convex sets using given convex sets Definition Let P R n, Q R n, α R. The scalar multiple αp of the set P is defined as αp = {x : x = αp, p P} The sum of two sets P and Q is the set, P + Q = {x : x = p + q, p P, q Q}.

Theorem If {C i }, i A, is any collection of convex sets, then i A C i is a convex set. Proof. Let C = i A C i. Case 1. If C is empty or singleton, then clearly C is a convex set. Case 2. To show that, for any x 1, x 2 C and λ [0, 1], λx 1 + (1 λ)x 2 C. Let x 1, x 2 C. Clearly, x 1, x 2 C i for every i A. Since every C i is convex, λx 1 + (1 λ)x 2 C i, i A, λ [0, 1] λx 1 + (1 λ)x 2 i A C i, λ [0, 1] λx 1 + (1 λ)x 2 C, λ [0, 1].

Theorem If C 1 and C 2 are convex sets, then C 1 + C 2 is a convex set. Proof. Let x 1, y 1 C 1 and x 2, y 2 C 2. So, x 1 + x 2, y 1 + y 2 C 1 + C 2. We need to show, λ(x 1 + x 2 ) + (1 λ)(y 1 + y 2 ) C 1 + C 2 λ [0, 1]. Since C 1 and C 2 are convex, z 1 = λx 1 + (1 λ)y 1 C 1 λ [0, 1], and z 2 = λx 2 + (1 λ)y 2 C 2 λ [0, 1]. Thus, z 1 + z 2 C 1 + C 2...... (1) Now, x 1 + x 2 C 1 + C 2 and y 1 + y 2 C 1 + C 2. Therefore, from (1), z 1 + z 2 = λ(x 1 + x 2 ) + (1 λ)(y 1 + y 2 ) C 1 + C 2 λ [0, 1]. Thus, C 1 + C 2 is a convex set.

Theorem If C is a convex set and α R, then αc is a convex set. Proof. Let x 1, x 2 C. Since C is convex, λx 1 + (1 λ)x 2 C λ [0, 1]. Also, αx 1, αx 2 αc. Therefore, λ(αx 1 ) + (1 λ)(αx 2 ) = α(λx 1 + (1 λ)x 2 ) αc for any λ [0, 1]. Therefore, αc is a convex set.

Hyperplane Definition Let b R and a R n, a 0. Then, the set is said to be a hyperplane in R n. H = {x : a T x = b} a denotes the normal to the hyperplane H. If a = 1, then b is the distance of H from the origin. In R 2, hyperplane is a line In R 3, hyperplane is a plane

Half-spaces The sets H + = {x : a T x b} and H = {x : a T x b} are called closed positive and negative half-spaces generated by H.

Some more examples of convex sets H = {x : a T x = b} is a convex set Let a 1, a 2,..., a m R n and b 1, b 2,..., b m R a T 1 b 1 a T 2 Define, A =. and b = b 2. a T m Then, {x : Ax = b} is a convex set. H + = {x : a T x b} and H = {x : a T x b} are convex sets. b m {x : Ax b} and {x : Ax b} are convex sets.

Convex Hull Definition The convex hull of a set S is the intersection of all the convex sets which contain S and is denoted by conv(s). Note: Convex Hull of a set S is the convex set. The smallest convex set that contains S is called the convex hull of the set S. Examples: conv({x, y}) = LS[x, y] where x and y are two points Let S = {(x, 0) : x R} {(0, y) : y R}. Then conv(s) = R 2.

S = {(x 1, x 2 ) : x 2 = x 2 1} conv(s) = {(x 1, x 2 ) : x 2 x 2 1}

Theorem CS1 Let S R n be a nonempty, closed convex set and y / S. Then there exists a unique point x 0 S with minimum distance from y. Further, x 0 is the minimizing point iff (y x 0 ) T (x x 0 ) 0 for all x S.

Proof Sketch Given S R n, a closed convex set and y / S. First show that there exists a minimizing point x 0 in S that is closest to y. Let δ = inf x S x y. Note that S is not bounded.

Proof Sketch (continued) Consider the set S B[y, 2δ]. By Weierstrass Theorem, there exists a minimizing point x 0 in S that is closest to y. x 0 = arg min x S B[y,2δ] x y

Proof Sketch (continued) Show uniqueness using triangle inequality Assume that there exists ˆx S such that y x 0 = y ˆx = δ. Since S is convex, 1 2 (x 0 + ˆx) S. Using triangle inequality, 2 y (x 0 + ˆx) y x 0 + y ˆx = 2δ 2 We have a contradiction if strict inequality holds.

Proof Sketch (continued) To show that x 0 is the unique minimizing point iff (y x 0 ) T (x x 0 ) 0 for all x S. Let x S. Assume (y x 0 ) T (x x 0 ) 0. y x 2 = y x 0 + x 0 x 2 = y x 0 2 + x 0 x 2 + 2(y x 0 ) T (x 0 x) Using the assumption, (y x 0 ) T (x 0 x) 0, we get y x 2 y x 0 2 This implies that x 0 is the minimizing point.

Proof Sketch (continued) Assume that x 0 is the minimizing point, that is, y x 0 2 y z 2 z S. Consider any x S. Since S is convex, λx + (1 λ)x 0 S λ [0, 1]. Therefore, y x 0 2 y x 0 λ(x x 0 ) 2 That is, y x 0 2 y x 0 2 +λ 2 x x 0 ) 2 2λ(y x 0 ) T (x x 0 ) 2(y x 0 ) T (x x 0 ) λ x x 0 ) 2 Letting λ 0 +, the result follows.

Theorem CS1 Let S R n be a nonempty, closed convex set and y / S. Then there exists a unique point x 0 S with minimum distance from y. Further, x 0 is the minimizing point iff (y x 0 ) T (x x 0 ) 0 for all x S.

Separating hyperplanes Definition Let S 1 and S 2 be nonempty subsets in R n and H = {x : a T x = b} be a hyperplane. 1 The hyperplane H is said to separate S 1 and S 2, if a T x b x S 1 and a T x b x S 2. 2 If a T x > b x S 1 and a T x < b x S 2, then the hyperplane H is said to strictly separate S 1 and S 2. 3 The hyperplane H is said to strongly separate S 1 and S 2 if a T x b + ɛ x S 1 and a T x b x S 2 where ɛ is a positive scalar.

H: Separating hyperplane The hyperplane H = {x : a T x = b} is said to separate S 1 and S 2, if a T x b x S 1 and a T x b x S 2.

H: Strictly separating hyperplane If a T x > b x S 1 and a T x < b x S 2, then the hyperplane H = {x : a T x = b} is said to strictly separate S 1 and S 2.

H: Strongly separating hyperplane The hyperplane H = {x : a T x = b} is said to strongly separate S 1 and S 2 if a T x b + ɛ x S 1 and a T x b x S 2 where ɛ is a positive scalar.

Separation of a closed convex set and a point Result Let S be a nonempty closed convex set in R n and y / S. Then there exists a nonzero vector a and a scalar b such that a T y > b and a T x b x S. Proof. By Theorem CS1, there exists a unique minimizing point x 0 S such that (x x 0 ) T (y x 0 ) 0 for each x S. Letting a = (y x 0 ) and b = a T x 0, we get a T x b for each x S and a T y b = (y x 0 ) T (y x 0 ) > 0 (since y x 0 ).

Cone Definition A set K R n is called a cone if for every x K and λ 0, we have λx K. K is a convex cone if it is convex and a cone.

Some examples of cones

Farkas Lemma Farkas Lemma Let A R m n and c R n. Then, exactly one of the following two systems has a solution: (I) Ax 0, c T x > 0 for some x R n (II) A T y = c, y 0 for some y R m.

Let a 1, a 2,..., a m R n. a T 1 a T 2 Define, A =.. a T m For example,

Farkas Lemma (Geometrical Interpretation) Exactly one of the following two systems has a solution: (I) Ax 0, c T x > 0 for some x R n (II) A T y = c, y 0 for some y R m

Farkas Lemma Let A R m n and c R n. Then, exactly one of the following two systems has a solution: (I) Ax 0, c T x > 0 (II) A T y = c, y 0 for some x R n for some y R m Proof of Farkas Lemma (a) Suppose, system (II) has a solution. Therefore, y 0 such that c = A T y. Consider x R n such that Ax 0. So, c T x = y T Ax 0. That is, System I has no solution.

Farkas Lemma Let A R m n and c R n. Then, exactly one of the following two systems has a solution: (I) Ax 0, c T x > 0 (II) A T y = c, y 0 for some x R n for some y R m Proof of Farkas Lemma (continued) (b) Suppose, system (II) has no solution. Let S = {x : x = A T y, y 0}, a closed convex set and c / S. Therefore, p R n and α R such that p T x α x S and c T p > α. This means, α 0 (Since 0 S) c T p > 0. Also, α p T A T y = y T Ap y 0. Since y 0, Ap 0 (as y can be made arbitrarily large). Thus, p R n such that Ap 0, c T p > 0 System (I) has a solution.

Farkas Lemma Let A R m n and c R n. Then, exactly one of the following two systems has a solution: (I) Ax 0, c T x > 0 for some x R n (II) A T y = c, y 0 for some y R m. Corollary Let A R m n. Then exactly one of the following systems has a solution: (I) Ax < 0 for some x R n (II) A T y = 0, y 0 for some nonzero y R m.

Corollary Let A R m n. Then exactly one of the following systems has a solution: (I) Ax < 0 for some x R n (II) A T y = 0, y 0 for some nonzero y R m. Proof. We can write system I as Ax + ze 0 for some x R n, z > 0 where e is a m-dimensional vector containing all 1 s. That is, ( ) ( ) x x (A e) 0, (0,..., 0, 1) > 0 z z for some (x z) T R n+1.

Proof. (continued) Recall Farkas Lemma. Therefore, System II is ( ) A T y = (0,..., 0, 1) T, y 0 for some y R m. e T That is A T y = 0, e T y = 1 and y 0 for some y R m. Using Farkas Lemma, the result follows.

Supporting Hyperplanes of sets at boundary points Definition Let S R n, S φ. Let x 0 be a boundary point of S. A hyperplane H = {x : a T (x x 0 ) = 0} is called a supporting hyperplane of S at x 0, if either S H + (that is, a T (x x 0 ) 0 x S), or S H (that is, a T (x x 0 ) 0 x S).

Theorem Let S be a nonempty convex set in R n and x 0 be a boundary point of S. Then, there exists a hyperplane that supports S at x 0.

Theorem Let S 1 and S 2 be two nonempty disjoint convex sets in R n. Then, there exists a hyperplane that separates S 1 and S 2. Proof. Consider the set S = S 1 S 2 = {x 1 x 2 : x 1 S 1, x 2 S 2 }. Note that S is a convex set and 0 / S. Thus, we have a convex set and a point not in that set. So, there exists a vector a such that a T x 0 x S. In other words, a T x 1 a T x 2 x 1 S 1, x 2 S 2.