Lagrangian Duality: Jos Ededon A
Lagrangian Duality: Jos Ededon A
José De Doná
September 2004
3 Weak Duality
4 Strong Duality
Example
Primal Problem P
minimise f (x ), (1)
subject to:
gi (x ) ≤ 0 for i = 1, . . . , m,
hi (x ) = 0 for i = 1, . . . , ,
x ∈ X.
where,
m
θ(u, v ) = inf{f (x ) + ui gi (x ) + vi hi (x ) : x ∈ X }, (3)
i =1 i =1
Primal Problem P z
X (g , xxxxxxxxxxxxxxxxxxxxxx
f)
xxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxx
minimise f (x ), x xxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxx
subject to: xxxxxxxxxxxxxxxxxxxxxx
G
xxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxx
[g (x ), f (x )]
g ( x ) ≤ 0, xxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxx
x ∈ X,
where f : Rn → R and
g : Rn → R. y
z
X (g , xxxxxxxxxxxxxxxxxxxxxx
f)
Primal Problem P xxxxxxxxxxxxxxxxxxxxxx
x xxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxx
G
minimise f (x ), xxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxx
[g (x ), f (x )]
subject to: xxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxx
(ȳ , z̄ ) xxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxx
g ( x ) ≤ 0, xxxxxxxxxxxxxxxxxxxxxx
x ∈ X.
y
Finally, to solve the dual problem, we have to find the line with
slope −u (u ≥ 0) such that the last intercept on the z -axis, θ(u), is
maximal. Such a line has slope −ū and supports the set G at the
point (ȳ , z̄ ). Thus, the solution to the dual problem is ū, and the
optimal dual objective value is z̄ .
Centre of Complex Dynamic
Systems and Control
Geometric Interpretation
z
X (g , xxxxxxxxxxxxxxxxxxxxxx
f)
xxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxx
x xxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxx
G
xxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxx
[g (x ), f (x )]
xxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxx
(ȳ , z̄ ) xxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxx
Slope −ū
z + uy = α
θ(u) Slope −u
y
The following result shows that the objective value of any feasible
solution to the dual problem constitutes a lower bound for the
objective value of any feasible solution to the primal problem.
f (x ) ≥ θ(u, v ).
◦
Proof.
We use the definition of θ given in (3), and the facts that x ∈ X ,
u ≥ 0, g (x ) ≤ 0 and h (x ) = 0. We then have
Corollary
Note from the corollary that the optimal objective value of the
primal problem is greater than or equal to the optimal objective
value of the dual problem.
X (g , f )xxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxx Notice that, in the
x xxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxx case shown in the
G
xxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxx figure, there exists
[g (x ), f (x )]
xxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxx a duality gap due to
Duality gap xxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxx
Optimal primal objective the nonconvexity of
xxxxxxxxxxxxxxxxxxxxxxxxxx
the set G .
Optimal dual objective
Example
h
Consider the functions: 3
{α(x ), h (x ) : x ∈ X }
α(x ) = (x − 1) −
2 1
4,
h (x ) = 2x − 1, α
X = {x ∈ R : |x | ≤ 2}. −5
Example (Ctd.)
r
3 xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxx
α(x ) = (x − 1) −
2 1
4, S
h (x ) = 2x − 1, p
X = {x ∈ R : |x | ≤ 2}.
−5 xxxxx
Continuing with the proof of the Lemma, we have the convex set:
Corollary
Let S be a nonempty convex set in Rn and x̄ int S. Then there is
a nonzero vector p such that p t (x − x̄ ) ≤ 0 for each x ∈ cl S.
◦
We then have, from the above corollary, that there exists a nonzero
vector (u0 , u, v ) such that
for each (p , q, r ) ∈ cl S .
Example (Ctd.)
r
3 xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxx
h (x ) = 2x − 1, p
X = {x ∈ R : |x | ≤ 2}.
−5 xxxxxx
for each (p , q, r ) ∈ cl S.
r
3 xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxx
h (x ) = 2x − 1, p
S
X = {x ∈ R : |x | ≤ 2}.
−5 xxxxxx
for each (p , q, r ) ∈ cl S.
Centre of Complex Dynamic
Systems and Control
Proof of the Lemma (Ctd.)
for each (p , q, r ) ∈ cl S .
Hence,
h
{α(x ), h (x ) : x ∈ X }
3
( u0 , v )
α(x ) = (x − 1)2 − 14 ,
h (x ) = 2x − 1, α
X = {x ∈ R : |x | ≤ 2}.
−5
If System 2: u0 α(x ) + v t h (x ) ≥ 0 for some (u0 , v ) (0, 0), u0 ≥ 0
and for all x ∈ X , has a solution such that u0 > 0, and x ∈ X is
such that h (x ) = 0, we can see that α(x ) must be α(x ) ≥ 0, and
hence System 1 has no solution.
Centre of Complex Dynamic
Systems and Control
Strong Duality
The following result, known as the strong duality theorem, shows
that, under suitable convexity assumptions and under a constraint
qualification, there is no duality gap between the primal and dual
optimal objective function values.
Theorem (Strong Duality Theorem)
Let X be a nonempty convex set in Rn . Let f : Rn → R and
g : Rn → Rm be convex, and h : Rn → R be affine. Suppose that
the following constraint qualification is satisfied. There exists an
x̂ ∈ X such that g (x̂ ) < 0 and h (x̂ ) = 0, and 0 ∈ int h (X ), where
h (X ) = {h (x ) : x ∈ X }. Then,
Example
Consider the following optimisation problem:
Primal Problem P
minimise (x − 1)2 ,
subject to:
2x − 1 = 0,
x ∈ X = {x ∈ R : |x | ≤ 2}.
Example (Ctd.)
Lagrangian Dual Problem D
maximise θ(v ),
where the Lagrangian dual function is,
θ(v ) = inf{(x − 1)2 + v (2x − 1) : |x | ≤ 2}.
Differentiating w.r.t. x and equating to zero, we get that the
optimiser of the dual Lagrangian subproblem is x ∗ = −v + 1 (if
−1 ≤ v ≤ 3).
Hence θ(v ) = (−v + 1 − 1)2 + v (−2v + 2 − 1) = −v 2 + v.
Differentiating w.r.t. v and equating to zero, we get that the
1
optimiser of the dual problem is v ∗ = and the optimal value of
2
∗ ∗ 1
the dual problem is −v + v = . Thus, there is no duality gap.
2
4
f
X
[ ]
-2 2
[h (x ), f (x )] = 2x − 1, (x − 1)2
-5 3 h
Optimal primal
Slope = −v ∗ = −1/2
and dual objectives