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.

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

Convex Optimization. Optimality conditions. (EE227BT: UC Berkeley) Lecture 9 (Optimality; Conic duality) 9/25/14. Laurent El Ghaoui.

Solutions to exam in SF1811 Optimization, Jan 14, 2015

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

APPENDIX A Some Linear Algebra

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

Lagrange Multipliers Kernel Trick

Lecture 10 Support Vector Machines II

princeton univ. F 17 cos 521: Advanced Algorithm Design Lecture 7: LP Duality Lecturer: Matt Weinberg

Yong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 )

6.854J / J Advanced Algorithms Fall 2008

Affine transformations and convexity

Assortment Optimization under MNL

p(z) = 1 a e z/a 1(z 0) yi a i x (1/a) exp y i a i x a i=1 n i=1 (y i a i x) inf 1 (y Ax) inf Ax y (1 ν) y if A (1 ν) = 0 otherwise

MMA and GCMMA two methods for nonlinear optimization

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

Appendix for Causal Interaction in Factorial Experiments: Application to Conjoint Analysis

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

College of Computer & Information Science Fall 2009 Northeastern University 20 October 2009

find (x): given element x, return the canonical element of the set containing x;

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

Errors for Linear Systems

Linear, affine, and convex sets and hulls In the sequel, unless otherwise specified, X will denote a real vector space.

Lecture 10 Support Vector Machines. Oct

Proseminar Optimierung II. Victor A. Kovtunenko SS 2012/2013: LV

1 Convex Optimization

PHYS 705: Classical Mechanics. Calculus of Variations II

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Perron Vectors of an Irreducible Nonnegative Interval Matrix

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

Research Article Global Sufficient Optimality Conditions for a Special Cubic Minimization Problem

Appendix B. Criterion of Riemann-Stieltjes Integrability

On a direct solver for linear least squares problems

Non-negative Matrices and Distributed Control

Lecture Notes on Linear Regression

Report on Image warping

Example: (13320, 22140) =? Solution #1: The divisors of are 1, 2, 3, 4, 5, 6, 9, 10, 12, 15, 18, 20, 27, 30, 36, 41,

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2

On the Global Linear Convergence of the ADMM with Multi-Block Variables

THE WEIGHTED WEAK TYPE INEQUALITY FOR THE STRONG MAXIMAL FUNCTION

Norms, Condition Numbers, Eigenvalues and Eigenvectors

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

8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS

REAL ANALYSIS I HOMEWORK 1

General viscosity iterative method for a sequence of quasi-nonexpansive mappings

MATH 241B FUNCTIONAL ANALYSIS - NOTES EXAMPLES OF C ALGEBRAS

BOUNDEDNESS OF THE RIESZ TRANSFORM WITH MATRIX A 2 WEIGHTS

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family

The Order Relation and Trace Inequalities for. Hermitian Operators

ADVANCED MACHINE LEARNING ADVANCED MACHINE LEARNING

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

Dynamic Systems on Graphs

Lecture 12: Discrete Laplacian

First day August 1, Problems and Solutions

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

Lecture 20: November 7

Support Vector Machines

Lecture 6: Support Vector Machines

Marginal Effects in Probit Models: Interpretation and Testing. 1. Interpreting Probit Coefficients

Kernel Methods and SVMs Extension

The proximal average for saddle functions and its symmetry properties with respect to partial and saddle conjugacy

A SURVEY OF PROPERTIES OF FINITE HORIZON DIFFERENTIAL GAMES UNDER ISAACS CONDITION. Contents

Notes on Analytical Dynamics

The Second Anti-Mathima on Game Theory

STEINHAUS PROPERTY IN BANACH LATTICES

Infinitely Split Nash Equilibrium Problems in Repeated Games

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

Support Vector Machines

Modelli Clamfim Equazioni differenziali 7 ottobre 2013

Economics 101. Lecture 4 - Equilibrium and Efficiency

Convergence rates of proximal gradient methods via the convex conjugate

Canonical transformations

Mixed Taxation and Production Efficiency

Bernoulli Numbers and Polynomials

Module 9. Lecture 6. Duality in Assignment Problems

Maximal Margin Classifier

MATH Homework #2

COS 521: Advanced Algorithms Game Theory and Linear Programming

Problem Set 9 Solutions

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION

Inexact Alternating Minimization Algorithm for Distributed Optimization with an Application to Distributed MPC

TAIL BOUNDS FOR SUMS OF GEOMETRIC AND EXPONENTIAL VARIABLES

Convexity preserving interpolation by splines of arbitrary degree

Lecture 17: Lee-Sidford Barrier

Game Theory. Lecture Notes By Y. Narahari. Department of Computer Science and Automation Indian Institute of Science Bangalore, India February 2008

Approximations of Set-Valued Functions Based on the Metric Average

Projective change between two Special (α, β)- Finsler Metrics

Topic 5: Non-Linear Regression

Math 217 Fall 2013 Homework 2 Solutions

Some modelling aspects for the Matlab implementation of MMA

ECE 534: Elements of Information Theory. Solutions to Midterm Exam (Spring 2006)

Lecture 3. Ax x i a i. i i

Perfect Competition and the Nash Bargaining Solution

Generalized Linear Methods

FIRST AND SECOND ORDER NECESSARY OPTIMALITY CONDITIONS FOR DISCRETE OPTIMAL CONTROL PROBLEMS

NP-Completeness : Proofs

Control of Uncertain Bilinear Systems using Linear Controllers: Stability Region Estimation and Controller Design

The Minimum Universal Cost Flow in an Infeasible Flow Network

ON THE EXTENDED HAAGERUP TENSOR PRODUCT IN OPERATOR SPACES. 1. Introduction

for Linear Systems With Strictly Diagonally Dominant Matrix

Transcription:

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 problem s fx : gx 0, x E }. Lx; λ = fx λ T gx, and the prmal and dual values d are defned, respectvely, as We have to show that p = sup Lx; λ, λ R m d = sup Lx; λ. λ R m p d. Wthout loss of generalty, we may assume that p < else ths nequalty s trval. Then the defnton of p mples that for every q > p there exsts some x E wth sup Lx; λ < q. λ R m In other words: For every q > p there exsts some x q E such that Lx q ; λ < q for all λ R m. d = sup Lx; λ sup Lx q ; λ < q. λ R m λ R m d < q for all q > p, whch mples that d p. 1 3. Slater and compactness Prove the Slater condton holds for problem 4.3.1 f and only f there s a pont ˆx E for whch the level sets λ R m : Lˆx; λ α } are compact. Note that the level sets of Lˆx; are always closed because L s contnuous n λ. the compactness of the level sets s equvalent to ther boundedness. Now recall that the Slater condton stated that there exsts ˆx dom f wth g ˆx < 0 for all. Assume that ths condton holds, and let α R. Let moreover λ R m satsfy Lˆx; λ α. Then, usng the fact that λ 0 for all and g ˆx < 0 for all, we obtan α fˆx λ g ˆx λ max g ˆx. Date: November 3, 015. 1

SELECTED SOLUTIONS, SECTION 4.3 As a consequence, λ α fˆx max g ˆx. the set gven n 1 s bounded. Conversely, assume that there exsts ˆx E such that the set gven n 1 s bounded for every α. Then, for each, the set of λ 0 such that fˆx λg ˆx α s bounded, whch mples that fˆx > ˆx dom f. Replacng now α by ˆα := α fˆx, we obtan the set of λ 0 such that λ g ˆx ˆα s bounded for all. Ths, however, s only possble f g ˆx < 0 for all, Slater s condton holds. 4. Examples of duals Lots and lots of them. The followng examples are orgnally stated wth condtons of the form a, x b for = 1,,..., m. I have rewrtten these nequaltes as Ax b, where A R n m denotes the matrx wth rows equal to the transposes of the vectors a. Also, the computatons soon become a bt repettve. a The lnear program We have b Another lnear program } c, x : Ax b. Lx; λ = c, x λ, Ax b Φλ = x Rnc, x λ, Ax b = λ, b c A λ, x f c A λ 0, = λ, b f c A λ = 0. Lx, λ = c, x δr n x : Ax b }. c, x λ, Ax b f x 0, f x 0. Φλ = c, x λ, Ax b x 0 = λ, b c x 0 A λ, x λ, b f c A λ 0, = f c A λ 0. c The quadratc program for C n } : Ax b. Lx; λ = xt Cx λ, Ax b.

SELECTED SOLUTIONS, SECTION 4.3 3 Φλ = λ, Ax b = λ, b A λ, x. Snce C S, n the mnmzaton problem to be solved here s convex and the mnmum s attaned at the pont where the gradent wth respect to x s zero. That s, at x satsfyng or and Cx A λ = 0, x = C 1 A λ. A λ, x = λt AC 1 A λ Φλ = λ, b λt AC 1 A λ. d The separable problem n } px j : Ax b. Lx; λ = Φλ = px j λ, Ax b = b, λ = b, λ = b, λ = b, λ = b, λ = b, λ px j λ, Ax b px j A λ, x pxj A λ j x j pt A λ j t t R sup t R A λ j t pt p A λ j p a j, λ, where a j denotes the j-th columns of the matrx A. e The penalzed lnear program } c, x ε lbx : Ax b.

4 SELECTED SOLUTIONS, SECTION 4.3 Recall that lbx = log x for x > 0. Φλ = c, x ε lbx λ, Ax b x>0 = λ, b c x>0 A λ, x ε lbx c A λ/ε, x lbx = λ, b ε sup x>0 = λ, b ε lb c A λ/ε Proposton 3.3.3 computaton of lb mples that λ, b ε lb c A λ/ε εn f c A λ > 0, Φλ = f c A λ 0. f The semdefnte program C, X δs n X : A, X b for = 1,,..., m }. Φλ = C, X = λ, b λ A, X b =1 C λ A, X. Now recall that S n s a self-dual cone, that s, N S n 0 = S n see Exercse.1.3., f C λ A S n, then C λ A, X 0 for all X S n C λ A, X = 0. On the other hand, f C λ A S n, then there exsts ˆX S n such that C λ A, ˆX > 0. C λ A, X C t>0 λ A, t ˆX =. Φλ = λ, b f C λ A S n, f C λ A S n.

SELECTED SOLUTIONS, SECTION 4.3 5 g The penalzed semdefnte program C, X ε ldx : A, X b for = 1,,..., m }. Φλ = C, X ε ldx = λ, b = λ, b ε sup = λ, b ε ld C C C λ A, X b λ A, X ε ldx λ A /ε, X ldx λ A /ε. Now Proposton 3.3.3 mples that Φλ = λ, b ε ld C λ A /ε εn. 9. Fenchel and Lagrangan dualty Let Y be a Eucldean space. By sutably rewrtng the prmal Fenchel problem } fx gax for gven functons f : E, ], g : Y, ], and lnear A: E Y, nterpret the dual Fenchel problem f A φ g φ } sup φ Y as a Lagrangan dual problem. We rewrte the problem as the constrant optmzaton problem fx gy} subject to ± Ax y 0. The Lagrangan of ths problem s Lx, y; λ = fx gy λ, Ax y µ, y Ax. Now the Lagrangan dual problem s sup fx gy λ µ, Ax y λ, µ Y = sup λ, µ Y fx µ λ, Ax gy λ µ, y = sup sup A µ λ, x fx sup λ µ, y gy λ, µ Y = sup f A µ λ g λ µ λ, µ Y = sup φ Y f A φ g φ.