理论教育 约束优化问题局部解的一阶必要条件

约束优化问题局部解的一阶必要条件

时间:2023-06-17 理论教育 版权反馈
【摘要】:对于一般的n维等式约束优化问题:其最优解的一阶必要条件为图2-7 二维等式约束优化问题4.GP型约束问题的最优解条件由上述IP型和EP型约束问题的最优解条件,可推出一般约束优化问题的一阶必要条件。 设二维约束优化问题的数学模型为它的当前迭代点为,试用K-T条件判断其是否为约束最优点。

约束优化问题局部解的一阶必要条件

1.起作用约束与不起作用约束

如前所述,一般的约束优化问题包含有不等式约束guX)≥0(u=1,2,…,p)和等式约束hvX)=0(v=1,2,…,q),在可行设计点Xk处,若有guX)=0,则该约束guX)称为可行点Xk的起作用约束;而若guXk)>0,则该约束guX)称为可行点Xk的不起作用约束。对于等式约束hvX)=0,显然在任一可行点处的等式约束都是起作用约束。

在可行点Xk处,起作用约束在Xk的邻域内起限制可行域的作用,而不起作用约束在Xk处的邻域内就不产生影响。因此应把注意力集中在起作用约束上。

2.IP型约束问题的最优解条件

图2-6所示为具有三个不等式约束的二维约束优化问题,分为以下三种情况。

1)最优点X*在可行域内(见图2-6a)。在此情况下,X*点的全部约束函数值guX*)>0(u=1,2,3),即这组约束条件对其最优点X*都不起作用。换言之,如果去掉全部约束,其最优点也仍是同一个X*点。因此,这种约束优化问题与无约束优化问题是等价的。

978-7-111-29617-1-Chapter02-75.jpg

图2-6 二维不等式约束优化问题

2)约束最优点X*位于g1X)的边界曲线与目标函数等值线的切点处(见图2-6b)。此时g1X*)=0,g2X*)>0,g3X*)>0。所以,起作用的约束仅为g1X)。既然X*fX)等值线与g1X)边界线的切点,则对应的梯度矢量ΔfX*)与Δg1X*)必共线,且方向一致。若取非负乘子λ1*≥0,则在X*处存在如下关系:

978-7-111-29617-1-Chapter02-76.jpg

3)当前迭代点Xk在两约束的交点上(见图2-6c)。这时,ΔfXk)夹于Δg1Xk)与Δg2Xk)之间。显然,在Xk点处,起作用约束为g1X)和g2X)。在这种情况下,由于Xk点邻近的可行域内部不存在目标函数值比fXk)更小的设计点,因而有:Xk=X*。由图2-6c可知,此时Xk点的目标函数的梯度ΔfXk)可表达为约束函数梯度Δg1Xk)和Δg2Xk)的线性组合。若用X*代替Xk,则有

978-7-111-29617-1-Chapter02-77.jpg

成立,且有λ1*≥0,λ2*≥0。

总结以上情况,对图2-6所示的二维约束优化问题,其最优解的一阶必要条件是

978-7-111-29617-1-Chapter02-78.jpg

根据上面的分析,可归纳出IP型约束问题最优解的一阶必要条件如下:

设最优点X*位于J个约束边界的汇交处,则J个约束条件组成一个起作用的约束集,对于X*点,下式必成立:

978-7-111-29617-1-Chapter02-79.jpg

然而,在实际求解过程中,并不能预先知道最优点X*位于哪一个或哪几个约束边界的汇交处。为此,把p个约束全部考虑进去,并取相应于不起作用约束的乘子λu*为零,这样,上式可改写成:

978-7-111-29617-1-Chapter02-80.jpg

此即为IP型约束优化问题最优解一阶必要条件的一般表达式,它实际上等价于式(2-30)。这是因为,在X*处,对于起作用约束,必有guX*)=0(u=1,2,…,J);对于不起作用约束,尽管guX*)>0,但λu=0。

3.EP型约束问题的最优解条件

图2-7所示为具有一个等式约束条件的二维优化问题,其数学模型

978-7-111-29617-1-Chapter02-81.jpg

该问题中,等式约束曲线hX)=0是它的可行域,而且目标函数等值线fX)=ChX)=0的切点X*是该约束问题的最优点。在X*点处,目标函数的梯度ΔfX)与约束函数的梯度ΔhX*)共线。因此,必定存在一个乘子μ*,使得ΔfX*)-μ*ΔhX*)=0成立。对于一般的n维等式约束优化问题:

978-7-111-29617-1-Chapter02-82.jpg

其最优解的一阶必要条件为

978-7-111-29617-1-Chapter02-83.jpg

978-7-111-29617-1-Chapter02-84.jpg

图2-7 二维等式约束优化问题

4.GP型约束问题的最优解条件

由上述IP型和EP型约束问题的最优解条件,可推出一般约束优化问题的一阶必要条件。设n维一般约束优化问题的数学模型为

978-7-111-29617-1-Chapter02-85.jpg

X*为其最优解的一阶必要条件应为

978-7-111-29617-1-Chapter02-86.jpg

若令 978-7-111-29617-1-Chapter02-87.jpg

则称LXλμ)为关于问题(2-32)的广义拉格朗日(Lagrange)函数,式中,978-7-111-29617-1-Chapter02-88.jpg978-7-111-29617-1-Chapter02-89.jpg拉格朗日乘子。(www.daowen.com)

于是,式(2-33)中第一式可写成:

978-7-111-29617-1-Chapter02-90.jpg

5.库恩-塔克(Kuhn-Tucker)条件

在实际优化计算中,常常需要判断某可行迭代点Xk是否可作为约束最优点X*输出而结束迭代,或者检验所输出的可行结果是否已满足约束最优解的必要条件。这种判断或检查通常是借助于库恩-塔克条件(简称K-T条件)进行的。

K-T条件指出,点Xk要成为约束优化问题最优解的必要条件是:目标函数fX)的负梯度-ΔfXk)可以表示为起作用约束函数梯度ΔguXk)和ΔhvXk)的线性组合。即

978-7-111-29617-1-Chapter02-91.jpg

式中其他符号的意义同前。当不等式约束guXk)≥0时,上式可写成:

978-7-111-29617-1-Chapter02-92.jpg

若引入广义拉格朗日函数及如下的向量:

978-7-111-29617-1-Chapter02-93.jpg

则式(2-37)可写成更紧凑形式:

978-7-111-29617-1-Chapter02-94.jpg

凡满足式(2-36)、(2-37)或(2-38)的迭代点978-7-111-29617-1-Chapter02-95.jpg即可以作为最优点输出。这时,978-7-111-29617-1-Chapter02-96.jpg称为K-T点。

【例2-5】 设二维约束优化问题的数学模型为

978-7-111-29617-1-Chapter02-97.jpg

它的当前迭代点为978-7-111-29617-1-Chapter02-98.jpg,试用K-T条件判断其是否为约束最优点。

解(1)计算Xk)点的诸约束函数值

978-7-111-29617-1-Chapter02-99.jpg

可知Xk点是可行点。

(2)Xk点起作用的约束

g1X)=1-x21-x2=0

g2X)=x2

(3)求Xk点的诸梯度

978-7-111-29617-1-Chapter02-100.jpg

(4)求拉格朗日乘子 把以上计算结果代入式(2-37),可得

978-7-111-29617-1-Chapter02-101.jpg

写成线性方程组

-2+2λ1=0

λ1-λ2=0

解得:λ1=λ2=1>0,满足K-T条件。所以,978-7-111-29617-1-Chapter02-102.jpg为约束局部最优点。而且,容易验证,fX)是凸函数,可行域为凸集,故978-7-111-29617-1-Chapter02-103.jpg也是该问题的全局最优解,如图2-8所示。

978-7-111-29617-1-Chapter02-104.jpg

图2-8 例2-5约束优化问题图解

6.K-T条件的几何意义

K-T条件的几何意义如图2-9所示,起作用约束函数的梯度向量,在设计空间内构成一个锥体。若满足K-T条件,则目标函数的负梯度向量-ΔfX)应包含在此锥体之内。显然,图2-9a表示X并非目标函数的最优点,因为它不满足K-T条件;而图2-9b所示的X才是目标函数的一个局部最优点。

978-7-111-29617-1-Chapter02-105.jpg

图2-9 K-T条件的几何解析

免责声明:以上内容源自网络,版权归原作者所有,如有侵犯您的原创版权请告知,我们将尽快删除相关内容。

我要反馈