约束非线性规划讲解

合集下载

第六讲线性规划与非线性规划

第六讲线性规划与非线性规划
f=f(x); •
(2)若有非线性约束条件:c1 x 0 或c2 x 0, 则建立M
文件c.m定义函数c1 x,c2 x, 一般形式为
function [c1,c2]=c(x)
c1=…
c2=… (3)建立主程序。求解非线性规划的函数是fmincon,
调用格式为 x=fmincon(‘fun’,x0,A1,b1);
故它属于一个整数线性规划问题,这里当成一个线 性规划求解,求得最优解刚好是整数x1=9,x2=0, 故它就是该整数规划的最优解.若用线性规划解法求 得的最优解不是整数,将其取整后不一定是相应整 数规划的最优解,这样的整数规划应用专门的方法 求解.
二、非线性规划
1、二次规划

标准形式:min
z
1
xT
x1 4x2 5

x1, x2 0

改写成标准形式:min z
x1 2x2
1 2
x12
1 2
x22
s.t.
2x1 3x2 x1 4x2
6 5
0 0
0 0
x1 x2
❖ 建立M文件fun1.m
❖ 建立主程序(见MATLAB程序(feixianxingguihua1))
工费用如下表.问怎样分配车床的加工任务,才能既满足加
工工件的要求,又使加工费用最低?
车床 类型


单位工件所需加工台时数 工件 1 工件 2 工件 3
0.4
1.1
1.0
0.5
1.2
1.3
单位工件的加工费用 工件 1 工件 2 工件 3
13
9
10
11
12
8
可用台 时数
800

第四章 约束非线性优化的理论与方法

第四章 约束非线性优化的理论与方法


h( x) x1 x2 6 0.
f
解:
(
x
)

(2
x1,1)T ,g1(
2x1 1
x) (1,0)T
21 x2 1
, g2 0,
(
x
)

(
2
x1
,2
x2
)T
,
h1
(
x
)

(1,1)
;
K T条件为

1 22 x2 1 0,
x* (1,0)T 是 最 优 点 , 起 作 用 约 束I ( x*) {1,3},f ( x*) (1,0)T
g1( x*) (0,1)T , g3 ( x*) (0,1)T , 显 然 找 不 到1 , 3 0, 使 f ( x*) 1g1( x*) 5g3 ( x*)成 立 。
F(x0)= f (x0) + g (x0)=0.
若 F ( x0 , h) hT 2F ( x0 )h 0, h 0
gi ( x0; h) 0, i 1,2,, m
则x0是局部极小点。 二,具不等式约束的问题
1,下降方向和可行方向
考虑一般非线性约束优化问题:
满足:
f
(
x

)

m

i 1
1 ,m T ;
p
i gi ( x ) jhj ( x
j1
1 )
,,
0
ቤተ መጻሕፍቲ ባይዱ

p
T

i gi ( x ) 0, i 0, i 1,, m
称上式为K-T条件,满足上式的点称K-T点。

高级运筹学-第9章:非线性规划

高级运筹学-第9章:非线性规划
第 九 章:非线性规划模型
9.1 基本概念和基本原理
一、什么是非线性规划:
目标函数和约束条件中有非线性函数的规划问题。
例9-1 某企业生产一种产品y需要生产资料x1和x2,用经济计量学方法根据统计 资料可写出生产函数为:
1/ 3 2/ 3 y 2 x1 x2
、 但是投入的资源有限,能源总共1O个单位,而每单位生产资料 x1要消耗1单位 能源,每单位生产资料x2要消耗2单位能源。问:应如何安排生产资料使产出最大?
y=f(x)
X2xBiblioteka 2、性质:fi(X) 为凸集R上的凸函数,则对ki≥0,i=1,2,…,m,有
k1f1(X)+k2f2(X)+…+kmfm(X)仍为凸函数。
3、凸函数的判定:f(X) 定义在凸集R上,若f(X)有连续的二阶导数, 则f(X)为凸函数 H为半正定。
f(X)为严格凸函数 H为正定。
则停止,得近似极值点 . 否则,继续缩短区间, 直至满足给定的精度为 止。
② f(x2)>f(x1),取[a1=x2,b1=b0] x’1=a1+(b1-a1) x’2=x1 ⑵ 若求f(x)的极大值点,则 ① f(x2)≥f(x1),取[a1=a0,b1=x1] x’1=x2 x’2=b1-(b1-a1) ② f(x2)<f(x1),取[a1=x2,b1=b0] x’1=a1+(b1-a1) 10 x’2=x1
例 求解 f(x)=-18x2+72x+28 的极大值点,≤0.1,起始搜索区间为[0,3]
解:①用间接法:令 f’(x)=-36x+72=0,得驻点 x=2 又因为f’’(x)=-36<0,故 x=2 为f(x)的极大值点,fmax=100 ②用直接法中的黄金分割法:令 n-1=,得n=1+(lg)/(lg)≈5.78442 约6步即可,计算结果见下表:

非线性规划作业

非线性规划作业

非线性规划作业非线性规划是运筹学中的一个重要分支,用于解决具有非线性目标函数和约束条件的优化问题。

本次作业将介绍非线性规划的基本概念、求解方法和应用,并提供一个实际问题供你进行求解和分析。

一、基本概念1. 非线性规划:非线性规划是指目标函数或者约束条件中至少有一个是非线性的优化问题。

它与线性规划相比,更具有灵便性和适合性。

2. 目标函数:非线性规划的目标函数是优化问题的目标,通常是最大化或者最小化的数学表达式。

它可能包含非线性项,如幂函数、指数函数等。

3. 约束条件:非线性规划的约束条件是对决策变量的限制条件,用于定义可行解的集合。

约束条件可以是等式或者不等式,也可以包含非线性项。

4. 局部最优解与全局最优解:非线性规划问题可能存在多个极值点,其中局部最优解是在某一特定区域内最优的解,而全局最优解是在整个可行域内最优的解。

二、求解方法1. 数学方法:非线性规划问题可以通过数学方法进行求解,如拉格朗日乘子法、KKT条件等。

这些方法基于数学推导和分析,可以得到问题的解析解。

2. 迭代方法:对于复杂的非线性规划问题,往往采用迭代方法进行求解。

典型的迭代方法包括牛顿法、拟牛顿法、单纯形法等。

这些方法通过不断迭代逼近最优解。

3. 优化软件:为了简化非线性规划问题的求解过程,可以使用专门的优化软件,如MATLAB、Gurobi、CPLEX等。

这些软件提供了丰富的求解算法和工具,可快速求解复杂问题。

三、应用案例假设你是一家创造公司的生产经理,你需要确定每一个产品的生产数量,以最大化公司的利润。

公司生产两种产品:A和B。

每一个产品的生产成本和利润如下:产品A:生产成本为1000元/件,利润为300元/件。

产品B:生产成本为1500元/件,利润为500元/件。

公司的生产能力有限,每天最多只能生产1000件产品。

此外,由于市场需求的限制,产品A和B的销售量之和不能超过800件。

你的任务是确定每一个产品的生产数量,以最大化公司的利润。

第6讲非线性规划基本知识

第6讲非线性规划基本知识
2 1
2 5 5 5 1 1 5 5 5
2 5 1 5
5 5
2 4 x1 x2 x2 | X1
26 2 5 0.7 5
f
X 3x
(0)
2 1
2 4 x1 x2 x2 | X ( 0 ) 1
第4讲 非线性规划及一维搜索(第3章)
非线性规划基本概念 一维搜索
(1学时) (1学时)
重 点:下降迭代算法、黄金分割法、二次插值法。 难 点:下降迭代算法构造 基本要求:了解非线性规划的分类,掌握梯度的计算和性质, 会用海赛阵判断凸规划,掌握用黄金分割法、二次插值法。
非线性规划基本概念(3.1)
1 2 n
Hesse矩阵
x 2 2 f (1 X ) x2x1 H ( X ) (f ( X )) 2 f ( X ) xnx1
2 f ( X )
2 f ( X ) x1x2 2 f ( X )
2 x2

2 f ( X ) xn x2
( )
f ( X ( k 1) )
( x1 k 1)
p
(k ) 1

f ( X ( k 1) ) x ( k 1) 2
p
(k ) 2
...
f ( X ( k 1) ) x ( k 1) n
( pnk )
f ( X ( k ) P ( k ) )T P ( k ) 0
令 k k 1 依次用(1)式计算得一个点列 X ( k ) 若有: f ( X (0) ) f ( X (1) ) f ( X (2) ) ... f ( X ( k ) ) ... 则称(1)为下降迭代算法

运筹学-非线性规划(四)(名校讲义)

运筹学-非线性规划(四)(名校讲义)

1.外点法(又称惩罚法)
其思路是在目标函数中增加一项使之变为无约束问题,同时 对破坏约束项需付出高昂代价,该法的起始点在可行域外, 一旦进入可行域内便得到最优解。
§1 多维有约束寻优方法 (9)
①思路 设原问题为min:f[X] 约束:XS S是En中一个约束子集,即X的可行域为S,
则将其变成无约束问题: min:f[X]+ p(X)
p(x) =1 =10 =100
=100
b
=10
=1
图4-19
x
a
§1 多维有约束寻优方法 (11)
②求解过程
令{k }(k=1,2,…,)是一无穷序列,且k≥0,k+1 >k,定义函数q(,X)=f(X)+Xk,若原问题有解,则当k→∞,
§1 多维有约束寻优方法 (2)
一、库恩-塔克(简称库塔)条件
1.可行方向和起作用约束 ①可行方向:
设X(0)是可行点,即X(0) R,若对于某一方向D,存在一 个数 0>0,使对于任意 (0≤≤0 )均有下式成立:X(0) +DR,则称方向D是点X(0)处的可行方向。 ②下降方向:对于f(X)的台劳级数展开,若▽f[X(0)]T· D<0, 则称D方向为f[X]的下降方向。
第二十四讲 非线性规划(四)
§1 多维有约束寻优方法
§1 多维有约束寻优方法 (1)
非线性规划的一般形式
min f(X) hi(X)=0 i=1,2,…,m (1)
gi(X)≥0 j=1,2,…,l
下面,先阐述非线性规划的重要理论成果——库恩-塔克 条件(Kuhn-Tucker),然后介绍比较重要的几种有约束的寻 优方法。
0 1 2 x1

《非线性规划》课件

《非线性规划》课件
非线性规划的优化目标是找到使目标函数达到最大值或最小值的最优解。这些目标可以是经济、社会或 科学领域中的实际问题。
非线性规划的约束条件
非线性规划的约束条件是指限制问题解的一组方程或不等式。这些约束条件可以包括物理限制、资源约 束和行为限制等。
非线性规划的求解方法
线性化方法
将非线性问题转化为等价的 线性问题,然后使用线性规 划方法求解。
牛顿法
使用牛顿迭代法逐步逼近最 优解。
拟牛顿法
使用近似Hessian矩阵的方法 优化牛顿法。
变尺度法、全局优化方法
1
变尺度法
通过改变尺度,将问题转化为更易求解的形式。
2
全局优化方法
使用启发式算法寻找全局最优解。
非线性规划的应用领域
生产计划问题
优化生产计划,提高效率和利润。
交通运输问题
优化交通网络和运输流程。
优化电力系统
使电力系统运行更加高效和可靠。
决策支持系统
为决策者提供优化建议和决策支持。
医资源分配和治疗方案。
非线性规划的挑战
复杂的问题结构和求解困难。
未来的研究方向
未来的研究方向包括改进算法性能、适用于大规模问题的方法和考虑不确定性的优化模型等。
《非线性规划》PPT课件
在这个《非线性规划》PPT课件中,我们将深入探讨非线性规划的各个方面, 并介绍其在不同领域的应用。让我们一起开启这个激动人心的学习之旅!
什么是非线性规划?
非线性规划是一种在优化问题中寻找最优解的数学方法。它处理的是有非线 性约束条件和目标函数的优化问题。
非线性规划的优化目标

非线性规划-KT条件

非线性规划-KT条件

第三讲非线性规划§4约束极值问题(1)问题min(),{|()0,1,}jf XR X g X j l⎧⎨=≥=⎩<1>思路:有约束→无约束; 非线性→线性; 复杂→简;一、最优性条件1. 可行下降方向(有用约束,可行方向,下降方向) (1) 有用(效)约束设<1>式的(),()jf Xg X有一阶连续偏导设(0)X是一个可行解, 下一步考察时,要讨论约束.分析: 应有(0)(0)(0)()0()0()0j j j g X g X g X ⎧>⎪≥→⎨=⎪⎩ 若(0)()0j g X>,则在(0)()U X 内, 有()0j g X >, 此时各个方向均可选. 若(0)()0j g X =,则(0)X∈()0j g X =形成的边界, 影响下一步选向.1x 2x {()0}R X g X =≥()f X ()0j g X =(0)X故称()0j g X =是(0)X 点的有效约束.(2) 可行方向(对可行域来说)设(0)X为可行点, P 为某方向,若存在00λ>, 使得(0)0,[0,]X P R λλλ+∈∈ 则称P 是(0)X点的一个可行方向.(a) 可行方向P 与有效约束(0)()0j g X =的梯度(0)()j g X ∇关系是:(0)()0T j g X P ∇≥.记有效约束下标集(0){|()0,1}j J j g X j l ==≤≤ 若P 为(0)X的可行方向, 则存在00λ>, 使得当0[0,]λλ∈,有(0)(0)()()0,j j g X P g X j J λ+≥=∈从而(0)(0)0d ()()0,d j T j g X P g X P j J λλλ=+=∇≥∈见下图.(b)反之, 若(0)()0Tj g XP ∇>, 则P 必为可行方向.(0)(0)(0)()()()()T j j j g X P g X g X P o λλλ+=+∇+<1>对有效约束(0)()0j g X=,只要λ充分小,得(0)1()0g X =(0)2()0g X =(0)2()g X ∇(0)1()g X ∇P(0)X ∙(0)()0j g X P λ+≥, 所以P 是可行方向;<2>对无效约束(0)()0j g X >,同样只要λ充分小, 就有(0)()0j g X P λ+≥,故P 也是可行方向; 事实上, 对无效(0)()0j g X>,P ∀都是可行方向.(3) 下降方向(对目标函数来说) 设(0)XR ∈, 对某P 方向, 若在00[0,],0λλλ''∈>内, 有(0)(0)()()f X P f X λ+<则称P 是一个下降方向. 下降方向判定:若(0)()0T f X P ∇<,则P 是(0)X 的一个下降方向.因为(0)()()f X f X P λ=+(0)(0)()()()T f X f X P o λλ=+∇+,只要λ充分小, 都有(0)()()f X f X <.(4) 可行下降方向 若(0)XR ∈的某方向P 是可行方向+下降方向,则称P 是(0)X的可行下降方向.即 存在00λ>,当0[0,]λλ∈时,有(0)()0j g X P λ+≥且(0)(0)()()f X P f X λ+<,是继续寻优方向. 讨论: (0)X非极小值点⇔存在可行下降方向P ; (0)X极小值点⇔无可行下降方向P ;(可行但不下降,或下降不可行)定理(局部极(最)小必要条件)设X *是min (),{()0}i f X X g X ∈≥局部极小点,(),(),j f X g X j J ∈(有效约束下标集)在X *处可微, (),j g X j J∉在X *处连续, 则在X *处无可行下降方向P ,即不存在P , 使**()0,,()0,T j Tg X P j J f X P ⎧∇>∈⎨∇<⎩(**) 证 否则由(**)及前面的分析, 可找出可行下降点 →X *非局部极小值点→矛盾.如图 所示问题:min (),{|()0,1,}j f X R X g X j l ⎧⎨=≥=⎩<1>2. 库恩—塔克条件(局部最小的必要条件) 是非线性规划中最重要成果之一 (1) Gordan 引理(不加证明) 设12,,...,l A A A 是l 个n 维向量, 则1x 2x ()f X *∇()f X 1()g X *∇*P ∃/,使0,1,2,...,T j A P j l <=⇔ 0j μ∃≥,不全为零, 使10lj j j A μ==∑.(不指向同侧的向量, 正组合为零)(如l =3,n =2)若同侧, 则有P (图a), 否则无P (图b),但可正组为0.3A 1A 2A PH ()a 3A 1A 2A P H()b(2) Fritz John 定理设X *是<1>极小值点, ()f X 和()j g X 有一阶连续偏导数, 则存在不全为零的01,,...,l μμμ, 使⎧⎪⎨⎪⎩01()()0()0,1,2,...,0,1,2,...,lj j j j j j f X g X g X j l j lμμμμ**=*∇-∇===≥=∑证明 因X *是问题<1>的解, 故由定理4, 不存在 可行下降方向P, 使()0()0,TTj f X P g X P j J **⎧∇<⎪⎨-∇<∈⎪⎩由Gordan 引理,存在不全为零非负数0,,j j J μμ∈使0()()0j j j Jf Xg X μμ**∈∇-∇=∑对无效约束j J ∉, 令0j μ=, 则()0j j g X μ*∇=从而有(对所有l )01()()0lj j j f X g X μμ**=∇-∇=∑且有()0,0,1,2,.j j j g X j l μμ*=≥=, 证毕.注1: 类似于条件极值的必要条件.注2 若00μ=,则有效约束的()j g X *∇正线性相关→同侧→有可行下降方向→X *非极值点. 故一般设()j g X *∇线性无关→00μ>. 以上条件称为 Fritz John 条件, *X 称为Fritz John 点.(3) 必要条件 (库恩-塔克条件)设X *是<1>极小值点, ()f X 和()j g X 有一阶连续偏导,且有效约束梯度线性无关,则1,...,l μμ**∃, 使⎧⎪⎨⎪⎩1()()0()0,1,2,...,0,1,2,...,lj j j j j j f X g X g X j l j lμμμ***=***∇-∇===≥=∑<2> 证明 由Fritz John 引理, ()j g X *∇j J ∈线性无关得00μ>, 作00/0j μμμ*=>, 即得<2>. 式<2>=库恩-塔克条件. 相应点=库恩-塔克点. 简称K-T 条件, K-T 点. 对一般非线性规划min (),()0,1,()0,1,i j f X h X i m g X j l⎧⎪==⇒⎨⎪≥=⎩ min (),()0,()0,1,()0,1,i ij f X h X h X i m g X j l⎧≥⎪⎨-≥=⎪≥=⎩<3> 它的K-T 条件如下设X *是<3>极小值点, 相应函数有一阶连续偏导,且有效约束的()i h X *∇和(),j g X j J *∇∈线性无关,则12(,,...,)Tm Γγγγ****∃=和1(,...,)T l M μμ***=,使11()()()0()0,1,2,...,0,1,2,...,mlii j j i j j j j f X h X g X g X j lj lγμμμ*****==***∇-∇-∇===≥=∑∑<4>其中12,,...,m γγγ***,1,...,l μμ**称为广义Lagrange 乘子.注1 对凸规划, K-T 条件也是充分的.设kX 为某可行解, 若kX 是极小点, 且1()0k g X =, 和2()0k g X =, 则()()k f X ∇必与,1()k g X ∇和2()k g X ∇同侧, 否则有可行下降方向.由1()k g X ∇和2()k g X ∇线性无关1122()()()k k f X g X g X μμ*∇=∇+∇即1x 2x ()k f X -∇()f X 1()0k g X =2()0k g X =2()k g X ∇1()k g X ∇1122()()()0k k f X g X g X μμ*∇-∇-∇=例10 用库恩-塔克条件解非线性规划2max ()(4)16f x x x ⎧=-⎨≤≤⎩. 1()g X ∇()k X ()()k f X -∇86()1()0k g X =()2()0k g X =2()g X ∇R164解 变为 212min ()(4)()10()60f x x g x x g x x ⎧=--⎪=-≥⎨⎪=-≥⎩, 12()2(4),()1,()1f x x g x g x ∇=--∇=∇=-,引入广义拉格朗日乘子12,μμ**, 则有 所以1212122(4)0(1)0(6)0,0x x x μμμμμμ*********⎧---+=⎪-=⎪⎨-=⎪⎪≥⎩, 具体分析如下.若120,0,μμ**>>引出矛盾, 无解;若120,0μμ**>=:1x *=,点; ()9f x *=(1μ*=6) 若120,0μμ**==:4x *=,()0f x *=;若120,0μμ**=>:6x *=,()4f x *=(2μ*=4) 所以最大值点1x *=, 最大值()9f x *=. 注: 2()(4)f x x =--非凸函数, 在[1,6]上有两个局部最小值点.还有一个”驻点”164附加例题(略)用K-T 条件解非线性规划2min ()(3)05f x x x ⎧=-⎨≤≤⎩. 解 212min ()(3),()0,()50f x x g x x g x x ⎧=-⎪=≥⎨⎪=-≥⎩,(是凸规划) 12()2(3),()1,()1f x x g x g x ∇=-∇=∇=-,所以1212122(3)00(5)0,0x x x μμμμμμ*********⎧--+=⎪=⎪⎨-=⎪⎪≥⎩, 具体分析如下. 若120,0,μμ**≠≠引出矛盾, 无解;若120,0,μμ**≠=解得10,6x μ**==-,非K-T 点; 若120,0,μμ**=≠解得15,4x μ**==-,非K-T 点;若120,0,μμ**==解得3x *=,()0f x *=全局最小.习题4.1 已知非线性规划131212max ()(1)0,0f X x x x x x =⎧⎪--≥⎨⎪≥⎩的极大点为(1,0), 试(1) 转化目标函数后, 写出其K-T 条件;(2) 求出K-T 点..4.2 试用K-T 条件求解问题2min ()(4)16f X x x =-≤≤.。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
g1 ( x ) x1 x2 4,
g1 ( x) [ 1 , 1 ]T
g2 ( x) x1 ,
g2 ( x) [ 1 , 0 ]T 。
g3 ( x) x2 ,
g3 ( x) [ 0 , 1 ]T 。
18
由K T条件得
x1 3 1 1 0 x 3 1 1 2 0 3 1 0 2
分析:
(1) 如果 I ( x*)中只有一个指标,不妨 设 g1 ( x)为积极约束。
则不存在向量d 使得 g1 ( x*)T d 0 T f ( x *) d0 成立。
12
则不存在向量d 使得 g1 ( x*)T d 0 成立。 T f ( x*) d 0
令 Q { x | h( x ) 0 , g ( x ) 0 } , 称 Q 为此约束极值问题的
可行域。
2
min f ( x ) hi ( x ) 0 i 1 , 2 , , m s.t. g j ( x ) 0 j 1 , 2 ,, l
hi ( x ) 0 hi ( x ) 0 hi ( x ) 0
gi ( x ) ( i I ( x*) ) 在 点 x * 处 连 续, { gi ( x*)| i I ( x*) } 线性无关。若 x *是约束极值问题 (1)的 局 部 极 小 点 , 则存在一组实数 i 使 其 满 足
l f ( x*) i gi ( x*) 0 i 1 () i gi ( x*) 0 , i 0, i 1 , 2 , , l
8
2 g3 ( x ) 0。 2
I ( x ) { 1 , 2 }。
x2 g2 ( x ) 0
g3 ( x ) 0
O
g1 ( x ) 0
x
x1
②如何判断一个方向是可行方向?
9
定理1:
给 定 点x Q , 记 点 x 的 积 极 约 束 指 标 集 为 I ( x )。 给 定 向 量 d , 如果对任意的 i I ( x ) 有 gi ( x )T d 0 , 则 d 是 点 x 的 可 行 方 向 。
min f ( x ) s.t. g( x ) 0
(1)
g1 ( x ) 0
x0 d2 d1
可行域为 Q { x | g( x ) 0 }。
①可行方向与积极约束: 可行方向:
x1 d 1 d2
g2 ( x ) 0
设 x 0 Q, d 为 一 个 向 量 。 如 果 存实 在 数 0, 使得对任意的 [ 0 , ] 有 x 0 d Q , 则 称d 为 x 0 处 的 一个可行方向。
这 里x *最 优 , f ( x * )与 h( x * )共 线 , 而x' 非 最 优
h( x* )
f ( x' )与h( x' )不 共 线 ,
h( x' )
最优性条件即:
f ( x*) * j h j ( x*)
j 1
6
l
(2)不等式约束极值问题的最优性条件
可行下降方向:
设点 x Q , 给定向量d ,如果d 既是点 x 处的可行方向, 又是该点的下降方向, 则称 d 为点 x 处的可行下降方向。
10
I ( x )。 给 定 定理2: 给 定 点x Q , 记 点 x 的 积 极 约 束 指 标 集 为 向 量d , 如 果d 满 足 g i ( x )T d 0 T f ( x ) d 0 i I ( x)
7
g i ( x ) 0, 如 果 积极约束: 设 点 x Q , 对 于 不 等 式 约 束 gi ( x ) 0 , 则 称 gi ( x ) 0 是 点 x 处 的 积 极 约 束 。
或 起作用约束(紧约束\积极约束\有效约束)。
记 I ( x ) { i | gi ( x ) 0 , 1 i l } , 称 I ( x )为点 x 处的积极 约束指标集。
() 式称为 K T条件(库恩 塔克条件),满足 () 式的点 称为K T点。
16
(4) 对于有等式约束的极值 问题
min f ( x ) h( x ) 0 s.t. g( x ) 0
K T条件可写为
m l f ( x*) u j h j ( x*) i g i ( x*) 0 j 1 i 1 i gi ( x*) 0 , i 0, i 1 , 2 , , l
以下分情况讨论:
19
(1) 若 x1 x 2 0 :
由1 ( x1 x2 4) 0 可得1 0。
1 2 3 2 3
这与 2 0 矛盾。 ( 2) 若 x1 0 , x 2 0 : 3 0
x1 1 2 3 x 3 1 3 2 1 ( x1 x 2 4) 0 2 x1 0 3 x2 0 x1 x 2 4 0 1 , 2 , 3 , x1 , x 2 0
证明: 令 x' x t d , t 0。 则对任意的i I ( x ) , 有
gi ( x' ) gi ( x ) t gi ( x )T d o( || td ||2 ) t gi ( x )T d o( || td ||2 ) 0
x' Q , 即 d 为可行方向。
约束极值问题也可记为
min f ( x ) s.t. g( x ) 0
3
2 约束极值及最优性条件——Kuhn-Tucker 条件
(1)等式约束性问题的最优性条件 考虑 min f(x) s.t. h(x)=0 回顾高等数学中所学的条件极值: 问题 即: 求 z = f(x,y)极值,在ф(x,y)=0的条件下。 min f(x,y)
记 h( x ) ( h1 ( x ) , h2 ( x ) ,, hm ( x ) )T , g( x ) ( g1 ( x ) , g2 ( x ) ,, gl ( x ) )T ,
则约束极值问题可记为 min f ( x )
h( x ) 0 s .t . g( x ) 0
由K T条件及约束条件得
x1 1 2 3 x 3 1 3 2 1 ( x1 x 2 4) 0 2 x1 0 3 x2 0 x1 x 2 4 0 1 , 2 , 3 , x1 , x 2 0
( 3)
15
i gi ( x*) 0 , i 1 , 2 , , l i 0, i 1 , 2 , , l
i 0 , gi ( x*) 0 ; i 0 , gi ( x*) 0 ;
定理4(K-T条件): 设 x* Q,f ( x ) 和 gi ( x ) ( i I ( x*) ) 在x * 处 可 微 ,
g2 ( x*) g1 ( x ) 0
g1 ( x*)
g2 ( x ) 0
x*
f ( x*)
则 存在1 , 2 0 , 使得 f ( x*) 1g1 ( x*) 2g2 ( x*)。
14

f ( x*) 1g1 ( x*) 2g2 ( x*) 0。
11
3. K T 条件(库恩 塔克条件)
设I ( x*) 是 其 积 极 约 束 指 标 , 集 点 x *是约束极值问题 (1) 的局部极小点 ,则在点 x*处不存在下降方向 d。
则由定理 2可 知 , 不 存 在 向 量 d, 使 下 式 成 立 g i ( x*)T d 0 T f ( x *) d0 i I ( x*) 。
若x*是其的最优解 , 则存在υ*∈ Rl 使
f ( x * )
* h ( x )0 * j j j 1
5
l
f ( x )
*
* * h ( x j j )0 j 1
l
几何意义:考虑一个约束的情况:
f ( x* )
f ( x' )
x'
h( x )
g1 ( x*) g1 ( x ) 0
x*
f ( x*)
则有
f ( x*) g1 ( x*) , 0。

f ( x*) g1 ( x*) 0。
13
( 2) 如果 I ( x*)中有两个指标,不妨设g1 ( x )和 g2 ( x )为积极约束。 并设g1 ( x*)和g2 ( x*)线性无关。
(3) 一般情况:设{gi ( x*) | i I ( x*) } 线性无关。 则存在非负实数i ( i I ( x*) ), 使得
f ( x*)
iI ( x*)
g ( x*) 0
i i
( 2)
( 2) 式可改写为
l f ( x*) i gi ( x*) 0 i 1 i gi ( x*) 0 , i 0, i 1 , 2 , , l
则 向 量d 是 点 x 处 的 可 行 下 降 方 向 。
证略
③极值点的必要条件: 定理3:
设 x* Q, I ( x*)是其积极约束指标集。
f ( x) 和 gi ( x) (i I ( x*)) 在点x * 处可微,
gi ( x ) (i I ( x*) ) 在点x * 处连续。
如 果 x * 是约 束极 值问 题 (1)的局 部极 小点 ,则在 点 x * 处没 有可 行下 降方向 。
()
17
相关文档
最新文档