拉格朗日乘数法
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
§4 条件极值
(一) 教学目的:了解拉格朗日乘数法,学会用拉格朗日乘数法求条件极值.
(二) 教学内容:条件极值;拉格朗日乘数法.
基本要求:
(1)了解拉格朗日乘数法的证明,掌握用拉格朗日乘数法求条件极值的方法.
(2) 较高要求:用条件极值的方法证明或构造不等式.
(三) 教学建议:
(1) 本节的重点是用拉格朗日乘数法求条件极值.要求学生熟练掌握.
(2) 多个条件的的条件极值问题,计算量较大,可布置少量习题.
(3) 在解决很多问题中,用条件极值的方法证明或构造不等式,是个好方法.可推荐给
较好学生.
——————————————————————
在许多极值问题中,函数的自变量往往要受到一些条件的限制,比如,要设计一个容积为V 的长方体形开口水箱,确定长、宽和高, 使水箱的表面积最小. 设水箱的长、宽、高分别为 z y x ,,, 则水箱容积 xyz V =
焊制水箱用去的钢板面积为 xy yz xz z y x S ++=)(2),,(
这实际上是求函数 ),,(z y x S 在 xyz V = 限制下的最小值问题。
这类附有条件限制的极值问题称为条件极值问题, 其一般形式是在条件
)(,,,2,1,0),,,(21n m m k x x x n k <== ϕ 限制下,求函数 ),,,(21n x x x f 的极值
条件极值与无条件极值的区别
条件极值是限制在一个子流形上的极值,条件极值存在时无条件极值不一定存在,即使存在二者也不一定相等。
例如,求马鞍面 12
2+-=y x z 被平面 XOZ 平面所截的曲线上的最低点。请看这个问题的几何图形(x31马鞍面)
从其几何图形可以看出整个马鞍面没有极值点,但限制在马鞍面被平面 XOZ 平面所截的曲线上,有极小值 1,这个极小值就称为条件极值。
二. 条件极值点的必要条件
设在约束条件0),(=y x ϕ之下求函数=z ),(y x f 的极值 . 当满足约束条件 的点),(00y x 是函数),(y x f 的条件极值点 , 且在该点函数),(y x ϕ满足隐函数存在条件时, 由方程0),(=y x ϕ决定隐函数)(x g y =, 于是点0x 就是一元函数())( , x g x f z =的极限点 , 有
0)(='+=x g f f dx dz y x . 代入 )
,(),()(00000y x y x x g y x ϕϕ-=', 就有 0)
,(),(),(),(00000000=-y x y x y x f y x f y x y x ϕϕ, ( 以下x f 、y f 、x ϕ、y ϕ均表示相应偏导数在点),(00y x 的值 . )
即 x f y ϕ—y f x ϕ0= , 亦即 (x f , y f ) (⋅y ϕ ,x ϕ-)0= .
可见向量(x f , y f )与向量(
y ϕ , x ϕ-)正交. 注意到向量(x ϕ , y ϕ)也与向量(y ϕ , x ϕ-)正交, 即得向量(x f , y f )与向量(x ϕ , y ϕ)线性相关, 即存在实数λ, 使
(x f ,y f ) + λ(x ϕ,y ϕ)0=.
亦即 ⎩⎨⎧=+=+. 0 , 0y y
x x f f λϕλϕ
Lagrange 乘数法 :
由上述讨论可见 , 函数=z ),(y x f 在约束条件0),(=y x ϕ之下的条件极值点应是方程组
⎪⎩⎪⎨⎧==+=+.
0),(, 0),(),(, 0),(),(y x y x y x f y x y x f y y x x ϕλϕλϕ
的解.
引进所谓Lagrange 函数
),(),(),,(y x y x f y x L λϕλ+=, ( 称其中的实数λ为Lagrange 乘数 )
则上述方程组即为方程组
⎪⎩⎪⎨⎧===.
0),,( , 0),,( , 0),,(λλλλy x L y x L y x L y x
因此,解决条件极值通常有两种方法
1)直接的方法是从方程组(1)
,
,,2,1,0),,,(21m k x x x n k ==ϕ
中解出 m x x x ,,,21 并将其表示为 m k x x x g x n m m k k ,,2,1,),,,(21 ==++
代入 ),,,(21n x x x f 消去 m x x x ,,,21 成为变量为 n m x x ,,1 +的函数 ),,(),,,,,(),,(1111n m n m m n x x F x x g g f x x f ++==
将问题化为函数 ),,(1n m x x F + 的无条件极值问题;
2)在一般情形下,要从方程组(1)中解出 m x x x ,,,21 来是困难的,甚至是不可能的,因此上面求解方法往往是行不通的。通常采用的拉格朗日乘数法,是免去解方程组(1)的困难,将求 ),(1n x x f 的条件极值问题化为求下面拉格朗日函数
∑=+=m k n k k
n m n x x x x f x x L 11111),,(),,(),,;,,( ϕλλλ
的稳定点问题,然后根据所讨论的实际问题的特性判断出哪些稳定点是所求的极值的。
一. 用Lagrange 乘数法解应用问题举例 :
例1 抛物面z y x =+22被平面1=++z y x 截成一个椭圆. 求该椭圆到坐标 原点的最长和最短距离.
例3求函数xyz z y x f =),,( 在条件
)0,0,0,0( 1111>>>>=++r z y x r
z y x 下的极小值. 并证明不等式 31
1113abc c b a ≤⎪⎭
⎫ ⎝⎛++-, 其中 c b a , , 为任意正常数 .
现在就以上面水箱设计为例, 看一看拉格朗日乘数法求解条件极值的过程 解: 这个问题的实质是求函数 xy yz xz z y x S ++=)(2),,(
在条件 0=-V xyz
下的最小值问题, 应用拉格朗日乘法,令
L='2*(x*z+y*z)+x*y+v*(x*y*z-V)';
dLdx=diff(L,'x')
dLdy=diff(L,'y')
dLdz=diff(L,'z')
dLdv=diff(L,'v')
dLdx =2*z+y+v*y*z
dLdy =2*z+x+v*x*z
dLdz =2*x+2*y+v*x*y
dLdv =x*y*z-V
令 L 的各偏导等零,解方程组求稳定点
s1='2*z+y+v*y*z';
s2='2*z+x+v*x*z';
s3='2*x+2*y+v*x*y';
s4='x*y*z-V';
[v,x0,y0,z0]=solve(s1,s2,s3,s4)
v =
[ -2*2^(2/3)/V^(1/3)]