漫步最优化十四——凸函数与凹函数




 




——

通常在实际中,最小化的函数有几个极值,所以最优化算法得出的极值不确实是否为全局的极值,对于一些特殊的函数,凸函数与凹函数,任何局部极值也是全局极致,因此如果目标函数是凸的或凹的,那么优化算法就能保证是全局的。

1集合RcEn是凸集,如果对每对点x1,x2Rc,每个实数α,0<α<1,点

x=αx1+(1α)x2

位于Rc,即xRc

效果上,如果任何两点x1,x2Rc用直线相连,x1,x2之间线上的每个点都在Rc中,那么Rc是凸的。如果存在点不在Rc中,那么该集合是非凸的,凸集合如图1所示。

凸的概念也可以用到函数上。

2

  1. 我们称定义在凸集Rc上的函数f(x)为凸的,如果对每对x1,x2Rc与每个实数α,0<α<1,不等式
    f[αx1+(1α)x2]αf(x1)+(1α)f(x2)

满足。如果x1x2

f[αx1+(1α)x2]<αf(x1)+(1α)f(x2)

满足,那么f(x)是严格凸的。

  1. 如果φ(x)定义在凸集Rc上且f(x)=φ(x)是凸的,那么φ(x)是凹的。如果f(x)是严格凸的,那么φ(x)是严格凹的。

上述定义中的不等式,左边是点x1,x2之间某处的f(x)值,而右边是基于线性插值得到的f(x)的近似,因此如果任何两点的线性插值大于函数的值,那么该函数就是凸的,图2a,b中的函数为凸的,2c为非凸的。


漫步最优化十四——凸函数与凹函数
图1

1如果

f(x)=af1(x)+bf2(x)

其中a,b0,f1(x),f2(x)是凸集Rc上的凸函数,那么f(x)是集合Rc上的凸函数。

因为f1(x),f2(x)是凸函数,a,b0,所以对于x=αx1+(1α)x2,我们有

af1(αx1+(1α)x2)a[αf1((x)1)+(1α)f1(x2)]bf2(αx1+(1α)x2)b[αf2((x)1)+(1α)f2(x2)]

其中0<α<1,因此

f(x)f(αx1+(1α)x2)=af1(x)+bf2(x)=af1(αx1+(1α)x2)+bf2(αx1+(1α)x2)α[af1(x1)+bf2(x1)]+(1α)[af1(x2)+bf2(x2)]

因为

af1(x1)+bf2(x1)=f(x1)af1(x2)+bf2(x2)=f(x2)

所以上面的不等式可以写成

f(αx1+(1α)x2)αf(x1)+(1α)f(x2)

f(x)是凸函数。


漫步最优化十四——凸函数与凹函数
图2

2如果f(x)是凸集Rc上的凸函数,那么对每个实数K而言,集合

Sc={x:xRc,f(x)K}

都是凸集。

如果x1,x2Sc,那么根据Sc的定义,f(x1)K,f(x2)K,因为f(x)是凸集,所以

f[αx1+(1α)x2]αf(x1)+(1α)f(x2)αK+(1α)K

或者

f(x)Kfor x=αx1+(1α)x2and0<α<1

所以

xSc

Sc是凸的。

定理2的图示如图3,其中集合Sc是凸集,如果f(x)在凸集Rc上是凸函数的话。


漫步最优化十四——凸函数与凹函数
图3

另一种考虑凸的角度是测试f(x)的梯度与海森矩阵。

3如果f(x)C1,那么f(x)在凸集Rc上是凸函数,当且仅当对所有x,x1Rc

f(x1)f(x)+g(x)T(x1x)

其中g(x)f(x)的梯度。

这个定理的证明由两部分组成。首先我们证明如果f(x)是凸函数,那么不等式成立。然后证明如果不等式成立,那么f(x)是凸函数。首先如果f(x)是凸函数,那么对于所有α,0<α<1

f[αx1+(1α)x]αf(x1)+(1α)f(x)

或者

f[x+α(x1x)]f(x)α[f(x1)f(x)]

α0,由f[x+α(x1x)]的泰勒级数可得

f(x)+g(x)Tα(x1x)f(x)α[f(x1)f(x)]

所以

f(x1)f(x)+g(x)T(x1x)

接下来,如果不等式在x,x2Rc处成立,那么

f(x2)f(x)+g(x)T(x2x)

从上面两式可得

αf(x1)+(1α)f(x2)αf(x)+αg(x)T(x1x)+(1α)f(x)+(1α)g(x)T(x2x)

或者

αf(x1)+(1α)f(x2)f(x)+gT(x)[αx1+(1α)x2x]

代入

x=αx1+(1α)x2

可得

f[αx1+(1α)x2]αf(x1)+(1α)f(x2)

其中0<α<1,因此f(x)是凸函数。

定理3说明f(x)在点x处基于f(x)导数的线性插值小于函数值,如图4所示。

4函数f(x)C2是凸集Rc上的凸函数,当且仅当f(x)的海森矩阵H(x)xRc是半正定的。

如果x1=x+d,其中x1,xRc中的任意点,那么由泰勒级数可得

f(x1)=f(x)+gT(x)(x1x)+12dTH(x+αd)d

其中0α1,接下来如果H(x)Rc中是半正定的,那么

12dTH(x+αd)d0

所以

f(x1)f(x)+gT(x)(x1x)

所以由定理3可知f(x)是凸函数。

如果H(x)Rc任何处都是半正定的,那么存在点x与方向d使得

dTH(x+αd)<0

所以

f(x1)<f(x)+gT(x)(x1x)

根据定理3可知f(x)是非凸的,所以当且仅当H(x)Rc任何地方是半正定时f(x)是凸函数。


漫步最优化十四——凸函数与凹函数
图4

对于严格凸函数,上面的定理修改如下:

5

  • 如果f(x)是凸集Rc上的严格凸函数,那么对每个实数K而言,集合
    Sc={x:xRc,f(x)<K}

都是凸集。

  • 如果f(x)C1,那么f(x)在凸集Rc上的严格凸函数,当且仅当对所有x,x1Rc
    f(x1)>f(x)+g(x)T(x1x)

其中g(x)f(x)的梯度。

  • 函数f(x)C2是凸集Rc上的凸函数,当且仅当f(x)的海森矩阵H(x)xRc是正定的。

如果φ(x)定义在凸集Rc上,且f(x)=φ(x)是严格凸函数,那么φ(x)是严格凹函数且其海森矩阵是负定的。反过来,如果φ(x)的海森矩阵是负定的,那么φ(x)是严格凹的。