Convex analysis

From HandWiki
Revision as of 19:58, 6 March 2023 by Sherlock (talk | contribs) (url)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
A 3-dimensional convex polytope. Convex analysis includes not only the study of convex subsets of Euclidean spaces but also the study of convex functions on abstract spaces.

Convex analysis is the branch of mathematics devoted to the study of properties of convex functions and convex sets, often with applications in convex minimization, a subdomain of optimization theory.

Convex sets

Main page: Convex set

A subset [math]\displaystyle{ C \subseteq X }[/math] of some vector space [math]\displaystyle{ X }[/math] is convex if it satisfies any of the following equivalent conditions:

  1. If [math]\displaystyle{ 0 \leq r \leq 1 }[/math] is real and [math]\displaystyle{ x, y \in C }[/math] then [math]\displaystyle{ r x + (1 - r) y \in C. }[/math][1]
  2. If [math]\displaystyle{ 0 \lt r \lt 1 }[/math] is real and [math]\displaystyle{ x, y \in C }[/math] with [math]\displaystyle{ x \neq y, }[/math] then [math]\displaystyle{ r x + (1 - r) y \in C. }[/math]
Convex function on an interval.
Main page: Convex function

Throughout, [math]\displaystyle{ f : X \to [-\infty, \infty] }[/math] will be a map valued in the extended real numbers [math]\displaystyle{ [-\infty, \infty] = \mathbb{R} \cup \{ \pm \infty \} }[/math] with a domain [math]\displaystyle{ \operatorname{domain} f = X }[/math] that is a convex subset of some vector space. The map [math]\displaystyle{ f : X \to [-\infty, \infty] }[/math] is a convex function if

[math]\displaystyle{ f( r x + (1 - r) y ) \leq r f(x) + (1 - r) f(y) }[/math]

 

 

 

 

(Convexity ≤)

holds for any real [math]\displaystyle{ 0 \lt r \lt 1 }[/math] and any [math]\displaystyle{ x, y \in X }[/math] with [math]\displaystyle{ x \neq y. }[/math] If this remains true of [math]\displaystyle{ f }[/math] when the defining inequality (Convexity ≤) is replaced by the strict inequality

[math]\displaystyle{ f( r x + (1 - r) y ) \lt r f(x) + (1 - r) f(y) }[/math]

 

 

 

 

(Convexity <)

then [math]\displaystyle{ f }[/math] is called strictly convex.[1]

Convex functions are related to convex sets. Specifically, the function [math]\displaystyle{ f }[/math] is convex if and only if its epigraph

A function (in black) is convex if and only if its epigraph, which is the region above its graph (in green), is a convex set.
A graph of the bivariate convex function [math]\displaystyle{ x^2 + x y + y^2. }[/math]

[math]\displaystyle{ \operatorname{epi} f := \left\{ (x,r) \in X \times \mathbb{R} ~:~ f(x) \leq r \right\} }[/math]

 

 

 

 

(Epigraph def.)

is a convex set.[2] The epigraphs of extended real-valued functions play a role in convex analysis that is analogous to the role played by graphs of real-valued function in real analysis. Specifically, the epigraph of an extended real-valued function provides geometric intuition that can be used to help formula or prove conjectures.

The domain of a function [math]\displaystyle{ f : X \to [-\infty, \infty] }[/math] is denoted by [math]\displaystyle{ \operatorname{domain} f }[/math] while its effective domain is the set[2]

[math]\displaystyle{ \operatorname{dom} f := \{ x \in X ~:~ f(x) \lt \infty \}. }[/math]

 

 

 

 

(dom f def.)

The function [math]\displaystyle{ f : X \to [-\infty, \infty] }[/math] is called proper if [math]\displaystyle{ \operatorname{dom} f \neq \varnothing }[/math] and [math]\displaystyle{ f(x) \gt -\infty }[/math] for all [math]\displaystyle{ x \in \operatorname{domain} f. }[/math][2] Alternatively, this means that there exists some [math]\displaystyle{ x }[/math] in the domain of [math]\displaystyle{ f }[/math] at which [math]\displaystyle{ f(x) \in \mathbb{R} }[/math] and [math]\displaystyle{ f }[/math] is also never equal to [math]\displaystyle{ -\infty. }[/math] In words, a function is proper if its domain is not empty, it never takes on the value [math]\displaystyle{ -\infty, }[/math] and it also is not identically equal to [math]\displaystyle{ +\infty. }[/math] If [math]\displaystyle{ f : \mathbb{R}^n \to [-\infty, \infty] }[/math] is a proper convex function then there exist some vector [math]\displaystyle{ b \in \mathbb{R}^n }[/math] and some [math]\displaystyle{ r \in \mathbb{R} }[/math] such that

[math]\displaystyle{ f(x) \geq x \cdot b - r }[/math]     for every [math]\displaystyle{ x }[/math]

where [math]\displaystyle{ x \cdot b }[/math] denotes the dot product of these vectors.

Convex conjugate

Main page: Convex conjugate

The convex conjugate of an extended real-valued function [math]\displaystyle{ f : X \to [-\infty, \infty] }[/math] (not necessarily convex) is the function [math]\displaystyle{ f^* : X^* \to [-\infty, \infty] }[/math] from the (continuous) dual space [math]\displaystyle{ X^* }[/math] of [math]\displaystyle{ X, }[/math] and[3]

[math]\displaystyle{ f^*\left(x^*\right) = \sup_{z \in X} \left\{ \left\langle x^*, z \right\rangle - f(z) \right\} }[/math]

where the brackets [math]\displaystyle{ \left\langle \cdot, \cdot \right\rangle }[/math] denote the canonical duality [math]\displaystyle{ \left\langle x^*, z \right\rangle := x^*(z). }[/math] The biconjugate of [math]\displaystyle{ f }[/math] is the map [math]\displaystyle{ f^{**} = \left( f^* \right)^* : X \to [-\infty, \infty] }[/math] defined by [math]\displaystyle{ f^{**}(x) := \sup_{z^* \in X^*} \left\{ \left\langle x, z^* \right\rangle - f\left( z^* \right) \right\} }[/math] for every [math]\displaystyle{ x \in X. }[/math] If [math]\displaystyle{ \operatorname{Func}(X; Y) }[/math] denotes the set of [math]\displaystyle{ Y }[/math]-valued functions on [math]\displaystyle{ X, }[/math] then the map [math]\displaystyle{ \operatorname{Func}(X; [-\infty, \infty]) \to \operatorname{Func}\left( X^*; [-\infty, \infty] \right) }[/math] defined by [math]\displaystyle{ f \mapsto f^* }[/math] is called the Legendre-Fenchel transform.

Subdifferential set and the Fenchel-Young inequality

If [math]\displaystyle{ f : X \to [-\infty, \infty] }[/math] and [math]\displaystyle{ x \in X }[/math] then the subdifferential set is

[math]\displaystyle{ \begin{alignat}{4} \partial f(x) :&= \left\{ x^* \in X^* ~:~ f(z) \geq f(x) + \left\langle x^*, z - x \right\rangle \text{ for all } z \in X \right\} && (\text{“} z \in X \text{''} \text{ can be replaced with: } \text{“} z \in X \text{ such that } z \neq x \text{''}) \\ &= \left\{ x^* \in X^* ~:~ \left\langle x^*, x \right\rangle - f(x) \geq \left\langle x^*, z \right\rangle - f(z) \text{ for all } z \in X \right\} && \\ &= \left\{ x^* \in X^* ~:~ \left\langle x^*, x \right\rangle - f(x) \geq \sup_{z \in X} \left\langle x^*, z \right\rangle - f(z) \right\} && \text{ The right hand side is } f^*\left( x^* \right) \\ &= \left\{ x^* \in X^* ~:~ \left\langle x^*, x \right\rangle - f(x) = f^*\left( x^* \right) \right\} && \text{ Taking } z := x \text{ in the } \sup{} \text{ gives the inequality } \leq. \\ \end{alignat} }[/math]

For example, in the important special case where [math]\displaystyle{ f = \| \cdot \| }[/math] is a norm on [math]\displaystyle{ X }[/math], it can be shown[proof 1] that if [math]\displaystyle{ 0 \neq x \in X }[/math] then this definition reduces down to:

[math]\displaystyle{ \partial f (x) = \left\{ x^* \in X^* ~:~ \left\langle x^*, x \right\rangle = \| x \| \text{ and } \left\| x^* \right\| = 1 \right\} }[/math]     and     [math]\displaystyle{ \partial f(0) = \left\{ x^* \in X^* ~:~ \left\| x^* \right\| \leq 1 \right\}. }[/math]

For any [math]\displaystyle{ x \in X }[/math] and [math]\displaystyle{ x^* \in X^*, }[/math] [math]\displaystyle{ f(x) + f^*\left(x^*\right) \geq \left\langle x^*, x \right\rangle, }[/math] which is called the Fenchel-Young inequality. This inequality is an equality (i.e. [math]\displaystyle{ f(x) + f^*\left(x^*\right) = \left\langle x^*, x \right\rangle }[/math]) if and only if [math]\displaystyle{ x^* \in \partial f(x). }[/math] It is in this way that the subdifferential set [math]\displaystyle{ \partial f (x) }[/math] is directly related to the convex conjugate [math]\displaystyle{ f^*\left( x^* \right). }[/math]

Biconjugate

The biconjugate of a function [math]\displaystyle{ f : X \to [-\infty, \infty] }[/math] is the conjugate of the conjugate, typically written as [math]\displaystyle{ f^{**} : X \to [-\infty, \infty]. }[/math] The biconjugate is useful for showing when strong or weak duality hold (via the perturbation function).

For any [math]\displaystyle{ x \in X, }[/math] the inequality [math]\displaystyle{ f^{**}(x) \leq f(x) }[/math] follows from the Fenchel–Young inequality. For proper functions, [math]\displaystyle{ f = f^{**} }[/math] if and only if [math]\displaystyle{ f }[/math] is convex and lower semi-continuous by Fenchel–Moreau theorem.[3][4]

Convex minimization

Main page: Convex optimization

A convex minimization (primal) problem is one of the form

find [math]\displaystyle{ \inf_{x \in M} f(x) }[/math] when given a convex function [math]\displaystyle{ f : X \to [-\infty, \infty] }[/math] and a convex subset [math]\displaystyle{ M \subseteq X. }[/math]

Dual problem

Main page: Duality (optimization)

In optimization theory, the duality principle states that optimization problems may be viewed from either of two perspectives, the primal problem or the dual problem.

In general given two dual pairs separated locally convex spaces [math]\displaystyle{ \left(X, X^*\right) }[/math] and [math]\displaystyle{ \left(Y, Y^*\right). }[/math] Then given the function [math]\displaystyle{ f : X \to [-\infty, \infty], }[/math] we can define the primal problem as finding [math]\displaystyle{ x }[/math] such that

[math]\displaystyle{ \inf_{x \in X} f(x). }[/math]

If there are constraint conditions, these can be built into the function [math]\displaystyle{ f }[/math] by letting [math]\displaystyle{ f = f + I_{\mathrm{constraints}} }[/math] where [math]\displaystyle{ I }[/math] is the indicator function. Then let [math]\displaystyle{ F : X \times Y \to [-\infty, \infty] }[/math] be a perturbation function such that [math]\displaystyle{ F(x, 0) = f(x). }[/math][5]

The dual problem with respect to the chosen perturbation function is given by

[math]\displaystyle{ \sup_{y^* \in Y^*} -F^*\left(0, y^*\right) }[/math]

where [math]\displaystyle{ F^* }[/math] is the convex conjugate in both variables of [math]\displaystyle{ F. }[/math]

The duality gap is the difference of the right and left hand sides of the inequality[6][5][7]

[math]\displaystyle{ \sup_{y^* \in Y^*} -F^*\left(0, y^*\right) \le \inf_{x \in X} F(x, 0). }[/math]

This principle is the same as weak duality. If the two sides are equal to each other, then the problem is said to satisfy strong duality.

There are many conditions for strong duality to hold such as:

Lagrange duality

For a convex minimization problem with inequality constraints,

[math]\displaystyle{ \min {}_{x} f(x) }[/math] subject to [math]\displaystyle{ g_i(x) \leq 0 }[/math] for [math]\displaystyle{ i = 1, \ldots, m. }[/math]

the Lagrangian dual problem is

[math]\displaystyle{ \sup {}_{u} \inf {}_{x} L(x, u) }[/math] subject to [math]\displaystyle{ u_i(x) \geq 0 }[/math] for [math]\displaystyle{ i = 1, \ldots, m. }[/math]

where the objective function [math]\displaystyle{ L(x, u) }[/math] is the Lagrange dual function defined as follows:

[math]\displaystyle{ L(x, u) = f(x) + \sum_{j=1}^m u_j g_j(x) }[/math]

See also

Notes

  1. 1.0 1.1 Rockafellar, R. Tyrrell (1997). Convex Analysis. Princeton, NJ: Princeton University Press. ISBN 978-0-691-01586-6. 
  2. 2.0 2.1 2.2 Rockafellar & Wets 2009, pp. 1-28.
  3. 3.0 3.1 Zălinescu 2002, pp. 75-79.
  4. Borwein, Jonathan; Lewis, Adrian (2006). Convex Analysis and Nonlinear Optimization: Theory and Examples (2 ed.). Springer. pp. 76–77. ISBN 978-0-387-29570-1. https://archive.org/details/convexanalysisno00borw_812. 
  5. 5.0 5.1 Boţ, Radu Ioan; Wanka, Gert; Grad, Sorin-Mihai (2009). Duality in Vector Optimization. Springer. ISBN 978-3-642-02885-4. 
  6. Zălinescu 2002, pp. 106-113.
  7. Csetnek, Ernö Robert (2010). Overcoming the failure of the classical generalized interior-point regularity conditions in convex optimization. Applications of the duality theory to enlargements of maximal monotone operators. Logos Verlag Berlin GmbH. ISBN 978-3-8325-2503-3. 
  8. Borwein, Jonathan; Lewis, Adrian (2006). Convex Analysis and Nonlinear Optimization: Theory and Examples (2 ed.). Springer. ISBN 978-0-387-29570-1. 
  9. Boyd, Stephen; Vandenberghe, Lieven (2004). Convex Optimization. Cambridge University Press. ISBN 978-0-521-83378-3. https://web.stanford.edu/~boyd/cvxbook/bv_cvxbook.pdf. Retrieved October 3, 2011. 
  1. The conclusion is immediate if [math]\displaystyle{ X = \{ 0 \} }[/math] so assume otherwise. Fix [math]\displaystyle{ x \in X. }[/math] Replacing [math]\displaystyle{ f }[/math] with the norm gives [math]\displaystyle{ \partial f(x) = \left\{ x^* \in X^* ~:~ \left\langle x^*, x \right\rangle - \| x \| \geq \left\langle x^*, z \right\rangle - \| z \| \text{ for all } z \in X \right\}. }[/math] If [math]\displaystyle{ x^* \in \partial f(x) }[/math] and [math]\displaystyle{ r \geq 0 }[/math] is real then using [math]\displaystyle{ z := r x }[/math] gives [math]\displaystyle{ \left\langle x^*, x \right\rangle - \| x \| \geq \left\langle x^*, r x \right\rangle - \| r x \| = r \left[ \left\langle x^*, x \right\rangle - \| x \| \right], }[/math] where in particular, taking [math]\displaystyle{ r := 2 }[/math] gives [math]\displaystyle{ x^*(x) \geq \| x \| }[/math] while taking [math]\displaystyle{ r := \frac1{2} }[/math] gives [math]\displaystyle{ x^*(x) \leq \| x \| }[/math] and thus [math]\displaystyle{ x^*(x) = \| x \| }[/math]; moreover, if in addition [math]\displaystyle{ x \neq 0 }[/math] then because [math]\displaystyle{ x^*\left(\frac{x}{\|x\|}\right) = 1, }[/math] it follows from the definition of the dual norm that [math]\displaystyle{ \left\| x^* \right\| \geq 1. }[/math] Because [math]\displaystyle{ \partial f(x) \subseteq \left\{ x^* \in X^* ~:~ x^*(x) = \| x \| \right\}, }[/math] which is equivalent to [math]\displaystyle{ \partial f(x) = \partial f(x) \cap \left\{ x^* \in X^* ~:~ x^*(x) = \| x \| \right\}, }[/math] it follows that [math]\displaystyle{ \partial f(x) = \left\{ x^* \in X^* ~:~ x^*(x) = \| x \| \text{ and } \| z \| \geq \left\langle x^*, z \right\rangle \text{ for all } z \in X \right\}, }[/math] which implies [math]\displaystyle{ \left\| x^* \right\| \leq 1 }[/math] for all [math]\displaystyle{ x^* \in \partial f(x). }[/math] From these facts, the conclusion can now be reached. ∎

References

External links