多元积分学及其应用

多元积分学及其应用

三重积分

定义

三重积分的定义通过分割、近似、求和、取极限给出:

\iiint_{\Omega} f(x,y,z) \, dV = \lim_{\lambda \to 0} \sum_{k=1}^{n} f(\xi_k, \eta_k, \zeta_k) \Delta v_k

其中,\Omega 为空间有界闭区域,\lambda 是各小闭区域直径的最大值,\Delta v_k 表示第 k 个小闭区域的体积。

性质

(三重积分的基本性质与二重积分类似,包括线性性、区域可加性、比较定理等,具体内容可参考二重积分性质拓展)

计算方法

  1. 直角坐标下的计算
    依据积分区域的特征,可选择不同的积分顺序:

    • 先一后二法(投影法):
      将空间区域 \Omega 投影到 xy 平面得平面区域 D_{xy},在 D_{xy} 内任取一点 (x,y),过该点作平行于 z 轴的直线,与 \Omega 的边界交于 z=z_1(x,y)z=z_2(x,y),则:
\iiint_{\Omega} f(x,y,z) \, dV = \iint_{D_{xy}} \left( \int_{z_1(x,y)}^{z_2(x,y)} f(x,y,z) \, dz \right) d\sigma
  • 先二后一法(截面法):
    若空间区域 \Omega 介于平面 z=c_1z=c_2 之间,过 z 轴上任一点 z \in [c_1, c_2] 作垂直于 z 轴的平面,与 \Omega 相交得平面区域 D_z,则:
\iiint_{\Omega} f(x,y,z) \, dV = \int_{c_1}^{c_2} \left( \iint_{D_z} f(x,y,z) \, dxdy \right) dz
  1. 柱坐标下的计算
    当积分区域 \Omega 具有旋转对称性(如圆柱、圆锥等)时,可采用柱坐标变换简化计算:
    • 坐标变换公式
\begin{cases} x = r \cos\theta, \\ y = r \sin\theta, \\ z = z \end{cases} \quad (0 \leq r < +\infty,\ 0 \leq \theta \leq 2\pi,\ -\infty < z < +\infty)
  • 体积元素
dV = r \, dr d\theta dz
  • 积分转换公式
\iiint_{\Omega} f(x,y,z) \, dV = \iiint_{\Omega} f(r \cos\theta, r \sin\theta, z) \cdot r \, dr d\theta dz

(注:适用场景包括积分区域由圆柱面、旋转抛物面等围成,或被积函数含 x^2 + y^2 项)

球坐标

坐标变换与体积元素

  • 变量替换公式
\begin{cases} x = r \sin\varphi \cos\theta, \\ y = r \sin\varphi \sin\theta, \\ z = r \cos\varphi \end{cases} \quad (0 \leq r < +\infty,\ 0 \leq \varphi \leq \pi,\ 0 \leq \theta \leq 2\pi)
  • 体积元素
dV = r^2 \sin\varphi \, dr d\varphi d\theta
  • 积分转换公式
\iiint_{\Omega} f(x,y,z) \, dV = \iiint_{\Omega} f(r \sin\varphi \cos\theta,\ r \sin\varphi \sin\theta,\ r \cos\varphi) \cdot r^2 \sin\varphi \, dr d\varphi d\theta

(适用场景:积分区域为球体、锥体等具有球对称性,或被积函数含 x^2+y^2+z^2 项)

积分简化技巧

利用奇偶性

  • 对称条件:若积分域 \Omega 关于 xoy 平面对称,则:
\iiint_{\Omega} f(x,y,z) \, dV = \begin{cases} 2 \iiint_{\Omega_{z \geq 0}} f(x,y,z) \, dV & \text{若 } f(x,y,-z) = f(x,y,z) \, (\text{偶函数}) \\ 0 & \text{若 } f(x,y,-z) = -f(x,y,z) \, (\text{奇函数}) \end{cases}

利用变量对称性

  • 示例:若积分域关于 x,y,z 具有轮换对称性(如球体 x^2+y^2+z^2 \leq R^2),则:
\iiint_{\Omega} (x^2 + y^2) \, dV = \frac{2}{3} \iiint_{\Omega} (x^2 + y^2 + z^2) \, dV

(通过对称性将 x^2,y^2,z^2 积分等价替换,简化计算)

对弧长的线积分(第一类线积分)

  • 定义
    设函数 f(x,y) 在曲线 L 上有定义,将 L 任意分成 n 个小弧段,各小弧段的长度为 \Delta s_i,在每个小弧段上任取一点 (\xi_i, \eta_i),作和式 \sum_{i=1}^n f(\xi_i, \eta_i) \Delta s_i。若当各小弧段长度的最大值 \lambda \to 0 时,该和式的极限存在,则称此极限为函数 f(x,y) 在曲线 L 上对弧长的线积分,记作:
\int_L f(x,y) \, ds = \lim_{\lambda \to 0} \sum_{i=1}^n f(\xi_i, \eta_i) \Delta s_i
  • 性质
    对弧长的线积分与积分路径的方向无关,即:
\int_{L(AB)} f(x,y) \, ds = \int_{L(BA)} f(x,y) \, ds
  • 计算方法

    • 直接法
      根据曲线 L 的不同方程形式,将线积分转化为定积分计算:
      1. 若曲线 L 的参数方程为 \begin{cases} x = x(t) \\ y = y(t) \end{cases}\alpha \leq t \leq \beta,且 x(t), y(t) 具有连续导数),则:
\int_L f(x,y) \, ds = \int_{\alpha}^{\beta} f(x(t), y(t)) \sqrt{[x'(t)]^2 + [y'(t)]^2} \, dt
2. 若曲线 $L$ 的方程为 $y = y(x)$($a \leq x \leq b$,且 $y(x)$ 具有连续导数),则:
\int_L f(x,y) \, ds = \int_a^b f(x, y(x)) \sqrt{1 + [y'(x)]^2} \, dx
  1. 若曲线 L 的极坐标方程为 \rho = \rho(\theta)\alpha \leq \theta \leq \beta,且 \rho(\theta) 具有连续导数),则通过 x = \rho(\theta) \cos\theta, y = \rho(\theta) \sin\theta 转化为参数方程,得:
\int_L f(x,y) \, ds = \int_{\alpha}^{\beta} f(\rho(\theta)\cos\theta, \rho(\theta)\sin\theta) \sqrt{[\rho(\theta)]^2 + [\rho'(\theta)]^2} \, d\theta
  • 利用奇偶性简化计算
    1. 若积分曲线 L 关于 y-轴对称:
      • f(-x, y) = f(x, y)(即 f 关于 x 为偶函数),则 \int_L f(x,y) \, ds = 2\int_{L_{x \geq 0}} f(x,y) \, ds
      • f(-x, y) = -f(x, y)(即 f 关于 x 为奇函数),则 \int_L f(x,y) \, ds = 0
    2. 若积分曲线 L 关于 x-轴对称:
      • f(x, -y) = f(x, y)(即 f 关于 y 为偶函数),则 \int_L f(x,y) \, ds = 2\int_{L_{y \geq 0}} f(x,y) \, ds
      • f(x, -y) = -f(x, y)(即 f 关于 y 为奇函数),则 \int_L f(x,y) \, ds = 0

对弧长的线积分(第一类线积分)

  • 利用对称性
    若积分曲线 C 关于直线 y = x 对称,则:
\int_C f(x,y) \, ds = \int_C f(y,x) \, ds

特别地,当 f(x,y) = f(x)(仅依赖于 x)时:

\int_C f(x) \, ds = \int_C f(y) \, ds
  • 空间曲线的计算
    设空间曲线 L 的参数方程为 \begin{cases} x = x(t) \\ y = y(t) \\ z = z(t) \end{cases}\alpha \leq t \leq \beta,且 x(t), y(t), z(t) 具有连续导数),则:
\int_L f(x,y,z) \, ds = \int_{\alpha}^{\beta} f(x(t), y(t), z(t)) \sqrt{[x'(t)]^2 + [y'(t)]^2 + [z'(t)]^2} \, dt

对坐标的线积分(第二类线积分)

  • 定义
    设函数 P(x,y), Q(x,y) 在有向曲线 L 上有定义,将 L 任意分成 n 个有向小弧段,各小弧段在 x-轴、y-轴上的投影分别为 \Delta x_i, \Delta y_i,在每个小弧段上任取一点 (\xi_i, \eta_i)。若当各小弧段长度的最大值 \lambda \to 0 时,和式 \sum_{i=1}^n [P(\xi_i, \eta_i) \Delta x_i + Q(\xi_i, \eta_i) \Delta y_i] 的极限存在,则称此极限为向量函数 (P, Q) 在有向曲线 L 上对坐标的线积分,记作:
\int_L P(x,y) \, dx + Q(x,y) \, dy = \lim_{\lambda \to 0} \sum_{i=1}^n [P(\xi_i, \eta_i) \Delta x_i + Q(\xi_i, \eta_i) \Delta y_i]
  • 性质
    对坐标的线积分与积分路径方向有关,即:
\int_{L(AB)} P \, dx + Q \, dy = -\int_{L(BA)} P \, dx + Q \, dy

(其中 L(AB) 表示从 AB 的路径,L(BA) 为反向路径)

  • 计算方法(平面曲线)

    • 直接法
      若有向曲线 L 的参数方程为 \begin{cases} x = x(t) \\ y = y(t) \end{cases}(当参数 t\alpha 单调变到 \beta 时,对应路径从起点到终点),则:
\int_L P \, dx + Q \, dy = \int_{\alpha}^{\beta} \left[ P(x(t), y(t)) x'(t) + Q(x(t), y(t)) y'(t) \right] dt
  • 格林公式
    设闭区域 D 由分段光滑的有向曲线 L(取正向,即沿 L 行走时 D 始终在左侧)围成,函数 P(x,y), Q(x,y)D 上具有一阶连续偏导数,则:
\int_L P \, dx + Q \, dy = \iint_D \left( \frac{\partial Q}{\partial x} - \frac{\partial P}{\partial y} \right) d\sigma

(注:若 L 非闭,可通过“补线”使其闭合后应用格林公式,再减去补线部分的积分)

  • 利用线积分与路径无关
    1. 判定条件:在单连通区域 D 内,若 \frac{\partial P}{\partial y} = \frac{\partial Q}{\partial x},则线积分 \int_L P \, dx + Q \, dy 与路径无关,仅依赖于起点和终点。
    2. 计算方法
      • 改换路径:选择简便路径(如折线)计算积分;
      • 原函数法:若存在函数 F(x,y) 使得 dF = P \, dx + Q \, dy(即 FP \, dx + Q \, dy 的原函数),则:
\int_{(x_1,y_1)}^{(x_2,y_2)} P \, dx + Q \, dy = F(x_2,y_2) - F(x_1,y_1)
     (求原函数的常用方法:偏积分法、凑微分法)
  • 两类线积分的联系
    平面上对坐标的线积分可转化为对弧长的线积分:
\int_L P \, dx + Q \, dy = \int_L (P \cos\alpha + Q \cos\beta) \, ds

其中 \cos\alpha, \cos\beta 为有向曲线 L 上点 (x,y) 处切向量的方向余弦。

  • 计算方法(空间曲线)

    • 直接法
      设空间有向曲线 L 的参数方程为 \begin{cases} x = x(t) \\ y = y(t) \\ z = z(t) \end{cases}(参数 t\alpha 变到 \beta 对应路径方向),则:
\int_L P \, dx + Q \, dy + R \, dz = \int_{\alpha}^{\beta} \left[ P x'(t) + Q y'(t) + R z'(t) \right] dt
(其中 $P = P(x(t), y(t), z(t))$, $Q, R$ 类似)
  • 斯托克斯公式
    \varSigma 是以空间有向闭曲线 L 为边界的分片光滑有向曲面(L 的正向与 \varSigma 的侧符合右手规则),函数 P, Q, R 在包含 \varSigma 的空间区域内具有一阶连续偏导数,则:
\int_L P \, dx + Q \, dy + R \, dz = \iint_{\varSigma} \left( \frac{\partial R}{\partial y} - \frac{\partial Q}{\partial z} \right) dy dz + \left( \frac{\partial P}{\partial z} - \frac{\partial R}{\partial x} \right) dz dx + \left( \frac{\partial Q}{\partial x} - \frac{\partial P}{\partial y} \right) dx dy

(或用行列式形式表示为:\iint_{\varSigma} \begin{vmatrix} dy dz & dz dx & dx dy \\ \frac{\partial}{\partial x} & \frac{\partial}{\partial y} & \frac{\partial}{\partial z} \\ P & Q & R \end{vmatrix}

对面积的面积分(第一类面积分)

  • 定义
    设函数 f(x,y,z) 在分片光滑曲面 \varSigma 上有定义,将 \varSigma 任意分成 n 个小曲面片,各小曲面片的面积为 \Delta S_i,在每个小曲面片上任取一点 (\xi_i, \eta_i, \zeta_i),作和式 \sum_{i=1}^n f(\xi_i, \eta_i, \zeta_i) \Delta S_i。若当各小曲面片直径的最大值 \lambda \to 0 时,该和式的极限存在,则称此极限为函数 f(x,y,z) 在曲面 \varSigma 上对面积的面积分,记作:
\iint_{\varSigma} f(x,y,z) \, dS = \lim_{\lambda \to 0} \sum_{i=1}^n f(\xi_i, \eta_i, \zeta_i) \Delta S_i
  • 性质
    对面积的面积分与积分曲面的方向(即曲面的“正侧”或“负侧”)无关,即:
\iint_{\varSigma} f(x,y,z) \, dS = \iint_{-\varSigma} f(x,y,z) \, dS

(其中 -\varSigma 表示与 \varSigma 取向相反的同一曲面)

  • 计算方法

    • 直接法
      若曲面 \varSigma 由方程 z = z(x,y) 给出,其在 xoy-平面上的投影区域为 D,且 z(x,y)D 上具有一阶连续偏导数,则:
\iint_{\varSigma} f(x,y,z) \, dS = \iint_{D} f(x,y,z(x,y)) \sqrt{1 + \left( \frac{\partial z}{\partial x} \right)^2 + \left( \frac{\partial z}{\partial y} \right)^2} \, d\sigma
(注:若曲面由 $x = x(y,z)$ 或 $y = y(x,z)$ 给出,可类似投影到 $yoz$-平面或 $xoz$-平面,替换相应的偏导数项)
  • 利用奇偶性
    若曲面 \varSigma 关于 xoy-面对称,则:
\iint_{\varSigma} f(x,y,z) \, dS = \begin{cases} 2\iint_{\varSigma_{z \geq 0}} f(x,y,z) \, dS, & \text{若 } f(x,y,-z) = f(x,y,z) \, (\text{即 } f \text{ 关于 } z \text{ 为偶函数}) \\ 0, & \text{若 } f(x,y,-z) = -f(x,y,z) \, (\text{即 } f \text{ 关于 } z \text{ 为奇函数}) \end{cases}
(类似结论可推广到关于 $yoz$-面或 $xoz$-面对称的曲面)
  • 利用对称性(轮换对称性)
    若曲面 \varSigma 具有轮换对称性(如球面 x^2 + y^2 + z^2 = 1),即对换变量 x, y, z 后曲面方程不变,则:
\iint_{\varSigma} f(x,y,z) \, dS = \iint_{\varSigma} f(y,z,x) \, dS = \iint_{\varSigma} f(z,x,y) \, dS

示例:对球面 \varSigma: x^2 + y^2 + z^2 = 1,有 \iint_{\varSigma} x^2 \, dS = \iint_{\varSigma} y^2 \, dS = \iint_{\varSigma} z^2 \, dS,因此:

\iint_{\varSigma} (x^2 + y^2) \, dS = \frac{2}{3} \iint_{\varSigma} (x^2 + y^2 + z^2) \, dS = \frac{2}{3} \iint_{\varSigma} 1 \, dS = \frac{2}{3} \cdot 4\pi = \frac{8\pi}{3}
(其中 $\iint_{\varSigma} 1 \, dS = 4\pi$ 为球面面积)

对坐标的面积分(第二类面积分)

  • 定义
    对坐标的面积分(第二类面积分)涉及向量场在有向曲面上的累积效应,典型形式为 \iint_{\varSigma} P(x,y,z) \, dy dz + Q(x,y,z) \, dz dx + R(x,y,z) \, dx dy。以 R(x,y,z) \, dx dy 分量为例,其定义为:
\iint_{\varSigma} R(x,y,z) \, dx dy = \lim_{\lambda \to 0} \sum_{i=1}^n R(\xi_i, \eta_i, \zeta_i) (\Delta S_i)_{xy}

其中 (\Delta S_i)_{xy} 表示小曲面片 \Delta S_ixoy- 平面上的投影(带符号,取决于曲面侧的取向)。

  • 性质
    对坐标的面积分与积分曲面的方向(“侧”)密切相关。若将曲面 \varSigma 的取向反转(记为 -\varSigma),则积分值变号:
\iint_{\varSigma} P \, dy dz + Q \, dz dx + R \, dx dy = -\iint_{-\varSigma} P \, dy dz + Q \, dz dx + R \, dx dy
  • 计算方法

    • 直接法(按投影方向转化为二重积分)
      根据曲面方程和投影方向选择合适的坐标面,需结合曲面取向确定积分符号(“上正下负、前正后负、右正左负”):
      1. 投影到 xoy-平面
        若曲面 \varSigma: z = z(x,y)(x,y) \in D_{xy}D_{xy} 为投影区域),取向“上侧”(法向量 z-分量为正)时:
\iint_{\varSigma} R(x,y,z) \, dx dy = \iint_{D_{xy}} R(x,y,z(x,y)) \, dx dy

取向“下侧”时,等式右侧加负号。

  1. 投影到 yoz-平面
    若曲面 \varSigma: x = x(y,z)(y,z) \in D_{yz}),取向“前侧”(法向量 x-分量为正)时:
\iint_{\varSigma} P(x,y,z) \, dy dz = \iint_{D_{yz}} P(x(y,z), y,z) \, dy dz

取向“后侧”时,等式右侧加负号。

  1. 投影到 xoz-平面
    若曲面 \varSigma: y = y(z,x)(z,x) \in D_{xz}),取向“右侧”(法向量 y-分量为正)时:
\iint_{\varSigma} Q(x,y,z) \, dz dx = \iint_{D_{xz}} Q(x,y(z,x), z) \, dz dx
   取向“左侧”时,等式右侧加负号。
  • 高斯公式(Gauss's Theorem)
    设空间闭区域 \varOmega 由分片光滑的有向闭曲面 \varSigma(取外侧,即法向量指向区域外部)围成,函数 P, Q, R\varOmega 上具有一阶连续偏导数,则:
\iint_{\varSigma_{\text{外}}} P \, dy dz + Q \, dz dx + R \, dx dy = \iiint_{\varOmega} \left( \frac{\partial P}{\partial x} + \frac{\partial Q}{\partial y} + \frac{\partial R}{\partial z} \right) dV

补面技巧:若 \varSigma 非闭,可通过添加辅助曲面 \varSigma' 使其闭合(形成 \varSigma_{\text{ 闭}} = \varSigma + \varSigma'),应用高斯公式后减去辅助曲面的积分:

\iint_{\varSigma} \cdots = \iint_{\varSigma_{\text{闭}}} \cdots - \iint_{\varSigma'} \cdots