非线性规划的概念和原理
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
第五章 非线性规划的概念和原理
非线性规划的理论是在线性规划的基础上发展起来的。1951年,库恩(H.W.Kuhn )和塔克(A.W.Tucker )等人提出了非线性规划的最优性条件,为它的发展奠定了基础。以后随着电子计算机的普遍使用,非线性规划的理论和方法有了很大的发展,其应用的领域也越来越广泛,特别是在军事,经济,管理,生产过程自动化,工程设计和产品优化设计等方面都有着重要的应用。
一般来说,解非线性规划问题要比求解线性规划问题困难得多,而且也不像线性规划那样有统一的数学模型及如单纯形法这一通用解法。非线性规划的各种算法大都有自己特定的适用范围。都有一定的局限性,到目前为止还没有适合于各种非线性规划问题的一般算法。这正是需要人们进一步研究的课题。
5.1 非线性规划的实例及数学模型
[例题6.1] 投资问题:
假定国家的下一个五年计划内用于发展某种工业的总投资为b 亿元,可供选择兴建的项目共有几个。已知第j 个项目的投资为j a 亿元,可得收益为j c 亿元,问应如何进行投资,才能使盈利率(即单位投资可得到的收益)为最高?
解:令决策变量为j x ,则j x 应满足条件()
10j j x x -= 同时j x 应满足约束条件
1
n
j j
j a x
b =≤∑
目标函数是要求盈利率()1121
,,,n
j
j
j n n
j j
j c x
f x x x a x
===
∑∑最大。
[例题6.2] 厂址选择问题:
设有n 个市场,第j 个市场位置为()
,j j p q ,它对某种货物的需要量为j b ()1,2,,j n =。
现计划建立m 个仓库,第i 个仓库的存储容量为i a ()1,2,,i m =。试确定仓库的位置,使
各仓库对各市场的运输量与路程乘积之和为最小。
解:设第i 个仓库的位置为(),i i x y ()1,2,,i m =,第i 个仓库到第j 个市场的货物供
应量为i j z ()1,2,
,,1,2,
,i m j n ==,则第i 个仓库到第j 个市场的距离为
i j d =
目标函数为
11
11
m
n
m
n
i j
i j i j
i j i j z
d z =====∑∑∑∑
约束条件为:
(1) 每个仓库向各市场提供的货物量之和不能超过它的存储容量; (2) 每个市场从各仓库得到的货物量之和应等于它的需要量; (3) 运输量不能为负数。
因此,问题的数学模型为:
11
min m n
i i j z ==∑∑s.t.
1
n
ij
i j z
a =≤∑,()1,2,,i m =
1
m
i j
j i z
b =≤∑,()1,2,
,j n =
0ij z ≥,()1,2,
,,1,2,
,i m j n ==
一般非线性规划的数学模型可表示为:
()min f x ;
s.t. ()0i g X ≥ ()1,2,,i m =, ()0j h X =,()1,2,
,j l =
式中()12,,
,T
n n X x x x R =∈是n 维向量,,i f g ()1,2,
,i m =,j h ()1,2,
,j l =都是
1n R R →的映射(即自变量是n 维向量,因变量是实数的函数关系),且其中至少存在一个
非线性映射。
与线性规划类似,把满足约束条件的解称为可行解。若记
()(){}
0,1,2,,,0,1,2,,i j X g X i m h X j l χ=≥===
则称χ为可行域。因此上述模型可简记为
()min f X
s.t. X χ∈
当一个非线性规划问题的自变量X 没有任何约束,或说可行域即是整个n 维向量空间,
即n
R χ=,则称这样的非线性规划问题为无约束问题:
()min f X 或()min n
X R
f X ∈ 有约束问题与无约束问题是非线性规划的两大类问题,它们在处理方法上有明显的不
同。
5.2 无约束非线性规划问题
5.2.1无约束极值条件
对于二阶可微的一元函数()f x ,如果x *
是局部极小点,则()
0f x *'=,并且
()0f x *''>;反之,如果()0f x *'=,()0f x *''<,则x *是局部极大点。关于多元函数,
也有与此类似的结果,这就是下述的各定理。
考虑无约束极值问题:
()min f x ,n x E ∈
定理6.1 (必要条件)设()f x 是n 元可微实函数,如果x *
是以上问题的局部极小解,
则()
0f x *∇=。
定理6.2 (充分条件)设()f x 是n 元二次可微实函数,如果x *
是上述问题的局部最小
解,则()
0f x *∇=,()
2f x *∇半正定;反之,如果在x *
点有()
0f x *∇=,()
2f x *∇正定,则x *
为严格局部最小解。
定理6.3 设()f x 是n 元可微凸函数,如果()
0f x *∇=,则x *
是上述问题的最小解。
[例题6.3] 试求二次函数()2
2
121122,282420f x x x x x x =-+-+的极小点。
解:由极值存在的必要条件求出稳定点:
1148f x x ∂=-∂, 22
44f x x ∂=-∂,则由()0f x ∇=得12x =,21x = 再用充分条件进行检验:
2214f x ∂=∂,2224f x ∂=∂,2212210f f x x x x ∂∂==∂∂∂∂,则由2
4004f ⎛⎫∇= ⎪⎝⎭
为正定矩阵得
极小点为()2,1T
x *
=。