Overview

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+x2≤0g2(x)=x15−x2≤0g3(x)=−x2≤0S = \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
Figure 1: Feasible region and constraints

Let x1k=1kx_1^k = \frac{1}{k}, by g1(x)g_1(\mathbf{x}), x2k≤1k3x_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,

lim⁡k→∞xk=xˉλk(xk−xˉ)=[λ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 k→∞k \to \infty.

But then, 0≤λkx2k≤λkk3→00 \leq \lambda^k x_2^k \leq \frac{\lambda^k}{k^3} \to 0 as k→∞k \to \infty.

So, TS(xˉ)={p∈R2∣p1≥0,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]T∇g2(x)=[5x14−1]T  ⟹  ∇g2(xˉ)=[0−1]T∇g3(x)=[0−1]T  ⟹  ∇g3(xˉ)=[0−1]T  ⟹  {∇g1(xˉ)Tp=0⋅p1+1⋅p2=p2≤0∇g2(xˉ)Tp=0⋅p1−1⋅p2=−p2≤0∇g3(xˉ)Tp=0⋅p1−1⋅p2=−p2≤0  ⟹  G(xˉ)={p∈R2∣p2≤0,−p2≤0}={p∈R2∣p2=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⋆),i∈I(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),i∈I(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 x∈S)\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 thatgi(x0)<0,i=1,…,m\exists\, \mathbf{x}_0 \, \text{such that} \quad g_i(\mathbf{x}_0) < 0, i = 1, \ldots, m.

Theorem 3Slater’s CQ  ⟹  Abadie’s CQ. (In all points x∈S)\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 togi(x)≤0,i=1,…,m hj(x)=0,j=1,…,l\begin{align*} \min & f(\mathbf{x}) \newline \text{subject to} \quad & 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 togi(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} \quad & 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μi∇gi(x⋆)+∑j=1lλ1,j∇hj(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,j≥0,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 λj∈R\lambda_j \in \mathbb{R}, then the KKT conditions can be written as,

{∇f(x⋆)+∑i=1mμi∇gi(x⋆)+∑j=1lλj∇hj(x⋆)=0μigi(x⋆)=0,i=1,…,mμi≥0,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)≔{p∈Rn∣∇hj(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 x⋆∈S\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 x⋆∈S\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),∀x∈Sg_i(\mathbf{x}), \forall \mathbf{x} \in S and i∈I(x⋆)i \in \mathcal{I}(\mathbf{x}^{\star}). By using the characterization of convexity (of C1C^1 functions),

∇gi(x⋆)T(x−x⋆)≤gi(x)⏟≤0−gi(x⋆)⏟= 0 since i∈I(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) 1,

∇hj(x⋆)T(x−x⋆)=hj(x)⏟= 0−hj(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}), ∀x∈S\forall \mathbf{x} \in S,

f(x)−f(x⋆)≥∇f(x⋆)T(x−x⋆)≥(−∑i=1mμi⏟≥0∇gi(x⋆)⏟≤0⏟≤0⏟≤0⏟≥0−∑j=1lλj∇hj(x⋆)⏟= 0)T(x−x⋆)≥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 togi(x)≤0,i=1,…,m\begin{align*} \min & f(\mathbf{x}) \newline \text{subject to} \quad & 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μi∇gi(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 tox∈S\begin{align*} (P) \quad & \begin{cases} \min & f(\mathbf{x}) \newline \text{subject to} \quad & \mathbf{x} \in S \end{cases} \end{align*}

The problem,

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

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

  1. fR(x)≤f(x),∀x∈Sf_R(\mathbf{x}) \leq f(\mathbf{x}), \forall \mathbf{x} \in S

  2. S⊆SRS \subseteq S_R

Lemma 1 (Relaxation Lemma)

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

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

is a relaxation of the problem,

min⁡f(x)subject togi(x)≤0,i=1,…,m\begin{align*} \min & f(\mathbf{x}) \newline \text{subject to} \quad & 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≔{x∈Rn∣gi(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 x∈S\mathbf{x} \in S,

L(x,μ)=f(x)+∑i=1mμi⏟≥0gi(x)⏟≤0⏟≤0≤f(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 fR⋆f_R^{\star} be the minimum value of (P′)(P^{\prime}), and f⋆f^{\star} be the minimum value of (P)(P). Then,

fR⋆≤f⋆f_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 xR⋆∈S\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).

Footnotes
  1. For clarity, we can use strict equality here because,

    hj(x)=ajTx+bj∇hj(x)=aj  ⟹  ∇hj(x⋆)T(x−x⋆)=ajTx+bj−(ajTx⋆+bj)ajT(x−x⋆)=ajTx−ajTx⋆ajT(x−x⋆)=ajT(x−x⋆)\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}↩

Search

Essays, lecture notes, and research

Type to search the site.