本周源文件是 Lee 的手写作业,题号为 2、4、5、8、9。下面根据扫描中的题面与作答过程还原,并用参考解核对候选点、乘子和二阶型。统一把不等式写成 gi(x)≥0,并取
L(x,w,v)=f(x)−i∑wigi(x)−j∑vjhj(x),wi≥0.
等式乘子的正负只取决于 Lagrange 函数的记号约定,不影响最优性结论。
第 2 题:圆盘与直线交集上的投影
mins.t.(x1−3)2+(x2−2)2−x12−x22+5≥0,x1+2x2−4=0,x1,x2≥0.
第一条是圆盘约束,不是线性约束。由等式消去 x1=4−2x2,圆盘与非负约束给 1≤x2≤2,目标变成
ϕ(x2)=5x22−8x2+5.
无约束极小点 x2=4/5 不在区间内,因此
x∗=(2,1)T.
在该点,圆盘约束和等式都活动。梯度为
∇f=(−2,−2)T,∇g=(−4,−2)T,∇h=(1,2)T.
驻点方程
∇f−w∇g−v∇h=0
给出
w=31,v=−32.
w≥0,KKT 成立。目标严格凸,可行域是圆盘、仿射直线与非负正交象限的凸交集,所以该点还是唯一全局最优解。
第 4 题:逐个检查候选点
根据作答过程还原问题:
mins.t.(x1−49)2+(x2−2)2−x12+x2≥0,−x1−x2+6≥0,x1,x2≥0.
题目给出三个候选点。
候选点 x(1)=(3/2,9/4)T
它满足全部约束,且只有 g1=−x12+x2 活动。此时
∇f=(−23,21)T,∇g1=(−3,1)T.
由 ∇f−w1∇g1=0 得
w1=21>0.
KKT 成立,目标值为
f(x(1))=(−43)2+(41)2=85.
目标严格凸,而可行域是凸函数 x12 的上图集与半空间的交,因此
x∗=(23,49)T,f∗=85
是唯一全局最优解。
候选点 x(2)=(9/4,2)T
代入第一条约束:
−(49)2+2=−1649<0,
所以它不可行,直接排除。
候选点 x(3)=(0,−2)T
扫描和参考解都印出了这个候选点,但它与题面中的 x2≥0 直接矛盾,也不满足 −x12+x2≥0。因此这里存在源材料笔误,不能把它当作可行 KKT 点;手写过程中得到负乘子,也会独立排除它。
第 5 题:等式与不等式混合
根据作答过程还原问题:
mins.t.x12−x2−3x3g(x)=−x1−x2−x3≥0,h(x)=x12+2x2−x3=0.
解 KKT 方程得到
x∗=(−27,−1235,1277)T,w∗=37,v∗=32.
两个约束在 x∗ 处都活动。按本页的记号
L=f−wg−vh,
其 Hessian 为
∇xx2L=2/300000000.
正乘子活动约束与等式约束给临界方向条件
−d1−d2−d3=0,−7d1+2d2−d3=0.
解得
d=(d1,2d1,−3d1)T.
因此对任意非零临界方向,
dT∇xx2Ld=32d12>0.
故 x∗ 是严格局部极小点,目标值为
f(x∗)=−1249.
第 8 题:两圆约束的三个候选点
根据作答过程还原问题:
mins.t.x2g(x)=−x12−(x2−4)2+16≥0,h(x)=(x1−2)2+(x2−3)2−13=0.
题目给出
x(1)=(0,0)T,x(2)=(516,532)T,x(3)=(2,3+13)T.
在 x(1)=(0,0)T
两条约束都活动,KKT 乘子为
w=81,v=0.
临界方向需满足
(0,8)d=0,(−4,−6)d=0,
只能得到 d=0。两个活动边界在该点横截,局部没有非零临界方向,因此二阶充分条件成立,x(1) 是局部极小点。
在 x(2)=(16/5,32/5)T
两条约束也都活动。解驻点方程得到
w=403,v=51.
两条活动约束梯度线性无关,临界方向仍只有 d=0,所以 x(2) 也是局部极小点。
在 x(3)=(2,3+13)T
圆盘约束不活动,故 w=0;等式约束给
v=2613.
等式切方向满足 d2=0。此时
∇xx2L=−131I,
所以对 d=(d1,0)T=0,
dT∇xx2Ld=−13d12<0.
它不是局部极小点。
第 9 题:参数 β 改变二阶性质
mins.t.21[(x1−1)2+x22]−x1+βx22=0.
考察 (0,0)T。驻点方程给等式乘子 v=1,切方向由
−d1=0
得到 d=(0,d2)T。Lagrange Hessian 在切方向上的二阶型是
dT∇xx2Ld=(1−2β)d22.
因此 β<1/2 时二阶充分条件直接给严格局部极小;β>1/2 时存在二阶下降方向,不是局部极小。边界 β=1/2 时二阶型为零,还不能停在必要条件上。
用约束精确消元 x1=βt2, x2=t:
f(t)=21[(βt2−1)2+t2]=21+(21−β)t2+2β2t4.
当 β=1/2 时,非零小 t 仍使四次项严格增加。因此完整结论是
β≤21 时 (0,0) 是严格局部极小点;β>21 时不是局部极小点。
作业自检
- 先验可行性,再列活动约束;不可行候选点无需解 KKT。
- 二阶条件检查 dT∇xx2Ld,不是只看 ∇2f。
- 二阶型等于零时必要条件没有给答案,要像参数题一样回到原函数或更高阶项。