考研数学——曲线曲面积分篇

曲线曲面积分

第一类曲线积分

物理背景:曲线的质量计算

假设平面上有一条曲线段 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ρ 为变量,需用“分割、近似、求和、取极限”的积分思想:
    1. 分割:将曲线 LLL 任意分成 nnn 个小弧段 Δs1,Δs2,…,Δsn\Delta s_1, \Delta s_2, \dots, \Delta s_nΔs1,Δs2,,ΔsnΔsi\Delta s_iΔsi 既表示第 iii 段的弧长,也表示该弧段);
    2. 近似:在每个小弧段 Δ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
    3. 求和:总质量的近似值为所有小弧段质量之和:m≈∑i=1nρ(ξi,ηi)Δsim \approx \sum_{i=1}^n \rho(\xi_i, \eta_i) \Delta s_imi=1nρ(ξi,ηi)Δsi
    4. 取极限:记所有小弧段的最大长度为 λ=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λ0i=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_iLf(x,y)ds=limλ0i=1nf(ξi,ηi)Δsi
存在,则称该极限为 函数 f(x,y)f(x,y)f(x,y) 在曲线 LLL 上对弧长的曲线积分(即第一类曲线积分)。

  • LLL 是闭合曲线(首尾相连),积分记为 ∮Lf(x,y) ds\oint_L f(x,y) \, dsLf(x,y)ds
  • f(x,y)=1f(x,y) = 1f(x,y)=1,则积分结果就是曲线 LLL 的长度:∫L1⋅ds=l\int_L 1 \cdot ds = lL1ds=llllLLL 的弧长)。

(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λ0i=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)} dtLf(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)
LLLy=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(α<β)

第一类曲线积分的应用场景

第一类曲线积分的核心应用是计算曲线的“分布量”,常见场景包括:

  1. 曲线的质量:若线密度为 ρ(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

  2. 曲线的重心(形心)

    • 平面曲线的重心 (xˉ,yˉ)(\bar{x}, \bar{y})(xˉ,yˉ)xˉ=1m∫Lxρds\bar{x} = \frac{1}{m} \int_L x\rho dsxˉ=m1Lxρdsyˉ=1m∫Lyρds\bar{y} = \frac{1}{m} \int_L y\rho dsyˉ=m1Lyρds
    • ρ=1\rho = 1ρ=1(均匀曲线),则为形心,此时 m=lm = lm=l(曲线长度)。
  3. 曲线的转动惯量

    • 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)dyLP(x,y)dx+Q(x,y)dy
    其中,dxdxdxdydydy分别是曲线LLLxxx轴和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)dzLP(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沿LLLBBB,则:
    ∫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] dtLPdx+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] dtLPdx+Qdy+Rdz=αβ[Px+Qy+Rz]dt

计算方法2:格林公式

什么是格林公式?

设区域 D⊂R2D\subset \mathbb R^2DR2 满足:

  • 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(xQyP)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(QxPy)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)=Cydx.


② 面积的对称形式

A(D)=12∮C(x dy−y dx) A(D)=\frac12\oint_C (x\,dy - y\,dx) A(D)=21C(xdyydx)


格林公式的物理意义(旋度与环流)

格林公式是二维的 Stokes 定理:

边界上的环流=区域内部旋度之和 \text{边界上的环流}=\text{区域内部旋度之和} 边界上的环流=区域内部旋度之和

  • Qx−Py=0Q_x - P_y = 0QxPy=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(QxPy)dA

要求:

  • CCC 逆时针
  • P,QP,QP,Q 一阶连续可偏导
  • CCC 为简单闭曲线
  • DDDCCC 围成的区域

对于一个积分与路径无关的场,是不是任意闭合曲线在该场中的第二类曲线积分都是0?

当然是的,在无旋场中,积分与路径无关,只与曲线的起点和终点位置有关,因为闭合曲线的起点和终点是重合的,所以无旋场中环路曲线积分始终为0。

既然如此,请你解释解释,在下面这个例子中,通过计算我们发现∂Q∂x=∂P∂y\frac{\partial Q}{\partial x} = \frac{\partial P}{\partial y}xQ=yP,为什么这个闭合环路积分不为0呢?
在这里插入图片描述
原因很简单,因为我们前面说的条件中明确说了 “对于一个积分与路径无关的场”,这就意味着场中的任意一点,都要满足∂Q∂x=∂P∂y\frac{\partial Q}{\partial x} = \frac{\partial P}{\partial y}xQ=yP,而在上面的例子中,(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}xQyP都不存在!于是我们得到的无旋场,是不能包含(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+QdyC2Pdx+Qdy=D(xQyP)dA.

由于在红色区域内,∂Q∂x=∂P∂y\frac{\partial Q}{\partial x} = \frac{\partial P}{\partial y}xQ=yP始终成立,因此

∮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+QdyC2Pdx+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_ii=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λ0i=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=xzzy=∂z∂yz_y = \frac{\partial z}{\partial y}zy=yz1+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面),也可类似地转化为相应的二重积分。

计算过程中的注意事项

  1. 第一类曲面积分是有对称性的! 我们一定要尽可能地利用积分区域的对称性,对曲面积分狠狠进行化简,如果最终里面的被积函数为常数,那么曲面积分就转化为积分曲面的表面积了
  2. 计算复杂的曲面积分时,一定要用好曲面表达式这一条件,通过将其代入,简化曲面积分的被积函数

第二类曲面积分

第二类曲面积分的定义

Σ\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ΔSixOyxOyxOy面上的投影为(Δ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λ0i=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面,区分前侧/后侧),也可类似转化,注意方向对符号的影响。

特殊计算方法:将第二类曲面积分 转化成 第一类曲面积分

  1. 求隐函数的偏导数:计算Fx=∂F∂xF_x = \frac{\partial F}{\partial x}Fx=xFFy=∂F∂yF_y = \frac{\partial F}{\partial y}Fy=yFFz=∂F∂zF_z = \frac{\partial F}{\partial z}Fz=zF
  2. 构造法向量:n⃗=(Fx,Fy,Fz)\vec{n} = (F_x, F_y, F_z)n=(Fx,Fy,Fz),定向决定符号(如取外侧时,法向量指向曲面外);
  3. 计算法向量模长:∣n⃗∣=Fx2+Fy2+Fz2|\vec{n}| = \sqrt{F_x^2 + F_y^2 + F_z^2}n=Fx2+Fy2+Fz2
  4. 求方向余弦:cos⁡α=Fx∣n⃗∣\cos\alpha = \frac{F_x}{|\vec{n}|}cosα=nFxcos⁡β=Fy∣n⃗∣\cos\beta = \frac{F_y}{|\vec{n}|}cosβ=nFycos⁡γ=Fz∣n⃗∣\cos\gamma = \frac{F_z}{|\vec{n}|}cosγ=nFz
    其中cos⁡α\cos\alphacosα是法向量与xxx轴夹角的余弦
  5. 利用两类曲面积分的联系转化
    第二类曲面积分可通过法向量方向余弦转化为第一类曲面积分:
    Σ\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αdSdzdx=cos⁡βdS\mathrm{d}z\mathrm{d}x = \cos\beta \mathrm{d}Sdzdx=cosβdSdxdy=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
    转化为第一类曲面积分后,可利用第一类的计算方法(如投影到坐标面)求解。
  6. 也可以用换底公式
    • 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=Ω(xP+yQ+zR)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}xP+yQ+zR称为向量场v⃗\vec{v}v散度,记为divv⃗\text{div} \vec{v}divv。高斯公式表明:闭曲面的通量等于区域内散度的三重积分,反映了“通量的来源与散度的分布”之间的关系。

高斯公式使用条件

  1. 曲面Σ\SigmaΣ闭曲面(包围一个空间区域Ω\OmegaΩ);
  2. Σ\SigmaΣ外侧(若取内侧,公式需加负号:∯Σ内侧⋯=−∭Ωdivv⃗ dxdydz\oiint_{\Sigma^{\text{内侧}}} \cdots = -\iiint_\Omega \text{div} \vec{v} \, dxdydzΣ内侧=Ωdivvdxdydz);
  3. 函数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=Σ(yRzQ)dydz+(zPxR)dxdz+(xQyP)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}Γvdl=Σ(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=ΣdydzxPdxdzyQdxdyzR
其中,行列式按第一行展开即得展开形式,旋度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=(yRzQ,zPxR,xQyP)

斯托克斯公式的物理意义(环流与旋度)

  • 左侧的第二类曲线积分∮Γv⃗⋅dl⃗\oint_\Gamma \vec{v} \cdot d\vec{l}Γvdl称为向量场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描述了“向量场在某点的局部旋转强度与方向”。
    斯托克斯公式表明:向量场沿闭曲线的环流,等于其旋度通过以该曲线为边界的任意曲面的通量,反映了“整体环流”与“局部旋度”的内在联系。

斯托克斯公式的使用条件

  1. 曲线与曲面的关联Γ\GammaΓΣ\SigmaΣ边界曲线Σ\SigmaΣΓ\GammaΓ为边界,且仅以Γ\GammaΓ为边界);
  2. 方向匹配Γ\GammaΓ的正向与Σ\SigmaΣ的侧必须满足右手规则(若方向相反,公式需加负号);
  3. 函数光滑性P,Q,RP, Q, RP,Q,R在包含Σ\SigmaΣ的空间区域内一阶偏导数连续

计算方法总结

写在最前面:计算复杂的曲面积分时,一定要用好曲面表达式这一条件,通过将其代入,简化曲面积分的被积函数

带入曲面表达式化简之后,依次开始考虑,能不能用高斯公式算、能不能用斯托克斯公式算,如果都不能,再考虑能不能转化成一型曲面积分来算,如果都不好办,再考虑投影法,老老实实算

【复现】基于改进秃鹰算法的微电网群经济优化调度研究(Matlab代码实现)内容概要:本文围绕“基于改进秃鹰算法的微电网群经济优化调度研究”展开,重点介绍了利用改进秃鹰算法(Bald Eagle Search Algorithm, BES)对微电网群进行经济优化调度的Matlab代码实现方法。该研究旨在解决微电网系统中多能源协调、运行成本最小化与供电可靠性之间的平衡问题,通过引入算法改进策略提升寻优能力与收敛速度,有效应对微电网中可再生能源出力波动性和负荷不确定性带来的调度挑战。文中详细阐述了微电网群的系统架构、目标函数构建、约束条件设定以及改进算法的应用流程,并通过仿真实验验证了其在降低运行成本、提高能源利用效率方面的优越性。; 适合人群:具备一定电力系统基础知识和Matlab编程能力的高校研究生、科研人员及从事微电网、智能电网优化调度相关工作的工程技术人员,尤其适合正在开展能源优化、智能算法应用研究的研究者。; 使用场景及目标:①用于微电网群多目标经济调度模型的构建与求解;②作为智能优化算法(如秃鹰算法)改进与应用的学习案例;③支撑科研论文复现、算法对比实验及实际项目中的能源管理系统开发。; 阅读建议:建议结合提供的Matlab代码进行同步调试与仿真,重点关注目标函数设计、约束处理方式与算法改进机制,同时可将其与其他智能算法(如粒子群、遗传算法)进行对比分析,深化对优化调度问题求解路径的理解。
基于分布式模型预测控制的多智能体点对点轨迹生成研究(Matlab代码实现)内容概要:本文围绕“基于分布式模型预测控制的多智能体点对点轨迹生成研究”展开,重点介绍利用分布式模型预测控制(DMPC)方法实现多智能体系统(如无人机、无人车等)在点对点转移过程中的协同轨迹规划。文中结合Matlab代码实现,展示了该方法在固定翼无人机一致性控制、异构车辆编队、多无人机协同路径规划等场景中的应用,强调系统间的通信拓扑、状态约束、避障机制与优化求解过程。同时,文档还列举了大量相关科研方向与代码资源,涵盖智能优化算法、路径规划、电力系统、信号处理等多个领域,旨在为科研人员提供全面的技术支持与复现基础。; 适合人群:具备一定自动化、控制理论或计算机编程基础,从事多智能体系统、路径规划、模型预测控制等相关方向研究的研究生、科研人员及工程技术人员。; 使用场景及目标:①研究多智能体系统在复杂环境下的协同轨迹生成与避障控制;②深入理解分布式模型预测控制(DMPC)的数学建模与优化求解过程;③借助提供的Matlab代码实现算法复现、仿真验证与二次开发,推动科研项目进展。; 阅读建议:建议读者结合文中提供的网盘资源下载完整代码与案例,按照目录顺序系统学习,并重点关注DMPC的构建流程、约束处理与多智能体通信机制,同时可拓展至其他智能优化与控制算法的应用实践中。
评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值