9.3 习题课讲解

9.3.1 Green公式

例 9.3.1 (习题课8·例4) C:|x|+|y|=2为正向闭曲线,用Green公式计算:(1)Caxdybydx|x|+|y|

由Green公式可得(2)I=12C(axdybydx)=a+b2Ddxdy=a+b28=4(a+b)

例 9.3.2 (习题课8·例6) fC1[1,4]f(1)=f(4)γxy直角坐标平面的第一象限中由直线y=xy=4x和曲线xy=1xy=4在所围成的平面有界区域D的正向边界(图8.3.3),用Green公式计算(3)γf(xy)ydy

由Green公式可得(4)Df(xy)ydy=Dxf(xy)ydxdy=Df(xy)dxdy=[1,4]2f(s)|det(x,y)(t,s)|dtds,(t=yx, s=xy)=1414f(s)1tdtds(t,s)(x,y)=|yx21xyx|=yx=t=(f(4)f(1))ln4=0 这个解法需要f连续可微。

例 9.3.3 (习题课8·例7) Dt={(x,y)R2x2+y2t2, t>0}f(x,y)Dt上连续,在Dt内存在连续偏导数。若f(x,y)Dt上满足方程2fx2+2fy2=12f(x,y)n为有向曲线Dt的外单位法向量,求极限limt01πt2Dtfndl

注意到(5)fn=fn,Δf=divf 由Green公式可得(6)Dtfndl=Dtfndl=Dtdivfdxdy=DtΔfdxdy=Dtf(x,y)2dxdy=(f(0,0)2+o(1))|Dt|,t0 所以所求极限为f(0,0)2

例 9.3.4 (习题课8·例8) f(x,y)C2(R)fxx(x,y)+fyy(x,y)=e(x2+y2),证明:(7)x2+y21(xfx+yfy)dxdy=π2e

CrDr分别表示圆心位于原点、半径为r、自然正向的圆周和圆盘,由Green公式可得(8)I=01rdr02πf(rn)dθ=01(Crfndl)rdr=01(DrΔfdxdy)rdr=01(Dreρ2ρdρdθ)rdr=01πeρ2|ρ=0rrdr=π01(1er2)rdr=π2(r2+er2)r=01=π2e

例 9.3.5 (习题课8·例9(1)) f(x)是正值连续函数,D为圆心在原点的单位圆,DD的正向边界,用Green公式证明:(9)D[xf(y)dyyf(x)dx]=D[yf(x)dx+xf(y)dy]

证明 由Green公式可得(10)D[xf(y)dyyf(x)dx]=D[(xf(y))x(yf(x))y]dxdy=D[f(y)+1f(x)]dxdy=D[f(x)+1f(y)]dxdy=D[yf(x)dx+xf(y)dy]

重积分中积分换元(x,y)(y,x)的Jacobi行列式为1,积分换元公式取行列式绝对值、为1,这就是重积分中的对称性。

例 9.3.6 (习题课8·例10) 设在上半平面D={(x,y)y>0}内,函数f:R2R2具有连续偏导数,且对任意的t>0都有f(tx,ty)=t2f(x,y)用Green公式证明:对D内的任意分段光滑的有向简单闭曲线L,都有(11)I=L[yf(x,y)dxxf(x,y)dy]=0

L围成的区域为D,由Green公式可得(12)L[yf(x,y)dxxf(x,y)dy]=D[2f(x,y)+xfx(x,y)+yfy(x,y)]dxdy 注意到(13)f(tx,ty)=t2f(x,y)xf1(tx,ty)+yf2(tx,ty)=2t3f(x,y)t=1,可得(14)xfx+yfy=2f 故有(15)L[yf(x,y)dxxf(x,y)dy]=0

例 9.3.7 Dδ={(x,y)δ2x2+y21},设f:DδR2C1满足x2+y2=1f(x,y)=0,证明:(16)limδ0+Dδxfx(x,y)+yfy(x,y)x2+y2dxdy=2πf(0,0)

证明 由题得(17)I(δ):=Dδxfx(x,y)+yfy(x,y)x2+y2dxdy=Dδfrr2dS=DδflnrdS 注意到2(18)(flnr)=flnr+f2lnr,2lnr=Δlnr=0Cr:={(x,y)x2+y2=r2}(自然正向),则有(19)I(δ)=Dδ(flnr)dS=Dδ+flnrndl=C1frr2rrdlCδfrr2rrdl=Cδfdlr 由积分中值定理可得θδ使得(20)I(δ)=f(δcosθδ,δsinθδ)Cδdlr=2πf(δcosθδ,δsinθδ)δ0+2πf(0,0)

例 9.3.8 u为有界开集DR2上的调和函数(记Δu:=uxx+uyy=0),证明:

(1)

(21)u(x0,y0)=12πD+(ulnrnlnrun)dl 其中(x0,y0)Dr=(x,y)(x0,y0)r=|r|nD的单位外法向量。

(2)

(22)u(x0,y0)=12πRL+u(x,y)dl 其中L为以(x0,y0)为圆心、R为半径的位于D中的任意圆周。

证明 (1) 注意到(23)ulnrnlnrun=(ulnrlnru)n 原积分在(x0,y0)处有瑕点。记Cδ:={(x,y)(xx0)2+(yy0)2=δ2}(自然正向),Dδ=D{(x,y)(xx0)2+(yy0)2δ2},由Green公式可得(24)RHS=12πD+Cδ+Cδ(ulnrlnru)ndl,(D+Cδ=Dδ+)=12πDδ(ulnrlnru)dS+12πCδ(ulnrlnru)ndl=12πDδ(u2lnrlnr2u)dS+12πCδ(urr2lnru)rrdl 注意到2lnr=2u=0,计算可得(25)RHS=0+12πCδudlr12πCδunlnrdl 由积分中值定理可得θδ使得(26)Cδudlr=u(x0+δcosθδ,y0+δsinθδ)Cδdlrδ0+2πu(x0,y0) 由积分的三角不等式可得(27)|Cδunlnrdl|Cδ|un||lnr|dl2πδ|lnδ|maxCδuδ0+0 因此(28)RHS=u(x0,y0)=LHS

(2) 由(1)的结论可得(29)u(x0,y0)=12πL+(ulnrnlnrun)dl=12πL+urr2rrdllnR2πL+undl=12πRL+u(x,y)dllnR2πDR2udS=12πRL+u(x,y)dl

例 9.3.9 φ:R2R2是微分同胚,(x,y)T=φ(u,v)D0为有界闭区域,其边界D0分段C1D1:=φ(D0)。用Green公式证明:(30)|D1|=D0|det(x,y)(u,v)|dudv

证明 首先,我们需要证明以下引理:

引理 9.3.10 φD0映射为D1,则映射前后曲线正向是否改变取决于det(x,y)(u,v)的符号。

引理证明 我们首先需要给“曲线正向”一个数学表述。

对于平面闭合曲线L,其前向单位切向量为τ^、外向单位法向量n^z轴的(正向)单位向量为k^,则k^(n^×τ^)的符号反映了L的定向。设τ,n分别为与τ^,n^同向的切向量、法向量,则k=n1τ2n2τ1的符号反映了L的定向。

D0f(u,v)0确定(否则考虑f),则n=f为外法向量(指向f增大的方向)。设D0+可以参数化为t(uv),则τ=(u(t)v(t))为前切向量。

现考虑变换φ:D0D1(u,v)(x,y),令f~=fφ1,则D1f~(x.y)0确定,n=(x,y)f~D1+可以参数化为t(uv)φ(xy),则τ=(x(t)y(t))。记J=Jφ=(x,y)(u,v),由链式法则可得(31)(xy)=(Jφ)(uv),(x,y)f~=(x,y)(fφ1)=(Jφ1)T(u,v)f 因此τ=Jτn=JTn

k=n1τ2n2τ1,等式两端分别同乘τ1,τ2,结合nτ=n1τ1+n2τ2=0可得(32){kτ1=n1τ1τ2n2τ12=n2(τ12+τ22)kτ2=n1τ22n2τ2τ1=n1(τ12+τ22){n1=kτ2τ12+τ22=k~τ2n2=kτ1τ12+τ22=k~τ1(33)J=(abcd)JT=1detJ(dcba) 则有(34)k=n1τ2n2τ1=1detJ[(dn1cn2)(cτ1+dτ2)(bn1+an2)(aτ1+bτ2)]=k~detJ[(dτ2+cτ1)(cτ1+dτ2)(bτ2aτ1)(aτ1+bτ2)]=k~detJ(τ12+τ22)=kdetJτ12+τ22τ12+τ22k的符号与kdetJ相同。

引理证明结束,我们回到原命题。注意到(35)I1=D1xdy=abx(t)y(t)dt=abx(t)[yuu(t)+yvv(t)]dt=D0xyudu+xyvdv=GreenD0[u(xuv)v(xuu)]dudv=D0(xuyvxvyu)dudv=D0det(x,y)(u,v)dudvdetJ>0,则映射前后曲线正向不变,故有(36)|D1|=I1=D0|det(x,y)(u,v)|dudvdetJ<0,则映射前后曲线正向改变,故有(37)|D1|=I1=D0|det(x,y)(u,v)|dudv 综上所述,原命题得证。

9.3.2 恰当方程与积分因子

例 9.3.11 (习题课9·例1) 求以下微分方程的通解:

(1)

(x2y)dx(x+sin2y)dy=0

(2)

eydx+(xey2y)dy=0

(3)

xdx+ydyx2+y2=ydxxdyx2

(4)

(cosx+1y)dx+(1yxy2)dy=0

(1) 注意到(38)x2dx+cos2y12ydy(ydx+xdy)=d(x33+sin2y4y2xy) 因此通解为x33+sin2y4y2xy=C。如不便凑出全微分,则可先验证无旋条件(39)d[(x2y)dx(x+sin2y)dy]=dydxdxdy=0 微分形式在R2上定义,R2是单连通集合,所以有原函数。计算可得(40) (0,0)(x,y)[(x2y)dx(x+sin2y)dy]=(0,0)(x,0)x2dx(x,0)(x,y)(x+sin2y)dy=x33xy+sin2y2y2

(2) 注意到(41)(eydx+xeydy)2ydy=d(xeyy2)=0 因此通解为xeyy2=C

(3) 注意到(42)d(x2+y2)=d(yx) 因此通解为x2+y2+yx=C

(4) 注意到(43)cosxdx+dyy+ydxxdyy2=d(sinx+lny+xy)=0 因此通解为sinx+lny+xy=C

例 9.3.12 (习题课9·例2) 求以下微分方程的通解:

(1)

(ycosxxsinx)dx+(ysinx+xcosx)dy=0

(2)

(x2sin2y)dx+xsin2ydy=0

(3)

(x+y)dx+(yx)dy=0

(4)

(x+y)(dxdy)=dx+dy

(5)

(3x2+y)dx+(2x2yx)dy=0

(1) 原方程不恰当,故需要引入积分因子μ,其满足(44)(μ(ycosxxsinx))y=(μ(ysinx+xcosx))x 亦即(45)μy(ycosxxsinx)=μ(ycosxxsinx)+μx(ysinx+xcosx)μx=0,即可由μy=μ解得μ=ey,此时有(46)ey(ycosxxsinx)dx+ey(ysinx+xcosx)dy=d(eyxcosx+yeysinxeysinx)=0 因此通解为eyxcosx+yeysinxeysinx=C

(2) 原方程不恰当,故需要引入积分因子μ,其满足(47)(μ(x2sin2y))y=(μxsin2y)x 亦即(48)μy(x2sin2y)=μxxsin2y+2μsin2yμy=0,即可由xμx+2μ=0解得μ=1x2,此时有(49)dx+xsin2ydysin2ydxx2=d(x+sin2yx) 因此通解为x2+sin2y=Cx

(3) 这不是恰当方程。注意到(50)xdx+ydyx2+y2+ydxxdyx2+y2=12d(x2+y2)+darctanxy=0 故通解为(51)12ln(x2+y2)+arctanxy=C 另一种做法是利用极坐标换元,计算可得原式为(52)ω=rdrr2dθ=0drrdθ=0 故通解为(53)lnrθ=C 写成直角坐标形式即为(54)x=Ceθcosθ,y=Ceθsinθ

(4) 这不是恰当方程。注意到(55)d(xy)=d(x+y)x+y 故通解为(56)xy=ln|x+y|+C 若引入积分因子μ,其需满足(57)(μ(x+y1))y=(μ(xy1))x 亦即(58)2μ+(x+y+1)μx+(x+y1)μy=0 考虑微分方程组的一组线性无关的特解(59){x(t)=x(t)+y(t)+1y(t)=x(t)+y(t)1ddt(xy)=2x(t)y(t)=2t 从而(60)x=x+y+1=2x2t+1(e2tx)=e2t(12t)=(e2tt) 因此(61)x=Ce2t+t,y=Ce2tt,C0 其中C0是为了保证x,y线性无关,不妨选C=1,代回可得(62)ddtμ(e2t+t,e2tt)=μxx(t)+μyy(t)=2μμ(e2t+t,e2tt)=e2t 于是选择μ(x,y)=eyx为积分因子,故有(63)eyx(x+y)=C

(5) 原方程不恰当,故需要引入积分因子μ,其满足(64)(μ(3x2+y))y=(μ(2x2yx))x 亦即(65)2(12xy)μ+x(12xy)μx+(3x2+y)μy=0μy=0,即可由xμx+2μ=0解得μ=1x2,此时有(66)3dx+2ydy+ydxxdyx2=d(3x+y2yx) 因此通解为3x2+xy2y=Cx

9.3.3 第二型曲面积分

例 9.3.13 (例3) Σ为锥面z=x2+y2被柱面x2+y2=2x所截的有限部分,规定曲面S的正向向下,所得的定向曲面记为Σ+。求下面两个积分的值:

(1)

ΣzdS

(2)

Σ+x2+y2+z2(xdydz+ydzdx+zdxdy)

(1) Σ在柱坐标系(r,θ,z)下的参数方程为(67)x=rcosθ,y=rsinθ,z=r 其中r22rcosθ,即(68)0r2cosθ,π2θπ2 计算可得(69)xr=(cosθ,sinθ,1)T,xθ=(rsinθ,rcosθ,0)T 从而(70)E=xrxr=2,F=xrxθ=0,G=xθxθ=r2 故有(71)dS=EGF2drdθ=2rdrdθ 因此(72)ΣzdS=π/2π/2dθ02cosθr2rdr=2π/2π/2r33|02cosθdθ=16230π/2cos3θdθ=3229

(2) 直接利用楔积计算可得(73) Σ+x2+y2+z2(xdydz+ydzdx+zdxdy)=0r2cosθ2r2(rcosθd(rsinθ)dr+rsinθdrd(rcosθ)+rd(rcosθ)d(rsinθ))=0r2cosθ2r2(r2cos2θdθdrr2sin2θdrdθ+r2drdθ)=0 或者利用其几何意义:在笛卡尔坐标系下,dydz对应于yOz平面中的投影面积,其系数对应于向量场的x坐标,所以xdydz+ydzdx+zdxdy的系数构成向量场(x,y,z)T。向量场V=(x,y,z)T在点(x,y,z)处与锥面Σ相切,从而与法向量正交。所以所求积分(V的通量)为零。

请注意:Σ锥面而不是柱面的一部分,故dS不为rdrdθ

例 9.3.14 (例4) Σ+为圆柱面Σ:x2+y2=1位于0z2的部分(图9.3.1),外法向为正,计算曲面积分:(74)I=Σ+[x(yz)dydz+(xy)dxdy]

PIC

图 9.3.1: 圆柱面Σ

1 记向量场V=(x(yz),0,xy)。根据题设,Σ+的单位正法向量n=(x,y,0),当(x,y,z)Σ+。曲面Σ在柱面坐标下的方程为(75)x=cosφ, y=sinφ, z=z,0φ2π, 0z2r=(x,y,z),则rφ=(sinφ,cosφ,0)rz=(0,0,1)。于是rφ×rz=(cosφ,sinφ,0)Σ+的单位正法向量n=(x,y,0)一致。因此(76)I=Σ+VndS=0φ2π, 0z2V(r(φ,z))(rφ×rz)dφdz=0φ2π, 0z2(cosφ(sinφz),0,cosφsinφ)(cosφ,sinφ,0)dφdz=0φ2π, 0z2(cos2φsinφzcos2φ)dφdz=2π

2 直接利用楔积计算可得(77)I=0φ2π,0z2[cosφ(sinφz)dsinφdz+(cosθsinφ)dcosθdsinθ]=0φ2π,0z2cos2φ(sinφz)dφdz=0φ2π,0z2zcos2φdφdz=2π. 在点(1,0,0)处,绕z轴逆时针旋转(φ增加)的切向量为(0,1,0),沿z轴上升(z增加)的切向量为(0,0,1)和曲面外法向(1,0,0)构成右手系,所以dφdz=dφdz

例 9.3.15 (例5) 求向量场V=xy i^+yz j^+zx k^从里向外穿过第一卦限中球面x2+y2+z2=1时的通量。

1 球面参数方程为(78)x=cosφsinθ, y=sinφsinθ, z=cosθ,(θ,φ)D=[0,π2]2 通量为(79)I=Σ+(xydydz+yzdzdx+zxdxdy) 向南-向东-向上构成球面的右手系,所以dθdφ符合球面从里向外。考虑到x,y,z 的轮转对称(这对称性是保方向的),所以(80)I=3Σ+zxdxdy=3Σ+zdx22dy=32Σ+zd(1y2z2)dy=3Σ+z2dzdy=3Dcos2θdcosθd(sinφsinθ)=3Dcos2θsinθcosφsinθdθdφ 所以(81)I=30π2cosφdφ0π21cos4θ8dθ=3π16

2 用Gauss公式,取(82)Ω={(x,y,z)x,y,z0,x2+y2+z21}(83) d(xydydz+yzdzdx+zxdxdy)=d(xy)dydz+d(yz)dzdx+d(zx)dxdy=ydxdydz+zdydx+xdzdxdy=(x+y+z)dxdydz 因此(84)Ω(xydydz+yzdzdx+zxdxdy)=Ω(y+z+x)dxdydz 由对称性可得(85)I=3Ωxdxdydz=301xdxy2+z21x2,y,z0dydz=301π(1x2)4xdx=3π16Ω的三个坐标平面上,由(86)(xydydz+yzdzdx+zxdxdy)=0 所以(87)Σ+ω=Ωω=3π16

例 9.3.16 (例6) Σx2+y2+z2=R2的上侧,求(88)Σ(x2dydz+y2dzdx+z2dxdy)

1 关于平面yz的反射:(x,y,z)(x,y,z)是一个反射变换,且使Σ改变方向,所以(89)Σx2dydz=Σ(x)2dydz=Σx2dydz=0 类似可得Σy2dzdx=0,所以(90)I=Σz2dxdy=x2+y2R2(R2x2y2)dxdy=πR42

2 注意到z=R2x2y2,计算可得(91)dydz=dydR2x2y2=dyd(R2x2y2)2R2x2y2=xdxdyR2x2y2,dzdx=dR2x2y2dx=d(R2x2y2)dx2R2x2y2=ydxdyR2x2y2, 代入积分中可得(92)I=Dx3+y3R2x2y2+(R2x2y2)dxdy=D(R2x2y2)dxdy=πR42 式中用到了函数奇偶性和重积分的对称性。

例 9.3.17 (例7) Σx2a2+y2b2+z2c2=1的外侧,计算:(93)Σ(x3dydz+y3dzdx+z3dxdy)

x=aX,y=bY,z=cZ,则Σ1:X2+Y2+Z2=1且为外侧,代入积分中可得(94)I=Σ(x3dydz+y3dzdx+z3dxdy)=Σ1(a3bcX3dYdZ+ab3cY3dZdX+abc3Z3dXdY) 变换(X,Y,Z)(Y,Z,X)保向,且保球面,所以(95)Σ1X3dYdZ=Σ1Y3dZdX=Σ1Z3dXdY 所以(96)I=(a3bc+ab3c+abc3)Σ1Z3dXdY=abc(a2+b2+c2)Σ1Z3dXdY=2abc(a2+b2+c2)X2+Y21(1X2Y2)32dXdY=abc(a2+b2+c2)2π01(1r2)32dr2=4πabc(a2+b2+c2)5

9.3.4 曲面积分中的坐标变换

引理 9.3.18 Σ为正则参数曲面,φ:R3R3是微分同胚(记φ,φ1都可微)且Jφ始终为正交矩阵。记Σ=φ(Σ)X=(x,y,z)TU=(u,v,w)T,则对任何Σ上的连续函数g(X),证明:(97)Σg(X)dσ=Σg(φ1(U))dσ

证明Σ有正则参数表示(98)X(s,t)=(x(s,t)y(s,t)z(s,t)),(s,t)D 其中DR2为平面有界闭区域。由此得到曲面Σ=φ(Σ)的一个参数表示(99)U(s,t)=φ(X(s,t)) 于是(100)dσ=det[((u,v,w)(s,t))T((u,v,w)(s,t))]dsdt=det[((x,y,z)(s,t))T((u,v,w)(x,y,z))T((u,v,w)(x,y,z))((x,y,z)(s,t))]dsdt=det[((x,y,z)(s,t))T((x,y,z)(s,t))]dsdt=dσ 故有(101)Σg(X)dσ=Dg(X(s,t))det[((x,y,z)(s,t))T((x,y,z)(s,t))]dsdt=Dg(φ1(U(s,t)))det[((u,v,w)(s,t))T((u,v,w)(s,t))]dsdt=Σg(φ1(U))dσ

例 9.3.19 (例8) 设一元函数f(u)于整个实轴上连续,S代表单位球面x2+y2+z2=1。证明Poisson公式:(102)Sf(ax+by+cz)dS=2π11f(ρu)du 这里ρ=a2+b2+c2

证明 显然ρ=0时结论成立,下设ρ>0。设P是正交矩阵,其第一行为(aρ,bρ,cρ)。在正交变换U=PX下,ax+by+cz=ρu,单位球面x2+y2+z2=1仍为单位球面u2+v2+w2=1。在以u轴为旋转轴的柱坐标系中,球面参数方程为(103)u=u,v=1u2cosθ,w=1u2sinθ 由此得到(104)dS=dudθ 根据上述引理可知(105) x2+y2+z2=1f(ax+by+cz)dS=u2+v2+w2=1f(ρu)dS=u[1,1],θ[0,2π]f(ρu)dudθ=2π11f(ρu)du.

2有关向量分析的更多细节,将在“向量场的旋度和散度”一节中详细介绍。