第十章 多元函数积分学中的基本公式及其应用
填空题
设 L 是区域
D:x2+y2≤−2x的正向边界,则
I=∮L(x3−y)dx+(x−y3)dy=.令
P(x,y)=x3−y,Q(x,y)=x−y3.曲线 L 是区域 D 的正向边界,由格林公式,
∮LPdx+Qdy=∬D(∂x∂Q−∂y∂P)dxdy.计算得
∂x∂Q=1,∂y∂P=−1,故
∂x∂Q−∂y∂P=2.又
x2+y2≤−2x⟺(x+1)2+y2≤1,所以 D 是半径为 1 的圆域,其面积为 π。因此
I=∬D2dxdy=2π.设 L 是平面上从圆周
x2+y2=a2上任一点到圆周
x2+y2=b2上任一点的一条光滑曲线 (a>0,b>0),r=x2+y2,则
I=∫Lr3(xdx+ydy)=.由
r2=x2+y2两边微分,得
2rdr=2xdx+2ydy,即
xdx+ydy=rdr.于是
r3(xdx+ydy)=r4dr=d(5r5).因此该曲线积分与路径无关,只取决于起点和终点的半径。起点位于半径为 a 的圆周上,终点位于半径为 b 的圆周上,所以
I=5r5r=ar=b=5b5−a5.设 r=x2+y2,常数 λ 使得曲线积分
∮Lyxrλdx−y2x2rλdy=0对上半平面的任何光滑闭曲线 L 成立,则 λ=。
上半平面 y>0 是单连通区域。令
P(x,y)=yxrλ,Q(x,y)=−y2x2rλ.要使曲线积分对上半平面的任意光滑闭曲线均为零,需满足
∂y∂P=∂x∂Q.注意
∂x∂rλ=λxrλ−2,∂y∂rλ=λyrλ−2.因此
∂y∂P=−y2xrλ+λxrλ−2,而
∂x∂Q=−y22xrλ−y2λx3rλ−2.令二者相等并整理,得
y2xrλ+λxrλ−2+y2λx3rλ−2=0.提取公因式后,
xrλ−2[y2r2+λ+y2λx2]=0.由于 r2=x2+y2,方括号内为
y2x2+y2+λy2x2+y2=(1+λ)y2r2.要使等式在上半平面恒成立,必须有
1+λ=0,故
λ=−1.求下列平面上曲线积分
其中 L 为椭圆
a2x2+b2y2=1的右半部分,从 A(0,−b) 到 B(0,b)。
注意到
d[ycos(x2)]=−2xysin(x2)dx+cos(x2)dy.因此原积分可写为
I=∫Ly2dx+∫Ld[ycos(x2)].第二项只与端点有关,故
∫Ld[ycos(x2)]=ycos(x2)AB=b−(−b)=2b.对于第一项,将椭圆右半部分参数化为
x=acost,y=bsint,−2π≤t≤2π.此时
dx=−asintdt,所以
∫Ly2dx=−ab2∫−2π2πsin3tdt=0,因为 sin3t 是奇函数,积分区间关于原点对称。
综上,
I=0+2b=2b.其中 A(0,−1),B(1,0),\wideparenAB 为单位圆在第四象限部分。
单位圆第四象限从 A(0,−1) 到 B(1,0) 可参数化为
x=cost,y=sint,−2π≤t≤0.于是
dx=−sintdt,dy=costdt.分子为
xdy−ydx=cos2tdt+sin2tdt=dt.分母为
(x−y)2=(cost−sint)2=2cos2(t+4π).因此
I=∫−2π0(cost−sint)2dt=21∫−2π0sec2(t+4π)dt=21tan(t+4π)−2π0=21(1−(−1))=1.其中 \wideparenAB 是沿椭圆
a2x2+b2y2=1正向从 A(a,0) 到 (0,b) 的一段弧,a=1。
注意到
d[21ln((x−1)2+y2)]=(x−1)2+y2(x−1)dx+ydy.由于 a=1,起点 A(a,0) 不是奇点 (1,0),所给第一象限椭圆弧上也不经过 (1,0)。因此该积分可直接利用原函数计算:
I=21ln((x−1)2+y2)(a,0)(0,b)=21ln(1+b2)−21ln(a−1)2=21ln(a−1)21+b2.也可写成
I=ln∣a−1∣1+b2.其中 L 是椭圆周
a2x2+b2y2=1,取逆时针方向。
将被积表达式拆分为
x2+y2x−y2dx+x2+y2x2+ydy=x2+y2xdx+ydy+x2+y2−y2dx+x2dy.第一部分为全微分:
x2+y2xdx+ydy=d(21ln(x2+y2)).由于 L 是闭曲线,故
∮Lx2+y2xdx+ydy=0.对第二部分,将椭圆参数化为
x=acost,y=bsint,0≤t≤2π.于是
dx=−asintdt,dy=bcostdt,从而
∮Lx2+y2−y2dx+x2dy=∫02πa2cos2t+b2sin2tab2sin3t+a2bcos3tdt.记被积函数为 F(t)。由于
F(t+π)=−F(t),故
∫02πF(t)dt=0.因此
I=0+0=0.其中
L:{x=a(t−sint),y=a(1−cost),t 从 0 到 π,a>0.由全微分公式,
d(exsiny)=exsinydx+excosydy,且
d(xy)=ydx+xdy.因此
(exsiny−my−y)dx+(excosy−mx)dy=d(exsiny)−md(xy)−ydx.曲线起点和终点分别为
A=(0,0),B=(aπ,2a).故
I=[exsiny−mxy]AB−∫Lydx=eaπsin(2a)−2mπa2−∫Lydx.由参数方程,
dx=a(1−cost)dt,y=a(1−cost),所以
∫Lydx=a2∫0π(1−cost)2dt=a2∫0π(1−2cost+cos2t)dt=a2(π+2π)=23πa2.因此
I=eaπsin(2a)−2mπa2−23πa2=eaπsin(2a)−(2m+23)πa2.求下列曲面积分
其中 Σ 为由曲面
y=x2+z2与平面 y=1、y=2 所围立体表面的外侧。
记
r=x2+z2.所围立体为
Ω:1≤y≤2,0≤r≤y.闭曲面 Σ 由三部分组成:
- 上底面 Σ1:y=2,外法向量指向 y 轴正向;
- 下底面 Σ2:y=1,外法向量指向 y 轴负向;
- 侧面 Σ3:y=x2+z2,外法向量指向立体外部。
由于被积函数在 r=0 处含有奇性,直接分片计算较为稳妥。
在上底面 Σ1 上,y=2,且投影区域为
x2+z2≤2.故
I1=∬x2+z2≤2x2+z2e2dxdz=∫02π∫02re2rdrdθ=22πe2.在下底面 Σ2 上,外法向量沿 y 轴负向,故
I2=−∬x2+z2≤1x2+z2edxdz=−∫02π∫01rerdrdθ=−2πe.侧面可写为
y=x2+z2.立体位于 x2+z2−y≤0 一侧,因此外法向量可取与
(2x,−1,2z)同向。其 y 分量为负,故
dzdx=−dxdz.侧面在 xOz 平面上的投影为
1≤x2+z2≤2.于是
I3=−∬1≤x2+z2≤2x2+z2ex2+z2dxdz=−∫02π∫12rerrdrdθ=−2π(e2−e).因此
I=I1+I2+I3=22πe2−2πe−2π(e2−e)=2π(2−1)e2.其中 Σ1 为圆柱面
x2+y2=a2上 x≥0、0≤z≤1 部分,法向量与 x 轴正向成锐角;Σ2 为 Oxy 平面上半圆域
x2+y2≤a2,x≥0部分,法向量与 z 轴正向相反。
令
F=((z+1)x,xy,0).则所求积分是 F 通过 Σ1∪Σ2 的通量。
在 Σ2 上,法向量为
n=(0,0,−1),而 F 的第三个分量为零,因此
∬Σ2F⋅ndS=0.对半圆柱面 Σ1,取参数方程
x=acosθ,y=asinθ,−2π≤θ≤2π,0≤z≤1.由于要求法向量与 x 轴正向成锐角,取单位外法向量
n=(cosθ,sinθ,0),dS=adθdz.于是
F⋅ndS=[(z+1)xcosθ+xysinθ]adθdz=[a2(z+1)cos2θ+a3cosθsin2θ]dθdz.因此
I=a2∫01(z+1)dz∫−2π2πcos2θdθ+a3∫01dz∫−2π2πcosθsin2θdθ.其中
∫01(z+1)dz=23,∫−2π2πcos2θdθ=2π,且令 u=sinθ,得
∫−2π2πcosθsin2θdθ=∫−11u2du=32.故
I=43πa2+32a3.S 是
a2x2+b2y2+c2z2=1(z≥0)的上侧。
令
F=(x2−y2,y2−z2,z2−x2).以椭圆域
S0:z=0,a2x2+b2y2≤1封闭上半椭球,并令 S0 取向下的外侧法向量。由高斯公式,
∬SF⋅ndS+∬S0F⋅ndS=∭ΩdivFdV,其中 Ω 为上半椭球体。
先求散度:
divF=2x+2y+2z.由于 Ω 关于 x=0、y=0 对称,故
∭Ω2xdV=0,∭Ω2ydV=0.于是
∭ΩdivFdV=∭Ω2zdV.当 z∈[0,c] 时,截面椭圆的面积为
πab(1−c2z2).因此
∭Ω2zdV=∫0c2z⋅πab(1−c2z2)dz=2πabc2.在底面 S0 上,向下的单位法向量为
n=(0,0,−1).又因 z=0,故
F⋅n=−(z2−x2)=x2.于是
∬S0F⋅ndS=∬Dx2dxdy,其中
D:a2x2+b2y2≤1.令
x=arcosθ,y=brsinθ,0≤r≤1,0≤θ≤2π,则
∬Dx2dxdy=a3b∫02πcos2θdθ∫01r3dr=4πa3b.所以
I=2πabc2−4πa3b=4πab(2c2−a2).求下列空间中的曲线积分
其中 L 是曲线
{x2+y2=4y,3y−z+1=0,且顺着 z 轴的正向看是沿逆时针方向。
令
F=(yz,3zx,−xy).曲线 L 位于平面
z=3y+1内,其在 Oxy 平面上的投影为圆
x2+(y−2)2=4.题目规定从 z 轴正向观察为逆时针,因此取曲面法向量的 z 分量为正。
计算旋度:
rotF=(∂y∂(−xy)−∂z∂(3zx),∂z∂(yz)−∂x∂(−xy),∂x∂(3zx)−∂y∂(yz))=(−4x,2y,2z).取 L 所围的平面区域为积分曲面 S。因为
z=3y+1,其向上的有向面积元为
ndS=(0,−3,1)dxdy.由斯托克斯公式,
I=∬SrotF⋅ndS=∬D(−4x,2y,2z)⋅(0,−3,1)dxdy=∬D(−6y+2z)dxdy.代入 z=3y+1,得
−6y+2z=−6y+2(3y+1)=2.投影区域 D 是半径为 2 的圆域,面积为 4π。故
I=2⋅4π=8π.其中 Γ 是沿螺旋线
x=acosθ,y=asinθ,z=2πhθ,从 A(a,0,0) 到 B(a,0,h) 的有向曲线。
记
P=x2−yz,Q=y2−xz,R=z2−xy.检验交叉偏导数:
Py=−z=Qx,Pz=−y=Rx,Qz=−x=Ry.因此微分形式
Pdx+Qdy+Rdz在整个空间内为全微分,曲线积分与路径无关。
求其原函数。由
ux=x2−yz积分得
u=3x3−xyz+g(y,z).再由
uy=−xz+gy=y2−xz可得
gy=y2,g=3y3+h(z).最后由
uz=−xy+h′(z)=z2−xy得
h′(z)=z2,h(z)=3z3.故可取原函数
u(x,y,z)=3x3+y3+z3−xyz.于是
I=u(B)−u(A)=(3a3+h3)−3a3=3h3.平面上第二类曲线积分与路径无关问题
判断下列曲线积分在指定区域 D 是否与路径无关,为什么?
(Ⅰ)
∫Lf(x2+y2)(xdx+ydy),其中 f(u) 为连续函数,D 为全平面;
(Ⅱ)
∫Lx2+y2−ydx+xdy,D={(x,y)∣全平面除去 −∞<x≤0, y=0}.(Ⅰ)与路径无关;
(Ⅱ)与路径无关。
(Ⅰ)令
u=x2+y2,则
du=2xdx+2ydy.因此
f(x2+y2)(xdx+ydy)=21f(u)du.取 f(u) 的一个原函数为 F(u),即
F′(u)=f(u),则
f(x2+y2)(xdx+ydy)=d[21F(x2+y2)].所以被积表达式在全平面上是某个函数的全微分,曲线积分只与起点、终点有关,而与路径无关。
也可取原函数
U(x,y)=21∫0x2+y2f(t)dt.(Ⅱ)记
P(x,y)=−x2+y2y,Q(x,y)=x2+y2x.在区域 D 内有
∂y∂P=(x2+y2)2y2−x2,∂x∂Q=(x2+y2)2y2−x2.故
Py=Qx.区域 D 是去掉非正 x 轴后的平面,它是单连通区域。可在该区域内取连续的辐角函数
θ=Arg(x+iy),−π<θ<π.由极坐标关系可知
dθ=x2+y2−ydx+xdy.所以该微分形式在 D 内为全微分,曲线积分与路径无关。
需要注意:若定义域是整个去心平面 R2∖{(0,0)},则区域不是单连通的,沿绕原点一周的闭曲线积分为 2π,此时曲线积分不与路径无关。题中删去了一条从原点延伸到无穷远的射线,正是为了使辐角函数能够单值连续地定义。
设 φ(x) 在 (0,+∞) 有连续导数,φ(π)=1,试确定 φ(x),使积分
I=∫AB[sinx−φ(x)]xydx+φ(x)dy在 x>0 与路径无关,并求当 A,B 分别为 (1,1)、(π,π) 时的积分值。
且
I=1+cos1.记
P(x,y)=[sinx−φ(x)]xy,Q(x,y)=φ(x).区域 x>0 为单连通区域。曲线积分与路径无关的充要条件为
Py=Qx.计算得
Py=xsinx−φ(x),Qx=φ′(x).故 φ(x) 应满足
φ′(x)=xsinx−φ(x),即
xφ′(x)+φ(x)=sinx.左端恰为乘积的导数:
[xφ(x)]′=sinx.积分得
xφ(x)=−cosx+C.由条件 φ(π)=1,有
π=−cosπ+C=1+C,从而
C=π−1.因此
φ(x)=xπ−1−cosx.由微分方程可知
φ′(x)=xsinx−φ(x),所以
P(x,y)=yφ′(x),Q(x,y)=φ(x).于是
Pdx+Qdy=d[yφ(x)].可取原函数
U(x,y)=yφ(x).故
I=U(π,π)−U(1,1)=πφ(π)−φ(1)=π−(π−1−cos1)=1+cos1.设
(P(x,y),Q(x,y))=((x2+y2)nx−y,(x2+y2)nx+y),n 为常数,问
∫LPdx+Qdy在区域
D={(x,y)∣(x,y)∈R2,(x,y)=(0,0)}是否与路径无关。
无论 n 取何值,该曲线积分在区域 D 内都不与路径无关。
令
s=x2+y2.则
P=(x−y)s−n,Q=(x+y)s−n.计算偏导数:
Py=−s−n−2ny(x−y)s−n−1,Qx=s−n−2nx(x+y)s−n−1.于是
Qx−Py=2s−n+2n[y(x−y)−x(x+y)]s−n−1=2s−n−2nss−n−1=2(1−n)s−n.当 n=1 时,
Qx−Py=0,因此微分形式不是闭形式,曲线积分不可能与路径无关。
下面考察 n=1。此时
\displaystyle \begin{aligned} P\,dx+Q\,dy &=\dfrac{x-y}{x^2+y^2}\,dx+rac{x+y}{x^2+y^2}\,dy\\ &=\dfrac{x\,dx+y\,dy}{x^2+y^2} +\dfrac{-y\,dx+x\,dy}{x^2+y^2}. \end{aligned}其中第一项为
x2+y2xdx+ydy=d(21ln(x2+y2)),第二项在局部等于 dθ,但在去心平面 D 内不存在单值的全局辐角函数。
取单位圆周
L:x=cost,y=sint,0≤t≤2π,沿逆时针方向,则
xdx+ydy=0,且
−ydx+xdy=dt.因此
∮LPdx+Qdy=∫02πdt=2π=0.若曲线积分与路径无关,则任意闭曲线积分都应为零,故当 n=1 时也不与路径无关。
综上,无论 n 取何值,该曲线积分在区域 D 内均不与路径无关。
原函数问题
求 Pdx+Qdy 在指定区域 D 上的原函数,其中
(P,Q)={1−x2y2cosxy, sinxy+xycosxy},D={(x,y)∣x>0}.原函数为
u(x,y)=x+ysinxy+C.要求函数 u(x,y),使
du=Pdx+Qdy,即
ux=1−x2y2cosxy,uy=sinxy+xycosxy.先对 Q 关于 y 积分。注意到
∂y∂(ysinxy)=sinxy+xycosxy.故
u(x,y)=ysinxy+φ(x).再对 x 求偏导:
ux=ycosxy(−x2y)+φ′(x)=−x2y2cosxy+φ′(x).与已知的 P 比较,得
φ′(x)=1,所以
φ(x)=x+C.因此原函数为
u(x,y)=x+ysinxy+C.直接验证:
d(x+ysinxy)=(1−x2y2cosxy)dx+(sinxy+xycosxy)dy.选择 a,b,使 Pdx+Qdy 在区域
D={(x,y)∣x2+y2=0}内为某函数 u(x,y) 的全微分,其中
P=(x2+y2)2y2+2xy+ax2,Q=−(x2+y2)2x2+2xy+by2.此时可取
u(x,y)=x2+y2x−y+C.若 Pdx+Qdy 是某函数的全微分,则必有
Py=Qx.计算并整理得
Py−Qx=−(x2+y2)34xy[(a+1)x+(b+1)y].要使上式对区域 D 中任意 (x,y) 恒为零,必须有
a+1=0,b+1=0.因此
a=b=−1.代入后,
P=(x2+y2)2y2+2xy−x2,Q=(x2+y2)2−x2−2xy+y2.观察函数
u(x,y)=x2+y2x−y.其偏导数为
ux=(x2+y2)2(x2+y2)−2x(x−y)=(x2+y2)2y2+2xy−x2=P,以及
uy=(x2+y2)2−(x2+y2)−2y(x−y)=(x2+y2)2−x2−2xy+y2=Q.因此
Pdx+Qdy=d(x2+y2x−y).虽然区域 D 不是单连通区域,但这里已经构造出了在整个 D 上单值可微的原函数,所以该微分形式确实是全微分。
求向量的散度与旋度
已知
E=r3qr,其中 r=(x,y,z),r=∣r∣,q 为常数,求 divE 与 rotE。
在 r=0 的区域内,
divE=0,rotE=0.原点处 E 无定义。
由
E=r3q(x,y,z),r=x2+y2+z2,可写成
E=(r3qx,r3qy,r3qz).先求散度。注意
∂x∂r−3=−3xr−5,∂y∂r−3=−3yr−5,∂z∂r−3=−3zr−5.因此
divE=q[∂x∂(xr−3)+∂y∂(yr−3)+∂z∂(zr−3)]=q[(r−3−3x2r−5)+(r−3−3y2r−5)+(r−3−3z2r−5)]=q[3r−3−3(x2+y2+z2)r−5]=q(3r−3−3r2r−5)=0.再求旋度。设
Ex=r3qx,Ey=r3qy,Ez=r3qz.例如,旋度的第一分量为
∂y∂Ez−∂z∂Ey=qz(−3yr−5)−qy(−3zr−5)=0.同理,
∂z∂Ex−∂x∂Ez=0,∂x∂Ey−∂y∂Ex=0.故
rotE=0.也可从势函数看出:当 r=0 时,
∇(r1)=−r3r,从而
E=−q∇(r1).因此 E 是梯度场,其旋度必为零。
需要注意,上述结论只在 r=0 时成立,因为原点处向量场无定义。