最优化

约束最优化方法 (一) 最优性条件

2019-12-07  本文已影响0人  小小何先生

  之前讨论的是无约束最优化方法,这一节主要介绍的是带有约束的非线性规划问题,所谓的非线性规划,就是约束项含有平方这种。解这类问题有两种方法,一个是容许方向法、它是一种直接处理约束的方法;另一个是罚函数法,它是将约束问题转变成一系列无约束问题,用无约束的极小点去逐渐逼近约束问题的极小点。但是在介绍这两种方法之前,要先介绍一些概念。

最优性条件

  本节主要讨论一般约束问题的最优性条件。我们将先从仅含等式约束或不等式约束的问题入手,然后自然过渡到一般约束问题。所以这一节主要介绍各种约束下的最优性条件,也就是各种约束下,什么样的条件能够推出这个点是最优点、另外一种,已知各种约束下的最优点,能够推出什么条件。整体目录结构如下:

等式约束问题的最优性条件

  考虑仅含等式约束的问题1
\begin{array}{l}{\min f(x)} \\ {\text {s.t. } h_{j}(x)=0, \quad j=1,2, \cdots, l .}\end{array}

  这个问题的最优性条件与求解方法在微积分中已从理论上得到了解决,这就是Lagrange定理和Lagrange乘子法。

定理1:假设

  这个定理的意义还在于,它把对等式约束问题的求解转化为对无约束问题的求解。

  上式是最优性一阶必要条件

定理2:在约束问题1中,假设:

v^{T} \bigtriangledown h_{j}(x^{*})=0,j=1,2,···,l

  的任意非零向量v \in R^{n},都有:
v^{T} \bigtriangledown_{x}^{2}L(x^{*},\lambda^{*})v > 0

  这个定理的几何意义是,在Lagrange函数的驻点[x^{*} \ \ \lambda^{*}]处,如果Lagrange函数关于x的Hesse矩阵在l个约束(超)曲面的切平面的交集上正定(注意,并不需要在原来的空间中正定),那么x^{*}就是严格局部极小点。

  这里就是直接给出两个定理,没办法,理解记忆吧。第一个定理相对来说比较重要一点。

不等式约束问题的最优性条件

几何最优性条件

  下面将给出约束问题2
\begin{array}{l}{\min f(x)} \\ {\text {s.t. } s_{i}(x) \geqslant 0, \quad i=1,2, \cdots, m .}\end{array}
  的最优性条件。

定义1 :对于约束问题,设\overline{x} \in DD=\{x|s_{i}(x)\geq 0 , i =1,2, \cdots, m \}\overline{x} 使得某个不等式约束有s_{i}(\overline{x})=0,则该不等式约束s_{i}(x) \geq 0称为是关于容许点\overline{x}起作用约束;否者,若s_{i}(\overline{x}) > 0,则该不等式约束称为是关于容许点的不起作用约束

定义2 :设CR中的非空集,且x \in C。对于\forall p \in R^{n},若当x+p \in C时,对于\forall t \geq 0,必有x+tp \in C,则集合C称为x为顶点的锥。若锥C是凸集,则称为凸锥

定义3 :设DR^{n}中的非空集,且x \in D。对于非零向量p \in R^{n},若存在\delta > 0,当t \in (0, \delta)时,必有x+tp \in D,则p称为点x容许方向向量,其方向称为点x容许方向。由点x的全部容许方向向量构成的集合称为点x容许方向集,或者说容许方向锥

如图所示
引理 :设,;并设当时,在点处可微,当时,在点处连续。若对于所有的,向量都使得,则是点的一个容许方向。

  约束曲面s_{i}(x)=0把整个空间分成两部分,梯度\bigtriangledown_{s_{i}}(\widetilde{x})总是指向包含容许集的那一侧。

  把由点x的所有下降方向向量构成的集合称为点x下降方向锥

定理:设f:R^{n} \rightarrow R^{1}在点x可微,则点x处的下降方向向量p必满足:
\bigtriangledown f(x)^{T}p < 0

  记S(x)=\{p|\bigtriangledown f(x)^{T}p < 0 \},则S(x)是点x处的下降方向集。显然S(x)R^{n}中的半空间。

  接下来就是几何最优性条件的定义:(因为这个条件是仅借助点集的概念给出的,所以称为几何最优性条件):

定理: 在约束问题2中,若x^{*}是局部最优点,则点x^{*}处的容许方向锥下降方向锥的交集是空集。

  上面这个定理仅仅是必要的,而不是充分的。也就是说知道这个点是最优点能够推断出容许方向锥和下降方向锥的交集是空集,但由容许方向锥和下降方向锥的交集是空集并不能推断出其是最优点。

Fritz John条件

  这里要介绍:引理(Farkas)、引理(Gordan)、定理:Fritz John

引理(Farkas):设a_{1}a_{2}\cdotsa_{m}bn维向量,则满足:
a_{i}^{T}p \geq 0, \ i=1,2,\cdots , m
  的向量p也满足
b^{T}p \geq 0
  的充要条件是,存在非负数\gamma_{1}\gamma_{2}\cdots\gamma_{n},使得:
b=\sum_{i=1}^{m}\gamma_{i}a_{i}

  这个依旧不需要证明,相信它就完事了,因为直观上感觉就是非常正确的。可以看课本图4-6。或者下面这张图理解(b理解为f(x^{*})):

image
引理(Gordan):设,,,是维向量,则不存在向量使得:

  成立的充要条件是,存在不全为零的非负数,,,,使得:

  这个怎么理解呢?不存在向量使得,所表示的几何意义就是,,,不会在超平面的一侧,因为要是在一侧的话,就会存在这样一个向量满足要求。既然不会在一侧,那么就一定会有一个非零的线性组合,使其最终结果为0。

定理:Fritz John: 在问题2中,设x^{*}是局部最优解,f(x)s_{1}(x)s_{2}(x)\cdotss_{m}(x)在点x^{*}处可微,那么,存在不全为零的实数\mu_{0}\mu_{1}\cdots\mu_{m}使得:
\left.\begin{array}{rl}{\mu_{0} \nabla f\left(x^{*}\right)-\sum_{i=1}^{m} \mu_{i} \nabla s_{i}\left(x^{*}\right)} {=0} \\ {\mu_{i} s_{i}\left(x^{*}\right)=0,} {\ \ \ i=1,2, \cdots, m} \\ {\mu_{i} \geq 0,} {\ \ \ i=0,1, \cdots, m}\end{array}\right\}

  上面这个式子我们来理解一下,因为这个x^{*}是最优点,所以这个容许集和下降方向集是空集。所以不存在向量p,使得:
\nabla f\left(x^{*}\right)^{T} p<0, \\ \left(-\nabla s_{i}\left(x^{*}\right)\right)^{T} p<0, i \in I\left(x^{*}\right)

  \mu_{i}s_{i}(x^{*})=0\ (i=1,2,\cdots ,m)称为互补松弛条件。它表明:若s_{i}(x^{*})>0,即i \overline{\in} I,则必有\mu_{i}=0;若\mu_{i}>0,则必有s_{i}(x^{*})=0,即i \in I

  这个定理给你了问题2的一个最优性必要条件。上式称为Fritz John条件,满足Fritz-John条件的点称为Fritz-John点\mu_{1}\mu_{2}\cdots\mu_{m}也称为Lagrange乘子。

  Fritz-John条件仅是判别某一容许点是否为最优点的必要条件,而不是充分条件。

Kuhn-Tucker条件

  如果Fritz-John条件中\mu_{0}=0时,Fritz-John条件失去实用价值。为了使得其大于0,就有了Kuhn-Tucker条件。

定理:Kuhn-Tucker:

在问题2中,假设
i)x^{*}是局部最优点;
ii) f(x),s_{1}(x),s_{2}(x),···,s_{m}(x)在点x^{*}处可微;
iii) 点x^{*}处的全部起作用约束的梯度线性无关。那么存在实数\mu_{1}、\mu_{2},···,\mu_{m}使得:
\begin{aligned} \nabla f\left(x^{*}\right)-\sum_{i=1}^{m} \mu_{i} \nabla s_{i}\left(x^{*}\right) =&0 \\ \mu_{i} S_{i}\left(x^{*}\right)=0, i=1,2, \cdots, &m, \\ \mu_{i} \geq 0, i=1, \cdots, &m \end{aligned}

  Kuhn-Tucker条件是Fritz John条件的特殊情况。

  Kuhn-Tucker条件有明显的几何意义。在Kuhn-Tucker定理的公式中,删去不起作用约束项(注意,它们的系数是\mu_{i}=0,当i \overline{\in}I,Kuhn-Tucker条件可简写成:存在\mu_{i} \geq 0i \in I)使得:
\nabla f(x^{*})=\sum_{i \in I}\mu_{i}\nabla_{s_{i}}(x^{*})

  此式在几何上表示:若x^{*}是问题地最优解,根据Farkas引理可知,在此点处地目标函数地梯度必处在由起作用约束函数地梯度张成地凸锥之中。

一般约束问题的最优性条件

  之前是不等式约束,现在考虑一般约束问题地最优性条件,既有等式约束还有不等式约束的情况。我们这节就介绍一般约束问题下的Fritz-John定理Kuhn-Tucker定理

Fritz-John定理

  考虑问题:

\left.\begin{array}{ll}{\min } & {f(x)} \\ {\text {s.t.}} & {s(x) \geq 0} \\ {} & {h(x)=0}\end{array}\right\}

  在上述问题中设x^{*}是局部最优点f(x)s_{1}(x), \cdots ,s_{m}(x), h(x), \cdots ,h_{l}(x)在点x^{*}处可微。那么,存在不全为零的实数\mu_{0},\mu_{1},\cdots .\mu_{m}, \lambda_{1},\lambda_{2},\cdots \lambda_{l},使得:
\begin{aligned} \mu_{0} \nabla f\left(x^{*}\right) &-\sum_{i=1}^{m} \mu_{i} \nabla s_{i}\left(x^{*}\right)-\sum_{j=1}^{l} \lambda_{j} \nabla h_{j}\left(x^{*}\right)=0 \\ \mu_{i} S_{i}\left(x^{*}\right) &=0, \quad i=1,2, \cdots, m \\ \mu_{i} & \geq 0, \quad i=0,1, \cdots, m \end{aligned}

  这个定理可以看成是Lagrange定理的结论与Fritz-John定理的结论的合并。相当于多了l个等式约束。

Kuhn-Tucker定理

  考虑问题:

\left.\begin{array}{ll}{\min } & {f(x)} \\ {\text {s.t.}} & {s(x) \geq 0} \\ {} & {h(x)=0}\end{array}\right\}

  假设:
  i) x^{*}是局部最优解;
  ii)f(x),s_{1}(x),\cdots , s_{m}(x),h_{1}(x),\cdots ,h_{l}(x)在点x^{*}处可微;
  iii)点x^{*}处的全部起作用约束的梯度线性无关。

那么存在实数\mu_{0},\mu_{1},\cdots .\mu_{m}, \lambda_{1},\lambda_{2},\cdots \lambda_{l}使得:

\begin{aligned} \nabla f\left(x^{*}\right)-\sum_{i=1}^{m} \mu_{i} \nabla s_{i}\left(x^{*}\right)-\sum_{j=1}^{l} \lambda_{j} \nabla h_{j}\left(x^{*}\right)=0 \\ \mu_{i} S_{i}\left(x^{*}\right)=0, \quad i=1,2, \cdots, m \\ \mu_{i} \geq 0, i=1, \cdots, m& \end{aligned}

定理:
  在以下凸规划问题中:
\left.\begin{array}{ll}{\min } & {f(x)} \\ {\text {s.t.}} & {s(x) \geq 0} \\ {} & {h(x)=0}\end{array}\right\}
  假设f是可微凸函数s_{1},\cdots ,s_{m}是可微凹函数,h_{1},\cdots ,h_{l}是线性函数。若x^{*}是上述问题的Kuhn-Tucker点,则x^{*}就是上述问题的全局最优点。(用Hesson矩阵即可证明是不是凸优化问题)。

我的微信公众号名称:深度学习与先进智能决策
微信公众号ID:MultiAgent1024
公众号介绍:主要研究分享深度学习、机器博弈、强化学习等相关内容!期待您的关注,欢迎一起学习交流进步!

上一篇 下一篇

猜你喜欢

热点阅读