最优性条件

等式约束问题

image-20231122194143520

image-20231122200155941

image-20231122200142522


不等式约束问题

image-20231122201051270

定义II:起作用约束标号集 ,将容许点带入s等于零的s的序号集合,比如{1,2},容许点带入s等于零的s被称作起作用约束


不等式的FJ条件

image-20231123204144358


不等式的KT条件

image-20231123205800495

KT条件在FJ条件的前提下要求 起作用约束条件的梯度s\nabla s 非线性相关


例题:

image-20231123210431617

image-20231123210602798

image-20231123210925116

image-20231123210946523


一般约束问题

image-20231123211028836

同时具有等式约束和不等式约束


与不等式约束相似,一般约束问题的FJ和KT条件在不等式约束的基础上添加了等式部分,即h和λ\lambda,不同的是等式部分一定为起作用约束

FJ条件

image-20231123211321498

KT条件

image-20231123211654194


例题:

image-20231123211807438

Zoutendijk 容许方向法

image-20231123212437354

容许方向法是求解约束最优化问题的一类基本方法。这类方法一般先从线性约束问题开始讨论,然后再推广到非线性约束问题。它的基本迭代格式是

image-20231123212734457

下降容许方向的确定

非零向量pp 为点xx容许方向向量的充要条件是pp满足 Ap>=0,Cp=0A^`p>=0,Cp=0

CC是等式约束的系数矩阵, AA^`是不等式约束的起作用约束部分系数矩阵

image-20231123213239415

若p还满足f(x)Tp<0\nabla f(x)^Tp<0,那么p是下降容许方向向量

确认p是否为下降容许向量,即求下式最优值,其中e为分量全为1的n维向量,即[1,1,1,1,1...1]T[1,1,1,1,1...1]^T

image-20231123214509812

该式最优值非正, 若最优值<0,则p是下降容许方向向量, 若最优值=0,则容许点x是KT点

直线搜索

AA^{``}是不等式约束的不起作用约束部分系数矩阵

image-20231123220145278

终止条件

image-20231123220617860

算法

image-20231123221415613

image-20231123221714365

例题见书p245

外部罚函数法

image-20231125145955593

其中uu是阶跃函数image-20231125150206437

image-20231125150445856

image-20231125150554816

若无约束问题(4.79)的极小点x 是约束问题(4.67)的容许点(即α(x)=0\alpha (x)=0),则x也是约束问题(4.67)的极小点,表明约束不起作用,这时只解一次无约束问题(4.80)就能求到约束问题(4.67)的极小点

若不是容许点(即α(x)0\alpha(x) ≠ 0),则x不是约束问题(4.67)的极小点,表明约束对极小点的求解有影响,这时就要解一系列无约束问题(4.80),得到的极小点序列{xkx_k }向容许集靠近,向(4.67)的极小点逼近


算法

image-20231125151834861

image-20231125152016338


例子

image-20231125163443026

image-20231125163459706

image-20231125163519936

乘子法

​ 乘子法解决了外部罚函数法因罚因子增大导致的数值计算不稳定问题,是外部罚函数法的一种改进方法.
​ 乘子法是在约束问题的 Lagrange 函数上加入惩罚,罚因子因此不必趋于无穷大,数值计算稳定性有保证.


等式约束问题

image-20231125161416397

image-20231125161432549

image-20231125161441781

image-20231125161522865

image-20231125162228477

如果无约束问题(4.121)的最优解x 是(4.109)的容许点,那么x 是(4.109)的最优解,且参数$\lambda =\lambda^, μ=μ^ $

image-20231125162553884

终止准则 h(x)=0h(x)=0


算法:


例子

image-20231125163849047

image-20231125163907469

不等式约束问题

image-20231125164124451

image-20231125164139093


例子

一般约束问题

image-20231125164917728

image-20231125164855857

image-20231125164937553


算法

image-20231125165030218