斯托克斯公式及其应用

文章目录

斯托克斯公式定理定理条件右手规则定理结论和公式👺公式形式的选用

斯托克斯公式和格林公式的关系例例例分析

斯托克斯公式

斯托克斯(Stokes)公式是格林公式的推广格林公式表达了平面闭区域上的二重积分与其边界曲线上的曲线积分之间的关系而斯托克斯公式把曲面

Σ

\Sigma

Σ上的曲面积分与沿曲面

Σ

\Sigma

Σ的边界曲线的曲线积分来联系起来

定理

定理条件

设

Γ

\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

Σ(连同边界

Γ

\Gamma

Γ)上具有一阶连续偏导数,则

∬

Σ

(

R

y

−

Q

z

)

d

y

d

z

+

(

P

z

−

R

x

)

d

z

d

x

+

(

Q

x

−

P

y

)

d

x

d

y

\iint_{\Sigma} (R_{y}-Q_{z})\mathrm{d}y\mathrm{d}z+(P_{z}-R_{x})\mathrm{d}z\mathrm{d}x+(Q_{x}-P_{y})\mathrm{d}x\mathrm{d}y

∬Σ​(Ry​−Qz​)dydz+(Pz​−Rx​)dzdx+(Qx​−Py​)dxdy=

∮

Γ

P

d

x

+

Q

d

y

+

R

d

z

\oint_{\Gamma}{P\mathrm{d}x+Q\mathrm{d}y+R\mathrm{d}z}

∮Γ​Pdx+Qdy+Rdz(1),将其等号左端记为(1-0)

公式(1)两端对调,得

∮

Γ

P

d

x

+

Q

d

y

+

R

d

z

\oint_{\Gamma}{P\mathrm{d}x+Q\mathrm{d}y+R\mathrm{d}z}

∮Γ​Pdx+Qdy+Rdz=

∬

Σ

(

R

y

−

Q

z

)

d

y

d

z

+

(

P

z

−

R

x

)

d

z

d

x

+

(

Q

x

−

P

y

)

d

x

d

y

\iint_{\Sigma} (R_{y}-Q_{z})\mathrm{d}y\mathrm{d}z+(P_{z}-R_{x})\mathrm{d}z\mathrm{d}x+(Q_{x}-P_{y})\mathrm{d}x\mathrm{d}y

∬Σ​(Ry​−Qz​)dydz+(Pz​−Rx​)dzdx+(Qx​−Py​)dxdy(1') 公式(1)左端(即公式(1’)的右端有规律可循

借助式

∮

Γ

P

d

x

+

Q

d

y

+

R

d

z

\oint_{\Gamma}{P\mathrm{d}x+Q\mathrm{d}y+R\mathrm{d}z}

∮Γ​Pdx+Qdy+Rdz的被积表达式:

P

d

x

+

Q

d

y

+

R

d

z

P\mathrm{d}x+Q\mathrm{d}y+R\mathrm{d}z

Pdx+Qdy+Rdz(a)先写出框架

(

)

d

y

d

z

+

(

)

d

z

d

x

+

(

)

d

x

d

y

()\mathrm{d}y\mathrm{d}z+()\mathrm{d}z\mathrm{d}x+()\mathrm{d}x\mathrm{d}y

()dydz+()dzdx+()dxdy(有向曲面元三个分量前加括号)

第一个括号内是式(a)中

d

z

\mathrm{d}z

dz的系数

P

P

P和

d

y

\mathrm{d}y

dy的系数

Q

Q

Q分别对

y

,

z

y,z

y,z求偏导做差.即

R

y

−

Q

z

R_{y}-Q_{z}

Ry​−Qz​类似的填充后面2个括号左后得

(

R

y

−

Q

z

)

d

y

d

z

+

(

P

z

−

R

x

)

d

z

d

x

+

(

Q

x

−

P

y

)

d

x

d

y

(R_{y}-Q_{z})\mathrm{d}y\mathrm{d}z+(P_{z}-R_{x})\mathrm{d}z\mathrm{d}x+(Q_{x}-P_{y})\mathrm{d}x\mathrm{d}y

(Ry​−Qz​)dydz+(Pz​−Rx​)dzdx+(Qx​−Py​)dxdy 但是仍不如下面的形式好记

利用行列式记号(以及微分算子)把斯托克斯公式(1)的左端写成

∬

Σ

∣

d

y

d

z

d

z

d

x

d

x

d

y

∂

∂

x

∂

∂

y

∂

∂

z

P

Q

R

∣

\iint_{\Sigma} \begin{vmatrix} \mathrm{d}y\mathrm{d}z&\mathrm{d}z\mathrm{d}x&\mathrm{d}x\mathrm{d}y\\ \frac{\partial}{\partial{x}}&\frac{\partial}{\partial{y}}&\frac{\partial}{\partial{z}}\\ P&Q&R \end{vmatrix}

∬Σ​

​dydz∂x∂​P​dzdx∂y∂​Q​dxdy∂z∂​R​

​​(1-1)或

∬

Σ

∣

cos

⁡

α

cos

⁡

β

cos

⁡

γ

∂

∂

x

∂

∂

y

∂

∂

z

P

Q

R

∣

⋅

d

S

\iint_{\Sigma} \begin{vmatrix} \cos\alpha&\cos\beta&\cos\gamma\\ \frac{\partial}{\partial{x}}&\frac{\partial}{\partial{y}}&\frac{\partial}{\partial{z}}\\ P&Q&R \end{vmatrix}\cdot\mathrm{d}S

∬Σ​

​cosα∂x∂​P​cosβ∂y∂​Q​cosγ∂z∂​R​

​⋅dS(1-2)应用行列式降阶展开规则(Laplace展开),即得到式(1-1)得到公式(1)的形式;展开时要注意偏微分符号和函数的乘积理解为偏导,例如,将

∂

∂

x

\frac{\partial}{\partial{x}}

∂x∂​和

R

R

R的乘积理解为

∂

R

∂

x

\frac{\partial{R}}{\partial{x}}

∂x∂R​,其他组合类似而公式(1-2)式利用(2)(两类曲面积分间的联系)

d

y

d

z

\mathrm{d}y\mathrm{d}z

dydz=

cos

⁡

α

d

S

\cos\alpha{\mathrm{d}S}

cosαdS

d

z

d

x

\mathrm{d}z\mathrm{d}x

dzdx=

cos

⁡

β

d

S

\cos\beta{\mathrm{d}S}

cosβdS

d

x

d

y

\mathrm{d}x\mathrm{d}y

dxdy=

cos

⁡

γ

d

S

\cos\gamma{\mathrm{d}S}

cosγdS 代换式(1-1),即可得到式(1-2),向量

n

\bold{n}

n=

(

cos

⁡

α

,

cos

⁡

β

,

cos

⁡

γ

)

(\cos\alpha,\cos\beta,\cos\gamma)

(cosα,cosβ,cosγ)为有向曲面

Σ

\Sigma

Σ在

(

x

,

y

,

z

)

(x,y,z)

(x,y,z)处的单位法向量 公式(1)称为斯托克斯公式

式(1-0),(1-1)都是曲面积分中的第二类积分;而式(1-2)则是第一类曲面积分

公式形式的选用

而公式(1-2)在

Γ

\Gamma

Γ张成的

Σ

\Sigma

Σ是平面的情形使用比较合适,这时候计算比较简单,

给出具体平面(方程)时的单位法向量(

cos

⁡

α

,

cos

⁡

β

,

cos

⁡

γ

\cos\alpha,\cos\beta,\cos\gamma

cosα,cosβ,cosγ)是常数 否则

Σ

\Sigma

Σ不是平面,那么求法向量可能就复杂,则应该使用(1-0)或(1-1)

斯托克斯公式和格林公式的关系

若

Σ

\Sigma

Σ是

x

O

y

xOy

xOy面上的一块平面闭区域,斯托克斯公式就变成格林公式因此格林公式是斯托克斯公式的一种特殊情形

例

利用斯托克斯公式计算

I

=

∮

Γ

z

d

x

+

x

d

y

+

y

d

z

I=\oint_{\Gamma}z\mathrm{d}x+x\mathrm{d}y+y\mathrm{d}z

I=∮Γ​zdx+xdy+ydz

其中

Γ

\Gamma

Γ为平面

x

+

y

+

z

=

1

x+y+z=1

x+y+z=1被三个坐标面面所截成的三角形的整个边界

Γ

\Gamma

Γ的正向与这个平面三角形

Σ

\Sigma

Σ上侧的法向量之间符合右手规则 按斯托克斯公式:

I

I

I=

∬

Σ

∣

d

y

d

z

d

z

d

x

d

x

d

y

∂

∂

x

∂

∂

y

∂

∂

z

z

x

y

∣

\iint_{\Sigma} \begin{vmatrix} \mathrm{d}y\mathrm{d}z&\mathrm{d}z\mathrm{d}x&\mathrm{d}x\mathrm{d}y\\ \frac{\partial}{\partial{x}}&\frac{\partial}{\partial{y}}&\frac{\partial}{\partial{z}}\\ z&x&y \end{vmatrix}

∬Σ​

​dydz∂x∂​z​dzdx∂y∂​x​dxdy∂z∂​y​

​=

∬

Σ

d

y

d

z

+

d

z

d

x

+

d

x

d

y

\iint_{\Sigma} \mathrm{d}y\mathrm{d}z+\mathrm{d}z\mathrm{d}x+\mathrm{d}x\mathrm{d}y

∬Σ​dydz+dzdx+dxdy

分项积分,得

I

I

I=

∬

D

y

z

d

σ

\iint_{D_{yz}}\mathrm{d}\sigma

∬Dyz​​dσ+

∬

D

z

x

d

σ

\iint_{D_{zx}}\mathrm{d}\sigma

∬Dzx​​dσ+

∬

D

x

y

d

σ

\iint_{D_{xy}}\mathrm{d}\sigma

∬Dxy​​dσ=

1

2

+

1

2

+

1

2

\frac{1}{2}+\frac{1}{2}+\frac{1}{2}

21​+21​+21​=

3

2

\frac{3}{2}

23​这里

D

y

z

,

D

y

z

,

D

x

y

D_{yz},D_{yz},D_{xy}

Dyz​,Dyz​,Dxy​分别为

Σ

\Sigma

Σ在

y

O

z

,

z

O

x

,

x

O

y

yOz,zOx,xOy

yOz,zOx,xOy面上的投影区域

例

利用斯托克斯公式计算

I

=

∮

Γ

(

y

2

−

z

2

)

d

x

+

(

z

2

−

x

2

)

d

y

+

(

x

2

−

y

2

)

d

z

I=\oint_{\Gamma}(y^2-z^2)\mathrm{d}x+(z^2-x^2)\mathrm{d}y+(x^2-y^2)\mathrm{d}z

I=∮Γ​(y2−z2)dx+(z2−x2)dy+(x2−y2)dz

其中

Γ

\Gamma

Γ是用平面

Π

:

x

+

y

+

z

=

3

2

\Pi:x+y+z=\frac{3}{2}

Π:x+y+z=23​截立方体

Ω

=

{

(

x

,

y

,

z

)

∣

x

∈

[

0

,

1

]

,

y

∈

[

0

,

1

]

,

z

∈

[

0

,

1

]

}

\Omega=\set{(x,y,z)|x\in[0,1],y\in[0,1],z\in[0,1]}

Ω={(x,y,z)∣x∈[0,1],y∈[0,1],z∈[0,1]}的表面所得的截痕若从

O

x

Ox

Ox轴的正向看去,取逆时针方向 分析

平面

Π

\Pi

Π和三个坐标轴的交点坐标分别是

A

(

3

2

,

0

,

0

)

A(\frac{3}{2},0,0)

A(23​,0,0),

B

(

0

,

3

2

,

0

)

B(0,\frac{3}{2},0)

B(0,23​,0),

C

(

0

,

0

,

3

2

)

C(0,0,\frac{3}{2})

C(0,0,23​)并且通过计算可知,原点

O

(

0

,

0

,

0

)

O(0,0,0)

O(0,0,0)到直线

A

B

,

B

C

,

A

C

AB,BC,AC

AB,BC,AC的距离都是

3

2

⋅

2

2

\frac{3}{2}\cdot{\frac{\sqrt2}{2}}

23​⋅22

​​=

3

2

4

\frac{3\sqrt{2}}{4}

432

​​而立方体在三个坐标面到

O

O

O的距离都是

2

\sqrt{2}

2

​;而

(

3

2

4

)

2

(\frac{3\sqrt{2}}{4})^2

(432

​​)2=

9

8

\frac{9}{8}

89​=,而

2

2

=

2

\sqrt{2}^2=2

2

​2=2,显然后者更大,因此

A

B

,

B

C

,

A

C

AB,BC,AC

AB,BC,AC分别穿过立方体在坐标面上的2条棱,将这些交点连结,得到一个六边形(属于平面

Π

\Pi

Π的一部分)

这些交点也可以求,例如求令

y

=

0

,

z

=

1

y=0,z=1

y=0,z=1两平面解平面

Π

\Pi

Π,可以求出

x

=

1

2

x=\frac{1}{2}

x=21​,即点

(

1

2

,

0

,

1

)

(\frac{1}{2},0,1)

(21​,0,1)就是其中一个交点 解

为了利用Stokes公式,我们要曲定一个有向面

这里取有向面

Σ

\Sigma

Σ为

Π

\Pi

Π的上侧被

Γ

\Gamma

Γ所围成的部分区域

Σ

\Sigma

Σ的一个法向量为

(

1

,

1

,

1

)

(1,1,1)

(1,1,1)单位化,得到一个单位法向量:

n

\bold{n}

n=

1

3

(

1

,

1

,

1

)

\frac{1}{\sqrt{3}}(1,1,1)

3

​1​(1,1,1)即

cos

⁡

α

=

cos

⁡

β

=

cos

⁡

γ

\cos\alpha=\cos\beta=\cos\gamma

cosα=cosβ=cosγ=

1

3

\frac{1}{\sqrt{3}}

3

​1​,按stokes公式

I

I

I=

∬

Σ

∣

1

3

1

3

1

3

∂

∂

x

∂

∂

y

∂

∂

z

y

2

−

z

2

z

2

−

x

2

x

2

−

y

2

∣

⋅

d

S

\iint_{\Sigma} \begin{vmatrix} \frac{1}{\sqrt{3}}&\frac{1}{\sqrt{3}}&\frac{1}{\sqrt{3}}\\ \frac{\partial}{\partial{x}}&\frac{\partial}{\partial{y}}&\frac{\partial}{\partial{z}}\\ y^2-z^2&z^2-x^2&x^2-y^2 \end{vmatrix}\cdot\mathrm{d}S

∬Σ​

​3

​1​∂x∂​y2−z2​3

​1​∂y∂​z2−x2​3

​1​∂z∂​x2−y2​

​⋅dS=

1

3

[

(

−

2

y

−

2

z

)

−

(

2

x

+

2

z

)

+

(

−

2

x

−

2

y

)

]

\frac{1}{\sqrt{3}}[(-2y-2z)-(2x+2z)+(-2x-2y)]

3

​1​[(−2y−2z)−(2x+2z)+(−2x−2y)]=

−

4

3

∬

Σ

(

x

+

y

+

z

)

d

S

-\frac{4}{\sqrt{3}}\iint_{\Sigma}(x+y+z)\mathrm{d}S

−3

​4​∬Σ​(x+y+z)dS 利用第一类曲面积分公式

而

Σ

\Sigma

Σ是

Π

\Pi

Π上的一部分,将

Π

\Pi

Π方程变形为

z

=

3

2

−

x

−

y

z=\frac{3}{2}-x-y

z=23​−x−y,所以被积函数为

x

+

y

+

z

x+y+z

x+y+z=

x

+

y

+

(

3

2

−

x

−

y

)

x+y+(\frac{3}{2}-x-y)

x+y+(23​−x−y)=

3

2

\frac{3}{2}

23​

且

G

=

1

+

(

−

1

)

2

+

(

−

1

)

2

G=\sqrt{1+(-1)^2+(-1)^2}

G=1+(−1)2+(−1)2

​=

3

\sqrt 3

3

​

所以:

I

I

I=

−

4

3

∬

D

x

y

3

2

3

d

x

d

y

-\frac{4}{\sqrt{3}}\iint_{D_{xy}}\frac{3}{2}\sqrt{3}\mathrm{d}x\mathrm{d}y

−3

​4​∬Dxy​​23​3

​dxdy=

−

6

∬

D

x

y

d

x

d

y

-6\iint_{D_{xy}}\mathrm{d}x\mathrm{d}y

−6∬Dxy​​dxdy=

−

6

σ

x

y

-6\sigma_{xy}

−6σxy​

或者分步计算

I

I

I=

−

4

3

3

2

∬

Σ

d

S

-\frac{4}{\sqrt{3}}\frac{3}{2}\iint_{\Sigma}\mathrm{d}S

−3

​4​23​∬Σ​dS=

−

2

3

∬

D

x

y

3

d

x

d

y

-2\sqrt{3}\iint_{D_{xy}}\sqrt{3}\mathrm{d}x\mathrm{d}y

−23

​∬Dxy​​3

​dxdy=

−

6

σ

x

y

-6\sigma_{xy}

−6σxy​

其中

D

x

y

D_{xy}

Dxy​为

x

O

y

xOy

xOy面上的投影区域,其面积

σ

x

y

\sigma_{xy}

σxy​=

1

−

2

×

1

8

1-2\times\frac{1}{8}

1−2×81​=

3

4

\frac{3}{4}

43​

将

D

x

y

D_{xy}

Dxy​单独分析:容易确定,

D

x

y

D_{xy}

Dxy​会落在

x

O

y

xOy

xOy面上以

(

0

,

0

,

0

)

,

(

1

,

0

,

0

)

(0,0,0),(1,0,0)

(0,0,0),(1,0,0),

(

1

,

1

,

0

)

,

(

0

,

1

,

0

)

(1,1,0),(0,1,0)

(1,1,0),(0,1,0)的连线而成的正方形

S

1

S_1

S1​内部的一个坐标面内六边形将6顶点分别求出:

(

1

2

,

0

)

(\frac{1}{2},0)

(21​,0),

(

1

,

0

)

(1,0)

(1,0),

(

1

,

1

2

)

(1,\frac{1}{2})

(1,21​),

(

1

2

,

1

)

(\frac{1}{2},1)

(21​,1),

(

0

,

1

)

(0,1)

(0,1),

(

0

,

1

2

)

(0,\frac{1}{2})

(0,21​) 分析其草图,可以用边长为1的正方形

S

1

S_1

S1​的面积减去两个相同的,直角边为

1

2

\frac{1}{2}

21​的等腰直角小三角形的面积得到 所以

I

I

I=

−

6

σ

x

y

{-6}{\sigma_{xy}}

−6σxy​=

−

9

2

-\frac{9}{2}

−29​

例

I

=

∮

Γ

(

y

−

z

)

d

x

+

(

z

−

x

)

d

y

+

(

x

−

y

)

d

z

I=\oint_{\Gamma}(y-z)\mathrm{d}x+(z-x)\mathrm{d}y+(x-y)\mathrm{d}z

I=∮Γ​(y−z)dx+(z−x)dy+(x−y)dz(0)

其中

Γ

\Gamma

Γ是曲线

x

2

+

y

2

=

1

x^2+y^2=1

x2+y2=1(0-1);

x

−

y

+

z

=

2

x-y+z=2

x−y+z=2(0-2)的交线从

z

z

z轴正向(往

z

z

z轴负向)看时,

Γ

\Gamma

Γ的方向是顺时针

分析

使用原始的曲线积分方法计算

需要得到空间有向曲线

Γ

\Gamma

Γ的参数方程形式(设其由两个曲面

Σ

1

,

Σ

2

\Sigma_1,\Sigma_2

Σ1​,Σ2​相交)为了得到参数方程,通常是做

Γ

\Gamma

Γ的投影(到坐标面,例如

x

O

y

xOy

xOy面),先得到投影(记为

C

C

C)的方程

曲线

C

C

C的参数方程通常是容易直接得到然后将得到的两个变量的参数方程代入到

Σ

1

,

Σ

2

\Sigma_1,\Sigma_2

Σ1​,Σ2​(其中的一个),得第三个变量的参数方程 对于参数方程而言,其三个变量的参数方程(设参数为

t

t

t)都是关于参数

t

t

t的式子,因此具有相对独立的特点,可以分批求参数方程对于本例:曲线

Γ

\Gamma

Γ的投影

C

C

C方程为

x

2

+

y

2

=

1

x^2+y^2=1

x2+y2=1其参数方程为

x

=

cos

⁡

θ

x=\cos\theta

x=cosθ,(1-1)

y

=

sin

⁡

θ

y=\sin\theta

y=sinθ;(1-2)

再将两者代入

x

−

y

+

z

=

2

x-y+z=2

x−y+z=2,可得

z

z

z=

2

−

x

+

y

2-x+y

2−x+y=

2

−

cos

⁡

θ

+

sin

⁡

θ

2-\cos\theta+\sin\theta

2−cosθ+sinθ(1-3) 将所以(1-1,1-2,1-3)代入式(0),可以算得

I

I

I=

∫

2

π

0

−

2

(

cos

⁡

θ

+

sin

⁡

θ

)

+

2

cos

⁡

2

θ

+

1

d

θ

\int_{2\pi}^{0}-2(\cos\theta+\sin\theta)+2\cos{2\theta}+1\mathrm{d}\theta

∫2π0​−2(cosθ+sinθ)+2cos2θ+1dθ=

−

2

π

-2\pi

−2π(2)

代入过程中以及化简有一定计算量,注意

∫

0

2

π

\int_{0}^{2\pi}

∫02π​是逆时针圆周的情形,而

∫

2

π

0

\int_{2\pi}^{0}

∫2π0​才是本例的情形掌握几个常用三角定积分特点可以加速某些计算过程,例如

∫

0

2

π

\int_{0}^{2\pi}

∫02π​或

∫

2

π

0

\int_{2\pi}^{0}

∫2π0​或积分区间宽度为被积函数的正数倍时,积分为0,这里被积函数为

sin

⁡

θ

,

cos

⁡

θ

\sin\theta,\cos\theta

sinθ,cosθ时结果总是为0,因为两个函数都是周期为

2

π

2\pi

2π的函数,它们的定积分区间由恰好是周期的整数倍,验证可知,确实如此 使用Stokes公式计算

首先要选择一个

Γ

\Gamma

Γ张成的有向曲面

Σ

\Sigma

Σ,选取尽量简单的面不妨取

Σ

\Sigma

Σ为

Γ

\Gamma

Γ在平面

x

−

y

+

z

=

2

x-y+z=2

x−y+z=2(即平面(0-2))的部分在根据右手法则,结合

Γ

\Gamma

Γ的方向,可定出

Σ

\Sigma

Σ的方向朝下因此化为二重积分时,取符号由stokes公式:

I

I

I=

∬

Σ

(

−

1

+

1

)

d

y

d

z

+

(

1

−

1

)

d

z

d

x

+

(

1

+

1

)

d

x

d

y

\iint_{\Sigma}(-1+1)\mathrm{d}y\mathrm{d}z+(1-1)\mathrm{d}z\mathrm{d}x+(1+1)\mathrm{d}x\mathrm{d}y

∬Σ​(−1+1)dydz+(1−1)dzdx+(1+1)dxdy=

∬

Σ

2

d

x

d

y

\iint_{\Sigma}2\mathrm{d}x\mathrm{d}y

∬Σ​2dxdy=

−

∬

D

x

y

2

d

x

d

y

-\iint_{D_{xy}}2\mathrm{d}x\mathrm{d}y

−∬Dxy​​2dxdy=

−

2

σ

x

y

-2\sigma_{xy}

−2σxy​=

−

2

π

-2\pi

−2π

其中

σ

x

y

\sigma_{xy}

σxy​时

D

x

y

D_{xy}

Dxy​的面积,而

D

x

y

D_{xy}

Dxy​是半径为

1

1

1的圆域,面积为

π

\pi

π 综上,使用stokes公式要简单的多