Essentials of Convex Optimization Max Welling Department of Computer Science University of Toronto ... Thus, the dual problem always provides lower bo...

0 downloads 0 Views 55KB Size

Loading...

Max Welling Department of Computer Science University of Toronto 10 King’s College Road Toronto, M5S 3G5 Canada [email protected]

Abstract This is a note to explain duality and convex optimization. It is based on Stephen Boyd’s book, chapter 5 (available online).

1 Lagrangians and all that Most kernel-based algorithms fall into two classes, either they use spectral techniques to solve the problem, or they use convex optimization techniques to solve the problem. Here we will discuss convex optimization. A constrained optimization problem can be expressed as follows, f0 (x) fi (x) ≤ 0 ∀i hj (x) = 0 ∀j

minimizex subject to

(1)

That is we have inequality constraints and equality constraints. We now write the primal Lagrangian of this problem, which will be helpful in the following development, X X λi fi (x) + νj hj (x) (2) LP (x, λ, ν) = f0 (x) + i

j

where we will assume in the following that λi ≥ 0 ∀i. From here we can define the dual problem by, LD (λ, ν) = infx LP (x, λ, ν) (3) This objective can actually become −∞ for certain values of its arguments. We will call parameters λ ≥ 0, ν for which LD > −∞ dual feasible. It is important to notice that the dual Lagrangian is a concave function of λ, ν because it is a pointwise infimum of a function. Hence, even if the primal is not convex, the dual is certainly concave! It is not hard to show that ∗

LD (λ, ν) ≤ p∗

where p is the primal optimal point. This simply follows because P ∗ ν h j j (x) ≤ 0 for a primal feasible point x . j

P i

(4) λi fi (x) +

Thus, the dual problem always provides lower bound to the primal problem. The optimal lower bound can be found by solving the dual problem, maximizeλ,ν subject to

LD (λ, ν) λi ≥ 0 ∀i

(5) ∗

which is therefore a convex optimization problem. If we call d the dual optimal point we always have: d∗ ≤ p∗ , which is called weak duality. p∗ − d∗ is called the duality gap. Strong duality holds when p∗ = d∗ . Strong duality is very nice, in particular if we can express the primal solution x∗ in terms of the dual solution λ∗ , ν ∗ , because then we can simply solve the dual problem and convert to the answer to the primal domain since we know that solution must then be optimal. Often the dual problem is easier to solve. So when does strong duality holds. Up to some mathematical details the answer is: if the primal problem is convex and the equality constraints are linear. This means that f0 (x) and {fi (x)} are convex functions and hj (x) = Ax − b. The primal problem can be written as follows, p∗ = inf sup LP (x, λ, ν) x λ≥0,ν

(6)

This can be seen as follows by noting that supλ≥0,ν LP (x, λ, ν) = f0 (x) when x is feasible but ∞ otherwise. To see this first check that by violating one of the constraints you can find a choice of λ, ν that makes the Lagrangian infinity. Also, when all the constraints are satisfied, the best we can do is maximize the additional terms to be zero, which is always possible. The dual problem by definition is given by, d∗ = sup inf LP (x, λ, ν) λ≥0,ν x

(7)

Hence, the “sup” and “inf” can be interchanged if strong duality holds, hence the optimal solution is a saddle-point. It is important to realize that the order of maximization and minimization matters for arbitrary functions (but not for convex functions). Try to imagine a “V” shapes valley which runs diagonally across the coordinate system. If we first maximize over one direction, keeping the other direction fixed, and then minimize the result we end up with the lowest point on the rim. If we reverse the order we end up with the highest point in the valley. There are a number of important necessary conditions that hold for problems with zero duality gap. These Karush-Kuhn-Tucker conditions turn out to be sufficient for convex optimization problems. They are given by, X X ∇f0 (x∗ ) + λ∗i ∇fi (x∗ ) + νj∗ ∇hj (x∗ ) = 0 (8) i

fi (x∗ ) ≤ 0 hj (x∗ ) = 0 λ∗i ≥ 0 λ∗i fi (x∗ ) = 0

j

(9) (10) (11) (12)

The first equation is easily derived because we already saw that p∗ = inf x LP (x, λ∗ , ν ∗ ) and hence all the derivatives must vanish. This condition has a nice interpretation as a “balancing of forces”. Imagine a ball rolling down a surface defined by f0 (x) (i.e. you are doing gradient descent to find the minimum). The ball gets blocked by a wall, which is the constraint. If the surface and constraint is convex then if the ball doesn’t move we have reached the optimal solution. At that point, the forces on the ball must balance. The

first term represent the force of the ball against the wall due to gravity (the ball is still on a slope). The second term represents the reaction force of the wall in the opposite direction. The λ represents the magnitude of the reaction force, which needs to be higher if the surface slopes more. We say that this constraint is “active”. Other constraints which do not exert a force are “inactive” and have λ = 0. The latter statement can be read of from the last KKT condition which we call “complementary slackness”. It says that either fi (x) = 0 (the constraint is saturated and hence active) in which case λ is free to take on a non-zero value. However, if the constraint is inactive: fi (x) ≤ 0, then λ must vanish. As we will see soon, the active constraints will correspond to the support vectors in SVMs! Complementary slackness is easily derived by, f0 (x∗ ) = LD (λ∗ , ν ∗ ) = inf f0 (x) +

X

x

∗

≤ f0 (x ) +

X i

≤ f0 (x∗ )

λ∗i fi (x∗ )

+

X

i ∗

νj hj (x )

λ∗i fi (x) +

X

νj hj (x)

j

(13)

j

(14)

where the first line follows from the definition the second because the inf is always smaller than any x∗ and the last because fi (x∗ ) ≤ 0, λ∗i ≥ 0 and hj (x∗ ) = 0. Hence all inequalities are equalities and each term is negative so each term must vanish separately.