# 变量代换

在引入积分问题之前,先思考如何改变曲面片的参数化的变量

CC^\infty 映射 f:U(Rm)V(Rn)f: U(\subset \mathbb R^m) \to V(\subset \mathbb R^n)微分同胚 (Diffeomorphism)「微分同相」,当且仅当 ff 是双射,并且 fff1f^{-1} 都是光滑映射

命题
EEsts-t 平面 R2\mathbb R^2 中的非空开集,DDuvu-v 平面 R2\mathbb R^2 中的非空开集,σ:DR3\boldsymbol \sigma: D \to \mathbb R^3 为曲面片的参数化。
Φ:ED\Phi: E \to D 是一个微分同胚,那么 σΦ:ER3\boldsymbol \sigma \circ \Phi: E \to \mathbb R^3 也是一个曲面片的参数化,并且 σ(E)=(σΦ)(E)\boldsymbol \sigma(E) = (\boldsymbol \sigma \circ \Phi)(E)

证明

τ:=σΦ\boldsymbol \tau := \boldsymbol \sigma \circ \Phi,作为复合映射,显然 τ\boldsymbol \tau 是光滑且单射的
利用复合映射的偏微分,有

τs(s,t)=usσu+vsσv,τt(s,t)=utσu+vtσv\boldsymbol \tau_s(s,t) = \frac{\partial u}{\partial s} \boldsymbol \sigma_u + \frac{\partial v}{\partial s} \boldsymbol \sigma_v,\quad \boldsymbol \tau_t(s,t) = \frac{\partial u}{\partial t} \boldsymbol \sigma_u + \frac{\partial v}{\partial t} \boldsymbol \sigma_v

计算向量积

τs(s,t)×τt(s,t)=(usσu+vsσv)×(utσu+vtσv)=(usvtutvs)σu×σv=JΦ(s,t)σu×σv\begin{aligned} \boldsymbol \tau_s(s,t) \times \boldsymbol \tau_t(s,t) &= \left(\frac{\partial u}{\partial s} \boldsymbol \sigma_u + \frac{\partial v}{\partial s} \boldsymbol \sigma_v\right) \times \left(\frac{\partial u}{\partial t} \boldsymbol \sigma_u + \frac{\partial v}{\partial t} \boldsymbol \sigma_v\right) \\ &= \left( \frac{\partial u}{\partial s} \frac{\partial v}{\partial t} - \frac{\partial u}{\partial t} \frac{\partial v}{\partial s} \right) \boldsymbol \sigma_u \times \boldsymbol \sigma_v \\ &= J_\Phi(s,t) \cdot \boldsymbol \sigma_u \times \boldsymbol \sigma_v \end{aligned}

由于 Φ\Phi 是微分同胚,所以 Jacobian 行列式 JΦ(s,t)0J_\Phi(s,t) \neq 0,又由于 σ\boldsymbol \sigma 是曲面片的参数化,所以 σu×σv0\boldsymbol \sigma_u \times \boldsymbol \sigma_v \neq \boldsymbol 0,因此 τs(s,t)×τt(s,t)0\boldsymbol \tau_s(s,t) \times \boldsymbol \tau_t(s,t) \neq \boldsymbol 0τ\boldsymbol \tau 也是一个曲面片的参数化
\square

通过这样的方法,参数化 σ\boldsymbol \sigma 变成了新的参数化 σΦ\boldsymbol \sigma \circ \Phi,称这样的操作为 变量代换 (Change of Variables)「変数変換」

  • 称该代换是保持方向的,当且仅当 JΦ(s,t)>0J_\Phi(s,t) \gt 0 对任意 (s,t)E(s,t) \in E 都成立
    • 此时 nσΦ(s,t)=nσ(u,v)\boldsymbol n_{\boldsymbol \sigma \circ \Phi}(s,t) = \boldsymbol n_{\boldsymbol \sigma}(u,v)
  • 称该代换是改变方向的,当且仅当 JΦ(s,t)<0J_\Phi(s,t) \lt 0 对任意 (s,t)E(s,t) \in E 都成立
    • 此时 nσΦ(s,t)=nσ(u,v)\boldsymbol n_{\boldsymbol \sigma \circ \Phi}(s,t) = -\boldsymbol n_{\boldsymbol \sigma}(u,v)

示例
一般地,由正则矩阵 AA 给出的线性映射

Φ:R2R2,Φ(s,t)=A(st)\Phi: \mathbb R^2 \to \mathbb R^2,\quad \Phi(s,t) = A \begin{pmatrix}s \\ t\end{pmatrix}

是一个微分同胚,且 JΦ(s,t)=detAJ_\Phi(s,t) = \det A 是常数

  • detA>0\det A \gt 0 时,Φ\Phi 是保持方向的变量代换
  • detA<0\det A \lt 0 时,Φ\Phi 是改变方向的变量代换

# 曲面积分

类似于在正则曲面之前引入了曲面片,曲面片作为曲面的基本单位,具备参数化的能力,这使得其法向量场可以由参数化的偏微分稳定给出

SR3S \subset \mathbb R^3 为以 σ:DR3\boldsymbol \sigma: D \to \mathbb R^3 参数化的曲面片
我们来分析一下如何计算曲面片 SS 上的面积:
考虑位于点 σ(u,v)\boldsymbol \sigma(u,v) 处的一个小片段
曲面积近似示例

考虑偏微分向量积的模长

σu(u,v)×σv(u,v)=limu,v0σ(u+u,v)σ(u,v)u×σ(u,v+v)σ(u,v)v=limu,v0σ(u+u,v)σ(u,v)u×σ(u,v+v)σ(u,v)v=limu,v0(σ(u+u,v)σ(u,v))×(σ(u,v+v)σ(u,v))uv\begin{aligned} \|\boldsymbol \sigma_u(u,v) \times \boldsymbol \sigma_v(u,v)\| &= \left\| \lim_{\triangle u, \triangle v \to 0} \frac{\boldsymbol \sigma(u + \triangle u, v) - \boldsymbol \sigma(u,v)}{\triangle u} \times \frac{\boldsymbol \sigma(u, v + \triangle v) - \boldsymbol \sigma(u,v)}{\triangle v} \right\| \\ &= \lim_{\triangle u, \triangle v \to 0} \left\| \frac{\boldsymbol \sigma(u + \triangle u, v) - \boldsymbol \sigma(u,v)}{\triangle u} \times \frac{\boldsymbol \sigma(u, v + \triangle v) - \boldsymbol \sigma(u,v)}{\triangle v} \right\| \\ &= \lim_{\triangle u, \triangle v \to 0} \frac{\|(\boldsymbol \sigma(u + \triangle u, v) - \boldsymbol \sigma(u,v)) \times (\boldsymbol \sigma(u, v + \triangle v) - \boldsymbol \sigma(u,v))\|}{|\triangle u| \cdot |\triangle v|} \end{aligned}

因此,当 u\triangle uv\triangle v 足够小时,可以近似

(σ(u+u,v)σ(u,v))×(σ(u,v+v)σ(u,v))σu(u,v)×σv(u,v)uv\|(\boldsymbol \sigma(u + \triangle u, v) - \boldsymbol \sigma(u,v)) \times (\boldsymbol \sigma(u, v + \triangle v) - \boldsymbol \sigma(u,v))\| \approx \|\boldsymbol \sigma_u(u,v) \times \boldsymbol \sigma_v(u,v)\| \cdot |\triangle u| \cdot |\triangle v|

也就是说,令矩形 R:=[u,u+u]×[v,v+v]R := [u, u + \triangle u] \times [v, v + \triangle v],则 σ(R)\boldsymbol \sigma(R) 的面积近似为这里的曲面面积。
那么统合所有小矩形的面积近似,就可以得到曲面片 SS 的面积近似

Rσu(u,v)×σv(u,v)uv\sum_R \|\boldsymbol \sigma_u(u,v) \times \boldsymbol \sigma_v(u,v)\| \cdot |\triangle u| \cdot |\triangle v|

定义
SR3S \subset \mathbb R^3 为以 σ:DR3\boldsymbol \sigma: D \to \mathbb R^3 参数化的曲面片
设曲面的一部分 TST \subset S 的原象可测,那么其 曲面积 (Surface Area)「曲面積」 定义为

Area(T)=σ1(T)σu(u,v)×σv(u,v)dudv\mathrm{Area}(T) = \iint_{\boldsymbol \sigma^{-1}(T)} \|\boldsymbol \sigma_u(u,v) \times \boldsymbol \sigma_v(u,v)\| du dv

从上述结果不难看出,曲面积分上的基本单位为 dA=σu(u,v)×σv(u,v)dudvdA = \|\boldsymbol \sigma_u(u,v) \times \boldsymbol \sigma_v(u,v)\| du dv,对于标量场的面积分,只需要给这个基本单位加上标量场的值。
但是对于向量场,类似向量场的曲线积分处理方式,需要将其与法向量场点积后再积分
一般地,曲面片上的积分按照如下给出

定义 曲面片的面积分
SR3S \subset \mathbb R^3 为以 σ:DR3\boldsymbol \sigma: D \to \mathbb R^3 参数化的曲面片,TST \subset S 的原象可测

  • 定义 TT 上连续标量场 f:TRf: T \to \mathbb R 的面积分为

TfdA=σ1(T)f(σ(u,v))σu(u,v)×σv(u,v)dudv\iint_T f dA = \iint_{\boldsymbol \sigma^{-1}(T)} f(\boldsymbol \sigma(u,v)) \|\boldsymbol \sigma_u(u,v) \times \boldsymbol \sigma_v(u,v)\| du dv

  • 定义 TT 上连续向量场 F:TR3\boldsymbol F: T \to \mathbb R^3 的面积分为

TFdA=σ1(T)F(σ(u,v))(σu(u,v)×σv(u,v))dudv\iint_T \boldsymbol F d\boldsymbol A = \iint_{\boldsymbol \sigma^{-1}(T)} \boldsymbol F(\boldsymbol \sigma(u,v)) \cdot (\boldsymbol \sigma_u(u,v) \times \boldsymbol \sigma_v(u,v)) du dv

分别称

  • dA=σu(u,v)×σv(u,v)dudvdA = \|\boldsymbol \sigma_u(u,v) \times \boldsymbol \sigma_v(u,v)\| du dv 为曲面片 SS 上的 面积元素 (Area Element)「面積要素」
  • dA=σu(u,v)×σv(u,v)dudvd\boldsymbol A = \boldsymbol \sigma_u(u,v) \times \boldsymbol \sigma_v(u,v) du dv 为曲面片 SS 上的 向量面积元素 (Vector Area Element)「ベクトル面積要素」

接下来是正则曲面的面积分问题。由于正则曲面上的点不一定都可以被一个曲面片参数化,因此无法直接套用上述定义来计算面积分。
SS 为正则曲面。我们想要考虑某一个标量场或者向量场在 TST \subset S 上的面积分

如果 TT 是一个可以被正向参数化的曲面片,那么标量场和向量场的面积微元可以分别由下式给出

dA=σu(u,v)×σv(u,v)dudv,dA=σu(u,v)×σv(u,v)dudvdA = \|\boldsymbol \sigma_u(u,v) \times \boldsymbol \sigma_v(u,v)\| du dv,\quad d\boldsymbol A = \boldsymbol \sigma_u(u,v) \times \boldsymbol \sigma_v(u,v) du dv

但是通常来说无法默认 TT 是一个曲面片。因此可以类似于曲线积分中分段光滑的概念,考虑 TT 是否可以被多个曲面片分割。

定义 正则曲面的面积分
令 S 为 R3\mathbb R^3 中的正则曲面,TST \subset S,设存在一系列子集 T1,T2,,TnST_1, T_2, \cdots, T_n \subset S 满足

  • T=i=1nTiT = \bigcup_{i=1}^n T_i
  • 每个 TiT_i 都可以被一个曲面片参数化,即存在 σi:DiR3\boldsymbol \sigma_i: D_i \to \mathbb R^3 使得 σi(Di)=Ti\boldsymbol \sigma_i(D_i) = T_i,并且 DiD_i 可测
  • iji \neq j 时,TiTjT_i \cap T_j 的面积为零

那么,定义沿曲面 TT

  • 标量场 f:TRf: T \to \mathbb R 的面积分为

TfdA=i=1nTifdA\iint_T f dA = \sum_{i=1}^n \iint_{T_i} f dA

  • 向量场 F:TR3\boldsymbol F: T \to \mathbb R^3 的面积分为

TFdA=i=1nTiFdA\iint_T \boldsymbol F d\boldsymbol A = \sum_{i=1}^n \iint_{T_i} \boldsymbol F d\boldsymbol A

这样一来,正则曲面的面积分问题转为多个曲面片的面积分问题。用求和符号可以改写为

TfdA=i=1nDif(σi(u,v))σi,u(u,v)×σi,v(u,v)dudvTFdA=i=1nDiF(σi(u,v))(σi,u(u,v)×σi,v(u,v))dudv\begin{aligned} \iint_T f dA &= \sum_{i=1}^n \iint_{D_i} f(\boldsymbol \sigma_i(u,v)) \|\boldsymbol \sigma_{i,u}(u,v) \times \boldsymbol \sigma_{i,v}(u,v)\| du dv \\ \iint_T \boldsymbol F d\boldsymbol A &= \sum_{i=1}^n \iint_{D_i} \boldsymbol F(\boldsymbol \sigma_i(u,v)) \cdot (\boldsymbol \sigma_{i,u}(u,v) \times \boldsymbol \sigma_{i,v}(u,v)) du dv \end{aligned}

# Stokes 定理

Stokes 定理一定程度上可以说是 Green 定理的曲面推广

定理 Stokes 定理
SS 为以 σ:DR3\boldsymbol \sigma: D \to \mathbb R^3 参数化的曲面片
领域 ΩD\Omega \subset D 的边界 Ω\partial \Omega 是由分段光滑的正则单纯闭曲线 C0,C1,,CC_0, C_1, \cdots, C_\ell 组成的,并且满足 ΩD\overline \Omega \subset D
那么经由 σ\boldsymbol \sigma 的映射,各个像 σ(Ci)\boldsymbol \sigma(C_i) 也是分段光滑的正则单纯闭曲线,且在 SS 中包围了一个区域 σ(Ω)\boldsymbol \sigma(\overline \Omega)
此时,对于定义在 Uσ(Ω)U \supset \boldsymbol \sigma(\overline \Omega) 上的连续向量场 F:UR3\boldsymbol F: U \to \mathbb R^3,有

σ(Ω)curlFdA=i=0σ(Ci)Fdx\iint_{\boldsymbol \sigma(\overline \Omega)} \mathrm{curl} \boldsymbol F \cdot d\boldsymbol A = \sum_{i=0}^\ell \int_{\boldsymbol \sigma(C_i)} \boldsymbol F \cdot d\boldsymbol x

证明

Stokes 定理的正式证明将在微分形式的章节给出,此处暂时省略

\square

作为一个特例,Stokes 在闭曲面上的结论广为人知

定义
称正则曲面 SS封闭曲面 (Closed Surface)「閉曲面」,当且仅当 SS 是有界的,且是闭集

在开始给出闭曲面上的 Stokes 定理之前,让我们接触三角剖分

定义 三角形
SS 为正则曲面,σ:DR3\boldsymbol \sigma:D \to \mathbb R^3SS 的局部参数化
ΩR2\Omega \subset \mathbb R^2 满足 ΩD\overline \Omega \subset D,称像 σ(Ω)\boldsymbol \sigma(\Omega)SS 上的一个 三角形 (Triangle)「三角形」,当且仅当 Ω\Omega 以三个正则曲线连接所得的单纯闭曲线作为边界的区域