12374 words
62 minutes
工科数学分析-第六章

第六章 多元函数积分学及其应用#

第一节 多元数量值函数积分的概念与性质#

一、 核心概念与统一数学定义#

1. 多元数量值函数积分的统一思想#

多元积分的建立均采用“分(分割)、匀(近似)、合(求和)、精(取极限)”的步骤。它不考虑方向性(即与方向无关),称为第一类多元积分

2. 统一数学定义(定义 1.1)#

(Ω)(\Omega) 是一个有界的、可度量的(可求长/宽/高/面积/体积等)几何形体,f(M)f(M) 是定义在 (Ω)(\Omega) 上的有界数量值函数。

  1. :将 (Ω)(\Omega) 任意划分为 nn 个小部分 (ΔΩk)(\Delta \Omega_k),其度量记为 ΔΩk\Delta \Omega_k

  2. :在每个 (ΔΩk)(\Delta \Omega_k) 上任取一点 MkM_k,作乘积 f(Mk)ΔΩkf(M_k) \Delta \Omega_k

  3. :作和式 k=1nf(Mk)ΔΩk\sum_{k=1}^n f(M_k) \Delta \Omega_k

  4. :设 dd 为这 nn 个子域直径的最大值。若当 d0d \to 0 时,该和式的极限存在,且与分割方式及点 MkM_k 的取法无关,则称 f(M)f(M)(Ω)(\Omega)可积,其极限值称为 f(M)f(M)(Ω)(\Omega) 上的多元数量值积分,记作:

   (Ω)f(M)dΩ=limd0k=1nf(Mk)ΔΩk\int_{(\Omega)} f(M) \mathrm{d}\Omega = \lim_{d \to 0} \sum_{k=1}^n f(M_k) \Delta \Omega_k


二、 五大积分形式对照表#

根据积分域 (Ω)(\Omega) 的维度和几何特征不同,该统一积分对应以下五种具体形式:

| 积分域 (Ω)(\Omega) 的类型 | 几何形体度量 ΔΩk\Delta\Omega_k | 积分名称与记号 | 积分微元 |

| :--- | :--- | :--- | :--- |

| 区间 [a,b][a,b] | 子区间长度 Δxk\Delta x_k | 定积分abf(x)dx\int_{a}^{b} f(x) \mathrm{d}x | 长度微元 dx\mathrm{d}x |

| 平面区域 (σ)(\sigma) | 子区域面积 Δσk\Delta \sigma_k | 二重积分(σ)f(x,y)dσ\iint_{(\sigma)} f(x,y) \mathrm{d}\sigma | 面积微元 dσ\mathrm{d}\sigma (或 dxdy\mathrm{d}x\mathrm{d}y) |

| 空间区域 (V)(V) | 子区域体积 ΔVk\Delta V_k | 三重积分(V)f(x,y,z)dV\iiint_{(V)} f(x,y,z) \mathrm{d}V | 体积微元 dV\mathrm{d}V (或 dxdydz\mathrm{d}x\mathrm{d}y\mathrm{d}z) |

| 曲线弧段 (C)(C) | 子弧段长度 Δsk\Delta s_k | 第一类曲线积分(C)f(x,y)ds\int_{(C)} f(x,y) \mathrm{d}s | 弧长微元 ds\mathrm{d}s |

| 曲面片 (S)(S) | 子曲面面积 ΔSk\Delta S_k | 第一类曲面积分(S)f(x,y,z)dS\iint_{(S)} f(x,y,z) \mathrm{d}S | 面积微元 dS\mathrm{d}S |


三、 积分的重要性质与定理#

f(M),g(M)f(M), g(M) 在可度量的闭区域 (Ω)(\Omega) 上连续,则积分一定存在,且满足以下性质:

1. 线性性质#

  • (Ω)kf(M)dΩ=k(Ω)f(M)dΩ(k 为常数)\int_{(\Omega)} k f(M) \mathrm{d}\Omega = k \int_{(\Omega)} f(M) \mathrm{d}\Omega \quad (k\text{ 为常数})

  • (Ω)[f(M)±g(M)]dΩ=(Ω)f(M)dΩ±(Ω)g(M)dΩ\int_{(\Omega)} [f(M) \pm g(M)] \mathrm{d}\Omega = \int_{(\Omega)} f(M) \mathrm{d}\Omega \pm \int_{(\Omega)} g(M) \mathrm{d}\Omega

2. 积分域的可加性#

(Ω)=(Ω1)(Ω2)(\Omega) = (\Omega_1) \cup (\Omega_2),且 (Ω1)(\Omega_1)(Ω2)(\Omega_2) 除边界外无公共部分,则:

(Ω)f(M)dΩ=(Ω1)f(M)dΩ+(Ω2)f(M)dΩ\int_{(\Omega)} f(M) \mathrm{d}\Omega = \int_{(\Omega_1)} f(M) \mathrm{d}\Omega + \int_{(\Omega_2)} f(M) \mathrm{d}\Omega

3. 积分不等式#

  • 保号/比较性:若 f(M)g(M)f(M) \le g(M),则 (Ω)f(M)dΩ(Ω)g(M)dΩ\int_{(\Omega)} f(M) \mathrm{d}\Omega \le \int_{(\Omega)} g(M) \mathrm{d}\Omega

  • 绝对值不等式(Ω)f(M)dΩ(Ω)f(M)dΩ\left| \int_{(\Omega)} f(M) \mathrm{d}\Omega \right| \le \int_{(\Omega)} |f(M)| \mathrm{d}\Omega

  • 估值定理:若 lf(M)Ll \le f(M) \le L,则 lΩ(Ω)f(M)dΩLΩl \cdot \Omega \le \int_{(\Omega)} f(M) \mathrm{d}\Omega \le L \cdot \Omega (其中 Ω\Omega 为几何形体的总度量)。

4. 积分中值定理#

f(M)f(M)有界闭连通集 (Ω)(\Omega) 上连续,则在 (Ω)(\Omega) 上至少存在一点 PP,使得:

(Ω)f(M)dΩ=f(P)Ω\int_{(\Omega)} f(M) \mathrm{d}\Omega = f(P) \cdot \Omega


四、 核心重难点解析#

1. 为什么极限条件必须是“直径 d0d \to 0”而非“测度 ΔΩk0\Delta \Omega_k \to 0”?(习题 A.2)#

  • 原理解析:几何形体的度量(如面积、体积)趋于 00,并不意味着其形状在所有方向上都收缩。例如,在二维平面中,一个长为 11、宽为 ϵ\epsilon 的极扁矩形,当 ϵ0\epsilon \to 0 时其面积趋于 00,但它在长度方向上并没有收缩成点。

  • 如果不要求直径 d0d \to 0(即每个方向上的尺寸都趋于 00),就无法保证每个子区域内的函数值可以用单一点 f(Mk)f(M_k) 的函数值来近似,从而导致积分和式的极限可能不唯一或不存在。因此,必须要求所有子区域的最大直径 d0d \to 0

2. 被积函数 f(M)f(M) 的有界性要求#

定义 1.1 的批注中指出,被积函数 ff 必须在积分域上有界。若 ff 无界,其 Riemann 积分的和式极限通常不存在(此时需要引入广义积分)。


五、 典型课后题(习题 6.1)精选解析#

1. 物理意义分析(习题 A.1)#

  • 题目:当 f(M)=1f(M)=1 时,积分 (Ω)f(M)dΩ\int_{(\Omega)} f(M) \mathrm{d}\Omega 表示什么意义?

  • 解答:根据定义,此时积分式为 (Ω)1dΩ=limd0ΔΩk=Ω\int_{(\Omega)} 1 \mathrm{d}\Omega = \lim_{d \to 0}\sum \Delta\Omega_k = \Omega。它表示积分域 (Ω)(\Omega) 的度量。具体地:

  * 在定积分中,表示区间长度 bab-a

  * 在二重积分中,表示平面区域的面积;

  * 在三重积分中,表示空间区域的体积;

  * 在第一类线/面积分中,分别表示曲线的弧长或曲面的面积。

2. 不用计算比较积分大小(习题 A.5 选讲)#

  • 题目 (1):比较 (σ)(x+y)dσ\iint_{(\sigma)} (x+y) \mathrm{d}\sigma(σ)(x+y)2dσ\iint_{(\sigma)} (x+y)^2 \mathrm{d}\sigma,其中 (σ)={(x,y)x0,y0,x+y1}(\sigma) = \{ (x,y) \mid x \ge 0, y \ge 0, x+y \le 1 \}

  * 解析:在区域 (σ)(\sigma) 内,由于 x0,y0x \ge 0, y \ge 0x+y1x+y \le 1,所以有 0x+y10 \le x+y \le 1

  * 在区间 [0,1][0,1] 上,显然有 (x+y)2x+y(x+y)^2 \le x+y

  * 根据积分的比较性质,可得:

    (σ)(x+y)dσ(σ)(x+y)2dσ\iint_{(\sigma)} (x+y) \mathrm{d}\sigma \ge \iint_{(\sigma)} (x+y)^2 \mathrm{d}\sigma

  • 题目 (4):比较 (σ1)x2ydσ\iint_{(\sigma_1)} x^2y \mathrm{d}\sigma(σ2)x2ydσ\iint_{(\sigma_2)} x^2y \mathrm{d}\sigma,其中 (σ1)(\sigma_1) 为上半单位圆,(σ2)(\sigma_2) 为整个单位圆。

  * 解析:根据积分域的可加性,整个单位圆积分可拆为上半圆和下半圆:

    (σ2)x2ydσ=(σ1)x2ydσ+(下半圆)x2ydσ\iint_{(\sigma_2)} x^2y \mathrm{d}\sigma = \iint_{(\sigma_1)} x^2y \mathrm{d}\sigma + \iint_{(\text{下半圆})} x^2y \mathrm{d}\sigma

    在上半圆 (σ1)(\sigma_1) 中, y0y \ge 0,所以被积函数 x2y0x^2y \ge 0,积分值大于 00

    在下半圆中, y<0y < 0,所以被积函数 x2y0x^2y \le 0,积分值小于 00

    根据对称性,上下半圆的积分大小相等符号相反,即整个单位圆 (σ2)(\sigma_2) 上的积分为 00

    因此:

    (σ1)x2ydσ>(σ2)x2ydσ\iint_{(\sigma_1)} x^2y \mathrm{d}\sigma > \iint_{(\sigma_2)} x^2y \mathrm{d}\sigma

3. 积分极限计算(习题 B.4 难点)#

  • 题目:设 f(x,y)f(x,y) 为连续函数,求 limr0+1πr2(σr)f(x,y)dσ\lim_{r \to 0^+} \frac{1}{\pi r^2} \iint_{(\sigma_r)} f(x,y) \mathrm{d}\sigma,其中 (σr)={(x,y)(xx0)2+(yy0)2r2}(\sigma_r) = \{ (x,y) \mid (x-x_0)^2 + (y-y_0)^2 \le r^2 \}

  • 解析

  1. 积分域 (σr)(\sigma_r) 是以 (x0,y0)(x_0, y_0) 为圆心, rr 为半径的闭圆盘,其面积为 πr2\pi r^2

  2. 由于 f(x,y)f(x,y) 在有界闭区域 (σr)(\sigma_r) 上连续,由二重积分中值定理,必存在圆盘内一点 (ξr,ηr)(σr)(\xi_r, \eta_r) \in (\sigma_r),使得:

     (σr)f(x,y)dσ=f(ξr,ηr)面积=f(ξr,ηr)πr2\iint_{(\sigma_r)} f(x,y) \mathrm{d}\sigma = f(\xi_r, \eta_r) \cdot \text{面积} = f(\xi_r, \eta_r) \pi r^2

  3. 代入极限式中:

     limr0+1πr2[f(ξr,ηr)πr2]=limr0+f(ξr,ηr)\lim_{r \to 0^+} \frac{1}{\pi r^2} \left[ f(\xi_r, \eta_r) \pi r^2 \right] = \lim_{r \to 0^+} f(\xi_r, \eta_r)

  4. 当 r0+r \to 0^+ 时,圆盘收缩至圆心,因此夹逼得出 (ξr,ηr)(x0,y0)(\xi_r, \eta_r) \to (x_0, y_0)。由于 f(x,y)f(x,y) 连续,有:

     limr0+f(ξr,ηr)=f(x0,y0)\lim_{r \to 0^+} f(\xi_r, \eta_r) = f(x_0, y_0)

  * 结论:该极限值为 f(x0,y0)f(x_0, y_0)。这一结论在物理和场论中常用于定义局域物理量的密度。


第二节 二重积分的计算#

一、 核心理论与计算公式#

1. 二重积分的几何意义#

若被积函数 f(x,y)0f(x,y) \ge 0,则二重积分 (σ)f(x,y)dσ\iint_{(\sigma)} f(x,y) \mathrm{d}\sigma 的值等于以平面闭区域 (σ)(\sigma) 为底,以曲面 z=f(x,y)z=f(x,y) 为顶的曲顶柱体的体积

  • 想一想拓展

  1. 若 f(x,y)0f(x,y) \le 0,则积分值等于该曲顶柱体体积的相反数(即 V-V)。

  2. 若 f(x,y)f(x,y)(σ)(\sigma) 上变号,则由曲面与 xOyxOy 面所围成立体的体积为 V=(σ)f(x,y)dσV = \iint_{(\sigma)} |f(x,y)| \mathrm{d}\sigma

2. 直角坐标系下的计算(累次积分法)#

根据积分区域 (σ)(\sigma) 的几何特征,主要分为以下两类:

  • XX 型区域(适合先对 yy 积分,再对 xx 积分):

  (σ)={(x,y)axb,y1(x)yy2(x)}(\sigma) = \{ (x,y) \mid a \le x \le b, \, y_1(x) \le y \le y_2(x) \}

  (σ)f(x,y)dσ=abdxy1(x)y2(x)f(x,y)dy\iint_{(\sigma)} f(x,y) \mathrm{d}\sigma = \int_{a}^{b} \mathrm{d}x \int_{y_1(x)}^{y_2(x)} f(x,y) \mathrm{d}y

  • YY 型区域(适合先对 xx 积分,再对 yy 积分):

  (σ)={(x,y)cyd,x1(y)xx2(y)}(\sigma) = \{ (x,y) \mid c \le y \le d, \, x_1(y) \le x \le x_2(y) \}

  (σ)f(x,y)dσ=cddyx1(y)x2(y)f(x,y)dx\iint_{(\sigma)} f(x,y) \mathrm{d}\sigma = \int_{c}^{d} \mathrm{d}y \int_{x_1(y)}^{x_2(y)} f(x,y) \mathrm{d}x

  • Fubini 定理(交换积分顺序)

  若 f(x,y)f(x,y) 在区域 (σ)(\sigma) 上连续,且 (σ)(\sigma) 既是 XX 型又是 YY 型区域,则两种累次积分结果相等:

  abdxy1(x)y2(x)f(x,y)dy=cddyx1(y)x2(y)f(x,y)dx\int_{a}^{b} \mathrm{d}x \int_{y_1(x)}^{y_2(x)} f(x,y) \mathrm{d}y = \int_{c}^{d} \mathrm{d}y \int_{x_1(y)}^{x_2(y)} f(x,y) \mathrm{d}x


二、 极坐标系下的计算#

当积分区域为圆、圆环、扇形,或被积函数含有 x2+y2x^2+y^2 项时,采用极坐标变换极为便利。

1. 坐标变换公式#

x=ρcosθ,y=ρsinθx = \rho \cos\theta, \quad y = \rho \sin\theta

此时,直角坐标系下的面积微元与极坐标面积微元的关系为:

dσ=ρdρdθ\mathrm{d}\sigma = \rho \mathrm{d}\rho \mathrm{d}\theta

2. 累次积分公式(以 θ\theta 一步分割为例)#

若积分域 (σ)(\sigma) 在极坐标下表示为 αθβ,ρ1(θ)ρρ2(θ)\alpha \le \theta \le \beta, \, \rho_1(\theta) \le \rho \le \rho_2(\theta),则:

(σ)f(x,y)dσ=αβdθρ1(θ)ρ2(θ)f(ρcosθ,ρsinθ)ρdρ\iint_{(\sigma)} f(x,y) \mathrm{d}\sigma = \int_{\alpha}^{\beta} \mathrm{d}\theta \int_{\rho_1(\theta)}^{\rho_2(\theta)} f(\rho\cos\theta, \rho\sin\theta) \rho \mathrm{d}\rho


三、 曲线坐标系下的变换(一般换元法)#

对于边界复杂、直角或极坐标不便处理的区域,需采用一般坐标变换 T:u=u(x,y),v=v(x,y)T: u=u(x,y), v=v(x,y)

1. 雅可比行列式(Jacobian)#

设逆变换为 x=x(u,v),y=y(u,v)x=x(u,v), y=y(u,v),其雅可比行列式为:

J=(x,y)(u,v)=xuxvyuyvJ = \frac{\partial(x,y)}{\partial(u,v)} = \begin{vmatrix} x_u & x_v \\ y_u & y_v \end{vmatrix}

变换后的面积元素为:dσ=(x,y)(u,v)dudv\mathrm{d}\sigma = \left| \frac{\partial(x,y)}{\partial(u,v)} \right| \mathrm{d}u \mathrm{d}v

2. 换元公式#

(σ)f(x,y)dσ=(σ)f[x(u,v),y(u,v)](x,y)(u,v)dudv\iint_{(\sigma)} f(x,y) \mathrm{d}\sigma = \iint_{(\sigma')} f[x(u,v), y(u,v)] \left| \frac{\partial(x,y)}{\partial(u,v)} \right| \mathrm{d}u \mathrm{d}v

  • 常用特例:广义极坐标变换(适用于椭圆区域 x2a2+y2b21\frac{x^2}{a^2} + \frac{y^2}{b^2} \le 1

  x=aρcosθ,y=bρsinθx = a\rho\cos\theta, \quad y = b\rho\sin\theta

  此时,其雅可比行列式计算得面积微元为:dσ=abρdρdθ\mathrm{d}\sigma = ab\rho \mathrm{d}\rho \mathrm{d}\theta


四、 本节核心重难点解析#

1. 积分顺序的合理选择(如例 2.4)#

计算二重积分时,若积分顺序选择不当,不仅会增加计算量,甚至可能导致原函数无法用初等函数表示。

  • 例 2.4 启示:计算 (σ)sinxxdσ\iint_{(\sigma)} \frac{\sin x}{x} \mathrm{d}\sigma,若先对 xx 积分,由于 sinxxdx\int \frac{\sin x}{x} \mathrm{d}x 无法用初等函数积出,计算直接受阻;而改用先对 yy 积分,再对 xx 积分,则可顺利解出。

2. 对称性的巧用(奇偶性与轮换对称性)#

充分利用区域和函数的对称性可极大地简化计算,甚至可以直接得出积分为 00 的结论(如例 2.5、例 2.11):

  • 奇偶对称性:若积分域关于某坐标轴对称,而被积函数关于该坐标轴对应的变量为奇函数,则积分为 00

  • 轮换对称性:若积分域关于直线 y=xy=x 对称,则被积函数中的 xxyy 对调,积分值保持不变。


五、 典型习题与例题精选分析#

1. 奇偶对称性应用(课后题 3.(9))#

  • 题目:计算 (σ)(xsiny+ycosx)dσ\iint_{(\sigma)} (x\sin y + y\cos x) \mathrm{d}\sigma,其中 (σ)(\sigma) 是以 (1,1),(1,1),(1,1)(1,1), (-1,1), (-1,-1) 为顶点的三角形区域。

  • 解析

  1. 描绘积分区域 (σ)(\sigma)。该三角形的边界线为 y=1y=1x=1x=-1,以及 y=xy=x

  2. 观察可知,该区域关于直线 y=xy = -x 对称。

  3. 我们将累次积分列出(作为 XX 型区域):

     I=11dxx1(xsiny+ycosx)dyI = \int_{-1}^{1} \mathrm{d}x \int_{x}^{1} (x\sin y + y\cos x) \mathrm{d}y

     先对 yy 积分:

     x1(xsiny+ycosx)dy=[xcosy+y22cosx]y=x1=xcos1+12cosx+xcosxx22cosx\int_{x}^{1} (x\sin y + y\cos x) \mathrm{d}y = \left[ -x\cos y + \frac{y^2}{2}\cos x \right]_{y=x}^{1} = -x\cos 1 + \frac{1}{2}\cos x + x\cos x - \frac{x^2}{2}\cos x

  4. 对 xx 进行 [1,1][-1, 1] 区间上的积分:

     由于 xcos1-x\cos 1xcosxx\cos x 均为关于 xx 的奇函数,在对称区间 [1,1][-1, 1] 上的积分为 00

     剩余部分为偶函数积分:

     I=11(12cosxx22cosx)dx=01(1x2)cosxdxI = \int_{-1}^{1} \left( \frac{1}{2}\cos x - \frac{x^2}{2}\cos x \right) \mathrm{d}x = \int_{0}^{1} (1-x^2)\cos x \mathrm{d}x

     利用分部积分法可求得:

     (1x2)cosxdx=(3x2)sinx2xcosx\int (1-x^2)\cos x \mathrm{d}x = (3-x^2)\sin x - 2x\cos x

     代入上下限 0011

     I=2sin12cos1I = 2\sin 1 - 2\cos 1

2. 极坐标下二重积分计算(课后题 6.(1))#

  • 题目:计算 (σ)ex2+y2dσ\iint_{(\sigma)} e^{x^2+y^2} \mathrm{d}\sigma,其中 (σ)={(x,y)a2x2+y2b2}(\sigma) = \{ (x,y) \mid a^2 \le x^2+y^2 \le b^2 \}

  • 解析

  积分区域为圆环域,被积函数含有 x2+y2x^2+y^2,极坐标范围为:0θ2π,aρb0 \le \theta \le 2\pi, \, a \le \rho \le b

  I=02πdθabeρ2ρdρ=2π[12eρ2]ab=π(eb2ea2)I = \int_{0}^{2\pi} \mathrm{d}\theta \int_{a}^{b} e^{\rho^2} \rho \mathrm{d}\rho = 2\pi \cdot \left[ \frac{1}{2} e^{\rho^2} \right]_{a}^{b} = \pi (e^{b^2} - e^{a^2})

3. 利用二重积分求平面图形面积(课后题 8.(1))#

  • 题目:求曲线 xy=a2xy=a^2 与直线 x+y=52ax+y=\frac{5}{2}aa>0a>0)所围成的平面图形面积。

  • 解析

  1. 求解交点横坐标:

     x+a2x=52a    2x25ax+2a2=0    x1=a2,x2=2ax + \frac{a^2}{x} = \frac{5}{2}a \implies 2x^2 - 5ax + 2a^2 = 0 \implies x_1 = \frac{a}{2}, \, x_2 = 2a

  2. 确定积分域:当 x[a2,2a]x \in [\frac{a}{2}, 2a] 时,上限为直线 y=52axy = \frac{5}{2}a - x,下限为双曲线 y=a2xy = \frac{a^2}{x}

  3. 计算面积 AA

     A=a/22a(52axa2x)dx=[52ax12x2a2lnx]a/22aA = \int_{a/2}^{2a} \left( \frac{5}{2}a - x - \frac{a^2}{x} \right) \mathrm{d}x = \left[ \frac{5}{2}ax - \frac{1}{2}x^2 - a^2\ln x \right]_{a/2}^{2a}

     代入计算得:

     A=(5a22a2a2ln(2a))(54a218a2a2lna2)=a2(1582ln2)A = \left( 5a^2 - 2a^2 - a^2\ln(2a) \right) - \left( \frac{5}{4}a^2 - \frac{1}{8}a^2 - a^2\ln\frac{a}{2} \right) = a^2 \left( \frac{15}{8} - 2\ln 2 \right)


第三节 三重积分的计算#

一、 核心计算方法(直角坐标系)#

在直角坐标系下,计算三重积分的核心思想是化三重积分为累次积分(单积分与二重积分的结合)。主要有两种切分路径:

1. 投影法(又称“先单后重”或“拉条法”)#

  • 适用场景:积分域 (V)(V) 的“顶”和“底”边界曲面易于用 z=z2(x,y)z=z_2(x,y)z=z1(x,y)z=z_1(x,y) 表示,且 (V)(V)xOyxOy 平面上的投影区域 (σ)(\sigma) 较简单。

  • 计算公式(公式 3.1):

  (V)f(x,y,z)dV=(σ)[z1(x,y)z2(x,y)f(x,y,z)dz]dσ\iiint_{(V)} f(x,y,z) \mathrm{d}V = \iint_{(\sigma)} \left[ \int_{z_1(x,y)}^{z_2(x,y)} f(x,y,z) \mathrm{d}z \right] \mathrm{d}\sigma

2. 截面法(又称“先重后单”或“切片法”)#

  • 适用场景:被积函数主要与 zz 相关,或者用平行于 xOyxOy 平面的平面 z=constz = \text{const} 去截积分域 (V)(V) 时,所得的截面 (σz)(\sigma_z) 的边界和面积易于表达。

  • 计算公式(公式 3.2):

  (V)f(x,y,z)dV=ab[(σz)f(x,y,z)dσ]dz\iiint_{(V)} f(x,y,z) \mathrm{d}V = \int_{a}^{b} \left[ \iint_{(\sigma_z)} f(x,y,z) \mathrm{d}\sigma \right] \mathrm{d}z

  其中 aabb 分别是积分域 (V)(V)zz 轴上的投影下限与上限。


二、 三维坐标变换与曲线坐标系#

对于边界含有圆柱面、锥面、球面等几何体的区域,直角坐标系往往计算繁琐,需要引入三维空间下的坐标变换。

1. 三维空间下的雅可比行列式(Jacobian)#

设变换关系为 x=x(u,v,w),y=y(u,v,w),z=z(u,v,w)x=x(u,v,w), y=y(u,v,w), z=z(u,v,w),其雅可比行列式为:

J=(x,y,z)(u,v,w)=xuxvxwyuyvywzuzvzwJ = \frac{\partial(x,y,z)}{\partial(u,v,w)} = \begin{vmatrix} x_u & x_v & x_w \\ y_u & y_v & y_w \\ z_u & z_v & z_w \end{vmatrix}

体积微元变换关系为:dV=Jdudvdw\mathrm{d}V = |J| \mathrm{d}u\mathrm{d}v\mathrm{d}w

2. 柱面坐标系(Cylindrical Coordinates)#

  • 变换公式x=ρcosθ,y=ρsinθ,z=zx = \rho \cos\theta, \, y = \rho \sin\theta, \, z = z

  • 体积微元(公式 3.11):dV=ρdρdθdz\mathrm{d}V = \rho \mathrm{d}\rho \mathrm{d}\theta \mathrm{d}z

  • 适用场景:积分域边界含有圆柱面(如 x2+y2=R2x^2+y^2=R^2)或旋转抛物面(如 z=x2+y2z=x^2+y^2)等。

3. 球面坐标系(Spherical Coordinates)#

  • 变换公式(公式 3.12):

  x=rsinφcosθ,y=rsinφsinθ,z=rcosφx = r \sin\varphi \cos\theta, \quad y = r \sin\varphi \sin\theta, \quad z = r \cos\varphi

  (其中 r0,0φπ,0θ2πr \ge 0, \, 0 \le \varphi \le \pi, \, 0 \le \theta \le 2\pi

  • 体积微元(公式 3.13):dV=r2sinφdrdφdθ\mathrm{d}V = r^2 \sin\varphi \mathrm{d}r \mathrm{d}\varphi \mathrm{d}\theta

  • 适用场景:积分域边界含有球面(如 x2+y2+z2=R2x^2+y^2+z^2=R^2)或圆锥面(如 z=x2+y2z=\sqrt{x^2+y^2})等。


三、 核心重难点:对称性的推广与应用#

合理利用对称性可以大幅缩减三重积分的计算量,常用于以下两个维度:

  1. 奇偶对称性:若积分域 (V)(V) 关于某坐标面(如 xOyxOy 平面)对称,而被积函数关于该面对应变量(如 zz)为奇函数,则积分为 00

   * 例 3.4 启示:在计算 (V)(x+y+z)dV\iiint_{(V)} (x+y+z) \mathrm{d}V 时,区域 (V)(V) 关于 yOzyOzxOzxOz 面对称,由于 xxyy 是奇函数,故 xdV=0\iiint x \mathrm{d}V = 0ydV=0\iiint y \mathrm{d}V = 0,积分直接简化为 zdV\iiint z \mathrm{d}V

  1. 轮换对称性:若积分域关于变量 x,y,zx, y, z 具有某种轮换对称,则可通过变量对调简化积分项。

四、 典型习题与例题精选分析#

1. 概念判断题(课后习题 2)#

  • 题目 (1):设 (V)(V) 是球体 x2+y2+z24x^2+y^2+z^2 \le 4(V1)(V_1) 为其上半球体,判断下列等式是否正确:

  (V)(x+y+z)2dV=2(V1)(x+y+z)2dV\iiint_{(V)} (x+y+z)^2 \mathrm{d}V = 2 \iiint_{(V_1)} (x+y+z)^2 \mathrm{d}V

  • 解析

  将积分项展开:(x+y+z)2=x2+y2+z2+2(xy+yz+zx)(x+y+z)^2 = x^2+y^2+z^2 + 2(xy+yz+zx)

  由于球体 (V)(V) 关于三个坐标面均对称:

  * 交叉项 xy,yz,zxxy, yz, zx 在整球 (V)(V) 上的积分为 00

  * 在上半球 (V1)(V_1) 上,利用关于 xOzxOzyOzyOz 面的对称性,交叉项的积分同样为 00

  * 剩余项 x2+y2+z2x^2+y^2+z^2 是关于 zz 的偶函数,由于上半球 (V1)(V_1)z0z \ge 0,在整球上的积分值确实等于上半球的 22 倍。

  * 结论:该式正确

  • 题目 (4):判断等式是否正确:

  (V)3(x2+y2+z2)dV=3(V)4dV\iiint_{(V)} 3(x^2+y^2+z^2) \mathrm{d}V = 3 \iiint_{(V)} 4 \mathrm{d}V

  • 解析

  在积分区域内,x2+y2+z2x^2+y^2+z^2 是自变量,其值在 0044 之间变化,不能在积分符号内部直接代入边界最大值 44

  * 结论:该式错误

2. 重难点题型:两球体交集的体积与积分(课后习题 4.(10) / 对应例 3.9)#

  • 题目:计算 (V)z2dV\iiint_{(V)} z^2 \mathrm{d}V,其中 (V)(V) 为两球体 x2+y2+z2R2x^2+y^2+z^2 \le R^2x2+y2+z22Rzx^2+y^2+z^2 \le 2Rz 的公共部分。

  • 解法对比与分析(本题有三种解法,其繁简对比极具教学意义):

  * 方法一:截面法(最简解法)

    两球面的交线处满足 R2z2=2Rzz2    z=R2R^2 - z^2 = 2Rz - z^2 \implies z = \frac{R}{2}

    * 当 0zR20 \le z \le \frac{R}{2} 时,横截面由下半部分球体边界限制,截面半径满足 r12=2Rzz2r_1^2 = 2Rz - z^2。其截面积为 π(2Rzz2)\pi(2Rz - z^2)

    * 当 R2zR\frac{R}{2} \le z \le R 时,横截面由上半部分球体边界限制,截面半径满足 r22=R2z2r_2^2 = R^2 - z^2。其截面积为 π(R2z2)\pi(R^2 - z^2)

    利用截面法列式:

    I=0R2z2π(2Rzz2)dz+R2Rz2π(R2z2)dzI = \int_{0}^{\frac{R}{2}} z^2 \cdot \pi(2Rz - z^2) \mathrm{d}z + \int_{\frac{R}{2}}^{R} z^2 \cdot \pi(R^2 - z^2) \mathrm{d}z

    通过简单的一元多项式积分即可求得:

    I=59480πR5I = \frac{59}{480}\pi R^5

  * 方法二:球面坐标法(运算量适中)

    在球面坐标系下,交线对应的夹角为 cosφ=12    φ=π3\cos\varphi = \frac{1}{2} \implies \varphi = \frac{\pi}{3}

    积分域需要以 φ=π3\varphi = \frac{\pi}{3} 为界拆分为两部分:

    I=02πdθ0π3dφ0R(rcosφ)2r2sinφdr+02πdθπ3π2dφ02Rcosφ(rcosφ)2r2sinφdrI = \int_{0}^{2\pi} \mathrm{d}\theta \int_{0}^{\frac{\pi}{3}} \mathrm{d}\varphi \int_{0}^{R} (r\cos\varphi)^2 \cdot r^2\sin\varphi \mathrm{d}r + \int_{0}^{2\pi} \mathrm{d}\theta \int_{\frac{\pi}{3}}^{\frac{\pi}{2}} \mathrm{d}\varphi \int_{0}^{2R\cos\varphi} (r\cos\varphi)^2 \cdot r^2\sin\varphi \mathrm{d}r

    计算后同样可得 59480πR5\frac{59}{480}\pi R^5。该方法展示了根据边界分区域积分的技巧。

3. 极坐标/球面坐标计算体积(课后习题 6.(3))#

  • 题目:求由曲面 (x2+y2+z2)2=a3z(x^2+y^2+z^2)^2 = a^3za>0a>0)所围成的立体的体积。

  • 解析

  由于边界方程含有 (x2+y2+z2)(x^2+y^2+z^2) 以及 zz,极适合采用球面坐标系。

  将方程化为球面坐标:

  (r2)2=a3(rcosφ)    r3=a3cosφ    r=acosφ3(r^2)^2 = a^3(r\cos\varphi) \implies r^3 = a^3\cos\varphi \implies r = a \sqrt[3]{\cos\varphi}

  因为 r0r \ge 0,要求 cosφ0    φ[0,π2]\cos\varphi \ge 0 \implies \varphi \in [0, \frac{\pi}{2}](立体位于上半空间)。

  立体的体积 VV 为:

  V=(V)dV=02πdθ0π2sinφdφ0acos1/3φr2drV = \iiint_{(V)} \mathrm{d}V = \int_{0}^{2\pi} \mathrm{d}\theta \int_{0}^{\frac{\pi}{2}} \sin\varphi \mathrm{d}\varphi \int_{0}^{a\cos^{1/3}\varphi} r^2 \mathrm{d}r

  先对 rr 积分:

  0acos1/3φr2dr=a33cosφ\int_{0}^{a\cos^{1/3}\varphi} r^2 \mathrm{d}r = \frac{a^3}{3}\cos\varphi

  代入外层积分:

  V=2πa330π2cosφsinφdφ=23πa3[12sin2φ]0π2=13πa3V = 2\pi \cdot \frac{a^3}{3} \int_{0}^{\frac{\pi}{2}} \cos\varphi \sin\varphi \mathrm{d}\varphi = \frac{2}{3}\pi a^3 \cdot \left[ \frac{1}{2}\sin^2\varphi \right]_{0}^{\frac{\pi}{2}} = \frac{1}{3}\pi a^3


第四节 含参变量的积分与反常重积分#

一、 含参变量积分的性质#

D=[a,b]×[c,d]D = [a,b] \times [c,d]f(x,y)f(x,y) 是定义在 DD 上的连续函数。定义含参变量 yy 的积分:

F(y)=abf(x,y)dx(y[c,d])F(y) = \int_{a}^{b} f(x,y) \mathrm{d}x \quad (y \in [c,d])

1. 基本性质定理(常数积分限)#

  • 定理 4.1(连续性):若 fC(D)f \in C(D),则 F(y)F(y)[c,d][c,d] 上连续。其重要极限形式为(极限与积分号可交换顺序):

  limyy0abf(x,y)dx=ab[limyy0f(x,y)]dx\lim_{y \to y_0} \int_{a}^{b} f(x,y) \mathrm{d}x = \int_{a}^{b} \left[ \lim_{y \to y_0} f(x,y) \right] \mathrm{d}x

  • 定理 4.2(可导性 / 积分号下求导):若 fC(D)f \in C(D) 且偏导数 fyC(D)f_y \in C(D),则 F(y)F(y)[c,d][c,d] 上可导,且:

  F(y)=ddyabf(x,y)dx=abf(x,y)ydxF'(y) = \frac{\mathrm{d}}{\mathrm{d}y} \int_{a}^{b} f(x,y) \mathrm{d}x = \int_{a}^{b} \frac{\partial f(x,y)}{\partial y} \mathrm{d}x

  • 定理 4.3(积分顺序交换性):若 fC(D)f \in C(D),则累次积分可以交换顺序:

  cd(abf(x,y)dx)dy=ab(cdf(x,y)dy)dx\int_{c}^{d} \left( \int_{a}^{b} f(x,y) \mathrm{d}x \right) \mathrm{d}y = \int_{a}^{b} \left( \int_{c}^{d} f(x,y) \mathrm{d}y \right) \mathrm{d}x

2. 变限参数积分求导(Leibniz 积分法则)#

  • 定理 4.5:若 f(x,y)f(x,y)fy(x,y)f_y(x,y)DD 上连续,且积分限函数 x1(y),x2(y)x_1(y), x_2(y)[c,d][c,d] 上可导,则 F(y)=x1(y)x2(y)f(x,y)dxF(y) = \int_{x_1(y)}^{x_2(y)} f(x,y) \mathrm{d}x 可导,且:

  F(y)=x1(y)x2(y)f(x,y)ydx+f[x2(y),y]x2(y)f[x1(y),y]x1(y)F'(y) = \int_{x_1(y)}^{x_2(y)} \frac{\partial f(x,y)}{\partial y} \mathrm{d}x + f[x_2(y), y] \cdot x_2'(y) - f[x_1(y), y] \cdot x_1'(y)


二、 反常重积分理论#

反常重积分主要分为两类:无界区域上的重积分无界函数的重积分

1. 无界区域的二重积分(定义 4.1)#

(σ)(\sigma) 为无界区域。若对于任何趋于整个区域的单调有界区域序列 (σn)(\sigma_n),极限 limn(σn)f(x,y)dσ\lim_{n \to \infty} \iint_{(\sigma_n)} f(x,y) \mathrm{d}\sigma 均存在且相等,则称该反常重积分收敛。

  • 定理 4.6(无穷远处的收敛判别法):若在无界区域上存在常数 M>0M > 0α>2\alpha > 2,当 ρ=x2+y2ρ0\rho = \sqrt{x^2+y^2} \ge \rho_0 时有:

  f(x,y)Mρα|f(x,y)| \le \frac{M}{\rho^\alpha}

  则反常重积分 (σ)f(x,y)dσ\iint_{(\sigma)} f(x,y) \mathrm{d}\sigma 收敛。(若 α2\alpha \le 2,在正值函数下通常发散)。

2. 无界函数的二重积分(奇点位于原点,定义 4.2)#

P0(x0,y0)P_0(x_0, y_0) 为函数的无界奇点。

  • 定理 4.7(奇点处的收敛判别法):若在奇点邻域内有 f(x,y)Mρα|f(x,y)| \le \frac{M}{\rho^\alpha},则:

  * 当 α<2\alpha < 2 时,反常重积分收敛

  * 当 α2\alpha \ge 2 时,反常重积分发散。


三、 重难点突破:两个经典数学工具的诞生#

1. 交换积分顺序简化一元积分(例 4.2)#

通过构造含参变量积分,并利用定理 4.3 交换积分顺序,可以避开直接积分的困境。

  • 典型方法:求解 01xbxalnxdx\int_{0}^{1} \frac{x^b - x^a}{\ln x} \mathrm{d}x

  注意到被积函数可以写为参数积分:xbxalnx=abxydy\frac{x^b - x^a}{\ln x} = \int_{a}^{b} x^y \mathrm{d}y

  从而原积分为:

  I=01dxabxydy=abdy01xydx=ab1y+1dy=ln1+b1+aI = \int_{0}^{1} \mathrm{d}x \int_{a}^{b} x^y \mathrm{d}y = \int_{a}^{b} \mathrm{d}y \int_{0}^{1} x^y \mathrm{d}x = \int_{a}^{b} \frac{1}{y+1} \mathrm{d}y = \ln \frac{1+b}{1+a}

2. 高斯概率积分的推导(例 4.5)#

利用二重反常积分可证出概率统计中极其重要的高斯积分值 π\sqrt{\pi}

I=+ex2dxI = \int_{-\infty}^{+\infty} e^{-x^2} \mathrm{d}x。通过构建全平面的二重反常积分:

R2e(x2+y2)dxdy=(+ex2dx)2\iint_{\mathbb{R}^2} e^{-(x^2+y^2)} \mathrm{d}x\mathrm{d}y = \left( \int_{-\infty}^{+\infty} e^{-x^2} \mathrm{d}x \right)^2

另一方面,在极坐标下计算该反常二重积分:

R2e(x2+y2)dxdy=02πdθ0+eρ2ρdρ=π\iint_{\mathbb{R}^2} e^{-(x^2+y^2)} \mathrm{d}x\mathrm{d}y = \int_{0}^{2\pi} \mathrm{d}\theta \int_{0}^{+\infty} e^{-\rho^2} \rho \mathrm{d}\rho = \pi

从而证得:+ex2dx=π\int_{-\infty}^{+\infty} e^{-x^2} \mathrm{d}x = \sqrt{\pi}


四、 重点习题精选分析#

1. 积分与极限顺序交换(课后题 A.1.(3))#

  • 题目:求极限 limα0011+α2x2dx\lim_{\alpha \to 0} \int_{0}^{1} \sqrt{1+\alpha^2-x^2} \mathrm{d}x

  • 解析

  被积函数 f(x,α)=1+α2x2f(x,\alpha) = \sqrt{1+\alpha^2-x^2} 在区域 [0,1]×[δ,δ][0,1] \times [-\delta, \delta] 上连续。由定理 4.1,极限与积分号可直接交换:

  limα0011+α2x2dx=01(limα01+α2x2)dx=011x2dx\lim_{\alpha \to 0} \int_{0}^{1} \sqrt{1+\alpha^2-x^2} \mathrm{d}x = \int_{0}^{1} \left( \lim_{\alpha \to 0} \sqrt{1+\alpha^2-x^2} \right) \mathrm{d}x = \int_{0}^{1} \sqrt{1-x^2} \mathrm{d}x

  此积分为四分之一单位圆面积,结果为 π4\frac{\pi}{4}

2. Feynman积分求导妙用(课后题 A.3.(3) - 经典难题)#

  • 题目:计算积分 I(a,b)=0π/2ln(a2sin2x+b2cos2x)dxI(a,b) = \int_{0}^{\pi/2} \ln(a^2\sin^2 x + b^2\cos^2 x) \mathrm{d}xa>0,b>0a>0, b>0)。

  • 解析

  1. 设 I(a,b)I(a,b),对参数 aa 求偏导:

     Ia=0π/22asin2xa2sin2x+b2cos2xdx\frac{\partial I}{\partial a} = \int_{0}^{\pi/2} \frac{2a\sin^2 x}{a^2\sin^2 x + b^2\cos^2 x} \mathrm{d}x

  2. 分子分母同除以 cos2x\cos^2 x,令 t=tanxt = \tan x 进行代换:

     Ia=0+2at2(a2t2+b2)(1+t2)dt\frac{\partial I}{\partial a} = \int_{0}^{+\infty} \frac{2a t^2}{(a^2t^2+b^2)(1+t^2)} \mathrm{d}t

  3. 利用部分分式拆分后积分,可算得极为简洁的偏导关系式:

     Ia=πa+b\frac{\partial I}{\partial a} = \frac{\pi}{a+b}

  4. 同理,由对称性可知:Ib=πa+b\frac{\partial I}{\partial b} = \frac{\pi}{a+b}

  5. 积分恢复原函数,得 I(a,b)=πln(a+b)+CI(a,b) = \pi \ln(a+b) + C

  6. 代入特殊值 a=b=1a=b=1 定常数:I(1,1)=0π/2ln(1)dx=0    πln(2)+C=0    C=πln2I(1,1) = \int_0^{\pi/2} \ln(1) \mathrm{d}x = 0 \implies \pi\ln(2) + C = 0 \implies C = -\pi\ln 2

  * 结论I(a,b)=πlna+b2I(a,b) = \pi \ln \frac{a+b}{2}

3. 积分号下交换积分计算(课后题 A.4 - Frullani积分)#

  • 题目:计算积分 0+eaxebxxdx\int_{0}^{+\infty} \frac{e^{-ax} - e^{-bx}}{x} \mathrm{d}xa>0,b>0a>0, b>0)。

  • 解析

  由于 eaxebxx=abeyxdy\frac{e^{-ax} - e^{-bx}}{x} = \int_{a}^{b} e^{-yx} \mathrm{d}y,可将原积分化为含参变量的反常重积分:

  I=0+dxabeyxdyI = \int_{0}^{+\infty} \mathrm{d}x \int_{a}^{b} e^{-yx} \mathrm{d}y

  交换积分顺序:

  I=abdy0+eyxdx=ab[1yeyx]0+dy=ab1ydy=lnbaI = \int_{a}^{b} \mathrm{d}y \int_{0}^{+\infty} e^{-yx} \mathrm{d}x = \int_{a}^{b} \left[ -\frac{1}{y} e^{-yx} \right]_{0}^{+\infty} \mathrm{d}y = \int_{a}^{b} \frac{1}{y} \mathrm{d}y = \ln \frac{b}{a}


第五节 重积分的应用#

一、 重积分的微元法基础(定义与定理)#

将一元函数定积分的“微元法”推广至多元重积分,其核心是引入区域函数对域的导数

1. 区域函数与对域的导数(定义 5.1)#

设平面区域 (σ)(\sigma),对其任一子域 (Δσ)(\Delta\sigma),若存在某一对应的物理或几何量 F((Δσ))F((\Delta\sigma)),则称 FF区域函数

若对于子域内的一点 M(x,y)M(x,y),当子域直径 d0d \to 0 时,极限:

dFdσ=limd0F((Δσ))Δσ=f(x,y)\frac{\mathrm{d}F}{\mathrm{d}\sigma} = \lim_{d \to 0} \frac{F((\Delta\sigma))}{\Delta \sigma} = f(x,y)

存在,则称 f(x,y)f(x,y) 为区域函数 FF 在点 MM 处对区域面积的导数。其对应的线性主部即为该物理量的微元(微分)

dF=f(x,y)dσ\mathrm{d}F = f(x,y) \mathrm{d}\sigma

2. 变限重积分对域的导数(定理 5.2)#

设连续函数 f(M)f(M) 在变域上的重积分定义为区域函数:Φ((Δσ))=(Δσ)f(M)dσ\Phi((\Delta\sigma)) = \iint_{(\Delta\sigma)} f(M) \mathrm{d}\sigma。由积分中值定理及极限过程可证:

dΦdσ=f(x,y)\frac{\mathrm{d}\Phi}{\mathrm{d}\sigma} = f(x,y)

该定理表明,连续函数在变域上的重积分作为区域函数,其对域的导数就等于被积函数在该点的值。这就为物理微元的合理性提供了严格的数学支撑。

3. 微元法的标准步骤#

若某一物理量 Φ\Phi 关于区域具有可加性,则:

  1. 建立微元:在微小区域 dσ\mathrm{d}\sigma(或 dV\mathrm{d}V)上,将分布不均匀的物理量看作均匀分布,写出微元表达式 dΦ=f(x,y)dσ\mathrm{d}\Phi = f(x,y)\mathrm{d}\sigma

  2. 累加求和(求积分):将微元在整个区域上进行重积分,得到总量 Φ=(σ)f(x,y)dσ\Phi = \iint_{(\sigma)} f(x,y)\mathrm{d}\sigma


二、 典型物理应用公式系统#

根据微元法,本节推导并归纳了以下重要的物理应用公式(以均匀分布或已知密度分布为例):

1. 质心与静矩(Center of Mass & Static Moments)#

对于密度分布为常数 KK(或变化密度为 f(x,y,z)f(x,y,z))的空间物体 (V)(V)

  • 静矩公式(对三个坐标面的静矩):

  Myz=(V)xdm,Mzx=(V)ydm,Mxy=(V)zdmM_{yz} = \iiint_{(V)} x \mathrm{d}m, \quad M_{zx} = \iiint_{(V)} y \mathrm{d}m, \quad M_{xy} = \iiint_{(V)} z \mathrm{d}m

  • 质心坐标公式

  xˉ=Myzm,yˉ=Mzxm,zˉ=Mxym\bar{x} = \frac{M_{yz}}{m}, \quad \bar{y} = \frac{M_{zx}}{m}, \quad \bar{z} = \frac{M_{xy}}{m}

  其中 m=(V)dmm = \iiint_{(V)} \mathrm{d}m 为物体的总质量。

2. 转动惯量(Moment of Inertia)#

物体的转动惯量微元满足 dI=r2dm\mathrm{d}I = r^2 \mathrm{d}m,其中 rr 为微元到转轴的垂直距离。

  • 平面薄板 (σ)(\sigma)xx 轴和 yy 轴的转动惯量

  Ix=(σ)y2dm,Iy=(σ)x2dmI_x = \iint_{(\sigma)} y^2 \mathrm{d}m, \quad I_y = \iint_{(\sigma)} x^2 \mathrm{d}m

  • 平面薄板对原点(极轴)的转动惯量

  I0=Ix+Iy=(σ)(x2+y2)dmI_0 = I_x + I_y = \iint_{(\sigma)} (x^2+y^2) \mathrm{d}m

3. 引力计算(Gravitational Attraction)#

设一质量为 mm' 的质点置于 P(x0,y0,z0)P'(x_0,y_0,z_0),连续分布物体 (V)(V) 密度为 μ\mu,则 (V)(V) 对该质点的万有引力矢量为:

F=Gm(V)rr3μdV\vec{F} = G m' \iiint_{(V)} \frac{\vec{r}}{r^3} \mu \mathrm{d}V

其中 r\vec{r} 是从 PP' 指向物体内部微元点的向量。


三、 本节重难点解析#

1. 引力场极限的物理学推广(例 5.4)#

例 5.4 求解了平面圆盘对轴线上单质点的万有引力:

Fz=2πμG(1hh2+R2)F_z = -2\pi\mu G \left( 1 - \frac{h}{\sqrt{h^2+R^2}} \right)

在物理理论分析中,若圆盘半径 RR \to \infty(即圆盘被视为无限大平面),极限结果为:

limRFz=2πμG\lim_{R \to \infty} F_z = -2\pi \mu G

此极限值表明:无限大均匀带质量平面在其附近空间产生的引力场是均匀场,其大小仅取决于面密度 μ\mu,与质点到平面的距离 hh 无关。

2. 质心计算中的截面法与柱坐标法的合理选择#

在求解旋转对称体(如半球体、半椭球体)的质心时:

  • 采用截面法(Slicing Method)通常比柱/球面坐标法更为简单。因为截面法能直接把三重积分转化为关于对称轴一元变量的简单积分。

四、 重点习题与例题精选分析#

1. 均匀半椭球体的质心坐标(课后习题 A.2.(1))#

  • 题目:求边界为上半椭球面 z=c1x2a2y2b2z = c \sqrt{1 - \frac{x^2}{a^2} - \frac{y^2}{b^2}} 及平面 z=0z=0a>0,b>0,z0a>0, b>0, z \ge 0)所围成的均匀物体的质心。

  • 解析

  1. 对称性分析:由于椭球体关于 yOzyOz 轴和 xOzxOz 轴对称且密度均匀,因此其质心必在 zz 轴上,即:

     xˉ=0,yˉ=0\bar{x} = 0, \quad \bar{y} = 0

  2. 总体积(质量 mm

     已知整椭球体积为 43πabc\frac{4}{3}\pi abc,均匀半椭球体的体积为:

     V=23πabcV = \frac{2}{3}\pi abc

  3. 利用截面法求静矩 MxyM_{xy}

     在高度 z[0,c]z \in [0,c] 处,平行于 xOyxOy 面的截面 (σz)(\sigma_z) 是一个椭圆:

     x2a2+y2b21z2c2\frac{x^2}{a^2} + \frac{y^2}{b^2} \le 1 - \frac{z^2}{c^2}

     该椭圆的半轴分别为 a1z2c2a \sqrt{1 - \frac{z^2}{c^2}}b1z2c2b \sqrt{1 - \frac{z^2}{c^2}}。因此其截面积为:

     S(z)=πa1z2c2b1z2c2=πab(1z2c2)S(z) = \pi \cdot a \sqrt{1 - \frac{z^2}{c^2}} \cdot b \sqrt{1 - \frac{z^2}{c^2}} = \pi ab \left( 1 - \frac{z^2}{c^2} \right)

     静矩为:

     Mxy=0czS(z)dz=πab0cz(1z2c2)dz=πab[z22z44c2]0c=14πabc2M_{xy} = \int_{0}^{c} z \cdot S(z) \mathrm{d}z = \pi ab \int_{0}^{c} z \left( 1 - \frac{z^2}{c^2} \right) \mathrm{d}z = \pi ab \left[ \frac{z^2}{2} - \frac{z^4}{4c^2} \right]_{0}^{c} = \frac{1}{4} \pi abc^2

  4. 质心高度 zˉ\bar{z}

     zˉ=MxyV=14πabc223πabc=38c\bar{z} = \frac{M_{xy}}{V} = \frac{\frac{1}{4} \pi abc^2}{\frac{2}{3}\pi abc} = \frac{3}{8}c

  * 结论:均匀半椭球体的质心坐标为 (0,0,38c)(0, 0, \frac{3}{8}c)。(此结论在 a=b=ca=b=c 时退化为均匀半球体质心高 38a\frac{3}{8}a 的经典结论)。

2. 平行轴定理的推导与应用(课后习题 B.4 / B.5 核心)#

  • 定理内容:若物体 (V)(V) 对通过其质心且平行于给定轴 LL 的轴 LcL_c 的转动惯量为 IcI_c,则物体对轴 LL 的转动惯量 ILI_L 为:

  IL=Ic+md2I_L = I_c + m d^2

  其中 mm 为物体总质量,dd 为轴 LLLcL_c 之间的距离。

  • 应用举例(习题 B.5):利用平行轴定理求半径为 RR 的均匀球体对于任一条切线 TT 的转动惯量 ITI_T

  1. 通过球心的轴 LcL_c 的转动惯量为经典的 Ic=25mR2I_c = \frac{2}{5}mR^2

  2. 球的切线 TT 与通过球心的平行轴之间的距离为 d=Rd = R

  3. 直接套用平行轴定理:

     IT=Ic+mR2=25mR2+mR2=75mR2I_T = I_c + m R^2 = \frac{2}{5}mR^2 + mR^2 = \frac{7}{5}mR^2


第六节 第一型线积分与面积分#

一、 第一型线积分(对弧长的曲线积分)#

1. 概念与定义(公式 6.1)#

设曲线弧段 (C)(C) 是有界的,其长度微元 ds>0\mathrm{d}s > 0。函数 f(x,y,z)f(x,y,z)(C)(C) 上的第一型线积分定义为:

(C)f(x,y,z)ds=limd0k=1nf(ξk,ηk,ζk)Δsk\int_{(C)} f(x,y,z) \mathrm{d}s = \lim_{d \to 0} \sum_{k=1}^n f(\xi_k, \eta_k, \zeta_k) \Delta s_k

  • 核心特征:由于弧长微元 Δsk>0\Delta s_k > 0 恒正,第一型线积分的值与积分路径的走向无关。

2. 参数计算公式(公式 6.3)#

(C)(C) 的参数方程为 x=x(t),y=y(t),z=z(t)x=x(t), y=y(t), z=z(t)αtβ\alpha \le t \le \beta,且 α<β\alpha < \beta),则:

(C)f(x,y,z)ds=αβf[x(t),y(t),z(t)]x˙2(t)+y˙2(t)+z˙2(t)dt\int_{(C)} f(x,y,z) \mathrm{d}s = \int_{\alpha}^{\beta} f[x(t), y(t), z(t)] \sqrt{\dot{x}^2(t) + \dot{y}^2(t) + \dot{z}^2(t)} \mathrm{d}t


二、 第一型面积分(对面积的曲面积分)#

1. 曲面面积与面积微元(公式 6.8 & 6.11)#

  • 参数式曲面 r(u,v)=(x,y,z)\vec{r}(u,v) = (x,y,z):面积微元为 dS=ru×rvdudv\mathrm{d}S = \| \vec{r}_u \times \vec{r}_v \| \mathrm{d}u\mathrm{d}v

  • 显式曲面 z=z(x,y)z = z(x,y)

  面积微元满足投影关系:dS=dxdycosγ\mathrm{d}S = \frac{\mathrm{d}x\mathrm{d}y}{\cos\gamma},其中 cosγ\cos\gamma 是曲面法向量与 zz 轴正向的夹角余弦,即:

  dS=1+zx2+zy2dxdy\mathrm{d}S = \sqrt{1 + z_x^2 + z_y^2} \mathrm{d}x\mathrm{d}y

2. 第一型面积分计算公式(公式 6.15)#

若曲面 (S)(S) 的方程为 z=z(x,y)z=z(x,y),在 xOyxOy 面投影为 (σ)(\sigma),则:

(S)f(x,y,z)dS=(σ)f[x,y,z(x,y)]1+zx2+zy2dxdy\iint_{(S)} f(x,y,z) \mathrm{d}S = \iint_{(\sigma)} f[x, y, z(x,y)] \sqrt{1 + z_x^2 + z_y^2} \mathrm{d}x\mathrm{d}y


三、 解题精髓:对称性化简技术#

第一型线积分与面积分在几何对称体(如圆、球)上的应用极具技巧性,合理利用对称性往往能免去繁琐的积分计算。

1. 奇偶对称性(消项技术)#

  • 若积分路径 (C)(C) 关于坐标原点(或某坐标轴)对称,而被积函数 f(x,y,z)f(x,y,z) 关于其变量为奇函数,则积分为 00

  • 例 6.2 启示:计算 (C)yds\int_{(C)} y \mathrm{d}s, parabolic 弧段 y2=2xy^2 = 2x[2,2][-2, 2] 内关于 xx 轴对称。由于被积函数 yy 是奇函数,可直接断定积分为 00

2. 轮换对称性(均分技术)#

  • 空间曲线圆周的轮换对称(例 6.3)

  若积分路径是球面与过原点平面的交线圆周 (C)(C),由于其对于 x,y,zx, y, z 的几何地位完全相同:

  (C)x2ds=(C)y2ds=(C)z2ds=13(C)(x2+y2+z2)ds\oint_{(C)} x^2 \mathrm{d}s = \oint_{(C)} y^2 \mathrm{d}s = \oint_{(C)} z^2 \mathrm{d}s = \frac{1}{3} \oint_{(C)} (x^2+y^2+z^2) \mathrm{d}s

  • 球面面积分的轮换对称(例 6.9)

  若积分域为球面 (S):x2+y2+z2=a2(S): x^2+y^2+z^2=a^2,则:

  (S)x2dS=(S)y2dS=(S)z2dS=13(S)(x2+y2+z2)dS=13a24πa2=43πa4\iint_{(S)} x^2 \mathrm{d}S = \iint_{(S)} y^2 \mathrm{d}S = \iint_{(S)} z^2 \mathrm{d}S = \frac{1}{3} \iint_{(S)} (x^2+y^2+z^2) \mathrm{d}S = \frac{1}{3} a^2 \cdot 4\pi a^2 = \frac{4}{3}\pi a^4


四、 精选课后习题详细解答(打勾/标记题目)#

1. 抛物线弧段的第一型线积分(习题 6.6 A.1.(1))#

  • 题目:计算 (C)yds\int_{(C)} y \mathrm{d}s,其中 (C)(C) 为抛物线 y2=2xy^2=2x 上自点 (0,0)(0,0) 到点 (2,2)(2,2) 的一段弧。

  • 解析

  1. 采用 yy 作为积分参数:x=12y2x = \frac{1}{2}y^2,其中 y[0,2]y \in [0,2]

  2. 计算弧长微元:

     ds=1+(dxdy)2dy=1+y2dy\mathrm{d}s = \sqrt{1 + \left(\frac{\mathrm{d}x}{\mathrm{d}y}\right)^2} \mathrm{d}y = \sqrt{1+y^2} \mathrm{d}y

  3. 列式并积分:

     (C)yds=02y1+y2dy=[13(1+y2)3/2]02=5513\int_{(C)} y \mathrm{d}s = \int_{0}^{2} y \sqrt{1+y^2} \mathrm{d}y = \left[ \frac{1}{3} (1+y^2)^{3/2} \right]_{0}^{2} = \frac{5\sqrt{5} - 1}{3}

2. 空间曲面交线的第一型线积分(习题 6.6 A.1.(5))#

  • 题目:计算 (C)x2ds\oint_{(C)} x^2 \mathrm{d}s,其中 (C)(C) 为圆周 x2+y2+z2=4x^2+y^2+z^2=4 与平面 z=3z=\sqrt{3} 的交线。

  • 解析

  1. 分析交线几何结构:将 z=3z=\sqrt{3} 代入球面方程,得:

     x2+y2=43=1x^2+y^2 = 4 - 3 = 1

     此交线 (C)(C) 为位于高度 z=3z = \sqrt{3} 处的水平单位圆,半径为 R=1R=1,其周长为 L=2πL = 2\pi

  2. 利用圆周的轮换对称性:

     (C)x2ds=(C)y2ds    (C)x2ds=12(C)(x2+y2)ds\oint_{(C)} x^2 \mathrm{d}s = \oint_{(C)} y^2 \mathrm{d}s \implies \oint_{(C)} x^2 \mathrm{d}s = \frac{1}{2} \oint_{(C)} (x^2+y^2) \mathrm{d}s

  3. 由于在圆周上 x2+y2=1x^2+y^2 = 1 恒成立,原积分直接化简为:

     (C)x2ds=12(C)1ds=12圆周长=122π=π\oint_{(C)} x^2 \mathrm{d}s = \frac{1}{2} \oint_{(C)} 1 \mathrm{d}s = \frac{1}{2} \cdot \text{圆周长} = \frac{1}{2} \cdot 2\pi = \pi

3. 半球面上的第一型面积分(习题 6.6 A.10.(3))#

  • 题目:计算 (S)(x+y+z)dS\iint_{(S)} (x+y+z) \mathrm{d}S,其中 (S)(S) 是上半球面 z=a2x2y2z = \sqrt{a^2-x^2-y^2}

  • 解析

  1. 对称性分析:上半球面 (S)(S) 关于 yOzyOz 面及 xOzxOz 面对称。

     * 由于被积函数中的 xxyy 是奇函数,对应的积分项为 00

       (S)xdS=0,(S)ydS=0\iint_{(S)} x \mathrm{d}S = 0, \quad \iint_{(S)} y \mathrm{d}S = 0

     * 积分式直接简化为求:(S)zdS\iint_{(S)} z \mathrm{d}S

  2. 显式公式计算:上半球面在 xOyxOy 面的投影区域 (σ)(\sigma) 为圆域 x2+y2a2x^2+y^2 \le a^2

     计算偏导数:zx=xzz_x = -\frac{x}{z}zy=yzz_y = -\frac{y}{z},则:

     dS=1+zx2+zy2dxdy=1+x2+y2z2dxdy=azdxdy\mathrm{d}S = \sqrt{1 + z_x^2 + z_y^2} \mathrm{d}x\mathrm{d}y = \sqrt{1 + \frac{x^2+y^2}{z^2}} \mathrm{d}x\mathrm{d}y = \frac{a}{z} \mathrm{d}x\mathrm{d}y

  3. 代入并消去 zz

     (S)zdS=(σ)zazdxdy=a(σ)dxdy\iint_{(S)} z \mathrm{d}S = \iint_{(\sigma)} z \cdot \frac{a}{z} \mathrm{d}x\mathrm{d}y = a \iint_{(\sigma)} \mathrm{d}x\mathrm{d}y

     上式右端的二重积分为投影圆域的面积 πa2\pi a^2,因此:

     (S)(x+y+z)dS=aπa2=πa3\iint_{(S)} (x+y+z) \mathrm{d}S = a \cdot \pi a^2 = \pi a^3


第七节 第二型线积分与面积分#

一、 第二型线积分(对坐标的曲线积分)#

1. 概念与定义(公式 7.2)#

A(M)=(P,Q,R)\vec{A}(M) = (P, Q, R) 为有界向量场,(C)(C) 为从起点 AA 到终点 BB 的有向曲线。第二型线积分定义为向量点积和的极限:

(C)A(M)ds=limd0k=1nA(Mˉk)Mk1Mk\int_{(C)} \vec{A}(M) \cdot \mathrm{d}\vec{s} = \lim_{d \to 0} \sum_{k=1}^n \vec{A}(\bar{M}_k) \cdot \mathbf{{M}_{k-1}M_k}

其直角坐标分量形式为(对坐标的曲线积分,公式 7.3):

(C)P(x,y,z)dx+Q(x,y,z)dy+R(x,y,z)dz\int_{(C)} P(x,y,z) \mathrm{d}x + Q(x,y,z) \mathrm{d}y + R(x,y,z) \mathrm{d}z

2. 参数计算公式(公式 7.4)#

设光滑有向曲线 (C)(C) 的参数方程为 x=x(t),y=y(t),z=z(t)x=x(t), y=y(t), z=z(t),其中参数 tt 从起点对应的 α\alpha 变化到终点对应的 β\beta注意:不要求 α<β\alpha < \beta),则:

(C)Pdx+Qdy+Rdz=αβ[P(x(t),y(t),z(t))x˙(t)+Q(x(t),y(t),z(t))y˙(t)+R(x(t),y(t),z(t))z˙(t)]dt\int_{(C)} P\mathrm{d}x + Q\mathrm{d}y + R\mathrm{d}z = \int_{\alpha}^{\beta} \left[ P(x(t), y(t), z(t))\dot{x}(t) + Q(x(t), y(t), z(t))\dot{y}(t) + R(x(t), y(t), z(t))\dot{z}(t) \right] \mathrm{d}t

3. 方向变号性质#

若将积分路径的方向反转(记为 C-C),积分值变号:

(C)Ads=(C)Ads\int_{(C)} \vec{A} \cdot \mathrm{d}\vec{s} = -\int_{(-C)} \vec{A} \cdot \mathrm{d}\vec{s}


二、 第二型面积分(对坐标的曲面积分)#

1. 有向曲面与流量模型#

为了定义第二型面积分,曲面必须指定侧(即确定法向量 en\vec{e}_n 的指向,如上侧/下侧、内侧/外侧)。

第二型面积分在物理上代表向量场(如流速场)穿过有向曲面的通量(Flux)。

2. 概念与定义(公式 7.8 & 7.11)#

(S)A(M)dS=(S)(Aen)dS\iint_{(S)} \vec{A}(M) \cdot \mathrm{d}\vec{S} = \iint_{(S)} (\vec{A} \cdot \vec{e}_n) \mathrm{d}S

A=(P,Q,R)\vec{A} = (P, Q, R),法单位向量 en=(cosα,cosβ,cosγ)\vec{e}_n = (\cos\alpha, \cos\beta, \cos\gamma)。利用有向投影元定义写为:

(S)Pdydz+Qdzdx+Rdxdy\iint_{(S)} P \mathrm{d}y \wedge \mathrm{d}z + Q \mathrm{d}z \wedge \mathrm{d}x + R \mathrm{d}x \wedge \mathrm{d}y

3. 显式方程下的计算法与“符号法则”(公式 7.12)#

以含有 RdxdyR \mathrm{d}x \wedge \mathrm{d}y 的单项积分为例,若曲面方程为 z=z(x,y)z=z(x,y),其在 xOyxOy 面投影为 (σxy)(\sigma_{xy}),则:

(S)R(x,y,z)dxdy=±(σxy)R[x,y,z(x,y)]dxdy\iint_{(S)} R(x,y,z) \mathrm{d}x \wedge \mathrm{d}y = \pm \iint_{(\sigma_{xy})} R[x, y, z(x,y)] \mathrm{d}x\mathrm{d}y

  • 正负号判定规则

  * 若有向曲面 (S)(S) 的法向量 n\vec{n}zz 轴正方向夹角为锐角(即取上侧),前号取 ++

  * 若夹角为钝角(即取下侧),前号取 -

  * 对 dydz\mathrm{d}y \wedge \mathrm{d}z(投影到 yOzyOz 面,看前侧 [+][+] / 后侧 [][-])和 dzdx\mathrm{d}z \wedge \mathrm{d}x(投影到 zOxzOx 面,看右侧 [+][+] / 左侧 [][-])同理适用。


三、 两类积分的桥梁与联系#

第二型积分本质上可以通过方向余弦投影,转化为第一型无方向的积分。

  • 曲线积分联系(公式 7.5 & 7.6):

  (C)Pdx+Qdy+Rdz=(C)(Pcosα+Qcosβ+Rcosγ)ds\int_{(C)} P\mathrm{d}x + Q\mathrm{d}y + R\mathrm{d}z = \int_{(C)} (P\cos\alpha + Q\cos\beta + R\cos\gamma)\mathrm{d}s

  • 曲面积分联系(公式 7.9):

  (S)Pdydz+Qdzdx+Rdxdy=(S)(Pcosα+Qcosβ+Rcosγ)dS\iint_{(S)} P \mathrm{d}y \wedge \mathrm{d}z + Q \mathrm{d}z \wedge \mathrm{d}x + R \mathrm{d}x \wedge \mathrm{d}y = \iint_{(S)} (P\cos\alpha + Q\cos\beta + R\cos\gamma)\mathrm{d}S


四、 重点习题详细解答(打勾/标记题目)#

1. 沿椭圆的第二型线积分(习题 6.7 A.2.(3))#

  • 题目:计算 (C)ydxxdy\oint_{(C)} y \mathrm{d}x - x \mathrm{d}y,其中 (C)(C) 为正向椭圆 x2a2+y2b2=1\frac{x^2}{a^2} + \frac{y^2}{b^2} = 1

  • 解析

  1. 引入参数方程:x=acost,y=bsintx = a\cos t, \, y = b\sin t。由于是正向(逆时针),参数 tt00 变化到 2π2\pi

  2. 微元代换:dx=asintdt,dy=bcostdt\mathrm{d}x = -a\sin t \mathrm{d}t, \, \mathrm{d}y = b\cos t \mathrm{d}t

  3. 代入计算:

     (C)ydxxdy=02π[bsint(asint)acost(bcost)]dt\oint_{(C)} y \mathrm{d}x - x \mathrm{d}y = \int_{0}^{2\pi} \left[ b\sin t (-a\sin t) - a\cos t (b\cos t) \right] \mathrm{d}t

     =02πab(sin2t+cos2t)dt=ab02πdt=2πab= \int_{0}^{2\pi} -ab (\sin^2 t + \cos^2 t) \mathrm{d}t = -ab \int_{0}^{2\pi} \mathrm{d}t = -2\pi ab

2. 半球面上的第二型面积分(习题 6.7 A.10.(1))#

  • 题目:计算 (S)(x+1)2dxdy\iint_{(S)} (x+1)^2 \mathrm{d}x \wedge \mathrm{d}y,其中 (S)(S) 为球面 x2+y2+z2=R2x^2+y^2+z^2=R^2 (z0z \ge 0) 的上侧。

  • 解析

  1. 符号判定:曲面为上半球面且取上侧,法向量与 zz 轴夹角为锐角,因此投影到 xOyxOy 面时,系数取 ++

  2. 投影区域:投影区域 (σxy)(\sigma_{xy}) 为圆域 x2+y2R2x^2+y^2 \le R^2

  3. 转换计算

     (S)(x+1)2dxdy=+(σxy)(x2+2x+1)dxdy\iint_{(S)} (x+1)^2 \mathrm{d}x \wedge \mathrm{d}y = + \iint_{(\sigma_{xy})} (x^2+2x+1) \mathrm{d}x\mathrm{d}y

  4. 利用对称性化简:由于积分域 (σxy)(\sigma_{xy}) 关于 yy 轴对称,奇次方项 2x2x 的积分为 00

     (σxy)(x2+2x+1)dxdy=(σxy)x2dxdy+(σxy)1dxdy\iint_{(\sigma_{xy})} (x^2+2x+1) \mathrm{d}x\mathrm{d}y = \iint_{(\sigma_{xy})} x^2 \mathrm{d}x\mathrm{d}y + \iint_{(\sigma_{xy})} 1 \mathrm{d}x\mathrm{d}y

     * 第一项:由极轴对称性,x2dxdy=12(x2+y2)dxdy=1202πdθ0Rρ3dρ=πR44\iint x^2 \mathrm{d}x\mathrm{d}y = \frac{1}{2}\iint (x^2+y^2) \mathrm{d}x\mathrm{d}y = \frac{1}{2} \int_0^{2\pi}\mathrm{d}\theta\int_0^R \rho^3\mathrm{d}\rho = \frac{\pi R^4}{4}

     * 第二项:等于圆域面积 πR2\pi R^2

  * 结论

     (S)(x+1)2dxdy=πR44+πR2=πR2(1+R24)\iint_{(S)} (x+1)^2 \mathrm{d}x \wedge \mathrm{d}y = \frac{\pi R^4}{4} + \pi R^2 = \pi R^2 \left( 1 + \frac{R^2}{4} \right)

3. 旋转抛物面与平面的通量流量(习题 6.7 A.13.(2) / 对应例 7.9)#

  • 题目:求向量 r=(x,y,z)\vec{r} = (x,y,z) 穿过圆柱体 x2+y2a2,0zhx^2+y^2 \le a^2, \, 0 \le z \le h 的全表面外侧的通量。

  • 解析

  圆柱全表面外侧 (S)(S) 可以分为三部分:侧面 (S1)(S_1),底面 (S2)(S_2),顶面 (S3)(S_3)

  * 侧面 (S1)(S_1):在柱面上,径向向量 r\vec{r} 的水平分量与柱面的外法线单位向量 en\vec{e}_n 平行且同向,而竖直分量 zkz\vec{k} 与侧法线垂直(点积为 00)。因此在侧面上:

    ren=a\vec{r} \cdot \vec{e}_n = a

    (S1)rdS=(S1)adS=a(侧面积)=a2πah=2πa2h\iint_{(S_1)} \vec{r} \cdot \mathrm{d}\vec{S} = \iint_{(S_1)} a \mathrm{d}S = a \cdot (\text{侧面积}) = a \cdot 2\pi a h = 2\pi a^2 h

  * 底面 (S2)(S_2):方程为 z=0z=0,外法线指向下,即 en=k\vec{e}_n = -\vec{k}

    ren=(x,y,0)(0,0,1)=0    (S2)rdS=0\vec{r} \cdot \vec{e}_n = (x, y, 0) \cdot (0, 0, -1) = 0 \implies \iint_{(S_2)} \vec{r} \cdot \mathrm{d}\vec{S} = 0

  * 顶面 (S3)(S_3):方程为 z=hz=h,外法线指向上,即 en=+k\vec{e}_n = +\vec{k}

    ren=(x,y,h)(0,0,1)=h\vec{r} \cdot \vec{e}_n = (x, y, h) \cdot (0, 0, 1) = h

    (S3)rdS=(S3)hdS=h(顶圆面积)=hπa2=πa2h\iint_{(S_3)} \vec{r} \cdot \mathrm{d}\vec{S} = \iint_{(S_3)} h \mathrm{d}S = h \cdot (\text{顶圆面积}) = h \cdot \pi a^2 = \pi a^2 h

  * 总通量 Φ\Phi

    Φ=2πa2h+0+πa2h=3πa2h\Phi = 2\pi a^2 h + 0 + \pi a^2 h = 3\pi a^2 h


第八节 各种积分的联系及其在场论中的应用#

一、 平面格林公式(Green’s Theorem)#

格林公式是连接平面闭区域上的二重积分与沿其边界的第二型线积分的桥梁。

1. 平面区域的连通性#

  • 单连通域:直观上,区域内部没有“洞”。即区域内任何简单闭曲线围成的部分都完全属于该区域。

  • 复连通域:区域内含有“洞”(如图 6.61 所示)。

2. 格林公式定理(定理 8.1)#

设平面有界闭区域 (σ)(\sigma) 由分段光滑的简单闭曲线 (C)(C) 围成。约定 (+C)(+C)(C)(C)正向(即当人沿此方向行走时,区域 (σ)(\sigma) 始终在人的左侧)。

若函数 P(x,y),Q(x,y)P(x,y), Q(x,y)(σ)(\sigma) 上具有连续的一阶偏导数,则:

(σ)(QxPy)dσ=(+C)Pdx+Qdy\iint_{(\sigma)} \left( \frac{\partial Q}{\partial x} - \frac{\partial P}{\partial y} \right) \mathrm{d}\sigma = \oint_{(+C)} P \mathrm{d}x + Q \mathrm{d}y (公式 8.1)

3. 复连通域下的推广#

格林公式对复连通域同样成立。此时,边界曲线 (+C)(+C) 包括:

  • 外边界:逆时针方向;

  • 内边界:顺时针方向。

通过引入辅助切割线(如图 6.64 ),切割线两侧的线积分因方向相反相互抵消,从而证明了格林公式在复连通域上的普适性。


二、 平面曲线积分与路径无关的等价条件#

在单连通域内,力场做功与路径无关,在数学上对应于全微分形式。

1. 四大等价命题(定理 8.2)#

(σ)(\sigma) 是平面单连通域, P,QP, Q(σ)(\sigma) 内连续。以下四个命题相互等价:

  1. 沿 (σ)(\sigma) 内任意分段光滑的闭曲线 (C)(C),其第二型线积分为 00

   (C)Pdx+Qdy=0(环量为 0,即无旋场)\oint_{(C)} P \mathrm{d}x + Q \mathrm{d}y = 0 \quad \text{(环量为 0,即无旋场)}

  1. 沿 (σ)(\sigma) 内任意曲线的积分值仅与起点 AA 和终点 BB 有关,与积分路径无关:

   (C)Pdx+Qdy=ABPdx+Qdy(保守力场做功)\int_{(C)} P \mathrm{d}x + Q \mathrm{d}y = \int_{A}^{B} P \mathrm{d}x + Q \mathrm{d}y \quad \text{(保守力场做功)}

  1. 被积表达式 Pdx+QdyP \mathrm{d}x + Q \mathrm{d}y 是某个二元函数 u(x,y)u(x,y) 的全微分:

   du=Pdx+Qdy(势函数存在,即有势场)\mathrm{d}u = P \mathrm{d}x + Q \mathrm{d}y \quad \text{(势函数存在,即有势场)}

  1. (σ)(\sigma) 内处处满足偏导数相等条件(定理 8.3):

   PyQx\frac{\partial P}{\partial y} \equiv \frac{\partial Q}{\partial x}

2. 势函数 u(x,y)u(x,y) 的求法#

若已知 Py=Qx\frac{\partial P}{\partial y} = \frac{\partial Q}{\partial x},求势函数(全微分原函数)主要有以下三种方法:

  • 折线积分法:选择基准点 (x0,y0)(x_0, y_0),沿平行于坐标轴的折线路径进行积分:

  u(x,y)=x0xP(t,y0)dt+y0yQ(x,t)dt+Cu(x,y) = \int_{x_0}^{x} P(t, y_0) \mathrm{d}t + \int_{y_0}^{y} Q(x, t) \mathrm{d}t + C

  • 偏积分法:先对 xx 求偏积分,再对 yy 求偏导数并与 Q(x,y)Q(x,y) 进行对比:

  u(x,y)=P(x,y)dx+φ(y)    uy=Q(x,y)    求出 φ(y)u(x,y) = \int P(x,y) \mathrm{d}x + \varphi(y) \implies \frac{\partial u}{\partial y} = Q(x,y) \implies \text{求出 } \varphi(y)

  • 凑微分法:利用微分运算法则,直接拼凑出 d()\mathrm{d}(\dots) 的形式(如例 8.4 方法三)。

三、 高斯公式与散度(Gauss’s Theorem & Divergence)#

高斯公式是格林公式在三维空间中的推广,它建立了空间闭区域上的三重积分与**其边界闭曲面上的第二型面积分(通量)**的对应关系。

1. 高斯公式定理(定理 8.4)#

设空间有界闭区域 (V)(V) 的边界为分片光滑的闭曲面 (S)(S)约定法向量朝外)。若 P,Q,RP, Q, R(V)(V) 上具有连续的一阶偏导数,则:

(V)(Px+Qy+Rz)dV=(S)Pdydz+Qdzdx+Rdxdy\iiint_{(V)} \left( \frac{\partial P}{\partial x} + \frac{\partial Q}{\partial y} + \frac{\partial R}{\partial z} \right) \mathrm{d}V = \oiint_{(S)} P \mathrm{d}y \wedge \mathrm{d}z + Q \mathrm{d}z \wedge \mathrm{d}x + R \mathrm{d}x \wedge \mathrm{d}y (公式 8.7)

  • 向量形式(利用哈密顿算子 \nabla,公式 8.10):

  (V)AdV=(S)AdS\iiint_{(V)} \nabla \cdot \vec{A} \mathrm{d}V = \oiint_{(S)} \vec{A} \cdot \mathrm{d}\vec{S}

2. 散度(Divergence)#

向量场 A(M)=(P,Q,R)\vec{A}(M) = (P, Q, R) 在点 MM 处的散度定义为:

divA=A=Px+Qy+Rz\operatorname{div} \vec{A} = \nabla \cdot \vec{A} = \frac{\partial P}{\partial x} + \frac{\partial Q}{\partial y} + \frac{\partial R}{\partial z}

  • 物理意义:散度代表单位体积内流体流出的净通量,即源强密度

  * divA>0\operatorname{div} \vec{A} > 0:该点有正源(有流体产生);

  * divA<0\operatorname{div} \vec{A} < 0:该点有负源(有流体汇集/消亡)。

四、 斯托克斯公式与旋度(Stokes’ Theorem & Curl)#

斯托克斯公式将空间闭合曲线上的第二型线积分以该曲线为边界的任意有向曲面上的第二型面积分联系起来,是格林公式在三维空间中的推广。

1. 斯托克斯公式定理(定理 8.5)#

(S)(S) 为分片光滑的有向曲面,其边界为分段光滑的简单闭曲线 (C)(C)(C)(C) 的方向与 (S)(S) 的法向量 en\vec{e}_n 符合右手螺旋法则(如图 6.72 )。

P,Q,RP, Q, R 在包含 (S)(S) 的区域上具有连续的一阶偏导数,则:

(C)Pdx+Qdy+Rdz=(S)(RyQz)dydz+(PzRx)dzdx+(QxPy)dxdy\oint_{(C)} P \mathrm{d}x + Q \mathrm{d}y + R \mathrm{d}z = \iint_{(S)} \left( \frac{\partial R}{\partial y} - \frac{\partial Q}{\partial z} \right) \mathrm{d}y \mathrm{d}z + \left( \frac{\partial P}{\partial z} - \frac{\partial R}{\partial x} \right) \mathrm{d}z \mathrm{d}x + \left( \frac{\partial Q}{\partial x} - \frac{\partial P}{\partial y} \right) \mathrm{d}x \mathrm{d}y (公式 8.18)

  • 向量形式(公式 8.19 & 8.26):

  (C)Ads=(S)(×A)dS=(S)rotAdS\oint_{(C)} \vec{A} \cdot \mathrm{d}\vec{s} = \iint_{(S)} (\nabla \times \vec{A}) \cdot \mathrm{d}\vec{S} = \iint_{(S)} \operatorname{rot} \vec{A} \cdot \mathrm{d}\vec{S}

2. 环量与环量密度#

  • 环量(Circulation):第二型闭合线积分 (C)Ads\oint_{(C)} \vec{A} \cdot \mathrm{d}\vec{s} 刻画了向量场沿曲线的旋转趋势。

  • 环量密度(Circulation Density):在点 MM 处沿法向 en\vec{e}_n 的单位面积环量极限。

旋度是一个向量,其方向代表环量密度最大的方向(即旋转轴方向),其模代表最大环量密度。

3. 旋度的坐标计算公式(公式 8.25)#

rotA=×A=ijkxyzPQR=(RyQz)i+(PzRx)j+(QxPy)k\operatorname{rot} \vec{A} = \nabla \times \vec{A} = \begin{vmatrix} \vec{i} & \vec{j} & \vec{k} \\ \frac{\partial}{\partial x} & \frac{\partial}{\partial y} & \frac{\partial}{\partial z} \\ P & Q & R \end{vmatrix} = \left( \frac{\partial R}{\partial y} - \frac{\partial Q}{\partial z} \right) \vec{i} + \left( \frac{\partial P}{\partial z} - \frac{\partial R}{\partial x} \right) \vec{j} + \left( \frac{\partial Q}{\partial x} - \frac{\partial P}{\partial y} \right) \vec{k}

4. 核心场论恒等式#

  • 梯度的旋度恒为零×(u)=0(即势场无旋)\nabla \times (\nabla u) = \vec{0} \quad \text{(即势场无旋)}

  • 旋度的散度恒为零(×A)=0(即旋度场无源)\nabla \cdot (\nabla \times \vec{A}) = 0 \quad \text{(即旋度场无源)}

  • 拉普拉斯恒等式×(×A)=(A)ΔA\nabla \times (\nabla \times \vec{A}) = \nabla(\nabla \cdot \vec{A}) - \Delta \vec{A}


五、 三种重要的特殊向量场#

根据散度和旋度的不同组合,空间向量场可以分为以下三类:

1. 无旋场(Irrotational Field / 保守场 / 有势场)#

  • 定义:满足 rotA0\operatorname{rot} \vec{A} \equiv \vec{0}

  • 等价性质(一维单连通域下):

  线积分与路径无关,环量恒为 00,且存在标量势函数 uu 使得 A=u\vec{A} = \nabla u

2. 无源场(Solenoidal Field / 螺线场 / 管状场)#

  • 定义:满足 divA0\operatorname{div} \vec{A} \equiv 0

  • 等价性质(二维单连通域下):

  任意闭合曲面通量恒为 00,且存在向量势 B\vec{B} 使得 A=×B\vec{A} = \nabla \times \vec{B}。在向量管中满足通量守恒(强度不变)。

3. 调和场(Harmonic Field)#

  • 定义:既无源又无旋的场,即满足:

  divA=0,rotA=0\operatorname{div} \vec{A} = 0, \quad \operatorname{rot} \vec{A} = \vec{0}

  • 数学特征:由于其无旋,必有 A=u\vec{A} = \nabla u;代入无源条件,其势函数 uu 必满足 Laplace(拉普拉斯)方程

  Δu=2ux2+2uy2+2uz2=0\Delta u = \frac{\partial^2 u}{\partial x^2} + \frac{\partial^2 u}{\partial y^2} + \frac{\partial^2 u}{\partial z^2} = 0


六、 重点习题与例题精选分析#

1. 奇点隔离法(挖洞法)—— 处理分母为 0 的奇点问题(例 8.3 & 例 8.9)#

当被积函数在区域内的某点(奇点)无定义时,无法直接套用格林公式或高斯公式。

  • 解决方法:以该奇点为圆心(或球心),取一极小半径 ϵ\epsilon,做一小圆周(或小球面)将其隔离。在排除奇点后的复连通区域上应用定理。

  • 经典结论

  * 平面涡旋源积分:若 (C)(C) 包围原点,则 (+C)xdyydxx2+y2=2π\oint_{(+C)} \frac{x\mathrm{d}y-y\mathrm{d}x}{x^2+y^2} = 2\pi

  * 库仑电场通量(高斯定理):若空间闭曲面 (S)(S) 包围点电荷 qq,则通过该闭曲面的位移通量恒为:

    (S)DdS=q\oiint_{(S)} \vec{D} \cdot \mathrm{d}\vec{S} = q

2. 非闭合曲线与曲面的“补线法”与“补面法”#

格林公式和高斯公式只适用于闭合曲线或曲面。若给定的积分路径不闭合,需通过“补线”或“补面”使其闭合,应用公式计算后,再扣除所补部分的积分值(需特别注意边界方向)。

3. 非闭合曲面积分应用高斯公式(例 8.7 (2))#

  • 题目:计算 I=(S)x3dydz+y3dzdx+(z3+x2+y2)dxdyI = \iint_{(S)} x^3 \mathrm{d}y \wedge \mathrm{d}z + y^3 \mathrm{d}z \wedge \mathrm{d}x + (z^3 + x^2 + y^2) \mathrm{d}x \wedge \mathrm{d}y,其中 (S)(S) 是上半球面 z=R2x2y2z = \sqrt{R^2-x^2-y^2} 的上侧。

  • 解析

  1. 上半球面 (S)(S) 开放,需补面。

  2. 补面:在 xOyxOy 平面上补平面圆盘 (S1):x2+y2R2,z=0(S_1): x^2+y^2 \le R^2, z=0。为了与上半球组成闭合曲面外侧, (S1)(S_1) 的法向量必须朝下(即 cosγ=1\cos\gamma = -1)。

  3. 对闭合曲面 (SS1)(S \cup S_1) 应用高斯公式,围成的半球体记为 (V1)(V_1)

     (SS1)=(V1)((x3)x+(y3)y+(z3+x2+y2)z)dV\oiint_{(S \cup S_1)} = \iiint_{(V_1)} \left( \frac{\partial(x^3)}{\partial x} + \frac{\partial(y^3)}{\partial y} + \frac{\partial(z^3+x^2+y^2)}{\partial z} \right) \mathrm{d}V

     =(V1)3(x2+y2+z2)dV=65πR5(利用球坐标易得)= \iiint_{(V_1)} 3(x^2+y^2+z^2) \mathrm{d}V = \frac{6}{5}\pi R^5 \quad \text{(利用球坐标易得)}

  4. 计算在补偿面 (S1)(S_1) 上的积分:

     在 (S1)(S_1) 上, z=0    dz=0z=0 \implies \mathrm{d}z = 0,前两项投影微元含 dz\mathrm{d}z 故积分为 00。第三项中由于曲面取下侧,需加负号:

     (S1)=(σ1)(0+x2+y2)dxdy=02πdθ0Rρ3dρ=πR42\iint_{(S_1)} = - \iint_{(\sigma_1)} (0 + x^2+y^2) \mathrm{d}x\mathrm{d}y = - \int_{0}^{2\pi}\mathrm{d}\theta\int_{0}^{R}\rho^3\mathrm{d}\rho = -\frac{\pi R^4}{2}

  5. 还原原非闭合曲面值:

     (S)=(SS1)(S1)=65πR5(πR42)=65πR5+πR42\iint_{(S)} = \oiint_{(S \cup S_1)} - \iint_{(S_1)} = \frac{6}{5}\pi R^5 - \left( -\frac{\pi R^4}{2} \right) = \frac{6}{5}\pi R^5 + \frac{\pi R^4}{2}

4. 利用路径无关性化简复杂积分(例 8.5)#

  • 题目:计算线积分 (C)cos(x+y2)dx+[2ycos(x+y2)11+y4]dy\int_{(C)} \cos(x+y^2)\mathrm{d}x + \left[ 2y\cos(x+y^2) - \frac{1}{\sqrt{1+y^4}} \right] \mathrm{d}y,其中 (C)(C) 为摆线自 O(0,0)O(0,0)A(2πa,0)A(2\pi a, 0) 的一段。

  • 解析

  1. 验证偏导数:

     Py=y[cos(x+y2)]=2ysin(x+y2)\frac{\partial P}{\partial y} = \frac{\partial}{\partial y}[\cos(x+y^2)] = -2y\sin(x+y^2)

     Qx=x[2ycos(x+y2)11+y4]=2ysin(x+y2)\frac{\partial Q}{\partial x} = \frac{\partial}{\partial x}\left[2y\cos(x+y^2) - \frac{1}{\sqrt{1+y^4}}\right] = -2y\sin(x+y^2)

     由于偏导数处处相等,该线积分与路径无关。

  2. 改道简化:用更简单的直线段 OAOA 代替复杂的摆线 (C)(C)

     在直线段 OAOA 上,y=0    dy=0y=0 \implies \mathrm{d}y=0,自变量 xx00 变到 2πa2\pi a

  3. 代入计算:

     (C)=OA=02πacos(x+0)dx=sin(2πa)\int_{(C)} = \int_{OA} = \int_{0}^{2\pi a} \cos(x+0) \mathrm{d}x = \sin(2\pi a)

5. 星形线围成图形的面积(习题 6.8 A.3)#

  • 题目:利用线积分计算星形线 x2/3+y2/3=a2/3x^{2/3} + y^{2/3} = a^{2/3} 所围图形的面积。

  • 解析

  1. 引入星形线的参数方程:x=acos3t,y=asin3t(t[0,2π])x = a\cos^3 t, \, y = a\sin^3 t \quad (t \in [0, 2\pi])

  2. 微元计算:

     dx=3acos2tsintdt,dy=3asin2tcostdt\mathrm{d}x = -3a\cos^2 t \sin t \mathrm{d}t, \quad \mathrm{d}y = 3a\sin^2 t \cos t \mathrm{d}t

  3. 应用平面面积的格林公式线积分表示形式:

     A=12(+C)xdyydxA = \frac{1}{2} \oint_{(+C)} x \mathrm{d}y - y \mathrm{d}x

     =1202π[acos3t(3asin2tcost)asin3t(3acos2tsint)]dt= \frac{1}{2} \int_{0}^{2\pi} \left[ a\cos^3 t (3a\sin^2 t \cos t) - a\sin^3 t (-3a\cos^2 t \sin t) \right] \mathrm{d}t

     =32a202πcos2tsin2t(cos2t+sin2t)dt=38a202πsin2(2t)dt=38πa2= \frac{3}{2} a^2 \int_{0}^{2\pi} \cos^2 t \sin^2 t (\cos^2 t + \sin^2 t) \mathrm{d}t = \frac{3}{8} a^2 \int_{0}^{2\pi} \sin^2(2t) \mathrm{d}t = \frac{3}{8}\pi a^2

6. 利用斯托克斯公式求空间闭曲线积分(习题 6.8 A.10 / 类似例题)#

  • 题目:应用 Stokes 公式计算线积分 (C)ydx+zdy+xdz\oint_{(C)} y \mathrm{d}x + z \mathrm{d}y + x \mathrm{d}z,其中 (C)(C) 为圆周 x2+y2+z2=a2,x+y+z=0x^2+y^2+z^2=a^2, \, x+y+z=0。从 zz 轴正向看去为逆时针方向。

  • 解析

  1. 计算旋度:设 A=(y,z,x)\vec{A} = (y, z, x)

     rotA=ijkxyzyzx=(1)i+(1)j+(1)k=(1,1,1)\operatorname{rot} \vec{A} = \begin{vmatrix} \vec{i} & \vec{j} & \vec{k} \\ \partial_x & \partial_y & \partial_z \\ y & z & x \end{vmatrix} = (-1)\vec{i} + (-1)\vec{j} + (-1)\vec{k} = (-1, -1, -1)

  2. 确定法向量:圆周 (C)(C) 位于平面 x+y+z=0x+y+z=0 上。根据右手螺旋法则,当曲线 (C)(C) 取逆时针方向时,平面的法单位向量指向第一卦限,即:

     en=13(1,1,1)\vec{e}_n = \frac{1}{\sqrt{3}}(1, 1, 1)

  3. 应用斯托克斯公式

     I=(S)rotAendS=(S)(1,1,1)13(1,1,1)dSI = \iint_{(S)} \operatorname{rot} \vec{A} \cdot \vec{e}_n \mathrm{d}S = \iint_{(S)} (-1, -1, -1) \cdot \frac{1}{\sqrt{3}}(1, 1, 1) \mathrm{d}S

     =(S)3dS=3S面积= \iint_{(S)} -\sqrt{3} \mathrm{d}S = -\sqrt{3} \cdot S_{\text{面积}}

     由于该圆周是大圆,其张成的圆盘面积为 πa2\pi a^2,因此:

     I=3πa2I = -\sqrt{3}\pi a^2