Part 7 - Optimality Conditions and Karush-Kuhn-Tucker Conditions II

Introduction

We have previously discussed Abadie’s constraint qualifications and the Karush-Kuhn-Tucker (KKT) conditions. However, Abadie’s constraint qualifications does not hold always (of course :)). Further, checking Abadie’s constraint qualifications requires us to compute the tangent cone TS(x)T_S(\mathbf{x}^*), which is not always easy.

Example 1

Let,

S={g1(x)=x13+x20g2(x)=x15x20g3(x)=x20S = \begin{cases} g_1(\mathbf{x}) = -x_1^3 + x_2 \leq 0 \newline g_2(\mathbf{x}) = x_1^5 - x_2 \leq 0 \newline g_3(\mathbf{x}) = -x_2 \leq 0 \end{cases}

Further, let xˉ=(00)T\bar{\mathbf{x}} = \begin{pmatrix} 0 & 0 \end{pmatrix}^T, compute TS(xˉ)T_S(\bar{\mathbf{x}}) and G(xˉ)G(\bar{\mathbf{x}}). Feasible region and constraints Let x1k=1kx_1^k = \frac{1}{k}, by g1(x)g_1(\mathbf{x}), x2k1k3x_2^k \leq \frac{1}{k^3}.

Further, let x2kx_2^k be any such sequence that is also 0,k\geq 0, \forall k (since g3(x)g_3(\mathbf{x})).

Then,

limkxk=xˉλk(xkxˉ)=[λkkλkx2k]\begin{align*} \lim_{k \to \infty} x^k = \bar{\mathbf{x}} \newline \lambda^k (\mathbf{x}^k - \bar{\mathbf{x}}) = \begin{bmatrix} \frac{\lambda^k}{k} \newline \lambda^k x_2^k \end{bmatrix} \end{align*}

Thus, we must select a λk\lambda^k such that λkk\frac{\lambda^k}{k} \to some finite value >0> 0, as kk \to \infty.

But then, 0λkx2kλkk300 \leq \lambda^k x_2^k \leq \frac{\lambda^k}{k^3} \to 0 as kk \to \infty.

So, TS(xˉ)={pR2p10,p2=0}T_S(\bar{\mathbf{x}}) = \{\mathbf{p} \in \mathbb{R}^2 \mid p_1 \geq 0, p_2 = 0\}.

For the gradient cone,

{g1(x)=[3x121]T    g1(xˉ)=[01]Tg2(x)=[5x141]T    g2(xˉ)=[01]Tg3(x)=[01]T    g3(xˉ)=[01]T    {g1(xˉ)Tp=0p1+1p2=p20g2(xˉ)Tp=0p11p2=p20g3(xˉ)Tp=0p11p2=p20    G(xˉ)={pR2p20,p20}={pR2p2=0}\begin{align*} & \begin{cases} \nabla g_1(\mathbf{x}) = \begin{bmatrix} -3x_1^2 & 1 \end{bmatrix}^T \implies \nabla g_1(\bar{\mathbf{x}}) = \begin{bmatrix} 0 & 1 \end{bmatrix}^T \newline \nabla g_2(\mathbf{x}) = \begin{bmatrix} 5x_1^4 & -1 \end{bmatrix}^T \implies \nabla g_2(\bar{\mathbf{x}}) = \begin{bmatrix} 0 & -1 \end{bmatrix}^T \newline \nabla g_3(\mathbf{x}) = \begin{bmatrix} 0 & -1 \end{bmatrix}^T \implies \nabla g_3(\bar{\mathbf{x}}) = \begin{bmatrix} 0 & -1 \end{bmatrix}^T \end{cases} \implies \newline & \begin{cases} \nabla g_1(\bar{\mathbf{x}})^T \mathbf{p} = 0 \cdot p_1 + 1 \cdot p_2 = p_2 \leq 0 \newline \nabla g_2(\bar{\mathbf{x}})^T \mathbf{p} = 0 \cdot p_1 - 1 \cdot p_2 = -p_2 \leq 0 \newline \nabla g_3(\bar{\mathbf{x}})^T \mathbf{p} = 0 \cdot p_1 - 1 \cdot p_2 = -p_2 \leq 0 \end{cases} \implies \newline & G(\bar{\mathbf{x}}) = \{\mathbf{p} \in \mathbb{R}^2 \mid p_2 \leq 0, -p_2 \leq 0\} = \{\mathbf{p} \in \mathbb{R}^2 \mid p_2 = 0\} \end{align*}

We can see that TS(xˉ)G(xˉ)T_S(\bar{\mathbf{x}}) \subseteq G(\bar{\mathbf{x}}), but TS(xˉ)G(xˉ)T_S(\bar{\mathbf{x}}) \neq G(\bar{\mathbf{x}}).

Constraint Qualifications

We will now introduce some more constraint qualifications that are easier to check than Abadie’s constraint qualifications, but imply that Abadie’s constraint qualifications hold.

Definition 1 (Linear Independence Constraint Qualification (LICQ))

The Linear Independence Constraint Qualification (LICQ) holds at a point x\mathbf{x}^{\star} if gi(x),iI(x)\nabla g_i(\mathbf{x}^{\star}), i \in \mathcal{I}(\mathbf{x}^{\star}) are linearly independent.

Theorem 1LICQ    Abadie’s CQ\text{LICQ} \implies \text{Abadie's CQ}
Definition 2 (Affine Independence Constraint Qualification (Affine CQ))

The Affine Independence Constraint Qualification (Affine CQ) holds if all gi(x),iI(x)g_i(\mathbf{x}), i \in \mathcal{I}(\mathbf{x}^{\star}) are affine functions and gi(x),i=1,,m\nabla g_i(\mathbf{x}^{\star}), i = 1, \ldots, m are affine.

Theorem 2Affine CQ    Abadie’s CQ. (In all points xS)\text{Affine CQ} \implies \text{Abadie's CQ. (In all points } \mathbf{x} \in S\text{)}
Definition 3 (Slater’s Constraint Qualification (Slater’s CQ))

Slater’s Constraint Qualification (Slater’s CQ) holds if all gi(x),i=1,,mg_i(\mathbf{x}), i = 1, \ldots, m are convex functions and if an interior point exists, i.e.,  x0 such that gi(x0)<0,i=1,,m\exists \ \mathbf{x}_0 \ \text{such that} \ g_i(\mathbf{x}_0) < 0, i = 1, \ldots, m.

Theorem 3Slater’s CQ    Abadie’s CQ. (In all points xS)\text{Slater's CQ} \implies \text{Abadie's CQ. (In all points } \mathbf{x} \in S\text{)}

Karush-Kuhn-Tucker (KKT) Conditions with Equality Constraints

So far, we have only considered inequality constraints when discussing the KKT conditions.

Let’s see how we can extend the KKT conditions to also include equality constraints.

Consider the following optimization problem with both equality and inequality constraints,

min f(x)subject to gi(x)0,i=1,,m hj(x)=0,j=1,,l\begin{align*} \min \ & f(\mathbf{x}) \newline \text{subject to} \ & g_i(\mathbf{x}) \leq 0, i = 1, \ldots, m \newline \ & h_j(\mathbf{x}) = 0, j = 1, \ldots, l \end{align*}

Let’s do a trick and reformulate the problem as,

min f(x)subject to gi(x)0,i=1,,m hj(x)0,j=1,,l hj(x)0,j=1,,l\begin{align*} \min \ & f(\mathbf{x}) \newline \text{subject to} \ & g_i(\mathbf{x}) \leq 0, i = 1, \ldots, m \newline \ & h_j(\mathbf{x}) \leq 0, j = 1, \ldots, l \newline \ & -h_j(\mathbf{x}) \leq 0, j = 1, \ldots, l \end{align*}

Thus, the KKT conditions for this problem are,

{f(x)+i=1mμigi(x)+j=1lλ1,jhj(x)+j=1lλ2,j(hj(x))=j=1l(λ1,jλ2,j)hj(x)=0μigi(x)=0,i=1,,mAlways fulfilled for feasible x,{λ1,jhj(x)=0,j=1,,lλ2,j(hj(x))=0,j=1,,lμi,λ1,j,λ2,j0,i=1,,m,j=1,,l\begin{cases} \nabla f(\mathbf{x}^{\star}) + \sum_{i=1}^{m} \mu_i \nabla g_i(\mathbf{x}^{\star}) + \underbrace{\sum_{j=1}^{l} \lambda_{1, j} \nabla h_j(\mathbf{x}^{\star}) + \sum_{j=1}^{l} \lambda_{2, j} \nabla (-h_j(\mathbf{x}^{\star}))}_{= \sum_{j=1}^{l} (\lambda_{1, j} - \lambda_{2, j}) \nabla h_j(\mathbf{x}^{\star})} = 0 \newline \mu_i g_i(\mathbf{x}^{\star}) = 0, i = 1, \ldots, m \newline \text{Always fulfilled for feasible } \mathbf{x}^{\star}, \begin{cases} \lambda_{1, j} h_j(\mathbf{x}^{\star}) = 0, j = 1, \ldots, l \newline \lambda_{2, j} (-h_j(\mathbf{x}^{\star})) = 0, j = 1, \ldots, l \newline \end{cases} \newline \mu_i, \lambda_{1, j}, \lambda_{2, j} \geq 0, i = 1, \ldots, m, j = 1, \ldots, l \newline \end{cases}

Thus, let λj=λ1,jλ2,j\lambda_j = \lambda_{1, j} - \lambda_{2, j}, which means that λjR\lambda_j \in \mathbb{R}, then the KKT conditions can be written as,

{f(x)+i=1mμigi(x)+j=1lλjhj(x)=0μigi(x)=0,i=1,,mμi0,i=1,,m(*)\begin{cases} \nabla f(\mathbf{x}^{\star}) + \sum_{i=1}^{m} \mu_i \nabla g_i(\mathbf{x}^{\star}) + \sum_{j=1}^{l} \lambda_j \nabla h_j(\mathbf{x}^{\star}) = 0 \newline \mu_i g_i(\mathbf{x}^{\star}) = 0, i = 1, \ldots, m \newline \mu_i \geq 0, i = 1, \ldots, m \newline \end{cases} \tag{*}

Now to fully generalize and capture the KKT conditions with equality constraints, we need to define a new cone.

Definition 4 (Tangent Space Cone for Equality Constraints)H(x){pRnhj(x)Tp=0,j=1,,l}H(\mathbf{x}) \coloneqq \{\mathbf{p} \in \mathbb{R}^n \mid \nabla h_j(\mathbf{x})^T \mathbf{p} = 0, j = 1, \ldots, l\}

With this we can define Abadie’s constraint qualifications for problems with equality constraints.

Definition 5 (Abadie’s Constraint Qualifications with Equality Constraints)

Abadie’s constraint qualifications hold at a point xS\mathbf{x}^{\star} \in S if,

TS(x)=G(x)H(x)T_S(\mathbf{x}^{\star}) = G(\mathbf{x}^{\star}) \cap H(\mathbf{x}^{\star})
Theorem 4 (KKT Conditions with Equality Constraints)

Assume Abadie’s constraint qualifications hold in xS\mathbf{x}^{\star} \in S. Then,

x is a local minimum     The system (*) is solvable for μ and λ.\mathbf{x}^{\star} \text{ is a local minimum } \implies \text{The system (*) is solvable for } \boldsymbol{\mu} \text{ and } \boldsymbol{\lambda}.

We can now define the sufficient conditions for optimality with equality constraints.

Theorem 5 (Sufficient Conditions for Optimality with Equality Constraints)

Assume that f(x)f(\mathbf{x}) is a convex function on SS, gi(x),i=1,,m\forall g_i(\mathbf{x}), i = 1, \ldots, m are convex functions on SS, and hj(x),j=1,,l\forall h_j(\mathbf{x}), j = 1, \ldots, l are affine functions on SS. Then,

x is a KKT point     x is globally optimal.\mathbf{x}^{\star} \text{ is a KKT point } \implies \mathbf{x}^{\star} \text{ is globally optimal.}
Proof (Sufficient Conditions for Optimality with Equality Constraints)

Let x\mathbf{x}^{\star} be a KKT point, by convexity of gi(x),xSg_i(\mathbf{x}), \forall \mathbf{x} \in S and iI(x)i \in \mathcal{I}(\mathbf{x}^{\star}). By using the characterization of convexity (of C1C^1 functions),

gi(x)T(xx)gi(x)0gi(x)= 0 since iI(x)0\begin{align*} \nabla g_i(\mathbf{x}^{\star})^T (\mathbf{x} - \mathbf{x}^{\star}) & \leq \underbrace{g_i(\mathbf{x})}_{\leq 0} - \underbrace{g_i(\mathbf{x}^{\star})}_{= \ 0 \text{ since } i \in \mathcal{I}(\mathbf{x}^{\star})} \newline & \leq 0 \end{align*}

Since hj(x),j=1,,lh_j(\mathbf{x}), j = 1, \ldots, l are affine functions, we have, again by the characterization of convexity (of C1C^1 functions) 1For clarity, we can use strict equality here because, hj(x)=ajTx+bjhj(x)=aj    hj(x)T(xx)=ajTx+bj(ajTx+bj)ajT(xx)=ajTxajTxajT(xx)=ajT(xx)\begin{aligned} h_j(\mathbf{x}) & = \mathbf{a}_j^T \mathbf{x} + b_j \newline \nabla h_j(\mathbf{x}) & = \mathbf{a}_j \newline \implies \newline \nabla h_j(\mathbf{x}^{\star})^T (\mathbf{x} - \mathbf{x}^{\star}) & = \mathbf{a}_j^T \mathbf{x} + b_j - (\mathbf{a}_j^T \mathbf{x}^{\star} + b_j) \newline \mathbf{a}_j^T(\mathbf{x} - \mathbf{x}^{\star}) & = \mathbf{a}_j^T \mathbf{x} - \mathbf{a}_j^T \mathbf{x}^{\star} \newline \mathbf{a}_j^T(\mathbf{x} - \mathbf{x}^{\star}) & = \mathbf{a}_j^T(\mathbf{x} - \mathbf{x}^{\star}) \end{aligned},

hj(x)T(xx)=hj(x)= 0hj(x)= 0=0\begin{align*} \nabla h_j(\mathbf{x}^{\star})^T (\mathbf{x} - \mathbf{x}^{\star}) & = \underbrace{h_j(\mathbf{x})}_{= \ 0} - \underbrace{h_j(\mathbf{x}^{\star})}_{= \ 0} \newline & = 0 \end{align*}

Further, by convexity of f(x)f(\mathbf{x}), xS\forall \mathbf{x} \in S,

f(x)f(x)f(x)T(xx)(i=1mμi0gi(x)0000j=1lλjhj(x)= 0)T(xx)0\begin{align*} f(\mathbf{x}) - f(\mathbf{x}^{\star}) & \geq \nabla f(\mathbf{x}^{\star})^T (\mathbf{x} - \mathbf{x}^{\star}) \newline & \geq \left(\underbrace{- \underbrace{\sum_{i=1}^{m} \underbrace{\underbrace{\mu_i}_{\geq 0} \underbrace{\nabla g_i(\mathbf{x}^{\star})}_{\leq 0}}_{\leq 0}}_{\leq 0}}_{\geq 0} - \underbrace{\sum_{j=1}^{l} \lambda_j \nabla h_j(\mathbf{x}^{\star})}_{= \ 0}\right)^T (\mathbf{x} - \mathbf{x}^{\star}) \newline & \geq 0 \end{align*}

Thus, by definition of global optimality, x\mathbf{x}^{\star} is globally optimal.

Connection to Lagrangian Relaxation

We can also connect the KKT conditions to the Lagrangian relaxation.

Let’s first define the Lagrangian.

Definition 6 (Lagrangian)

For the problem,

min f(x)subject to gi(x)0,i=1,,m\begin{align*} \min \ & f(\mathbf{x}) \newline \text{subject to} \ & g_i(\mathbf{x}) \leq 0, i = 1, \ldots, m \newline \end{align*}

The Lagrangian is defined as,

L(x,μ)=f(x)+i=1mμigi(x)\mathcal{L}(\mathbf{x}, \boldsymbol{\mu}) = f(\mathbf{x}) + \sum_{i=1}^{m} \mu_i g_i(\mathbf{x})
Note

It is not hard to see that,

xL(x,μ)=f(x)+i=1mμigi(x)    The first KKT condition is xL(x,μ)=0\begin{align*} \nabla_{\mathbf{x}} \mathcal{L}(\mathbf{x}, \boldsymbol{\mu}) & = \nabla f(\mathbf{x}) + \sum_{i=1}^{m} \mu_i \nabla g_i(\mathbf{x}) \newline & \implies \text{The first KKT condition is } \nabla_{\mathbf{x}} \mathcal{L}(\mathbf{x}^{\star}, \boldsymbol{\mu}) = 0 \end{align*}

We will now define what relaxtions means.

Definition 7 (Relaxation)

For a problem,

(P){min f(x)subject to xS\begin{align*} (P) \quad & \begin{cases} \min \ & f(\mathbf{x}) \newline \text{subject to} \ & \mathbf{x} \in S \end{cases} \end{align*}

The problem,

(P){min fR(x)subject to xSR\begin{align*} (P^{\prime}) \quad & \begin{cases} \min \ & f_R(\mathbf{x}) \newline \text{subject to} \ & \mathbf{x} \in S_R \end{cases} \end{align*}

The problem (P)(P^{\prime}) is called a relaxation if,

  1. fR(x)f(x),xSf_R(\mathbf{x}) \leq f(\mathbf{x}), \forall \mathbf{x} \in S

  2. SSRS \subseteq S_R

Lemma 1 (Relaxation Lemma)

For μ0\boldsymbol{\mu} \geq \mathbf{0}, the problem,

min L(x,μ)subject to xRn\begin{align*} \min \ & \mathcal{L}(\mathbf{x}, \boldsymbol{\mu}) \newline \text{subject to} \ & \mathbf{x} \in \mathbb{R^n} \end{align*}

is a relaxation of the problem,

min f(x)subject to gi(x)0,i=1,,m\begin{align*} \min \ & f(\mathbf{x}) \newline \text{subject to} \ & g_i(\mathbf{x}) \leq 0, i = 1, \ldots, m \newline \end{align*}
Proof (Relaxation Lemma)

By definition, the set SS is a subset of Rn\mathbb{R}^n,

S{xRngi(x)0,i=1,,m}RnS \coloneqq \{\mathbf{x} \in \mathbb{R}^n \mid g_i(\mathbf{x}) \leq 0, i = 1, \ldots, m\} \subseteq \mathbb{R}^n

Further, for μ0\boldsymbol{\mu} \geq \mathbf{0} and xS\mathbf{x} \in S,

L(x,μ)=f(x)+i=1mμi0gi(x)00f(x)\begin{align*} \mathcal{L}(\mathbf{x}, \boldsymbol{\mu}) & = f(\mathbf{x}) + \underbrace{\sum_{i=1}^{m} \underbrace{\mu_i}_{\geq 0} \underbrace{g_i(\mathbf{x})}_{\leq 0}}_{\leq 0} \newline & \leq f(\mathbf{x}) \end{align*}
Theorem 6 (Relaxtion Theorem)

a) Let fRf_R^{\star} be the minimum value of (P)(P^{\prime}), and ff^{\star} be the minimum value of (P)(P). Then,

fRff_R^{\star} \leq f^{\star}

also called the lower bound property.

b) If (P)(P^{\prime}) is infeasible, then (P)(P) is infeasible.

c) If (P)(P^{\prime}) has an optimal solution xRS\mathbf{x}_R^{\star} \in S, and fR(xR)=f(xR)f_R(\mathbf{x}_R^{\star}) = f(\mathbf{x}_R^{\star}), then xR\mathbf{x}_R^{\star} is an optimal solution of (P)(P).