Lecture 9 Monotone VIs/CPs Properties of cones and some existence results. October 6, 2008

Similar documents
Lecture 8 Plus properties, merit functions and gap functions. September 28, 2008

Lecture 7 Monotonicity. September 21, 2008

Course Notes for EE227C (Spring 2018): Convex Optimization and Approximation

Lecture 1. 1 Conic programming. MA 796S: Convex Optimization and Interior Point Methods October 8, Consider the conic program. min.

Lecture 8. Strong Duality Results. September 22, 2008

Lecture 5. Theorems of Alternatives and Self-Dual Embedding

Convex Sets. Prof. Dan A. Simovici UMB

GENERALIZED CONVEXITY AND OPTIMALITY CONDITIONS IN SCALAR AND VECTOR OPTIMIZATION

SOME STABILITY RESULTS FOR THE SEMI-AFFINE VARIATIONAL INEQUALITY PROBLEM. 1. Introduction

1 Review of last lecture and introduction

Chapter 1. Preliminaries

Convex Optimization Theory. Chapter 5 Exercises and Solutions: Extended Version

Game theory: Models, Algorithms and Applications Lecture 4 Part II Geometry of the LCP. September 10, 2008

MAT-INF4110/MAT-INF9110 Mathematical optimization

Some Properties of Convex Hulls of Integer Points Contained in General Convex Sets

Summer School: Semidefinite Optimization

Optimization and Optimal Control in Banach Spaces

Week 3: Faces of convex sets

On duality gap in linear conic problems

LECTURE 25: REVIEW/EPILOGUE LECTURE OUTLINE

Optimality Conditions for Nonsmooth Convex Optimization

Supplement: Hoffman s Error Bounds

Optimality Conditions for Constrained Optimization

Assignment 1: From the Definition of Convexity to Helley Theorem

Some Properties of the Augmented Lagrangian in Cone Constrained Optimization

c 2000 Society for Industrial and Applied Mathematics

Lecture: Cone programming. Approximating the Lorentz cone.

Lecture 13 Newton-type Methods A Newton Method for VIs. October 20, 2008

Convex Optimization M2

On the projection onto a finitely generated cone

On duality theory of conic linear problems

Approximate Farkas Lemmas in Convex Optimization

Mathematics 530. Practice Problems. n + 1 }

Appendix A: Separation theorems in IR n

Constraint qualifications for convex inequality systems with applications in constrained optimization

Appendix PRELIMINARIES 1. THEOREMS OF ALTERNATIVES FOR SYSTEMS OF LINEAR CONSTRAINTS

LECTURE 4 LECTURE OUTLINE

GEORGIA INSTITUTE OF TECHNOLOGY H. MILTON STEWART SCHOOL OF INDUSTRIAL AND SYSTEMS ENGINEERING LECTURE NOTES OPTIMIZATION III

LECTURE 9 LECTURE OUTLINE. Min Common/Max Crossing for Min-Max

Asteroide Santana, Santanu S. Dey. December 4, School of Industrial and Systems Engineering, Georgia Institute of Technology

Convex Analysis and Optimization Chapter 2 Solutions

Lecture 5. The Dual Cone and Dual Problem

ON A CLASS OF NONSMOOTH COMPOSITE FUNCTIONS

Extreme Abridgment of Boyd and Vandenberghe s Convex Optimization

Key words. Complementarity set, Lyapunov rank, Bishop-Phelps cone, Irreducible cone

Copositive Plus Matrices

CONSTRAINT QUALIFICATIONS, LAGRANGIAN DUALITY & SADDLE POINT OPTIMALITY CONDITIONS

Locally convex spaces, the hyperplane separation theorem, and the Krein-Milman theorem

Course 212: Academic Year Section 1: Metric Spaces

Maximal Monotone Inclusions and Fitzpatrick Functions

1 Maximal Lattice-free Convex Sets

Normal Fans of Polyhedral Convex Sets

E5295/5B5749 Convex optimization with engineering applications. Lecture 5. Convex programming and semidefinite programming

Convex Analysis and Optimization Chapter 4 Solutions

Boundary Behavior of Excess Demand Functions without the Strong Monotonicity Assumption

Research Division. Computer and Automation Institute, Hungarian Academy of Sciences. H-1518 Budapest, P.O.Box 63. Ujvári, M. WP August, 2007

Polynomial complementarity problems

Inequality Constraints

Lecture 19 Algorithms for VIs KKT Conditions-based Ideas. November 16, 2008

5. Duality. Lagrangian

UTILITY OPTIMIZATION IN A FINITE SCENARIO SETTING

IE 521 Convex Optimization Homework #1 Solution

Barrier Method. Javier Peña Convex Optimization /36-725

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

Lecture 3. Optimization Problems and Iterative Algorithms

Integer Programming, Part 1

A Geometric Framework for Nonconvex Optimization Duality using Augmented Lagrangian Functions

Convex Optimization & Lagrange Duality

CONVERGENCE PROPERTIES OF COMBINED RELAXATION METHODS

Chapter 2: Preliminaries and elements of convex analysis

Implications of the Constant Rank Constraint Qualification

Lagrange Duality. Daniel P. Palomar. Hong Kong University of Science and Technology (HKUST)

Math 350 Fall 2011 Notes about inner product spaces. In this notes we state and prove some important properties of inner product spaces.

Analysis III Theorems, Propositions & Lemmas... Oh My!

Applied Mathematics Letters

Introduction to Convex and Quasiconvex Analysis

A projection algorithm for strictly monotone linear complementarity problems.

Two-Term Disjunctions on the Second-Order Cone

3. Linear Programming and Polyhedral Combinatorics

Real Analysis Math 131AH Rudin, Chapter #1. Dominique Abdi

Lecture 1: Basic Concepts

Conic Linear Optimization and its Dual. yyye

A Parametric Simplex Algorithm for Linear Vector Optimization Problems

Convex Geometry. Carsten Schütt

CHAPTER 2: CONVEX SETS AND CONCAVE FUNCTIONS. W. Erwin Diewert January 31, 2008.

Lecture: Duality.

Mathematical Optimisation, Chpt 2: Linear Equations and inequalities

UNDERGROUND LECTURE NOTES 1: Optimality Conditions for Constrained Optimization Problems

Optimality, Duality, Complementarity for Constrained Optimization

Linear Algebra Review: Linear Independence. IE418 Integer Programming. Linear Algebra Review: Subspaces. Linear Algebra Review: Affine Independence

POLARS AND DUAL CONES

Introduction to Mathematical Programming IE406. Lecture 3. Dr. Ted Ralphs

Convex Sets with Applications to Economics

On John type ellipsoids

Chap 2. Optimality conditions

NORMS ON SPACE OF MATRICES

LP Duality: outline. Duality theory for Linear Programming. alternatives. optimization I Idea: polyhedra

Closedness of Integer Hulls of Simple Conic Sets

Primal/Dual Decomposition Methods

Extended Monotropic Programming and Duality 1

Transcription:

Lecture 9 Monotone VIs/CPs Properties of cones and some existence results October 6, 2008

Outline Properties of cones Existence results for monotone CPs/VIs Polyhedrality of solution sets Game theory: Models, Algorithms and Applications 1

Properties of cones Motivation: Recall that in the earlier section, existence statements required that the range of x ref satisfied a certain interioricity property: namely if F (x ref ) int((k ) ), then SOL(K.F ) is nonempty and compact. Essentially, this reduces to characterizing the interior of the dual cone Given an arbitrary closed convex cone, we provide results for a vector to be in the relative interior of C. Game theory: Models, Algorithms and Applications 2

Let C be a closed convex cone in R n Introduction The set C ( C) is a linear subspace in R n called the lineality space of C; This is the largest subspace contained in C and is denoted by lin C = C ( C). For S 1, S 2 R n, we have (S 1 + S 2 ) = S 1 S 2. If C 1 and C 2 are closed convex cones such that the sum C 1 + C 2 is closed then we have (C 1 C 2 ) = C 1 + C 2. We may then show that the linear hull of C denoted by linc satisfies linc = (C ( C)), where A is the orthogonal complement of A. The set v is a subspace containing vectors that are orthogonal to v Game theory: Models, Algorithms and Applications 3

Relative interior Definition 1 Let S be a subset of R n. The relative interior of S is the interior of S considered as a subset of the affine hull Aff(S) and is denoted by ri(s). The formal definition is ri(s) = {x S : ɛ > 0, (N ɛ aff(s)) S}. Consider S = {(x, y, 0) : 0 x, y 1} R 3. Its interior int(s) =. Its relative interior is nonempty since ri(s) = {(x, y, 0) : 0 < x, y < 1}. This follows from noting that aff(s) = {(x, y, 0) : x, y R} Consider S = {(x, y, z) : 0 x, y, z 1} R 3. Game theory: Models, Algorithms and Applications 4

Its interior and relative interior are the same and are given by int(s) = ri(s) = {(x, y, z) : 0 < x, y, z < 1}. If S = ri(s), S is said to be relatively open The relative boundary, rbd(s) = S ri(s). Prop: Let C be a closed convex cone in R n. Then the following are equivalent: (a) v ri(c ) (b) For all x 0, x C lin C, v T x > 0 (c) For every scalar η > 0, the set S(v, η) {x C lin C : v T x η} is bounded (d) v C and the intersection C v is a linear subspace If any of the above holds, then C v is equal to the lineality space of C. Game theory: Models, Algorithms and Applications 5

Solid and pointed cones Definition 2 A cone C is pointed if C ( C) = {0}. A set S is solid if int (S). R n + is pointed and solid pos(a) is also pointed and solid Lemma 1 Let C be a closed convex cone in R n. Then the following hold: 1. If C is solid, then its dual C is pointed. 2. If C is pointed, then its dual C is solid. Thus C is solid(pointed) if and only if C is pointed(solid). Proof: Game theory: Models, Algorithms and Applications 6

1. If C is solid, then to show that C is pointed, we need to show that C ( C ) is {0}. Let d (C ( C )). Therefore, we have d T x = 0 for all x C. Since C is closed and convex, we have C = C. Moreover, from C being solid, we have int (C). Let y be an interior point of C = C. Then, by definition y T d > 0 if d 0. But this contradicts d T x = 0, x C. Therefore, d = 0 and C is pointed. 2. Let C be pointed. To show that C has a nonempty interior, let k be the maximum number of linearly independent vectors in C and let {y 1,..., y k } represent a collection of such vectors: (a) If k = n, then C contains the simplicial (pos) cone generated by these vectors and is nonempty. Since C contains this cone, the result follows. (b) If k < n, then the system p T y i = 0, i = 1,..., k for some 0 p R n. But these k vectors span the space in C implying that Game theory: Models, Algorithms and Applications 7

p T y = 0, y C. Consequently, p C ( C ) or p C ( C) since C = C. Therefore C ( C) {0}, contradicting the pointedness of C. (c) To establish the last assertion, assume that C is pointed. Then by (b), we have that C is solid. But since C is closed and convex, then C = C. Thus C is solid. C is solid implies C is pointed in a similar fashion. Game theory: Models, Algorithms and Applications 8

Existence results Theorem 1 Let K be a closed convex cone in R n and let F be a pseudomonotone continuous map from K into R n. Then the following three statements are equivalent: (a) The CP(K,F) is strictly feasible. (b) The dual cone K has a nonempty interior and SOL(K,F) is a nonempty compact set. (c) The dual cone K has a nonempty interior and K [ (F (K) )] = {0}. Proof: Game theory: Models, Algorithms and Applications 9

(a) = (b): If the CP(K,F) is strictly feasible, then clearly int K is nonempty. From theorem 2.3.5, we have that SOL(K,F) is nonempty and compact. Theorem 2 Suppose F is pseudo-monotone on K. If there exists an x ref K satisfying F (x ref ) (int(k ) ), then SOL(K, F ) is nonempty, convex and compact. (b) = (c): This follows from the next result (provided without proof) and by noting that K = K : Theorem 3 Let K be a closed and convex set and F be a pseudomonotone continuous mapping. Then SOL(K,F) is nonempty and bounded if and only if K [ (F (K) )] = {0}. Game theory: Models, Algorithms and Applications 10

(c) = (a): By prop 2.3.17, deg(f nat K, Ω) = 1 for every bounded open set Ω containing SOL(K,F). Let q be an arbitrary vector in int K. For a fixed but arbitrary scalar ɛ > 0, we have (F ɛ,nat K (x) = x Π K (x F (x) + ɛq) be the natural map of the perturbed CP: K x F (x) ɛq K. For sufficiently small ɛ we have deg(f ɛ,nat K, Ω) = deg(fnat, Ω) = 1. K Thus the perturbed CP has a solution, say x (from degree-theoretic results proved earlier). Since q intk, x must belong to K and F (x) to K. Hence, the CP(K,F) is strictly feasible. Game theory: Models, Algorithms and Applications 11

Special case: an affine map Corollary 4 Let K be a closed convex cone in R n and let F (x) = q + Mx be an affine pseudo-monotone map from K into R n. Then the following three statements are equivalent: 1. There exists a vector x K such that Mx + q is in K. 2. The dual cone K has a nonempty interior and SOL(K,F) is a nonempty compact set. 3. The dual cone K has a nonempty interior and [d K, M T d ( K), q T d 0] = d = 0. Game theory: Models, Algorithms and Applications 12

It suffices to show that K [ (F (K) )] = {0} is equivalent to the assertion in (c) above. In fact the former may be stated as [d K, M T d ( K), q T d 0] =[d K, (q) T d 0, (Mx) T d 0 x K] =[d K, d T (q + Mx) 0, x K] =K ( F (K)), F (K) = Mx + q, ={0} = d = 0, giving us the required result. Game theory: Models, Algorithms and Applications 13

Example 1 Consider the NCP(F): Feasibility = Solvability? - No!! F (x) ( ) 2x1 x 2 2x 2 + 1 x 2, (x 1, x 2 ) R 2. 1 + 2x 1 1 We have F = ( 2x2 2x 1 2 2 2x 1 0 ) 0 implying that F is monotone on R 2 +. The feasible region F EA(F ) = {(x 1, x 2 ) R 2 : x 1 = 1, x 2 0}. However, this NCP has no solution. Game theory: Models, Algorithms and Applications 14

However it does have ɛ exact solutions. For a scalar ɛ (0, 1), the vector x ɛ = (1 ɛ, 1/(2ɛ)). It may be verified that and there exists an x satisfying lim min(x ɛ, F (x ɛ )) = 0 ɛ 0 min( x, F ( x)) ɛ for every ɛ > 0. NCP(F). Such an x is called an ɛ approximate solution of This leads to two important questions: When does feasibility of CP(K,F) imply its solvability? K is polyhedral and F is affine Game theory: Models, Algorithms and Applications 15

F is strictly monotone on K and K is pointed but non-polyhedral Does every feasible monotone CP have ɛ solutions? - yes! Theorem 5 Let K be R n + and F be a monotone affine map from Rn onto itself. Then the CP(K,F) is solvable if and only if it is feasible. Proof: Since K is a polyhedral cone, we have that K = K. Moreover, let F (x) = q + Mx, where M 0. Consider the gap program given by min x T (Mx + q) subject to x K Mx + q K. This is a convex QP with an objective that is bounded below by zero on K. We then use the Frank-Wolfe theorem - which states that a quadratic Game theory: Models, Algorithms and Applications 16

function, bounded below on a polyhedral set, attains its minimum on that set. Therefore the dual gap program attains a minimum at x. By the necessary conditions of optimality, there exists a λ R n such that K x v q + (M + M T )x M T λ K K q + Mx λ K = K. Since (x ) T v = 0 and λ T v 0, we have 0 (x λ) T v = (x λ) T (q + (M + M T )x M T λ) = (x λ) T (q + Mx ) + (x λ) T M T (x λ) (x λ) T (q + Mx ) = (x ) T (Mx + q). Game theory: Models, Algorithms and Applications 17

But (x ) T (Mx + q) 0 implying that (x ) T (Mx + q) = 0 and x lies in SOL(K,F). Game theory: Models, Algorithms and Applications 18

General strictly monotone F Theorem 6 (Th. 2.4.8) Let K be a pointed closed and convex cone in R n and F : K R n be a coninuous map. Consider the following statements (a) F is strictly monotone on K and FEA(K,F) is nonempty; (b) F is strictly monotone on K and CP(K,F) is strictly feasible; (c) the CP(K,F) has a unique solution. It holds that (a) (b) (c). Proof: Clearly (b) = (a) (a) = (b) (omitted) - see Theorem 2.4.8 in FP-I Converse is proved in Th. 2.4.8 (b) = (c) follows from Th 2.3.5 Game theory: Models, Algorithms and Applications 19

F is affine monotone Under the assumptions of affineness and monotonicity of F, we obtain the following: Lemma 2 Let K be a convex cone in R n and F (x) = q + Mx where M 0. For any y SOL(K, q, M) it holds that SOL(K, q, M) = {x K : q+mx K, (M T +M)(x y) = 0, q T (x y) = 0.} Proof: By proposition 2.3.6, we have from monotonicity that x, y SOL(K, q, M): (x y) T M(x y) = 0. Game theory: Models, Algorithms and Applications 20

For x y, we have that (M + M T )(x y) = 0 which may be further simplified as x T (M + M T )(x y) = 0 x T (M + M T )x = x T (M + M T )y x T Mx = 1 2 xt (M + M T )y similarly y T My = 1 2 xt (M + M T )y Since x SOL(K, F ), we have that x T (q + Mx) = 0 = x T q = x T Mx. Similarly, y T (q + My) = 0 = y T q = y T My. But y T My = x T Mx = = x T q = x T Mx = y T My = y T q. Game theory: Models, Algorithms and Applications 21

Therefore x lies in the set on the right hand side of the specified set equation. To prove the reverse inclusion, we need to show that if x lies in {x K : q + Mx K, (M T + M)(x y) = 0, q T (x y) = 0} then x T (Mx + q) = 0. Since q T (x y) = 0, (M + M T )(x y) = 0, we have q T x = q T y, (M + M T )x = (M + M T )y. The latter expression implies that x T Mx = y T My implying that x ( Mx + q) = y T (My + y) = 0. The result follows. Game theory: Models, Algorithms and Applications 22

Implication: The aforementioned lemma implies that q T x is a constant scalar and (M + M T )x is a constant vector for all x SOL(K, q, M). Therefore if K is polyhedral, then the set given by SOL(K,F) is also polyhedral. Game theory: Models, Algorithms and Applications 23

Polyhedrality of Monotone AVI Consider the affine variational inequality denoted by AVI(K,q,M) which requires a vector x such that (y x) T (Mx + q) 0, y K, where K is a polyhedral set. Note that if F is non necessarily affine, we have a linearly constrained VI. Before proving a result pertaining to the polyhedrality of the SOL(K,q,M), we provide a result that gives a KKT characterization of the VI. Proposition 1 Let K be given by K {x R n : Ax b, Cx = d} Game theory: Models, Algorithms and Applications 24

for some matrices A R m n and B R l n. A vector x solves the VI(K,F) if and only if there exists λ R m and µ R l such that F (x) + C T µ + A T λ = 0 Cx d = 0 0 b Ax λ 0. Proof: For a fixed x and K polyhedral, we may formulate the following LP: min y T F (x) subject to y K, Game theory: Models, Algorithms and Applications 25

by noting that the VI is equivalent to finding an x such that y T F (x) x T F (x) y K. Therefore if x SOL(K, F ), then x is a solution to this LP. But by LP duality, we have the existence of (λ, µ) such that the following holds: F (x) + C T µ + A T λ = 0 Cx d = 0 0 b Ax λ 0. Reverse holds in a similar fashion. We may now prove the polyhedrality of the solution set of AVI(K,q,M). Game theory: Models, Algorithms and Applications 26

Clearly, when K is a cone, this follows from the result for monotone CP(K,q,M). Theorem 7 Let K be polyhedral and F (x) Mx + q. The solution set of AVI(K,q,M) is polyhedral. Proof: Let K be given as K {x R n : Ax b, Cx = d} for some matrices A R m n and B R l n. However, x is a solution to AVI(K,q,m) if and only if (x, µ, λ) is a solution to Mx + q + C T µ + A T λ = 0 Cx d = 0 0 b Ax λ 0. Game theory: Models, Algorithms and Applications 27

But this is an mixed LCP (a special case of a monotone affine CP) and has a polyhedral solution set. However, the solution set of AVI(K,q,M) is a projection of this solution set given by the mapping: x µ Rn+l+m x R n. λ Therefore the solution set of AVI(K,q,M) is also polyhedral. Game theory: Models, Algorithms and Applications 28