看了不少大佬的博客,对拉格朗日乘子法和kkt条件有一点点认识,记录一下
拉格朗日乘子法
数学上我们经常遇到很多地方求极值问题,在没有约束条件下,相对比较容易。但有约束条件下就不太好解决。
(1)无约束
minf(x)minf(x)minf(x)
(2)等式约束
minf(x)minf(x)minf(x) s.t.g(x)==0s.t. g(x)==0s.t.g(x)==0
此类问题,就可以直接应用拉格朗日乘子法,引入系数变量λ,构造新的函数F(x,λ)=f(x)+∑i=1kλigi(x)F(x,λ)=f(x)+\sum_{i=1}^{k}λ_{i}g_{i}(x)F(x,λ)=f(x)+i=1∑kλigi(x)
对各个系数求导,并且令导数为0,∂F∂x=f′(x)+∑i=1kλigi′(x)=0\frac{\partial F}{\partial x}=f^{'} (x)+\sum_{i=1}^{k}λ_{i}g_{i}^{'}(x)=0∂x∂F=f′(x)+i=1∑kλigi′(x)=0
∂F∂λi=gi(x)=0\frac{\partial F}{\partial \lambda_{i}}=g_{i}(x)=0∂λi∂F=gi(x)=0
可以看到极值的时候要满足的条件是,gi(x)=0,原来的约束条件体现在这,并且F(x)取极值的时候与f(x)相等,F(x)的极值也等于f(x)的极值。因此拉格朗日乘子法可以将等式约束问题转成无约束问题。
看经典的图来直观解释一下,假设在三维坐标上,令z=f(x,y),我们要求z的极值,约束条件是g(x,y)=0,在xoy平面上,将函数f(x,y)的等高线投影到xoy平面上,交点满足约束条件,并且是f(x,y)的一个值,如果相交的话,一定存在更大或者更小的等高线,那么在相切的位置才可能取到极值。
(3) 不等式约束
minf(x)minf(x)minf(x) s.t.g(x)==0,h(x)<=0s.t. g(x)==0,h(x)<=0s.t.g(x)==0,h(x)<=0
构造函数F(x,λ)=f(x)+∑i=1kλigi(x)+∑i=1kμihi(x)F(x,λ)=f(x)+\sum_{i=1}^{k}λ_{i}g_{i}(x)+\sum_{i=1}^{k}\mu_{i}h_{i}(x)F(x,λ)=f(x)+i=1∑kλigi(x)+i=1∑kμihi(x)
kkt条件:
- μihi(x)=0\mu_{i}h_{i}(x)=0μihi(x)=0
- μi>=0\mu_{i}>=0μi>=0
- hi(x)<=0h_{i}(x)<=0hi(x)<=0
- f′(x)+∑i=1kλigi′(x)+∑i=1kμihi′(x)=0f^{'} (x)+\sum_{i=1}^{k}λ_{i}g_{i}^{'}(x)+\sum_{i=1}^{k}\mu_{i}h_{i}^{'}(x)=0f′(x)+i=1∑kλigi′(x)+i=1∑kμihi′(x)=0
要满足kkt条件,然后应用拉格朗日乘子法求解。
谈谈我对这些条件的认识:
前面已经讲到过在等式约束的时候,F取极值的时候一定也是f的极值,那么2,3条件可以得到μigi(x)<=0\mu_{i}g_{i}(x)<=0μigi(x)<=0
可以得到maxF(x)=f(x),只有μihi(x)=0,才能取到极小值,那么minxf(x)=minxmaxμ,λF(x)maxF(x)=f(x),只有\mu_{i}h_{i}(x)=0,才能取到极小值,那么min_{x}f(x)=min_{x} max_{\mu , \lambda}F(x)maxF(x)=f(x),只有μihi(x)=0,才能取到极小值,那么minxf(x)=minxmaxμ,λF(x)
转化为对偶问题,minxf(x)=maxμ,λminxF(x),假设x0处取得极小值.转化为对偶问题,min_{x}f(x)=max_{\mu , \lambda} min_{x}F(x),假设x_{0}处取得极小值.转化为对偶问题,minxf(x)=maxμ,λminxF(x),假设x0处取得极小值.
minxf(x)=f(x0)=maxμ,λ(f(x0)+∑i=1kλigi(x0)+∑i=1kμihi(x0))=f(x0)min_{x}f(x)=f(x_{0})=max_{\mu , \lambda} (f(x_{0})+\sum_{i=1}^{k}λ_{i}g_{i}(x_{0})+\sum_{i=1}^{k}\mu_{i}h_{i}(x_{0}))=f(x_{0})minxf(x)=f(x0)=maxμ,λ(f(x0)+i=1∑kλigi(x0)+i=1∑kμihi(x0))=f(x0)
再来看条件1:
满足条件1的时候:
- μi=0,此时约束条件不起作用,不出现在函数F中,极值应该是落在可行域内\mu_{i}=0,此时约束条件不起作用,不出现在函数F中,极值应该是落在可行域内μi=0,此时约束条件不起作用,不出现在函数F中,极值应该是落在可行域内
- gi(x)=0,此时约束条件起作用,极值应该是落在gi(x)=0上,不等式约束变成等式约束g_{i}(x)=0,此时约束条件起作用,极值应该是落在g_{i}(x)=0上,不等式约束变成等式约束gi(x)=0,此时约束条件起作用,极值应该是落在gi(x)=0上,不等式约束变成等式约束
看图理解一下
如有错误请麻烦指正