日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問 生活随笔!

生活随笔

當前位置: 首頁 > 编程资源 > 编程问答 >内容正文

编程问答

从拉格朗日乘数法到KKT条件

發布時間:2023/12/2 编程问答 46 豆豆
生活随笔 收集整理的這篇文章主要介紹了 从拉格朗日乘数法到KKT条件 小編覺得挺不錯的,現在分享給大家,幫大家做個參考.

從拉格朗日乘數法到KKT條件

最近看論文遇到了Karush–Kuhn–Tucker (KKT)條件,想搞清楚這是個什么東東,因此就把這個東西認真學習一下并且分享出來,希望對大家有用。學習KKT就不得不先學習一下拉格朗日乘數法,于是不得不重新翻出被記憶塵封的高數~~

1.拉格朗日乘數法

在數學最優問題中,拉格朗日乘數法是一種尋找變量受一個或多個條件所限制的多元函數的極值的方法。這種方法將一個有n 個變量與k 個約束條件的最優化問題轉換為一個有n + k個變量的方程組的極值問題,其變量不受任何約束。這種方法引入了一種新的標量未知數,即拉格朗日乘數:約束方程的梯度(gradient)的線性組合里每個向量的系數。

以二元函數為例:

設給定二元函數

maxz=f(x,y)s.t.φ(x,y)=0max z=f(x,y)\\s.t. φ(x,y)=0maxz=f(x,y)s.t.φ(x,y)=0

為尋找z=f(x,y)在附加條件下的極值點,先做拉格朗日函數

F(x,y,λ)=f(x,y)+λφ(x,y)F(x,y,\lambda)=f(x,y)+\lambda φ(x,y)F(x,y,λ)=f(x,y)+λφ(x,y)

其中λ\lambdaλ為參數。令F(x,y,λ)F(x,y,λ)F(x,y,λ)xxxyyyλλλ的一階偏導數等于零,即

Fx′=?x′(x,y)+λφx′(x,y)=0F'_x=?'_x(x,y)+λφ'_x(x,y)=0 Fx?=?x?(x,y)+λφx?(x,y)=0

Fy′=?y′(x,y)+λφy′(x,y)=0F'_y=?'_y(x,y)+λφ'_y(x,y)=0Fy?=?y?(x,y)+λφy?(x,y)=0

Fλ′=φ(x,y)=0F'_λ=φ(x,y)=0Fλ?=φ(x,y)=0

由上述方程組解出x,yx,yx,yλλλ,如此求得的(x,y)(x,y)(x,y),就是函數z=?(x,y)z=?(x,y)z=?(x,y)在附加條件φ(x,y)=0φ(x,y)=0φ(x,y)=0下的可能極值點。
若這樣的點只有一個,由實際問題可直接確定此即所求的點。

下面這篇博客可以很好的幫助理解:

支持向量機(SVM)課前準備(一)–拉格朗日乘子法 - be·freedom - 博客園 (cnblogs.com)

2.KKT條件

先給出一個KKT條件的實例,我們優化的目標是:

minimizef(x)s.t.ki=0gj≤0,i,j=1,2,3,...minimize f(x)\\ s.t. k_i = 0\\ g_j\leq0,i,j=1,2,3,... minimizef(x)s.t.ki?=0gj?0,i,j=1,2,3,...
其中,KKT條件如下:

?f+∑λiki+ωjgj=0.......(1)ki=0......(2)gj=0......(3)uj≥0......(4)ujgj=0......(5)\nabla f+\sum\lambda_ik_i+\omega_jg_j=0.......(1)\\ k_i=0......(2)\\ g_j=0......(3)\\ u_j\geq0......(4)\\ u_jg_j=0......(5)\\ ?f+λi?ki?+ωj?gj?=0.......(1)ki?=0......(2)gj?=0......(3)uj?0......(4)uj?gj?=0......(5)
公式1、2、3容易理解。

公式4、5,通過一個簡單例子說明:

minimizef(x)s.t.g1(x)=a?x≤0g2(x)=x?b≤0minimize f(x)\\s.t. g_1(x)=a-x\leq0\\g_2(x)=x-b\leq0minimizef(x)s.t.g1?(x)=a?x0g2?(x)=x?b0

gig_igi?添加一個 ≥0 的松弛變量a12,b12a_1^2,b_1^2a12?,b12?。得到
g1(x)=a?x+a12g2(x)=x?b+b12g_1(x)=a-x+a_1^2\\g_2(x)=x-b+b_1^2g1?(x)=a?x+a12?g2?(x)=x?b+b12?

由此,我們將不等式轉化為等式約束,應用拉格朗日乘子法。

拉格朗日方程如下:

L(x1,a1,b1,u1,u2)=f(x)+u1(a?x+a12)+u2(x?b+b12)...(1)L(x_1,a_1,b_1,u_1,u_2)=f(x)+u_1(a-x+a_1^2)+u_2(x-b+b_1^2)...(1)L(x1?,a1?,b1?,u1?,u2?)=f(x)+u1?(a?x+a12?)+u2?(x?b+b12?)...(1)
對方程求偏導如下:

{?F?x=?f?x+u1dg1dx+u2dg2dx=?f?x?u1+u2=0...(2)?F?u1=a?x+a12=g1+a12=0...(3)?F?u2=x?b+a12=g2+b12=0...(4)?F?a1=2u1a1=0...(5)?F?b1=2u2b1=0...(6),u1≥0,u2≥0\begin{cases}\frac {\partial F}{\partial x} =\frac{\partial f}{\partial x}+u_1\frac{ozvdkddzhkzdg_1}{ozvdkddzhkzdx}+u_2\frac {ozvdkddzhkzdg_2}{ozvdkddzhkzdx}=\frac{\partial f}{\partial x}-u_1+u_2=0...(2)\\\frac {\partial F}{\partial u_1}=a-x+a_1^2=g_1+a_1^2=0...(3)\\\frac {\partial F}{\partial u_2}=x-b+a_1^2=g_2+b_1^2=0...(4)\\\frac {\partial F}{\partial a_1}=2u_1a_1=0...(5)\\\frac {\partial F}{\partial b_1}=2u_2b_1=0...(6),u_1\geq0,u_2\geq0\end{cases}?????????????????x?F?=?x?f?+u1?dxdg1??+u2?dxdg2??=?x?f??u1?+u2?=0...(2)?u1??F?=a?x+a12?=g1?+a12?=0...(3)?u2??F?=x?b+a12?=g2?+b12?=0...(4)?a1??F?=2u1?a1?=0...(5)?b1??F?=2u2?b1?=0...(6),u1?0,u2?0?

那么現在開始解方程組

首先考慮 式5,

u1=0,a1≠0u1=0,a1≠0u1=0a1?=0時,即g1g_1g1?f(x)f(x)f(x)無約束

u1≠0,a1=0u1≠0,a1=0u1?=0a1=0時,即 g1g_1g1?f(x)f(x)f(x)有約束,且根據式3可知,g1g_1g1? 也等于0。

式6同理。注意,不等式對f(x)有約束效果時,不等式等于零。

此時,方程組簡化成

{?f?x+u1dg1dx+u2dg2dx=0u1g1(x)=0,u2g2(x)=0,μ1≥0,μ1≥0.\begin{cases}\frac{\partial f}{\partial x}+u_1\frac{ozvdkddzhkzdg_1}{ozvdkddzhkzdx}+u_2\frac {ozvdkddzhkzdg_2}{ozvdkddzhkzdx}=0\\u_1g_1(x)=0,u_2g_2(x)=0,\\\mu_1\geq0,\mu_1\geq0.\end{cases}???????x?f?+u1?dxdg1??+u2?dxdg2??=0u1?g1?(x)=0,u2?g2?(x)=0,μ1?0,μ1?0.?

推廣,對于多個不等式約束,有

minf(x)s.t.gj(x)≤0(j=1,2,...,m)min f(x)\\s.t. g_j(x)\leq0(j=1,2,...,m)minf(x)s.t.gj?(x)0(j=1,2,...,m)

我們有

{?f(x?)?x+∑j=1mujdgj(x?)dx?=0ujgj(x?)=0(j=1,2,...,m),μj≥0(j=1,2,...m).\begin{cases}\frac{\partial f(x^*)}{\partial x}+\sum_{j=1}^{m}u_j\frac{ozvdkddzhkzdg_j(x^*)}{ozvdkddzhkzdx^*}=0\\u_jg_j(x^*)=0(j=1,2,...,m),\\\mu_j\geq0(j=1,2,...m).\end{cases}???????x?f(x?)?+j=1m?uj?dx?dgj?(x?)?=0uj?gj?(x?)=0(j=1,2,...,m),μj?0(j=1,2,...m).?

上述不等式稱之為KKT條件,與本文一開始略有不同的是沒有了等式的約束(不用在此糾結,我們關心的是不等式的約束問題),其中uju_juj?稱之為KKT乘子

總結

以上是生活随笔為你收集整理的从拉格朗日乘数法到KKT条件的全部內容,希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯,歡迎將生活随笔推薦給好友。

主站蜘蛛池模板: 色综合久久久久无码专区 | 91精品国产一区二区在线观看 | 久久中文字幕一区二区 | av男人的天堂在线观看 | 性大片潘金莲裸体 | 免费观看黄色的网站 | 国产黄 | 欧美精品日韩少妇 | 中国少妇av | 欧美在线观看免费高清 | 精品啪啪 | 久久机热这里只有精品 | 成人免费激情视频 | 午夜寂寞视频 | 一卡二卡在线视频 | 大陆极品少妇内射aaaaa | 99久久久无码国产精品免费麻豆 | 有码一区二区 | 国产欧美日| www亚洲成人 | 99精品久久毛片a片 成人网一区 | 国内精品99 | 国产精品国产三级国产aⅴ原创 | 少妇太紧太爽又黄又硬又爽 | 久久久久久久美女 | 欧美xxx性 | 色在线播放 | www.黄色网| 青青草原国产在线观看 | 不卡的中文字幕 | 大尺码肥胖女系列av | 国产伦精品一区二区三区免费迷 | 草久av| 国模在线视频 | 欧美一区二区三区免费观看 | wwwxxx日本人 | 四虎永久免费在线观看 | 日本少妇高潮 | 边吃奶边添下面好爽 | 两个小y头稚嫩紧窄h文 | 性欧美长视频 | 成人在线免费视频播放 | 免费色片 | 天啪| 免费看一级 | 国产亚洲精品久久久久四川人 | 天堂视频免费 | 成人免费一区二区三区在线观看 | 国产调教一区 | 日本加勒比一区 | 成人在线播放网站 | 中文在线播放 | 久久久久女人精品毛片九一 | 久久特黄视频 | 少妇又紧又爽视频 | 国产精成人品 | 熊猫电影yy8y全部免费观看 | 国产污污在线观看 | 男男做性免费视频网 | 国产原创视频在线观看 | 亚洲熟妇无码另类久久久 | 女人高潮娇喘1分47秒 | 国产高清毛片 | yy77777丰满少妇影院 | 91手机在线观看 | 蝌蚪网在线视频 | 超碰爱爱 | 福利一二三区 | 国产色诱视频 | 亚洲黄色免费 | 欧美啊v| 国产欧美在线视频 | 青娱乐免费在线视频 | 亚洲精品国产精品国自产观看浪潮 | 在线观看日批 | 99色播 | 精品无码av一区二区三区四区 | av撸撸在线 | 麻豆视频在线播放 | 国产情侣啪啪 | 欧美日韩在线观看视频 | 成人a√| 最近国语视频在线观看免费播放 | 韩国伦理片观看 | 少妇高清精品毛片在线视频 | 欧美真人性野外做爰 | 夜夜嗨av一区二区三区免费区 | 国产成人高清视频 | 女人喂男人奶水做爰视频 | 男人猛吃奶女人爽视频 | 黄视频网站在线 | 在线观看污污视频 | av在线影音 | 国产69精品久久久久久 | 99国产精品白浆在线观看免费 | 扒开伸进免费视频 | jizz视频| 香港台湾日本三级大全 | 正在播放木下凛凛xv99 |