3.3 习题课讲解

3.3.1 高阶偏导数

例 3.3.1 (例1) 已知函数f(u,v)具有2阶连续偏导数,且函数g(x,y)=f(2x+y,3x−y)满足∂2g∂x2+∂2g∂x∂y−6∂2g∂y2=1。

(1)

求∂2f∂u∂v;

(2)

若∂f(u,0)∂u=ue−u,f(0,v)=150v2−1,求f(u,v)的表达式。

解 (1) 由链式法则得到(1){g1=2f1+3f2g2=f1−f2 从而(2){g11=4f11+12f12+9f22g12=2f11+f12−3f22g22=f11−2f12+f22 因此(3)1=g11+g12−6g22=25f12⟹∂2f∂u∂v=125

(2) 参考图3.3.1,我们有(4)∂f∂u(u,v)=∂f∂u(u,0)+∫0v∂2f∂u∂v(u,t)dt=ue−u+v25 进而(5)f(u,v)=f(0,v)+∫0u∂f∂u(t,v)dt=150v2−1+uv25+(−1−u)e−u ◻

PIC

图 3.3.1: 例1(2)计算图示

另解(1)1 由变量替换{u=2x+yv=3x−y以及Jacobi矩阵形式的链式法则,得到偏微分算子的关系(6)(∂∂x∂∂y)=(∂∂u∂∂v)(uxuyvxvy)=(∂∂u∂∂v)(213−1) 故有(7)∂2∂x2+∂2∂x∂y−6∂2∂y2=(∂∂x+3∂∂y)(∂∂x−2∂∂y)=(13)(∂∂x∂∂y)(∂∂x∂∂y)(1−2)=(13)(231−1)(∂∂u∂∂v)(∂∂u∂∂v)(231−1)(1−2)=(50)(∂2∂u2∂2∂u∂v∂2∂u∂v∂2∂v2)(05)=25∂2∂u∂v 因此(8)∂2f∂u∂v=125 ◻

另解(1)2 由另解(1)1第一步解得(9)(∂∂u∂∂v)=(∂∂x∂∂y)(213−1)−1=(∂∂x∂∂y)(151535−25) 从而(10)∂2∂x2+∂2∂x∂y−6∂2∂y2=(∂∂x+3∂∂y)(∂∂x−2∂∂y)=(5∂∂u)(5∂∂v)=25∂2∂u∂v ◻

例 3.3.2 (例2,star) 设A是m阶正交矩阵,函数f:Rm→R∈C2、F(x)=f(Ax)。令y=Ax,证明:(11)∑k=1m(∂F∂xk)2=∑k=1m(∂f∂yk)2,∑k=1m∂2F∂xk2=∑k=1m∂2f∂yk2

证明 (1) 可直接通过全微分算子的协变变换证明:(12)∂∂xk∂∂xk=∂∂yiAik∂∂yjAjk=∂∂yi(AAT)ij∂∂yj=∂∂yiδij∂∂yj=∂∂yi∂∂yi 或者利用梯度算子的逆变特性。坐标向量(逆变)的变换矩阵为A,故基向量(协变)的变换矩阵B=A−1,因此(13)∇y=B−1∇x=A∇x⟹∇xF(x)=AT∇yf(y)⟹∥∇xF(x)∥2=∥∇yf(y)∥2

(2) 可借助Hesse矩阵3证明4:(14)HF(x)=∂x(∇xF)(x)=A−1∂y(∇yf)(y)A=A−1Hf(y)A 故(15)∑k=1m∂2F∂xk2=trHF(x)=trHf(y)=∑k=1m∂2f∂yk2 ◻

例 3.3.3 (例3,star) 设x=x(u,v),y=y(u,v)满足(16)∂x∂u=∂y∂v,∂x∂v=−∂y∂u 证明:

(1)

若w满足∂2w∂x2+∂2w∂y2=0,则w~(u,v):=w(x(u,v),y(u,v))满足∂2w~∂u2+∂2w~∂v2=0。

(2)

∂2(xy)∂u2+∂2(xy)∂v2=0。

证明 可微函数z(u,v)=x(u,v)+iy(u,v)在单连通域C上满足Cauchy-Riemann方程,从而解析。

(1) w在单连通域R2上调和,故存在解析函数f使得w(x,y)=ℜf(x+iy),因此f∘z亦为解析函数,从而w~(u,v)=ℜf(z(u+iv))调和。

(2) 由(1)知w~=xy调和。 ◻

例 3.3.4 (star) 设f:R2→R调和,令u(x,y)=f(xx2+y2,yx2+y2)。证明:u调和。

证明 令ξ(x.y)=xx2+y2、η(x,y)=−yx2+y2,其满足CR方程,由上例可知u~=f(ξ,η)调和。由于(ξ,η)↦(ξ,−η)是正交变换,故u调和。 ◻

注 有关以上两题的更多信息,大家可以参考 王兆臻学长的习题课笔记5。

3.3.2 Taylor展开

例 3.3.5 (例4,star) 求cos⁡xcos⁡y在(0,0)的带Peano余项的4阶Taylor公式和带Lagrange余项的1阶Taylor公式。

解 f(x,y)=cos⁡xcos⁡y在y=0的邻域中是C∞函数。(17)f(x,y)=1−x22+x424+o(x4)1−y22+y424+o(y4)=(1−x22+x424+o(x4))(1+y22+5y424+o(x4))=1+y2−x22+x4−6x2y2+5y424+o(r4),r=x2+y2→0 由Peano余项Taylor公式的唯一性,这就是在(0,0)处的Peano余项4阶Taylor公式。

计算f(x,y)的Hesse矩阵(18)fxx=−cos⁡xcos⁡y,fxy=−sin⁡xsin⁡ycos2⁡yfyy=(cos⁡xsin⁡ycos2⁡y)y′=cos⁡xcos3⁡y+2cos⁡xcos⁡ysin2⁡ycos4⁡y 所以(19)f(x,y)=1+fxx(θx,θy)2x2+fxy(θx,θy)xy+fyy(θx,θy)2y2=1−cos⁡(θx)2cos⁡(θy)x2−sin⁡(θx)sin⁡(θy)cos2⁡(θy)xy+cos⁡(θx)cos⁡(θy)(1+sin2⁡(θy))2cos4⁡(θy)y2 其中0<θ<1。 ◻

注 以下是带r阶Lagrange余项的r−1阶Taylor公式:(20)f(x0+v)=f(x0)+∑k=1r−11k!∑i1,⋯,ik=1n∂kf∂xik⋯∂xi1(x0)vi1⋯vik+1r!∑i1,⋯,ir=1n∂rf∂xir⋯∂xi1(x0+θv)vi1⋯vir 其中θ∈(0,1)。具体来说,设f:R2→R∈C2,对于带2阶Lagrange余项的1阶Taylor公式,我们有(21)f(x,y)=f(x0,y0)+∂f∂x(x0,y0)(x−x0)+∂f∂y(x0,y0)(y−y0)+12[∂2f∂x2(x1,y1)(x−x0)2+2∂2f∂x∂y(x1,y1)(x−x0)(y−y0)+∂2f∂y2(x1,y1)(y−y0)2] 其中x1=x0+θ(x−x0),y1=y0+θ(y−y0),θ∈(0,1)。所以本题需要计算到函数的所有二阶偏导数,以得到带Lagrange余项的1阶Taylor公式。

例 3.3.6 (例5) 求(1−x2−y2)−1/2在(0,0)的带Peano余项和带Lagrange余项的Taylor公式。

解 由广义二项式定理可得(22)(1−r2)−1/2=1−12(−r2)+(−12)(−12−1)2(−r2)2+⋯+(2n−1)!!2nr2n+o(r2n) 故有(23)(1−x2−y2)−1/2=1+∑k=1n(2k−1)!!2k(x2+y2)k+o(r2n),r=x2+y2→0 注意到(24)((1+x)−1/2)(n)=(−1)n(2n−1)!!2n(1+x)−12−n 故有(25)(1+x)−1/2=1−12x+1⋅323x2+⋯+(−1)n(2n−1)!!2n(1+θx)−12−nxn 因此(26)(1−x2−y2)−1/2=1+12(x2+y2)+1⋅323(x2+y2)2+⋯ +(2n−1)!!2n(1−θ(x2+y2))−12−n(x2+y2)n ◻

例 3.3.7 (例6,star) 求xy在(1,1)处的带Peano余项的2阶Taylor展开。

解 令x=1+u、y=1+v,则当(u,v)→(0,0)时,有(27)xy=e(1+v)ln⁡(1+u)=e(1+v)(u−u22+o(u2))=eu+uv−u22+o(r2),r=u2+v2=1+(u+uv−u22)+12u2+o(r2)=1+u+uv+o(r2)=1+(x−1)+(x−1)(y−1)+o(r2),r=(x−1)2+(y−1)2 ◻

注 以下是带Peano余项的r阶Taylor公式:(28)f(x0+v)=∑k=0r1k!∑i1,⋯,ik=1n∂kf∂xik⋯∂xi1(x0)vi1⋯vik+o(∥v∥r) 具体来说,设f:R2→R∈C2,对于带2阶Lagrange余项的1阶Taylor公式,我们有(29)f(x,y)=f(x0,y0)+∂f∂x(x0,y0)(x−x0)+∂f∂y(x0,y0)(y−y0)+12[∂2f∂x2(x0,y0)(x−x0)2+2∂2f∂x∂y(x0,y0)(x−x0)(y−y0)+∂2f∂y2(x0,y0)(y−y0)2]+o((x−x0)2+(y−y0)2) 所以本题需要计算到函数的所有二阶偏导数,以得到带Peano余项的2阶Taylor公式。

特别提醒 Taylor公式中的易错点:

(1)

上文中提到的Taylor公式的阶。

(2)

Taylor多项式不能化简,需要保持xi−x0i因子。

(3)

在展开到第n阶时需要除以n!。

(4)

Taylor公式的唯一性是指Taylor多项式的唯一性,余项的形式不一定唯一。如f(x,y)=ex2−y2,利用上式可得其在原点附近的1阶Taylor公式为(30)f(x,y)=1+x2−y2+[(x2−y2)+2θ2(x2−y2)2]eθ2(x2−y2) 如果借助一元Taylor公式,则有(31)f(x,y)=1+x2−y2+12(x2−y2)2eθ(x2−y2)

3.3.3 极值与函数凹凸性

例 3.3.8 (例7) 求z=(x2+2y2)e−x2−y2的极值、最值、值域。

解1 只求一阶导数,用梯度向量场分析函数单调性。由z≥0=z(0,0)知z的最小值为0。求偏导可知(32)zx=2e−x2−y2x[1−(x2+2y2)]zy=2e−x2−y2y[2−(x2+2y2)] 函数可微,极值点必然是驻点,从而zx=zy=0,解得(33){x1=0y1=0,{x2=0y2=±1,{x3=±1y3=0 由梯度向量场分布(图3.3.2)知,(0,0)是最小值点,(0,1)和(0,−1)是最大值点,(1,0)和(−1,0)是鞍点,不是极值点。

现在讨论z的值域。由对称性,只需讨论x≥0,y≥0部分。显然z(0,0)=0是最小值。

1∘ 当x>1时,x2+2y2>1,zx<0,故z(x,y)关于x>0严格减,即z(x,y)<z(1,y)。由(34)zy(1,y)=2e−1−y2y(1−2y2) 知z(1,y)关于y在区间(0,12)上单增,在区间(12,+∞)上单减,从而z(1,y)≤z(1,12)。由(35)zx(x,12)=−2e−x2−12x3<0 知z(x,12)关于x>0是减函数,从而z(x,12)≤z(0,12)。再由(36)zy(0,y)=2e−y2y(2−2y2) 知z(0,y)在y=1取得最大值z(0,1)=2e−1。所以对于任意x≥1以及任意y≥0,都有z(x,y)≤z(0,1)。

2∘ 当0<x≤1时,则zy(x,y)的符号与2−(x2+2y2)相同,故z(x,y)关于y在(0,1−x22)上单增、在(1−x22,+∞)上单减,从而z(x,y)≤z(x,1−x22)=2e−12x2−1≤z(0,1)=2e−1。

由连续函数的介值性质,z的值域为区间[0,2e−1]。 ◻

PIC
(a) z=z(x,y)的等高线图和梯度向量场分布
PIC
(b) 函数值z(x,y)的变化(沿绿色箭头函数值增大)
图 3.3.2: 例7图示

解2 求驻点,计算Hesse矩阵,判断正、负定,分析极值类型。由z≥0=z(0,0)知z的最小值为0。利用(37)lim(x,y)→∞f(x,y)=0<z(1,0) 知,存在R>1使得当x2+y2≥R2时z(x,y)<z(1,0)。于是z(x,y)在有界闭集x2+y2≤R2中的最大值点就是函数在整个定义域上的最大值点,因此是驻点,即函数值最大的驻点就是最大值点。 ◻

解3 利用极坐标(38)0≤z=r2e−r23−cos⁡2θ2≤2e−1 ◻

例 3.3.9 (例10,star) 讨论函数u(x,y)=(1+ey)cos⁡x−yey的极值。

解 计算可知(39)∇u=(−(1+ey)sin⁡xey(cos⁡x−1−y)) 解得驻点(40)(2kπ0),((2k+1)π−2),k∈Z 随后分情况讨论。继续计算可得(41)uxx=−(1+ey)cos⁡x,uxy=−eysin⁡x,uyy=ey(cos⁡x−2−y) 代入第一类驻点可得(42)Hu(2kπ,0)=(−200−1) 负定,所以(2kπ,0)是极大值点,极大值为2。代入第二类驻点可得(43)Hu((2k+1)π,−2)=(1+e−200−e−2) 特征值一正一负,所以((2k+1)π,0)是鞍点。 ◻

PIC

图 3.3.3: z=u(x,y)的三维图像的透视图

例 3.3.10 (例11) 证明对任意x>0、y>0,都有(44)ex−1+xln⁡x+ey−1+yln⁡y≥2xy

证明 记(45)F(x,y)=ex−1+xln⁡x+ey−1+yln⁡y−2xy 计算可得(46)Fx=ex−1+1+ln⁡x−2y,Fy=ey−1+1+ln⁡y−2x 进一步计算可得(47)Fxx=ex−1+1x,Fxy=−2,Fyy=ey−1+1y 其中(48)Fxxx=x2ex−1−1x2,(x2ex−1−1)′=ex−1(x2+2x)>0,∀x>0 所以u(x)=x2ex−1−1关于x是增函数。注意到u(1)=0,故Fxx关于x在x=1时取最小值Fxx(1,y)=2,因此Fxx≥2。

同理Fyy≥2,故HF(x,y)总是半正定,从而F(x,y)是凸函数,在R2上存在最小值。再注意到Fx(1,1)=Fy(1,1)=0,所以F(1,1)=0是最小值。因此对任意x>0,y>0,总有F(x,y)≥0。 ◻

例 3.3.11 (例8) 设正数a1,…,an满足a1a2⋯an=1,求a1+a2+⋯+an的取值范围。

解1 记(49)f(x1,x2,…,xn−1)=x1+x2+⋯+xn−1+1x1x2⋯xn−1,x1,x2,…,xn−1>0 计算可得(50)fk(x)=1−1x1x2⋯xn−1xk,fkj(x)=1+δkjx1x2⋯xn−1xkxj 所以(51)Hf(x)=1x1x2⋯xn−1(2x121x1x2⋯1x1xn−11x1x22x22⋯1x2xn−1⋮⋮⋱⋮1x1xn−11x2xn−1⋯2xn−12)=1x1x2⋯xn−1(1x11x2⋮1xn−1)(1x11x2⋯1xn−1)+diag{1x12,1x22,…,1xn−12}x1x2⋯xn−1 总是半正定,且定义域[0,+∞)n−1是凸集,所以f是凸函数。

由f1=f2=⋯=fn−1=0得到x1=x2=⋯=xn−1=1,所以f(1,1,…,1)=n是最小值。又f(N,1,…,1)>N,所以f的值域为[n,+∞)。 ◻

解2 类似解1构造f,计算可得(52)fn−1=1−1x1x2⋯xn−1xn−1 所以对给定的x1,x2,…,xn−2,f关于xn−1在xn−1=1x1x2⋯xn−2处取得最小值。继续构造(53)g(x1,x2,…,xn−2)=f(x1,x2,…,xn−2,1x1x2⋯xn−2)=x1+x2+⋯+xn−2+2x1x2⋯xn−2 计算可得(54)gn−2=1−1x1x2⋯xn−2xn−22 所以g关于xn−2在xn−2=1x1x2⋯xn−33处取得最小值。如此类推,直到(55)h(x1,x2)=x1+x2+n−2x1x2n−2 由(56)h2=1−1x1x2n−1n−2 知h关于x2在x2=1x1n−1处取得最小值。(57)u(x1)=h(x1,1x1n−1)=x1+n−1x1n−1,u′=1−1x1nn−1 u在x1=1处取得最小值u(1)=n。因此f在x1=1,x2=1,…,xn−1=1时取得最小值n。又f(N,1,…,1)>N,所以f的值域为[n,+∞)。 ◻

3.3.4 与极值有关的证明题(1)

例 3.3.12 (例9) 设u(x,y)在x2+y2≤1上连续,在x2+y2<1内C2,且(58){uxx+uyy=u,x2+y2<1u(x,y)=0,x2+y2=1 证明:当x2+y2<1时,u(x,y)≤0。

证明 反证法。假设存在(x0,y0)满足x02+y02<1、u(x0,y0)>0,因为u(x,y)在x2+y2≤1上连续,所以u(x,y)有最大值u(x1,y1),于是u(x1,y1)≥ u(x0,y0)>0。因为沿边界x2+y2=1,u(x,y)=0,所以x12+y12<1,因此(x1,y1)是极大值,从而Hu(x1,y1)半正定,故有(59)uxx(x1,y1)+uyy(x1,y1)=trHu(x1,y1)≤0 这与(60)uxx(x1,y1)+uyy(x1,y1)=u(x1,y1)>0 矛盾。所以在x2+y2<1内,恒有u(x,y)≤0。 ◻

注 同理可得当x2+y2<1时,u(x,y)≥0,所以u(x,y)=0对所有(x,y)∈D―成立。

例 3.3.13 (例9拓展1,star) 设D⊆R2为有界开集,函数f:D―→R∈C(D―)满足f∈C2(D),且(61){Δf=f,(x,y)∈Df=φ(x,y),(x,y)∈∂D 证明:

(1)

若φ≥0对所有(x,y)∈∂D成立,则f≥0对所有(x,y)∈D成立。

(2)

若φ>0对所有(x,y)∈∂D成立,则f>0对所有(x,y)∈D成立。

证明 (1) 设(x0,y0)为f在D―上的最小值点。假设f(x0,y0)<0,则(x0,y0)∈D;又f∈C2(D),由Fermat引理可知(x0,y0)为驻点,且0>f(x0,y0)=Δf(x0,y0)=trHf(x0,y0)≥0,矛盾!故f(x0,y0)≥0,即f≥0对所有(x,y)∈D成立。

(2) 记m:=min(x,y)∈∂Df(x,y)>0、X:=max(x,y)∈∂Dx<+∞。构造函数(62)g(x,y):=f(x,y)−mex−X,(x,y)∈D 计算可得(63)Δg=Δf−mΔex−X=f−mex−X=g,(x,y)∈Dg=f−mex−X≥f−m≥0,(x,y)∈∂D 由(1)可得g≥0对所有(x,y)∈D成立,因此(64)f(x,y)=g(x,y)+mex−X≥mex−X>0,(x,y)∈D 故f>0对所有(x,y)∈D成立。 ◻

例 3.3.14 (例9拓展2,star) 设区域D⊆Rn,u,f,φ∈C(D―)、u∈C2(D)、c∈R,证明:以下边值问题的解唯一。(65){Δu=cu+f,r∈Du=φ,r∈∂D

证明 设u1,u2均满足该方程和边界条件,令u=u1−u2,则u满足(66){Δu=cu,r∈Du=0,r∈∂D

当c≠0时,令v=uc,则v满足(67){Δv=v,r∈Dv=0,r∈∂D v在有界闭集D―上存在最小值点x0,显然v(x0)<0。若v(x0)<0,则有x0∉∂D,故x0为极小值点,从而Δv(x0)=trHv(x0)≥0>v(x0),矛盾!故v(x)≥0对任意x∈D―成立。同理可得v(x)≤0,故v≡0,即u1≡u2。

当c=0时,注意到(68)∫D∥∇u∥2dV=∫D[∇⋅(u∇u)−∇2u]dV=∮∂Du∇u⋅ndS−∫D∇2udV=0 由u∈C(D―)可得u≡0,即u1≡u2。

综上所述,原边值问题具有唯一解。 ◻

3如前所述,trHf=gik∂i(gkj∂jf),其等于∂i∂i当且仅当G为单位矩阵。不要想当然地认为此式成立。

4坐标变换y↦x是逆变,协变变换矩阵为A;∇是逆变算子,变换矩阵A−1左乘;∂是协变算子,变换矩阵A右乘。

5./figure/laplace_wzz.pdf。