数学 旧 .com 迁移

数学分析:相关资料 / 曲线积分

从旧 .com 全量搬运的历史内容,来源路径:/math/课程/数学分析/resources/曲线积分/

迁移来源

曲线积分

曲线的弧长 #

definition

对于空间中的参数方程

{x=x(t),y=y(t),t[a,b]z=z(t),\begin{equation} \begin{cases} x=x(t), & \\ y=y(t), &\quad t\in[a,b] \\ z=z(t), & \end{cases} \end{equation}

所定义的曲线段 CC, 如果对任意的 at1<t2ba\leqslant t_1<t_2\leqslant b, 当 t1=at_1=at2=bt_2=b 不同时成立时有

(x(t1),y(t1),z(t1))(x(t2),y(t2),z(t2)),(x(t_1),y(t_1),z(t_1))\neq(x(t_2),y(t_2),z(t_2)),

则称 CC简单曲线. 更进一步的, 如果有 (x(a),y(a),z(a))=(x(b),y(b),z(b))(x(a),y(a),z(a))=(x(b),y(b),z(b)) 则称 CC简单闭曲线.

definition

设曲线段 CC 由 (\ref{空间曲线参数方程}) 所定义. 若存在 sRs\in\mathbb{R}, 使得对任意的 ε>0\varepsilon>0 而言, 存在 δ>0\delta>0, 对由区间 [a,b][a,b] 的任意一组满足 maxiΔti<δ\max\limits_{i}\Delta t_i<\delta 的分点

a=t0<t1<<tn=ba=t_0<t_1<\cdots<t_n=b

所定义的曲线上的点 Mi(x(ti),y(ti),z(ti))M_i(x(t_i),y(t_i),z(t_i)) 均有

1inMi1Mis<ε,\left|\sum\limits_{1\leqslant i\leqslant n}\overline{M_{i-1}M_i}-s\right|<\varepsilon,

那么就称曲线段 CC可求长的, 并称 ssCC弧长.

类似也可以给出由参数方程

{x=x(t),y=y(t),t[a,b]\begin{equation} \begin{cases} x=x(t), \\ y=y(t), \end{cases} \quad t\in[a,b] \end{equation}

所定义的平面上的曲线段及其弧长定义.

tip

CC 是由 (\ref{空间曲线参数方程}) 给出的可求长的曲线段, φ:[c,d][a,b]\varphi:[c,d]\longrightarrow[a,b] 是严格单调的满射, 并记

C1:{x=x(φ(u)),y=y(φ(u)),u[c,d]z=z(φ(u)),C_1:\begin{cases} x=x(\varphi(u)), & \\ y=y(\varphi(u)), &\quad u\in[c,d] \\ z=z(\varphi(u)), & \end{cases}

那么 C1C_1 也是可求长的曲线, 且其弧长等于 CC 的弧长. 简而言之, 曲线的弧长与参数方程的选取无关.

tip

如果 x(t),y(t),z(t)x(t),y(t),z(t) 均在区间 [a,b][a,b] 上连续可导, 则由 (\ref{空间曲线参数方程}) 所定义的曲线段 CC 是可求长的, 且弧长为

s=ab[x(t)]2+[y(t)]2+[z(t)]2dt.s=\int_a^b \sqrt{[x'(t)]^2+[y'(t)]^2+[z'(t)]^2}\t d t.
tip

如果 x(t),y(t)x(t),y(t) 均在区间 [a,b][a,b] 上连续可导, 那么平面上由 (\ref{平面曲线参数方程}) 所定义的曲线段 CC 是可求长的, 且弧长为

s=ab[x(t)]2+[y(t)]2dt.s=\int_a^b \sqrt{[x'(t)]^2+[y'(t)]^2}\t d t.
tip

对于定义在平面上的极坐标方程 r=r(θ) (θ[α,β])r=r(\theta)\ (\theta\in[\alpha,\beta]) 可以将其视作由参数方程

{x=r(θ)cosθ,y=r(θ)sinθ,θ[α,β]\begin{equation} \begin{cases} x=r(\theta)\cos\theta, \\ y=r(\theta)\sin\theta, \\ \end{cases} \quad \theta\in[\alpha,\beta] \end{equation}

那么此时就有

s=αβ[r(θ)]2+[r(θ)]2dθ.s=\int_\alpha^\beta\sqrt{[r'(\theta)]^2+[r(\theta)]^2}\t d\theta.
example

a>0a>0. 对于星形线 (astroid) {x=acos3t,y=asin3t,(t[0,2π])\begin{cases} x=a\cos^3 t, \\ y=a\sin^3 t, \end{cases} (t\in[0,2\pi]) 而言, 其弧长为

\mint[0]2π[x(t)]2+[y(t)]2dt=3a\mint[0]2πcos4tsin2t+sin4t+cos2tdt=3a\mint[0]2πsintcostdt=6a.\begin{aligned} & \mint[0]^{2\pi} \sqrt{[x'(t)]^2+[y'(t)]^2}\t d t \\ &\\ = & 3a\mint[0]^{2\pi} \sqrt{\cos^4 t\sin^2 t+\sin^4 t+\cos^2 t}\t d t \\ & \\ = & 3a\mint[0]^{2\pi}|\sin t\cos t|\t d t = 6a. \end{aligned}
example
tip

简单曲线 CC 的弧长在正交变换下保持不变.

第一型曲线积分 #

definition

CC 是一条可求长的曲线, 其两端点是 AABB (若是闭曲线则 AABB 是一个点), ff 是定义在 CC 上的一个函数. 如果存在实数 II, 使得对任意的 ε>0\varepsilon>0, 均存在 δ>0\delta>0, 当我们依次取分点

A=M0,M1,,Mn=BA=M_0,M_1,\ldots,M_n=B

时, 只要 max1inΔsi<δ\max\limits_{1\leqslant i\leqslant n}\Delta s_i<\delta (其中 Δsi\Delta s_i 表示曲线段 Mi1Mi^\widehat{M_{i-1}M_i} 的弧长), 就对任意的 ξiMi1Mi^\bm\xi_i\in\widehat{M_{i-1}M_i}

i=1nf(ξi)ΔsiI<ε,\left|\sum\limits_{i=1}^n f(\bm\xi_i)\Delta s_i-I\right|<\varepsilon,

那么就称 IIffCC 上的第一型曲线积分(line integral of the first kind), 记作

I=Cf ds.I=\int_C f\ \text{d} s.

特别地, 当 CC 是闭曲线时, 我们也采用记号

I=Cf ds.I=\oint_C f\ \text{d} s.
info

当第一型曲线积分存在时, 积分值与曲线的定向无关.

tip

CC 时一条可求长曲线, ffgg 是定义在 CC 上的两个函数, [leftmargin=1.5cm]

  • (1) 如果 ffggCC 上的第一型曲线积分都存在, 那么对任意的 α,βR\alpha,\beta\in\mathbb{R}, αf+βg\alpha f+\beta gCC 上的第一型曲线积分存在并且,
C(αf+βg)ds=αCfds+βCgds.\int_C (\alpha f + \beta g)\text{d} s=\alpha\int_C f\text{d} s + \beta\int_C g\text{d} s.
  • (2) 如果 C=C1C2C=C_1\cup C_2, C1,C2C_1,C_2 均是可求长曲线, 且公共点为端点, 那么当 C1,C2C_1,C_2 的第一型曲线积分都存在时, ffCC 上的第一型曲线积分也存在, 且
Cfds=C1fds+C2fds.\int_C f\text{d} s=\int_{C_1} f\text{d} s+\int_{C_2} f\text{d} s.
definition

CCR3\mathbb{R}^3 中的光滑曲线段, 即存在参数方程

{x=x(t),y=y(t),z=z(t),t[a,b]\begin{cases} x=x(t),\\ y=y(t),\\ z=z(t), \end{cases} t\in [a,b]

表示 CC, 且 x(t),y(t),z(t)x(t),y(t),z(t) 均在 [a,b][a,b] 上连续可微.

取分点 ,求黎曼和, 用积分第一中值定理及闵可夫斯基不等式进行等价, 可得上述光滑曲线段的第一型曲线积分为

Cf(x,y,z)ds=abf(x(t),y(t),z(t))[x(t)]2+[y(t)]2+[z(t)]2dt.\begin{equation} \int_C f(x,y,z)\text{d} s=\int_a^b f(x(t),y(t),z(t)) \sqrt{[x'(t)]^2+[y'(t)]^2+[z'(t)]^2} \text{d} t. \end{equation}

类似地, 如果是平面上的曲线, 则有

Cf(x,y)ds=abf(x(t),y(t))[x(t)]2+[y(t)]2dt.\begin{equation} \int_C f(x,y)\text{d} s=\int_a^b f(x(t),y(t)) \sqrt{[x'(t)]^2+[y'(t)]^2} \text{d} t. \end{equation}

第二型曲线积分 #

definition

CCR3\mathbb{R}^3 中的一条定向的可求长的曲线, 起点AA, 终点BB, 在 CC 上定义映射 f=(P,Q,R)T:CR3f=(P,Q,R)^T:C\longrightarrow \mathbb{R}^3. 若存在实数 II, 使得对任意的 ε>0\varepsilon>0, 均存在 δ>0\delta>0, 当我们在 CC 上从 AABB 依次取分点

A=M0,M1,,Mn=BA=M_0,M_1,\ldots,M_n=B

时, 只要 max1inMi1Mi<δ\max\limits_{1\leqslant i\leqslant n}\overline{M_{i-1}M_i}<\delta, 就对任意的 ξiMi1Mi^\bm\xi_i\in\widehat{M_{i-1}M_i}

i=1n<f(ξi),Mi1Mi>I<ε,\left|\sum\limits_{i=1}^n\left<f(\bm\xi_i),\overrightarrow{M_{i-1}M_i}\right>-I\right|<\varepsilon,

则称 IIf=(P,Q,R)Tf=(P,Q,R)^T 沿定向曲线 CC第二型曲线积分 (line integral of the second kind). 也称作 ff 沿道路 AB^\widehat{AB}第二型曲线积分, 记作

I=CPdx+Qdy+Rdz=AB^Pdx+Qdy+Rdz.I=\int_C P\text{d} x+Q\text{d} y+R\text{d} z=\int_{\widehat{AB}}P\text{d} x+Q\text{d} y+R\text{d} z.

特别地, 当 CC 是闭曲线时, 我们也采用记号

I=CPdx+Qdy+Rdz.I=\oint_C P\text{d} x+Q\text{d} y+R\text{d} z.

类似可定义 R2\mathbb{R}^2 中定向曲线 CC 的第二型曲线积分

CPdx+Qdy.\int_C P\text{d} x+Q\text{d} y.
info

在计算第二型曲线积分时, 需注意曲线的定向, 因为对于以 A,BA,B 为端点的曲线

AB^Pdx+Qdy+Rdz=BA^Pdx+Qdy+Rdz\int_{\widehat{AB}}P\text{d} x+Q\text{d} y+R\text{d} z=-\int_{\widehat{BA}}P\text{d} x+Q\text{d} y+R\text{d} z
tip

AB^\widehat{AB}R3\mathbb{R}^3 中的一条可求长的定向曲线, f=(P1,Q1,R1)Tf=(P_1,Q_1,R_1)^Tg=(P2,Q2,R2)Tg=(P_2,Q_2,R_2)^T 均是从 AB^\widehat{AB}R3\mathbb{R}^3 的映射. [leftmargin=1.5cm]

  • (1) 若 ff,gg 沿 AB^\widehat{AB} 的第二型曲线积分均存在, 则对任意的 α,βR, αf+βg\alpha,\beta\in \mathbb{R},\ \alpha f+\beta g 沿 AB^\widehat{AB} 的第二型曲线积分也存在, 并且等于
α(AB^P1dx+Q1dy+R1dz)+β(AB^P2dxQ2dyR2dz).\alpha\left(\int_{\widehat{AB}} P_1\text{d} x+Q_1\text{d} y+R_1\text{d} z\right)+\beta\left(\int_{\widehat{AB}}P_2\text{d} xQ_2\text{d} yR_2\text{d} z\right).
  • (2) 设 DDAB^\widehat{AB} 上一点, 如果 ff 沿 AD^\widehat{AD}DB^\widehat{DB} 的第二型曲线积分均存在, 则 ff 沿 AB^\widehat{AB} 的第二型曲线积分也存在, 并且等于
AD^P1dx+Q1dy+R1dz+DB^P1dx+Q1dy+R1dz.\int_{\widehat{AD}} P_1\text{d} x+Q_1\text{d} y+R_1\text{d} z+\int_{\widehat{DB}} P_1\text{d} x+Q_1\text{d} y+R_1\text{d} z.

AB^\widehat{AB}R3\mathbb{R}^3 中的定向光滑曲线段, 再设

f(P,Q,R)T:AB^R3.f(P,Q,R)^T:\widehat{AB}\longrightarrow \mathbb{R}^3.

则有

\mint[AB^]Pdx+Qdy+Rdz=\mint[a]b[P(x(t),y(t),z(t))x(t)+Q(x(t),y(t),z(t))y(t)+R(x(t),y(t),z(t))z(t)]dt\begin{aligned} \mint[\widehat{AB}] P\text{d} x+Q\text{d} y+R\text{d} z=&\mint[a]^b[P(x(t),y(t),z(t))x'(t)+Q(x(t),y(t),z(t))y'(t)\\ &+R(x(t),y(t),z(t))z'(t)]\text{d} t \end{aligned}

格林公式 #

definition

对于 R2\mathbb{R}^2 平面上的有界闭区域 DD, 其边界 D\partial D, 是由有限条光滑曲线组成. 当在边界上行走时, 如果与之相邻的区域的内部总是在左侧, 则称这个方向是正向

tip

SSR2\mathbb{R}^2 中的有界闭区域, S\partial S 由有限多条分段光滑曲线组成, 若 P,QC1(S)P,Q\in C^1(S), 则

SPdx+Qdy=S(QxPy)dxdy\begin{equation} \int_{\partial S}P\text{d} x+Q\text{d} y = \iint\limits_S\left(\frac{\partial Q}{\partial x}-\frac{\partial P}{\partial y}\right)\text{d} x\text{d} y \end{equation}

其中 S\partial S 的定向为正向.

在定理 \ref{格林公式} 条件下, 再设 u(x,y)u(x,y)SS 上连续可微, 那么将 (\ref{格林公式1}) 中的 PP 换为 uPuP, 并取 Q=0Q=0 可得

SuPdx=S(uP)ydxdy=S(Puy+uPy)dxdy,\begin{equation*} \int_{\partial S} uP\text{d} x = -\iint\limits_{S}\frac{\partial(uP)}{\partial y}\text{d} x\text{d} y = -\iint\limits_{S}\left(P\frac{\partial u}{\partial y}+u\frac{\partial P}{\partial y}\right) \text{d} x\text{d} y, \end{equation*}

也即

SuPydxdy=SuPdx+SPuydxdy.\begin{equation} -\iint\limits_{S}u\frac{\partial P}{\partial y}\text{d} x\text{d} y=\int_{\partial S}uP\text{d} x+\iint\limits_{S} P\frac{\partial u}{\partial y}\text{d} x\text{d} y. \end{equation}

同理, 将 QQ 换为 uQuQ 可得,

SuQxdxdy=SuQdySQuxdxdy.\begin{equation} \iint\limits_{S}u\frac{\partial Q}{\partial x}\text{d} x\text{d} y=\int_{\partial S}uQ\text{d} y-\iint\limits_{S} Q\frac{\partial u}{\partial x}\text{d} x\text{d} y. \end{equation}

相加后可得,

Su(QxPy)dxdy=(SuPdx+uQdy)S(QuxPuy)dxdy.\begin{equation} \iint\limits_{S}u\left(\frac{\partial Q}{\partial x}-\frac{\partial P}{\partial y}\right)\text{d} x\text{d} y=\left(\int_{\partial S}uP\text{d} x + uQ\text{d} y\right)-\iint\limits_{S}\left( Q\frac{\partial u}{\partial x}-P\frac{\partial u}{\partial y}\right)\text{d} x\text{d} y. \end{equation}

以上三式均被称作平面上的分部积分公式.

definition

对于 R2\mathbb{R}^2 中的一个区域 DD, 若 DD 中任意一条简单闭曲线所围成的区域均包含于 DD, 则称 DD单连通的 (simply connected), 否则称 DD多连通的 (multiply connected) 或者称作复连通的.

tip

利用格林公式计算闭曲线围成的面积. 设 SSR2\mathbb{R}^2 中的一个有界闭区域, 且 S\partial S 由有限多条光滑曲线组成, 那么由格林公式

μ(S)=Sdxdy=Sxdy=Sydx.\begin{equation} \mu(S)=\iint\limits_{S}\text{d} x\text{d} y = \int_{\partial S}x\text{d} y=-\int_{\partial S}y \text{d} x. \end{equation}

更进一步的, 有

μ(S)=12Sxdyydx.\begin{equation} \mu(S)=\frac 1 2\int_{\partial S}x\text{d} y-y\text{d} x. \end{equation}

虽然看上去 (\ref{格林公式面积2}) 和 (\ref{格林公式面积1}) 没有实质上的差异. 但在实际计算中, 如果曲线有一定的对称性 (\ref{格林公式面积2}) 能带来很大的便利.

tip

DDR2\mathbb{R}^2 中的一个单连通区域, P,QC1(D)P,Q\in C^1(D), 则下列命题等价: [leftmargin=1.5cm]

  • (1) 对 DD 中任意两点 A,BA,B 以及 DD 中从 AABB 的任意两条分段光滑曲线 C1,C2C_1,C_2
C1Pdx+Qdy=C2Pdx+dy.\int_{C_1} P\text{d} x+Q\text{d} y = \int_{C_2}P\text{d} x+\text{d} y.

第二型曲线积分与路径无关.

  • (2) 对于 DD 中由有限多条光滑曲线组成的任一闭曲线 CC
CPdx+Qdy=0.\int_C P\text{d} x+Q\text{d} y = 0.
  • (3) 在 DD 上有 Py=Qx\dfrac{\partial P}{\partial y} = \dfrac{\partial Q}{\partial x}.

应用: 调和函数 #

definition

DD 是一个平面 (闭) 区域, ff 是定义在 DD 上的具有二阶偏导数的函数, 若在 DD 上有

2fx2+2fy2=0,\frac{\partial^2 f}{\partial x^2}+\frac{\partial^2 f}{\partial y^2} = 0,

则称 ffDD 上的调和函数 (harmonic function).

通常记

Δf=2fx2+2fy2,\Delta f = \frac{\partial^2 f}{\partial x^2}+\frac{\partial^2 f}{\partial y^2},

并称 Δ=2x2+2y2\Delta = \dfrac{\partial^2}{\partial x^2}+\dfrac{\partial^2}{\partial y^2}拉普拉斯算子 (Laplace operator).

abstract

ffC2C^2 类的调和函数, A=[abcd]\\ \\A=\left[\begin{aligned} a & b \\ c & d \end{aligned}\right] 是一个\hr{正交矩阵}, g(x,y)=f(ax+by,cx+dy)g(x,y)=f(ax+by,cx+dy). 则有 Δf=Δg\Delta f = \Delta g.

note

x=ax+by,y=cx+dyx'=ax+by,y'=cx+dy, 利用偏导数的链式法则可得

\begin{equation*} \begin{aligned} \Delta g &=& \dfrac{\partial^2 g}{\partial x^2} + \dfrac{\partial^2 g}{\partial y^2} \

10pt] &=& \dfrac{\partial }{\partial x}\left(\dfrac{\partial g}{\partial x’}\cdot\dfrac{\partial x’}{\partial x}+\dfrac{\partial g}{\partial y’}\cdot\dfrac{\partial y’}{\partial x}\right)+\dfrac{\partial }{\partial y}\left(\dfrac{\partial g}{\partial x’}\cdot\dfrac{\partial x’}{\partial y}+\dfrac{\partial g}{\partial y’}\cdot\dfrac{\partial y’}{\partial y}\right) \

10pt] &=&\dfrac{\partial }{\partial x}\left(\dfrac{\partial f}{\partial x'}\cdot a+\dfrac{\partial f}{\partial y'}\cdot c\right)+\dfrac{\partial }{\partial y}\left(\dfrac{\partial f}{\partial x'}\cdot b+\dfrac{\partial f}{\partial y'}\cdot d\right) \

10pt] &=&\dfrac{\partial }{\partial x’}\left(\dfrac{\partial f}{\partial x’}\cdot a+\dfrac{\partial f}{\partial y’}\cdot c\right)\dfrac{\partial x’}{\partial x}+\dfrac{\partial }{\partial y’}\left(\dfrac{\partial f}{\partial x’}\cdot a+\dfrac{\partial f}{\partial y’}\cdot c\right)\dfrac{\partial y’}{\partial x} \

10pt] &&+\dfrac{\partial }{\partial x'}\left(\dfrac{\partial f}{\partial x'}\cdot b+\dfrac{\partial f}{\partial y'}\cdot d\right)\dfrac{\partial x'}{\partial y}+\dfrac{\partial }{\partial y'}\left(\dfrac{\partial f}{\partial x'}\cdot b+\dfrac{\partial f}{\partial y'}\cdot d\right)\dfrac{\partial y'}{\partial y} \

10pt] &=& (a^2+b^2)\dfrac{\partial^2 f}{\partial {x’}^2}+(c^2+d^2)\dfrac{\partial^2 f}{\partial {y’}^2}+(ac+ac+bd+bd)\dfrac{\partial^2 f}{\partial x’\partial y’}\

10pt] &=& \dfrac{\partial^2 f}{\partial {x'}^2}+\dfrac{\partial^2 f}{\partial {y'}^2} = \Delta f. \end{aligned} \end{equation*}

上述最后一行利用了正交矩阵的性质, 任意两行向量点积是 00, 即 ac+bd=0ac+bd=0.

更进一步的, 如果是 nn 元调和函数 g(x1,x2,,xn)=f(x1,x2,,xn)g(x_1,x_2,\ldots,x_n)=f(x_1',x_2',\ldots,x_n') 其中, (x1,,xn)T=A(x1,,xn)T(x_1',\ldots,x_n')^T=A(x_1,\ldots,x_n)^T, 且 AAnn 阶正交矩阵.

那么有 xixj=ai,j\dfrac{\partial x_i'}{\partial x_j}=a_{i,j}.

Δg=i=1n2gxi2[10pt]=i=1nxi(j=1ngxjxjxi)[10pt]=i=1nk=1nxk(j=1nfxjaj,i)xkxi[10pt]=i=1nk=1nxk(j=1nfxjaj,i)ak,i[10pt]=k=1nj=1ni=1nak,iaj,i2fxkxj[10pt]=k=1nfxk2=Δf.\begin{equation*} \begin{aligned} \Delta g &=&\sum\limits_{i=1}^n \dfrac{\partial^2 g}{\partial x_i^2}\\ [10pt] &=& \sum\limits_{i=1}^n \dfrac{\partial}{\partial x_i}(\sum\limits_{j=1}^n\dfrac{\partial g}{\partial x_j'}\cdot\dfrac{\partial x_j'}{\partial x_i}) \\ [10pt] &=& \sum\limits_{i=1}^n \sum\limits_{k=1}^n \dfrac{\partial}{\partial x_k'}(\sum\limits_{j=1}^n\dfrac{\partial f}{\partial x_j'}\cdot a_{j,i})\dfrac{\partial x_k'}{\partial x_i} \\ [10pt] &=& \sum\limits_{i=1}^n \sum\limits_{k=1}^n \dfrac{\partial}{\partial x_k'}(\sum\limits_{j=1}^n\dfrac{\partial f}{\partial x_j'}\cdot a_{j,i}) a_{k,i} \\ [10pt] &=&\sum\limits_{k=1}^n\sum\limits_{j=1}^n\sum\limits_{i=1}^n a_{k,i}a_{j,i}\dfrac{\partial^2 f}{\partial x_k'\partial x_j'} \\ [10pt] &=&\sum\limits_{k=1}^n \dfrac{\partial^f}{\partial x_k'^2}=\Delta f. \end{aligned} \end{equation*}

利用到了 i=1naj,iak,i={1,j=k0,jk\sum\limits_{i=1}^n a_{j,i}a_{k,i}=\begin{cases} 1,\quad j=k \\ 0,\quad j\neq k \end{cases}

tip

DD 是平面上由有限多条光滑曲线所围城的有界闭区域, uuvv 是定义在 DD 上的两个函数, 且 uC2(D), vC1(D)u\in C^2(D),\ v \in C^1(D), 则

DvΔudxdy=D(ux vx+uy vy)\iint\limits_{D}v\Delta u\text{d} x\text{d} y = -\int\limits_{D}\left(\frac{\partial u}{\partial x}\ \frac{\partial v}{\partial x}+\frac{\partial u}{\partial y}\ \frac{\partial v}{\partial y}\right)

讨论

评论

正在加载评论...