Proximal Operator
   HOME

TheInfoList



OR:

In
mathematical optimization Mathematical optimization (alternatively spelled ''optimisation'') or mathematical programming is the selection of a best element, with regard to some criterion, from some set of available alternatives. It is generally divided into two subfi ...
, the proximal operator is an operator associated with a proper,An (extended) real-valued function ''f'' on a
Hilbert space In mathematics, Hilbert spaces (named after David Hilbert) allow generalizing the methods of linear algebra and calculus from (finite-dimensional) Euclidean vector spaces to spaces that may be infinite-dimensional. Hilbert spaces arise natural ...
is said to be ''proper'' if it is not identically equal to +\infty, and -\infty is not in its image.
lower semi-continuous In mathematical analysis, semicontinuity (or semi-continuity) is a property of extended real-valued functions that is weaker than continuity. An extended real-valued function f is upper (respectively, lower) semicontinuous at a point x_0 if, rou ...
convex function In mathematics, a real-valued function is called convex if the line segment between any two points on the graph of a function, graph of the function lies above the graph between the two points. Equivalently, a function is convex if its epigra ...
f from a
Hilbert space In mathematics, Hilbert spaces (named after David Hilbert) allow generalizing the methods of linear algebra and calculus from (finite-dimensional) Euclidean vector spaces to spaces that may be infinite-dimensional. Hilbert spaces arise natural ...
\mathcal to \infty,+\infty/math>, and is defined by: ::\operatorname_f(v) = \arg \min_ \left(f(x) + \frac 1 2 \, x - v\, _\mathcal^2\right). For any function in this class, the minimizer of the right-hand side above is unique, hence making the proximal operator well-defined. The proximal operator is used in proximal gradient methods, which is frequently used in optimization algorithms associated with non-
differentiable In mathematics, a differentiable function of one real variable is a function whose derivative exists at each point in its domain. In other words, the graph of a differentiable function has a non-vertical tangent line at each interior point in its ...
optimization problems such as
total variation denoising In signal processing, particularly image processing, total variation denoising, also known as total variation regularization or total variation filtering, is a noise removal process (filter). It is based on the principle that signals with excessi ...
.


Properties

The \text of a proper, lower semi-continuous convex function f enjoys several useful properties for optimization. * Fixed points of \text_f are minimizers of f: \ = \arg \min f. * Global convergence to a minimizer is defined as follows: If \arg \min f \neq \varnothing, then for any initial point x_0 \in \mathcal, the recursion (\forall n \in \mathbb)\quad x_ = \text_f x_n yields convergence x_n \to x \in \arg \min f as n \to +\infty. This convergence may be weak if \mathcal is infinite dimensional. * The proximal operator can be seen as a generalization of the
projection operator In linear algebra and functional analysis, a projection is a linear transformation P from a vector space to itself (an endomorphism) such that P\circ P=P. That is, whenever P is applied twice to any vector, it gives the same result as if it wer ...
. Indeed, in the specific case where f is the 0-\infty indicator function \iota_C of a nonempty, closed, convex set C we have that : \begin \operatorname_(x) &= \operatorname\limits_y \begin \frac \left\, x-y \right\, _2^2 & \text y \in C \\ + \infty & \text y \notin C \end \\ &=\operatorname\limits_ \frac \left\, x-y \right\, _2^2 \end : showing that the proximity operator is indeed a generalisation of the projection operator. * A function is firmly non-expansive if (\forall (x,y) \in \mathcal^2) \quad \, \text_fx - \text_fy\, ^2 \leq \langle x-y\ , \ \text_fx - \text_fy\rangle. * The proximal operator of a function is related to the gradient of the
Moreau envelope The Moreau envelope (or the Moreau-Yosida regularization) M_f of a proper lower semi-continuous convex function f is a smoothed version of f. It was proposed by Jean-Jacques Moreau in 1965. The Moreau envelope has important applications in mathemat ...
M_ of a function \lambda f by the following identity: \nabla M_(x) = \frac (x - \mathrm_(x)). * The proximity operator of f is characterized by inclusion p=\operatorname_f(x) \Leftrightarrow x-p \in \partial f(p) , where \partial f is the
subdifferential In mathematics, the subderivative, subgradient, and subdifferential generalize the derivative to convex functions which are not necessarily differentiable. Subderivatives arise in convex analysis, the study of convex functions, often in connection ...
of f, given by : \partial f(x) = \ In particular, If f is differentiable then the above equation reduces to p=\operatorname_f(x) \Leftrightarrow x-p = \nabla f(p) .


Notes


References

{{Reflist Mathematical optimization


See also

*
Proximal gradient method Proximal gradient methods are a generalized form of projection used to solve non-differentiable convex optimization problems. Many interesting problems can be formulated as convex optimization problems of the form \operatorname\limits_ \sum_^n ...


External links

* Th
Proximity Operator repository
a collection of proximity operators implemented in
Matlab MATLAB (an abbreviation of "MATrix LABoratory") is a proprietary multi-paradigm programming language and numeric computing environment developed by MathWorks. MATLAB allows matrix manipulations, plotting of functions and data, implementation ...
and
Python Python may refer to: Snakes * Pythonidae, a family of nonvenomous snakes found in Africa, Asia, and Australia ** ''Python'' (genus), a genus of Pythonidae found in Africa and Asia * Python (mythology), a mythical serpent Computing * Python (pro ...
.
ProximalOperators.jl
a
Julia Julia is usually a feminine given name. It is a Latinate feminine form of the name Julio and Julius. (For further details on etymology, see the Wiktionary entry "Julius".) The given name ''Julia'' had been in use throughout Late Antiquity (e.g. ...
package implementing proximal operators.
ODL
a Python library for
inverse problems An inverse problem in science is the process of calculating from a set of observations the causal factors that produced them: for example, calculating an image in X-ray computed tomography, source reconstruction in acoustics, or calculating the ...
that utilizes proximal operators.