HOME

TheInfoList



OR:

In mathematical optimization, the Rosenbrock function is a non- convex function, introduced by Howard H. Rosenbrock in 1960, which is used as a performance test problem for optimization
algorithm In mathematics and computer science, an algorithm () is a finite sequence of rigorous instructions, typically used to solve a class of specific problems or to perform a computation. Algorithms are used as specifications for performing ...
s. It is also known as Rosenbrock's valley or Rosenbrock's banana function. The global minimum is inside a long, narrow, parabolic shaped flat valley. To find the valley is trivial. To converge to the global
minimum In mathematical analysis, the maxima and minima (the respective plurals of maximum and minimum) of a function, known collectively as extrema (the plural of extremum), are the largest and smallest value of the function, either within a given r ...
, however, is difficult. The function is defined by f(x, y) = (a-x)^2 + b(y-x^2)^2 It has a global minimum at (x, y)=(a, a^2), where f(x, y)=0. Usually, these parameters are set such that a = 1 and b = 100. Only in the trivial case where a=0 the function is symmetric and the minimum is at the origin.


Multidimensional generalisations

Two variants are commonly encountered. One is the sum of N/2 uncoupled 2D Rosenbrock problems, and is defined only for even Ns: : f(\mathbf) = f(x_1, x_2, \dots, x_N) = \sum_^ \left 00(x_^2 - x_)^2 + (x_ - 1)^2 \right This variant has predictably simple solutions. A second, more involved variant is : f(\mathbf) = \sum_^ 00 (x_ - x_i^2 )^2 + (1-x_i)^2\quad \mbox \quad \mathbf = (x_1, \ldots, x_N) \in \mathbb^N. has exactly one minimum for N=3 (at (1, 1, 1)) and exactly two minima for 4 \le N \le 7—the global minimum at (1, 1, ..., 1) and a local minimum near \hat = (-1, 1, \dots, 1). This result is obtained by setting the gradient of the function equal to zero, noticing that the resulting equation is a rational function of x. For small N the polynomials can be determined exactly and Sturm's theorem can be used to determine the number of real roots, while the roots can be bounded in the region of , x_i, < 2.4. For larger N this method breaks down due to the size of the coefficients involved.


Stationary points

Many of the stationary points of the function exhibit a regular pattern when plotted. This structure can be exploited to locate them.


Optimization examples

The Rosenbrock function can be efficiently optimized by adapting appropriate coordinate system without using any gradient information and without building local approximation models (in contrast to many derivate-free optimizers). The following figure illustrates an example of 2-dimensional Rosenbrock function optimization by
adaptive coordinate descent Adaptive coordinate descent is an improvement of the coordinate descent algorithm to non-separable optimization by the use of adaptive encoding. The adaptive coordinate descent approach gradually builds a transformation of the coordinate system suc ...
from starting point x_0=(-3,-4). The solution with the function value 10^ can be found after 325 function evaluations. Using the
Nelder–Mead method The Nelder–Mead method (also downhill simplex method, amoeba method, or polytope method) is a numerical method used to find the minimum or maximum of an objective function in a multidimensional space. It is a direct search method (based on ...
from starting point x_0=(-1,1) with a regular initial simplex a minimum is found with function value 1.36 \cdot 10^ after 185 function evaluations. The figure below visualizes the evolution of the algorithm.


See also

*
Test functions for optimization In applied mathematics, test functions, known as artificial landscapes, are useful to evaluate characteristics of optimization algorithms, such as: * Convergence rate. * Precision. * Robustness. * General performance. Here some test functions are ...


References


External links


Rosenbrock function plot in 3D
* {{MathWorld , title=Rosenbrock Function , urlname=RosenbrockFunction Mathematical optimization Polynomials Functions and mappings