中科大-凸优化 笔记(lec25)-等价变换
全部筆記的匯總貼(視頻也有傳送門):中科大-凸優化
?fT(x?)(y?x)≥0\nabla f^T(x^*)(y-x)\ge0?fT(x?)(y?x)≥0線性規劃的解在邊界上
一、等價變換
例:食譜問題
選擇食譜,使得mmm種營養元素分別不小于b1,?,bmb_1,\cdots,b_mb1?,?,bm?
有nnn種食物,單位食物營養為a1j,a2j,?,amj,?j=1,?,na_{1j},a_{2j},\cdots,a_{mj},\forall j=1,\cdots,na1j?,a2j?,?,amj?,?j=1,?,n
目標:單位食物價格CjC_jCj?,找出總價最小的食譜
設食物量分別為x1,?,xnx_1,\cdots,x_nx1?,?,xn?
min∑j=1nCjxjs.t.∑j=1mCijxj≥bj,i=1,?,m,xj≥0,j=1,?,nmin\;\;\sum_{j=1}^nC_jx_j\\s.t.\;\;\sum_{j=1}^m C_{ij}x_j\ge b_j,i=1,\cdots,m,x_j\ge0,j=1,\cdots,nminj=1∑n?Cj?xj?s.t.j=1∑m?Cij?xj?≥bj?,i=1,?,m,xj?≥0,j=1,?,n
min[C1,?,Cn][x1?xn]s.t.[a11?a1n??am1?amn][x1?xn]=[b1?bm],[x1?xn]≥0min\;\;[C_1,\cdots,C_n]\left[ \begin{array}{c}x_1\\\vdots\\x_n\end{array} \right] \\s.t.\left[ \begin{array}{c} a_{11}\cdots a_{1n}\\\vdots\;\;\;\;\;\;\;\;\;\vdots \\a_{m1}\cdots a_{mn} \end{array} \right]\left[ \begin{array}{c} x_1\\\vdots \\x_n \end{array} \right]=\left[ \begin{array}{c} b_1\\\vdots \\b_m \end{array} \right],\left[ \begin{array}{c} x_1\\\vdots \\x_n \end{array} \right]\ge0min[C1?,?,Cn?]????x1??xn??????s.t.????a11??a1n???am1??amn??????????x1??xn??????=????b1??bm??????,????x1??xn??????≥0
例:線性分數規劃(linear fractional programing)
(P0P_0P0?)min?f0(x)\min f_0(x)minf0?(x)線性分數函數(凸×擬凸√)f0(x)=CTx+deTx+fdomf0={x∣eTx+f>0}f_0(x)=\frac{C^Tx+d}{e^Tx+f}\;dom\;f_0=\{x|e^Tx+f>0\}f0?(x)=eTx+fCTx+d?domf0?={x∣eTx+f>0}
s.t.Gx≤h,Ax=bs.t.\;\;Gx\le h,Ax=bs.t.Gx≤h,Ax=b
不是標準的凸問題,需要轉化
?(P1)\Rightarrow(P_1)?(P1?)
min?CTy+dz\min C^Ty+dzminCTy+dz
s.t.Gy?hz≤0,eTx+fz=1,Ay?bz=0,z≥0(z>0等價)s.t.\;Gy-hz\le0,e^Tx+fz=1,Ay-bz=0,z\ge0(z>0等價)s.t.Gy?hz≤0,eTx+fz=1,Ay?bz=0,z≥0(z>0等價)
證明
(1)若xxx在P0P_0P0?可行,y=xeTx+f,z==1eTx+fy=\frac x{e^Tx+f},z==\frac 1{e^Tx+f}y=eTx+fx?,z==eTx+f1?
{Gx≤hAx=beTx+f>0\left\{ \begin{array}{l} Gx\le h \\ Ax=b \\ e^Tx+f>0 \end{array} \right. ????Gx≤hAx=beTx+f>0?
Gy?hz=Gx?heT+f≤0,eTx+fz=eTx+feTx+f=1,Ay?bz=Ax?beTx+f=0z≤0Gy-hz=\frac{Gx-h}{e^T+f}\le0,e^Tx+fz=\frac{e^Tx+f}{e^Tx+f}=1,Ay-bz=\frac{Ax-b}{e^Tx+f}=0\;\;z\le0Gy?hz=eT+fGx?h?≤0,eTx+fz=eTx+feTx+f?=1,Ay?bz=eTx+fAx?b?=0z≤0
cTy+dz=cTx+deTx+f=f0(x)c^Ty+dz=\frac{c^Tx+d}{e^Tx+f}=f_0(x)cTy+dz=eTx+fcTx+d?=f0?(x)
(2)若y,zy,zy,z在P1P_1P1?種可行,若z>0z>0z>0,則x=yzx=\frac yzx=zy?,則xxx在P0P_0P0?中可行且兩問題目標函數值相同。
(3)(在下一小結講,筆記也做在這里算了)
若y,zy,zy,z在P1P_1P1?中可行,若z=0z=0z=0,設x0x_0x0?為P0P_0P0?的可行解
則x=x0+tyx=x_0+tyx=x0?+ty對P0P_0P0?可行,?t≥0\forall t\ge0?t≥0
Gy≤0,Ay=0,eTy=1Gx=Gx0+Gy≤hAx=Ax0+tAh=beTx+f=eTx0+f+teTy>0f0(x)=f0(x0+ty)=CTx0+CTty+deTx0+eTy+f=(t→∞)CTyGy\le0,Ay=0,e^Ty=1\\Gx=Gx_0+Gy\le h\\Ax=Ax_0+tAh=b\\e^Tx+f=e^Tx_0+f+te^Ty>0\\f_0(x)=f_0(x_0+ty)=\frac{C^Tx_0+C^Tty+d}{e^Tx_0+e^Ty+f}=(t\rightarrow \infty)C^TyGy≤0,Ay=0,eTy=1Gx=Gx0?+Gy≤hAx=Ax0?+tAh=beTx+f=eTx0?+f+teTy>0f0?(x)=f0?(x0?+ty)=eTx0?+eTy+fCTx0?+CTty+d?=(t→∞)CTy
下一章傳送門:中科大-凸優化 筆記(lec26)-二次規劃
總結
以上是生活随笔為你收集整理的中科大-凸优化 笔记(lec25)-等价变换的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: ker矩阵是什么意思_“拨开迷雾”,如何
- 下一篇: 敬告青年工学家