罗森布罗克函数

双变数的罗森布罗克函数

在数学最佳化中，罗森布罗克函数是一个用来测试最佳化演算法性能的非凸函数，由霍华德·哈里·罗森布罗克】在1960年提出。也称为罗森布罗克山谷或罗森布罗克香蕉函数，也简称为香蕉函数。

罗森布罗克函数的定义如下：

$f(x,y)=(1-x)^{2}+100(y-x^{2})^{2}.\quad$

罗森布罗克函数的每个等高线大致呈抛物线形，其全域最小值也位在抛物线形的山谷中（香蕉型山谷）。很容易找到这个山谷，但由于山谷内的值变化不大，要找到全域的最小值相当困难。

其全域最小值位于$(x,y)=(1,1)$点，数值为$f(x,y)=0$。有时第二项的系数不同，但不会影响全域最小值的位置。

用罗森布罗克函数测试梯度下降法的结果，约1000次才到达全域最小值

多变数下的扩展

多变数的罗森布罗克k函数有以下二种形式。一种是$N/2$个独立二维罗森布罗克函数的和：

$f(\mathbf {x} )=f(x_{1},x_{2},\dots ,x_{N})=\sum _{i=1}^{N/2}\left[100(x_{2i-1}^{2}-x_{2i})^{2}+(x_{2i-1}-1)^{2}\right].$

此形式只在$N$为偶数时有定义，而且其解较简单。

另一个较复杂的形式为：

$f(\mathbf {x} )=\sum _{i=1}^{N-1}\left[(1-x_{i})^{2}+100(x_{i+1}-x_{i}^{2})^{2}\right]\quad \forall \mathbf {x} \in \mathbb {R} ^{N}.$

可证明当$N=3$时，此形式的罗森布罗克函数只有一个最小值（位置在$(1,1,1)$），在 $4\leq N\leq 7$时只有二个最小值，所有变数均为1时有全域最小值，而在$(x_{1},x_{2},\dots ,x_{N})=(-1,1,\dots ,1)$附近有局部最小值。此结果是将令函数的梯度为0后求得，罗森布罗克函数的梯度仍为一个$x$的多项式，在$N$较小时，可以精确的列出多项式，再求出实根的个数，而其根限制在$|x_{i}|<2.4$的范围内。若$N$较大时因为相关的系数太多，无法用以上方式进行。

随机函数

有许多方式可以将罗森布罗克函数延伸到随机（stochastic）函数，以下是一种例子：

$f(\mathbf {x} )=\sum _{i=1}^{n-1}{\Big [}(1-x_{i})^{2}+100\epsilon _{i}(x_{i+1}-x_{i}^{2})^{2}{\Big ]},$

其中随机变数$\epsilon _{i}(i=1,2,...,n-1)$服从均匀分布 Unif(0,1)。原则上，此随机函数的全域最小值仍在(1,1,...,1)，不过因为其随机的特性，任何以梯度下降法为基础的最佳化演算法均无法用来求得此随机函数的最小值。

可适用的最佳化演算法

经若经过适当的坐标系调整，可以在没有梯度资讯及不建立局部近似模型的情形下（和其他不使用梯度资讯的最佳化演算法相反），用最佳化演算法求得罗森布罗克函数的最小值。以下的例子说明如何用适应坐标下降法对二维的罗森布罗克函数进行最佳化，启始点$x_{0}=(-3,-4)$。在325次函数的运算后可找到最小值的位置，函数的数值为$10^{-10}$。

相关条目

• Himmelblau函数
• Rastrigin函数
• 格里旺克函数

外部链接

• Rosenbrock function plot in 3D（页面存档备份，存于互联网档案馆）
• Minimizing the Rosenbrock Function（页面存档备份，存于互联网档案馆） by Michael Croucher, The Wolfram Demonstrations Project.
• 埃里克·韦斯坦因. Rosenbrock Function. MathWorld.