第一型曲面积分&第二型曲面积分

第一型曲面积分

定义1(第一型曲面积分)

SR3S\subset \mathbb R^3 为光滑曲面,f:SRf:S\rightarrow \mathbb R,设 r:[a,b]×[c,d]S\vec{r}:[a,b]\times[c,d]\rightarrow SSS 的参数方程,设

π:a=s0<s1<<sN1=bc=t0<t1<<tN2=d\begin{aligned} \pi : a=&s_0<s_1<\cdots<s_{N_1} = b\\ c=&t_0<t_1<\cdots<t_{N_2} = d \end{aligned}

[a,b]×[c,d][a,b]\times[c,d] 的一个分划,考虑和式

=i,jf(ξij)σ(Sij)\sum = \sum_{i,j}f(\xi_{ij})\sigma(S_{ij})

其中 Sij=r([si1,si]×[tj1,tj]), ξijSijS_{ij} = \vec{r}([s_{i-1},s_i]\times[t_{j-1},t_j]),\ \xi_{ij}\in S_{ij}
如果 limπ0\lim\limits_{\triangle\pi\rightarrow 0}\sum 存在,记

Sfdσ=limπ0\int_Sf\,d\sigma = \lim_{\triangle\pi\rightarrow 0}\sum

Sfdσ\int_Sf\,d\sigma 称为 ff 沿曲面 SS 的积分。


如果令 f(x)=1f(x) = 1,结合定理2和上节定义的 曲面面积,则 σ(S)=S1dσ\sigma(S) = \int_S 1\, d\sigma

定理2(第一型曲面积分的计算公式)

SR3S\subset\mathbb R^3 为光滑曲面,r:[a,b]×[c,d]S\vec{r}:[a,b]\times[c,d]\rightarrow SSS 的参数方程,若 f:SRf:S\rightarrow \mathbb R 连续,则

Sfdσ=[a,b]×[c,d]fr rs×rtdsdt\int_Sf\,d\sigma = \int_{[a,b]\times[c,d]}f\circ \vec{r}\ |\vec{r_s}\times\vec{r_t}|\,ds\,dt


证明方法和 第一型曲线积分的证明方法 类似(去估计两者的差值,建议先搞懂曲线积分的证明方法),但会更加复杂。

证明: (其中 \sum 和分划 π\pi 的定义与 定义1 中一致)

[a,b]×[c,d]f(r(s,t)) rs×rtdsdti=1N1j=1N2f(r(si,tj))σ(Sij)i=1N1j=1N2si1sitj1tjf(r(s,t))rs×rtdsdti=1N1j=1N2f(r(si,tj))si1sitj1tjrs×rtdsdtsi1sitj1tjf(r(s,t))rs×rtdsdt=i=1N1j=1N2si1sitj1tj(f(r(si,tj))f(r(s,t)))rs×rtdsdti=1N1j=1N2si1sitj1tjf(r(si,tj))f(r(s,t))rs×rtdsdt σ(S)ω(Δπ)0\begin{aligned} &\left|\sum - \int_{[a,b]\times[c,d]}f(\vec{r}(s,t))\ |\vec{r_s}\times\vec{r_t}|\,ds\,dt\right|\\ \leqslant &\left|\sum_{i=1}^{N_1}\sum_{j=1}^{N_2}f(\vec{r}(s_{i},t_{j}))\sigma(S_{ij})- \sum_{i=1}^{N_1}\sum_{j=1}^{N_2}\int_{s_{i-1}}^{s_i}\int_{t_{j-1}}^{t_j}f(\vec{r}(s,t))|\vec{r}_s\times\vec{r}_t|\,ds\,dt\right|\\ \leqslant &\sum_{i=1}^{N_1}\sum_{j=1}^{N_2}\left|f(\vec{r}(s_{i},t_{j}))\int_{s_{i-1}}^{s_i}\int_{t_{j-1}}^{t_j}|\vec{r}_s\times\vec{r}_t|\,ds\,dt-\int_{s_{i-1}}^{s_i}\int_{t_{j-1}}^{t_j}f(\vec{r}(s,t))|\vec{r}_s\times\vec{r}_t|\,ds\,dt\right|\\ = &\sum_{i=1}^{N_1}\sum_{j=1}^{N_2}\left|\int_{s_{i-1}}^{s_i}\int_{t_{j-1}}^{t_j}(f(\vec{r}(s_{i},t_{j}))-f(\vec{r}(s,t)))|\vec{r}_s\times\vec{r}_t|\,ds\,dt\right|\\ \leqslant &\sum_{i=1}^{N_1}\sum_{j=1}^{N_2}\int_{s_{i-1}}^{s_i}\int_{t_{j-1}}^{t_j}|f(\vec{r}(s_{i},t_{j}))-f(\vec{r}(s,t))|\cdot|\vec{r}_s\times\vec{r}_t|\,ds\,dt\\ \leqslant &\ \sigma(S)\cdot \omega(\Delta \pi)\rightarrow 0 \end{aligned}

其中,ω(Δπ)=sup(u,v)(s,t)Δπf(r(u,v))f(r(s,t))\omega(\Delta\pi) = \sup\limits_{|(u,v)-(s,t)|\leqslant\Delta\pi}|f(\vec{r}(u,v))-f(\vec{r}(s,t))|,当 Δπ0\Delta\pi\rightarrow 0 时,ω(Δπ)0\omega(\Delta\pi)\rightarrow 0

QED

定义3(非闭方体的积分域)

SR3S\subset \mathbb R^3 为光滑曲线,r:DS\vec{r}:\overline{D}\rightarrow SSS 的参数方程,如果 f:SRf:S\rightarrow \mathbb R 连续,定义

Sfdσ=Dfr rs×rtdsdt\int_Sf\,d\sigma = \int_Df\circ \vec{r}\ |\vec{r}_s\times\vec{r}_t|\,ds\,dt


要确保这里的定义是良定义,还需要证明对于不同的参数方程,该积分值都相同(然而我不会证😢)

定理4(第一型曲面积分的性质)

SR3S\subset \mathbb R^3 为光滑曲面,f,g:SRf, g:S\rightarrow \mathbb R 连续。

  1. 如果 f0f\geqslant 0,则 Sfdσ0\int_Sf\,d\sigma \geqslant 0,“==”成立当且仅当 f0f\equiv 0

  2. 如果 k,lRk, l\in \mathbb R,则 S{kf+lg}dσ=kSfdσ+lSdσ\int_S\{kf+lg\}\,d\sigma = k\int_Sf\,d\sigma+l\int_S\,d\sigma

  3. 如果 fgf \leqslant g,则 SfdσSgdσ\int_Sf\,d\sigma\leqslant \int_Sg\,d\sigma

  4. SfdσSfdσ|\int_Sf\,d\sigma|\leqslant\int_S|f|\,d\sigma


这几个性质都可以根据定义直接得出。

BR(x0)={xR3:xx0R}B_R(x_0) = \{x\in \mathbb R^3:|x-x_0|\leqslant R\},特别的 BR={xR3:xR}B_R = \{x\in\mathbb R^3:|x|\leqslant R\}

定理5(球面上的积分)

fC(BR)f\in C(\partial B_R),则

BRfdσ=02π0πf(Rsinφcosθ,Rsinφsinθ,Rcosφ)R2sinφdφdθ\int_{\partial B_R}f\,d\sigma = \int_0^{2\pi}\int_0^{\pi}f(R\sin\varphi\cos\theta,R\sin\varphi\sin\theta,R\cos\varphi)R^2\sin\varphi\,d\varphi\,d\theta


思路: BR\partial B_R 的参数方程为 r=(Rsinφcosθ,Rsinφsinθ,Rcosφ)\vec{r} = (R\sin\varphi\cos\theta,R\sin\varphi\sin\theta,R\cos\varphi),且 rθ×rφ=R2sinφ|\vec{r}_{\theta}\times\vec{r}_{\varphi}| = R^2\sin\varphi,所以代入 计算公式 即可得出结论。

推论1

fC(BR)f\in C(\partial B_{R}),则

BRf(x)dσ=R2B1f(Rx)dσ\int_{\partial B_R} f(x)\,d\sigma = R^2\int_{\partial B_1}f(Rx)\,d\sigma

注: 这里的缩放只能是对球进行,而不能对椭球缩放。(因为球具有很好的对称性?)

推论2

fC(BR)f\in C(\overline{B}_R)

BRf=0R{Brfdσ}dr\int_{B_R}f = \int_0^R\left\{\int_{\partial B_r}f\,d\sigma\right\}\,dr

推论3(球体表面积公式)

f=1f = 1,则

σ(BR)=BR1dσ=02π0πR2sinφdφdθ=4πR2\sigma(\partial B_R) = \int_{\partial B_R}1\,d\sigma = \int_0^{2\pi}\int_0^{\pi}R^2\sin\varphi\,d\varphi\,d\theta=4\pi R^2

定义6(分片曲面积分)

SR3, S=i=1NSiS\subset \mathbb R^3,\ S = \bigcup\limits_{i=1}^NS_i,其中 Si (i=1N)S_i \ (i=1\sim N) 为内部互不相交的光滑曲面,f:SRf:S\rightarrow \mathbb R 连续,定义

Sfdσ=i=1NSifdσ\int_Sf\,d\sigma = \sum_{i=1}^N\int_{S_i}f\,d\sigma


注意现在定义的第一型曲面积分和曲线积分都是和坐标系无关的定义,所以可以在不同位置建立坐标系,相当于对原来的点做了平移变换,正交变换,比如

BR(x0)f(x)dσ=BRf(x+x0)dσ\int_{\partial B_R(x_0)}f(x)\,d\sigma = \int_{\partial B_R}f(x+x_0)\,d\sigma

PoissonPoisson 公式:

B1f(ax+by+cz)dσ=2π11f(ua2+b2+c2)du\int_{ \partial B_1}f(ax+by+cz)\,d\sigma = 2\pi\int_{-1}^1f(u\sqrt{a^2+b^2+c^2})\,du

例题

S={(x,y,z):(xa)2+(yb)2+(zc)2=R2 且 zc}S = \{(x, y, z): (x-a)^2+(y-b)^2+(z-c)^2 = R^2\text{ 且 }z\geqslant c\},求 S(x+y+z)dσ\int_S(x+y+z)\,d\sigma

解:(BR)+={(x,y,z):x2+y2+z2=R 且 z0}(\partial B_R)^+ = \{(x, y, z):x^2+y^2+z^2 = R\text{ 且 }z\geqslant 0\}

S(x+y+z)dσ=(BR)+(x+y+z+a+b+c)dσ= 2πR2(a+b+c)+(BR)+(x+y+z)dσ=故可以消去由于x,y在曲面上的对称性 2πR2(a+b+c)+(BR)+zdσ= 2πR2(a+b+c)+12BRzdσ= 2πR2(a+b+c)+1202π0πRcosφR2sinφdφdθ= 2πR2(a+b+c)+πR320π2cosφsinφdφ= 2πR2(a+b+c)+πR3\begin{aligned} \int_S(x+y+z)\,d\sigma =& \int_{(\partial B_R)^+}(x+y+z+a+b+c)\,d\sigma\\ = &\ 2\pi R^2(a+b+c)+\int_{(\partial B_R)^+}(x+y+z)\,d\sigma\\ \xlongequal[\text{故可以消去}]{\text{由于}x,y\text{在曲面上的对称性}}&\ 2\pi R^2(a+b+c)+\int_{(\partial B_R)^+}z\,d\sigma\\ = &\ 2\pi R^2(a+b+c)+\frac{1}{2}\int_{\partial B_R}|z|\,d\sigma\\ = &\ 2\pi R^2(a+b+c)+\frac{1}{2}\int_0^{2\pi}\int_0^{\pi}|R\cos\varphi|R^2\sin\varphi\,d\varphi\,d\theta\\ = &\ 2\pi R^2(a+b+c)+\pi R^3\cdot2\int_0^{\frac{\pi}{2}}\cos\varphi\sin\varphi\,d\varphi\\ = &\ 2\pi R^2(a+b+c)+\pi R^3 \end{aligned}

第二型曲面积分

SR3S\subset \mathbb R^3 是一个光滑曲面,r:DS\vec{r}:\overline{D}\rightarrow SSS 的参数方程。

定义1(曲面的定向)

SS 上的一个连续单位法向量场,称为 SS 的一个定向

命题2(定向只有两个)

SS 有且仅有两个定向,它们分别为

n=rs×rtrs×rt(r)1n\vec{n} = \frac{\vec{r}_s\times\vec{r}_t}{|\vec{r_s}\times\vec{r}_t|}\circ(\vec{r})^{-1}\quad\text{和}\quad -\vec{n}

定义3(定向曲面)

n\vec{n}SS 的一个定向{S,n}\{S,\vec{n}\} 为定向曲面,简记为 SS,称 SS 为定向曲面,n\vec{n}SS 的正向,S={S,n}-S=\{S,-\vec{n}\}

定义4(第二型曲面积分)

SR3S\subset \mathbb R^3 为定向曲面,n\vec{n}SS 的正向,F:SR3\vec{F}:S\rightarrow \mathbb R^3,设 r:[a,b]×[c,d]\vec{r}:[a, b]\times[c,d]SS 的参数方程,设

π:a=s0<s1<<sN1=bc=t0<t1<<tN2=d\begin{aligned} \pi : a =& s_0<s_1<\cdots<s_{N_1} = b\\ c =& t_0<t_1<\cdots<t_{N_2} = d \end{aligned}

[a,b]×[c,d][a, b]\times[c, d] 的分划,考虑和式

=i,j(Fn)(ξij)σ(Sij)\sum = \sum_{i,j}(\vec{F}\cdot\vec{n})(\xi_{ij})\sigma(S_{ij})

其中 Sij=r([si1,si]×[tj1,tj]), ξijSijS_{ij} = \vec{r}([s_{i-1},s_i]\times[t_{j-1},t_j]),\ \xi_{ij}\in S_{ij}
如果 limπ0\lim\limits_{\triangle\pi\rightarrow 0}\sum 存在,记

SFdσ=limπ0\int_S\vec{F}\cdot d\vec{\sigma} = \lim_{\triangle\pi\rightarrow 0}\sum

SFdσ\int_S\vec{F}\cdot d\vec{\sigma} 称为向量场 F\vec{F} 沿曲面 SS 的积分。


不难发现,第二型曲面积分的定义式和 第一型曲面积分 的定义式区别只有 ffFn\vec{F}\cdot \vec{n} 这里,于是有如下定理。

定理5(第二型和第一型曲面积分的转化)

SR3S\subset \mathbb R^3 为定向曲面,r:[a,b]×[c,d]S\vec{r}:[a,b]\times [c,d]\rightarrow SSS 的参数方程,n\vec{n}SS 的正向,设 F:SR3\vec{F}:S\rightarrow \mathbb R^3 连续,则

SFdσ=SFndσ\int_S\vec{F}\cdot d\vec{\sigma} = \int_S\vec{F}\cdot\vec{n}\,d\sigma

定义6(非闭方体的积分域)

SR3S\subset \mathbb R^3 为定向曲面,n\vec{n}SS 的正向,设 F:SR3\vec{F}:S\rightarrow \mathbb R^3 连续,定义

SFdσ=SFndσ\int_S\vec{F}\cdot d\vec{\sigma} = \int_S\vec{F}\cdot\vec{n}\,d\sigma

定理7(第二型曲面积分的性质)

SR3S\subset \mathbb R^3 为定向曲面,F,G:SR3\vec{F},\vec{G}:S\rightarrow \mathbb R^3 连续,

  1. k,lRk, l\in \mathbb R,则 S(kF+lG)dσ=kSFdσ+lSGdσ\int_S(k\vec{F}+l\vec{G})\cdot d\vec{\sigma} = k\int_S\vec{F}\cdot d\vec{\sigma}+l\int_S\vec{G}\cdot d\vec{\sigma}

  2. SFdσ=SFdσ\int_S \vec{F}\cdot d\vec{\sigma} = -\int_{-S}\vec{F}\cdot d\vec{\sigma}

  3. r:DS\vec{r}:\overline{D}\rightarrow SSS 的参数方程,且 n=rs×rtrs×rt(r)1\vec{n} = \dfrac{\vec{r}_s\times\vec{r}_t}{|\vec{r_s}\times\vec{r}_t|}\circ(\vec{r})^{-1},则

SFdσ=D(Fr)(rs×rt)dsdt\int_S\vec{F}\,d\vec{\sigma} = \int_D(\vec{F}\circ\vec{r})\cdot(\vec{r}_s\times\vec{r}_t)\,ds\,dt


证明下第三条,利用 定义6 即可,

SFdσ=SFnds=D(Fn)rrs×rtdsdt=D(Fr)(nr) rs×rtdsdt=代入n的表达式D(Fr)(rs×rt)dsdt\begin{aligned} \int_S\vec{F}\,d\vec{\sigma} =& \int_S\vec{F}\cdot\vec{n}\,ds = \int_D(\vec{F}\cdot\vec{n})\circ\vec{r}\,|\vec{r}_s\times\vec{r}_t|\,ds\,dt \\ =& \int_D(\vec{F}\circ\vec{r})\cdot(\vec{n}\circ\vec{r})\ |\vec{r}_s\times\vec{r}_t|\,ds\,dt\\ \xlongequal{\text{代入}\vec{n}\text{的表达式}}&\int_D(\vec{F}\circ\vec{r})\cdot(\vec{r}_s\times\vec{r}_t)\,ds\,dt \end{aligned}

定义8(分片曲面积分)

S=i=1NSiS=\bigcup\limits_{i=1}^NS_iSi (i=1N)S_i\ (i=1\sim N) 为内部互不相交的定向曲面,设 F:SR3\vec{F}:S\rightarrow \mathbb R^3 连续,定义

SFdσ=i=1NSiFdσ\int_S\vec{F}\cdot d\vec{\sigma} = \sum_{i=1}^N\int_{S_i}\vec{F}\cdot\,d\vec{\sigma}


则有 dσ=ndσ=(rs×rt)dsdtd\vec{\sigma} = \vec{n}\,d\sigma = (\vec{r}_s\times\vec{r}_t)\,ds\,dt

在具体坐标系下第二型曲面积分的计算公式

建立空间直角坐标系 OxyzOxyz,则 n=(cosα,cosβ,cosγ)\vec{n} = (\cos\alpha, \cos\beta, \cos\gamma)

dσ=ndσ=(cosαdσ,cosβdσ,cosγdσ)=(dydz,dzdx,dxdy)d\vec{\sigma} = \vec{n}d\sigma = (\cos\alpha\,d\sigma,\cos\beta\,d\sigma,\cos\gamma\,d\sigma) \xlongequal{\text{令}} (dy\,dz,dz\,dx,dx\,dy)。(可以理解为曲面在法向量方向上的投影)

F=(F1(x,y,z),F2(x,y,z),F3(x,y,z))\vec{F} = (F_1(x, y, z),F_2(x, y, z), F_3(x, y, z)),则

SFdσ=SF1(x,y,z)dydz+F2(x,y,z)dzdx+F3(x,y,z)dxdy\int_S\vec{F}\cdot d\vec{\sigma} = \int_SF_1(x, y, z)\,dy\,dz+F_2(x, y, z)\,dz\,dx+F_3(x, y, z)\,dx\,dy

r=(x(s,t),y(s,t),z(s,t))\vec{r} = (x(s, t), y(s, t), z(s, t)),曲面的正向:n=rs×rtrs×rt(r)1\vec{n} = \dfrac{\vec{r}_s\times\vec{r}_t}{|\vec{r}_s\times\vec{r}_t|}\circ (\vec{r})^{-1},则

(Fr)(rs×rt)= [Fr,rs,rt](混合积)=F1rF2rF3rxsyszsxtytzt= (F1r)(y,z)(s,t)+(F2r)(z,x)(s,t)+(F3r)(x,y)(s,t)\begin{aligned} (\vec{F}\circ \vec{r})\cdot (\vec{r}_s\times\vec{r}_t) =&\ [\vec{F}\circ \vec{r},\vec{r}_s,\vec{r}_t]&(\text{混合积})\\ =& \left | \begin{matrix} F_1\circ\vec{r}&F_2\circ\vec{r}&F_3\circ\vec{r}\\ x_s&y_s&z_s\\ x_t&y_t&z_t \end{matrix}\right |\\ =&\ (F_1\circ\vec{r})\frac{\partial(y,z)}{\partial(s,t)}+(F_2\circ\vec{r})\frac{\partial(z,x)}{\partial(s,t)}+(F_3\circ\vec{r})\frac{\partial(x,y)}{\partial(s,t)}\\ \end{aligned}

其中,(y,z)(s,t)=ysytzszt=yszsytzt\dfrac{\partial (y,z)}{\partial (s,t)} = \left|\begin{matrix}y_s&y_t\\z_s&z_t\end{matrix}\right| = \left|\begin{matrix}y_s&z_s\\y_t&z_t\end{matrix}\right|JacobianJacobian 行列式。

则(由 定理7 计算公式展开)

SFdσ=SF1dydz+F2dzdx+F3dxdy=D((F1r)(y,z)(s,t)+(F2r)(z,x)(s,t)+(F3r)(x,y)(s,t))dsdt\begin{aligned} \int_S\vec{F}\cdot d\vec{\sigma} =& \int_SF_1\,dy\,dz + F_2\,dz\,dx + F_3\,dx\,dy\\ =& \int_D\left((F_1\circ\vec{r})\frac{\partial(y,z)}{\partial(s,t)}+(F_2\circ\vec{r})\frac{\partial(z,x)}{\partial(s,t)}+(F_3\circ\vec{r})\frac{\partial(x,y)}{\partial(s,t)}\right)\,ds\,dt \end{aligned}

所以有如下对应关系

SRdxdy=D(Rr)(x,y)(s,t)dsdtSQdzdx=D(Qr)(z,x)(s,t)dsdtSPdydz=D(Pr)(y,z)(s,t)dsdt\begin{aligned} \int_SR\,dx\,dy =& \int_D(R\circ\vec{r})\frac{\partial(x,y)}{\partial(s,t)}\,ds\,dt\\ \int_SQ\,dz\,dx =& \int_D(Q\circ\vec{r})\frac{\partial(z,x)}{\partial(s,t)}\,ds\,dt\\ \int_SP\,dy\,dz =& \int_D(P\circ\vec{r})\frac{\partial(y,z)}{\partial(s,t)}\,ds\,dt \end{aligned}


第一型曲面积分&第二型曲面积分
https://wty-yy.github.io/posts/29775/
作者
wty
发布于
2021年11月14日
许可协议