曲线曲面积分
第一类曲线积分
物理背景:曲线的质量计算
假设平面上有一条曲线段 LLL,其线密度(单位长度的质量)不是常数,而是曲线上点的函数 ρ(x,y)\rho(x,y)ρ(x,y)(若为空间曲线则为 ρ(x,y,z)\rho(x,y,z)ρ(x,y,z))。如何计算这条曲线的总质量?
- 若线密度 ρ\rhoρ 为常数,总质量 = 线密度 × 曲线长度(即 m=ρ⋅lm = \rho \cdot lm=ρ⋅l);
- 若线密度 ρ\rhoρ 为变量,需用“分割、近似、求和、取极限”的积分思想:
- 分割:将曲线 LLL 任意分成 nnn 个小弧段 Δs1,Δs2,…,Δsn\Delta s_1, \Delta s_2, \dots, \Delta s_nΔs1,Δs2,…,Δsn(Δsi\Delta s_iΔsi 既表示第 iii 段的弧长,也表示该弧段);
- 近似:在每个小弧段 Δsi\Delta s_iΔsi 上取任意一点 (ξi,ηi)(\xi_i, \eta_i)(ξi,ηi),用该点的线密度 ρ(ξi,ηi)\rho(\xi_i, \eta_i)ρ(ξi,ηi) 近似整个小弧段的密度,則小弧段质量近似为 ρ(ξi,ηi)⋅Δsi\rho(\xi_i, \eta_i) \cdot \Delta s_iρ(ξi,ηi)⋅Δsi;
- 求和:总质量的近似值为所有小弧段质量之和:m≈∑i=1nρ(ξi,ηi)Δsim \approx \sum_{i=1}^n \rho(\xi_i, \eta_i) \Delta s_im≈∑i=1nρ(ξi,ηi)Δsi;
- 取极限:记所有小弧段的最大长度为 λ=max{Δs1,…,Δsn}\lambda = \max\{\Delta s_1, \dots, \Delta s_n\}λ=max{Δs1,…,Δsn},当 λ→0\lambda \to 0λ→0 时,若和式的极限存在,则该极限即为曲线的总质量:
m=limλ→0∑i=1nρ(ξi,ηi)Δsim = \lim_{\lambda \to 0} \sum_{i=1}^n \rho(\xi_i, \eta_i) \Delta s_im=limλ→0∑i=1nρ(ξi,ηi)Δsi
第一类曲线积分的定义
将上述物理问题抽象为数学定义:
(1)平面曲线的第一类曲线积分
设 LLL 是平面上的光滑曲线段(切线连续变化),f(x,y)f(x,y)f(x,y) 是定义在 LLL 上的有界函数。若极限
∫Lf(x,y) ds=limλ→0∑i=1nf(ξi,ηi)Δsi\int_L f(x,y) \, ds = \lim_{\lambda \to 0} \sum_{i=1}^n f(\xi_i, \eta_i) \Delta s_i∫Lf(x,y)ds=limλ→0∑i=1nf(ξi,ηi)Δsi
存在,则称该极限为 函数 f(x,y)f(x,y)f(x,y) 在曲线 LLL 上对弧长的曲线积分(即第一类曲线积分)。
- 若 LLL 是闭合曲线(首尾相连),积分记为 ∮Lf(x,y) ds\oint_L f(x,y) \, ds∮Lf(x,y)ds;
- 若 f(x,y)=1f(x,y) = 1f(x,y)=1,则积分结果就是曲线 LLL 的长度:∫L1⋅ds=l\int_L 1 \cdot ds = l∫L1⋅ds=l(lll 为 LLL 的弧长)。
(2)空间曲线的第一类曲线积分
类似地,若 Γ\GammaΓ 是空间中的光滑曲线段,f(x,y,z)f(x,y,z)f(x,y,z) 定义在 Γ\GammaΓ 上,则其第一类曲线积分为:
∫Γf(x,y,z) ds=limλ→0∑i=1nf(ξi,ηi,ζi)Δsi\int_\Gamma f(x,y,z) \, ds = \lim_{\lambda \to 0} \sum_{i=1}^n f(\xi_i, \eta_i, \zeta_i) \Delta s_i∫Γf(x,y,z)ds=limλ→0∑i=1nf(ξi,ηi,ζi)Δsi
计算方法:转化为定积分
第一类曲线积分的核心计算思路是参数化曲线,将其转化为关于参数的定积分(本质是“弧长元素 dsdsds 的参数化”)。
1. 平面曲线 LLL 的三种参数形式
(1)参数方程形式(最通用)
设 LLL 的参数方程为:
{x=φ(t)y=ψ(t)(t∈[α,β])\begin{cases}
x = \varphi(t) \\
y = \psi(t)
\end{cases} \quad (t \in [\alpha, \beta]){x=φ(t)y=ψ(t)(t∈[α,β])
其中 φ(t),ψ(t)\varphi(t), \psi(t)φ(t),ψ(t) 在 [α,β][\alpha, \beta][α,β] 上具有连续导数,且 φ′2(t)+ψ′2(t)≠0\varphi'^2(t) + \psi'^2(t) \neq 0φ′2(t)+ψ′2(t)=0(保证曲线光滑)
-
弧长元素 dsdsds 的参数化公式:
ds=(dxdt)2+(dydt)2dt=φ′2(t)+ψ′2(t)dtds = \sqrt{\left( \frac{dx}{dt} \right)^2 + \left( \frac{dy}{dt} \right)^2} dt = \sqrt{\varphi'^2(t) + \psi'^2(t)} dtds=(dtdx)2+(dtdy)2dt=φ′2(t)+ψ′2(t)dt -
积分转化公式:
∫Lf(x,y)ds=∫αβf(φ(t),ψ(t))⋅φ′2(t)+ψ′2(t)dt\int_L f(x,y) ds = \int_\alpha^\beta f\left( \varphi(t), \psi(t) \right) \cdot \sqrt{\varphi'^2(t) + \psi'^2(t)} dt∫Lf(x,y)ds=∫αβf(φ(t),ψ(t))⋅φ′2(t)+ψ′2(t)dt
注意:定积分的上下限需满足 α<β\alpha < \betaα<β(因 ds>0ds > 0ds>0,参数 ttt 需从“小”到“大”遍历曲线)。
(2)直角坐标形式(y=y(x)y = y(x)y=y(x))
若 LLL 由 y=y(x)y = y(x)y=y(x)(x∈[a,b]x \in [a, b]x∈[a,b])表示,且 y(x)y(x)y(x) 在 [a,b][a, b][a,b] 上连续可导,则可将 xxx 视为参数(t=xt = xt=x),参数方程为:
{x=xy=y(x)(x∈[a,b])\begin{cases}
x = x \\
y = y(x)
\end{cases} \quad (x \in [a, b]){x=xy=y(x)(x∈[a,b])
-
弧长元素 ds=1+(dydx)2dxds = \sqrt{1 + \left( \frac{dy}{dx} \right)^2} dxds=1+(dxdy)2dx
-
积分公式:
∫Lf(x,y)ds=∫abf(x,y(x))⋅1+y′2(x)dx(a<b)\int_L f(x,y) ds = \int_a^b f\left( x, y(x) \right) \cdot \sqrt{1 + y'^2(x)} dx \quad (a < b)∫Lf(x,y)ds=∫abf(x,y(x))⋅1+y′2(x)dx(a<b)
(3)极坐标形式(r=r(θ)r = r(\theta)r=r(θ))
若 LLL 由极坐标方程 r=r(θ)r = r(\theta)r=r(θ)(θ∈[θ1,θ2]\theta \in [\theta_1, \theta_2]θ∈[θ1,θ2])表示,且 r(θ)r(\theta)r(θ) 在 [θ1,θ2][\theta_1, \theta_2][θ1,θ2] 上连续可导,利用极坐标与直角坐标的关系 x=r(θ)cosθx = r(\theta)\cos\thetax=r(θ)cosθ,y=r(θ)sinθy = r(\theta)\sin\thetay=r(θ)sinθ,参数方程为:
{x=r(θ)cosθy=r(θ)sinθ(θ∈[θ1,θ2])\begin{cases}
x = r(\theta)\cos\theta \\
y = r(\theta)\sin\theta
\end{cases} \quad (\theta \in [\theta_1, \theta_2]){x=r(θ)cosθy=r(θ)sinθ(θ∈[θ1,θ2])
-
弧长元素 ds=(dxdθ)2+(dydθ)2dθ=r2(θ)+r′2(θ)dθds = \sqrt{\left( \frac{dx}{d\theta} \right)^2 + \left( \frac{dy}{d\theta} \right)^2} d\theta = \sqrt{r^2(\theta) + r'^2(\theta)} d\thetads=(dθdx)2+(dθdy)2dθ=r2(θ)+r′2(θ)dθ
-
积分公式:
∫Lf(x,y)ds=∫θ1θ2f(r(θ)cosθ,r(θ)sinθ)⋅r2(θ)+r′2(θ)dθ(θ1<θ2)\int_L f(x,y) ds = \int_{\theta_1}^{\theta_2} f\left( r(\theta)\cos\theta, r(\theta)\sin\theta \right) \cdot \sqrt{r^2(\theta) + r'^2(\theta)} d\theta \quad (\theta_1 < \theta_2)∫Lf(x,y)ds=∫θ1θ2f(r(θ)cosθ,r(θ)sinθ)⋅r2(θ)+r′2(θ)dθ(θ1<θ2)
2. 空间曲线 Γ\GammaΓ 的参数形式
设 Γ\GammaΓ 的参数方程为:
{x=φ(t)y=ψ(t)z=ω(t)(t∈[α,β])\begin{cases}
x = \varphi(t) \\
y = \psi(t) \\
z = \omega(t)
\end{cases} \quad (t \in [\alpha, \beta])⎩⎨⎧x=φ(t)y=ψ(t)z=ω(t)(t∈[α,β])
其中 φ(t),ψ(t),ω(t)\varphi(t), \psi(t), \omega(t)φ(t),ψ(t),ω(t) 在 [α,β][\alpha, \beta][α,β] 上连续可导,且 φ′2(t)+ψ′2(t)+ω′2(t)≠0\varphi'^2(t) + \psi'^2(t) + \omega'^2(t) \neq 0φ′2(t)+ψ′2(t)+ω′2(t)=0。
-
弧长元素 ds=φ′2(t)+ψ′2(t)+ω′2(t)dtds = \sqrt{\varphi'^2(t) + \psi'^2(t) + \omega'^2(t)} dtds=φ′2(t)+ψ′2(t)+ω′2(t)dt
-
积分公式:
∫Γf(x,y,z)ds=∫αβf(φ(t),ψ(t),ω(t))⋅φ′2(t)+ψ′2(t)+ω′2(t)dt(α<β)\int_\Gamma f(x,y,z) ds = \int_\alpha^\beta f\left( \varphi(t), \psi(t), \omega(t) \right) \cdot \sqrt{\varphi'^2(t) + \psi'^2(t) + \omega'^2(t)} dt \quad (\alpha < \beta)∫Γf(x,y,z)ds=∫αβf(φ(t),ψ(t),ω(t))⋅φ′2(t)+ψ′2(t)+ω′2(t)dt(α<β)
第一类曲线积分的应用场景
第一类曲线积分的核心应用是计算曲线的“分布量”,常见场景包括:
-
曲线的质量:若线密度为 ρ(x,y)\rho(x,y)ρ(x,y)(或 ρ(x,y,z)\rho(x,y,z)ρ(x,y,z)),则质量 m=∫Lρdsm = \int_L \rho dsm=∫Lρds;
-
曲线的重心(形心):
- 平面曲线的重心 (xˉ,yˉ)(\bar{x}, \bar{y})(xˉ,yˉ):xˉ=1m∫Lxρds\bar{x} = \frac{1}{m} \int_L x\rho dsxˉ=m1∫Lxρds,yˉ=1m∫Lyρds\bar{y} = \frac{1}{m} \int_L y\rho dsyˉ=m1∫Lyρds;
- 若 ρ=1\rho = 1ρ=1(均匀曲线),则为形心,此时 m=lm = lm=l(曲线长度)。
-
曲线的转动惯量:
- 对 xxx 轴的转动惯量 Ix=∫Ly2ρdsI_x = \int_L y^2 \rho dsIx=∫Ly2ρds;
- 对原点的转动惯量 IO=∫L(x2+y2)ρdsI_O = \int_L (x^2 + y^2) \rho dsIO=∫L(x2+y2)ρds。
第二类曲线积分
物理背景
第二类曲线积分(也称为“对坐标的曲线积分”)起源于对变力沿曲线做功的计算。假设一个质点在变力F⃗(x,y,z)=P(x,y,z)i⃗+Q(x,y,z)j⃗+R(x,y,z)k⃗\vec{F}(x, y, z) = P(x, y, z)\vec{i} + Q(x, y, z)\vec{j} + R(x, y, z)\vec{k}F(x,y,z)=P(x,y,z)i+Q(x,y,z)j+R(x,y,z)k的作用下,沿空间曲线LLL从点AAA移动到点BBB,则变力所做的功可通过第二类曲线积分表示。
数学定义和表达式
设曲线LLL是从点AAA到点BBB的有向曲线段,函数P(x,y)P(x,y)P(x,y)、Q(x,y)Q(x,y)Q(x,y)在LLL上有界,则第二类曲线积分定义为:
-
平面曲线积分:
∫LP(x,y)dx+Q(x,y)dy\int_L P(x,y)dx + Q(x,y)dy∫LP(x,y)dx+Q(x,y)dy
其中,dxdxdx和dydydy分别是曲线LLL在xxx轴和yyy轴上的有向投影。 -
空间曲线积分:
∫LP(x,y,z)dx+Q(x,y,z)dy+R(x,y,z)dz\int_L P(x,y,z)dx + Q(x,y,z)dy + R(x,y,z)dz∫LP(x,y,z)dx+Q(x,y,z)dy+R(x,y,z)dz
计算方法1:参数方程法
设曲线LLL的参数方程为:
-
平面情形:x=x(t)x = x(t)x=x(t),y=y(t)y = y(t)y=y(t),当ttt从α\alphaα变到β\betaβ时,点(x,y)(x,y)(x,y)从AAA沿LLL到BBB,则:
∫LPdx+Qdy=∫αβ[P(x(t),y(t))x′(t)+Q(x(t),y(t))y′(t)]dt\int_L Pdx + Qdy = \int_\alpha^\beta \left[ P(x(t), y(t))x'(t) + Q(x(t), y(t))y'(t) \right] dt∫LPdx+Qdy=∫αβ[P(x(t),y(t))x′(t)+Q(x(t),y(t))y′(t)]dt -
空间情形:x=x(t)x = x(t)x=x(t),y=y(t)y = y(t)y=y(t),z=z(t)z = z(t)z=z(t),则:
∫LPdx+Qdy+Rdz=∫αβ[Px′+Qy′+Rz′]dt\int_L Pdx + Qdy + Rdz = \int_\alpha^\beta \left[ Px' + Qy' + Rz' \right] dt∫LPdx+Qdy+Rdz=∫αβ[Px′+Qy′+Rz′]dt
计算方法2:格林公式
什么是格林公式?
设区域 D⊂R2D\subset \mathbb R^2D⊂R2 满足:
- DDD 为有界闭区域
- 边界是一条分段光滑的简单闭曲线 C=∂DC=\partial DC=∂D
- P(x,y),Q(x,y)P(x,y), Q(x,y)P(x,y),Q(x,y) 在包含 DDD 的开区域上连续可偏导
且曲线 CCC 的方向为 沿边界的逆时针方向(正向),
则:
∮C(P dx+Q dy)=∬D(∂Q∂x−∂P∂y) dA \boxed{\oint_{C} \left(P\,dx + Q\,dy\right)= \iint_{D} \left(\frac{\partial Q}{\partial x}-\frac{\partial P}{\partial y}\right)\, dA } ∮C(Pdx+Qdy)=∬D(∂x∂Q−∂y∂P)dA
格林公式的几何含义
- 左边:围着边界曲线 CCC 转一圈所做的“环量”
- 右边:区域内部的“整体旋度”的面积积分
也就是说:
边界上的环流 = 区域内的旋度之和
格林公式中曲线的正向方向(非常重要)
格林公式对方向有严格要求:
- 我们规定格林公式中曲线的正向为逆时针,也就是区域在行进方向的 左侧
- 如果实际计算的时候发现曲线的方向是 顺时针,则化为二重积分之后需要在前面多加一个负号:
∮C顺时针(Pdx+Qdy)=−∬D(Qx−Py)dA \oint_{C}^{\text{顺时针}}(Pdx+Qdy) = -\iint_D\left(Q_x - P_y\right)dA ∮C顺时针(Pdx+Qdy)=−∬D(Qx−Py)dA
用第二类曲线积分描述平面图形的面积公式
① 基本公式
取:
P=0,Q=x P=0,\quad Q=x P=0,Q=x
得:
A(D)=∬D1 dA=∮Cx dy. A(D)=\iint_D 1\, dA = \oint_C x\, dy. A(D)=∬D1dA=∮Cxdy.
或取 P=−y,Q=0P=-y,Q=0P=−y,Q=0:
A(D)=∮C−y dx. A(D)=\oint_C -y\, dx. A(D)=∮C−ydx.
② 面积的对称形式
A(D)=12∮C(x dy−y dx) A(D)=\frac12\oint_C (x\,dy - y\,dx) A(D)=21∮C(xdy−ydx)
格林公式的物理意义(旋度与环流)
格林公式是二维的 Stokes 定理:
边界上的环流=区域内部旋度之和 \text{边界上的环流}=\text{区域内部旋度之和} 边界上的环流=区域内部旋度之和
-
若 Qx−Py=0Q_x - P_y = 0Qx−Py=0,则线积分与路径无关
-
若区域单连通,则存在势函数 ϕ\phiϕ,使得
P=ϕx,Q=ϕy. P=\phi_x, \quad Q=\phi_y. P=ϕx,Q=ϕy.
⭐ 小结
格林公式:
∮C(Pdx+Qdy)=∬D(Qx−Py)dA\oint_C(Pdx+Qdy)= \iint_D(Q_x-P_y)dA ∮C(Pdx+Qdy)=∬D(Qx−Py)dA
要求:
- CCC 逆时针
- P,QP,QP,Q 一阶连续可偏导
- CCC 为简单闭曲线
- DDD 是 CCC 围成的区域
对于一个积分与路径无关的场,是不是任意闭合曲线在该场中的第二类曲线积分都是0?
当然是的,在无旋场中,积分与路径无关,只与曲线的起点和终点位置有关,因为闭合曲线的起点和终点是重合的,所以无旋场中环路曲线积分始终为0。
既然如此,请你解释解释,在下面这个例子中,通过计算我们发现∂Q∂x=∂P∂y\frac{\partial Q}{\partial x} = \frac{\partial P}{\partial y}∂x∂Q=∂y∂P,为什么这个闭合环路积分不为0呢?

原因很简单,因为我们前面说的条件中明确说了 “对于一个积分与路径无关的场”,这就意味着场中的任意一点,都要满足∂Q∂x=∂P∂y\frac{\partial Q}{\partial x} = \frac{\partial P}{\partial y}∂x∂Q=∂y∂P,而在上面的例子中,(0,0)(0,0)(0,0)这一点就不满足,为啥呢?因为在(0,0)(0,0)(0,0)这一点,∂Q∂x和∂P∂y\frac{\partial Q}{\partial x} 和 \frac{\partial P}{\partial y}∂x∂Q和∂y∂P都不存在!于是我们得到的无旋场,是不能包含(0,0)(0,0)(0,0)这一点的!
如何对一个中间有洞的区域使用格林公式?
例如下图红色区域,就是一个无旋场,下面我们就对这个区域使用格林公式

格林公式可以用于有多个边界的区域(含“空洞”的区域)
如果区域 DDD 是一个“外圈”和一个“内圈(洞)”组成的环状区域,例如:
- 外边界 C1C_1C1:上图中的星型曲线
- 内边界 C2C_2C2:上图中的小椭圆
那么格林公式仍然成立,只需在 方向上进行特殊处理。
注意这时候的曲线方向!⭐⭐ 方向规则:
| 边界类型 | 方向 |
|---|---|
| 外边界 C1C_1C1 | 逆时针(正向) |
| 内边界(洞) C2C_2C2 | 顺时针(反向) |
原因:格林公式要求“区域永远在曲线左侧”。
- 沿外边界逆时针走 → 区域在左侧
- 沿内边界必须顺时针走 → 才能让区域仍在左侧
所以,当区域有洞时,只有外圈和内圈的方向相反,才能应用格林公式。
题目给的外圈方向是顺时针,与正方向相反,为了方便处理,我们可以先算外圈方向是正方向时候的曲线积分,最后再将结果取反,就可以得到顺时针方向的曲线积分结果
为什么在无源场中可以换路径积分?
假设外圈曲线C1C_1C1的方向是逆时针,我们选取内圈C2C_2C2的曲线方向是顺时针,则∮C1P dx+Q dy − ∮C2P dx+Q dy=∬D(∂Q∂x−∂P∂y)dA. \oint_{C_1} P\,dx + Q\,dy \;-\; \oint_{C_2} P\,dx + Q\,dy= \iint_D \left(\frac{\partial Q}{\partial x}-\frac{\partial P}{\partial y}\right)dA. ∮C1Pdx+Qdy−∮C2Pdx+Qdy=∬D(∂x∂Q−∂y∂P)dA.
由于在红色区域内,∂Q∂x=∂P∂y\frac{\partial Q}{\partial x} = \frac{\partial P}{\partial y}∂x∂Q=∂y∂P始终成立,因此
∮C1P dx+Q dy − ∮C2P dx+Q dy=0
\oint_{C_1} P\,dx + Q\,dy \;-\; \oint_{C_2} P\,dx + Q\,dy=
0
∮C1Pdx+Qdy−∮C2Pdx+Qdy=0
即
∮C1P dx+Q dy = ∮C2P dx+Q dy
\oint_{C_1} P\,dx + Q\,dy \;=\; \oint_{C_2} P\,dx + Q\,dy
∮C1Pdx+Qdy=∮C2Pdx+Qdy
这就是保守场(无源场)可以换路径积分的原因!
计算方法3:通过斯托克斯公式转化成第二类曲面积分
第一类曲面积分
第一类曲面积分定义
设曲面Σ\SigmaΣ是光滑的,函数f(x,y,z)f(x,y,z)f(x,y,z)在Σ\SigmaΣ上有界。将Σ\SigmaΣ任意分成nnn个小块ΔSi\Delta S_iΔSi(ΔSi\Delta S_iΔSi同时表示第iii个小块的面积),在每个ΔSi\Delta S_iΔSi上任取一点(ξi,ηi,ζi)(\xi_i, \eta_i, \zeta_i)(ξi,ηi,ζi),作和∑i=1nf(ξi,ηi,ζi)ΔSi\sum_{i=1}^n f(\xi_i, \eta_i, \zeta_i) \Delta S_i∑i=1nf(ξi,ηi,ζi)ΔSi。如果当各小块曲面的直径的最大值λ→0\lambda \to 0λ→0时,这个和的极限存在,则称此极限为函数f(x,y,z)f(x,y,z)f(x,y,z)在曲面Σ\SigmaΣ上的第一类曲面积分(或对面积的曲面积分),记作∬Σf(x,y,z)dS\iint_\Sigma f(x,y,z) dS∬Σf(x,y,z)dS,即:
∬Σf(x,y,z)dS=limλ→0∑i=1nf(ξi,ηi,ζi)ΔSi\iint_\Sigma f(x,y,z) dS = \lim_{\lambda \to 0} \sum_{i=1}^n f(\xi_i, \eta_i, \zeta_i) \Delta S_i∬Σf(x,y,z)dS=limλ→0∑i=1nf(ξi,ηi,ζi)ΔSi
其中,f(x,y,z)f(x,y,z)f(x,y,z)称为被积函数,Σ\SigmaΣ称为积分曲面,dSdSdS称为曲面的面积元素。
第一类曲面积分物理意义
若f(x,y,z)f(x,y,z)f(x,y,z)是曲面Σ\SigmaΣ的面密度(单位面积的质量),则第一类曲面积分∬Σf(x,y,z)dS\iint_\Sigma f(x,y,z) dS∬Σf(x,y,z)dS表示曲面Σ\SigmaΣ的总质量。
第一类曲面积分的计算方法(投影转化为二重积分)
设曲面Σ\SigmaΣ的方程为z=z(x,y)z = z(x,y)z=z(x,y),Σ\SigmaΣ在xOyxOyxOy面上的投影区域为DxyD_{xy}Dxy,函数z=z(x,y)z = z(x,y)z=z(x,y)在DxyD_{xy}Dxy上具有连续偏导数,则第一类曲面积分可转化为以下二重积分:
∬Σf(x,y,z)dS=∬Dxyf(x,y,z(x,y))1+zx2(x,y)+zy2(x,y) dxdy\iint_\Sigma f(x,y,z) dS = \iint_{D_{xy}} f(x,y,z(x,y)) \sqrt{1 + z_x^2(x,y) + z_y^2(x,y)} \, dxdy∬Σf(x,y,z)dS=∬Dxyf(x,y,z(x,y))1+zx2(x,y)+zy2(x,y)dxdy
其中,zx=∂z∂xz_x = \frac{\partial z}{\partial x}zx=∂x∂z,zy=∂z∂yz_y = \frac{\partial z}{\partial y}zy=∂y∂z,1+zx2+zy2 dxdy\sqrt{1 + z_x^2 + z_y^2} \, dxdy1+zx2+zy2dxdy是曲面Σ\SigmaΣ的**面积元素dSdSdS**的表达式。
同理,若曲面Σ\SigmaΣ的方程为x=x(y,z)x = x(y,z)x=x(y,z)(投影到yOzyOzyOz面)或y=y(x,z)y = y(x,z)y=y(x,z)(投影到xOzxOzxOz面),也可类似地转化为相应的二重积分。
计算过程中的注意事项
- 第一类曲面积分是有对称性的! 我们一定要尽可能地利用积分区域的对称性,对曲面积分狠狠进行化简,如果最终里面的被积函数为常数,那么曲面积分就转化为积分曲面的表面积了
- 计算复杂的曲面积分时,一定要用好曲面表达式这一条件,通过将其代入,简化曲面积分的被积函数
第二类曲面积分
第二类曲面积分的定义
设Σ\SigmaΣ是有向光滑曲面(规定了侧,如左侧、右侧、上侧、下侧等),函数 R(x,y,z)R(x,y,z)R(x,y,z) 在Σ\SigmaΣ上有界。将Σ\SigmaΣ任意分成nnn个小块ΔSi\Delta S_iΔSi(ΔSi\Delta S_iΔSi同时表示第iii个小块的面积),ΔSi\Delta S_iΔSi在xOyxOyxOy面上的投影为(ΔSi)xy(\Delta S_i)_{xy}(ΔSi)xy(投影有正负,取决于Σ\SigmaΣ的侧),在ΔSi\Delta S_iΔSi上任取一点(ξi,ηi,ζi)(\xi_i, \eta_i, \zeta_i)(ξi,ηi,ζi),作和∑i=1nR(ξi,ηi,ζi)(ΔSi)xy\sum_{i=1}^n R(\xi_i, \eta_i, \zeta_i) (\Delta S_i)_{xy}∑i=1nR(ξi,ηi,ζi)(ΔSi)xy。若当各小块曲面的直径的最大值λ→0\lambda \to 0λ→0时,该和的极限存在,则称此极限为函数R(x,y,z)R(x,y,z)R(x,y,z)在Σ\SigmaΣ上对**x,yx,yx,y坐标的曲面积分**,记作∬ΣR(x,y,z)dxdy\iint_\Sigma R(x,y,z) dxdy∬ΣR(x,y,z)dxdy,即:
∬ΣR(x,y,z)dxdy=limλ→0∑i=1nR(ξi,ηi,ζi)(ΔSi)xy\iint_\Sigma R(x,y,z) dxdy = \lim_{\lambda \to 0} \sum_{i=1}^n R(\xi_i, \eta_i, \zeta_i) (\Delta S_i)_{xy}∬ΣR(x,y,z)dxdy=limλ→0∑i=1nR(ξi,ηi,ζi)(ΔSi)xy
同理,可定义对y,zy,zy,z坐标的曲面积分∬ΣP(x,y,z)dydz\iint_\Sigma P(x,y,z) dydz∬ΣP(x,y,z)dydz和对x,zx,zx,z坐标的曲面积分∬ΣQ(x,y,z)dxdz\iint_\Sigma Q(x,y,z) dxdz∬ΣQ(x,y,z)dxdz。
第二类曲面积分的一般形式为这三个积分的和:
∬ΣPdydz+Qdxdz+Rdxdy\iint_\Sigma P dydz + Q dxdz + R dxdy∬ΣPdydz+Qdxdz+Rdxdy
第二类曲面积分的物理意义
若v⃗(x,y,z)=Pi⃗+Qj⃗+Rk⃗\vec{v}(x,y,z) = P\vec{i} + Q\vec{j} + R\vec{k}v(x,y,z)=Pi+Qj+Rk是流速场(单位时间内流体通过单位面积的流量),则第二类曲面积分∬ΣPdydz+Qdxdz+Rdxdy\iint_\Sigma P dydz + Q dxdz + R dxdy∬ΣPdydz+Qdxdz+Rdxdy表示单位时间内流体通过有向曲面Σ\SigmaΣ的流量(通量)。
基本计算方法:投影法
设曲面Σ\SigmaΣ的方程为z=z(x,y)z = z(x,y)z=z(x,y),Σ\SigmaΣ在xOyxOyxOy面上的投影区域为DxyD_{xy}Dxy,函数z=z(x,y)z = z(x,y)z=z(x,y)在DxyD_{xy}Dxy上具有连续偏导数:
- 若Σ\SigmaΣ取上侧(法向量与zzz轴正方向夹角为锐角),则
∬ΣR(x,y,z)dxdy=∬DxyR(x,y,z(x,y))dxdy\iint_\Sigma R(x,y,z) dxdy = \iint_{D_{xy}} R(x,y,z(x,y)) dxdy∬ΣR(x,y,z)dxdy=∬DxyR(x,y,z(x,y))dxdy - 若Σ\SigmaΣ取下侧(法向量与zzz轴正方向夹角为钝角),则
∬ΣR(x,y,z)dxdy=−∬DxyR(x,y,z(x,y))dxdy\iint_\Sigma R(x,y,z) dxdy = -\iint_{D_{xy}} R(x,y,z(x,y)) dxdy∬ΣR(x,y,z)dxdy=−∬DxyR(x,y,z(x,y))dxdy
同理,对于∬ΣPdydz\iint_\Sigma P dydz∬ΣPdydz(曲面Σ:x=x(y,z)\Sigma: x = x(y,z)Σ:x=x(y,z),投影到yOzyOzyOz面,区分左侧/右侧)和∬ΣQdxdz\iint_\Sigma Q dxdz∬ΣQdxdz(曲面Σ:y=y(x,z)\Sigma: y = y(x,z)Σ:y=y(x,z),投影到xOzxOzxOz面,区分前侧/后侧),也可类似转化,注意方向对符号的影响。
特殊计算方法:将第二类曲面积分 转化成 第一类曲面积分
- 求隐函数的偏导数:计算Fx=∂F∂xF_x = \frac{\partial F}{\partial x}Fx=∂x∂F,Fy=∂F∂yF_y = \frac{\partial F}{\partial y}Fy=∂y∂F,Fz=∂F∂zF_z = \frac{\partial F}{\partial z}Fz=∂z∂F;
- 构造法向量:n⃗=(Fx,Fy,Fz)\vec{n} = (F_x, F_y, F_z)n=(Fx,Fy,Fz),定向决定符号(如取外侧时,法向量指向曲面外);
- 计算法向量模长:∣n⃗∣=Fx2+Fy2+Fz2|\vec{n}| = \sqrt{F_x^2 + F_y^2 + F_z^2}∣n∣=Fx2+Fy2+Fz2;
- 求方向余弦:cosα=Fx∣n⃗∣\cos\alpha = \frac{F_x}{|\vec{n}|}cosα=∣n∣Fx,cosβ=Fy∣n⃗∣\cos\beta = \frac{F_y}{|\vec{n}|}cosβ=∣n∣Fy,cosγ=Fz∣n⃗∣\cos\gamma = \frac{F_z}{|\vec{n}|}cosγ=∣n∣Fz;
其中cosα\cos\alphacosα是法向量与xxx轴夹角的余弦 - 利用两类曲面积分的联系转化
第二类曲面积分可通过法向量方向余弦转化为第一类曲面积分:
设Σ\SigmaΣ的单位法向量为n=(cosα,cosβ,cosγ)\boxed{n} = (\cos\alpha, \cos\beta, \cos\gamma)n=(cosα,cosβ,cosγ),则
dydz=cosαdS\mathrm{d}y\mathrm{d}z = \cos\alpha \mathrm{d}Sdydz=cosαdS,dzdx=cosβdS\mathrm{d}z\mathrm{d}x = \cos\beta \mathrm{d}Sdzdx=cosβdS,dxdy=cosγdS\mathrm{d}x\mathrm{d}y = \cos\gamma \mathrm{d}Sdxdy=cosγdS
因此:∬ΣPdydz+Qdzdx+Rdxdy=∬Σ(Pcosα+Qcosβ+Rcosγ)dS\iint_{\Sigma} P \mathrm{d}y\mathrm{d}z + Q \mathrm{d}z\mathrm{d}x + R \mathrm{d}x\mathrm{d}y = \iint_{\Sigma} (P \cos\alpha + Q \cos\beta + R \cos\gamma) \mathrm{d}S∬ΣPdydz+Qdzdx+Rdxdy=∬Σ(Pcosα+Qcosβ+Rcosγ)dS
转化为第一类曲面积分后,可利用第一类的计算方法(如投影到坐标面)求解。 - 也可以用换底公式
-
dxdzdxdy=cosβcosγ=FyFz\frac{\mathrm{d}x\mathrm{d}z}{\mathrm{d}x\mathrm{d}y} = \frac{\cos\beta}{\cos\gamma} = \frac{F_y}{F_z}dxdydxdz=cosγcosβ=FzFy
-
dydzdxdy=cosαcosγ=FxFz\frac{\mathrm{d}y\mathrm{d}z}{\mathrm{d}x\mathrm{d}y} = \frac{\cos\alpha}{\cos\gamma} = \frac{F_x}{F_z}dxdydydz=cosγcosα=FzFx
-
特殊计算方法:高斯公式——将 空间闭曲面上的第二类曲面积分 转化成 闭曲面所围区域内的三重积分
高斯公式内容
设空间闭区域Ω\OmegaΩ是由分片光滑的闭曲面Σ\SigmaΣ围成,Σ\SigmaΣ取外侧(即闭曲面的“外法线方向”一侧);函数P(x,y,z)P(x,y,z)P(x,y,z)、Q(x,y,z)Q(x,y,z)Q(x,y,z)、R(x,y,z)R(x,y,z)R(x,y,z)在Ω\OmegaΩ上具有一阶连续偏导数,则:
∯ΣP dydz+Q dxdz+R dxdy=∭Ω(∂P∂x+∂Q∂y+∂R∂z)dxdydz\oiint_\Sigma P \, dydz + Q \, dxdz + R \, dxdy = \iiint_\Omega \left( \frac{\partial P}{\partial x} + \frac{\partial Q}{\partial y} + \frac{\partial R}{\partial z} \right) dxdydz∬ΣPdydz+Qdxdz+Rdxdy=∭Ω(∂x∂P+∂y∂Q+∂z∂R)dxdydz
高斯公式的物理意义(通量与散度)
- 若v⃗=(P,Q,R)\vec{v} = (P, Q, R)v=(P,Q,R)是流速场,则左侧的第二类曲面积分表示单位时间内流体通过闭曲面Σ\SigmaΣ的总流量(通量)。
- 右侧被积函数∂P∂x+∂Q∂y+∂R∂z\frac{\partial P}{\partial x} + \frac{\partial Q}{\partial y} + \frac{\partial R}{\partial z}∂x∂P+∂y∂Q+∂z∂R称为向量场v⃗\vec{v}v的散度,记为divv⃗\text{div} \vec{v}divv。高斯公式表明:闭曲面的通量等于区域内散度的三重积分,反映了“通量的来源与散度的分布”之间的关系。
高斯公式使用条件
- 曲面Σ\SigmaΣ是闭曲面(包围一个空间区域Ω\OmegaΩ);
- Σ\SigmaΣ取外侧(若取内侧,公式需加负号:∯Σ内侧⋯=−∭Ωdivv⃗ dxdydz\oiint_{\Sigma^{\text{内侧}}} \cdots = -\iiint_\Omega \text{div} \vec{v} \, dxdydz∬Σ内侧⋯=−∭Ωdivvdxdydz);
- 函数P,Q,RP, Q, RP,Q,R在Ω\OmegaΩ上一阶偏导数连续。
特殊计算方法:斯托克斯公式——将 空间有向闭曲线上的第二类曲线积分 转化成 以该曲线为边界的有向曲面上的第二类曲面积分
设Γ\GammaΓ是分段光滑的空间有向闭曲线,Σ\SigmaΣ是以Γ\GammaΓ为边界的分片光滑的有向曲面(Γ\GammaΓ称为Σ\SigmaΣ的“边界曲线”),且Γ\GammaΓ的正向与Σ\SigmaΣ的侧满足右手规则(右手四指沿Γ\GammaΓ正向弯曲,拇指指向即为Σ\SigmaΣ的法向量方向);函数P(x,y,z)P(x,y,z)P(x,y,z)、Q(x,y,z)Q(x,y,z)Q(x,y,z)、R(x,y,z)R(x,y,z)R(x,y,z)在包含Σ\SigmaΣ在内的一个空间区域内具有一阶连续偏导数,则:
1. 展开形式
∮ΓP dx+Q dy+R dz=∬Σ(∂R∂y−∂Q∂z)dydz+(∂P∂z−∂R∂x)dxdz+(∂Q∂x−∂P∂y)dxdy\oint_\Gamma P \, dx + Q \, dy + R \, dz = \iint_\Sigma \left( \frac{\partial R}{\partial y} - \frac{\partial Q}{\partial z} \right) dydz + \left( \frac{\partial P}{\partial z} - \frac{\partial R}{\partial x} \right) dxdz + \left( \frac{\partial Q}{\partial x} - \frac{\partial P}{\partial y} \right) dxdy∮ΓPdx+Qdy+Rdz=∬Σ(∂y∂R−∂z∂Q)dydz+(∂z∂P−∂x∂R)dxdz+(∂x∂Q−∂y∂P)dxdy
2. 行列式记忆形式(更直观)
为便于记忆,公式可写成“旋度的曲面积分”形式,其中rotv⃗\text{rot} \vec{v}rotv(或∇×v⃗\nabla \times \vec{v}∇×v)表示向量场v⃗=(P,Q,R)\vec{v} = (P, Q, R)v=(P,Q,R)的旋度:
∮Γv⃗⋅dl⃗=∬Σ(rotv⃗)⋅dS⃗\oint_\Gamma \vec{v} \cdot d\vec{l} = \iint_\Sigma (\text{rot} \vec{v}) \cdot d\vec{S}∮Γv⋅dl=∬Σ(rotv)⋅dS
对应的行列式表达式为:
∮ΓPdx+Qdy+Rdz=∬Σ∣dydzdxdzdxdy∂∂x∂∂y∂∂zPQR∣\oint_\Gamma P dx + Q dy + R dz = \iint_\Sigma \begin{vmatrix}
dydz & dxdz & dxdy \\
\frac{\partial}{\partial x} & \frac{\partial}{\partial y} & \frac{\partial}{\partial z} \\
P & Q & R
\end{vmatrix}∮ΓPdx+Qdy+Rdz=∬Σdydz∂x∂Pdxdz∂y∂Qdxdy∂z∂R
其中,行列式按第一行展开即得展开形式,旋度rotv⃗=(∂R∂y−∂Q∂z,∂P∂z−∂R∂x,∂Q∂x−∂P∂y)\text{rot} \vec{v} = \left( \frac{\partial R}{\partial y} - \frac{\partial Q}{\partial z}, \frac{\partial P}{\partial z} - \frac{\partial R}{\partial x}, \frac{\partial Q}{\partial x} - \frac{\partial P}{\partial y} \right)rotv=(∂y∂R−∂z∂Q,∂z∂P−∂x∂R,∂x∂Q−∂y∂P)。
斯托克斯公式的物理意义(环流与旋度)
- 左侧的第二类曲线积分∮Γv⃗⋅dl⃗\oint_\Gamma \vec{v} \cdot d\vec{l}∮Γv⋅dl称为向量场v⃗\vec{v}v沿闭曲线Γ\GammaΓ的环流,表示“向量场推动质点沿Γ\GammaΓ运动一周所做的功”(若v⃗\vec{v}v是力场),或“流体沿Γ\GammaΓ的旋转趋势强度”(若v⃗\vec{v}v是流速场)。
- 右侧的曲面积分是旋度rotv⃗\text{rot} \vec{v}rotv在曲面Σ\SigmaΣ上的通量,旋度rotv⃗\text{rot} \vec{v}rotv描述了“向量场在某点的局部旋转强度与方向”。
斯托克斯公式表明:向量场沿闭曲线的环流,等于其旋度通过以该曲线为边界的任意曲面的通量,反映了“整体环流”与“局部旋度”的内在联系。
斯托克斯公式的使用条件
- 曲线与曲面的关联:Γ\GammaΓ是Σ\SigmaΣ的边界曲线(Σ\SigmaΣ以Γ\GammaΓ为边界,且仅以Γ\GammaΓ为边界);
- 方向匹配:Γ\GammaΓ的正向与Σ\SigmaΣ的侧必须满足右手规则(若方向相反,公式需加负号);
- 函数光滑性:P,Q,RP, Q, RP,Q,R在包含Σ\SigmaΣ的空间区域内一阶偏导数连续。
计算方法总结
写在最前面:计算复杂的曲面积分时,一定要用好曲面表达式这一条件,通过将其代入,简化曲面积分的被积函数
带入曲面表达式化简之后,依次开始考虑,能不能用高斯公式算、能不能用斯托克斯公式算,如果都不能,再考虑能不能转化成一型曲面积分来算,如果都不好办,再考虑投影法,老老实实算
1366

被折叠的 条评论
为什么被折叠?



