21-第二十一章 偏导数的应用
依赖于
被以下题目直接调用
正文部分
21.1.1 曲线的切向量、切线与法平面
设 f(x,y) 为可微函数,曲线方程为
z=f(x,y),cosαx−x0=sinαy−y0.
求曲线上点 M0(x0,y0,z0) 的切线与 xOy 平面所成角的正切。
解
令
F(x,y,z)=z−f(x,y),G(x,y,z)=sinα(x−x0)−cosα(y−y0).
则
∂(y,z)∂(F,G)=−fy−cosα10=cosα,∂(z,x)∂(F,G)=10−fxsinα=sinα,
∂(x,y)∂(F,G)=−fxsinα−fy−cosα=fxcosα+fysinα.
于是曲线上点 M0(x0,y0,z0) 的切线与 z 轴正方向所成夹角 φ 的余弦为
cosφ=1+[fx(x0,y0)cosα+fy(x0,y0)sinα]2fx(x0,y0)cosα+fy(x0,y0)sinα.
由此可计算出切线与 xOy 平面所成角的正切为
fx(x0,y0)cosα+fy(x0,y0)sinα.
21.1.2 曲面的法向量、法线与切平面
求曲面
x=u+v,y=u2+v2,z=u3+v3
的切平面当切点 M(u,v),u=v,趋于曲面的边界 u=v 上的点 M0(u0,u0) 时的极限位置。
解
先求 u=v 时的切平面方程。由于
∂(u,v)∂(y,z)=2u3u22v3v2=6uv(v−u),
∂(u,v)∂(z,x)=3u213v21=3(u2−v2),∂(u,v)∂(x,y)=12u12v=2(v−u),
于是过 (x,y,z) 点的切平面方程为
6uv(v−u)(X−(u+v))+3(u2−v2)(Y−(u2+v2))+2(v−u)(Z−(u3+v3))=0.
消去 v−u,并整理得到
6uvX−3(u+v)Y+2Z=3uv(u+v)−u3−v3.
令 (u,v)→(u0,u0),则有
6u02X−6u0Y+2Z=4u03.
将 (X,Y,Z) 改写为 (x,y,z),可见所求的极限位置为
3u02x−3u0y+z=2u03.
21.1.3 曲线的夹角、曲面的夹角
设
zxy=u,x2+z2+y2+z2=v,x2+z2−y2+z2=w
是三个分别以 u,v,w 为参数的单参数曲面族,证明:过同一点的三曲面族的三个曲面是两两正交的。
证
设三曲面交于点 (x,y,z),其法向量分别为
(zy,zx,−z2xy),
(x2+z2x,y2+z2y,x2+z2z+y2+z2z),
(x2+z2x,−y2+z2y,x2+z2z−y2+z2z).
由于这三个向量的两两内积均为 0,所以过 (x,y,z) 点的这三个曲面族的三个曲面是两两正交的。
21.1.4 练习题
在曲线 x=t, y=t2, z=t3 上求出一点,使该点的切线平行于平面 x+2y+z=4。
解
曲线的切向量为 (1,2t,3t2)。切线平行于平面 x+2y+z=4,等价于此向量与平面的法向量 (1,2,1) 正交,于是
(1,2,1)⋅(1,2t,3t2)=1+4t+3t2=0.
故 t=−1,−31,所求两点为 (−1,1,−1) 与 (−31,91,−271)。
证明:斜驶线
tan(4π+2ψ)=ekφ,k=常数,
与地球的每一条子午线相交成定角,其中 φ 为地球上点的经度,ψ 为地球上点的纬度。
解
取半径为 R 的球面参数方程
x=Rcosψcosφ,y=Rcosψsinφ,z=Rsinψ.
球面第一基本形式为 ds2=R2(dψ2+cos2ψdφ2)。由题设
lntan(4π+2ψ)=kφ,
对 φ 求导得 secψdφdψ=k,即
dψdφ=kcosψ1.
若 α 为曲线与子午线所成角,则
tanα=R∣dψ∣Rcosψ∣dφ∣=cosψdψdφ=∣k∣1,
为常数,所以它与每一条子午线都交成定角。
求曲线 x2+y2+z2=6, x+y+z=0 在点 (1,−2,1) 处的切线与法平面方程。
解
设
F=x2+y2+z2−6,G=x+y+z.
在点 (1,−2,1) 处,∇F=(2,−4,2),∇G=(1,1,1)。切向量可取
∇F×∇G=(−6,0,6)∼(1,0,−1).
故切线为
x=1+t,y=−2,z=1−t.
法平面与切向量正交,因此其方程为
(1,0,−1)⋅(x−1,y+2,z−1)=0,
即 x−z=0。
求曲面 z=arctanxy 在点 (1,1,4π) 处的法线与切平面方程。
解
设 F(x,y,z)=z−arctanxy,则
Fx=x2+y2y,Fy=−x2+y2x,Fz=1.
在 (1,1,4π) 处有 ∇F=(21,−21,1)∼(1,−1,2)。故法线为
1x−1=−1y−1=2z−π/4,
切平面为
21(x−1)−21(y−1)+z−4π=0,
即 x−y+2z−2π=0。
证明:曲面 xyz=a3 (a>0) 的每一个切平面与坐标面形成体积相同的四面体。
解
设 (x0,y0,z0) 为曲面 xyz=a3 上一点,则切平面方程为
y0z0(x−x0)+x0z0(y−y0)+x0y0(z−z0)=0.
由 x0y0z0=a3 化简得
x0x+y0y+z0z=3.
它与三坐标轴的截距分别是 3x0,3y0,3z0,因而所成四面体体积为
V=61⋅3x0⋅3y0⋅3z0=29a3,
与点的选取无关。
21.2.2 梯度
求 u=x+y+z 在沿 x2+y2+z2=1 上点的外法向的方向导数,并问在球面上何点该方向导数取:(1) 最大值;(2) 最小值;(3) 等于 0。
解
球面 x2+y2+z2=1 上点的外法向方向为 (x,y,z),函数 u=x+y+z 沿 (x,y,z) 方向的方向导数为
∂l∂u=∇u⋅(x,y,z)=∣∇u∣cosθ=3cosθ,
其中 θ 是向量 ∇u=(1,1,1) 与 (x,y,z) 之间的夹角。于是当 x=y=z>0 时,方向导数取最大值 3,此时
x=y=z=33;
当 x=y=z<0 时,方向导数取最小值 −3,此时
x=y=z=−33;
当 x+y+z=0 时,方向导数为 0。
若 R2 上的可微函数 f(x,y) 满足
xfx(x,y)+yfy(x,y)=0,
则 f(x,y) 恒为常数。
证 1
设 x=rcosθ, y=rsinθ,取 l=(cosθ,sinθ),由题设得
∂l∂f=cosθfx(rcosθ,rsinθ)+sinθfy(rcosθ,rsinθ)=0.
于是任意一从原点出发的射线上的方向导数为 0,因此 f(x,y) 在该射线上为常数,即
f(x,y)=f(0,0).
证 2
令 x=rcosθ, y=rsinθ,由题设得
∂r∂f=fxcosθ+fysinθ=0.
即 f(rcosθ,rsinθ) 与 r 无关。对任意固定的 θ,令 r→0+,得
f(x,y)=r→0+limf(rcosθ,rsinθ)=f(0,0).
证 3
由
xfx(x,y)+yfy(x,y)=0=0⋅f(x,y)
及齐次函数的充分必要条件知 f 为零次齐次函数,于是 ∀t∈R,
f(tx,ty)=t0f(x,y)=f(x,y).
于其中令 t→0,得
f(x,y)=f(0,0).
证 4
对任意固定的 (x,y)∈R2,定义一元函数
φ(t)=f(tx,ty),t≥0.
则
φ′(t)=xf1(tx,ty)+yf2(tx,ty).
于是当 t>0 时,由题设得
φ′(t)=t1{txf1(tx,ty)+tyf2(tx,ty)}=0.
所以 φ(t) 为常数。取 t=1,则 φ(t)≡φ(1),即
f(tx,ty)=f(x,y).
令 t→0,得 f(x,y)=f(0,0)。
21.2.3 练习题
求函数 u=x2+y2+z2 在沿椭球面
a2x2+b2y2+c2z2=1
上点 (x,y,z) 处的外法线方向的方向导数。
解
设椭球面为 F(x,y,z)=a2x2+b2y2+c2z2−1=0,则外法向单位向量为
n=∣∇F∣∇F=a4x2+b4y2+c4z2(a2x,b2y,c2z).
又 ∇u=(2x,2y,2z),故所求方向导数为
∂n∂u=∇u⋅n=a4x2+b4y2+c4z22(a2x2+b2y2+c2z2)=a4x2+b4y2+c4z22.
设 ∂n∂u 和 ∂n∂v 为函数 u 和 v 在沿曲面 F(x,y,z)=0 上的点 (x,y,z) 的法线方向 n 的方向导数,证明:
∂n∂(uv)=u∂n∂v+v∂n∂u.
解
法向方向导数就是算子 ∂n∂=n⋅∇。于是
∂n∂(uv)=n⋅∇(uv)=n⋅(u∇v+v∇u)=u∂n∂v+v∂n∂u.
设 f(t), u=u(x,y,z) 为可微函数,证明:∇f(u)=f′(u)∇u。
解
由链式法则,
∂x∂f(u)=f′(u)ux,∂y∂f(u)=f′(u)uy,∂z∂f(u)=f′(u)uz.
因而
∇f(u)=f′(u)(ux,uy,uz)=f′(u)∇u.
设 f(x,y) 在点 p0(x0,y0) 处可微,l1,l2,⋯,ln 为 n 个单位向量,相邻的两向量夹角为 n2π。证明:
i=1∑n∂li∂f(x0,y0)=0.
解
设 ∇f(x0,y0)=(A,B),并令 li=(cosθi,sinθi),其中 θi=θ0+n2(i−1)π。则
∂li∂f(x0,y0)=Acosθi+Bsinθi.
求和得
i=1∑n∂li∂f(x0,y0)=Ai=1∑ncosθi+Bi=1∑nsinθi=0,
因为这些方向在单位圆上等分一周,余弦和正弦之和都为零。
设
u=a2x2+b2y2−c2z2,a>b>c>0.
求在点 (0,0,0) 处函数 u 增加最快的方向。
解
在原点处 ∇u(0,0,0)=0,故需比较二次主部。沿单位向量 l=(α,β,γ) 取点 (x,y,z)=t(α,β,γ),则
u=t2(a2α2+b2β2−c2γ2),α2+β2+γ2=1.
由于 a>b>c>0,有 b21>a21>0>−c21,故括号中最大值为 b21,当且仅当 α=γ=0, β=±1 时取得。故增加最快的方向是 y 轴方向,即 l=(0,±1,0)。
设
f(x,y)=⎩⎨⎧x−y+x2+y2xy2,0,(x,y)=(0,0),(x,y)=(0,0).
证明:f(x,y) 在原点处连续,沿任何方向的方向导数存在,但不可微。
解
先证连续性。对 (x,y)=(0,0),有
x2+y2xy2≤∣x∣,
故
∣f(x,y)∣≤∣x∣+∣y∣+∣x∣→0,
所以 f 在原点连续。
对任意单位向量 l=(cosθ,sinθ),
tf(tcosθ,tsinθ)−f(0,0)=cosθ−sinθ+1cosθsin2θ
与 t 无关,故沿任意方向的方向导数都存在。
但若 f 在原点可微,则由 fx(0,0)=1、fy(0,0)=−1 可知其微分应为 x−y。于是余项应满足
f(x,y)−(x−y)=x2+y2xy2=o(x2+y2).
取 y=x,得
2∣x∣f(x,x)−(x−x)=2∣x∣x/2=22sgnx,
极限不存在,更不趋于 0。故 f 在原点不可微。
21.3.1 Taylor 公式
设
f(x,y)=⎩⎨⎧x2+y21−ex(x2+y2),0,(x,y)=(0,0),(x,y)=(0,0).
求 f(x,y) 在 (0,0) 的 4 阶 Taylor 多项式,并求出
∂x∂y∂2f(0,0),∂x4∂4f(0,0).
解
由于
ex(x2+y2)=1+x(x2+y2)+21x2(x2+y2)2+o([x(x2+y2)]2),
于是
x2+y21−ex(x2+y2)=−x−21x2(x2+y2)+o(x2(x2+y2)).
由 Taylor 展式的唯一性知 f(x,y) 的 4 阶 Taylor 展开式为
−x−21x4−21x2y2.
由此得
∂x∂y∂2f(0,0)=0,∂x4∂4f(0,0)=4!(−21)=−12.
设
u=a2x2+b2y2+c2z2,a>b>c>0.
求在 (0,0,0) 处函数增加最快的方向。
解
由于 ∇u(0,0,0) 为零向量,故不能应用梯度方向是函数增长最快的方向的性质,需另想办法。考虑沿某单位方向 l=(α,β,γ) 函数 u 的变化
u(tα,tβ,tγ)−u(0,0,0),
其中 t 是参数。令 φ(t)=u(tα,tβ,tγ),由 Taylor 展式
φ(t)−φ(0)=φ′(0)t+2φ′′(0)t2+o(t2)=(a2α2+b2β2+c2γ2)t2+o(t2).
由于 a>b>c 及 α2+β2+γ2=1,当 t>0 充分小时,沿方向 l=(0,0,±1) 时 φ(t)−φ(0) 最大,即函数 u 增加最快。
注
设二次连续可微函数 f 在点 p0 处的梯度为零向量,但 Hesse 矩阵 Q 是非零矩阵,则 f 增长最快的可能方向在 Q 的最大特征值对应的特征子空间中。如果最大特征值对应的特征子空间为一维空间,则相应的方向就是 f 增长最快的方向。
设 f(x,y) 在单位圆盘
D={(x,y)∣x2+y2≤1}
上具有连续的一阶偏导数,且满足 ∣f(x,y)∣≤1, ∀(x,y)∈D。证明:存在点
(x0,y0)∈intD={(x,y)∣x2+y2<1},
使得在点 (x0,y0) 处有不等式 fx2+fy2≤16。
证
设
g(x,y)=f(x,y)+2(x2+y2).
在单位圆周上显然有 g(x,y)≥1,而在原点 g(0,0)≤1。所以或者 g 在 D 上恒等于 1,或者 g(x,y) 必在 D 的某个内点 (x0,y0) 处取到极小值。因此总存在 (x0,y0)∈intD,使得
∂x∂g(x0,y0)=∂y∂g(x0,y0)=0.
于是
fx(x0,y0)=−4x0,fy(x0,y0)=−4y0,
从而
(fx2+fy2)(x0,y0)=16(x02+y02)≤16.
证明:当 ∣x∣ 和 ∣y∣ 充分小时,有近似式
cosycosx≈1−21x2+21y2.
证
设 f(x,y)=cosycosx,则 f(x,y) 在原点附近无穷次可微,且
f(0,0)=1,fx(0,0)=0,fy(0,0)=0,
fxx(0,0)=−1,fyy(0,0)=1,fxy(0,0)=0.
由 Taylor 公式,当 ∣x∣ 和 ∣y∣ 充分小时,有
cosycosx=1−21x2+21y2+o(∣x∣2+∣y∣2),
此即所证近似式。
21.3.2 极值问题
求
z=21x2+xy+21y2−2x−2y+5
的全部极值点与极值。
解
由于
zx=zy=x+y−2,
则直线 x+y−2=0 上的点都是驻点。又
A=zxx=1,B=zxy=1,C=zyy=1,
于是 Δ=AC−B2=0,故不能用法则来判别这些驻点是否为极值点。从函数本身来看有
z=21(x+y)2−2(x+y)+5=21[(x+y)−2]2+3.
故 x+y=2 上的全部点均为 z 的极小值点,极小值为 3。由于 z 在 R2 上可微,故无其他极值点。
注
本题表明,对于二元函数,即使它不是常值函数,其驻点仍可以有无穷多个,且可以构成曲线。
设 f(x,y)=(y−x2)(y−2x2)。证明:沿着经过点 (0,0) 的每一条直线,点 (0,0) 均是 f(x,y) 在该直线上的极小值点。但点 (0,0) 不是 f(x,y) 在整体上的极小值点。
证
在 xOy 平面上画出 S1:y=x2 和 S2:y=2x2 表示的两条抛物线,可见 S1 和 S2 把平面分成四个区域。在每个区域 f 取确定的符号,而过 (0,0) 的每一条直线在 (0,0) 附近都位于 f 取 “+” 的区域内,f(0,0)=0,所以 (0,0) 是 f 在该直线上取到极小值的点。另一方面,(0,0) 在 R2 中的任一小邻域内都包含着上述四个区域,即 f 可以在这个小邻域取到正值,也可以取到负值。故 f 在 (0,0) 处不可能取到极值。
注
根据需要在不同的区域取不同的值,是构造二元函数反例的常用手法。
21.3.3 最大最小值问题
求
f(x,y)=x2−xy+y2−2x+y
在全平面上的最大、最小值。
解 1
令
fx=2x−y−2=0,fy=−x+2y+1=0,
得驻点 (1,0)。又
A=fxx(1,0)=2>0,B=fxy(1,0)=−1,C=fyy(1,0)=2,
且 Δ=AC−B2>0,故 (1,0) 是极小值点,极小值为
f(1,0)=−1.
又有
f(ρcosθ,ρsinθ)=ρ2(1−sinθcosθ)−ρ(2cosθ−sinθ)≥21ρ2−3ρ→+∞(ρ→+∞).
可见 f(x,y) 在全平面上无最大值。又可知存在 ρ0,当 ρ≥ρ0 时,f>−1。于是 f 在 R2 上的最小值与 f 在 D={x2+y2≤ρ02} 上的最小值相同,又 f(x,y) 在 D 内无不可导点,所以
R2minf=Dminf=min{f(1,0),f∣∂D}=−1.
解 2
由
f(x,y)=y2+(1−x)y+(x2−2x)=(y+21−x)2+(43x2−23x−41),
先固定 x,求 miny∈R{f(x,y)},得
y∈Rmin{f(x,y)}=43x2−23x−41=43(x−1)2−1.
从而
x∈Rmin{y∈Rminf(x,y)}=−1,
所以
R2minf(x,y)=−1.
又 f(0,y)=y2+y 无最大值,因此 f(x,y) 也无最大值。
解 3
将函数改写为
f(x,y)=(x−2y−1)2+43y2−1≥−1,
且 f(1,0)=−1,于是 f(x,y) 在 (1,0) 点取得最小值 −1,无最大值。
注
在上例的解 1 中,能否像一元函数那样,根据 f(x,y) 在 R2 上只有唯一的极值点,就断言该极值点是最值点呢?回答是否定的,例如可考虑 f(x,y)=x3−4x2+2xy−y2 在 R2 上的情况。
设 D 为有界闭区域,u(x,y) 在 D 上连续,存在偏导数,且
∂x∂u+∂y∂u=u,u∣∂D=0.
则 u 在 D 上恒为零。
证
用反证法。假设 u 在 D 上有正的最大值或负的最小值。由条件 u∣∂D=0 知这样的最大值或最小值只能在 D 的内部达到。不妨设 (x0,y0)∈intD,且满足
u(x0,y0)>0,u(x0,y0)≥u(x,y),∀(x,y)∈D.
于是 (x0,y0) 也是极大值点,从而
∂x∂u(x0,y0)=∂y∂u(x0,y0)=0.
另一方面,由已知条件
∂x∂u(x0,y0)+∂y∂u(x0,y0)=u(x0,y0),
此与上式矛盾,从而 u≡0。
在平面上给定不在同一直线上的三点 Mi(ai,bi),i=1,2,3。求平面内的一点,使它到三点的距离之和为最小。
解
任取点 M(x,y),令
ρi=(x−ai)2+(y−bi)2,i=1,2,3.
于是所给问题为研究函数
u(x,y)=i=1∑3ρi=i=1∑3(x−ai)2+(y−bi)2
的最小值。除了在三个给定点以外,它处处存在着偏导数
∂x∂u=i=1∑3ρix−ai=i=1∑3cosθi,∂y∂u=i=1∑3ρiy−bi=i=1∑3sinθi,
其中 θi 表示 x 轴正向与以 Mi 为起点的射线 MiM 的夹角。首先找驻点 M0,令两个偏导数为 0 得
cosθ1+cosθ2+cosθ3=0,sinθ1+sinθ2+sinθ3=0.
第一式乘以 sinθ2,第二式乘以 cosθ2,相减得
sin(θ2−θ1)=sin(θ3−θ2).
同样可得
sin(θ3−θ2)=sin(θ1−θ3).
由图 21.1 可知
∠M1M0M2=θ2−θ1,∠M2M0M3=θ3−θ2,∠M3M0M1=θ1−θ3+2π.
由上述等式得
∠M1M0M2=∠M2M0M3=∠M3M0M1=32π.
于是点 M0 可由下列方法求得:往三角形 M1M2M3 的三边上各作一含圆周角为 32π 的弧,三弧的公共点为 M0,称为 Fermat 点。
若三角形没有大于或等于 32π 的内角,则此弧确能在三角形之内相交而确定 M0,这时各边显然都对着顶点在 M0 的等于 32π 的角。在这种情形,就必须比较 u(x,y) 在 M0,M1,M2,M3 这四个点的值。事实上,由余弦定理以及
∠M1M0M2=32π,
有
(M1M2)2=(M0M2)2+(M0M1)2+M0M2⋅M0M1(M0M2+21M0M1)2.
于是
M1M2>M0M2+21M0M1.
同理,
M1M3>M0M3+21M0M1.
两式相加得
M1M2+M1M3>M0M1+M0M2+M0M3,
即 u(M1)>u(M0)。显然此处的 M1 可以换成 M2 或 M3。
若三角形 M1M2M3 有一个内角大于或等于 32π,则三条圆弧没有公共点,驻点也就不存在了。而函数 u(x,y) 在 M1,M2,M3 中之一处,也就是在钝角的顶点处达到其最小值。
注
这一问题再一次说明,在探求函数的最大最小值时,除了驻点以外,导数不存在的点也必须考虑在内。
设通过观测或实验得到一列数据 (xi,yi),i=1,2,⋯,n。它们大体上满足线性关系,即大体上可以用直线方程来反映变量 x 与变量 y 之间的对应关系。确定直线方程使这些数据代入后的偏差的平方和最小。
解
不妨设 xi (i=1,2,⋯,n) 不全相等,即直线非垂直。设方程为
y=ax+b.
偏差的平方和为
f(a,b)=i=1∑n(axi+b−yi)2.
现在要确定 a,b,使得 f(a,b) 为最小。为此,令
fa=2i=1∑nxi(axi+b−yi)=0,fb=2i=1∑n(axi+b−yi)=0.
按 a,b 合并整理得
ai=1∑nxi2+bi=1∑nxi=i=1∑nxiyi,ai=1∑nxi+bn=i=1∑nyi.
解此方程组得唯一的驻点 (a,b)=(aˉ,bˉ)。为进一步确定该点是否取到极小值,计算得到
A=faa=2i=1∑nxi2,B=fab=2i=1∑nxi,C=fbb=2n,
D=AC−B2=4ni=1∑nxi2−4(i=1∑nxi)2>0.
由前述定理知 f(a,b) 在点 (aˉ,bˉ) 取到唯一极小值。并且由 xi 不全相等的假定可知对任意 (a,b) 有
1≤i≤nmax∣axi+b∣>0
在 a2+b2=1 上有正下界,故当 ∣a∣+∣b∣→+∞ 时 f(a,b)→+∞。因此 f 的最小值在有界闭区域上取到,所得的唯一极小值就是最小值。
注 1
在具体问题中,主要指应用问题,如果问题确有最值,而边界值明显不是最值,在区域内部驻点又唯一,则此驻点必是最值点。
注 2
如例题 21.3.7 的注所指出的,多元函数的唯一极值点并不一定是最值点。但如果 f 是二次函数,则 f 的极值点一定是最值点。设 p0 是 f 的极值点,则 p0 是 f 的驻点,在 p0 处作 Taylor 展开,有
f(x)=f(p0)+21(x−p0)TQ(x−p0).
由于 f 是二次函数,故 Q 是常值矩阵。并且 f 在 p0 附近取到极值,因此 Q 是半正定或半负定矩阵。如果 Q 是半正定矩阵,则 f(x)−f(p0)≥0, ∀x,所以 f(p0) 是最小值;如果 Q 是半负定矩阵,则 f(x)−f(p0)≤0, ∀x,所以 f(p0) 是最大值。
21.3.4 练习题
设 u(x,y) 在 (x0,y0) 的邻域上具有连续的二阶偏导数,证明:当 h→0 时有
∂x2∂2u(x0,y0)+∂y2∂2u(x0,y0)=h21{u(x0+h,y0)+u(x0−h,y0)+u(x0,y0+h)+u(x0,y0−h)−4u(x0,y0)}+o(1).
解
由一元 Taylor 公式,
u(x0±h,y0)=u(x0,y0)±hux(x0,y0)+2h2uxx(x0,y0)+o(h2),
u(x0,y0±h)=u(x0,y0)±huy(x0,y0)+2h2uyy(x0,y0)+o(h2).
四式相加后一次项消去,得
u(x0+h,y0)+u(x0−h,y0)+u(x0,y0+h)+u(x0,y0−h)−4u(x0,y0)=h2(uxx(x0,y0)+uyy(x0,y0))+o(h2),
两边除以 h2 即得所证公式。
对于下列函数,点 (0,0) 是否为驻点?是否为极值点?
(1)(2)(3)f(x,y)=x2−4xy+5y2−1;f(x,y)=x2+y2;f(x,y)=(x+y)2−y2.
解
- ∇f=(2x−4y,−4x+10y),故 (0,0) 是驻点。Hesse 矩阵为
H=(2−4−410),
满足 2>0,detH=4>0,故 (0,0) 为极小值点,极小值为 −1。
-
f(x,y)=x2+y2 在原点不可微,故不是驻点;但 f≥0,且 f(0,0)=0,故 (0,0) 是极小值点。
-
f(x,y)=x2+2xy,∇f=(2x+2y,2x),故 (0,0) 是驻点。Hesse 矩阵
H=(2220)
的行列式为 −4<0,故 (0,0) 为鞍点,不是极值点。
求下列函数的极值点及相应的极值:
(1)(2)(3)u(x,y)=x2(y−1)2;u(x,y)=3x2y−x4−2y2;u(x,y)=(1+ey)cosx−yey.
解
-
u=x2(y−1)2≥0,且当 x=0 或 y=1 时取零。因此全部极小值点恰为两条直线 x=0 与 y=1,极小值为 0;无极大值。
-
ux=2x(3y−2x2),uy=3x2−4y。联立得唯一驻点 (0,0)。但沿 y=0 有 u=−x4<0,沿 y=43x2 有 u=81x4>0,故 (0,0) 为鞍点,没有极值。
-
ux=−(1+ey)sinx,uy=ey(cosx−y−1)。驻点满足 sinx=0 且 y=cosx−1,即
(x,y)=(2kπ,0),((2k+1)π,−2),k∈Z.
在这些点上 uxy=0,且 uyy=−ey<0。再有
uxx=−(1+ey)cosx.
因而在 (2kπ,0) 处,uxx=−2<0,故为极大值点,极大值为 2;在 ((2k+1)π,−2) 处,uxx=1+e−2>0,故为鞍点。函数无极小值。
求 f(x,y)=sinxsinysin(x+y) 在
D={(x,y)∣x≥0, y≥0, x+y≤π}
上的最大、最小值。
解
在区域 D 内三因子 sinx,siny,sin(x+y) 都非负,故最小值为 0,在边界 x=0、y=0 或 x+y=π 上取到。
对内点,令
F(x,y)=lnsinx+lnsiny+lnsin(x+y).
则极值点满足
cotx+cot(x+y)=0,coty+cot(x+y)=0.
两式相减得 cotx=coty,在 D 内故 x=y。代回得 cotx+cot2x=0,解得 x=y=3π。于是
fmax=sin3πsin3πsin32π=833.
证明:函数 f(x,y)=x3−4x2+2xy−y2 在 R2 上有唯一的极大值点,但该极大值点不是最大值点。
解
∇f=(3x2−8x+2y,2x−2y),驻点满足 y=x 及 3x2−6x=0,故有 (0,0) 与 (2,2) 两点。Hesse 矩阵为
H=(6x−822−2).
在 (0,0) 处,detH=12>0 且 fxx=−8<0,故为极大值点;在 (2,2) 处,detH=−12<0,故为鞍点。因此极大值点唯一。
但沿直线 y=x,有
f(x,x)=x3−3x2→+∞(x→+∞),
所以该极大值点不是最大值点。
21.4.1 条件极值
求 f(x,y)=ax2+2bxy+cy2 (b=0) 在条件 x2+y2=1 之下的极值。
解 1
作 Lagrange 函数
L(x,y,λ)=ax2+2bxy+cy2−λ(x2+y2−1).
解方程组
Lx=2ax+2by−2λx=0,Ly=2bx+2cy−2λy=0,x2+y2=1.
由 x2+y2=1 知 x,y 不可能同时为零,于是由前两式得
a−λbbc−λ=λ2−(a+c)λ+ac−b2=0.
因判别式
Δ=(a+c)2−4(ac−b2)=(a−c)2+4b2>0,
从而解得两个相异实根
λ±=21[a+c±(a−c)2+4b2].
对每个 λ,由方程与 x2+y2=1 解出两组解。又
LxxLxyLxyLyy=2(a−λ)2b2b2(c−λ)=0,
从而在驻点,矩阵
(LxxLxyLxyLyy)
既不是正定的,也不是负定的。此时,对约束式 x2+y2=1 求微分得
2xdx+2ydy=0,dy=−yxdx.
代入
d2L=2(a−λ)dx2+4bdxdy+2(c−λ)dy2,
得
d2L=y22{(a−λ)y2−2bxy+(c−λ)x2}dx2=y22(a+c−2λ)dx2{<0,0,λ=λ+,λ=λ−.
因此,对应于 λ+ 的点 f 取得极大值,对应于 λ− 的点 f 取得极小值。极大值、极小值分别为 λ+、λ−。
解 2
用隐函数方法。容易证明 b=0 时,f 的极值不在 (±1,0) 取得,而除去这两点外,x2+y2−1=0 在单位圆上任意一点的某邻域内有连续可导的隐函数 y=h(x),且
y′=−yx.
记 F(x)=f(x,h(x)),则
F′(x)=2ax+2b(y+xy′)+2cyy′=y2{(a−c)xy+b(y2−x2)}.
令 F′(x)=0,记 u=yx,由上式得
bu2−(a−c)u−b=0.
它的判别式 (a−c)2+4b2>0,故方程有相异实根
u±=2b1[a−c±(a−c)2+4b2],
且异号。将它们代入 x2+y2=1 可解得对应于 u+ 的解 x1,x2 和对应于 u− 的解 x3,x4。在这些驻点上
F′′(x)=y22{(a−c)(y2−x2)−4bxy}=−b2u{(a−c)2+4b2}{<0,0,u=u+,u=u−.
故 x1,x2 是 f 的极大值点,x3,x4 是 f 的极小值点。其极大值和极小值分别为
u+2+1au+2+2bu++c,u−2+1au−2+2bu−+c.
解 3
注意到
f(x,y)=(x,y)(abbc)(yx)
是一个实的二元二次型。利用高等代数中的有关定理:任意一个实二次型
i,j=1∑naijxixj,aij=aji,
都可以经过正交的线性变换变成平方和
λ1y12+λ2y22+⋯+λnyn2,
其中平方项的系数 λ1,λ2,⋯,λn 就是矩阵 (aij)n×n 的特征多项式全部的根。因此 f(x,y) 可化为
f(x,y)=λ1u2+λ2v2,u2+v2=1,
其中 u,v 是 x,y 的线性组合,然后再用 Lagrange 乘子法或初等方法求解。
注
有些极值问题会有初等解法,但这些初等解法一般都有一定的局限性。
21.4.2 条件最值
椭球面
3x2+y2+2z2=1
被通过原点的平面 2x+y+z=0 截成一个椭圆 l。求此椭圆的面积。
解
只要求出椭圆 l 的长短半轴即可。于是问题转化为在约束条件
F(x,y,z)=3x2+y2+2z2−1=0,G(x,y,z)=2x+y+z=0
下求 r2=x2+y2+z2 的最大、最小值。为此定义
L(x,y,z,λ,μ)=x2+y2+z2+λ(3x2+y2+2z2−1)+μ(2x+y+z).
令
Lx=2x+32λx+2μ=0,Ly=2y+2λy+μ=0,Lz=2z+λz+μ=0.
将这三个式子分别乘以 x,y,z,然后相加,利用约束条件得
λ=−r2.
将 λ 代入 Lx=0 得
(r2−3)x=3μ.
若 r2=3,则 μ=0,代入其余方程得 y=z=0,显然不满足两个约束条件。因此 r2−3=0,所以
x=r2−33μ.
同理
y=2(r2−1)μ,z=r2−2μ.
将它们代入 G=0,得
r2−36μ+2(r2−1)μ+r2−2μ=0.
由于 μ=0,消去 μ,得
15(r2)2−49r2+36=0.
由此可以解出两个根 r12 和 r22,即为条件驻点对应的函数值。由于约束集合是有界闭集,故 r2 的最值存在;又因为在约束集合上
∂(y,z)∂(F,G),∂(z,x)∂(F,G),∂(x,y)∂(F,G)
不同时为 0,故最值必定在条件驻点达到。因而这两个值也即为长、短半轴的平方。根据 Vieta 定理有
r12r22=1536=512.
于是所求的面积为
s=πr1r2=2π53.
注 1
本题所用解法也是条件极值或最值问题中常用的技巧,其特点是不必求出条件驻点,而直接求出所要求的条件最值。
注 2
在得到 λ=−r2 之后可用下面的解法来求 λ。由约束式可知 x,y,z 不能同时为 0,于是把方程组视 x,y,z,μ 为未知量,λ 为参数,系数矩阵的行列式为 0,即
2+32λ00202+2λ01002+λ12110=0.
展开整理得
15λ2−49λ+36=0.
设 p≥1。求
f(x,y)=21(xp+yp)
在 x+y=C (C 为正常数),x≥0, y≥0 的条件下的最小值。
解
设
L(x,y,λ)=21(xp+yp)+λ(x+y−C).
令
Lx=2pxp−1+λ=0,Ly=2pyp−1+λ=0,x+y=C.
解得驻点 x=y=C/2。由于约束集是有界闭集,比较驻点值与附加约束集合边界点的值
f(2C,2C)=(2C)p,f(0,C)=f(C,0)=21Cp,
可知 (2C)p 为 f 在条件 x,y≥0 下的最小值。
注
由本题结论可得
21(xp+yp)≥(2x+y)p.
一般地,可以证明:若 p≥1, xi≥0, i=1,2,⋯,n,则有
n1i=1∑nxip≥(n1i=1∑nxi)p.
很多不等式都可以类似地化为条件最值问题来证明。
求
f(x,y,z)=lnx+2lny+3lnz
的最大值,其中 x2+y2+z2=6r2 (r>0),x>0, y>0, z>0。且证明对任意正数 a,b,c,有
ab2c3≤108(6a+b+c)6.
分析
原问题等价于求 u=xy2z3 的最大值。这样做是因为对数函数在零点无定义,而幂函数则不然,因而可将约束条件化为 x≥0, y≥0, z≥0。注意到在边界上 u 取最小值,故在
D1={(x,y,z)∣x2+y2+z2=6r2, x≥0, y≥0, z≥0}
的内部一定有 u 的最大值点。
解
构造 Lagrange 函数
L(x,y,z,λ)=xy2z3+λ(x2+y2+z2−6r2).
令
Lx=y2z3+2λx=0,Ly=2xyz3+2λy=0,
Lz=3xy2z2+2λz=0,x2+y2+z2=6r2.
往前三式两边分别乘以 x,y,z,然后相加,利用约束条件得
xy2z3=−2λr2.
再将上式代入前三式中,明显地 λ=0,则
x=r,y=2r,z=3r,λ=−33r4.
对应了唯一驻点,于是 xy2z3 的最大值为 63r6,即
xy2z3≤63r6.
两边平方得
x2y4z6≤108r12.
令 a=x2, b=y2, c=z2,则 r2=6a+b+c。于是
ab2c3≤108(6a+b+c)6.
注
也可以用算术平均值—几何平均值不等式证明:
ab2c3=2233a(2b)2(3c)3≤108(6a+b+c)6.
21.4.3 隐函数的极值
求由方程
2x2+y2+z2+2xy−2x−2y−4z+4=0
所确定的函数 z=z(x,y) 的极值。
解 1(直接法)
设
g(x,y,z)=2x2+y2+z2+2xy−2x−2y−4z+4.
由于
∂z∂g=2z−4,
由 ∂z∂g=0 得 z=2。在满足方程的 (x,y,z),z=2 附近确定了隐函数 z=z(x,y)。由原方程对隐函数 z=z(x,y) 求导,得
4x+2z∂x∂z+2y−2−4∂x∂z=0,
2y+2z∂y∂z+2x−2−4∂y∂z=0.
令 ∂x∂z=∂y∂z=0,则
2x+y−1=0,y+x−1=0.
由此得到驻点为 (0,1),将之代入原方程得
z2−4z+3=0.
于是 z1=1, z2=3。
在上述两式两边求导,得
4+2(∂x∂z)2+2z∂x2∂2z−4∂x2∂2z=0,
2∂x∂z∂y∂z+2z∂x∂y∂2z+2−4∂x∂y∂2z=0,
2+2(∂y∂z)2+2z∂y2∂2z−4∂y2∂2z=0.
在 (x,y,z)=(0,1,1),∂x∂z=∂y∂z=0,得
A=zxx=2,B=zxy=1,C=zyy=1,
于是 D=AC−B2>0,故 z=1 为极小值。
在 (x,y,z)=(0,1,3),∂x∂z=∂y∂z=0,得
A=zxx=−2,B=zxy=−1,C=zyy=−1,
于是 D=AC−B2>0,由此可见 z=3 为极大值。
解 2
以原方程为约束条件,取目标函数 f(x,y,z)=z,则 Lagrange 函数为
L(x,y,z,λ)=z+λ(2x2+y2+z2+2xy−2x−2y−4z+4).
令
Lx=4λx+2λy−2λ=0,Ly=2λx+2λy−2λ=0,
Lz=1+2λz−4λ=0.
显然 λ=0,于是由前两式得驻点为 (0,1)。代入原方程得 z1=1, z2=3。再由 Lz=0 得
λ1=−21,λ2=21.
由于
Lxx=4λ,Lyy=Lzz=2λ,Lxy=2λ,Lxz=Lyz=0,
于是 L 在 (0,1,1) 与 (0,1,3) 的 Hesse 矩阵分别为
210110001,−2−10−1−1000−1.
前者正定,后者负定,所以 z=1 为极小值,z=3 为极大值。
注
在上述条件极值问题中约束条件给出了一个有界闭集,所以条件最值一定取到。而且可以证明条件最值一定是条件极值,所以 z=1 和 z=3 分别是该方程确定的隐函数 z=z(x,y) 的最小值和最大值。
21.4.4 练习题
求 f(x,y,z)=x4+y4+z4 在条件 xyz=1 下的极值。
解
由 AM—GM 不等式,
x4+y4+z4≥3(x4y4z4)1/3=3.
当且仅当 ∣x∣=∣y∣=∣z∣=1 时取等;再由 xyz=1 得极小值为 3,在 (1,1,1) 及恰有两个分量为 −1 的三点处取到。令 x=t,y=t,z=t−2,则函数值趋于 +∞,故无极大值。
求 u=x−2y+2z 在条件 x2+y2+z2=1 下的极值。
解
在线性函数 u=(1,−2,2)⋅(x,y,z) 在单位球面上的最值为
±12+(−2)2+22=±3.
最大值在 (31,−32,32) 处取得,最小值在其相反点处取得。
求 z=cos2x+cos2y 满足条件 x−y=4π 时的极值。
解
由约束 y=x−4π,得
z=cos2x+cos2(x−4π)=1+2cos2x+sin2x.
故
zmax=1+22,zmin=1−22.
当 x=8π+kπ, y=−8π+kπ 时取最大值;当 x=85π+kπ, y=83π+kπ 时取最小值。
求 f(x,y,z)=xlymzn 满足条件 ax+by+cz=k 的极值,常数和变量均正。
解
令 S=l+m+n,考虑 lnf=llnx+mlny+nlnz。Lagrange 方程为
xl=λa,ym=λb,zn=λc.
故
x=λal,y=λbm,z=λcn.
代入约束 ax+by+cz=k,得 λ=kS。于是极大值点为
x=aSkl,y=bSkm,z=cSkn,
极大值为
fmax=(Sk)Salbmcnllmmnn.
当某一变量趋于 0+ 时 f→0,故无最小值,只以下确界 0。
求
u=a2x2+b2y2+c2z2
在条件
x2+y2+z2=1,xcosα+ycosβ+zcosγ=0
下的极值,其中 a>b>c>0,cos2α+cos2β+cos2γ=1。
解
记
q1=a21,q2=b21,q3=c21,
以及 n=(cosα,cosβ,cosγ)。在约束 ∣x∣=1、x⋅n=0 下求二次型 u=xTQx 的最值,其中 Q=diag(q1,q2,q3)。Lagrange 方程给出
(Q−λI)x=μn,
从而
xi=qi−λμni.
再由 x⋅n=0 得
q1−λcos2α+q2−λcos2β+q3−λcos2γ=0.
即
λ2−Aλ+B=0,
其中
A=q1+q2+q3−q1cos2α−q2cos2β−q3cos2γ,
B=q2q3cos2α+q1q3cos2β+q1q2cos2γ.
此二次方程的两根就是所求最小值与最大值。
设约束条件为
x11+x21+⋯+xn1=a1,xi>0,i=1,2,⋯,n,
而 a>0 为定值。求 u=x1x2⋯xn 的极值。
解
设 g=∑i=1nxi−1−a−1=0。取对数后作 Lagrange 乘子法:
xi1−λxi21=0(i=1,…,n),
故 x1=⋯=xn=λ。代入约束得 λ=na,于是
x1=⋯=xn=na,
且
umin=(na)n.
另一方面让某个 xi→∞ 而其余变量调整保持约束,可使 u→∞,故无最大值。
求满足约束条件
x+y+z=2π,x>0, y>0, z>0
时,u=sinxsinysinz 的极值。
解
令 F=lnu=lnsinx+lnsiny+lnsinz,约束为 x+y+z=2π。Lagrange 方程给出
cotx=coty=cotz,
由 x,y,z∈(0,π/2) 得 x=y=z=6π。故
umax=sin36π=81.
当某个变量趋于 0+ 时 u→0,故无最小值,只有下确界 0。
若
a1x1+a2x2+⋯+anxn=1,
求 u=x12+x22+⋯+xn2 的极值,其中 ai>0 (i=1,2,⋯,n) 是常数,变量 x1,x2,⋯,xn 均大于 0。
解
由 Cauchy 不等式,
1=(i=1∑naixi)2≤(i=1∑nxi2)(i=1∑nai21).
因而
u≥∑i=1nai−21.
当且仅当 xi=λ/ai 时取等。由约束求得
xi=∑j=1naj−2ai−1,
所以最小值为
umin=∑i=1nai−21.
若要求 xi>0,则最大值不存在;其上确界为 maxiai2。
求 u=x3+y3+z3−2xyz 在 x2+y2+z2≤1 中的最大、最小值。
解
由于 u 是三次齐次函数,闭球上的最值必在边界球面上取得。由 Lagrange 方程
3x2−2yz=2λx,3y2−2zx=2λy,3z2−2xy=2λz,
相减可得
(x−y)(3(x+y)+2z−2λ)=0,
以及循环式。于是驻点只可能是三类:x=y=z=±31;某一坐标为 0、其余两坐标中一个为 ±1;或两坐标相等而第三个与之成固定比。分别代入可得函数值只可能是
±1,±331,±54196.
因而最大值为 1,最小值为 −1,分别在 (±1,0,0)、(0,±1,0)、(0,0,±1) 中符号同号与异号的点上取到。
求 z=x2−xy+y2 在 ∣x∣+∣y∣≤1 中的最大、最小值。
解
有
x2−xy+y2≥0,
故最小值为 0,在 (0,0) 处取得。又因
x2−xy+y2≤x2+∣xy∣+y2≤(∣x∣+∣y∣)2≤1,
所以最大值不超过 1,且在 (±1,0),(0,±1) 处取到。故最大值为 1。
在球面 x2+y2+z2=1 上求一点,使其到 n 个已知点
Mi(xi,yi,zi)(i=1,2,⋯,n)
距离的平方和最小。
解
设 P=(x,y,z),S=∑i=1n(xi,yi,zi)。则
i=1∑n∣P−Mi∣2=n+i=1∑n∣Mi∣2−2P⋅S,
因为 ∣P∣=1。故该和最小时 P⋅S 最大。若 S=0,最优点为
P=∣S∣S.
若 S=0,则上式与 P 无关,球面上任意点都可。
求点 (−1,0) 到半立方抛物线 y2=x3 的最短距离。
解
取参数方程 x=t2, y=t3。到点 (−1,0) 的距离平方为
d2(t)=(t2+1)2+t6=t6+t4+2t2+1.
其导数为 2t(3t4+2t2+2),只在 t=0 处为零,故最小距离为 1,最近点是原点。
解
设三边为 a,b,c,半周长为 s。由 Heron 公式
Δ=s(s−a)(s−b)(s−c).
当周长固定时 s 固定,于是由 AM—GM 不等式,(s−a)(s−b)(s−c) 在三者相等时最大,即 a=b=c。故面积最大的三角形是正三角形。
解
设 x1x2⋯xn=a,则
x11+⋯+xn1≥n(x1x2⋯xn1)1/n=a1/nn.
当且仅当 x1=⋯=xn=a1/n 时取等。故倒数和最小时各因子应相等。
长为 a 的铁丝截成两段,一段围成一个正方形,另一段围成一个圆。这两段各为多少时,正方形与圆面积之和达最大值。
解
设长度 x 围成正方形,则余长 a−x 围成圆,总面积为
A(x)=16x2+4π(a−x)2.
这是开口向上的二次函数,因此在闭区间 [0,a] 上最大值取于端点。由于 4π1>161,故最大值在 x=0 处取得,即全部铁丝围成圆。若要求两段都严格为正,则最大值不存在,其上确界仍为 4πa2。
过椭圆 3x2+2xy+3y2=1 上任意点作此椭圆的切线,求切线与两坐标轴围成的三角形面积最小者。
解
在椭圆上一点 (x,y) 处的切线为
(6x+2y)X+(2x+6y)Y=2.
它与两轴围成的三角形面积为
S=216x+2y22x+6y2=2∣(3x+y)(x+3y)∣1.
令 u=x+y, v=x−y,则约束变为 2u2+4v2=1,而
(3x+y)(x+3y)=3x2+10xy+3y2=4u2−21v2.
在约束下其绝对值最大为 4,当且仅当 v=0,即 x=y=±221。于是
Smin=81.
求原点到曲面 (x−y)2+z2=1 的最短距离。
解
令 u=x−y, v=x+y,则约束为 u2+z2=1。而
x2+y2=2u2+v2≥2u2,
等号当且仅当 v=0,即 y=−x。故距离平方满足
r2=x2+y2+z2≥2u2+z2=1−2u2≥21.
所以最短距离为 21,在 (±21,∓21,0) 处取得。
求 x2−xy+y2−z2=1 与 x2+y2=1 的交线到原点的最短距离。
解
由 x2+y2=1 与 x2−xy+y2−z2=1 得
z2=−xy.
因而 xy≤0,且到原点距离平方为
r2=x2+y2+z2=1−xy≥1.
当且仅当 xy=0 时取等,于是最短距离为 1,在 (±1,0,0) 与 (0,±1,0) 处取得。
解
设所给圆半径为 r,三角形内角为 A,B,C。则
S=r2(cot2A+cot2B+cot2C),2A+2B+2C=2π.
由 Jensen 不等式(或 AM—GM)知,上式在 A=B=C=3π 时最小。故面积最小的外切三角形是正三角形。
解
设四边长为 a,b,c,d,面积为 S,则由 Bretschneider 公式
S2=(s−a)(s−b)(s−c)(s−d)−abcdcos22B+D,
其中 s=2a+b+c+d。四边长固定时第一项固定,故当且仅当 cos2B+D=0,即 B+D=π 时面积最大。于是面积最大的四边形必为圆内接四边形。
在一个已知凸四边形内求一点 C,使其到四个顶点距离的平方和最小。
解
设四顶点位置向量分别为 A1,A2,A3,A4,点 C 的位置向量为 x。则
i=1∑4∣x−Ai∣2=4x−4A1+A2+A3+A42+常数.
因而最小值在
x=4A1+A2+A3+A4
处取得。这个点正是两条对角线中点连线的中点,也即 Varignon 平行四边形的中心,因此位于凸四边形内部。
证明:椭圆的内接三角形中,面积达最大的三角形的一个顶点处的法线,必定与此三角形的该顶点所对的边正交,并求面积最大的内接三角形。
解
设椭圆由线性变换 T(u,v)=(au,bv) 将单位圆 u2+v2=1 映成。单位圆内接三角形面积最大时必为正三角形,其面积为 433;经 T 后面积放大 ab 倍,故椭圆内接三角形的最大面积为
433ab.
极值三角形正是单位圆内接正三角形的像,可写成
(acost,bsint),(acos(t+32π),bsin(t+32π)),
(acos(t+34π),bsin(t+34π)).
在单位圆上,正三角形一顶点处的法线就是半径,并与对边垂直。该对边经 T 后的法向量变为 (au0,bv0),正是椭圆在点 (au0,bv0) 处的法向量,所以椭圆上极值三角形同样满足“一个顶点处的法线与其对边垂直”。
在椭圆 x2+4y2=4 上求一点,使到直线 3x+4y=12 的距离最短。
解
到直线 3x+4y=12 的距离为
d=512−3x−4y,
因为椭圆 x2+4y2=4 完全位于该直线的同侧。故只需求 3x+4y 在椭圆上的最大值。令 x=2cost, y=sint,则
3x+4y=6cost+4sint≤213,
等号在 cost=133、sint=132 时取到。故最近点为
(136,132),
最短距离为
dmin=512−213.
设椭圆
a2x2+b2y2=1
的任一切线与 x,y 轴分别交于 A,B 两点,试求线段 AB 长度的最小值。
解
设切点为 (x0,y0),则切线方程为
a2xx0+b2yy0=1.
与坐标轴交于 A(x0a2,0)、B(0,y0b2),从而
∣AB∣2=x02a4+y02b4.
再由 a2x02+b2y02=1,令 x0=acost, y0=bsint,得
∣AB∣2=cos2ta2+sin2tb2≥(a+b)2.
故 ∣AB∣min=a+b。等号条件为 cos2t=a+ba、sin2t=a+bb。
在抛物线 y=x2 的所有与法线重合的弦中求长度最短的弦。
解
设法线经过抛物线上点 (t,t2),则法线方程为
y−t2=−2tx−t(t=0).
与抛物线再交于 x=s,代入 y=x2 得
s=−t−2t1.
弦长平方为
L2=(s−t)2+(s2−t2)2.
利用 s+t=−2t1 化简得
L2=16t4(4t2+1)3.
令 u=4t2>0,则 L2=u2(u+1)3。求导得极小值在 u=2 处取得,即 t2=21。故最短弦长为
Lmin=233.
设椭圆
a2x2+b2y2=1
将圆 (x−1)2+y2=1 包围于其内。试求 a,b 的值,使符合上述条件的椭圆的面积为最小。
解
椭圆包住圆 (x−1)2+y2=1,等价于对一切 t 都有
a2(1+cost)2+b2sin2t≤1.
令 u=cos22t∈[0,1],并记 p=a24, q=b24,则上式化为
(p−q)u2+qu≤1(0≤u≤1).
为使面积 πab=pq4π 最小,只需在上述约束下使 pq 最大。最优情形必有 q>p,此时二次式在顶点处取最大值,故需
4(q−p)q2≤1,
即 p≤q−4q2。因此
pq≤q2−4q3.
右端在 q=38 时取最大值,且此时 p=98。于是
a2=p4=29,b2=q4=23.
因而面积最小的椭圆为
9/2x2+3/2y2=1.
设方程 F(x,y)=0 满足隐函数定理的条件,并由此确定了隐函数 y=f(x),又设 F(x,y) 具有连续的二阶偏导数。
- 求 f′′(x);
- 若 F(x0,y0)=0,y0=f(x0) 为 f(x) 的一个极值,证明:当 Fy(x0,y0) 与 Fxx(x0,y0) 同号时,f(x0) 为极大值;当 Fy(x0,y0) 与 Fxx(x0,y0) 异号时,f(x0) 为极小值;
- 对方程 x2+xy+y2=27 在隐函数形式下(不解出 y),求 y=f(x) 的极值,并用 (2) 的结论判别极大或极小。
解
- 由 F(x,f(x))≡0 得
Fx+Fyf′=0,
再求导一次,得
Fxx+2Fxyf′+Fyy(f′)2+Fyf′′=0.
因而
f′′=−FyFxx+2Fxyf′+Fyy(f′)2=−Fy3FxxFy2−2FxFxyFy+Fx2Fyy.
- 若 x0 是极值点,则 f′(x0)=0,故 Fx(x0,y0)=0,且
f′′(x0)=−Fy(x0,y0)Fxx(x0,y0).
因此当 Fy 与 Fxx 同号时,f′′(x0)<0,故为极大值;异号时 f′′(x0)>0,故为极小值。
- 取 F(x,y)=x2+xy+y2−27。极值点满足 F=0 且 Fx=2x+y=0。解得
(x,y)=(3,−6),(−3,6).
由 Fy=x+2y、Fxx=2,在 (3,−6) 处有 Fy=−9,与 Fxx 异号,故 y=−6 为极小值;在 (−3,6) 处有 Fy=9,与 Fxx 同号,故 y=6 为极大值。
设 D 是由两抛物线 y=x2−1, y=−x2+1 所围成的闭域,试在 D 内求一椭圆
a2x2+b2y2=1,
使其面积为最大。
解
记上边界为 y=1−x2。若椭圆内接于 D,由对称性只需考虑中心在原点的椭圆
a2x2+b2y2=1.
它位于 D 内当且仅当
b1−a2x2≤1−x2(∣x∣≤a≤1).
令 s=a2、t=x2,则
b2≤s−ts(1−t)2=:ψs(t),0≤t≤s.
求导得
ψs′(t)=(s−t)2s(t−1)(2s−t−1).
当 s≤21 时,minψs=s,故 ab≤s;当 s≥21 时,最小值在 t=2s−1 处取得,且
b2≤4s(1−s),ab≤2s1−s.
右端在 s=32 时最大,于是
a2=32,b2=98.
故最大面积椭圆为
2/3x2+8/9y2=1,
最大面积为
πab=334π.
21.6.2 第一组参考题
曲面 z=z(x,y) 由方程
F(z−cx−a,z−cy−b)=0
确定,其中 a,b,c 为常数。证明:
- 曲面的切平面经过定点 (a,b,c);
- 函数 z=z(x,y) 满足方程
∂x2∂2z⋅∂y2∂2z−(∂x∂y∂2z)2=0.
解
设 P=(x,y,z) 在曲面上,而 λ∈R。则
c+λ(z−c)−ca+λ(x−a)−a=z−cx−a,c+λ(z−c)−cb+λ(y−b)−b=z−cy−b,
故整条直线 (a,b,c)+λ(x−a,y−b,z−c) 都在曲面上。于是该曲面是以 (a,b,c) 为顶点的锥面,故每一点的切平面都经过定点 (a,b,c)。
若局部写成 z=z(x,y),则由上面的齐次性得
z(a+t(x−a),b+t(y−b))−c=t(z(x,y)−c).
对 t 在 t=1 处求导,得
(x−a)zx+(y−b)zy=z−c.
再对 x,y 求导,得
(x−a)zxx+(y−b)zxy=0,(x−a)zxy+(y−b)zyy=0.
消去 x−a,y−b 即得
zxxzyy−zxy2=0.
在平面上给定一个边长为 a,b,c 的三角形,在它上面可作无数个等体积的给定高 h 的三棱锥,在其中求出有最小侧面积 S 的那一个。
解
设三棱锥顶点在底面内的正射影为 P,到三边的距离分别为 da,db,dc。其侧面积为
S=21(ah2+da2+bh2+db2+ch2+dc2).
而底面积固定,满足 ada+bdb+cdc=2Δ。作 Lagrange 乘子法可得
h2+da2da=h2+db2db=h2+dc2dc,
从而 da=db=dc。故 P 必为底面三角形的内心,此时侧面积最小。
设 u=f(x,y,z) 为二次可微函数,记 cosα,cosβ,cosγ 为方向 l 的方向余弦。求
∂l2∂2u=∂l∂(∂l∂u).
解
设 l=(cosα,cosβ,cosγ),则
∂l∂u=uxcosα+uycosβ+uzcosγ.
再沿 l 求导,得
∂l2∂2u=uxxcos2α+uyycos2β+uzzcos2γ+2uxycosαcosβ+2uyzcosβcosγ+2uzxcosγcosα.
设 u=f(x,y,z) 为二次可微函数,li=(cosαi,cosβi,cosγi),i=1,2,3 为三个互相垂直的单位向量,证明:
1.
(∂l1∂u)2+(∂l2∂u)2+(∂l3∂u)2=(∂x∂u)2+(∂y∂u)2+(∂z∂u)2;
∂l12∂2u+∂l22∂2u+∂l32∂2u=∂x2∂2u+∂y2∂2u+∂z2∂2u;
- 若 ∂li∂f=0, i=1,2,3,则 f(x,y,z) 在 R3 中恒为常数。
解
设三个单位向量作为正交矩阵 A 的三行。则 A 为正交阵。
- 因为
(∂l1∂u,∂l2∂u,∂l3∂u)T=A(ux,uy,uz)T,
正交变换保持 Euclid 范数,故
i=1∑3(∂li∂u)2=ux2+uy2+uz2.
- 将上一题中的公式对 i=1,2,3 求和,并利用
i=1∑3cos2αi=i=1∑3cos2βi=i=1∑3cos2γi=1,
以及交叉项和为 0,便得
i=1∑3∂li2∂2u=uxx+uyy+uzz.
- 若 ∂li∂f=0 对 i=1,2,3 都成立,则由 (1) 得 ∣∇f∣2=0,故 ∇f≡0。由于 R3 连通,f 恒为常数。
设 f 在 Oδ(x0)⊂Rn 上连续,在 Oδ(x0)∖{x0} 可微,证明:
- 如果 ∀x∈Oδ(x0)∖{x0},有 (x−x0)⋅∇f(x)<0,则 x0 是 f 的一个极大值点;
- 如果 ∀x∈Oδ(x0)∖{x0},有 (x−x0)⋅∇f(x)>0,则 x0 是 f 的一个极小值点。
解
对任意 x=x0,令 ϕ(t)=f(x0+t(x−x0)),0≤t≤1。则
ϕ′(t)=(x−x0)⋅∇f(x0+t(x−x0)).
由题设可知,当 t>0 时有 ϕ′(t)<0,故 ϕ 在 (0,1] 上严格递减,从而
f(x)=ϕ(1)<ϕ(0)=f(x0).
这表明 x0 是严格极大值点。极小值情形同理。
设 f 在 Oδ0(x0)⊂Rn 上二次连续可微,Hesse 矩阵
(∂xi∂xj∂2f(x0))
为正定阵,∇f(x0)=0。证明:∃δ∈(0,δ0),使
(x−x0)⋅∇f(x)>0,∀x∈Oδ(x0)∖{x0}.
解
由 Taylor 公式,
∇f(x)=∇f(x0)+H(x0)(x−x0)+o(∣x−x0∣)=H(x0)(x−x0)+o(∣x−x0∣).
因而
(x−x0)⋅∇f(x)=(x−x0)TH(x0)(x−x0)+o(∣x−x0∣2).
由正定性,存在 m>0 使 hTH(x0)h≥m∣h∣2。取 δ 足够小,使余项绝对值不超过 2m∣x−x0∣2,便有
(x−x0)⋅∇f(x)≥2m∣x−x0∣2>0(0<∣x−x0∣<δ).
用条件极值的方法证明 H”older 不等式:
i=1∑naixi≤(i=1∑naip)1/p(i=1∑nxiq)1/q,
其中 ai≥0, xi≥0, i=1,2,⋯,n;p,q>1,p1+q1=1。
解
先在约束 ∑xiq=1 下最大化 ∑aixi。Lagrange 方程为
ai=λqxiq−1,
故 xi=caip−1。由约束求得
c=(i=1∑naip)−1/q,
从而最大值为 (∑aip)1/p。将 xi 按任意比例伸缩即得 H”older 不等式。
证明:二次型
f(x,y,z)=Ax2+By2+Cz2+2Dyz+2Ezx+2Fxy
在单位球面 x2+y2+z2=1 上的最大值与最小值恰好是矩阵
AFEFBDEDC
的最大特征值与最小特征值。
解
设矩阵 M 为题中对称矩阵,二次型可写为 f(x)=xTMx。在约束 ∣x∣=1 下,Lagrange 方程为
Mx=λx,
即驻点正是 M 的单位特征向量,而相应函数值为特征值。故球面上的最大值与最小值分别等于 M 的最大、最小特征值。
设 A=(aij) 是 n 阶方阵,且设
j=1∑naij2=Hi,i=1,2,⋯,n,
其中 H1,H2,⋯,Hn 是 n 个确定的非负实数。
- 证明:如果矩阵 A 的行向量是 Rn 中两两正交的向量,则 (detA)2 取得极值;
- 根据等式 (detA)2=detA⋅detAT,其中 AT 是矩阵 A 的转置矩阵,证明:
max{(detA)2}=H1×H2×⋯×Hn;
- 证明:对任意的矩阵 B=(bij) 有 Hadamard 不等式
(detB)2≤j=1∏n(i=1∑nbij2);
- 给 Hadamard 不等式以直观的几何解释。
解
记 A 的行向量为 r1,…,rn。则
(detA)2=det(AAT),
其中 G=AAT 的对角元恰是 H1,…,Hn。
- 若 ri 两两正交,则 G=diag(H1,…,Hn),故
(detA)2=H1H2⋯Hn.
这给出极值候选。
- 一般地,由 Hadamard 不等式对正定矩阵 G 有
detG≤i=1∏nGii=H1H2⋯Hn,
且等号当且仅当 G 对角,即行向量两两正交。因此
max(detA)2=H1H2⋯Hn.
- 对任意矩阵 B,把它的列向量看成行向量应用上式,即得
(detB)2≤j=1∏ni=1∑nbij2.
- 几何上,这就是说由列向量张成的平行多面体体积的平方,不超过各边长平方之积;等号当且仅当这些边互相正交。
设 u(x,y) 在 x2+y2≤1 上连续,在 x2+y2<1 上满足
∂x2∂2u+∂y2∂2u=u.
- 若在 x2+y2=1 上 u(x,y)≥0,证明:当 x2+y2≤1 时,u(x,y)≥0;
- 若在 x2+y2=1 上 u(x,y)>0,证明:当 x2+y2≤1 时,u(x,y)>0。
解
设 u 在闭圆盘上取得最小值。若边界上 u≥0 而某内点 p 上有 u(p)<0,则 p 是内极小点,故 uxx(p)+uyy(p)≥0;但方程给出 uxx(p)+uyy(p)=u(p)<0,矛盾。故结论 (1) 成立。
若边界上 u>0 而内部有点使 u=0,则 u 在内部取得最小值 0。由强极值原理可知 u 应恒为常数 0,这与边界严格为正矛盾,所以 (2) 也成立。
解
设 p0 为驻点,H 为 Hesse 矩阵。若 H 正定,则存在 m>0 使 hTHh≥m∣h∣2。Taylor 公式给出
f(x0+h)−f(x0)=21hTHh+o(∣h∣2)≥4m∣h∣2
对充分小的 h=0 成立,故 x0 为严格极小值点。负定情形同理给出严格极大值。若 H 不定,则存在 h1,h2 使二次项一正一负,再由 Taylor 展开可知 x0 不是极值点。
仿照 21.4.2 小节对具有多个约束条件的条件最值的确定给出详细的讨论。
解
设约束为 g1(x)=⋯=gm(x)=0,并假定在候选点 x0 处雅可比矩阵秩为 m。必要条件是存在 λ1,…,λm 使
∇f(x0)=j=1∑mλj∇gj(x0).
令 L=f−∑λjgj。若在切空间
Tx0M={h∈Rn∣∇gj(x0)⋅h=0, j=1,…,m}
上,二次型 d2L(x0)[h,h] 正定,则 x0 为条件极小值点;负定则为条件极大值点;不定则不是条件极值。这就是多重约束情形的判别原则。
21.6.2 第二组参考题
设二次连续可微函数 f 在点 p0 处的梯度为零向量,但 Hesse 矩阵 Q 为非零矩阵。证明:f 在 p0 处的函数值增长最快的方向位于 Q 的最大特征值对应的特征子空间中。如果最大特征值对应的特征子空间为一维空间,则相应的方向就是 f 增长最快的方向。
解
设 l 为单位向量,则
f(p0+tl)−f(p0)=2t2lTQl+o(t2).
当 t>0 很小时,增长最快的方向正是使 Rayleigh 商 lTQl 最大的方向。由谱定理,该最大值为 Q 的最大特征值,且达到最大值的单位向量恰是对应特征子空间中的单位向量。若该特征子空间一维,则方向唯一。
设 R3 中的光滑曲线段 L 的参数表达式为
x=x(t),y=y(t),z=z(t),t∈[a,b].
P(x(a),y(a),z(a)),Q(x(b),y(b),z(b)) 为 L 的起点与终点。Π 是通过直线段 PQ 的任一平面。证明:∃ξ∈(a,b),使得向量 (x′(ξ),y′(ξ),z′(ξ)) 与 Π 平行。
解
取平面 Π 的法向量 n。因为线段 PQ 平行于 Π,故 n⋅(Q−P)=0。令
ϕ(t)=n⋅(x(t),y(t),z(t)).
则 ϕ(a)=ϕ(b)。由 Rolle 定理,存在 ξ∈(a,b) 使 ϕ′(ξ)=0,即
n⋅(x′(ξ),y′(ξ),z′(ξ))=0.
这说明切向量与 Π 平行。
(Sanderson(桑德森)中值定理)设 f:[a,b]→Rm 为 k 次可微向量值函数,并且非零向量 v 与 f(a),f(b) 及 f(j)(a),j=1,2,⋯,k−1 正交。证明:∃ξ∈(a,b),使 v 与 f(k)(ξ) 正交。其中 f(j) 表示 f 的各个分量函数的 j 阶导数组成的向量。
解
令 ϕ(t)=v⋅f(t)。题设给出
ϕ(a)=ϕ(b)=0,ϕ(j)(a)=0(j=1,2,…,k−1).
于是对 ϕ 连续应用 k 次 Rolle 定理,可得某个 ξ∈(a,b) 使
ϕ(k)(ξ)=v⋅f(k)(ξ)=0.
证明:与曲面 ax2+by2+cz2=1 (abc=0) 相切的三个互相垂直的平面的交点在球面
x2+y2+z2=a1+b1+c1
上。
解
设三个切平面的法向量为 ni=(axi,byi,czi),其中切点 (xi,yi,zi) 满足 axi2+byi2+czi2=1。由切平面互相垂直知 ni 两两正交。设三平面的交点为 X,则
ni⋅X=1(i=1,2,3).
于是把 X 在正交基 mi=ni/∣ni∣ 下展开,得
∣X∣2=i=1∑3∣ni∣21.
又由
1=niTdiag(a1,b1,c1)ni=∣ni∣2miTdiag(a1,b1,c1)mi
得
∣ni∣21=miTdiag(a1,b1,c1)mi.
求和并利用 {mi} 正交归一,便有
∣X∣2=trdiag(a1,b1,c1)=a1+b1+c1.
设曲面 Σ:z=f(x,y),(x,y)∈D 的每一点的法线与 z 轴相交,其中 D 是 R2 中的圆环或圆盘,f 是 C1 函数。证明:Σ 是一个旋转曲面。
解
曲面法线在点 (x,y,f(x,y)) 处的方向为 (−fx,−fy,1)。法线与 z 轴相交,等价于存在 t 使
x−tfx=0,y−tfy=0.
消去 t 得
yfx=xfy.
而极坐标下 fθ=−yfx+xfy,故 fθ≡0。于是 f(x,y)=φ(x2+y2),曲面是旋转曲面。
证明:不等式
2ab≤ea−1+alna+eb−1+blnb
对所有 a>0, b>0 成立,并求出等式成立的充分必要条件。
解
先证基本不等式
uv≤eu−1+vlnv(u>0,v>0).
事实上令 t=eu−1/v,则
eu−1+vlnv−uv=v(t−lnt−1)≥0.
将它分别应用于 (u,v)=(a,b) 与 (b,a),相加即得
2ab≤ea−1+alna+eb−1+blnb.
等号成立当且仅当同时有 b=ea−1 与 a=eb−1。由 lnx≤x−1 知其唯一解为 a=b=1。
设 T:u=φ(x,y), v=ψ(x,y) 是平面上的一个变换。如果在该变换下任何两条相交曲线的交角保持不变,则称该变换是保角的。
- 如果 φ 与 ψ 满足 φx=ψy, φy=−ψx,则 T 是保角变换;
- 证明:下面的反演变换是保角的
ξ=x2+y2x,η=x2+y2y;
- 证明:任一圆的反演像是另一圆或直线;
- 求出反演变换的 Jacobi 行列式;
- 证明:连续可微变换 T 是保角的充分必要条件是它满足 Cauchy—Riemann 方程
φx=ψy,φy=−ψx
或
φx=−ψy,φy=ψx.
第一种情况保持角度方向,第二种情况角度反向。
解
设 J=(φxψxφyψy)。
- 若满足 Cauchy—Riemann 方程 φx=ψy, φy=−ψx,则
J=(ab−ba),
因而 JTJ=(a2+b2)I。故任意切向量经 J 作用后长度同乘以同一因子,夹角不变,所以 T 保角。
- 对反演变换
ξ=x2+y2x,η=x2+y2y,
直接计算得
ξx=−ηy=(x2+y2)2y2−x2,ξy=ηx=−(x2+y2)22xy,
它满足反向的 Cauchy—Riemann 方程,故也是保角变换。
- 任意圆或直线都可写成
α(x2+y2)+βx+γy+δ=0.
代入反演公式 x=ξ2+η2ξ、y=ξ2+η2η,化简为
δ(ξ2+η2)+βξ+γη+α=0.
若 δ=0,像是圆;若 δ=0,像是直线。故任一圆的反演像是另一圆或直线。
- Jacobi 行列式为
∂(x,y)∂(ξ,η)=−(x2+y2)21.
- 连续可微变换保角,当且仅当其微分在每一点都是“伸缩 × 正交变换”。对 2×2 矩阵,这等价于
JTJ=λI,
即两列向量等长且正交。解此代数条件恰得两类矩阵
(ab−ba),(abb−a),
这正对应于两组 Cauchy—Riemann 方程。前者保持角度方向,后者反向。
由公式
ξ=x2+y2+z2x,η=x2+y2+z2y,ζ=x2+y2+z2z
定义三维反演变换。
- 证明:任何两个曲面之间的夹角在反演变换下不变;
- 证明:球面被变换为球面或平面;
- 求变换的 Jacobi 行列式。
解
三维反演写作 T(x)=∣x∣2x。其微分为
dTx=∣x∣21(I−2Px),
其中 Px 是沿径向的正交投影。因此 dTx 是“按 ∣x∣−2 缩放,再关于过原点的切平面反射”,故保持任意两向量的夹角,于是任意两曲面的夹角也保持不变。
任意球面或平面都可写成
α(x2+y2+z2)+βx+γy+δz+ε=0.
代入反演公式后化为
ε(ξ2+η2+ζ2)+βξ+γη+δζ+α=0.
故像仍是球面或平面。
由 dTx 的特征值为 ∣x∣−2,∣x∣−2,−∣x∣−2,可得 Jacobi 行列式
∂(x,y,z)∂(ξ,η,ζ)=−(x2+y2+z2)31.
评论
支持 Markdown 和 LaTeX 数学公式。