第一型曲线积分与曲面积分
第一型曲线积分
曲线的弧长
在定积分的几何应用一节我们已经介绍了曲线弧长公式,现在,我们对曲线的弧长作一个完整的论述,以引出第一型曲线积分的定义。
定义17.1 对1元
n维向量函数
ϕ(t)=(ϕ1(t),⋯,ϕn(t)),t∈[a,b],如果
ϕ(t)是连续的,则称
γ:{ϕ(t):t∈[a,b]}是
Rn上的连续曲线,如果对于
t1,t2∈[a,b],t1=t2,t1=a或t2=b,都有
ϕ(t1)=ϕ(t2),则称
γ为若当曲线或简单曲线,如果
γ是简单曲线同时
ϕ(a)=ϕ(b),则称
γ为若当闭曲线,如果
ϕ(t)在
[a,b]上有连续的导数,则称
γ为光滑的曲线
假设曲线
γ的起点和终点分别为
A,B,在中间取
n个点
A1,⋯,An,可将
γ分为
n+1段。令
A0=A,An+1=B,连接
Ak−1,Ak,k=1,⋯,n+1,得到一段内接折线。以内接折线的长度作为曲线弧长的估计。
L≈k=1∑n+1∣∣Ak−Ak−1∣∣这里
Ak(k=0,⋯,n+1)视为向量。由于两点之间线段最短,如果取另外一个分划
P2:A0′=A,A1′,⋯,Am′,Am+1′=B,如果
P1:A0,⋯,An+1都在
P2内,则称
P2是
P1的加细,则
k=1∑m+1∣∣Ak′−Ak−1′∣∣≥k=1∑n+1∣∣Ak−Ak−1∣∣不断加细,如果有一个上确界,那么这个上确界就称为是曲线的弧长,如下图。
如果
A1′在
A0,A1之间,那么必有
∣∣A1′−A0∣∣+∣∣A1−A1′∣∣≥∣∣A1−A0∣∣。这就是为什么
k=1∑m+1∣∣Ak′−Ak−1′∣∣≥k=1∑n+1∣∣Ak−Ak−1∣∣
定义17.2
γ是一段
Rn上的连续曲线,如果上确界
L=Δ:A0,⋯,An,⋯,An+1是γ的分划supk=1∑n+1∣∣Ak−Ak−1∣∣存在,则称
γ为可求长曲线,其称
L为
γ的弧长
在定积分一节中,我们求弧长的办法是,对连续曲线
γ:ϕ(t)=(ϕ1(t),⋯,ϕn(t)),t∈[a,b]对
[a,b]的分划
Δ:a=t0<⋯<tn<tn+1=b,则
L=λ(Δ)→0limk=1∑n+1∣∣ϕ(tk)−ϕ(tk−1)∣∣这个公式和定义17.2形式上不是完全一致的,但下面我们将证明这两个定义的完全相同的。
定理17.1
γ:ϕ(t),a≤t≤b是
Rn的一段可求长曲线,
L是其弧长,则
L=λ(Δ)→0limk=1∑n+1∣∣ϕ(tk)−ϕ(tk−1)∣∣
该定理的过程与达布定理及其类似。
证:
由
L的定义,对任意的
ε>0,存在分划
Δ0:a=t0(0)<t1(0)<⋯<tm(0)<tm+1(0)=b,满足
L≥k=1∑m+1∣∣ϕ(tk(0))−ϕ(tk−1(0))∣∣>L−2ε由于
ϕ(t)在
[a,b]上连续,则
ϕ(t)在
[a,b]上一致连续,存在
δ>0,当
x1,x2∈[a,b],
∣x1−x2∣<δ时,就有
∣∣ϕ(x2)−ϕ(x1)∣∣<2m(m+1)ε对任意的分划
Δ:a=t0<t1<⋯<ts=b,再设
Δ′=Δ0∪Δ:a=t0′<t1′<⋯<tp′<tp+1′=b,则
L≥≥k=1∑p+1∣∣ϕ(tk′)−ϕ(tk−1′)∣∣k=1∑m+1∣∣ϕ(tk(0))−ϕ(tk−1(0))∣∣>L−2ε从而
∣∣∣∣∣k=1∑p+1∣∣ϕ(tk′)−ϕ(tk−1′)∣∣−L∣∣∣∣∣<2ε只要
λ(Δ)<δ,
Δ中至多有
m个小区间插入了
Δ0的分点,一个小区间至多插入
m个
Δ0的分点,从而
=<∣∣∣∣∣k=1∑p+1∣∣ϕ(tk′)−ϕ(tk−1′)∣∣−k=1∑s+1∣∣ϕ(tk)−ϕ(tk−1)∣∣∣∣∣∣∣k=1∑p+1∣∣ϕ(tk′)−ϕ(tk−1′)∣∣−k=1∑s+1∣∣ϕ(tk)−ϕ(tk−1)∣∣m(m+1)2m(m+1)ε=2ε从而
≤+∣∣∣∣∣k=1∑s+1∣∣ϕ(tk)−ϕ(tk−1)∣∣−L∣∣∣∣∣∣∣∣∣∣k=1∑p+1∣∣ϕ(tk′)−ϕ(tk−1′)∣∣−k=1∑s+1∣∣ϕ(tk)−ϕ(tk−1)∣∣∣∣∣∣∣∣∣∣∣∣k=1∑p+1∣∣ϕ(tk′)−ϕ(tk−1′)∣∣−L∣∣∣∣∣<ε因此
L=λ(Δ)→0limk=1∑n+1∣∣ϕ(tk)−ϕ(tk−1)∣∣
模仿定积分几何应用中
R2光滑曲线弧长的求法,可以证明
Rn中的光滑曲线
γ:γ(t),t∈[a,b]都是可求长曲线,并且
L=∫ab∣∣γ′(t)∣∣dt于是,如果连续曲线由有限段光滑曲线拼接而成,那么该连续曲线也是可求长曲线。
第一型曲线积分的物理背景及定义
对于
R3上的一条细钢丝
γ,在其上定义了一个密度函数
ρ(x,y,z),怎么求其质量呢?如果钢丝是均匀的,那么其质量应该是
ρ.L(γ),其中,
ρ为钢丝的密度,
L(γ)是钢丝的长度。如果钢丝不是均匀的,可以采取微元法:将钢丝分解为若干段小钢丝
γ1,⋯,γn,只要每段钢丝的弧长足够小,如果
ρ是连续的,在每段小钢丝就可以近似地视为均匀的小钢丝,任取
ξk∈γk(k=1,⋯,n),则估计其质量为
m(γ)≈k=1∑nρ(ξk)L(γk)当钢丝越分越细时,误差越来越小,则小钢丝的质量就为
m(γ)=λ(Δ)→0limk=1∑nρ(ξk)L(γk)其中,
λ(Δ)=1≤i≤nmaxL(γk),求解的思路和定积分是相同的,不同的是现在是对曲线的积分。对以上物理背景进行抽象,就得到第一型曲线积分的定义。首先,如果
γ是可求长曲线,则连续曲线
γk⊆γ也是可求长的,这由可求长曲线的定义是容易看出的。
定义17.3
γ是
Rn上一段可求长曲线,起点和终点分别为
A,B,
ρ(x)是
γ上的函数,如果存在实数
I,对任意的
ε>0,存在
δ>0,从
A到
B任取
A=A0,A1,⋯,An,An+1=B,只要
λ=1≤i≤n+1maxL(Ai−1Ai)<δ,任取
ξk∈Ak−1Ak(k=1,⋯,n+1),就有
∣∣∣∣∣k=1∑n+1ρ(ξk)L(Ak−1Ak)−I∣∣∣∣∣<ε则称
ρ在
γ上可积,
I为
ρ在
γ上的第一型曲线积分,记为
∫γf(x)ds
第一型曲线积分和定积分、重积分一样,有线性性质,不等式性质,区间可加性。
定理17.2
γ是
Rn上的可求长曲线,
f1,f2在
γ上可积,则对于任意的实数
α,β,
αf1+βf2在
γ上也可积,并且
∫γ(αf1(x)+βf2(x))ds=α∫γf1(x)ds+β∫γf2(x)ds
定理17.3
γ是
Rn上的可求长曲线,
f(x)=1,x∈γ在
Rn上可积,并且
∫γf(x)ds=L(γ)
定理17.4
γ是
Rn上的可求长曲线,
f1,f2在
γ上可积,并且
f1(x)≤f2(x)x∈γ,则
∫γf1(x)ds≤∫γf2(x)ds
定理17.5
γ是
Rn上的可求长曲线,起点和终点分别为
A,B,取一个分点
C,
f在
γ上可积,则
f在
AC段和
CB段都可积,并且
∫γf(x)ds=∫ACf(x)ds+∫CBf(x)ds
第一型曲线积分的计算公式
第一型曲线积分的计算也是通过定积分进行。
定理17.6 光滑曲线
γ:ϕ(t),t∈[a,b],设
ϕ′(t)=0,
f(x)在
γ上连续,则
f(x)在
γ上可积,并且
∫γf(x)ds=∫abf(ϕ(t))∣∣ϕ′(t)∣∣dt
证:
令
S(x1,x2)=∫x1x2∣∣ϕ′(t)∣∣dt,则设
m=t∈[a,b]min∣∣ϕ′(t)∣∣,M=t∈[a,b]max∣∣ϕ′(t)∣∣,有
m>0,因此
m(x2−x1)≤S(x1,x2)≤M(x2−x1)对
Δ:a=t0<t1<⋯<tn=b,对应
γ的分划
Δ′:A0=ϕ(a),A1=ϕ(t1),⋯,An=ϕ(tn)=B,由上式可以看出,
λ(Δ′)→0的充要条件是
λ(Δ)→0。对任意的分划
Δ:a=t0<⋯<tn=b,任取
ξk∈[tk−1,tk](k=1,⋯,n),由积分中值定理,存在
ζk∈[tk−1,tk](k=1,⋯,n),有
=k=1∑nf(ξk)[S(tk)−S(tk−1)]=k=1∑nf(ξk)∫tk−1tk∣∣ϕ′(t)∣∣dtk=1∑nf(ξk)∣∣ϕ′(ζk)∣∣Δtk
∀ε>0,设
M=t∈[a,b]max∣f(ϕ(t))∣,由于
∣∣ϕ′(t)∣∣在
[a,b]上一致连续,存在
δ1>0,对任意的
x1,x2∈[a,b],只要
∣x1−x2∣<δ1,就有
∣∣∣ϕ′(x1)∣∣−∣∣ϕ′(x2)∣∣∣<2M(b−a)ε令
I=∫abf(ϕ(t))∣∣ϕ′(t)∣∣dt,存在
δ2>0,就有
∣∣∣∣∣k=1∑nf(ϕ(ξk))∣∣ϕ′(ξk)∣∣Δtk−I∣∣∣∣∣<2ε只要
λ(Δ)<δ0=min(δ1,δ2),就有
≤≤∣∣∣∣∣k=1∑nf(ϕ(ξk))∣∣ϕ′(ξk)∣∣Δtk−k=1∑nf(ϕ(ξk))∣∣ϕ′(ζk)∣∣Δtk∣∣∣∣∣k=1∑n∣f(ϕ(ξk))∣∣∣∣ϕ′(ξk)∣∣−∣∣ϕ′(ζk)∣∣∣ΔtkM(b−a)2M(b−a)ε=2ε则
≤<∣∣∣∣∣k=1∑nf(ϕ(ξk))∣∣ϕ′(ζk)∣∣Δtk−I∣∣∣∣∣∣∣∣∣∣k=1∑nf(ϕ(ξk))∣∣ϕ′(ξk)∣∣Δtk−I∣∣∣∣∣+∣∣∣∣∣k=1∑nf(ϕ(ξk))∣∣ϕ′(ξk)∣∣Δtk−k=1∑nf(ϕ(ξk))∣∣ϕ′(ζk)∣∣Δtk∣∣∣∣∣ε因此
f(x)在
γ上可积,并且
∫γf(x)ds=∫abf(ϕ(t))∣∣ϕ′(t)∣∣dt
由第一型曲线积分的区间可加性,我们就可以计算逐段光滑的曲线上的曲线积分。计算第一型曲线积分的第一步,首先是要写出曲线的参数方程,然后再套用以上公式即可。
例17.1 计算
∫γ(x34+y34)ds,其中
γ为曲线
x32+y32=a32(a>0)
解:
将
γ写成参数方程形式
{x=acos3ty=asin3t其中
t∈[−π,π]则
{x′=−3asintcos2ty′=3acostsin2t则
x′2+y′2
=3a∣sintcost∣则
=∫γ(x34+y34)ds=∫−ππa34(sin4t+cos4t)(3a∣sintcost∣)dt3a37∫−ππ(sin4t+cos4t)∣sintcost∣dt=4a37
例17.2 求
∫Lxyds,其中
L为球面
x2+y2+z2=a2与平面
x+y+z=0的交线
解:
为了写出
L的参数方程,我们首先要求出
L在
Oxy平面上的投影,联立
{x2+y2+z2=a2x+y+z=0消去
z,得到
x2+y2+(−x−y)2=2x2+2xy+2y2=a2通过配方,得到
2(x+2y)2+23y2=a2通过变数替换后投影可以变换为椭圆,无论如何,从以上方程我们可以令
⎩⎪⎪⎨⎪⎪⎧x+2y=22
acosty=32
asintz=−x−y由此就可以得到
L的参数方程
⎩⎪⎨⎪⎧x=22
acost−66
asinty=36
asintz=−22
acost−66
asint其中
t∈[−π,π],则
x′2+y′2+z′2
=a从而
∫Lxyds==a∫−ππ(22
acost−66
asint)(36
asint)dt−3a3∫−ππsin2tdt=−3a3π
例17.3 计算
∫L(xy+xz+yz)ds,
L同例17.2
解:
注意到
L上
x,y,z的地位是相同,比如,求解
∫Lxzds,将参数方程写成
⎩⎪⎨⎪⎧x=22
acost−66
asintz=36
asinty=−22
acost−66
asint再代入就可以得到
∫Lxzds=−3a3π,从而
∫L(xy+xz+yz)ds=−a3π
例17.4 计算
∫L2y2+z2
ds,其中
L是
x2+y2+z2=a2(a>0)与
x=y的交线
解:
在
L上,有
2y2+z2=a2,有
∫L2y2+z2
ds=a2∫Lds=πa2
第一型曲面积分
第一型曲面积分的物理背景及定义
同第一型曲面积分的物理背景类似,第一型曲面积分的物理背景是空间曲面的质量。如果
S是可求面积的均匀的空间曲面,设其密度为
ρ,那么
S的质量为
ρ∣S∣,其中
∣S∣是
S的面积。但如果
S不是均匀的空间区间,在
S上定义了密度函数
ρ(x,y,z),可以将
S划分为可求面积的小区面块
S1,⋯,Sn,只要
1≤i≤nmaxdiam(Si)足够小,由一致连续性,
S1,⋯,Sn都可以视为均匀曲面,任取
ξk∈Sk(k=1,⋯,n),则估计其质量为
m(S)≈k=1∑nρ(ξk)∣Sk∣,当
1≤i≤nmaxdiam(Si)→0时,如果该和数有极限,则该极限为
S的质量。
定义17.4 设
S是一张可求面积的曲面,将
S分割为可求面积的小曲面块
S1,⋯,Sn,记该分划为
Δ,定义
λ(Δ)=1≤i≤nmaxdiam(Si)
f(x)是
S上的函数,如果存在实数
I,对任意的
ε>0,存在
δ>0,对任意的分划
Δ:S1,⋯,Sn,只要
λ(Δ)<δ,任取
ξk∈Sk(k=1,⋯,n),都有
∣∣∣∣∣k=1∑nf(ξk)∣Sk∣−I∣∣∣∣∣<ε则称
f在
S上可积,
I称为
f在
S上的第一型曲面积分,记为
∬Sf(x)ds
同样可以写出曲面积分的性质:线性,不等式,可加性等,这与曲线积分比较类似。
第一型曲面积分的计算公式
关于第一型曲面积分的存在性,可以模仿定积分和重积分的可积性理论,建立起曲面积分的可积性理论,就可以证明如下命题:
命题17.1
f(x,y,z)在光滑曲面
S上连续,则
f(x,y,z)在
S上可积
下面我们给出第一型曲面积分的计算公式:
定理17.6
f(x,y,z)在光滑曲面
⎩⎪⎨⎪⎧x=x(u,v)y=y(u,v)z=z(u,v)上连续,其中
(u,v)∈D,
D为可求面积的有界闭区域,其中
ru′×rv′=0,则
∬Sf(x,y,z)dS=∫Df(x(u,v),y(u,v),z(u,v))∣∣ru′×rv′∣∣dudv
该定理的证明和曲线积分的证明是类似的,需要用到重积分的积分中值定理,这里就不再赘述了。下面给出第一型曲面积分的几个算例。
例17.5 (利用曲面方程化简被积函数) 求解曲面积分
∬Sx2+y2dS,其中
S为柱面
x2+y2=R2被平面
z=0和
z=H所截的部分
解:
由曲面方程为
x2+y2=R2,得到
∬Sx2+y2dS=R21∬SdS=R22RHπ=R2Hπ
例17.6 求解曲面积分
∬Sz2dS,其中
S为
x=ucosv,y=usinv,z=v(0≤u≤a,0≤v≤2π)
解:
ru′=(cosv,sinv,0),rv′=(−usinv,ucosv,1),则曲面积分化为重积分为
∬z2dS==∬[0,a]×[0,2π]v2u2+1
dudv34π3(ln(a+1+a2
)+a1+a2
)
例17.7 (求解曲面积分时注意完整考虑整个曲面,不要遗漏某一两面)求解第一型曲面积分
∬Sx2+y2dS,其中
S为立体
x2+y2
≤z≤1的边界曲面
解:
需要注意的是,这个立体是一个倒圆锥,有一个底面和一个侧面,不要遗漏掉底面。设底面为
S1,侧面为
S2,则
S1:z=1,x2+y2≤1,从而
=∬S1x2+y2dS=∬x2+y2≤1x2+y2dS∫02πdθ∫01r3dr=2π侧面写成参数方程形式为
⎩⎪⎨⎪⎧x=rcosθy=rsinθz=r取值范围为
0≤r≤1,0≤θ≤2π,此时
A2+B2+C2
=2
r,则将曲面积分化为重积分即为
∬S2x2+y2dS=2
∫02πdθ∫01r3dr=22
π则
∬Sx2+y2dS=∬S1x2+y2dS+∬S2x2+y2dS=2(1+2
)π
第二型曲线积分与曲面积分
第二型曲线积分
第二型曲线积分的物理背景及定义
第二型曲线积分的物理背景是变力做功。
如上图,如果在牵引力
F的作用下,箱子移动的位移
s,则在力学中,力
F对箱子所作的功为
F.s,这是恒力对一个质点的作用。如果是变力,该如何求解力对质点所作的功呢。假设在平面上有一个力场
F(x,y)=(P(x,y),Q(x,y)),在力的作用下,质点的运动轨迹为
s:s(t)=(x(t),y(t)),t∈[α,β],这里假设运动轨迹是光滑的。则我们可以将运动曲线划分为
n段小曲线
s1,⋯,sn,每段运动都可以视为恒力做功,任取
(ξk,ζk)∈sk,在该段发生的位移为
Δsk,则估计该力所做的功为
W≈k=1∑nF(ξk,ζk).Δsk当曲线段最大直径趋于0时,如果以上和式右极限,即为力场对质点所作的功。将以上物理背景进行抽象,就得到第二型曲线积分。
定义17.5
L为
Rn上的连续曲线,起点为
A,终点为
B,
F(x)=(F1(x),⋯,Fn(x))是定义在
L上的
n维向量函数,如果存在实数
I,对任意的
ε>0,存在
δ>0,对任意的
L的分划
Δ:A=A0,A1,⋯,Am−1,Am=B,只要
λ(Δ)<δ,不论取何种
ξk∈Ak−1Ak,其中
Ak−1Ak表示
Ak−1到
Ak的弧段(
k=1,⋯,m),都有
∣∣∣∣∣k=1∑mF(ξk).(Ak−Ak−1)−I∣∣∣∣∣<ε在上式中,
Ak表示向量
(k=1,⋯,n),则称
F在
L上的第二型曲线积分存在,
I记为
∫Lk=1∑nFk(x)dxk
第二型曲线积分是有向的,计算第二型曲线积分时必须指定起点和终点,容易证明第二型曲线积分有线性性质,区间可加性,并且起终点互换时,第二型曲线积分的值为原来的值的相反数,具体的定理和证明就不一一写出了。
第二型曲线积分的计算方法
同样地,我们仅给出光滑曲线上的第二型曲线积分的计算公式。对光滑曲线
L:x(t)=(x1(t),⋯,xn(t))t∈[a,b]
F(x)=(F1(x),⋯,Fn(x))是定义在
L上的连续函数,对分划
Δ:a=t0<t1<⋯<tm=b,任取
ξk∈[tk−1,tk](k=1,⋯,m),得到一个和式
k=1∑mF(x(ξk)).(x(tk)−x(tk−1))=k=1∑mi=1∑nFi(x(ξk))(xi(tk)−xi(tk−1))由拉格朗日中值定理,对任意的
i=1,⋯,n,k=1,⋯,m,存在
ζik∈[tk−1,tk],满足:
xi(tk)−xi(tk−1)=xi′(ζik)Δtk代入,得到
k=1∑mF(x(ξk)).(x(tk)−x(tk−1))=k=1∑mi=1∑nFi(x(ξk))xi′(ζik)Δtk对任意的
ε>0,由于
F是连续的,设
M1=1≤k≤nmaxt∈[a,b]max∣Fi(x(t))∣,由一致连续性,存在
δ>0,当
x1,x2∈[a,b],∣x1−x2∣<δ1时,对任意的
i=1,⋯,n,都有
∣xi′(x1)−xi′(x2)∣<2M(b−a)ε又存在
δ2>0,当
λ(Δ)<δ2时,有
∣∣∣∣∣k=1∑mFi(x(ξk))xi′(ξk)Δtk−∫abFi(x(t))xi′(t)dt∣∣∣∣∣<2ε则当
λ(Δ)<min(δ1,δ2)时,有
≤<<∣∣∣∣∣∫abFi(x(t))xi′(t)dt−k=1∑mFi(x(ξk))xi′(ζik)Δtk∣∣∣∣∣∣∣∣∣∣k=1∑mFi(x(ξk))xi′(ζik)Δtk−k=1∑mFi(x(ξk))xi′(ξk)Δtk∣∣∣∣∣+∣∣∣∣∣∫abFi(x(t))xi′(t)dt−k=1∑mFi(x(ξk))xi′(ξk)Δtk∣∣∣∣∣2ε+k=1∑m∣Fi(x(ξk))∣∣xi′(ξk)−xi′(ζik)∣Δtk2ε+M2M(b−a)εk=1∑mΔtk=ε得出结论,对
i=1,⋯,n,有
λ(Δ)→0limk=1∑mFi(x(ξk))(xi(tk)−xi(tk−1))=∫abFi(x(t))xi′(t)dt故
λ(Δ)→0limi=1∑nk=1∑mFi(x(ξk))(xi(tk)−xi(tk−1))=i=1∑n∫abFi(x(t))xi′(t)dt而
1≤k≤mmaxdiam(sk)→0的充要条件就是
λ(Δ)→0。于是就有如下定理
定理17.7 对光滑曲线
L:x(t)=(x1(t),⋯,xn(t))t∈[a,b]
F(x)=(F1(x),⋯,Fn(x))是定义在
L上的连续函数,
F(x)在
L上的第二型曲线积分存在,并且
∫Lk=1∑nFk(x)dxk=k=1∑n∫abFi(x(t))xi′(t)dt
由定理17.7,实际上我们可以分开定义
∫LFi(x)dxi(i=1,⋯,n),但如果这样定义,就不能反映出第二型曲线积分的物理背景。但计算第二型曲线积分时,需要逐个计算
∫LFi(x)dxi(i=1,⋯,n),再求和。以下是第二型曲线积分的若干算例。
例17.8(确定曲线的方向) 计算第二型曲线积分
∫L(y+z)dx+(x+z)dy+(x+y)dz,其中
L为曲线
x2+y2=1,z=0,从
z轴正向看去取逆时针
解:
计算改积分的难点是确定曲线的方向,将曲线写成参数方程的形式
⎩⎪⎨⎪⎧x=costy=sintz=0方向如上图所示,则参数的取值范围是
[0,2π],起点为
t=0,终点为
t=2π,则
=∫L(y+z)dx+(x+z)dy+(x+y)dz∫02π(−sin2t+cos2t)dt=0
注:如何理解例17.8中的从z轴正向看去成逆时针方向:实际上,可以手比划出一个圈,从下网上看,在该视角下确定逆时针方向,标上箭头即可。当然,我们可以用右手螺旋法则确定该方向,伸出右手,拇指向上,另外四指所指的方向即为题目所指定的方向。
例17.9(利用投影法确定曲线参数方程) 计算第二型曲线积分
∫Lydx+zdy+xdz,其中
L为曲线
{x2+y2+z2=2zx+z=1从
z轴正向看去,
L取逆时针
解:
该题的难点在于确定曲线的参数方程,联立两个曲线,消去
z,得到
2x2+y2=1可见曲线在
Oxy平面上的投影是一个椭圆,由此,就可以很容易的写出该曲线的参数方程
⎩⎪⎨⎪⎧x=22
cosθy=sinθz=1−22
cosθ接下来,我们需要确定曲线的方向,由右手螺旋法则,
θ的取值范围是
[0,2π],起点是
θ=0,终点是
θ=2π,因此
==∫Lydx+zdy+xdz∫02π(−22
sin2θ+(1−22
cosθ)cosθ+21sinθcosθ)dθ−2
π
对于闭曲线,设其围成的图形为
S,指定其正向为向该方向移动,
S始终在其左手边,则对于闭曲线
L,记其正向的曲线积分为
∮LPdx+Qdy
例17.10 (利用曲线来简化被积函数)计算下列第二型曲线积分
∮Lx2+y2xdy−ydx,其中
L为
x2+y2=1
解:
在该曲线上,有
x2+y2=1,则
∮Lx2+y2xdy−ydx=∮Lxdy−ydx曲线的参数方程为
{x=cosθy=sinθ取值范围
[0,2π],起点为
θ=0,终点为
θ=2π,则
∮Lxdy−ydx=∫02π(cos2θ+sin2θ)dθ=2π
两类曲线积分的联系
第一型曲线积分和第二型曲线积分是有联系的。若光滑曲线
L:x(t),t∈[a,b],满足
x′(t)=0,∀t∈[a,b]则对
t∈[a,b],在该点处曲线的切向量为
(x1′(t),⋯,xn′(t)),单位向量为
(∑i=1nxi′2(t)
x1′(t),∑i=1nxi′2(t)
x2′(t),⋯,∑i=1nxi′2(t)
xn′(t))记
cosαi=∑i=1nxi′2(t)
xi′(t),若
fi(x1,⋯,xn)(i=1,⋯,n)在
L上连续,则
∫Li=1∑nfidxi=∫abi=1∑nficosαii=1∑nxi′2(t)
=∫Li=1∑nficosαids在这里
cosαi(i=1,⋯,n)是
L上的连续函数,这个公式就建立了第二型曲线积分和第一型曲线积分的联系,当然,第一型曲线积分是没有方向的,之所以第二型曲线积分有方向,是因为
cosαi(i=1,⋯,n)是由曲线的方向确定的。
第二型曲面积分
曲面的定向
在介绍第二型曲面积分的物理背景以及其定义之前,我们首先需要介绍曲面的定向。首先我们要给出可定向曲面和不可定向曲面的概念。所谓可定向曲面,又称双侧曲面,也就是对曲面上任意一点,都有两个方向相反的单位法向量,确定一个正方向后,绕着曲面环绕一周,单位法向量连续变化,回到原点后单位法向量应当同原来相同。举一个例子:圆柱面。下图是一个圆柱面
从
A点开始,绕着图中的闭曲面,环绕一周,
A处法向量确定为向外的法向量,则环绕一周后,法向量连续变化,回到
A点时,法向量和原来相同,此时称该曲面为可定向曲面或双侧曲面。不可定向曲面则相反,莫比乌斯带就是一个典型的单侧曲面或不可定向曲面。莫比乌斯带的构造如下,取一个细长的矩形纸带
ABB′A′,
A对应
B′,
B对应
A′将纸袋
AB边和
A′B′边拼接起来
沿着原来的直线
CC′运动,在
C处选定一个单位法向量,且连续变化,到达
C′点时,法向量则会变成原来的反方向。我们的讨论仅限于双侧曲面,对于双侧曲面,每一点都取一个正方向(要求这个单位向量作为曲面的向量函数在曲面上是连续的),这样就规定的曲面正侧,设此时单位向量为
(cosα,cosβ,cosγ)。则另一侧称为负侧,对应的法向量为
−(cosα,cosβ,cosγ)。当
cosγ>0,则称此时的正侧为上侧,反之称下侧,当
cosβ>0,则称此时的正侧为右侧,反之称左侧,当
cosα>0,则称此时的正侧为前侧,反之为后侧。对于封闭曲面,设其围成的立体为
V,若单位向量指向
V,称该侧为内侧,否则称外侧。
第二型曲面积分的物理背景及定义
第二型曲面积分的物理背景是通过一个空间曲面的流体的流量,设流体流动的速度向量为
v(x,y,z)=(P(x,y,z),Q(x,y,z),R(x,y,z))。如果流速是恒定的,通过一个平面区域
ΔS,该平面的法向量为
n,则流量应该为
v.n∣ΔS∣,其中
∣ΔS∣为
ΔS的面积,如果流速在空间上不是恒定的,则我们可以将曲面划分为若干小区块,只要最大直径足够小,小区块可以近似的认为是平直的,任取其中一点的单位法向量
n(ξk)=(cosαk,cosβk,cosγk),估计该曲面块上的流量为
n(ξk).v(ξk)∣ΔSk∣,则估计整个曲面的流量为
k=1∑mn(ξk).v(ξk)∣ΔSk∣,当
1≤i≤mmaxdiam(ΔSk)→0时,如果以上和式有极限存在,该极限就是单位时间内该流量场通过该空间曲面的流量。再考察以上和式
=k=1∑mn(ξk).v(ξk)∣ΔSk∣k=1∑mP(ξk)cosαk∣ΔSk∣+k=1∑mQ(ξk)cosβk∣ΔSk∣+k=1∑mR(ξk)cosγk∣ΔSk∣当
v连续时,三个曲面积分都存在
∬SPcosαdS,∬SQcosβdS,∬SRcosγdS并且
λ(Δ)→0limk=1∑mP(ξk)cosαk∣ΔSk∣=∬SPcosαdSλ(Δ)→0limk=1∑mQ(ξk)cosβk∣ΔSk∣=∬SQcosβdSλ(Δ)→0limk=1∑mR(ξk)cosγk∣ΔSk∣=∬SRcosγdS记
∬SPcosαdS=∬SPdydz,∬SQcosβdS=∬SQdzdx,∬SRcosγdS=∬SRdxdy。我们称这个积分是第二型曲面积分,从定义来看,第二型曲面积分是有方向的,方向就体现在正侧的选择上,第二型曲面积分和第一型曲面积分的联系是
∬SPdydz+Qdzdx+Rdxdy=∬SPcosα+Qcosβ+RcosγdS这与第一型曲线积分和第二型曲线积分的联系是一致的。
第二型曲面积分的计算方法
如果光滑曲面
⎩⎪⎨⎪⎧x=x(u,v)y=y(u,v)z=z(u,v)其中
(u,v)∈D,
D是可求面积的有界闭区域,满足
ru′×rv′=0,则对任意的
(u,v)∈D,在该点处的法向量为
±(A,B,C)A=∂(u,v)∂(y,z),B=∂(u,v)∂(z,x),C=∂(u,v)∂(x,y)选定其中一侧为正侧,设正侧的法向量为
(A,B,C),则第二型曲面积分的计算公式为
∬SPdydz+Qdzdx+Rdxdy=∬DPA+QB+RCdudv证明方法和第一型曲面积分公式的证明方法是类似的(利用积分中值定理和一致连续性),这里就不再具体写出,下面给出一个算例。
例17.11 (利用曲面特征化简被积函数)计算第二型曲面积分
I=∬S(x2+y2+z2)23xdydz+ydzdx+zdxdy,其中
S为
z=a2−x2−y2
的上侧
解:
由于在曲面上有
x2+y2+z2=a2,则
∬S(x2+y2+z2)23xdydz+ydzdx+zdxdy=a31∬Sxdydz+ydzdx+zdxdy曲面的参数方程为
⎩⎪⎨⎪⎧x=acosφcosθy=acosφsinθz=asinφ取值范围为
0≤φ≤2π,0≤θ≤2π,并且
rφ′=(−asinφcosθ,−asinφsinθ,acosφ),
rθ′=(−acosφsinθ,acosφcosθ,0),则
A=a2cos2φcosθB=a2cos2φsinθC=a2sinφcosφ代入
=∬Sxdydz+ydzdx+zdxdya3∫02πdφ∫02πdθ=2πa3因此
∬S(x2+y2+z2)23xdydz+ydzdx+zdxdy=2π
积分之间的联系
格林公式
定理17.8
P,Q在区域
Ω⊆R2上连续可导,
D⊆Ω是由有限条逐段光滑的简单闭曲线
L围成的有界闭区域,则有
∫LPdx+Qdy=∬D(∂x∂Q−∂y∂P)dxdy其中
L取正向
我们先来阐述几个概念,对于一条若当闭曲线
L,存在一个有界闭区域和一个无界闭区域,均以
L为边界,这就是若尔当定理。其具体的证明需要用到拓扑上的概念,这里我们就不给出证明了。但是这个定理的几何解释是相当直观的,见下图。
L是一条简单闭区间,
D1是一个有界区域,
D2是一个无界区域,但是两者均以
L为边界,现在我们引入单连通、多连通的概念,
D是一个区域或闭区域,如果在
D内任作若尔当闭曲线
L′,以
L′为边界的有界闭区域
D′满足
D′⊆D,则称
D是单连通区域,否则称为多连通区域,我们可以给出一个直观的几何解释。
如果区域
D内有一个“洞”,在这个“洞”外围画一条若尔当闭曲线,则其围成的有界闭区域不全在
D内,如上图。如果
D内没有“洞”,则无论在
D内如何画若当闭曲线,其围成的有界闭区域一定全部包含在
D内,记为单连通区域。所谓单连通区域,从几何直观讲,即是内部没有"洞"的区域。如果含有一个洞,上面左边的那幅图所示的区域,其边界是两条若当闭曲线,称为二连通区域,如果有两个洞,则边界由三条若当闭曲线构成,以此类推,如果边界由
k条若当闭曲线构成,则称该区域为
k连通区域。我们先就单连通区域讨论格林公式。先考察X型区域
D={(x,y):a≤x≤b,f(x)≤y≤g(x)}:
=∬D∂y∂Pdxdy=∫abdx∫f(x)g(x)∂y∂Pdy∫abP(x,g(x))dx−∫abP(x,f(x)dx其边界由四条曲线构成,如下图
设
f(x),g(x)都连续可微,则
L1的参数方程可表为
{x=ty=f(t)取值范围
[a,b],起点为
a,终点为
b,则
∫L1Pdx=∫abP(t,f(t))dt
L3的参数方程可表为
{x=ty=g(t)取值范围
[a,b],起点
b,终点
a,则
∫L3Pdx=∫baP(t,g(t))dt而
∫L2Pdx=∫L4Pdx=0于是
∫LPdx=∫abP(x,f(x))dx−∫abP(x,g(x))dx=−∬D∂y∂Pdxdy同样,对
Y型区域,成立
∫LQdy=∬D∂x∂Qdxdy如果
D即是
X型区域,又是
Y型区域,则
∫LPdx+Qdy=∬D(∂x∂Q−∂y∂P)dxdy于是格林公式成立,对于一般的由逐段光滑的闭曲线围成单连通区域,可以加若干光滑曲线将其划分为有限个由逐段光滑的闭曲线围成的单连通区域之并,每个小区域即是
X型区域,又是
Y型区域,具体操作超出数学分析课程的范畴。从下图可以看出,如果两个区域有公共边界,则公共边界曲线段上方向恰好相反,正负相抵消,最后只留下
D的边界曲线上的曲线积分。
最后,再由二重积分的区域可加性,可以证得一般单连通区域上的格林公式,对于多连通区域,也可以加若干条光滑曲线划分为若干个单连通区域,所加的曲线在应用格林公式时方向相反,正负相消,最后只留下边界的曲线。注意,一般而言,正向是逆时针的,但是如果是洞的边界则不一定,如上图,按照运动过程中
D始终在左手边,内部的光滑曲线正向应当取顺时针方向。我们之前曾经给出参数方程形式下求图形面积的公式
S=21∫xdy−ydx我们现在从格林公式的角度再推导这个公式,若
D是由有限条逐段光滑的若当闭曲线围成的区域,那么
Q(x,y)=2xP(x,y)=−2y则
Qx′−Py′=1由格林公式
21∫Lxdy−ydx=∬Ddxdy=S(D)这就是我们把面积公式写成这种形式的原因。在求解曲线积分时,如果使用常规方法,则对于若干段光滑曲线围成的闭曲线,则需要分段计算再相加,但是如果使用格林公式,可以化为二重积分,某些情况下可以简化计算。
例17.12 计算第二型曲线积分
∮L2xydx+y2dy,其中
L是由两条连接点
(0,0),(4,2)的曲线
y=2x与
y=x
组成的封闭曲线
解:
P(x,y)=2xy,Q(x,y)=y2,则
Qx′−Py′=−2x则由格林公式
∮L2xydx+y2dy=−2∫04dx∫2xx
xdy=−1564
例17.13 计算第二型曲线积分
∮L(x2+4xy)dx+(2x2+3y)dy,中
L为椭圆周
16x2+9y2=1
解:
P(x,y)=x2+4xy,Q(x,y)=2x2+3y,则
Qx′−Py′=4x−4x=0因此,由格林公式
∮L(x2+4xy)dx+(2x2+3y)dy=0
实际上,例17.13中无论取任何分段光滑的闭曲线,其曲线积分都是0。在这种条件下,积分和路径是无关的,只要确定一条
A到
B的逐段光滑的若当曲线
L1,再补充一条
B到
A的逐段光滑的若当曲线
L2,构成一条逐段光滑的若当闭曲线,如果
∫L1Pdx+Qdy+∫L2Pdx+Qdy=0则
−∫L2Pdx+Qdy=∫L1Pdx+Qdy
L2的反方向即是从
A到
B的一条逐段光滑的若当曲线,则两条路径的曲线积分是相同的,即所谓的积分与路径无关。就可以选择一条方便计算的路径来计算曲线积分,这是格林公式的小技巧之一。
例17.14 计算第二型曲线积分
∫L(exsiny−x−y)dx+(excosy−x)dy,其中
L是曲线
y=sinx从
(0,0)到
(π,0)的部分
解:
P(x,y)=exsiny−x−y,Q(x,y)=excosy−x,则
Qx′−Py′=0,设
L2为
(π,0)到
(0,0)的线段,则由格林公式
+∫L(exsiny−x−y)dx+(excosy−x)dy∫L2(exsiny−x−y)dx+(excosy−x)dy=0而
=∫L2(exsiny−x−y)dx+(excosy−x)dy∫π0−xdx=2π2因此
∫L(exsiny−x−y)dx+(excosy−x)dy=−2π2
即使
Qx′−Py′=0,应用格林公式改变路径某些情况下也可以简化计算。
例17.15 计算第二型曲线积分
∫L(2x2y−y2cosx)dx+(1−2ysinx+3x2y2)dy,其中
L为抛物线
x=2πy2从点
(0,0)到
(2π,1)的部分
解:
P(x,y)=2x2y−y2cosx,Q(x,y)=1−2sinx+2x2y2,则
Qx′−Py′=6xy2−2x2令
L1为
(0,0)到
(2π,0)之间的线段,
L2为
(2π,0)到
(2π,1)之间的线段,则
∫L1Pdx+Qdy=∫02π0dx=0∫L2Pdx+Qdy=∫011−2y+43π2y2dy=4π2由格林公式
4π2−∫LPdx+Qdy=∫01dy∫2πy22π6xy2−2x2dx=7π2−14π3因此
∫LPdx+Qdy=28π2(2π+3)
例17.16(利用格林公式转换积分曲线) 计算第二型曲线积分
∮Lx2+y2(ax−by)dx+(bx+ay)dy,其中
L为椭圆周
c2x2+d2y2=1
解:
如果用参数方程法直接计算,分母让我们的被积函数变得十分复杂,此时,我们可以试探性地用格林公式:
P(x,y)=x2+y2ax−by,Q(x,y)=x2+y2bx+ay,当
x2+y2=0时,有
∂x∂Q=(x2+y2)2by2−bx2−2axy∂y∂P=(x2+y2)2by2−bx2−2axy很巧合的是
∂x∂Q−∂y∂P=0对任意的
δ>0,令曲线
Lδ:x2+y2=δ2,方向呈逆时针,则
∮Lδx2+y2(ax−by)dx+(bx+ay)dy=2πb由格林公式,取
δ足够小,使得
Lδ完全包含在椭圆
c2x2+d2y2≤1内,则
∮LPdx+Qdy−∮LδPdx+Qdy=0故
∮LPdx+Qdy=2πb
例17.17 (逆向使用格林公式)
f(x,y)二阶连续可微,满足
∂x2∂2f+∂y2∂2f=e−(x2+y2),计算积分
∬x2+y2≤1(x∂x∂f+y∂y∂f)dxdy
解:
令
P(x,y)=−(x2+y2−1)∂y∂f,Q(x,y)=(x2+y2−1)∂x∂f,则
∂x∂Q=∂x2∂2f(x2+y2−1)+2x∂x∂f∂y∂P=−∂y2∂2f(x2+y2−1)−2y∂y∂f由格林公式有
+∮x2+y2=1Pdx+Qdy=∬x2+y2≤1(x2+y2−1)(∂x2∂2f+∂y2∂2f)dxdy2∬x2+y2≤1(x∂x∂f+y∂y∂f)dxdy则
===∬x2+y2≤1(x∂x∂f+y∂y∂f)dxdy−21∬x2+y2≤1(x2+y2−1)(∂x2∂2f+∂y2∂2f)dxdy−21∬x2+y2≤1(x2+y2−1)e−(x2+y2)dxdy−21∫02πdθ∫01r(r2−1)e−r2dr=2eπ
高斯公式
定理17.9
V⊆R3是由有限块光滑双侧曲面所围成的有界闭区域,
P,Q,R在包含
V的一个区域上连续可微,则
∭V(∂x∂P+∂y∂Q+∂z∂R)dxdydz=∬∂VPdydz+Qdzdx+Rdxdy其中曲面积分的方向向外
高斯公式的证明和格林公式是类似的,同样先证明高斯公式在一些特殊的区域上成立,再用这些区域逼近一般的区域。我们这里仅作简单的推导:
设
V={(x,y,z):(x,y)∈D,φ(x,y)≤z≤ψ(x,y)},其中
D是由逐段光滑的简单闭曲线围成的
R2上的闭区域,
φ,ψ连续可微,我们称这类型的区域为
Z型区域。下面我们证明,在
Z型区域上,成立
∭V∂z∂Rdxdydz=∬∂VRdxdy首先
=∭V∂z∂Rdxdydz=∬Ddxdy∫φ(x,y)ψ(x,y)∂z∂Rdz∬DR(x,y,ψ(x,y))dxdy−∬DR(x,y,φ(x,y))dxdy
V由一个下底面
S1,一个上底面
S2以及以下侧面
S3构成,设
D的边界曲线
L是光滑的,设
L:x=x(t),y=y(t),t∈[a,b]。则
S1:⎩⎪⎨⎪⎧x=xy=yz=φ(x,y)(x,y)∈D则法向量为
(φx′,φy′,−1),则
∬S1Rdxdy=−∬DR(x,y,φ(x,y))dxdy同理就有
∬S2Rdxdy=∬DR(x,y,ψ(x,y)dxdy侧面的参数方程可表为
⎩⎪⎨⎪⎧x=x(t)y=y(t)z=z其中
a≤t≤b,φ(x(t),y(t))≤z≤ψ(x(t),y(t)),法向量为
(x′(t),−y′(t),0),很显然
∬S3Rdxdy=0因此
∬∂VRdxdy=∬DR(x,y,ψ(x,y)dxdy−∬DR(x,y,φ(x,y)dxdy同样地可以定义
Y型区域和
X型区域,如果一个区域同时是
X,Y,Z型区域,高斯公式成立,对于一般区域,划分为
XYZ型区域的乘积即可,证明的细节超出数学分析的范畴,这里就不给出详细证明了。
例17.18 (利用高斯公式将第二型曲面积分化为三重积分)
∬Sx2dydz+y2dzdx+z2dxdy,
S为锥面
x2+y2=z2(0≤z≤h),下侧
解:
设
V为圆锥
x2+y2≤z2(0≤z≤h),设
S1为
z=h,x2+y2≤h2,上侧,则由高斯公式
=∬S1x2dydz+y2dzdx+z2dxdy+∬Sx2dydz+y2dzdx+z2dxdy2∭V(x+y+z)dxdydz求解三重积分
∭V(x+y+z)dxdydz,作变换
⎩⎪⎨⎪⎧x=rcosθy=rsinθz=u0≤u≤h,0≤θ≤2π,0≤r≤u
∣det(J)∣=r,有
=∭V(x+y+z)dxdydz∫0hdu∫0udr∫02π(r2(cosθ+sinθ)+ru)dθ=4h4π再求解
∬S1x2dydz+y2dzdx+z2dxdy=h2∬x2+y2≤h2dxdy=h4π因此
∬Sx2dydz+y2dzdx+z2dxdy=2h4π−h4π=−2h4π
对由逐片光滑的曲面围成的立体
V,由高斯公式,有
∬∂Vxdydz+ydzdx+zdxdy=3∭Vdxdydz因此,该立体的体积为
31∬∂Vxdydz+ydzdx+zdxdy,就能通过曲面积分求解立体的体积。
例17.19 求椭球
a2x2+b2y2+c2z2≤1的体积
解:曲面的参数方程为
⎩⎪⎨⎪⎧x=acosφcosθy=bcosφsinθz=csinφ−2π≤φ≤2π,0≤θ≤2π则外法向量为
n=(bccos2φcosθ,accos2φsinθ,absinφcosφ)
V=3abc∫−2π2πdφ∫02πcosφdθ=34abcπ
例17.20 计算第二型曲面积分
∬Sxdydz+ydzdx+zdxdy,其中
S为
z=4−(x2+y2)与平面
z=0所围立体的外侧
解:
由该曲面积分的几何意义,就有
∬Sxdydz+ydzdx+zdxdy=3∬D4−(x2+y2)dxdy其中
D={(x,y):x2+y2≤4},则
∬Sxdydz+ydzdx+zdxdy=3∫02dr∫02πr(4−r2)dθ=24π
例17.21 计算第二型曲面积分
∬Sxdydz+ydzdx+zdxdy,其中
S是上半球面
z=a2−x2−y2
的上侧
解:设
S2为
z=0,x2+y2≤a,下侧,则
∬Sxdydz+ydzdx+zdxdy+∬S2xdydz+ydzdx+zdxdy=2a3π而
∬S2xdydz+ydzdx+zdxdy=0因此
∬Sxdydz+ydzdx+zdxdy=2a3π
斯托克斯公式
定理17.10 设
S⊂R3是光滑双侧曲面,其边界由有限段逐段光滑曲线组成,给定
S的一侧,边界曲线取正向,
P,Q,R在包含
S的某个区域上连续可微,则
=∮∂SPdx+Qdy+Rdz∬S(∂y∂R−∂z∂Q)dydz+(∂z∂P−∂x∂R)dzdx+(∂x∂Q−∂y∂P)dxdy
斯托克斯公式可以视为格林公式在
R3上的推广,这里,我们不作详细证明,只给出一种特殊情况下的推导,假设
S:x(u,v),y(u,v),z(u,v),(u,v)∈D是二阶光滑曲面,有界闭区域
D的边界曲线是光滑的(逐段光滑情形的证明也是类似的),设为
u=u(t),v=v(t),t∈[a,b]。则
S的边界曲线为
⎩⎪⎨⎪⎧x=x(u(t),v(t))y=y(u(t),v(t))z=z(u(t),v(t))于是
∮LPdx=∫abP(xu′u′+xv′v′)dt=∮L′Pxu′du+xv′dv其中
L′为
D的边界曲线,由格林公式
==∮L′Pxu′du+Pxv′dv∬DPy′∂(u,v)∂(y,z)dudv+∬DPz′∂(u,v)∂(z,x)dudv−∬SPy′dxdy+∬SPz′dzdx同理可推得
∮LQdy=∬SQx′dxdy−∬SQz′dydz∮LRdz=∬SRy′dydz−∬SRx′dzdx这样,在这种特殊的曲面上就证得斯托克斯公式,实际上,斯托克斯公式不需要二阶光滑这种严苛的条件也能成立,这里就不提供详细的证明细节了。
应用斯托克斯公式的难点在于确定曲面
S,见下例:
例17.22 计算第二型曲线积分
∮Lydx+zdy+xdz,其中
L是
x2+y2+z2=a2,
x+y+z=0,从
x轴正向看去是逆时针方向
解:联立两个方程
{x2+y2+z2=a2x+y+z=0消去
x,所截曲线在
Oyz平面的投影为
(y+z)2+y2+z2=a2即
23(y+z)2+21(y−z)2=a2则平面取该截线所围成的在平面
x+y+z=0的曲面,取左侧,由斯托克斯公式
=∮Lydx+zdy+xdz=−3∬Sdydz+dzdx+dxdy−3∬23(y+z)2+21(y−z)2≤a2dydz作正交变换
{u=2
y+zv=2
y−z有
−3∬23(y+z)2+21(y−z)2≤a2dydz=−3∬3u2+v2≤a2dudv=−3
a2π
积分与路径无关
所谓积分与路径无关,即曲线积分只与起止点有关,而与积分路径无关。其物理背景是重力做功,重力做功只与高度的变化有关,而与下落的路径无关,那么,满足什么条件下曲线积分与路径无关呢?
定理17.11(积分与路径无关) 设
D是
R2上的一个区域,
P,Q在
D上连续可微,则以下三个命题等价:
(1)
D内任意逐段光滑的闭曲线
L,都有
∮LPdx+Qdy=0(2)
D内任意逐段光滑曲线
L的曲线积分
∫LPdx+Qdy只与
L的起止点有关,与
L的路径无关
(3)积分式
Pdx+Qdy是某函数
u(x,y)的全微分
证:
(1)→(2)是显然的
(2)→(3):
对任意的
(x0,y0)∈D,令
u(x,y)=∫(x0,y0)(x,y)Pdx+Qdy,存在
δ>0,当
∣∣(x,y)−(x0,y0)∣∣<δ时,有
(x,y)∈D,则当
Δx<δ时
u(x0+Δx,y0)−u(x0,y0)=∫x0x0+ΔxP(t,y0)dt由积分中值定理,存在
ξ介于
x0和
x0+Δx之间,使得
u(x0+Δx,y0)−u(x0,y0)=P(ξ,y0)Δx故
Δx→0limΔxu(x0+Δx,y0)−u(x0,y0)=P(x0,y0)同理可证
Δy→0limΔxu(x0,y0+Δx)−u(x0,y0)=Q(x0,y0)而
P,Q是连续的,因此,
u可微,并且
du=Pdx+Qdy
(3)→(1):假设存在
D内的连续可微函数
U,满足
∂x∂U=P,∂y∂U=Q,设
L是连接
A到
B的光滑曲线,
A∈D,B∈D,设参数方程为
x=x(t),y=y(t),t∈[a,b],其中
A=(x(a),y(a)),B=(x(b),y(b)),于是
===∫LPdx+Qdy=∫abP(x(t),y(t))x′(t)+Q(x(t),y(t))y′(t)dt∫abUx′(x(t),y(t))x′(t)+Uy′(x(t),y(t))y′(t)dt∫abdtdU(x(t),y(t))dtU(x(b),y(b))−U(x(a),y(a))因此,如果
L是逐段光滑的闭曲线,作分划
A=A0,A1,⋯,An=A,使得
L在
Ak−1Ak段是光滑曲线,记
Ak−1Ak段为
L(Ak−1Ak),则
=∫LPdx+Qdy=k=1∑n∫L(Ak−1Ak)Pdx+Qdyk=1∑n(U(Ak)−U(Ak−1))=U(An)−U(A0)=0
同样地方法可以证明:
定理17.12
V是
R3上区域,
P,Q,R在
V上连续可微,则以下命题等价:
(1)
V内任意逐段光滑的闭曲线
L,都有
∮LPdx+Qdy+Rdz=0(2)
V内任意逐段光滑曲线
L的曲线积分
∫LPdx+Qdy+Rdz只与
L的起止点有关,与
L的路径无关
(3)积分式
Pdx+Qdy+Rdz是某函数
u(x,y,z)的全微分
如果在区域
D内积分与路径无关,由定理17.11,
Pdx+Qdy在
D内是某个函数
U的全微分。则
∂y∂P=∂x∂y∂2U∂x∂Q=∂y∂x∂2U由于
P,Q都是连续可微的,则
∂x∂y∂2U=∂y∂x∂2U于是
∂y∂P=∂x∂Q这让我们想起格林公式,如果在
D内作一条逐段光滑的简单闭曲线
L,如果
L围成的有界区域
D0⊆D内,由格林公式
∫LPdx+Qdy=∫D0(∂x∂Q−∂y∂P)dxdy=0当然,定理17.11中 (1)要求的是任意的闭曲线,而非简单闭曲线,其次,应用格林公式需要
D0⊆D,因此,我们首先需要假设
D是单连通区域,那么,下一个问题是,如果在单连通区域
D内,都有
∂y∂P=∂x∂Q能否推出
D内积分与路径无关呢?答案是肯定的。
定理17.13
D⊆R2是平面上的单连通区域,并且
P,Q在
D内连续可微,则
D内积分与路径无关的充要条件是
∂y∂P=∂x∂Q
我们作一个简单的分析,如果逐段可微的闭曲线
L是简单闭曲线,那么毫无疑问,由格林公式,有
∫LPdx+Qdy=0但
L不一定是简单闭曲线,即
L可能存在自交点,如下图所示
上图所示的闭曲线就不是简单闭曲线,但是
L可以分为4段简单闭曲线,每段简单闭曲线都是逐段光滑的,也有
∫LPdx+Qdy=0。我们给自交点下个定义,对于逐段光滑的闭曲线
{x=x(t)y=y(t)t∈[a,b]很显然,即满足
a<t0<b,存在
a≤t0′<t0,有
x(t0)=x(t0′),y(t0)=y(t0′),就称
(x(t0),y(t0))是
L的一个自交点。如果
L只有有限个自交点,对应的参数从小到大排列为
a<t1<⋯<tn<b,对于
t1,存在
t1′∈[a,t1),满足
(x(t1),y(t1))=(x(t1′),y(t1′))。据此,我们可以分出一段闭曲线
L1:(x(t),y(t)),t1′≤t≤t1。
L去掉
L1段,再拼接起来,就形成一段新的闭曲线,新的闭曲线的自交点一定比原来的曲线自交点少,由数学归纳法,如果逐段可微的闭曲线
L只有有限个自交点,则
∫LPdx+Qdy=0实际上,对任意的闭曲线,我们可以以闭折线替换之,而闭折线一定满足
∫LPdx+Qdy=0,这就证得了定理17.13。问题是如何用闭折线取代一条闭曲线呢?我们引入一种区域——星形区域,对区域
G,如果存在
A∈G,对任意的
M∈G,直线段
AM⊆G,则称
G是关于
A的星形区域。我们下面证明,在星形区域上,如果满足
∂y∂P=∂x∂Q则积分与路径无关。假设
G是关于
A的星形区域,并且在
G上满足:
∂y∂P=∂x∂Q。定义
G上的函数
U(M)=∫AMPdx+Qdy,其中
AM是直线段
AM,从
A到
M。设
M(x0,y0),存在
δ>0,使得
B(M,δ)⊆G,当
Δx<δ时,取
M′(x0+Δx,y0),则
MM′的直线段都在
G内闭折线
AMM′A围成的三角形区域都在
G内,从而由格林公式
∫MM′Pdx+Qdy+U(M)−U(M′)=0而
∫MM′Pdx+Qdy=Δx∫01P(x0+tΔx,y0)dt故由积分中值定理,存在
ξ∈[0,1],有
ΔxU(M′)−U(M)=P(x0+ξΔx,y0)令
Δx→0,可知
∂x∂U=P,同理可证
∂y∂U=Q,因此,在
G上积分与路径无关。显然任何邻域都是星形区域。假设
L是单连通区域
D上的任意一条逐段光滑的闭曲线,则
Dc为闭集,
L∩D=∅,则
L与
Dc有一个正距离
δ0>0。取
L的一个分划:
Δ:A=A0,A1,⋯,An=A0,使得
λ(Δ)=1≤i≤nmaxdiamL(Ai−1Ai)<2δ0,则对于点
Ai−1,其邻域
B(Ai−1,2δ0)⊆G,并且直线段
Ai−1Ai及弧段
Ai−1Ai都在这个邻域内,而
B(Ai−1,2δ0)显然是星形区域,从而
∫L(Ai−1Ai)Pdx+Qdy=∫Ai−1AiPdx+Qdy因此,弧段
L(Ai−1Ai)可用直线段
Ai−1Ai取代,从而曲线段
L可由闭折线
L′:AA1⋯An−1A取代,故
∫LPdx+Qdy=∫L′Pdx+Qdy=0这就证明了在单连通区域上,积分与路径无关。定理17.13成立。对于
R3中积分与路径无关的条件,我们同样想到应用斯托克斯公式,所谓空间单连通区域,即
G内任意闭曲面围成的区域内部都在
G内;为了应用斯托克斯公式,我们需要引入曲面单连通区域的概念,即
G内任意逐段光滑的简单闭曲线都是
G内某片分片光滑的闭曲面的边界曲线,再应用斯托克斯公式,就有
定理17.14
G⊆R3是曲面单连通区域,并且
P,Q,R在
G上连续可微,则曲面积分
∫LPdx+Qdy+Rdz在
G内积分路径无关的充要条件为
∂z∂Q=∂y∂R∂x∂R=∂z∂P∂x∂Q=∂y∂P
该定理的证明同样要引入星形区域的概念,与
R2情形是类似,这里就不赘述了。
定理17.15
G⊆R3是星形区域,并且
P,Q,R在
G上连续可微,则曲面积分
∫LPdx+Qdy+Rdz在
G内积分路径无关的充要条件为
∂z∂Q=∂y∂R∂x∂R=∂z∂P∂x∂Q=∂y∂P