第二型曲线积分 Green公式

第七周定义了第二型曲线积分(物理含义是变力做功)及其计算方法,GreenGreen 公式定义基本完成。

第二型曲线积分

AB\mathop{AB}\limits^{-\rightarrow}nn 维空间中的向量,则称它的单位向量为

AB^=ABAB\widehat{\mathop{AB}\limits^{\rightarrow}}=\frac{\mathop{AB}\limits^{-\rightarrow}}{|\mathop{AB}\limits^{-\rightarrow}|}

F:ARnRn\vec{F}:A\subset \mathbb R^n\rightarrow \mathbb R^n 为向量场(类似于物理中的重力场,磁场)。

CRnC\subset \mathbb R^n 是光滑曲线,τ:CRn, τ=1\vec{\tau}:C\rightarrow \mathbb R^n,\ |\vec{\tau}| = 1τ(P)\vec{\tau}(P)CCPP 处的单位切向量,则称 τ\vec{\tau}CC 的单位切向量场。

定义1(定向)

如果 τ\vec{\tau}CC连续单位切向量场,则称 τ\vec{\tau}CC 的一个定向


因为 CC 上每一点的切向量都有两个方向,连续性要求不能两个相邻点方向是完全相反的,也就是说一个领域中的点都是指向一个方向的。

命题2(定向只有两个)

CRnC\subset \mathbb R^n 为光滑曲线,则 CC 有且仅有两个定向。


思路: 使用 CC 的参数表达式,构造出两个定向,注意定义域是否一样(一个是参数,一个是坐标),使用内积和函数的连续性证明唯二性。

证明:

  1. 存在性,设 α:[a,b]C\alpha:[a, b]\rightarrow CCC 的参数方程,则 CC 有两个定向,分别是:τ=α^α1\vec{\tau} = \widehat{\alpha'}\circ\alpha^{-1}τ-\vec{\tau}

  2. 唯二性,设 λ\vec{\lambda} 也是 CC 的定向,则 PC\forall P\in C,对两个定向 λ,τ\vec{\lambda}, \vec{\tau}PP 点处做内积,有

(λτ)(P)=±1=======则 λτ 也是连续的由于 λ,τ 都是连续的λτ1 或 λτ1\begin{aligned} (\vec{\lambda}\cdot\vec{\tau})(P) = \pm 1\mathop{=======\Rightarrow}\limits_{\text{则 }\vec{\lambda}\cdot\vec{\tau}\text{ 也是连续的}}^{\text{由于 }\vec{\lambda},\vec{\tau}\text{ 都是连续的}}\vec{\lambda}\cdot\vec{\tau}\equiv 1\text{ 或 }\vec{\lambda}\cdot\vec{\tau}\equiv -1 \end{aligned}

命题3(定向曲线)

τ\vec{\tau}CC 的一个定向,称 {C,τ}\{C,\vec{\tau}\} 是一个定向曲线

简记为 CCCC 是一个定向曲线,τ\vec{\tau}CC 的正向,如果是负向则记为 C={C,τ}-C = \{C, -\vec{\tau}\}


定义4(第二型曲线积分)

CRnC\subset \mathbb R^n 为定向曲线,τ\vec{\tau}CC 的正向,设 F:CRn\vec{F}:C\rightarrow \mathbb R^n

α:[a,b]C\alpha : [a, b]\rightarrow CCC 的参数方程,τ=α^α1\vec{\tau} = \widehat{\alpha'}\circ\alpha^{-1},方向为 α\alpha 参数增加的方向。

[a,b][a, b] 的一个分划 π:a=t0<t1<<tN=b\pi : a =t_0 < t_1 < \cdots < t_N = b,考虑和式:

S=i=1NF(ξi)τ(ξi)L(Ci)S = \sum_{i=1}^N\vec{F}(\xi_i)\cdot\vec{\tau}(\xi_i)L(C_i)

其中 Ci=α([ti1,ti]), ξiCiC_i = \alpha([t_{i-1}, t_i]),\ \xi_i\in C_i

如果存在 IRI\in \mathbb R,使得 ε>0\forall \varepsilon > 0,有 ISε|I - S| \leqslant \varepsilon,则称极限 limΔπ0S\lim\limits_{\Delta\pi\rightarrow 0}S 存在,且 limΔπ0S=I\lim\limits_{\Delta\pi\rightarrow 0} S = I

如果 limΔπ0S\lim\limits_{\Delta\pi\rightarrow 0} S 存在,定义

CFds=limΔπ0S\int_C\vec{F}\cdot d\vec{s} = \lim_{\Delta\pi\rightarrow 0} S

称为向量场 F\vec{F} 沿定向曲线 CC 的积分。

定理5(性质及计算方法)

CRnC\subset \mathbb R^n 为定向曲线,τ\vec{\tau}CC 的正向,α:[a,b]C\alpha:[a, b]\rightarrow CCC 的参数方程,且 CC 的正向是 α\alpha 参数增加的方向,设 F,G:CRn\vec{F},\vec{G}:C\rightarrow \mathbb R^n 连续,则

  1. CFds=CFτds\displaystyle \int_C\vec{F}\cdot d\vec{s} = \int_C\vec{F}\cdot\vec{\tau}\,ds(转化为第一型曲线积分)

  2. CFds=abF(α(t))α(t)dt\displaystyle \int_C\vec{F}\cdot d\vec{s} = \int_a^b\vec{F}(\alpha(t))\cdot\alpha'(t)\,dt

  3. k,lRk, l\in \mathbb R,则 C{kF+lG}ds=kCFds+lCGds\displaystyle \int_C\{k\vec{F} + l\vec{G}\}\cdot d\vec{s} = k\int_C\vec{F}\cdot d\vec{s}+l\int_C\vec{G}\cdot d\vec{s}(线性性)

  4. CFds=CFds\displaystyle \int_{-C}\vec{F}\cdot d\vec{s} = -\int_C\vec{F}\cdot d\vec{s}


思路: 通过第一型曲线积分的和形式即可转化为第二型曲线积分的和形式。

通过 τ\vec{\tau} 的参数化表达式,转化为参数形式,转化为一维积分形式以后,利用一维积分的性质,③④③④ 就不难证明了。

证明: (就证明 吧)

CFds=CFτds=CF(ααα1)ds=abF(α(t))αα(t)α(t)dt=abF(α(t))α(t)dt\begin{aligned} \int_C\vec{F}\cdot d\vec{s} = \int_C\vec{F}\cdot \vec{\tau}\,ds &= \int_C\vec{F}\cdot (\frac{\alpha'}{|\alpha'|}\circ\alpha^{-1})\,ds\\ &= \int_a^b\vec{F}(\alpha(t))\cdot\frac{\alpha'}{|\alpha'|}(t)\cdot|\alpha'(t)|\,dt\\ &= \int_a^b\vec{F}(\alpha(t))\cdot\alpha'(t)\,dt \end{aligned}


对于 dsd\vec{s} 的理解:

ds=τdsτ:曲线切线方向,ds:弧长的微分=α(t)dt位移=速度×时间\begin{aligned} d\vec{s} &= \vec{\tau}\,ds &\vec{\tau}:\text{曲线切线方向},ds:\text{弧长的微分}\\ &= \alpha'(t)\,dt &\text{位移=速度}\times\text{时间} \end{aligned}


更加一般的,在具体的坐标系中计算第二型曲线积分:

OXOX 为直角坐标系,x=(x1,x2,,xn)x=(x_1,x_2,\cdots, x_n),记 ds=(dx1,dx2,,dxn)d\vec{s} = (dx_1,dx_2,\cdots,dx_n)

向量场 F=(F1(x),F2(x),,Fn(x))\vec{F} = (F_1(x), F_2(x),\cdots, F_n(x)),则

CFds=CF1(x)dx1++Fn(x)dxn\int_C\vec{F}\cdot d\vec{s} = \int_CF_1(x)\,dx_1+\cdots+F_n(x)\,dx_n

光滑曲线 CC 的参数方程 α(t)=(x1(t),x2(t),,xn(t))=x(t)\alpha(t) = (x_1(t),x_2(t),\cdots,x_n(t)) = x(t),则

CF1(x)dx1++Fn(x)dxn=abF(α(t))α(t)dt=ab(F1(x(t))x1(t)++Fn(x(t))xn(t))dt\begin{aligned} \int_CF_1(x)\,dx_1+\cdots+F_n(x)\,dx_n &= \int_a^b\vec{F}(\alpha(t))\cdot\vec{\alpha}'(t)\,dt \\ &= \int_a^b(F_1(x(t))x_1'(t)+\cdots+F_n(x(t))x_n'(t))\,dt\\ \end{aligned}

从形式上看,就像是做了变量代换 dx1(t)=x1(t)dtdx_1(t) = x_1'(t)\,dt

定义6(分段光滑曲线的第二型曲线积分)

C=i=1NCiC = \bigcup\limits_{i=1}^NC_iCiC_i 为定向曲线,CiC_i 内部互不相交,F:CRnF:C\rightarrow \mathbb R^n 连续,则

CFds=i=1nCiFds\int_C\vec{F}\cdot d\vec{s} = \sum_{i=1}^n\int_{C_i}\vec{F}\cdot d\vec{s}

CC 分段光滑,设 α:[a,b]C\alpha:[a, b]\rightarrow CCC 的参数方程,设 F:CRn\vec{F}:C\rightarrow \mathbb R^n 连续,C=C1+C2++CNC = C_1+C_2+\cdots+C_N,且 CiC_i 的正向都是 α\alpha 的参数增加的方向,则

CFds=abF(α(t))α(t)dt\int_C\vec{F}\cdot d\vec{s} = \int_a^b\vec{F}(\alpha(t))\cdot\vec{\alpha}'(t)\,dt

例一

Γ={(x,y,z)R3:x2+y2+z2=a2,x+y+z=0}\Gamma = \{(x, y, z)\in \mathbb R^3:x^2+y^2+z^2 = a^2, x+y+z = 0\}Γ\Gamma 的正向为逆时针方向。

求:(1). Γxdx+ydy+zdz\displaystyle \int_\Gamma x\,dx+y\,dy+z\,dz,(2). Γzdx+xdy+ydz\displaystyle \int_\Gamma z\,dx+x\,dy+y\,dz


思路: 该题要结合图形解决,(1). (x,y,z)(x, y, z)τ\vec{\tau} 正交,所以 Γ(x,y,z)τdt=Γ0dt=0\int_\Gamma (x, y, z)\cdot \vec{\tau}\,dt = \int_\Gamma 0\,dt = 0

(2). 利用外积(叉乘)求出 τ=13a(1,1,1)×(x,y,z)\vec{\tau} = \frac{1}{\sqrt{3}a}(1,1,1)\times(x, y, z),代入 Γ(z,x,y)τds\int_\Gamma(z, x, y)\cdot\vec{\tau}\,ds,使用混合积和 xy+xz+yz=12((x+y+z)2x2y2z2)=a22xy+xz+yz = \frac{1}{2}((x+y+z)^2-x^2-y^2-z^2) = -\frac{a^2}{2},即可计算出结果为 3πa2\sqrt{3}\pi a^2

参考作图:

例一

P.S. 这道题还有 StokesStokes 公式的做法,详见 Stokes 公式 - 例一

Green公式

GreenGreen 公式能够将二维的第二型曲线积分(封闭曲线)转化为二重积分(由封闭曲线围成的区域),所以要能从封闭曲线围成的区域判断封闭曲线的定向,封闭曲线相当于区域的边界,判断方法是根据外法向来确定边界的定向(对外法向进行逆时针旋转),下面对这个过程进行详细定义。

定义1(某点处光滑,单位外法向量)

ΩRn\Omega\subset \mathbb R^n 为开集,PΩP\in \partial \Omega,如果 P\exists P 的领域 UU
开集 DRn1D\subset \mathbb R^{n-1}f:DR, fCk, (1k+)f:D\rightarrow \mathbb R,\ f\in C^k,\ (1\leqslant k\leqslant +\infty),使得在适当的坐标系下
P=(x0,f(x0)), x0DP = (x_0, f(x_0)),\ x_0\in D,且

ΩU={(x,y):xD,y=f(x)}ΩU={(x,y):xD,y>f(x)}U\begin{aligned} \partial\Omega\cap U &= \{(x, y):x\in D, y = f(x)\}\\ \Omega\cap U&=\{(x, y):x\in D, y > f(x)\}\cap U \end{aligned}

则称 Ω\partial \OmegaPP 点为 Ck光滑C^k-\text{光滑} (简称光滑),与 ff 的光滑性一致。

x0x_0 处的单位外法向量(也成单位外法向):n=(f(x0),1)1+f(x0)2\vec{n} = \dfrac{(\nabla f(x_0), -1)}{\sqrt{1+|\nabla f(x_0)|^2}},其中 f(x0)\nabla f(x_0)ffx0x_0 处的梯度。

内法向量就是 n=(f(x0),1)1+f(x0)2\vec{n} = \dfrac{(-\nabla f(x_0), 1)}{\sqrt{1+|\nabla f(x_0)|^2}},推导方法:

P=(x0,f(x0))P = (x_0,f(x_0)),则在该坐标系下 PP 点的切线方程为

l:0=f(x0)+f(x0)(xx0)yl:0 = f(x_0)+\triangledown f(x_0)(x-x_0)-y

(x1,y1),(x2,y2)l(x_1,y_1),(x_2,y_2)\in l,则

0= f(x0)+f(x0)(x1x0)y10= f(x0)+f(x0)(x2x0)y2}(f(x0),1)(x1x2,y1y2)=0\left.\begin{aligned} 0=&\ f(x_0)+\triangledown f(x_0)(x_1-x_0)-y_1\\ 0=&\ f(x_0)+\triangledown f(x_0)(x_2-x_0)-y_2\\ \end{aligned}\right\}\Rightarrow (\triangledown f(x_0), -1)\cdot(x_1-x_2,y_1-y_2)=0

看最后一个变量,1-1 则方向向下,是外法向,反之 (f(x0),1)(-\triangledown f(x_0),1) 则是内法向。

单位外法向量场


利用第一维的函数图像来刻画一个点附近的曲线图像,则该点的光滑性就可以用该函数的光滑性来表示。

定义2(平面分段光滑)

ΩR2\Omega\subset \mathbb R^2 为有界区域,Ω\partial \Omega 为有限条内部互不相交的分段光滑的闭曲线,设 Ω=i=1NCi\partial \Omega = \bigcup\limits_{i=1}^NC_iCiC_i 满足:

  1. CiC_i 为光滑曲线。

  2. CiC_i 的内部 CiC_i^\circ 互不相交。

  3. pCi\forall p\in C_i^\circ,有 Ω\partial \OmegaPP 点光滑。

则称 Ω\partial \Omega 分段光滑。

定义3(平面边界的单位内外法向量场,正向)

ΩR2\Omega\subset \mathbb R^2 为有界区域,Ω\partial \Omega 分段光滑,PΩ\forall P\in\partial \OmegaΩ\partial\OmegaPP 点光滑,则称 Ω\partial\Omega 光滑。

n=ΩR2\vec{n}=\partial \Omega\rightarrow \mathbb R^2n(P):=(n1(P),n2(P))\vec{n}(P):=(n_1(P),n_2(P))Ω\partial \OmegaPP 的单位外法向,称为 Ω\partial\Omega 的单位外法向(连续)。

ΩR2\Omega\subset \mathbb R^2 为有界区域,Ω\partial\Omega 分段光滑,Ω=i=1NCi\partial \Omega=\bigcup\limits_{i=1}^NC_i

τ\vec{\tau}CiC_i 的一个定向,PCiP\in C_i^\circ,则

τ(P)=(n2(P),n1(P)), τ=(n2,n1), 逆时针τ(P)=(n2(P),n1(P)), τ=(n2,n1), 顺时针\begin{aligned} \vec{\tau}(P) &= (-n_2(P), n_1(P)),\ \vec{\tau} = (-n_2, n_1),\ \text{逆时针}\\ \vec{\tau}(P) &= (n_2(P), -n_1(P)),\ \vec{\tau} = (n_2, -n_1),\ \text{顺时针} \end{aligned}

τ=(n2,n1)\vec{\tau} = (-n_2, n_1)CiC_i 的正向。

定理4(Green公式)

ΩR2\Omega\subset\mathbb R^2 为有界区域,Ω\partial \Omega 分段光滑,设 P,QC1(Ωˉ)P, Q\in C^1(\bar{\Omega}),则

ΩPdx+Qdy=Ω{QxPy}dxdy\int_{\partial\Omega}P\,dx+Q\,dy = \int_{\Omega}\left\{\frac{\partial Q}{\partial x} - \frac{\partial P}{\partial y}\right\}\,dxdy


证明留到下个 note 了。


第二型曲线积分 Green公式
https://wty-yy.github.io/posts/51465/
作者
wty
发布于
2021年10月30日
许可协议