1.一个约束条件#
二元函数z=f(x,y),目标函数,g(x,y)=0约束条件
设方程g(x,y)=0确定的隐函数y=y(x),gy(P0)=0,P0(x0,y0)
若极值点在点P0(x0,y0)取到,则z=f(x,y(x))在x=x0取到极值,于是
dxdz=fx+fydxdy=0
dxdg(x,y(x))=gx+gydxdy=0
∇f,∇g⊥(1,dxdy=T),∇f//∇g
因此,存在常数λ0,使得∇f(P0)=λ0∇g(P0),g(p0)=0(∗2)
引入变量λ与函数L(x,y,λ)=f(x,y)−λg(x,y)(∗3)
因为
LxLyLλ=fx−λgx=fy−λgy=−g所以(∗2)等价于∇L(P0)=0(∗4)
若P0(x0,y0)为条件极值,则(x0,y0,λ0)为L(x,y,λ)的驻点,该方法称为拉格朗日乘数法
若约束集{(x,y)∣g(x,y)=0}是有限闭集,则条件极值存在
设gx,gy在D上不全为零,求出拉格朗日函数(∗3)的全部驻点以及相应的函数值,再取最大或最小
ex.4 设x2+xy+y2=1,求2x+3y的最大值#
解
令L(x,y,λ)=2x+3y−λ(x2+xy+y2−1),由∇L=0,得
LxLyLλ=2−λ(2x+y)=0=3−λ(x+2y)=0=−(x2+xy+y2−1)=0x=21±1,y=21±4,2x+3y=±2114
三元函数,u=f(x,y,z),g(x,y,z)=0
若极值在点P0(x0,y0,z0)处取到,设gz=0,确定z=z(x,y),z0=z(x0,y0)
gx+gzzx=0,∇g⊥(1,0,zx)=v1
gy+gyzy=0,∇g⊥(0,1,zy)=v2
∇g//v1×v2=(−zx,−zy,1)=v3
因为u=f(x,y,z(x,y))在Q0(x0,y0)取极值,所以
uxuy=fx+fzzx=0,∇f⊥v1=fy+fzzy=0,∇f⊥v2说明∇f(p0)//∇g(P0)
令L(x,y,z,λ)=f(x,y,z)−λg(x,y,z)(6)
则条件极值=>,∇L=0即∇f−λ∇g=0,g=0
例4.8 设x,y,z>0,xyz=v,u=xy+2yz+2yz,求umin#
解:
令L(x,y,z,λ)=xy+2yz+2xz−λ(xyz−v),由∇L=0
LxLyLzLλ=y+2z−λ(yz)=0=x+2z−λ(xz)=0=2y+2x−λ(xy)=0=−(xyz−v)=0z=34v,u=12z2=12(4v)32
ex x,y,z>0,证明xyz≤(3x+y+z)3#
证明:
u=xyz,x+y+z=a>0
令L(x,y,z,λ)=xyz−λ(x+y+z−a),由∇L=0
LxLyLzLλ=yz−λ=0=zx−λ=0=xy−λ=0=−(x+y+z−a)=0x=y=z=3a,u=(3a)3
约束集S={(x,y,z)∣x+y+z=0,x,y,z≥0}
为△ABC,它是有界闭域u=xyz在S上连续,可取到最大值与最小值,在∂S上,u=0,
于是在大致在内部取到,因此最大值为(3a)3
2.两个约束条件#
u=f(x,y,z),g(z,y,z)=0,h(x,y,z)=0(8)
设方程组g(x,y,z)=0,h(x,y,z)=0(9)
确定隐函数组x=x(z),y=y(z)(10)
若u=f(x(z),y(z),z)取极值,则
dzdu=fxdzdx+fydzdy+fz=0=(fx,fy,fz)⋅(dzdx,dzdy,1),∇f⊥T将x=x(z),y=y(z)带入结束条件
gxdzdx+gydzdy+gz=0,∇g⊥T同理
hxdzdx+hydzdy+hz=0,∇h⊥T故存在λ,μ,使得
∇f+λ∇g+μ∇h=0
令L(x,y,z,λ,μ)=f(x,y,z)+λg(x,y,z)+μh(x,y,z)
ex.5 求u=x2+y2+z2在约束条件4x2+5y2+25z2=1,z=x+y 的极值#
解:
令L(x,y,z,λ,μ)=x2+y2+z2−λ(4x2+5y2+25z2−1)−μ(x+y−z)
由∇L=0
LxLyLzLλLμ=2x−21λx−μ=0=2y−52λy−μ=0=2z−252λz+μ=0=−(4x2+5y2+25z2−1)=0=−(x+y−z)=0(1)(2)(3)(4)(5)(1)×x+(2)×y+(3)×z:
2(x2+y2+z2)−λ(2x2+52y2+252z2)−μ(x+y−z0
结合(4)(5),x2+y2+z2=λ
由(1)(2)(3)
x=2−21λμ,y=2−52λμ,z=2−252λμ
带入(5),则λ=10,1775
Method 2#
把x,y,z,μ当做未知数,λ为系数,线性方程组有解
2−21λ00102−52λ01002−252λ−1−1−110=(2−21λ)(2−52λ+2−252λ)−(−1)(2−52λ)(2−252λ)=12−2598λ+12534λ23.一般形式#
uφk(X)=f(X)=0X=(x1,x2,...,xn)∈Rn1≤k≤m<nL(X,λ)=f(X)+λ1φ1(X)+...+λmφm(X),λ=(λ1,...,λm)∈Rm
由∇L=0,求出全部驻点(X0,λ0)
X0=(x0,1,x0,2,...,x0,n),λ0=(λ0,1,...,λ0,m)
并要求rank(∂xi∂φk)m×n=m
∂x1∂φ1dx1+...+∂xm∂φ1dxm+∂xm+1∂φ1dxm+1+...+∂xn∂φ1dxn=0...∂x1∂φmdx1+...+∂xm∂φmdxm+∂xm+1∂φmdxm+1+...+∂xn∂φmdxn=0设X=(xm+1,...,xn)且优越数据确定的隐函数组
x1=x1(X),...,xm=xm(X)
条件极值的充分条件#
对L的每个驻点(X0,λ0)
计算Φ(X)=L(X,λ)在点X0的Hesse矩阵,
H(X0)=(∂xi∂yj∂2L)n×n若H(X0)正(负)定在X0为条件极小(大)值点
证明:
注意到
∇Φ(X)=∇f(X)+λ0,1∇φ1(X)+...+λ0,m∇φm(X),∇Φ(X0)=0
定义约束集S={X∈Rn∣φk(X)=0,1≤k≤m}
在S上f(X)=Φ(X),设X,X+h∈S,则
f(X+h)−f(X)=Φ(X+h)−Φ(X)=∇Φ(X)h+21hTH(X)h+o(∣∣h∣∣2)=21hTH(X)h+o(∣∣h∣∣2)当H(X0)正定时,a=min∣∣h∣∣=1hTH(X0)h>0
此时hTH(X0)h≥a∣∣h∣∣2,于是
f(X0+h)−f(X0)>∣∣h∣∣2(fraca2+o(1))≥0
ex.6 g(x,y,z)=2x2+y2+z2+2xy−2x−2y−4z+4=0确定的隐函数z=z(x,y)的极值#
解:
uL(x,y,z,λ)LxLyLz=f(x,y,z)=z,g(x,y,z)=0=z+λg(x,y,z),∇L=0=λ(4x+2y−2)=0=λ(2y+2x−2)=0=1+λ(2z−4)=0x=0,y=0则
z1=1,λ1=21或z2=3,λ2=−21
Lxx=4λ
Lyy=2λ
Lzz=2λ
Lxy=2λ
Lxz=0
Lyz=0
H1=220120001正定,z1(0,1)=1为极小值
H2=−2−20−1−2000−1负定,z2(0,1)=3为极大值
蜂房问题#
容积相同,一个六棱柱用怎样的三个全等菱形为底,表面积最小