§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 的极值条件极值与无条件极值的区别条件极值是限制在一个子流形上的极值,条件极值存在时无条件极值不一定存在,即使存在二者也不一定相等。
例如,求马鞍面 122+-=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 yx 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 kn 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 rz y x 下的极小值. 并证明不等式 311113abc 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*zdLdy =2*z+x+v*x*zdLdz =2*x+2*y+v*x*ydLdv =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)][ -8*(-1/4*2^(1/3)*V^(1/3)+1/4*i*3^(1/2)*2^(1/3)*V^(1/3))^2/V][ -8*(-1/4*2^(1/3)*V^(1/3)-1/4*i*3^(1/2)*2^(1/3)*V^(1/3))^2/V] x0 =[ 2^(1/3)*V^(1/3)]y0 =[ 2^(1/3)*V^(1/3)]z0 =[ 1/2*2^(1/3)*V^(1/3)]这里显然只有实数解才有意义, 所以 L 的稳定点只有下面一个33221,2V z V y x === 又已知所求的问题确实存在最小值,从而解出的稳定点就是最小值点, 即水箱长宽与为高的2倍时用钢板最省。
下面再看一个条件极值求解问题例2 抛物面 z y x =+22 被平面 1=++z y x 截成一个椭圆,求这个椭圆到坐标原点的最长最短距离。
(x73)解: 这个问题的实质是求函数 222),,(z y x z y x f ++=在条件 022=-+z y x 与 01=-++z y x 下的最大、最小值问题, 应用拉格朗日乘法,令L='x^2+y^2+z^2+v*(x^2+y^2-z)+h*(x+y+z-1)';dLdx=diff(L,'x')dLdy=diff(L,'y')dLdz=diff(L,'z')dLdv=diff(L,'v')dLdh=diff(L,'h')dLdx =2*x+2*v*x+hdLdy =2*y+2*v*y+hdLdz =2*z-v+hdLdv =x^2+y^2-zdLdh =x+y+z-1s1='2*x+2*v*x+h';s2='2*y+2*v*y+h';s3='2*z-v+h';s4='x^2+y^2-z';s5='x+y+z-1';[h,v,x0,y0,z0]=solve(s1,s2,s3,s4,s5);x0,y0,z0x0 =[ 3/4-1/4*i*13^(1/2)][ 3/4+1/4*i*13^(1/2)][ -1/2+1/2*3^(1/2)][ -1/2-1/2*3^(1/2)]y0 =[ 3/4+1/4*i*13^(1/2)][ 3/4-1/4*i*13^(1/2)][ -1/2+1/2*3^(1/2)][ -1/2-1/2*3^(1/2)]z0 = -1/2, -1/2, 2-3^(1/2), 2+3^(1/2)即 L 的稳定点有两个 32,23132,231222111+=--==-=+-==z y x z y x 因为函数 ),,(z y x f 在有界闭集 }1,|),,({22=++=+z y x z y x z y x 上连续,必有最大值和最小值,而求得的稳定点又恰是两个,所以它们一个是最大点, 另一个是最小,其最大最小值为。
(x73)x1=-1/2+1/2*3^(1/2);x2=-1/2-1/2*3^(1/2);y1=-1/2+1/2*3^(1/2);y2=-1/2-1/2*3^(1/2);z1=2-3^(1/2);z2=2+3^(1/2);f1=(x1^2+y1^2+z1^2)^(1/2) f2=(x2^2+y2^2+z2^2)^(1/2)f1 = 0.5829 ; f2 = 4.2024。