過去研究的 $f(x)$ 或 $f(x,y)$ 把點映成一個數值;而向量值函數 $\mathbf{r}(t)$ 把參數映成一個向量。本節要介紹的向量場同時保留兩者最有用的部分:把平面或空間中的每一個點對應到一個向量,因此可以表示像力場、速度場、電場、磁場這類「每個位置都有方向 + 大小」的物理量。
在平面區域 $R$ 上的向量場是一個函數 $\mathbf{F}$,對 $R$ 內每個點都對應一個向量 $\mathbf{F}(x,y)$。
在空間區域 $Q$ 上的向量場是一個函數 $\mathbf{F}$,對 $Q$ 內每個點都對應一個向量 $\mathbf{F}(x,y,z)$。
用分量寫成 $\mathbf{F}(x,y)=M(x,y)\,\mathbf{i}+N(x,y)\,\mathbf{j}$ 或 $\mathbf{F}(x,y,z)=M\,\mathbf{i}+N\,\mathbf{j}+P\,\mathbf{k}$,當每個分量函數連續時稱 $\mathbf{F}$ 連續。
給定純量函數 $f(x,y)$,
$\nabla f = f_x\,\mathbf{i}+f_y\,\mathbf{j}$ 在每點指向「上升最快」的方向。例如 $f=x^2y+3xy^3$,則 $\nabla f=(2xy+3y^3)\,\mathbf{i}+(x^2+9xy^2)\,\mathbf{j}$。
$\mathbf{F}=\dfrac{k}{\|\mathbf{r}\|^2}\mathbf{u}$,方向始終指向(或背離)原點,且大小只與距離有關。Newton 重力場 $-\dfrac{Gm_1m_2}{\|\mathbf{r}\|^2}\mathbf{u}$、Coulomb 電場都屬此類。
想像繞 $z$ 軸轉動的輪子:每個點的速度方向切於圓周,且距軸越遠速度越大。典型表達式 $\mathbf{F}=-y\,\mathbf{i}+x\,\mathbf{j}$,整體呈逆時針渦流。
設 $\mathbf{r}(t)=x(t)\,\mathbf{i}+y(t)\,\mathbf{j}+z(t)\,\mathbf{k}$ 為位置向量,若向量場可寫成 $$\mathbf{F}(x,y,z)=\dfrac{k}{\|\mathbf{r}\|^2}\,\mathbf{u},\quad \mathbf{u}=\dfrac{\mathbf{r}}{\|\mathbf{r}\|}$$ 則稱 $\mathbf{F}$ 為反平方場;重力場與電場都是反平方場的例子。
向量場由無窮多個向量組成,圖示時的目標是挑代表性向量畫出來,看出方向與大小的分布規律。下面互動圖把六種典型向量場放進同一個介面,可以即時切換、調整格點密度、看 $|\mathbf{F}|$ 熱度圖,並比較 $\partial N/\partial x$ 與 $\partial M/\partial y$ 是否相等(即定理 15.1 的保守判定)。
令 $\|\mathbf{F}\|=c$,得 $\sqrt{4x^2+y^2}=c$,即橢圓 $\dfrac{x^2}{(c/2)^2}+\dfrac{y^2}{c^2}=1$。
同一條等大線上的向量「長度相同」,方向才會變。
取 $c=1, 2, 4$ 三圈橢圓上的點代入 $\mathbf{F}=(2x,\,y)$:
| 點 (x,y) | F = (2x, y) | ‖F‖ |
|---|---|---|
| (1/2, 0) | (1, 0) | 1 |
| (0, 1) | (0, 1) | 1 |
| (1, 0) | (2, 0) | 2 |
| (0, 2) | (0, 2) | 2 |
| (2, 0) | (4, 0) | 4 |
| (1, 2√3) | (2, 2√3) | 4 |
例題 2 畫出的箭頭看起來「總是垂直於某條等高線」,這正是梯度的特徵。並非所有向量場都能寫成某個 $f$ 的梯度 —— 可以的叫做保守,那個 $f$ 叫做位能函數。
若存在可微函數 $f$ 使得 $\mathbf{F}=\nabla f$,則稱向量場 $\mathbf{F}$ 為保守,並稱 $f$ 為 $\mathbf{F}$ 的位能函數。
設 $M$、$N$ 在開圓盤 $R$ 上具連續一階偏導,則 $\mathbf{F}(x,y)=M\,\mathbf{i}+N\,\mathbf{j}$ 在 $R$ 上為保守的充要條件是 $$\dfrac{\partial N}{\partial x}=\dfrac{\partial M}{\partial y}.$$
取 $f(x,y)=x^2+\tfrac{1}{2}y^2$,則 $f_x=2x$、$f_y=y$,正好等於 $\mathbf{F}$ 的兩個分量。所以 $\mathbf{F}=\nabla f$,由定義 $\mathbf{F}$ 是保守的。
對 $\mathbf{F}=\dfrac{k}{\|\mathbf{r}\|^2}\mathbf{u}=\dfrac{k(x\mathbf{i}+y\mathbf{j}+z\mathbf{k})}{(x^2+y^2+z^2)^{3/2}}$,取 $f(x,y,z)=\dfrac{-k}{\sqrt{x^2+y^2+z^2}}=\dfrac{-k}{\|\mathbf{r}\|}$。
驗算:$f_x=\dfrac{kx}{\|\mathbf{r}\|^3}$,同樣對 $y$、$z$,所以 $\nabla f = \mathbf{F}$。
$M=x^2y$、$N=xy$。
$\dfrac{\partial M}{\partial y}=x^2$、$\dfrac{\partial N}{\partial x}=y$。
兩者不等(在 $x^2=y$ 的曲線上才相等,但不是整個 $R$ 上),故 $\mathbf{F}$ 不保守。
$M=2x$、$N=y$。
$\dfrac{\partial M}{\partial y}=0$、$\dfrac{\partial N}{\partial x}=0$,兩者相等於整個 $\mathbb{R}^2$,故 $\mathbf{F}$ 保守。
$\dfrac{\partial M}{\partial y}=2x=\dfrac{\partial N}{\partial x}$,確認保守。
$f(x,y)=\displaystyle\int 2xy\,dx = x^2 y + g(y)$,其中 $g(y)$ 為僅與 $y$ 有關的待定函數。
對 $x^2 y + g(y)$ 對 $y$ 偏導:$f_y = x^2 + g'(y)$。
令 $x^2 + g'(y) = x^2 - y$ ⇒ $g'(y) = -y$ ⇒ $g(y) = -\dfrac{y^2}{2}+K$。
把定理 15.1 推到三維,需要先定義旋度,它本身是一個向量場;以及純量的散度。直觀上 curl 描述「在該點繞著轉動」的強度與軸向,div 描述「在該點向外發散」的速率。
設 $\mathbf{F}(x,y,z)=M\,\mathbf{i}+N\,\mathbf{j}+P\,\mathbf{k}$,則 $$\mathrm{curl}\,\mathbf{F}=\nabla\times\mathbf{F}=\left(\dfrac{\partial P}{\partial y}-\dfrac{\partial N}{\partial z}\right)\mathbf{i}-\left(\dfrac{\partial P}{\partial x}-\dfrac{\partial M}{\partial z}\right)\mathbf{j}+\left(\dfrac{\partial N}{\partial x}-\dfrac{\partial M}{\partial y}\right)\mathbf{k}.$$
記憶要點:i, j, k 三個分量裡,j 那一項要記得「先 P 後 M」並乘上負號,等價於下方行列式展開。若 $\mathrm{curl}\,\mathbf{F}=\mathbf{0}$,稱 $\mathbf{F}$ 為無旋的。
若 $M$、$N$、$P$ 在開球 $Q$ 上具連續一階偏導,則 $\mathbf{F}=M\,\mathbf{i}+N\,\mathbf{j}+P\,\mathbf{k}$ 保守 ⇔ $\mathrm{curl}\,\mathbf{F}=\mathbf{0}$,即 $$\dfrac{\partial P}{\partial y}=\dfrac{\partial N}{\partial z},\quad \dfrac{\partial P}{\partial x}=\dfrac{\partial M}{\partial z},\quad \dfrac{\partial N}{\partial x}=\dfrac{\partial M}{\partial y}.$$
$M=2xy,\ N=x^2+z^2,\ P=2yz$。
偏導:$M_y=2x,\ M_z=0;\ N_x=2x,\ N_z=2z;\ P_x=0,\ P_y=2z$。
$\mathbf{i}$ 分量:$P_y-N_z = 2z - 2z = 0$。
$\mathbf{j}$ 分量:$-(P_x-M_z) = -(0-0) = 0$。
$\mathbf{k}$ 分量:$N_x-M_y = 2x - 2x = 0$。
$f = \int 2xy\,dx = x^2 y + g(y,z)$。
$f_y = x^2 + g_y(y,z) = x^2 + z^2$ ⇒ $g_y = z^2$ ⇒ $g(y,z) = y z^2 + h(z)$。
所以 $f = x^2 y + y z^2 + h(z)$。
$f_z = 2 y z + h'(z) = 2 y z$ ⇒ $h'(z) = 0$ ⇒ $h(z) = K$。
$\mathrm{div}\,\mathbf{F}(x,y) = \nabla\cdot\mathbf{F} = \dfrac{\partial M}{\partial x}+\dfrac{\partial N}{\partial y}$(平面)
$\mathrm{div}\,\mathbf{F}(x,y,z) = \nabla\cdot\mathbf{F} = \dfrac{\partial M}{\partial x}+\dfrac{\partial N}{\partial y}+\dfrac{\partial P}{\partial z}$(空間)
若 $\mathrm{div}\,\mathbf{F}=0$,稱 $\mathbf{F}$ 為無散度的。流體力學裡叫不可壓縮,電磁學裡叫螺線型場。
$\dfrac{\partial M}{\partial x}=3 x^2 y^2 z,\quad \dfrac{\partial N}{\partial y}=0,\quad \dfrac{\partial P}{\partial z}=0$。
$\mathrm{div}\,\mathbf{F} = 3 x^2 y^2 z = 3 \cdot 4 \cdot 1 \cdot (-1) = -12$。
若 $\mathbf{F}=M\,\mathbf{i}+N\,\mathbf{j}+P\,\mathbf{k}$ 的三個分量都具有連續二階偏導,則 $$\mathrm{div}(\mathrm{curl}\,\mathbf{F})=0.$$ 證明來自 Clairaut 定理(混合偏導相等),兩兩抵消。
到目前為止接觸的積分都「沿一段區間 $[a,b]$」或「沿一塊平面區域 $R$」。這一節要把積分的領域換成一條曲線 $C$,得到所謂的線積分。它有兩種類型:純量被積函數 $\int_C f\,ds$(例如求帶質量密度的細線總質量)與向量被積函數 $\int_C \mathbf{F}\cdot d\mathbf{r}$(例如算力場做的功)。
平面曲線 $C:\mathbf{r}(t)=x(t)\,\mathbf{i}+y(t)\,\mathbf{j}$,$a\le t\le b$,光滑表示 $dx/dt,\,dy/dt$ 連續且不同時為零。空間曲線多一條 $dz/dt$。
若 $[a,b]$ 可被分割成有限段,每段都光滑,則 $C$ 為分段光滑。
$C_1:\ (0,0,0)\to(1,0,0)$ $\mathbf{r}_1(t) = t\,\mathbf{i},\ 0\le t\le 1$
$C_2:\ (1,0,0)\to(1,2,0)$ $\mathbf{r}_2(t) = \mathbf{i} + (t-1)\,2\,\mathbf{j},\ 1\le t\le 2$
$C_3:\ (1,2,0)\to(1,2,1)$ $\mathbf{r}_3(t) = \mathbf{i} + 2\,\mathbf{j} + (t-2)\,\mathbf{k},\ 2\le t\le 3$
設 $f$ 在含光滑、有限長度曲線 $C$ 的區域上有定義,把 $C$ 分割為 $n$ 段,每段弧長為 $\Delta s_i$,取樣點 $(x_i, y_i)$(平面)或 $(x_i, y_i, z_i)$(空間)。當 $\|\Delta\|\to 0$:
若極限存在則稱 $f$ 沿 $C$ 可積。可證:當 $f$ 連續時,極限對所有合理的光滑參數化都得到相同的值。
若 $f$ 連續、$C:\mathbf{r}(t)=x(t)\mathbf{i}+y(t)\mathbf{j}$($a\le t\le b$)光滑,則 $$\int_C f(x,y)\,ds=\int_a^b f(x(t),y(t))\sqrt{[x'(t)]^2+[y'(t)]^2}\,dt.$$ 空間版本 $\mathbf{r}(t)=x(t)\mathbf{i}+y(t)\mathbf{j}+z(t)\mathbf{k}$ 為 $$\int_C f(x,y,z)\,ds=\int_a^b f(x(t),y(t),z(t))\sqrt{[x'(t)]^2+[y'(t)]^2+[z'(t)]^2}\,dt.$$
特例:若 $f\equiv 1$,線積分給出曲線的弧長 $\int_C ds = \int_a^b \|\mathbf{r}'(t)\|\,dt$。
下面的動畫把定理 15.4 視覺化:選曲線 $C$ 與被積函數 $f$,按 ▶ 播放,紅色標記沿曲線移動;右側即時列出目前位置 $\mathbf{r}(t)$、弧長因子 $\|\mathbf{r}'(t)\|$、被積函數值 $f(\mathbf{r}(t))$ 與累積量 $\int_{t_{\min}}^t f(\mathbf{r}(s))\|\mathbf{r}'(s)\|\,ds$。進度條只表示動畫時間進度,真正的積分大小請看下方數值與目標值。
$\mathbf{r}(t) = t\,\mathbf{i} + 2t\,\mathbf{j} + t\,\mathbf{k},\quad 0\le t\le 1$。
$\mathbf{r}'(t) = (1, 2, 1)$,$\|\mathbf{r}'(t)\| = \sqrt{1+4+1} = \sqrt{6}$。
$f(x,y,z) = x^2 - y + 3z$ 沿 $C$:$f(t, 2t, t) = t^2 - 2t + 3t = t^2 + t$。
$\displaystyle\int_C f\,ds = \int_0^1 (t^2 + t)\sqrt{6}\,dt = \sqrt{6}\left[\dfrac{t^3}{3}+\dfrac{t^2}{2}\right]_0^1 = \sqrt{6}\left(\dfrac{1}{3}+\dfrac{1}{2}\right) = \dfrac{5\sqrt{6}}{6}.$
$\mathbf{r}_1(t) = t\,\mathbf{i},\ 0\le t\le 1$;$\|\mathbf{r}_1'\| = 1$。
$\displaystyle\int_{C_1} x\,ds = \int_0^1 t\cdot 1\,dt = \dfrac{1}{2}$。
$\mathbf{r}_2(t) = \cos t\,\mathbf{i} + \sin t\,\mathbf{j},\ 0\le t\le \pi/2$;$\|\mathbf{r}_2'\| = 1$。
$\displaystyle\int_{C_2} x\,ds = \int_0^{\pi/2}\cos t\,dt = 1$。
$\mathbf{r}'(t) = (1,\,2t^{1/2},\,t)$; $$\|\mathbf{r}'(t)\| = \sqrt{1+4t+t^2}.$$ 這個根號不需要先化成完全平方;因為被積函數是 $x+2=t+2$,正好可以和 $1+4t+t^2$ 的導數配合做變數變換。
$x = t$,所以 $f = t+2$。
$\displaystyle\int_C (x+2)\,ds=\int_0^2 (t+2)\sqrt{1+4t+t^2}\,dt.$
令 $u = 1+4t+t^2$,$du = (4+2t)\,dt = 2(t+2)\,dt$,故 $(t+2)\,dt = \dfrac{1}{2}\,du$。
$t=0$ 時 $u=1$;$t=2$ 時 $u=13$。
積分變為 $\displaystyle\int_1^{13} \dfrac{1}{2}\sqrt{u}\,du = \dfrac{1}{2}\cdot\dfrac{2}{3}u^{3/2}\Big|_1^{13} = \dfrac{1}{3}(13\sqrt{13} - 1).$
力 $\mathbf{F}$ 沿路徑 $C$ 對物體做的功,只與 $\mathbf{F}$ 在切方向上的投影 $\mathbf{F}\cdot\mathbf{T}$ 有關。一小段弧長 $\Delta s_i$ 的元功 $\approx [\mathbf{F}\cdot\mathbf{T}]\,\Delta s_i$;取極限即為下面的定義。
設 $\mathbf{F}$ 在光滑曲線 $C:\mathbf{r}(t),\ a\le t\le b$ 上連續,定義 $$\int_C \mathbf{F}\cdot d\mathbf{r}=\int_C \mathbf{F}\cdot \mathbf{T}\,ds = \int_a^b \mathbf{F}(\mathbf{r}(t))\cdot \mathbf{r}'(t)\,dt.$$
由 $\mathbf{T} = \mathbf{r}'(t)/\|\mathbf{r}'(t)\|$ 與 $ds = \|\mathbf{r}'(t)\|\,dt$ 即可消去分母;得到的形式直接適合計算。
直接視覺化 例題 6:力場 $\mathbf{F}=-\tfrac{1}{2}x\,\mathbf{i}-\tfrac{1}{2}y\,\mathbf{j}+\tfrac{1}{4}\,\mathbf{k}$(始終指向 $z$ 軸並有一點向上抬)作用在沿螺旋線 $\mathbf{r}(t)=(\cos t,\sin t,t)$ 移動的粒子上,從 $(1,0,0)$ 到 $(-1,0,3\pi)$ 做的功。右側的 $F\cdot r'$ 是當下瞬間功率,累積數值是 $\int_0^t \mathbf{F}(\mathbf{r}(s))\cdot\mathbf{r}'(s)\,ds$;進度條只表示螺旋路徑已播放到哪裡。
$\mathbf{F}(\mathbf{r}(t)) = (-\tfrac{1}{2}\cos t,\,-\tfrac{1}{2}\sin t,\,\tfrac{1}{4})$。
$\mathbf{r}'(t) = (-\sin t,\,\cos t,\,1)$。
$\mathbf{F}\cdot \mathbf{r}' = \tfrac{1}{2}\sin t\cos t - \tfrac{1}{2}\sin t\cos t + \tfrac{1}{4} = \dfrac{1}{4}$。
前兩項剛好抵消(這是因為力的 $xy$ 分量正交於螺旋切向的 $xy$ 部分)。
$W = \displaystyle\int_0^{3\pi} \dfrac{1}{4}\,dt = \dfrac{3\pi}{4}$。
$\mathbf{F}=(y,\,x^2)=(4t-t^2,\,(4-t)^2)$。$\mathbf{r}_1'(t)=(-1,\,4-2t)$。
$\mathbf{F}\cdot\mathbf{r}_1' = -(4t-t^2)+(4-t)^2(4-2t)$。
展開並對 $t\in[0,3]$ 積分得 $\displaystyle\int_{C_1}\mathbf{F}\cdot d\mathbf{r} = \dfrac{69}{2}$。
$\mathbf{F}=(4t-t^2,\,t^2)$,$\mathbf{r}_2'(t)=(1,\,4-2t)$。
$\mathbf{F}\cdot\mathbf{r}_2' = (4t-t^2)+t^2(4-2t)$。
對 $t\in[1,4]$ 積分得 $\displaystyle\int_{C_2}\mathbf{F}\cdot d\mathbf{r} = -\dfrac{69}{2}$。
把 $\mathbf{F}\cdot d\mathbf{r}$ 展開: $$\int_C \mathbf{F}\cdot d\mathbf{r} = \int_C (M\,\mathbf{i}+N\,\mathbf{j})\cdot(dx\,\mathbf{i}+dy\,\mathbf{j}) = \int_C M\,dx + N\,dy.$$ 在三維是 $\int_C M\,dx + N\,dy + P\,dz$。Green 定理 (§15.4) 就使用這種微分形式。
$dx = -3\sin t\,dt$,$dy = 3\cos t\,dt$。
$y^3 = 27\sin^3 t$;$x^3 + 3xy^2 = 27\cos^3 t + 81\cos t\sin^2 t = 27\cos t\,(1+2\sin^2 t)$。
$y^3\,dx = -81\sin^4 t\,dt$。
$(x^3+3xy^2)\,dy = 27\cos t(1+2\sin^2 t)\cdot 3\cos t\,dt = 81\cos^2 t(1+2\sin^2 t)\,dt$。
合計:$81[\cos^2 t + 2\sin^2 t\cos^2 t - \sin^4 t]\,dt$。
在 $[0,2\pi]$ 上:$\int\cos^2 t\,dt = \pi$,$\int 2\sin^2 t\cos^2 t\,dt = \pi/2$,$\int\sin^4 t\,dt = 3\pi/4$。
合計 $\displaystyle 81\!\left[\pi + \dfrac{\pi}{2} - \dfrac{3\pi}{4}\right] = 81\cdot\dfrac{3\pi}{4} = \dfrac{243\pi}{4}$。
在 §15.4 例題 2 用 Green 定理可一步算出同樣的 $\dfrac{243\pi}{4}$,比這裡直接展開簡潔得多。
重力場做的功有個漂亮性質:在物體從 A 移到 B 的過程中,不論走哪條路徑,功都一樣。本節要把這個現象推廣 —— 它等價於「向量場是保守的」,並且也等價於「沿任何閉曲線積分為零」。三件事互相蘊含,合稱線積分基本定理。
$M=\tfrac{1}{2}xy,\ N=\tfrac{1}{4}x^2$,$\partial M/\partial y=\tfrac{1}{2}x=\partial N/\partial x$ ✓ 保守。位能 $f=\tfrac{1}{4}x^2 y + K$(驗 $f_x=\tfrac{1}{2}xy$、$f_y=\tfrac{1}{4}x^2$)。
$\mathbf{F}=(\tfrac{1}{2}t^2,\tfrac{1}{4}t^2)$,$\mathbf{r}'=(1,1)$,$\mathbf{F}\cdot\mathbf{r}'=\tfrac{3}{4}t^2$。
$\int_0^1 \tfrac{3}{4}t^2\,dt = \tfrac{1}{4}.$
$\mathbf{F}=(\tfrac{1}{2}t^3,\tfrac{1}{4}t^4)$,$\mathbf{r}'=(2t,1)$,$\mathbf{F}\cdot\mathbf{r}'=t^4 + \tfrac{1}{4}t^4=\tfrac{5}{4}t^4$。
$\int_0^1 \tfrac{5}{4}t^4\,dt=\tfrac{1}{4}.$
$\mathbf{F}=(\tfrac{1}{2}t\cdot t^3,\tfrac{1}{4}t^2)=(\tfrac{1}{2}t^4,\tfrac{1}{4}t^2)$,$\mathbf{r}'=(1,3t^2)$。
$\mathbf{F}\cdot\mathbf{r}'=\tfrac{1}{2}t^4+\tfrac{3}{4}t^4=\tfrac{5}{4}t^4$。$\int_0^1 \tfrac{5}{4}t^4\,dt=\tfrac{1}{4}.$
設分段光滑曲線 $C:\mathbf{r}(t),\ a\le t\le b$ 位於開區域 $R$ 內。若 $\mathbf{F}=M\,\mathbf{i}+N\,\mathbf{j}$ 在 $R$ 上保守(即 $\mathbf{F}=\nabla f$)且 $M,N$ 連續,則 $$\int_C \mathbf{F}\cdot d\mathbf{r} = \int_C \nabla f\cdot d\mathbf{r} = f(\mathbf{r}(b)) - f(\mathbf{r}(a)).$$
空間版本完全類似:$\int_C \nabla f\cdot d\mathbf{r}=f(B)-f(A)$。直觀:保守場做的功=位能差,跟走哪條路無關。
$\partial M/\partial y = 2x = \partial N/\partial x$ ✓ 保守。由 §15.1 例題 6 已得位能 $f(x,y) = x^2 y - \dfrac{y^2}{2}$。
$f(1,2) = 1\cdot 2 - \dfrac{4}{2} = 2 - 2 = 0$。
$f(-1,4) = (-1)^2 \cdot 4 - \dfrac{16}{2} = 4 - 8 = -4$。
$\displaystyle\int_C \mathbf{F}\cdot d\mathbf{r} = f(1,2) - f(-1,4) = 0 - (-4) = 4$。
若 $\mathbf{F}$ 在開連通區域上連續,則 $\displaystyle\int_C \mathbf{F}\cdot d\mathbf{r}$ 路徑無關 ⇔ $\mathbf{F}$ 保守。
「連通」表示任兩點可由區域內的分段光滑曲線連起來。「路徑無關」表示固定起點與終點時,所有合法路徑的積分值相同。
下面互動圖固定起點 $A=(-1,4)$ 與終點 $B=(1,2)$(也就是例題 2 的端點);提供三條外型不同的路徑:直線、上凸圓弧、下凸圓弧。同時可切換兩種向量場:保守場 $\mathbf{F}=2xy\,\mathbf{i}+(x^2-y)\,\mathbf{j}$(背景顯示位能 $f=x^2 y - y^2/2$ 的等高線)與非保守場 $\mathbf{F}=-y\,\mathbf{i}+x\,\mathbf{j}$(旋轉場)。觀察三條路徑的積分值是否一致。
$M=e^x\cos y,\ N=-e^x\sin y,\ P=2$。
$P_y - N_z = 0 - 0 = 0$;
$M_z - P_x = 0 - 0 = 0$;
$N_x - M_y = -e^x\sin y - (-e^x\sin y) = 0$。
三項皆零 ⇒ $\mathrm{curl}\,\mathbf{F}=\mathbf{0}$ ⇒ 保守。
由 $f_x = e^x \cos y$:$f = e^x \cos y + g(y,z)$。
$f_y = -e^x \sin y + g_y(y,z) = -e^x \sin y$ ⇒ $g_y = 0$ ⇒ $g = h(z)$。
$f_z = h'(z) = 2$ ⇒ $h(z) = 2z + K$。
$\boxed{f(x,y,z) = e^x \cos y + 2z + K}$。
$f(1,\pi,3) = e\cos\pi + 6 = -e + 6$。
$f(0,\pi/2,1) = e^0 \cos(\pi/2) + 2 = 0 + 2 = 2$。
$W = f(1,\pi,3) - f(0,\pi/2,1) = (-e + 6) - 2 = 4 - e$。
設 $\mathbf{F}=M\,\mathbf{i}+N\,\mathbf{j}+P\,\mathbf{k}$ 在開連通區域 $R$ 上具連續一階偏導,$C$ 為 $R$ 內分段光滑曲線。下列三條互相等價:
令 $\mathbf{F}_1 = y^3\,\mathbf{i} + 3xy^2\,\mathbf{j}$ 與 $\mathbf{F}_2 = \mathbf{i} + \mathbf{j}$。則 $\mathbf{F}=\mathbf{F}_1+\mathbf{F}_2$。
$\partial M_1/\partial y = 3y^2 = \partial N_1/\partial x$ ✓。位能 $f_1 = xy^3$。
$\displaystyle\int_{C_1}\mathbf{F}_1\cdot d\mathbf{r} = f_1(2,0)-f_1(0,0) = 0 - 0 = 0$。
$C_1:\mathbf{r}(t)=(1-\cos t,\,\sin t),\ 0\le t\le \pi$(順著 $(0,0)\to(2,0)$ 走上半圓)。
$\mathbf{r}'(t)=(\sin t,\cos t)$;$\mathbf{F}_2\cdot\mathbf{r}'=\sin t + \cos t$。
$\displaystyle\int_0^\pi (\sin t+\cos t)\,dt = [-\cos t+\sin t]_0^\pi = (1+0)-(-1+0)=2$。
Green 定理把邊界線積分和區域雙重積分畫上等號 —— 它說:在一個簡單連通的平面區域 $R$ 上,沿其邊界 $C$ 的線積分(逆時針)等於某個雙重積分。這是把1D 微積分基本定理推廣到平面的版本,也是後面 Stokes 定理(§15.8)與散度定理(§15.7)的二維特例。
簡單曲線:$\mathbf{r}(c)\ne\mathbf{r}(d)$,其中 $c,d\in(a,b)$ — 不自交。
簡單連通區域:任意簡單閉曲線所圍區域都仍在 $R$ 內 — 直觀說「沒有洞」。圓盤是;圓環不是。
設 $R$ 為簡單連通區域,其邊界 $C$ 為分段光滑、逆時針定向(即 $R$ 永遠在 $C$ 的左側)。若 $M,N$ 在含 $R$ 的開區域上具連續一階偏導,則 $$\oint_C M\,dx + N\,dy = \iint_R\!\left(\dfrac{\partial N}{\partial x} - \dfrac{\partial M}{\partial y}\right) dA.$$
當 $\mathbf{F}=M\,\mathbf{i}+N\,\mathbf{j}$ 為保守場,$\partial N/\partial x = \partial M/\partial y$ ⇒ 右邊為 0 ⇒ 與定理 15.7 第三條(閉曲線積分為 0)一致。
$M=y^3,\ N=x^3+3xy^2$。
$\partial M/\partial y = 3y^2$;$\partial N/\partial x = 3x^2 + 3y^2$。
差為 $3x^2 + 3y^2 - 3y^2 = 3x^2$。
$y=x^3$ 與 $y=x$ 在 $0\le x\le 1$ 之間,下方為 $y=x^3$、上方為 $y=x$(因為 $x>x^3$ 當 $0
$\displaystyle\iint_R 3x^2\,dA = \int_0^1 \int_{x^3}^{x} 3x^2 \,dy\,dx = \int_0^1 3x^2(x - x^3)\,dx = \int_0^1 (3x^3 - 3x^5)\,dx = \dfrac{3}{4} - \dfrac{1}{2} = \dfrac{1}{4}$。
$\partial N/\partial x - \partial M/\partial y = 3x^2 + 3y^2 - 3y^2 = 3x^2$。
Green 給:$W = \displaystyle\iint_R 3x^2\,dA$,其中 $R$ 為以原點為中心、半徑 3 的圓盤。
$x = r\cos\theta$,$dA = r\,dr\,d\theta$:
$\displaystyle W = \int_0^{2\pi}\!\int_0^3 3r^2\cos^2\theta\cdot r\,dr\,d\theta = 3\!\left(\int_0^{2\pi}\cos^2\theta\,d\theta\right)\!\left(\int_0^3 r^3\,dr\right).$
$\int_0^{2\pi}\cos^2\theta\,d\theta=\pi$;$\int_0^3 r^3\,dr=\dfrac{81}{4}$。
$W = 3 \cdot \pi \cdot \dfrac{81}{4} = \dfrac{243\pi}{4}$。
$M = y^3$,$N = 3xy^2$。
$\partial M/\partial y = 3y^2 = \partial N/\partial x$ ✓ 保守。位能 $f = xy^3$(驗 $f_x = y^3$、$f_y = 3xy^2$)。
因為 $C$ 是閉曲線,且 $\mathbf{F}$ 在含 $R$ 的開區域上保守,由定理 15.7 第三條:
$\displaystyle\oint_C \mathbf{F}\cdot d\mathbf{r} = 0$。
或用 Green:$\partial N/\partial x - \partial M/\partial y = 0$,雙重積分自然為 0。
左右兩個圖同步播放動畫:左邊 marker 沿邊界 $C$ 走,累積 $\oint M\,dx + N\,dy$;右邊用固定中點格點在區域 $R$ 上累積 $\iint(\partial N/\partial x-\partial M/\partial y)\,dA$。動畫點的顯示會降採樣以保持流暢,但右側數值使用完整高解析度格點與 prefix sum 計算,因此 reset 後會得到相同掃描順序與相同終點誤差。進度條只表示動畫時間進度;Green 定理的驗證看 banner 中的左右數值差。
把 Green 定理「反過來用」:選擇 $M, N$ 使 $\partial N/\partial x - \partial M/\partial y = 1$,則 $\oint_C M\,dx + N\,dy = \iint_R dA = $ 區域面積。最常用的選擇是 $M = -y/2,\ N = x/2$。
若 $R$ 是由分段光滑簡單閉曲線 $C$(逆時針)圍成的區域,則 $$\text{面積}(R) = \dfrac{1}{2}\oint_C (x\,dy - y\,dx).$$
$\mathbf{r}(t) = a\cos t\,\mathbf{i} + b\sin t\,\mathbf{j},\ 0\le t\le 2\pi$ — 逆時針一週。
$dx = -a\sin t\,dt$,$dy = b\cos t\,dt$。
$x\,dy - y\,dx = (a\cos t)(b\cos t\,dt) - (b\sin t)(-a\sin t\,dt) = ab(\cos^2 t + \sin^2 t)\,dt = ab\,dt$。
$\text{面積} = \dfrac{1}{2}\int_0^{2\pi} ab\,dt = \dfrac{1}{2}\cdot 2\pi ab = \pi ab$。
設 $\mathbf{F}=M\,\mathbf{i}+N\,\mathbf{j}$,把 $\mathbf{F}$ 推廣到 3D($P=0$),則:
這兩個寫法分別會在 §15.8 Stokes 定理(推到 3D 曲面)與 §15.7 散度定理(推到 3D 區域)變得更一般。
把本章前四節要記住的公式、定理、運算子整理成一頁。第一次看可能不熟,多回來查幾次就熟了。
純量函數 → 向量場。在每點指向 $f$ 上升最快的方向。
向量場 → 向量場。curl F = 0 ⇔ F 保守(在開球上,§15.1 定理 15.2)。
向量場 → 純量。div F = 0 ⇔ F 不可壓縮 / 無散度。恆等式:div(curl F) = 0(定理 15.3)。
| 對象 | 輸入 | 輸出 |
|---|---|---|
| 梯度 ∇f | 純量場 f | 向量場 |
| 旋度 ∇×F | 向量場 F | 向量場 |
| 散度 ∇·F | 向量場 F | 純量場 |
對方向無關。$f\equiv 1$ 時給弧長。$\|\mathbf{r}'\| = \sqrt{x'^2 + y'^2 + z'^2}$。
方向會影響符號:$\int_{-C} = -\int_C$。物理意義:力 F 沿 C 做的功。
F 保守 ⇒ 線積分只取決於兩端點。
前提:F 在開連通區域上具連續一階偏導。
區域必須開且簡單連通(否則只有「F 保守 ⇒ curl = 0」一方向)。
前提:R 簡單連通,C 分段光滑、逆時針(R 在左側),M、N 一階偏導連續。
取 M = −y/2、N = x/2,∂N/∂x − ∂M/∂y = 1。橢圓 (a,b):$\pi ab$。
旋度形式是 Stokes 定理(§15.8)在平面上的特例;散度形式是散度定理(§15.7)的二維版本。
| 主題 | 關鍵物件 | 核心關係 |
|---|---|---|
| 15.1 向量場 | F、∇f、curl F、div F | F = ∇f ⇔ F 保守 |
| 15.2 線積分 | ∫f ds、∫F·dr | 把曲線參數化轉成 1D 定積分 |
| 15.3 路徑無關 | 位能 f(B) − f(A) | 保守 ⇔ 路徑無關 ⇔ ∮ = 0 |
| 15.4 Green | 邊界 C、區域 R | ∮_C = ∬_R(N_x − M_y) |