高等数学 第十一章 曲线积分与曲面积分（上）

第十一章 曲线积分与曲面积分

——made by njtech Melody

曲线积分 $\{\begin{array}{l}\text{ 对弧长的曲线积分 }\\ \text{ 对坐标的曲线积分 }\end{array}$ 曲面积分 $\{\begin{array}{l}\text{ 对面积的曲面积分 }\\ \text{ 对坐标的曲面积分 }\end{array}$

*.引言

• 与一元积分类似，都有“大化小”，“常代变”，“近似和”，”取极限“

“大化小”：分为n个小区域

“常代变”：每个小区域中算出题目所需的变量

“近似和”：整体取近似

”取极限“：n取无穷，区域越来越小，结果越来越逼近真实值

理解曲线积分与曲面积分

目前所学的所有积分类型：

• 一重积分：对线段 $x$ 进行微分

$$\underset{x\to 0}{\lim }\sum _{i=1}^{n}f\left({\xi }_t\right)\Delta x_i={\int }_a^{b}f(x)dx$$

• 二重积分：对面积 $xy$ 进行微分

$$\underset{\begin{array}{c}m\to \infty \\ n\to \infty \end{array}}{\lim }\frac{1}{mn}\sum _{i=1}^m\sum _{j=1}^{n}f\left(\frac{i}{m},\frac{j}{n}\right)={\iint }_{D}f(x,y)\mathrm{d}x\mathrm{d}y$$

• 三重积分：对体积 $xyz$ 进行微分

$$\underset{\lambda \to 0}{\lim }\sum _{i=1}^{n}f\left({\xi }_i,{\eta }_i,{\zeta }_i\right)\Delta v_i={\iiint }_{\Omega }f(x,y,z)\mathrm{d}v$$

​ 从积分的几何意义上说，每重积分都是用现有的维度向高纬度宣战（好中二），试图用当前维度的数值来求出高纬度的数值。

• 第一类曲线积分（对弧长的曲线积分）

对三维空间的弧长 $s$ 进行积分，$f(x,y)$ 代表“高度情况”，用二维数值求出三维弧长与底面之间的面积
$$\underset{\lambda \to 0}{\lim }\sum _{i=1}^{n}f\left({\xi }_i,{\eta }_i\right)\Delta S_i={\int }_{L}f(x,y)ds$$

• 第二类曲线积分（对坐标的曲线积分）

由于他是对坐标积分，在几何意义上我暂时理解不了有什么几何意义= =，下面会解释物理意义
$$\underset{l\to 0}{\lim }\sum _{i=1}^{n}p\left({\xi }_i,{\eta }_i\right)\Delta x_i={\int }_{L}Q(x,y)dx$$

• 第一类曲面积分（对面积的曲面积分）

对于被积函数如 $f(x,y,z)$ 意味着我们不能直接理解他在四维空间到底指的是什么，降到三维空间而言，他代表这一个曲面的“密度”，进而求出来的是曲面的质量
$$\underset{\lambda \to 0}{\lim }\sum _{i=1}^m\rho \left({\xi }_i,{\eta }_i,{\zeta }_i\right)\Delta S_i={\iint }_z\rho (x,y,z)dS$$

• 第二类曲面积分（对坐标的曲面积分）

同曲线积分，我也不会解释捏，下面会解释物理意义

一、对弧长的曲线积分

1.定义

设 $L$ 为 $xOy$ 平面内有限的光滑或逐段光滑的曲线段, 函数 $f(x,y)$ 在曲线 $L$ 上有界, 将 $L$ 分为 $n$ 个小的曲线段, 记为 $\Delta s_1,\Delta s_2,\cdots ,\Delta s_n$, 任取 $\left({\xi }_i,{\eta }_i\right)\in \Delta s_i$, 作 ${\sum }_{i=1}^{n}f\left({\xi }_i,{\eta }_i\right)\Delta s_i$, 令 $\lambda =$ max ⁡ 1 ⩽ i ⩽ n { Δ s i } \max _{1 \leqslant i \leqslant n}\left\{\Delta s_{i}\right\} max1⩽i⩽n​{Δsi​}, 若 ${\lim }_{\lambda \to 0}{\sum }_{i=1}^{n}f\left({\xi }_i,{\eta }_i\right)\Delta s_i$ 存在, 称此极限为函数 $f(x,y)$ 在曲线段 $L$ 上对弧长的曲线 积分, 记为 ${\int }_{L}f(x,y)\mathrm{d}s$, 即 ${\int }_{L}f(x,y)\mathrm{d}s={\lim }_{\lambda \to 0}{\sum }_{i=1}^{n}f\left({\xi }_i,{\eta }_i\right)\Delta s_i$.

2.物理意义

设 $L$ 为 $xOy$ 平面内有限的曲线段, 其线密度为 $\rho (x,y)$, 求其质量 $m$ 的过程如下:

(1) 将 $L$ 分为 $n$ 个小的曲线段, 记为 $\Delta s_1,\Delta s_2,\cdots ,\Delta s_n$;

(2) 任取 $\left({\xi }_i,{\eta }_i\right)\in \Delta s_i$, 则 $\Delta m_i\approx \rho \left({\xi }_i,{\eta }_i\right)\Delta s_i(i=1,2,\cdots ,n)$, 于是 $m\approx {\sum }_{i=1}^n\rho \left({\xi }_i,{\eta }_i\right)\Delta s_i$;

(3) 令 λ = max ⁡ 1 ⩽ i ⩽ n { Δ s i } \lambda=\max _{1 \leqslant i \leqslant n}\left\{\Delta s_{i}\right\} λ=max1⩽i⩽n​{Δsi​}, 则 $m={\lim }_{\lambda \to 0}{\sum }_{i=1}^n\rho \left({\xi }_i,{\eta }_i\right)\Delta s_i$.

3.基本性质

• 积分区域可加性

$${\int }_L[af(x,y)+bg(x,y)]\mathrm{d}s=a{\int }_{L}f(x,y)\mathrm{d}s+b{\int }_{L}g(x,y)\mathrm{d}s$$

• 积分限可加性

$${\int }_{L}f(x,y)\mathrm{d}s={\int }_{L_1}f(x,y)\mathrm{d}s+{\int }_{L_2}f(x,y)\mathrm{d}s,其中L=L_1+L_2$$

• 积分区域对称性（奇偶性）

  （1）若曲线 $L$ 关于 $y$ 轴 (即关于变量 $x$ ) 对称, $L_1$ 为位于 $y$ 轴右侧的部分,

  若 $f(-x,y)=-f(x,y)$, 则 ${\int }_{L}f(x,y)\mathrm{d}s=0$;

  若 $f(-x,y)=f(x,y)$, 则 ${\int }_{L}f(x,y)\mathrm{d}s=2{\int }_{L_1}f(x,y)\mathrm{d}s$

  （2）若曲线 $L$ 关于 $x$ 轴 (即关于变量 $y$ ) 对称, $L_1$ 为位于 $x$ 轴上侧的部分,

  若 $f(x,-y)=-f(x,y)$, 则 ${\int }_{L}f(x,y)\mathrm{d}s=0$;

  若 $f(x,-y)=f(x,y)$, 则 ${\int }_{L}f(x,y)\mathrm{d}s=2{\int }_{L_1}f(x,y)\mathrm{d}s$.

  （3）若 $L$ 关于直线 $y=x$ 对称, 则 ${\int }_{L}f(x,y)\mathrm{d}s={\int }_{L}f(y,x)\mathrm{d}s$.

• ${\int }_L\mathrm{d}s=l$, 其中 $l$ 为曲线 $L$ 的长度.

4.计算方法

①特殊代替法+奇偶性判断法

• 由于是对弧长进行微分，往不严谨的角度可以说微分的东西跟 $x,y$ 无关，因此可以直接把题目给的线段式子代进积分函数，再结合奇偶性计算。

【例】计算 ${\int }_L\left(x^2-2xy+3y^2\right)\mathrm{d}s$, 其中 $L:x^2+y^2=1$.

【解】由对称性与奇偶性得 ${\int }_L-2xy\mathrm{d}s=0$ ，所以 $I={\int }_L\left(x^2-2xy+3y^2\right)\mathrm{d}s={\int }_L\left(x^2+3y^2\right)\mathrm{d}s$.

​ 因为 $L$ 关于直线 $y=x$ 对称, 所以 $I={\int }_L\left(y^2+3x^2\right)\mathrm{d}s$.

​ 于是 $2I={\int }_L\left(x^2+3y^2\right)\mathrm{d}s+{\int }_L\left(y^2+3x^2\right)\mathrm{d}s=4{\int }_L\left(x^2+y^2\right)\mathrm{d}s=4{\int }_L\mathrm{d}s=8\pi$

​ 故 $I={\int }_L\left(x^2-2xy+3y^2\right)\mathrm{d}s=4\pi$.

②参数法（定积分法）

情形一、设 $L:y=\phi (x)(a⩽x⩽b)$, 则 ${\int }_{L}f(x,y)\mathrm{d}s={\int }_a^{b}f[x,\phi (x)]\sqrt{1+{\phi }^{\prime 2}(x)}\mathrm{d}x$.

情形二、设 $L:\{\begin{array}{l}x=\phi (t),\\ y=\psi (t),\end{array}(\alpha ⩽t⩽\beta )$, 则
$${\int }_{L}f(x,y)\mathrm{d}s={\int }_a^{\beta }f[\phi (t),\psi (t)]\sqrt{{\phi }^{\prime 2}(t)+{\psi }^{\prime 2}(t)}\mathrm{d}t.$$

• 极坐标法其实就是情形二、只不过变量由 $t$ 变成了 $\theta$ 计算方法同理（适用式子带 $x^2+y^2$）

  设 $L:\{\begin{array}{l}x=\rho cos\theta ,\\ y=\rho sin\theta ,\end{array}(\alpha ⩽\theta ⩽\beta )$, 则
  $${\int }_{L}f(x,y)\mathrm{d}s={\int }_a^{\beta }f[\rho cos\theta ,\rho sin\theta ]\rho \mathrm{d}\rho .$$

【例】 ${\int }_L(x+y)ds$ ，其中 $L$ 为连接（1,0）及（0,1）两点的直线段

【解】线段式子：$y=1-x$ ，$\sqrt{1+y^{\prime 2}(x)}\mathrm{d}x=\sqrt{2}$ ，积分上下限 $(0⩽x⩽1)$
$$原式={\int }_0^{1}1\sqrt{2}dx=\sqrt{2}$$

【例】计算 ${\int }_{L}x{\mathrm{e}}^{\sqrt{x^2+y^2}}\mathrm{d}s$, 其中 $L$ : $x^2+y^2=4,(0⩽t⩽\frac{\pi }{4})$

【解】令 $L_{}:\{\begin{array}{l}x=2\cos t,\\ y=2\sin t\end{array}\left(0⩽t⩽\frac{\pi }{4}\right)$,
$${\int }_{L_{}}x{\mathrm{e}}^{\sqrt{x^2+y^2}}\mathrm{d}s={\mathrm{e}}^2{\int }_{L_{}}x\mathrm{d}s={\mathrm{e}}^2{\int }_0^{\frac{\pi }{4}}2\cos t\cdot 2\mathrm{d}t=2\sqrt{2}{\mathrm{e}}^2,$$

二、对坐标的曲线积分

1.定义

• 二维空间

$L$ 为 $xOy$ 平面上的有向曲线段, 函数 $P(x,y),Q(x,y)$ 在 $L$ 上有界,

若 ${\lim }_{\lambda \to 0}{\sum }_{i=1}^{n}P\left({\xi }_i,{\eta }_i\right)\Delta x_i$ 存在, 称此极限为函数 $P(x,y)$ 沿有向曲线段 $L$ 对坐标 $x$ 的曲线积分, 记 为 ${\int }_{L}P(x,y)\mathrm{d}x$ ，即 ${\int }_{L}P(x,y)\mathrm{d}x={\lim }_{\lambda \to 0}{\sum }_{i=1}^{n}P\left({\xi }_i,{\eta }_i\right)\Delta x_i$ ；

若 ${\lim }_{\lambda \to 0}{\sum }_{i=1}^{n}Q\left({\xi }_i,{\eta }_i\right)\Delta y_i$ 存在, 称此极限为函数 $Q(x,y)$ 沿有向曲线段 $L$ 对坐标 $y$ 的曲线积 分, 记为 ${\int }_{L}Q(x,y)\mathrm{d}y$, 即 ${\int }_{L}Q(x,y)\mathrm{d}y={\lim }_{\lambda \to 0}{\sum }_{i=1}^{n}Q\left({\xi }_i,{\eta }_i\right)\Delta y_i$

• 三维空间

$L$ 为三维空间的有向曲线段,函数 $P(x,y,z),Q(x,y,z),R(x,y,z)$ 在$L$ 上有界,

若 ${\lim }_{\lambda \to 0}{\sum }_{i=1}^{n}P\left({\xi }_i,{\eta }_i,{\zeta }_i\right)\Delta x_i$ 存在, 称此极限为函数 $P(x,y,z)$ 沿有向曲线段 $L$ 对坐 标 $x$ 的曲线积分, 记为 ${\int }_{L}P(x,y,z)\mathrm{d}x$, 即 ${\int }_{L}P(x,y,z)\mathrm{d}x={\lim }_{\lambda \to 0}{\sum }_{i=1}^{n}P\left({\xi }_i,{\eta }_i,{\zeta }_i\right)\Delta x_i$;

若 ${\lim }_{\lambda \to 0}{\sum }_{i=1}^{n}Q\left({\xi }_i,{\eta }_i,{\zeta }_i\right)\Delta y_i$ 存在, 称此极限为函数 $Q(x,y,z)$ 沿有向曲线段 $L$ 对坐标 $y$ 的 曲线积分, 记为 ${\int }_{L}Q(x,y,z)\mathrm{d}y$, 即 ${\int }_{L}Q(x,y,z)\mathrm{d}y={\lim }_{\lambda \to 0}{\sum }_{i=1}^{n}Q\left({\xi }_i,{\eta }_i,{\zeta }_i\right)\Delta y_i$ ；

若 ${\lim }_{\lambda \to 0}{\sum }_{i=1}^{n}R\left({\xi }_i,{\eta }_i,{\zeta }_i\right)\Delta z_i$ 存在, 称此极限为函数 $R(x,y,z)$ 沿有向曲线段 $L$ 对坐标 $z$ 的 曲线积分, 记为 ${\int }_{L}R(x,y,z)\mathrm{d}z$, 即 ${\int }_{L}R(x,y,z)\mathrm{d}z={\lim }_{\lambda \to 0}{\sum }_{i=1}^{n}R\left({\xi }_i,{\eta }_i,{\zeta }_i\right)\Delta z_i$.

2.物理意义（以二维空间为例）

• ：物体在变力作用下沿有向曲线段做功

设 F ⃗ ( x , y ) = { P ( x , y ) , Q ( x , y ) } , L \vec{F}(x, y)=\{P(x, y), Q(x, y)\}, L F (x,y)={P(x,y),Q(x,y)},L 为 $xOy$ 平面内的有向曲线段, 物体在 $\vec{F}(x,y)$ 的作用下沿 $L$ 从起点 $A$ 到达终点 $B$, 元素法求力所做的功的方法如下:
取 $\mathrm{d}\vec{s}\subset L$, 其中 d s = { d x ,   d y } \mathrm{d} s=\{\mathrm{d} x, \mathrm{~d} y\} ds={dx, dy}, 则 $\mathrm{d}W=\vec{F}(x,y)\cdot \mathrm{d}\vec{s}=P(x,y)\mathrm{d}x+$ $Q(x,y)\mathrm{d}y$,于是$W={\int }_{L}P(x,y)\mathrm{d}x+Q(x,y)\mathrm{d}y$.

3.基本性质（二维空间为例）

• ${\int }_{L^{-}}P(x,y)\mathrm{d}x+Q(x,y)\mathrm{d}y=-{\int }_{L}P(x,y)\mathrm{d}x+Q(x,y)\mathrm{d}y.$
• ${\int }_{L}P(x,y)\mathrm{d}x+Q(x,y)\mathrm{d}y={\int }_{L_1}P(x,y)\mathrm{d}x+Q(x,y)\mathrm{d}y+{\int }_{L_2}P(x,y)\mathrm{d}x+Q(x,y)\mathrm{d}y$,
  其中 $L=L_1+L_2$.

4.计算方法

①参数法（定积分）

• 参数法适用于二维、三维空间

情形一、设 $L:y=\phi (x)(a⩽x⩽b)$, 则
∫ L P ( x , y ) d x + Q ( x , y ) d y = ∫ a b { P [ x , φ ( x ) ] + Q [ x , φ ( x ) ] φ ′ ( x ) } d x . \int_{L} P(x, y) \mathrm{d} x+Q(x, y) \mathrm{d} y=\int_{a}^{b}\left\{P[x, \varphi(x)]+Q[x, \varphi(x)] \varphi^{\prime}(x)\right\} \mathrm{d} x . ∫L​P(x,y)dx+Q(x,y)dy=∫ab​{P[x,φ(x)]+Q[x,φ(x)]φ′(x)}dx.

情形二、设 $L:\{\begin{array}{l}x=\phi (t),\\ y=\psi (t),\end{array}(\alpha ⩽t⩽\beta )$, 则
∫ L P ( x , y ) d x + Q ( x , y ) d y = ∫ a β { P [ φ ( t ) , ψ ( t ) ] φ ′ ( t ) + Q [ φ ( t ) , ψ ( t ) ] ψ ′ ( t ) } d t . \int_{L} P(x, y) \mathrm{d} x+Q(x, y) \mathrm{d} y=\int_{a}^{\beta}\left\{P[\varphi(t), \psi(t)] \varphi^{\prime}(t)+Q[\varphi(t), \psi(t)] \psi^{\prime}(t)\right\} \mathrm{d} t . ∫L​P(x,y)dx+Q(x,y)dy=∫aβ​{P[φ(t),ψ(t)]φ′(t)+Q[φ(t),ψ(t)]ψ′(t)}dt.

• 极坐标法其实就是情形二、只不过变量由 $t$ 变成了 $\theta$ 计算方法同理（适用式子带 $x^2+y^2$）

  设 $L:\{\begin{array}{l}x=\rho cos\theta ,\\ y=\rho sin\theta ,\end{array}(\alpha ⩽\theta ⩽\beta )$, 则
  ∫ L P ( x , y ) d x + Q ( x , y ) d y = ∫ a β { P [ ρ c o s θ , ρ s i n θ ] ( − ρ s i n θ ) + Q [ ρ c o s θ , ρ s i n θ ] ρ c o s θ } d t . \int_{L} P(x, y) \mathrm{d} x+Q(x, y) \mathrm{d} y=\int_{a}^{\beta}\left\{P[\rho cos \theta, \rho sin \theta] (-\rho sin \theta)+Q[\rho cos \theta, \rho sin \theta] \rho cos \theta\right\} \mathrm{d} t . ∫L​P(x,y)dx+Q(x,y)dy=∫aβ​{P[ρcosθ,ρsinθ](−ρsinθ)+Q[ρcosθ,ρsinθ]ρcosθ}dt.

【例】${\int }_{\Gamma }x\mathrm{d}x+y\mathrm{d}y+(x+y-1)\mathrm{d}z$, 其中 $\Gamma$ 是从点 $(1,1,1)$ 到点 $(2,3,4)$ 的一段直线.

【解】
$$\begin{array}{l}\{\begin{array}{l}x=1+t\\ y=1+2t\\ z=1+3t\end{array}\Rightarrow \{\begin{array}{l}dx=dt\\ dy=2dt\\ dz=3dt\end{array}\\ {\int }_0^1[(1+t)\times 1+(1+2t)\times z+(1+3t)\times 3]dt=13\end{array}$$

【例】 ${\oint }_{L}xy\mathrm{d}x$, 其中 $L$ 为圆周 $(x-R{)}^2+y^2=R^2(R>0)$ 及 $x$ 轴所围成的在第一象限内的区域的整个边界(按逆时针方向绕行).

【解】 ${\oint }_{L}xy\mathrm{d}x=L1(圆弧部分)+L2(坐标轴部分)$

​ L2：y=0，∴${\oint }_{L_2}xy\mathrm{d}x=0$

​ L1：令 { x = R + R c o s θ , y = R s i n θ } \{ x=R+Rcos\theta ,y=Rsin\theta \} {x=R+Rcosθ,y=Rsinθ} $(0\le \theta \le \pi )$
$${\oint }_{L_1}xy\mathrm{d}x={\int }_0^{\pi }{R}^2(1+cos\theta )Rsin\theta \cdot (-Rsin\theta )d\theta =-\frac{\pi }{2}{R}^3$$

②格林公式（二维空间）

• 格林公式是牛顿-莱布尼茨公式的推广。

牛顿-莱布尼茨公式：如果函数 $f(x)$在区间 $[a,b]$上连续，并且存在原函数 $F(x)$

$${\int }_a^{b}f(x)dx=F(b)-F(a)$$

我的理解： $F(x)$ 对 $x$ 微分（求导），实现将“边界计算”转为“长度计算”

格林公式：设 $D$ 为 $xOy$ 平面上连通的有限闭区域, $L$ 为闭区域 $D$ 的正向边界, 函数 $P(x,y)$, $Q(x,y)$ 在 $D$ 上连续可偏导, 则

$${\oint }_{L}P(x,y)dx+Q(x,y)dy={\iint }_D\left[\frac{\partial Q}{\partial x}-\frac{\partial P}{\partial y}\right]dxdy$$

• 正向边界：指闭区域 $D$ 始终在边界的“左边”，即外逆里顺

我的理解：同上，对 $P(x,y)$ 求 $x,y$ 微分（求全导），实现将”边界计算“转化为“面积计算”，只不过 $P(x)$已经对x求偏导了，那就让他对 $y$ 求一次偏导，同理对 $Q(x,y)$ 也是如此，让他对 $x$ 求一次偏导

格林公式的计算

单连通区域（无洞）

【例】${\oint }_L{x}^{2}ydx-x{y}^{2}dy，其中L:x^2+y^2=a^2\text{ 逆时针 }$

【解】令 $P=x^{2}y,Q=-x{y}^2$ $\frac{\partial p}{\partial y}=x^2\frac{\partial Q}{\partial x}=-y^2$ $\frac{\partial Q}{\partial x}-\frac{\partial P}{\partial y}=-(x^2+y^2)$

​ $原式={\iint }_D\left(-y^2-x^2\right)dxdy=-{\int }_0^{2\pi }d\theta {\int }_0^a{\rho }^2\cdot \rho d\rho =-\frac{\pi a^4}{2}$

复连通区域（有洞）

【例】$I={\oint }_L\frac{ydx-xdy}{2\left(x^2+y^2\right)}$ ，其中 $L:(x-1{)}^2+y^2=2$ （逆时针）

【解】显然在（0,0）处不可取，故是一个复连通区域，要挖一个洞。$I=D$，$D=D_1+D_2$

对于$D1$：（第一次格林公式：边界→面积）

令 $P=\frac{y}{2(x^2+y^2)}$，$Q=\frac{-x}{2(x^2+y^2)}$

$p_y^{\prime }=\frac{2\left(x^2+y^2\right)-y\cdot 4y}{4\left(x^2+y^2\right)}=\frac{2x^2-2y^2}{4{\left(x^2+y^2\right)}^2}$ $Q^{\prime }x=\frac{-2\left(x^2+y^2\right)+x\cdot 4x}{4{\left(x^2+y^2\right)}^2}=\frac{2x^2-2y^2}{4{\left(x^2+y^2\right)}^2}$ ∴ $Q^{\prime }x-P^{\prime }y=0$
$$l_1+l_2={\oint }_{l_1}\frac{ydx-xdy}{2\left(x^2+y^2\right)}+{\oint }_{l_2}\frac{ydx-xdy}{2\left(x^2+y^2\right)}={\iint }_{D_1}0dxdy=D_1=0$$
对于D2：（第二次格林公式：面积→边界）
$$D_2=-l_2=-{\oint }_{l_2}\frac{ydx-xdy}{2\left(x^2+y^2\right)}=\frac{1}{2r^2}{\int }_{l_2}-ydx+xdy$$
令 $P=-y$，$Q=x$ $p_y^{\prime }=-1$ $Q_x^{\prime }=1$ ∴ $Q^{\prime }x-P^{\prime }y=-2$
$$\frac{1}{2S^2}{\int }_c-ydx+xdy=-\frac{1}{r^2}{\iint }_{D2}dxdy=-\frac{1}{r^2}\cdot \pi r^2=-\pi$$
综上：$I={\oint }_L\frac{ydx-xdy}{2\left(x^2+y^2\right)}=-\pi$

格林公式的应用

求面积：
$$S_D={\iint }_{D}dxdy=\frac{1}{2}{\oint }_{L}xdy-ydx$$
【例】计算星形线 $x=a{\cos }^{3}t,y=a{\sin }^{3}t$ 所围图形的面积.

【解】
$$\begin{array}{rl}A& =\frac{1}{2}{\oint }_{L}x\mathrm{d}y-y\mathrm{d}x=\frac{1}{2}{\int }_0^{2\pi }\left[a{\cos }^{3}t{\left(a{\sin }^{3}t\right)}^{\prime }-a{\sin }^{3}t{\left(a{\cos }^{3}t\right)}^{\prime }\right]\mathrm{d}t\\ & =\frac{1}{2}{\int }_0^{2\pi }3a^2{\sin }^{2}t{\cos }^{2}t\mathrm{d}t=\frac{3a^2}{2}{\int }_0^{2\pi }\left(\frac{1-\cos 4t}{8}\right)\mathrm{d}t=\frac{3}{8}\pi a^2.\end{array}$$

曲线积分与路径无关的条件

解释：对于满足一些条件的曲线，起点和终点的位置固定，沿不同的路线积分，其积分值相同，即曲线积分只与起点和终点有关，与路线的选取无关。

定理 ：设 $D$ 为单连通区域,函数 $P(x,y),Q(x,y)$ 在区域 $D$ 内连续可偏导, 则下列四个命题等价：

• 曲线积分 ${\int }_{L}P(x,y)\mathrm{d}x+Q(x,y)\mathrm{d}y$ 与路径无关.

• 对区域 $D$ 内任意闭曲线 $C$, 有 ${\oint }_{C}P(x,y)\mathrm{d}x+Q(x,y)\mathrm{d}y=0$.

• 区域 $D$ 内恒有 $\frac{\partial Q}{\partial x}=\frac{\partial P}{\partial y}$ (柯西 $-$ 黎曼条件).

• 在区域 $D$ 内存在二元函数 $u(x,y)$, 使得 $\mathrm{d}u(x,y)=P(x,y)\mathrm{d}x+Q(x,y)\mathrm{d}y$.

利用该条件计算曲线积分

【例】已知点 $O(0,0)$ 及点 $A(1,1)$, 且曲线积分 $I={\int }_{\text{C }}\left(ax\cos y-y^2\sin x\right)\mathrm{d}x+\left(by\cos x-x^2\sin y\right)\mathrm{d}y$ 与路径无关,试确定常数 $a,b$, 并求 $I$.
【思路探索】由积分与路径无关, 可得 $\frac{\partial P}{\partial y}=\frac{\partial Q}{\partial x}$, 利用恒等式求 $a,b$; 再选择简单路径 $(0,0)$ 到 $(0,1)$ 和 $(0,1)$ 到 $(1,1)$ 求曲线积分.
解: 令 $P(x,y)=ax\cos y-y^2\sin x,Q(x,y)=by\cos x-x^2\sin y$, $\frac{\partial P}{\partial y}=-ax\sin y-2y\sin x,\frac{\partial Q}{\partial x}=-by\sin x-2x\sin y$
由题意知 $\frac{\partial P}{\partial y}=\frac{\partial Q}{\partial x}$,得 $a=b=2$,所以
$$\begin{array}{rl}I& ={\int }_{(0,0)}^{(1,1)}P\mathrm{d}x+Q\mathrm{d}y={\int }_{(0,0)}^{(0,1)}P\mathrm{d}x+Q\mathrm{d}y+{\int }_{(0,1)}^{(1,1)}P\mathrm{d}x+Q\mathrm{d}y\\ & ={\int }_0^{1}Q(0,y)\mathrm{d}y+{\int }_0^{1}P(x,1)\mathrm{d}x={\int }_0^{1}2y\mathrm{d}y+{\int }_0^1(2x\cos 1-\sin x)\mathrm{d}x=2\cos 1\end{array}$$

求全微分方程

定理：若一阶微分方程为 $P(x,y)dx+Q(x,y)dy $ 若曲线积分与路径无关，即$\frac{\partial Q}{\partial x}=\frac{\partial P}{\partial y}$ ，则称此方程为全微分方程，其隐式通解为： $u(x,y)={\int }_{x_0,y_0}^{(x,y)}P(x,y)dx+Q(x,y)dy=C$

【例】设 $P\mathrm{d}x+Q\mathrm{d}y=\left(x^2+2xy-y^2\right)\mathrm{d}x+\left(x^2-2xy-y^2\right)\mathrm{d}y$, 求 $u(x,y)$, 使得 $\mathrm{d}u(x,y)=P\mathrm{d}x$ $+Q\mathrm{d}y$.
【解】
由 $\frac{\partial u}{\partial x}=x^2+2xy-y^2$, 积分得 $u=\int \left(x^2+2xy-y^2\right)\mathrm{d}x+\phi (y)=\frac{1}{3}{x}^3+x^{2}y-x{y}^2+\phi (y).$
对 $y$ 求导得 $\frac{\partial u}{\partial y}=x^2-2xy+{\phi }^{\prime }(y)=Q=x^2-2xy-y^2$,
所以 ${\phi }^{\prime }(y)=-y^2$, 从而 $\phi (y)=-\frac{1}{3}{y}^3+C$,故 $u(x,y)=\frac{1}{3}{x}^3+x^{2}y-x{y}^2-\frac{1}{3}{y}^3+C$.

③斯托克斯公式