一、无约束优化问题
- 首先考虑一个不带任何约束的优化问题:minxf(x)\min_xf(x)xminf(x) 其中x∈Rdx∈R^dx∈Rd
- 根据Fermat定理,直接求解梯度等于0的方程:∇xf(x)=0∇_xf(x)=0∇xf(x)=0
- 这个方程的解叫做函数f(x)f(x)f(x)的无条件极值,可能有多个,最小的就是f(x)f(x)f(x)的最小值。
- 下面的推导中,我们把上面方程的解称为:f(x)f(x)f(x)的全局最优解
- 如果没有解析解,可以使用梯度下降或牛顿方法等迭代的手段来使 x 沿负梯度方向逐步逼近极小值点。
二、等式约束优化问题
- 问题描述:minxf(x)\min_xf(x)xminf(x)s.t.g(x)=0s.t. \quad g(x)=0s.t.g(x)=0其中x∈R^d.
- 二维平面如下图所示,约束条件在红线上:
- 从几何角度看,该问题的目标是在由方程g(x)=0g(x)=0g(x)=0确定的d-1维曲面上寻找能使目标函数f(x)f(x)f(x)最小化的点.
- 约束条件会将解的范围限定在一个可行域,此时∇f(x)=0∇f(x)=0∇f(x)=0的解不一定在这个可行域内,我们只需找到在可行域内使得f(x)f(x)f(x)最小的值即可,即:
- 如果∇f(x∗)=0∇f(x^* )=0∇f(x∗)=0的解x∗x^*x∗在约束区域内,那么x∗x^*x∗就是最优解 ,即全局最优解一定是约束条件下的最优解。
- 如果约束条件下,∇f(x)!=0∇f(x)!=0∇f(x)!=0 ,即全局最优解没用落在约束区域内,那么就在约束区间中寻找最接近全局最优解的点作为束条件下的最优解。
- 为了更好的理解拉格朗日乘子法推导,这里先讨论一些梯度的含义:
- 假设二维条件下的约束曲线为:g(x,y)=0g(x,y)=0g(x,y)=0,如上图的红线。其梯度为:∇g(x)=[∂g(x,y)∂x,∂g(x,y)∂y]∇g(x)=[\frac{∂g(x,y)}{∂x},\frac{∂g(x,y)}{∂y}]∇g(x)=[∂x∂g(x,y),∂y∂g(x,y)]其含义为g(x,y)=0g(x,y)=0g(x,y)=0的法向量,正交于约束曲线,并且其指向g(x,y)>0g(x,y)>0g(x,y)>0的方向。
- 也可以这么理解:z=g(x,y)z=g(x,y)z=g(x,y),为一个二元函数,其图像是三维空间总的有个立体图形,则约束曲面(g(x,y)=0g(x,y)=0g(x,y)=0)为:z=g(x,y)z=g(x,y)z=g(x,y)与xoyxoyxoy坐标面的交线。那么zzz的梯度为:∇z=∇g(x)=[∂z∂x,∂z∂y]=[∂g(x,y)∂x,∂g(x,y)∂y]∇z=∇g(x)=[\frac{∂z}{∂x},\frac{∂z}{∂y}]=[\frac{∂g(x,y)}{∂x},\frac{∂g(x,y)}{∂y}]∇z=∇g(x)=[∂x∂z,∂y∂z]=[∂x∂g(x,y),∂y∂g(x,y)]其含义为函数值zzz增长最快的方向,即沿着∇z∇z∇z改变(x,y)(x,y)(x,y),zzz增长最快。所以在g(x,y)=0g(x,y)=0g(x,y)=0的点沿着梯度方向走,zzz必定变大,所以梯度指向g(x,y)>0g(x,y)>0g(x,y)>0的方向。
- 所以约束条件g(x,y)=0g(x,y)=0g(x,y)=0其实就是目标函数目标函数(f(x)f(x)f(x))定义域中的一个区域。目标函数可以写做: z=f(x)z=f(x)z=f(x),即目标函数图像比约束条件多一个维度。
- 根据上面梯度含义,如果xxx(向量)沿着一个与−∇f(x)-∇f(x)−∇f(x)方向夹角小于90的方向移动f(x)f(x)f(x)一定会减小。
- 根据上面性质我们得出:在最优点x∗x^*x∗,目标函数在该点的梯度∇f(x∗)∇f(x^*)∇f(x∗)正交于约束曲面,即∇g(x∗)//∇f(x∗)∇g(x^*)//∇f(x^*)∇g(x∗)//∇f(x∗)。
- 证明:反证法,如果在最优点x∗x^*x∗,∇g(x∗)∇g(x^*)∇g(x∗)与∇f(x∗)∇f(x^*)∇f(x∗)不平行,那么g(x)g(x)g(x)在x∗x^*x∗点的切线方向(记作Δg(x∗)Δg(x^*)Δg(x∗))与∇f(x∗)∇f(x^*)∇f(x∗)的夹角不等于90度,
- 当夹角小于90度,那么−Δg(x∗)-Δg(x^*)−Δg(x∗)与−∇f(x∗)-∇f(x^*)−∇f(x∗)夹角小于90度,那只要我们沿着切线反方向走f(x)f(x)f(x)函数值会下降。
- 当夹角大于90度,那么Δg(x∗)Δg(x^*)Δg(x∗)与−∇f(x∗)-∇f(x^*)−∇f(x∗)夹角小于90度,那只要我们沿着切线方向走f(x)f(x)f(x)函数值会下降。
- 总结:只要∇f(x∗)∇f(x^*)∇f(x∗)与Δg(x∗)Δg(x^*)Δg(x∗)不正交,沿着切线方向或方向f(x)f(x)f(x)函数值都会下降。
- 证明:反证法,如果在最优点x∗x^*x∗,∇g(x∗)∇g(x^*)∇g(x∗)与∇f(x∗)∇f(x^*)∇f(x∗)不平行,那么g(x)g(x)g(x)在x∗x^*x∗点的切线方向(记作Δg(x∗)Δg(x^*)Δg(x∗))与∇f(x∗)∇f(x^*)∇f(x∗)的夹角不等于90度,
- 由此可知,在最优点x∗x^*x∗,梯度∇g(x∗)∇g(x^*)∇g(x∗)和∇f(x∗)∇f(x^*)∇f(x∗)的方向必相同或相反,即存在常数λλλ使得:∇f(x∗)+λ∇g(x∗)=0∇f(x^* )+λ∇g(x^* )=0∇f(x∗)+λ∇g(x∗)=0λ=0λ=0λ=0是为∇f(x∗)=0∇f(x^* )=0∇f(x∗)=0的情况,即全局最优解在约束区域中。
- λλλ称为拉格朗日乘子,定义拉格朗日函数:L(x,λ)=f(x)+λg(x)L(x,λ)=f(x)+λg(x)L(x,λ)=f(x)+λg(x)
- 求拉格朗日函数的无条件极值,即求下面方程组: {∂L(x,λ)∂x=∇f(x)+λ∇g(x)=0∂L(x,λ)∂w=g(x)=0\begin{cases} \frac{∂L(x,λ)}{∂x}=∇f(x)+λ∇g(x)=0\\\\ \frac{∂L(x,λ)}{∂w}=g(x)=0 \end{cases}⎩⎪⎨⎪⎧∂x∂L(x,λ)=∇f(x)+λ∇g(x)=0∂w∂L(x,λ)=g(x)=0
- 由此可以看出原约束优化问题可转化为对拉格朗日函数L(x,λ)L(x,λ)L(x,λ)的无约束优化问题.
三、不等式约束问题
- 问题描述:minxf(x)\min_xf(x)xminf(x)s.t.g(x)≤0s.t. \quad g(x)≤0s.t.g(x)≤0其中x∈R^d.
- 二维情况,如下图:
- 此时有两种情况:(这里和其他人写的可能不一样,如果有错误欢迎指正)
- 第一种情况:就是目标函数f(x)f(x)f(x)的无条件极值(全局最优解)落在约束区域(g(x)≤0g(x)≤0g(x)≤0)内,即下面公式有解:{∇f(x)=0g(x)≤0\begin{cases} ∇f(x)=0 \\ g(x)≤0 \end{cases}{∇f(x)=0g(x)≤0那么这个解就是我们要求的约束区域上的最优解。
- 第二种情况:就是目标函数f(x)f(x)f(x)的无条件极值落在约束区域外面,即在整个约束区域内∇f(x)!=0∇f(x) != 0∇f(x)!=0:
- 这时约束条件的下最优解应尽量靠近无约束时的全局最优解,所以这种情况的约束条件下的最优解一定在边界上。
- 那么怎样才算个最靠近全局最优解呢?假设x∗x^*x∗为这个要求的约束条件下的最优解,那么必定有:∇f(x∗)∇f(x^* )∇f(x∗)的方向必与∇g(x∗)∇g(x^* )∇g(x∗)相反,即存在常数λ>0λ>0λ>0使得:∇f(x∗)+λ∇g(x∗)=0∇f(x^* )+λ∇g(x^* )=0∇f(x∗)+λ∇g(x∗)=0为啥λ>0λ>0λ>0呢?我们先看一张具体的图:图中红色向量为∇g(x)∇g(x)∇g(x),蓝色向量为−∇f(x)-∇f(x)−∇f(x),因为−∇f(x)-∇f(x)−∇f(x)为目标函数减少最快的方向,那么它一定指向全局最优解,所以只有与−∇f(x)-∇f(x)−∇f(x)同∇g(x)∇g(x)∇g(x)向时,无论沿着边界往那个方向移动走f(x)f(x)f(x)都不会变小,f(x)f(x)f(x)要减小只有走出约束区域。
- 两种情况的图:
- 现在用拉格朗日乘子法将上面两种情况统一起来:
- 首先定义拉格朗日函数: L(x,λ)=f(x)+λg(x)L(x,λ)=f(x)+λg(x)L(x,λ)=f(x)+λg(x)
- 将两种情况统一如下:{∂L(x,λ)∂x=∇f(x)+λ∇g(x)=0g(x)≤0λ≥0g(x)λ=0 \begin{cases} \frac{∂L(x,λ)}{∂x}=∇f(x)+λ∇g(x)=0 \\ g(x)≤0 \\ λ≥0\\ g(x)λ=0 \end{cases}⎩⎪⎪⎪⎨⎪⎪⎪⎧∂x∂L(x,λ)=∇f(x)+λ∇g(x)=0g(x)≤0λ≥0g(x)λ=0
- 分析:
- λ=0λ=0λ=0时:对应第一种情况: {∇f(x)=0g(x)≤0 \begin{cases} ∇f(x)=0 \\ g(x)≤0 \end{cases}{∇f(x)=0g(x)≤0
- λ>0λ>0λ>0时,则根据最后一个约束条件可知g(x)=0g(x)=0g(x)=0,对应第二种情况:{∇f(x)+λ∇g(x)=0g(x)=0λ>0 \begin{cases} ∇f(x)+λ∇g(x)=0 \\ g(x)=0\\ λ>0 \end{cases}⎩⎪⎨⎪⎧∇f(x)+λ∇g(x)=0g(x)=0λ>0
- 一般很少求解下面方程来求最优解:{∂L(x,λ)∂x=∇f(x)+λ∇g(x)=0g(x)≤0λ≥0g(x)λ=0 \begin{cases} \frac{∂L(x,λ)}{∂x}=∇f(x)+λ∇g(x)=0 \\ g(x)≤0 \\ λ≥0\\ g(x)λ=0 \end{cases}⎩⎪⎪⎪⎨⎪⎪⎪⎧∂x∂L(x,λ)=∇f(x)+λ∇g(x)=0g(x)≤0λ≥0g(x)λ=0一般作为判断最优解的条件,即若x∗x^*x∗为最优解,那么必定满足上面方程,该方程组也叫做KKT条件。
四、多个约束问题
- 问题描述:minxf(x)\min_xf(x)xminf(x)s.t.hi(x)=0,i=1,2…,mgj(x)≤0,j=1,2….,ns.t. \quad h_i (x)=0,i=1,2…,m\\ g_j (x)≤0,j=1,2….,n s.t.hi(x)=0,i=1,2…,mgj(x)≤0,j=1,2….,n其中x∈R^d.
- 引入拉格朗日乘子λ=(λ1,λ2,…,λm)λ=(λ_1,λ_2,…,λ_m)λ=(λ1,λ2,…,λm)和μ=(μ1,μ2,…,μn)μ=(μ_1,μ_2,…,μ_n)μ=(μ1,μ2,…,μn),相应的拉格朗日函数为:L(x,λ,μ)=f(x)+∑i=1mλihi(x)+∑j=1nμjgj(x)L(x,λ,μ)=f(x)+∑_{i=1}^m λ_i h_i (x) +∑_{j=1}^n μ_j g_j (x)L(x,λ,μ)=f(x)+i=1∑mλihi(x)+j=1∑nμjgj(x)KKT条件为:{∂L(x,λ,μ)∂x=∇f(x)+∑i=1mλi∇hi(x)+∑j=1nμj∇gj(x)=0gj(x)≤0μj≥0gj(x)μj=0 \begin{cases} \frac{∂L(x,λ,μ)}{∂x}=∇f(x)+∑_{i=1}^m λ_i ∇h_i (x) +∑_{j=1}^n μ_j ∇g_j (x)=0 \\ g_j(x)≤0 \\ μ_j≥0\\ g_j(x)μ_j=0 \end{cases}⎩⎪⎪⎪⎨⎪⎪⎪⎧∂x∂L(x,λ,μ)=∇f(x)+∑i=1mλi∇hi(x)+∑j=1nμj∇gj(x)=0gj(x)≤0μj≥0gj(x)μj=0