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

歡迎訪問 生活随笔!

生活随笔

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

编程问答

UA SIE545 优化理论基础3 Fritz-John与Kuhn-Tucker理论总结 带等式约束与不等式约束的极值问题

發布時間:2025/4/14 编程问答 31 豆豆
生活随笔 收集整理的這篇文章主要介紹了 UA SIE545 优化理论基础3 Fritz-John与Kuhn-Tucker理论总结 带等式约束与不等式约束的极值问题 小編覺得挺不錯的,現在分享給大家,幫大家做個參考.

UA SIE545 優化理論基礎3 Fritz-John與Kuhn-Tucker理論總結 帶等式約束與不等式約束的極值問題

對于函數f:X→Yf:X \to Yf:XY,我們希望XXX是一個凸的度量空間(或者拓撲空間),YYY是一個定義了良序的內積空間。N?(x)N_{\epsilon}(x)N??(x)表示中心為xxx,半徑為?\epsilon?的領域,下面先列出幾個會用到的概念。

Global minimum If xˉ∈X\bar x \in XxˉX, ?x∈X\forall x \in X?xX, f(xˉ)≤f(x)f(\bar x) \le f(x)f(xˉ)f(x), then xˉ\bar xxˉ is a global minimum. If equality in f(xˉ)≤f(x)f(\bar x) \le f(x)f(xˉ)f(x) will never hold, then xˉ\bar xxˉ is strict global minimum.

Local minimum If xˉ∈X\bar x \in XxˉX, ??>0\exists \epsilon>0??>0, ?x∈N?(xˉ)\forall x \in N_{\epsilon}(\bar x)?xN??(xˉ), f(xˉ)≤f(x)f(\bar x) \le f(x)f(xˉ)f(x), then xˉ\bar xxˉ is a local minimum. If equality in f(xˉ)≤f(x)f(\bar x) \le f(x)f(xˉ)f(x) will never hold, then xˉ\bar xxˉ is strict local minimum.

Pseudoconvex ?x1,x2∈X\forall x_1,x_2 \in X?x1?,x2?X such that ?f(x1)T(x2?x1)≥0\nabla f(x_1)^T(x_2-x_1) \ge 0?f(x1?)T(x2??x1?)0, we have f(x2)≥f(x1)f(x_2) \ge f(x_1)f(x2?)f(x1?).

Quasiconvex ?x1,x2∈X,λ∈[0,1]\forall x_1,x_2 \in X,\lambda \in [0,1]?x1?,x2?X,λ[0,1] such that f(λx1+(1?λ)x2)≤max?(f(x),f(y))f(\lambda x_1+ (1-\lambda)x_2) \le \max(f(x),f(y))f(λx1?+(1?λ)x2?)max(f(x),f(y))


min?x∈Sf(x)\min_{x \in S} f(x)minxS?f(x),其中SSSXXX的子集。

Cone of feasible directions at xˉ∈clS\bar x \in clSxˉclS,
D={d≠0:?δ>0,?λ∈(0,δ),xˉ+λd∈S}D=\{d \ne 0:\exists \delta>0, \forall \lambda \in (0,\delta),\bar x+\lambda d \in S\}D={d?=0:?δ>0,?λ(0,δ),xˉ+λdS}

Each d∈Dd \in DdD is called feasible direction.

Cone of improving direction at xˉ∈clS\bar x \in clSxˉclS,
F={d:?δ>0,?λ∈(0,δ),f(xˉ+λd)<f(xˉ)}F=\{d:\exists \delta>0,\forall \lambda \in (0,\delta),f(\bar x+\lambda d)<f(\bar x)\}F={d:?δ>0,?λ(0,δ),f(xˉ+λd)<f(xˉ)}Each d∈Fd \in FdF is called improving direction.


無約束極值問題
min?x∈Xf(x)\min_{x \in X}f(x)xXmin?f(x)

Let xˉ∈X\bar x \in XxˉX be a candidate point. PSD = positive semi-definite, PD = positive definite.

First order necessary condition:
?f(xˉ)=0\nabla f(\bar x) = 0?f(xˉ)=0

Second order necessary condition:
H(xˉ)PSDH(\bar x)\ PSDH(xˉ)?PSD

Second order sufficient condition:
H(xˉ)PDH(\bar x)\ PDH(xˉ)?PD

For convex optimization, sufficient and necessary condition:
?f(xˉ)=0\nabla f(\bar x) = 0?f(xˉ)=0


含等式約束與不等式約束的極值問題
min?x∈Xf(x)s.t.g(x)≤0h(x)=0\min_{x \in X} f(x) \\ s.t. g(x) \le 0 \\ h(x)=0xXmin?f(x)s.t.g(x)0h(x)=0

Let LLL be Lagrange function:
L(x)=f(x)+uTg(x)+vTh(x)L(x) = f(x)+u^Tg(x)+v^Th(x)L(x)=f(x)+uTg(x)+vTh(x)

Fritz-John (FJ First order necessary condition):
u0T?f(xˉ)+uT?g(xˉ)+vT?h(xˉ)=0uTg(xˉ)=0(u0,u)≥0,(u0,u,v)≠0u_0^T\nabla f(\bar x)+u^T\nabla g(\bar x)+v^T\nabla h(\bar x)=0 \\ u^T g(\bar x) = 0 \\ (u_0,u) \ge 0, (u_0,u,v) \ne 0u0T??f(xˉ)+uT?g(xˉ)+vT?h(xˉ)=0uTg(xˉ)=0(u0?,u)0,(u0?,u,v)?=0

Karush-Kuhn-Tucker (KKT first order necessary condition):
g(xˉ)≤0?f(xˉ)+uT?g(xˉ)+vT?h(xˉ)=0,u≥0uTg(xˉ)=0g(\bar x) \le 0\\ \nabla f(\bar x)+u^T\nabla g(\bar x)+v^T\nabla h(\bar x)=0,u \ge 0 \\ u^T g(\bar x) = 0 g(xˉ)0?f(xˉ)+uT?g(xˉ)+vT?h(xˉ)=0,u0uTg(xˉ)=0

and ?g(xˉ)\nabla g(\bar x)?g(xˉ) and ?h(xˉ)\nabla h(\bar x)?h(xˉ) are linearly independent. 第一個條件叫primal feasibility(PF)、第二個條件dual feasibility(DF),第三個條件約束的是complementary slacks (CS),最后一個條件叫constraint qualification,保證KKT與FJ等價。

KKT second-order necessary condition:
HL(xˉ)PSDHL(\bar x)\ PSDHL(xˉ)?PSD

KKT second-order sufficient condition:
HL(xˉ)PDHL(\bar x)\ PDHL(xˉ)?PD


考慮下面的優化問題
min?f(x1,x2)=(x1?1)2+x22s.t.g(x1,x2)=2kx1?x22≤0,k>0\min f(x_1,x_2)=(x_1-1)^2+x_2^2 \\ s.t.g(x_1,x_2)=2kx_1-x_2^2 \le 0, k > 0minf(x1?,x2?)=(x1??1)2+x22?s.t.g(x1?,x2?)=2kx1??x22?0,k>0

因為
Hf(x1,x2)=[2002],Hg(x1,x2)=[000?2]Hf(x_1,x_2) = \left[ \begin{matrix} 2 & 0 \\ 0 & 2 \end{matrix}\right], \ Hg(x_1,x_2) = \left[ \begin{matrix} 0 & 0 \\ 0 & -2 \end{matrix}\right]Hf(x1?,x2?)=[20?02?],?Hg(x1?,x2?)=[00?0?2?]

顯然fff是凸函數但ggg不是,因此這個優化不是凸優化,那么如何計算這個優化的最優解呢?


注意到目標函數是圓,約束的邊界是拋物線,所以可以用作圖法求解,非常簡單讀者可以自行嘗試。我們來討論一下代數解法,考慮KKT條件,但需要注意這個問題非凸,所以KKT條件并不是充分條件。因為
?g(x1,x2)=2[k?x2],k>0\nabla g(x_1,x_2) = 2 \left[ \begin{matrix} k \\ -x_2 \end{matrix} \right],k>0?g(x1?,x2?)=2[k?x2??],k>0

也就是說?g(x1,x2)\nabla g(x_1,x_2)?g(x1?,x2?)非零,所以?g(xˉ)\nabla g(\bar x)?g(xˉ) and ?h(xˉ)\nabla h(\bar x)?h(xˉ) are linearly independent這個條件成立,因此KKT是必要條件。

定義Lagrange函數:
L(x1,x2,u)=f(x1,x2)+ug(x1,x2)=(x1?1)2+x22+u(2kx1?x22)L(x_1,x_2,u)=f(x_1,x_2)+ug(x_1,x_2) \\ =(x_1-1)^2+x_2^2 +u(2kx_1-x_2^2 ) L(x1?,x2?,u)=f(x1?,x2?)+ug(x1?,x2?)=(x1??1)2+x22?+u(2kx1??x22?)

寫出KKT條件:
{2(x1?1)+2ku=02x2?2ux2=0ug(x1,x2)=u(2kx1?x22)=0g(x1,x2)=2kx1?x22≤0u≥0\begin{cases} 2(x_1-1)+2ku = 0 \\ 2x_2 -2ux_2 = 0 \\ug(x_1,x_2) = u(2kx_1-x_2^2) = 0 \\ g(x_1,x_2)=2kx_1-x_2^2 \le 0 \\ u \ge 0 \end{cases}????????????????2(x1??1)+2ku=02x2??2ux2?=0ug(x1?,x2?)=u(2kx1??x22?)=0g(x1?,x2?)=2kx1??x22?0u0?

接下來做分類討論:

  • u=0u=0u=0, 則x2=0x_2=0x2?=0, x1=1x_1=1x1?=1, not feasible (x1=1,x2=0?g(x1,x2)>0x_1=1,x_2=0 \Rightarrow g(x_1,x_2) > 0x1?=1,x2?=0?g(x1?,x2?)>0)
  • u>0u>0u>0, 第二個方程說明2x2(u?1)=02x_2(u-1)=02x2?(u?1)=0,如果x2=0x_2 =0x2?=0,則x1=0,u=1/kx_1=0,u=1/kx1?=0,u=1/k;如果x2≠0x_2 \ne 0x2??=0, 則u=1u=1u=1,那么x1=1?kx_1=1-kx1?=1?k, x2=±2k(1?k)x_2=\pm \sqrt{2k(1-k)}x2?=±2k(1?k)?
  • 現在我們有了三組可能的解,但還需要根據kkk的值做分類討論:

  • k>1k>1k>1,只有(0,0)(0,0)(0,0)是KKT point,因此(0,0)(0,0)(0,0)是最優解;
  • k≤1k \le 1k1, (0,0),(1?k,2k(1?k)),(1?k,?2k(1?k))(0,0),(1-k,\sqrt{2k(1-k)}),(1-k,-\sqrt{2k(1-k)})(0,0),(1?k,2k(1?k)?),(1?k,?2k(1?k)?)都是KKT point,我們需要通過比較它們對應的目標函數的取值來確定誰是最優解
  • {f(0,0)=1f(1?k,±2k(1?k))=2k?k2\begin{cases} f(0,0) = 1 \\ f(1-k,\pm \sqrt{2k(1-k)})=2k-k^2 \end{cases}{f(0,0)=1f(1?k,±2k(1?k)?)=2k?k2?

    因為
    2k?k2?1=?(k?1)2≤02k-k^2-1 = -(k-1)^2 \le 02k?k2?1=?(k?1)20

    所以(1?k,2k(1?k)),(1?k,?2k(1?k))(1-k,\sqrt{2k(1-k)}),(1-k,-\sqrt{2k(1-k)})(1?k,2k(1?k)?),(1?k,?2k(1?k)?)是最優點。

    總結

    以上是生活随笔為你收集整理的UA SIE545 优化理论基础3 Fritz-John与Kuhn-Tucker理论总结 带等式约束与不等式约束的极值问题的全部內容,希望文章能夠幫你解決所遇到的問題。

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

    主站蜘蛛池模板: 上床视频在线观看 | 国产欧美视频一区 | 国产精品久久av无码一区二区 | 老司机精品视频在线 | 天天干,夜夜爽 | 亚洲一区二区色图 | 亚洲欧美色图视频 | 国产精品久久久久久久久久久久久久久 | 无码人妻精品一区二区蜜桃色欲 | 韩国黄色一级片 | 日本美女黄网站 | 欧美破处女 | 日韩成人av毛片 | 欧美美女性视频 | 日本人体一区 | 国产精品视频观看 | 亚洲精品中文无码AV在线播放 | 熟睡侵犯の奶水授乳在线 | 久久99国产精品成人 | www三级免费 | 日本精品在线播放 | 99久久久无码国产 | av动漫网 | 欧美激情国产在线 | 尹人av| 日本精品一区在线观看 | 欧美丰满艳妇bbwbbw | 绿帽h啪肉np辣文 | 精品一区二区三区四区视频 | 亚洲综合精品在线 | 麻豆av在线免费观看 | 免费av观看 | 男女深夜福利 | 91桃色视频在线观看 | 国产91在线免费 | 免费处女在线破视频 | 在线一区二区三区四区 | 中文字幕精 | 老头吃奶性行交 | 久久久无码精品亚洲国产 | 丰满少妇一级片 | 亚洲最新中文字幕 | 九色视频偷拍少妇的秘密 | 久久国产精品久久 | 亚洲国产免费av | www.97视频| 欧美yyy | 亚洲中文无码久久 | 免费a网 | 国产精品久久久久久久久久久久久久久久久 | 黄色片日韩| 91久久综合精品国产丝袜蜜芽 | 国产精品午夜在线观看 | av777777| 婷婷色亚洲 | 国模大尺度视频 | 国产区网址 | 天天爽天天搞 | www污网站| 大尺度做爰呻吟舌吻情头 | 麻豆激情网 | 香蕉视频色 | 国产欧美日韩在线观看 | 日日夜夜精品视频免费 | 免费看日批视频 | 91看片黄 | 高清视频一区二区三区 | 亚洲黄色小说网 | 日韩国产综合 | 五月婷婷免费视频 | 久久久久久亚洲av毛片大全 | 亚洲国产精品久久精品怡红院 | 久久亚洲精选 | 国产主播精品 | 男女午夜免费视频 | 久久牛牛 | 午夜裸体性播放 | 欧美黑人啪啪 | 久久在线视频 | 狠狠操中文字幕 | 插插看看 | 91丝袜呻吟高潮美腿白嫩 | xxxx少妇| 亚洲乱码国产乱码精品精剪 | 国产av无码专区亚洲av毛片搜 | 国产精品99精品 | 琪琪女色窝窝777777 | 爱豆国产剧免费观看大全剧集 | 色图自拍偷拍 | 天堂国产一区二区三区 | 色婷婷www | 99mav| 黑丝一区| 久久精品无码一区二区三区免费 | 婷婷精品一区二区三区 | 国产美女无遮挡永久免费 | 超碰97色| 欧美理论在线 | 黄色免费91 |