拉格朗日乘数和KTT条件

拉格朗日乘数和KTT条件

通常来说最优化问题是指:对于给定的某一函数,求其在指定作用域上的全局最优解。在求解数学的最优化问题中,Lagrange Multiplier (拉格朗日乘子法)和 Karush-Kuhn-Tucker Conditions(KTT条件)是两种常用于寻找多元函数在其变量受到一个或多个条件的约束时的极值的方法。

拉格朗日乘子法常用于等式约束的优化问题,求取出最优值;KKT条件常用于受不等式约束的优化问题。其中KKT条件是拉格朗日乘数的推广。

常见的最优化问题有三种:

  1. 无约束优化问题

min ⁡ f ( x ) \min f\left( x \right) minf(x)

  1. 等式约束问题 ( s . t . s.t. s.t. 表示subject to,受限于)

min ⁡ f ( x ) \min f\left( x \right) minf(x)
s . t .        h i ( x ) = 0                i = 1 , 2 , … , l s.t.{\rm{ }}\;\;\;{h_i}\left( x \right) = 0\;\;\;\;\;\;\;i = 1,2, \ldots ,l s.t.hi(x)=0i=1,2,,l

  1. 不等式约束问题

min ⁡ f ( x ) \min f\left( x \right) minf(x)
s . t .        h i ( x ) = 0                i = 1 , 2 , … , l                  g j ( x ) ≤ 0            j = 1 , 2 , … , m \begin{array}{l} s.t.{\rm{ }}\;\;\;{h_i}\left( x \right) = 0\;\;\;\;\;\;\;i = 1,2, \ldots ,l\\ \;\;\;\;\;\;\;\;{g_j}\left( x \right) \le 0\;\;\;\;\;j = 1,2, \ldots ,m \end{array} s.t.hi(x)=0i=1,2,,lgj(x)0j=1,2,,m

对于第一类优化问题,解决方法通常是函数对变量求导,令求导函数等于0的点可能是极值点,最后将结果带回原函数进行验证。

拉格朗日乘数——求解等式约束问题

首先我们通过维基百科的一个二维实例来进行介绍。

设目标函数为 f ( x , y ) f(x,y) f(x,y) ,约束等式为 g ( x , y ) = c g(x,y)=c g(x,y)=c

min ⁡ f ( x , y ) \min f\left( {x,y} \right) minf(x,y)
s . t .        g ( x , y ) = c        ( c    i s    c o n s t a n t ) s.t.\;\;\;g\left( {x,y} \right) = c\;\;\;\left( {c\;is\;constant} \right) s.t.g(x,y)=c(cisconstant)

我们可以想象一下, f ( x ) f\left( x \right) f(x) 在三维空间中是一个曲面,当它与一个面 z = d n z={{d}_{n}} z=dn 相交构成等高线。对于不同的 d n d_n dn 值, f ( x ) f\left( x \right) f(x) 有不同的等高线。

想象一下我们沿着约束条件 g ( x , y ) = c g\left( x,y \right)=c g(x,y)=c 构成的可行集走,在大部分情况下, f f f 的等高线和 g g g 的可行集不会重合,但是如果有解的情况下,两条线会有交点。想像此时我们移动 g ( x , y ) = c g\left( x,y \right)=c g(x,y)=c 上的点,因为 f f f 是连续的方程,我们因此能走到 f ( x , y ) = d n f\left( x,y \right)={{d}_{n}} f(x,y)=dn 更高或更低的等高线上,也就是说 d n d_n dn 可以变大或变小。只有当 g ( x , y ) = c g\left( x,y \right)=c g(x,y)=c f ( x , y ) = d n f\left( x,y \right)={{d}_{n}} f(x,y)=dn 相切,也就是说,此时,我们正同时沿着 g ( x , y ) = c g\left( x,y \right)=c g(x,y)=c f ( x , y ) = d n f\left( x,y \right)={{d}_{n}} f(x,y)=dn 走。这种情况下,会出现极值。

用向量的形式来表达的话,我们说相切的性质在此意味着 f f f g g g 的切线在某点上平行,同时也意味着两者的梯度平行。此时引入一个未知标量 λ λ λ ,并求解:

∇ [ f ( x , y ) + λ ( g ( x , y ) − c ) ] = 0 \nabla \left[ f\left( x,y \right)+\lambda \left( g\left( x,y \right)-c \right) \right]=0 [f(x,y)+λ(g(x,y)c)]=0

λ ≠ 0 \lambda \ne 0 λ=0 .
一旦求出 λ λ λ 的值,将其套入下式,易求在无约束条件下的极值和对应的极值点。

F ( x , y , λ ) = f ( x , y ) + λ ( g ( x , y ) − c ) F\left( x,y,\lambda \right)=f\left( x,y \right)+\lambda \left( g\left( x,y \right)-c \right) F(x,y,λ)=f(x,y)+λ(g(x,y)c)

新方程 F ( x , y , λ ) F\left( x,y,\lambda \right) F(x,y,λ)在达到极值时与 f ( x , y ) f\left( x,y \right) f(x,y) 相等,因为 F ( x , y , λ ) F\left( x,y,\lambda \right) F(x,y,λ) 达到极值时 g ( x , y ) − c g\left( x,y \right)-c g(x,y)c 总等于零。

上图蓝线为 $f\left( x,y \right)={{d}_{n}} $ 的等高线,绿线为 $g\left( x,y \right)=c $ 的点的轨迹。从梯度的方向上来看,显然有 d 1 > d 2 {{d}_{1}}>{{d}_{2}} d1>d2 。绿色的线是约束,也就是说,只要正好落在这条绿线上的点才可能是满足要求的点。如果没有这条约束, f ( x , y ) f\left( x,y \right) f(x,y) 的最小值应该会落在最小那圈等高线内部的某一点上。而现在加上了约束,最小值点应该在哪里呢?显然应该是在 f ( x , y ) f\left( x,y \right) f(x,y) 的等高线正好和约束线相切的位置,因为如果只是相交意味着肯定还存在其它的等高线在该条等高线的内部或者外部,使得新的等高线与目标函数的交点的值更大或者更小,只有到等高线与目标函数的曲线相切的时候,可能取得最优值。

KTT条件——求解不等式约束问题

设目标函数为 f ( x , y ) f(x,y) f(x,y) ,约束等式为 h i ( x ) {h_i}(x) hi(x),约束不等式为 g j ( x ) {g_j}(x) gj(x)

min ⁡ f ( x ) \min f\left( x \right) minf(x)
s . t .        h i ( x ) = 0                i = 1 , 2 , … , l                  g j ( x ) ≤ 0            j = 1 , 2 , … , m \begin{array}{l} s.t.{\rm{ }}\;\;\;{h_i}\left( x \right) = 0\;\;\;\;\;\;\;i = 1,2, \ldots ,l\\ \;\;\;\;\;\;\;\;{g_j}\left( x \right) \le 0\;\;\;\;\;j = 1,2, \ldots ,m \end{array} s.t.hi(x)=0i=1,2,,lgj(x)0j=1,2,,m

则我们定义不等式约束下的拉格朗日函数 F F F ,则 F F F 的表达式为

F ( x , λ , μ ) = f ( x ) + ∑ i = 1 l λ i h i ( x ) + ∑ j = 1 m μ j g j ( x ) F(x,\lambda ,\mu ) = f(x) + \sum\limits_{i = 1}^l {{\lambda _i}} {h_i}(x) + \sum\limits_{j = 1}^m {{\mu _j}} {g_j}(x) F(x,λ,μ)=f(x)+i=1lλihi(x)+j=1mμjgj(x)

未完待续。。。

参看文献

  • 2
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值