Duality and KKT
Duality and KKT
• sensitivity analysis
• generalized inequalities
1
ECEN 629 Shuguang Cui
Lagrangian
standard form problem, (for now) we don’t assume convexity
minimize f0 (x)
subject to fi (x) ≤ 0, i = 1, . . . , m
hi(x) = 0, i = 1, . . . , p
Lagrangian L : Rn+m → R
m
X p
X
L(x, λ, ν) = f0 (x) + λifi (x) + νihi(x)
i=1 i=1
2
ECEN 629 Shuguang Cui
example: LP
minimize cT x
subject to aTi x − bi ≤ 0, i = 1, . . . , m
m
T
λi(aTi x−bi ) = −bT λ+ (AT λ+ c)T x
X
Note that L(x, λ) = c x+
i=1
−bT λ T
if A λ + c = 0
hence g(λ) =
−∞ otherwise
3
ECEN 629 Shuguang Cui
4
ECEN 629 Shuguang Cui
maximize g(λ, ν)
subject to λ0
5
ECEN 629 Shuguang Cui
Strong duality
d ⋆ = p⋆
when strong duality holds, dual optimal λ⋆ serves as certificate of optimality for primal
optimal point x⋆
many conditions or constraint qualifications guarantee strong duality for convex problems
Slater’s condition: if primal problem is strictly feasible (and convex), i.e., there exists
x ∈ relint D with
fi (x) < 0, i = 1, . . . , m
hi(x) = 0, i = 1, . . . , p
then we have p⋆ = d⋆
6
ECEN 629 Shuguang Cui
maximize −bT λ
subject to AT λ + c = 0
λ0
for LP we have strong duality except in one (pathological) case: primal and dual both
infeasible (p⋆ = +∞, d⋆ = −∞)
7
ECEN 629 Shuguang Cui
(primal) QP
minimize xT P x
subject to Ax b
T T
we assume P ≻ 0 for simplicity Lagrangian is L(x, λ) = x P x + λ (Ax − b)
∇xL(x, λ) = 0 yields x = −(1/2)P −1 AT λ, hence dual function is
T −1 T T
g(λ) = −(1/4)λ AP A λ−b λ
dual of QP is
maximize −(1/4)λT AP −1AT λ − bT λ
subject to λ0
. . . another QP
8
ECEN 629 Shuguang Cui
minimize xT x
subject to Ax = b
A is fat, full rank (solution is x⋆ = AT (AAT )−1 b)
dual function is
“ ” 1
g(ν) = inf x x + ν (Ax − b) = − ν T AAT ν − bT ν
T T
x 4
dual problem is
maximize − 14 ν T AAT ν − bT ν
solution: ν ⋆ = −2(AAT )−1b
can check d⋆ = p⋆
9
ECEN 629 Shuguang Cui
idea: simple transformation of primal problem can lead to very different dual
primal problem:
m
T
X
minimize log exp(ai x − bi )
i=1
⋆
dual function is constant g = p (we have strong duality, but it’s useless)
10
ECEN 629 Shuguang Cui
let us introduce
• m new variables y1, . . . , ym
• m new equality constraints y = Ax − b
dual function
m
!
exp yi + ν T (Ax − b − y)
X
g(ν) = inf log
x,y
i=1
• infimum is −∞ if AT ν 6= 0
• assuming AT ν = 0, let’s minimize over y :
eyi
Pm yj = νi
j=1 e
νi log νi − bT ν
X
g(ν) = −
i
11
ECEN 629 Shuguang Cui
dual problem
T
X
maximize −b ν − νi log νi
i
T
subject to 1 ν = 1, (ν 0)
AT ν = 0
moral: trivial reformulation can yield different dual
12
ECEN 629 Shuguang Cui
Duality in algorithms
13
ECEN 629 Shuguang Cui
f0 (x(k) ) − p⋆
relative error = ≤ǫ
|p⋆|
stopping criterion:
„ « „ «
f0 (x(k) )−g(λ(k) ,ν (k) ) f0 (x(k) )−g(λ(k) ,ν (k) )
until g(λ(k) , ν (k)) > 0 & ≤ǫ or f0 (x(k) ) < 0 & ≤ǫ
g(λ(k) ,ν (k) ) −f0 (x(k) )
14
ECEN 629 Shuguang Cui
Complementary slackness
suppose x⋆, λ⋆, and ν ⋆ are primal, dual feasible with zero duality gap (hence, they are
primal, dual optimal)
m p m
!
⋆ ⋆ ⋆ ⋆ ⋆ ⋆ ⋆ ⋆
X X X
f0 (x ) = g(λ , ν ) = inf f0(x) + λi fi (x) + νi hi (x) ≤ f0(x )+ λi fi (x )
x
i=1 i=1 i=1
Pm
hence we have i=1 λ⋆ifi (x⋆ ) = 0, and so
⋆ ⋆
λi fi (x ) = 0, i = 1, . . . , m
15
ECEN 629 Shuguang Cui
suppose
• fi are differentiable
• x⋆, λ⋆ are (primal, dual) optimal, with zero duality gap
therefore
⋆ ⋆ ⋆ ⋆ ⋆
X X
∇f0 (x ) + λi ∇fi (x ) + νi ∇hi(x ) = 0
i i
16
ECEN 629 Shuguang Cui
so if x⋆ , λ⋆, and ν ⋆ are (primal, dual) optimal, with zero duality gap, they satisfy
fi (x⋆ ) ≤ 0
hi (x⋆ ) = 0
λ⋆i ≥ 0
λ⋆ifi (x⋆ ) = 0
∇f0(x⋆ ) + ⋆ ⋆
νi⋆∇hi(x⋆ ) = 0
P P
i λi ∇fi (x ) + i
conversely, if the problem is convex and x⋆, λ⋆ satisfy KKT, then they are (primal, dual)
optimal
17
ECEN 629 Shuguang Cui
consider set
m+1
A = { (u, t) ∈ R | ∃x fi (x) ≤ ui , f0 (x) ≤ t }
• A is convex if fi are
• for λ 0, ( » –T » – ˛˛ » – )
λ u u
g(λ) = inf ∈A
˛
1 t t
˛
˛
t
A
T
t + λ u = g(λ)
g(λ)
» –
λ
1
u
18
ECEN 629 Shuguang Cui
p⋆ A
λ⋆
» –
1 u
19
ECEN 629 Shuguang Cui
epi p⋆
p⋆ (u)
p⋆ (0) − λ⋆T u
0
u
Pm
λ⋆ gives lower bound on p⋆(u): p⋆(u) ≥ p⋆ − i=1 λ⋆iui
⋆ ⋆ ∂p⋆(0) ⋆
if p (u) is differentiable, λi =− , λi is sensitivity of p⋆ w.r.t. ith constraint
∂ui
20
ECEN 629 Shuguang Cui
Generalized inequalities
minimize f0 (x)
subject to fi (x) Ki 0, i = 1, . . . , L
dual function
“ ”
g(λ1 , . . . , λL) = inf f0 (x) + λT1 f1 (x) + ··· + λTL fL(x)
x
21
ECEN 629 Shuguang Cui
dual problem
maximize g(λ1 , . . . , λL)
subject to λi K ⋆ 0, i = 1, . . . , L
i
then d⋆ = p⋆
22
ECEN 629 Shuguang Cui
dual problem
maximize Tr F0Z
subject to Tr FiZ + ci = 0, i = 1, . . . , n
Z = ZT 0
strong duality holds if there exists x with F0 + x1 F1 + · · · + xnFn ≺ 0
23