39203 圖解梯度、散度與旋度

終極密碼

遊戲規則:本遊戲為猜密碼的遊戲。密碼為0到100之間的其中1個整數,電腦會提示密碼的所在範圍,玩家必須在6次之內猜到密碼才能過關。
★ 終極密碼為0到100之間 ★
您共有六次機會

『如果 $\sigma$ 代表向量函數, $\nabla \sigma$ 就包含有向量與純量兩部分, 它們可以寫成 $S\nabla \sigma$ 和 $V\nabla \sigma$, 我建議純量部分為 $\sigma$ 的散度 $($divergence$)$, 這是向量函數效應的好名字。 但是 $\nabla \sigma$, 一般說來也是有向量的部分, 可以稱這種向量為原向量函數的旋度 $($curl $)$ 或變體 $($version$)$1 1 version 是翻譯、譯本、敘述或異見, 但這個字應分解為versio-n-其拉丁文是verterevers- 的意思是轉動(to turn)。所以它也翻為a special form or variant of something $($變體$)$ , 它代表了由向量$\sigma$ 所承載之內容旋轉的量與方向。』
                                        --- James Clerk Maxwell (1831$\sim$1879) ---

一、前言

正如微積分是為了解決牛頓力學而發明的, 向量分析則是應運電磁學而產生的。 學向量分析心中常有電磁學讓你永遠不孤單也不會迷失。我願意再強調一次

『數學沒有物理是瞎子, 物理沒有數學是跛子。』

一個數學系特別是應用數學系把『普通物理』列為選修(等價於不用修), 等於告訴世人它是一個爛系。歷史告訴我們的還有時代的潮流就是跨領域! 一個只懂數學(本質上不是真的懂), 而對其它領域無知的人是絕對沒有競爭力。 就算是數學專業的人也必須同時對於非主修的科目有所涉獵, 如此才可以培養寬廣的胸襟與開放的心靈, 對於別人的研究也有興趣並懂得欣賞。否則會越來越孤立以至於枯竭, 片面內省的數學玄想只會導致貧瘠。 百年前德國哥廷根偉大數學傳統的建立者 Felix Klein (1849$\sim$1925) 就已經有這樣的呼籲:

『我們強烈感受到, 在現代科學思想的快速發展下, 我們的學科正瀕臨孤立的危險。 自從現代分析興起以來, 讓數學與理論自然科學雙方皆獲益的親密互動關係, 已經面臨崩潰的威脅。 \cdots 今日數學界的當務之急, 是將其純科學部門和其具有最重要應用價值的物理科學部門再次緊密結合, 這種關係以往所得的豐碩成果可以從 Lagrange 和 Gauss 的研究裡看到。』

他的吶喊現在看來仍然是先知的警語, 令人深思。 看來數學的發展與文明或民主一樣, 不見得都是微分大於0(增函數), 永遠是進步的! 如果你不守護它, 深化其內涵並隨時進步, 一旦獨裁的敗家子出現, 整個數學社會被迫走入歧途, 則幾代人的努力馬上就煙消雲散。

1865年左右, 馬克士威(James Clerk Maxwell, 1831$\sim$1879)綜合電磁現象的安培定律及法拉第定律, 建立電磁場的數學理論, 他提出了變化的磁場激發電場的概念之後, 發現還缺少某些重要的東西。 假如對這些方程式進行數學變換, 那麼就得出一種與 電磁學方程組體系的其它公式矛盾的關係。假如這個方程式藉著一個附加項, 即所謂的位移電流來補充, 那麼這個矛盾就可以克服。 所以他提出 位移電流來克服電磁學方程組之矛盾, 而且整理成完整的電磁場理論, 因此馬克士威就從法拉第的解釋者, 成為獨創的學者。 馬克士威方程包含了完整的電學、磁學和光學的理論, 這是馬克士威對物理學所作的 一項偉大的統一, 將電、磁、光三者合而為一。但是馬克士威的成果在當時(1865)並沒有立即引起數學界的注意與興趣, 所以並不像牛頓與 Leibniz (1646$\sim$1716)那個年代承襲古希臘的信念:

『大自然是根據數學設計的。』

他們將數學看成是通往大自然的真理之路。 這個信仰吸引了 Bernoulli 家族、 Euler (1707 $\sim$1783)、 Lagrange (1736$\sim$1813)、 Laplace (1749$\sim$1827)、 $\cdots\cdots$、 一直到 Fourier (1768 $\sim$1830), 這些同時具有科學家身份的大數學家為之獻身, 進而開創了 17$\sim$18世紀數學創造最偉大的年代。 更重要的是除了科學之外他們也將理性的精神注入到人文、社會領域。 如果 19 世紀後半葉的數學家們能夠像 Euler 等人一樣將牛頓的工作視為研究的核心, 也將 Maxwell 方程當為最根本重要的問題來關心, 則他們極有可能早就發現愛因斯坦的狹義相對論 (special relativity)、 拓撲群及其表現理論還有雙曲偏微分程的大部分理論以及泛函分析。 但是數學界卻失去了這大好的機會與 Maxwell 方程攜手開創數學與科學史的第二個黃金年代。 事實上那時候數學界並不缺乏有見識的人, 例如, 牛津的數學家 Henry Smith (1826$\sim$1883) 就已經體會到需要新的數學以了解 Maxwell 方程, 並在 1873 年左右呼籲數學家研究馬克士威的著作。 但 Smith 已經 46 歲而且他的學養不足以成為新領域的開創者, 有能力的年輕人並不聽從他的建議。 在 20 世紀完成 Smith 先知呼聲的 Hermann Minkowski (1864$\sim$1909)與 Jacoques Hadamard (1865$\sim$1963)那時分別是 9 歲與 8 歲, Elie Cartan (1869$\sim$1951) 只有 3 歲, 而 Hermann Weyl (1885$\sim$1955), Jean Leray (1906$\sim$1998) 與 Harish Chandra (1923$\sim$1983)則還沒有出生。 我想最主要的原因是19世紀後半葉數學家的興趣是複變函數論、 解析數論與不變量理論, 而這些豐富的成果已經形成數學界美學與價值觀的標準, 以至於形式上複雜的 Maxwell 方程被拒絕了。

雖然我們可以說馬克士威是向量分析的創始人, 但他那個時候還沒有向量的符號, 以至於方程式必須以各個分量寫出所以非常的繁複, 這使得他偉大的創見並沒有馬上為當時的人所認可。 如果沒有專門且精準的符號, 則現代數學根本不可能產生。 向量分析最重要的進展主要是歸功於 J. W. Gibbs (1839$\sim$1903)與 O. Heaviside (1850$\sim$1925), 甚至連符號也是他們兩人的貢獻。 其中 Gibbs 出版了第一本向量分析課本(Elements of Vector Analysis), 向量分析這個專有名詞正式在此出現。 實際上基於非原創性的理由 Gibbs 並不情願出版這本教材, 後來是由他的學生 E. B. Wilson 整理之後才成為正式的教科書。 Gibbs 與 Heaviside 向量符號之比較如下:

內積 外積 梯度 散度 旋度
Gibbs $\bf a\cdot\bf b$ $\bf a\times\bf b$ $\nabla$ $\nabla\cdot$ $\nabla\times$
Heaviside $\bf a\bf b$ $V\bf a\bf b$ $\nabla$ div curl

梯度、散度與旋度形式上來看都是一次微分, 但本質上卻有極大的差異。要瞭解它們就必須從 Gauss 散度定理與 Stokes 定理著手, 雖然這兩個定理的內涵都是微積分基本定理, 但是一個關心的是法向量 (散度)另一個則是切向量(旋度)。 在二維平面給定切向量 $(a,b)$ 則朝外法向量為 $(b,-a)$, 反過來亦然, 知道其中一個向量立即可得另一個向量。 所以二維的 Green 定理同時具有 Gauss 散度定理與 Stokes 定理兩個身分。 讀者如果會怕三維空間, 那麼就先把 Green 定理弄清楚, 但不能只學一般微積分課本教你的切向量之情形。 這就是主動與被動學習的差異, 你(妳)是要當自己生命的主人還是等著讓別人來決定你(妳)的未來! 這兩個積分定理的主要貢獻者為

(1) G. Green (1793$\sim$1841) : 1828年提出 Green 定理, Green 等式。

(2) 高斯 (Karl Friedrich Gauss, 1777$\sim$1855) : 提出散度定理。

(3) M. Ostogradsky (1801$\sim$1861) : 1831年研究熱傳導方程(heat equation) 時也導出散度定理。

(4) George Gabriel Stokes (1819$\sim$1903) : 1850 年給 Lord Kelvin (1824$\sim$1907)的信中提到 Stokes定理, Stokes 1854 年在劍橋的數學競試以這定理為題目, 其中得獎的學生有 J. C. Maxwell。 由此可見學生 Maxwell 遠比老師 Stokes 厲害的多。

本文主要的目的是想要以直觀、 圖解的角度來介紹曲線座標系統的梯度、 散度與旋度這三個向量分析最重要的概念。 大部分數學系只活在一維空間, 把向量分析視為是微積分的推廣。 這是極大的誤解, 同時也帶來錯誤的數學觀。 通常學生對於非垂直座標的梯度、散度與旋度都非常之陌生, 頂多就是把它當成公式, 需要的時候再查表。對此我們並不滿意, 公式並不僅僅是公式, 它會講話, 因此我們特別介紹量綱(因次)分析 (dimensional analysis) 這個簡單 但重要的觀念、直觀且有感覺地認識所有的方程式與(不)等式。 讀者也應嘗試能從公式中解釋各項之意義, 並從中體會公式(或方程式) 本身之物理或幾何意義。 想更深入瞭解梯度、散度與旋度, 我們需要微分形式 (differential form) 這個微分幾何的基本知識, 我們將在下一篇文章介紹。 關於量綱(或因次)分析, 請讀者參閱我在數學傳播的文章。 這個觀念是我最有心得也最引以為豪的方法, 它幫助我直觀地了解方程式, 讓我所學的數學變得有血有淚, 不再是冰冷無情的怪咖。

二、曲線座標系統

回顧一下直角座標的散度 (divergence)、 旋度 (curl) 與梯度 (gradient), 這三個算子 正好可以對應到向量的三個重要運算 : 內積 (inner product)、 外積 (cross product) 與 直積 (direct product) (已知 ${\bf F} = (F_1, F_2, F_3)$ 是向量函數、 $f$ 是純量函數) \begin{equation} \begin{array}{ccl} \displaystyle {內積}:\quad~&{向量} \mapsto {純量}\quad~& {\rm div}\, {\bf F} =\nabla\cdot {\bf F}= \dfrac{\partial F_1}{\partial x} + \dfrac{\partial F_2}{\partial y} + \dfrac{\partial F_3}{\partial z}\\[3mm] \displaystyle {外積}:\quad~ &{向量} \mapsto {向量}\quad~& {\rm curl}\, {\bf F} =\nabla\times {\bf F} =\left|\begin{array}{ccc} {\bf i} & {\bf j} & {\bf k}\\ ~\frac{\partial }{\partial x}~ & ~\frac{\partial }{\partial y}~ & ~\frac{\partial }{\partial z}~\\ F_1 & F_2 & F_3 \end{array}\right| \\ \displaystyle {直積}:\quad~&{純量} \mapsto{向量}\quad~& {\rm grad}f=\nabla f= \bigg(\dfrac{\partial f}{\partial x},\dfrac{\partial f}{\partial y}, \dfrac{\partial f}{\partial z}\bigg) \end{array} \label{2.1} \tag{2.1} \end{equation} 直角座標之弧長元素與體積元素 \begin{equation} ds^2 =d{\bf r}\cdot d {\bf r}= dx^2 + dy^2 +dz^2,\qquad dV = dx dy dz \label{2.2} \tag{2.2} \end{equation} 還有 ${\bf i}$, ${\bf j}$, ${\bf k}$ 三個方向之面積元素為 \begin{equation} \begin{array}{c} dS_x=dydz,\qquad {x=c\, 之平面 } \\[3pt] dS_y=dzdx,\qquad {y=c\, 之平面} \\[3pt] dS_z=dxdy,\qquad {z=c\, 之平面} \end{array} \label{2.3} \tag{2.3} \end{equation} 在彎曲空間梯度、散度與旋度, 這三個算子與弧長元素、體積元素與面積元素要如何重新定義呢? 這是我們所關心的問題。

空間上任一位置向量(position vector)表示為 \begin{equation} {\bf r}=(x,y,z)=x{\bf{i}}+y{\bf{j}}+z{\bf{k}} \label{2.4} \tag{2.4} \end{equation} 其中 ${\bf{i}} =(1,0,0)$, ${\bf{j}} =(0,1,0)$, ${\bf{k}} =(0,0,1)$ 就是直角座標的標準基底, 而且滿足右手法則 (${\bf i}\to{\bf j}\to{\bf k}$) $$ {\bf i}\cdot {\bf i}={\bf j}\cdot {\bf j}={\bf k}\cdot {\bf k}=1 $$ $$ {\bf i}\times {\bf j}={\bf k},\qquad {\bf j}\times {\bf k}={\bf i},\qquad {\bf k}\times {\bf i}={\bf j} $$ $$ {\bf i}\times {\bf j}=-{\bf j}\times {\bf i},\qquad {\bf j}\times {\bf k}=-{\bf k}\times {\bf j},\qquad {\bf k}\times {\bf i}=-{\bf i}\times {\bf k} $$ 如果在線性空間上, 我們可以透過線性變換(即矩陣)轉換成其它基底。 但是在彎曲空間, 這種想法已無法適用, 那麼座標變換要如何重新定義呢? 仍然回到\eqref{2.4}, 我們觀察到如果把 $x,y,z$ 視為變數, 那麼對位置向量 ${\bf r}$ 取偏導數, 座標系統可表示為 \begin{equation} \{{\bf{i}},{\bf{j}},{\bf{k}}\} =\{{\bf e}_x,{\bf e}_y,{\bf e}_z\} =\Big\{\frac{\partial {\bf r}}{\partial x}, \frac{\partial {\bf r}}{\partial y}, \frac{\partial {\bf r}}{\partial z}\Big\} \label{2.5} \tag{2.5} \end{equation} 這個事實由全微分來看是自然的 \begin{equation} d{\bf r} =\frac{\partial {\bf r}}{\partial x}dx+ \frac{\partial {\bf r}}{\partial y}dy+ \frac{\partial {\bf r}}{\partial z}dz = dx {\bf i} + dy {\bf j} +dz {\bf k} \label{2.6} \tag{2.6} \end{equation} 換句話說, $\{{\bf i},{\bf j},{\bf k}\}$ 或 $\{{\bf e}_x,{\bf e}_y,{\bf e}_z\}$ 可由位置向量${\bf r}$之一階偏導數來表示。 全微分並不受座標之限制, \eqref{2.6} 可以換成任意的座標, 也因此提供我們推廣至彎曲空間的思考方式。

考慮一般的曲線座標系統 (curvilinear system)。 已知 $x_1,x_2,x_3$ 是 $u_1,u_2,u_3$ 的函數 \begin{equation} x_i = x_i(u_1,u_2,u_3),\quad u_i=u_i(x_1,x_2,x_3),\quad i=1,2,3 \label{2.7} \tag{2.7} \end{equation} 我們要求 $(x_1,x_2,x_3)\longmapsto (u_1,u_2,u_3)$ 這個變換是一對一, 即 Jacobian \begin{equation} \frac{\partial (x_1,x_2,x_3)}{\partial (u_1,u_2,u_3)} =\left|\begin{array}{ccc} ~\frac{\partial x_1}{\partial u_1}~ & ~\frac{\partial x_1}{\partial u_2}~ & ~\frac{\partial x_1}{\partial u_3}~\cr \frac{\partial x_2}{\partial u_1} & \frac{\partial x_2}{\partial u_2} & \frac{\partial x_2}{\partial u_3}\cr \frac{\partial x_3}{\partial u_1} & \frac{\partial x_3}{\partial u_2} & \frac{\partial x_3}{\partial u_3}\cr \end{array}\right| =\frac{\partial {\bf r}}{\partial u_1}\cdot \Big(\frac{\partial {\bf r}}{\partial u_2}\times \frac{\partial {\bf r}}{\partial u_3}\Big)\not= 0 \label{2.8} \tag{2.8} \end{equation} 向量 $\{\frac{\partial {\bf r}}{\partial u_1}, \frac{\partial {\bf r}}{\partial u_2}, \frac{\partial {\bf r}}{\partial u_3}\}$ 正是 $\{u_1, u_2,u_3\}$ 三個曲線之切向量。 條件 \eqref{2.8} 保證這組向量是線性獨立 (linear independent), 仿 \eqref{2.5} 我們可以選取這組來做為基底, 再化為單位向量(這是學數學的好習慣) \begin{eqnarray} {\bf e}_1&=&\frac1{h_1}\frac{\partial {\bf r}}{\partial u_1},\quad {\bf e}_2=\frac1{h_2}\frac{\partial {\bf r}}{\partial u_2},\quad {\bf e}_3=\frac1{h_3}\frac{\partial {\bf r}}{\partial u_3} \label{2.9} \tag{2.9}\\ h_1&=&\Big|\frac{\partial {\bf r}}{\partial u_1}\Big|,\quad\hskip6pt h_2=\Big|\frac{\partial {\bf r}}{\partial u_2}\Big|,\quad\hskip8pt h_3=\Big|\frac{\partial {\bf r}}{\partial u_3}\Big| \label{2.10} \tag{2.10} \end{eqnarray} 因此 $\{{\bf e}_1,{\bf e}_2,{\bf e}_3\}$ 是 $u_1,u_2,u_3$-座標曲線之單位向量, 就是我們所要的彎曲空間之基底, 同時也是 $u_1, u_2, u_3$-座標曲線的單位切向量。 若 $\{ {\bf e}_1,{\bf e}_2,{\bf e}_3 \}$ 互相垂直則稱曲線座標系為 正交曲線座標系。因為我們的目的是詮釋梯度、散度與旋度的物理意義, 為避免節外生枝, 在這篇文章我們只討論正交曲線座標系。

圖1.正交曲線座標系統

由全微分得 \begin{equation} d{\bf r} = {\partial {\bf r} \over\partial u_1} du_1 + {\partial {\bf r} \over\partial u_2} du_2 + {\partial {\bf r} \over\partial u_3} du_3 = {h_1 du_1} {\bf e}_1 + {h_2 du_2} {\bf e}_2 + {h_3 du_3} {\bf e}_3 \label{2.11} \tag{2.11} \end{equation} 值得一提的是$du_i$永遠要配合$h_i$, 因為$h_i du_i$才是長度, 其量綱為$[h_i du_i]={\rm L}$, 這個事實我們一直會用到。由畢氏定理容易看出來弧長元素 \begin{equation} ds^2 = d{\bf r}\cdot d{\bf r} = h^2_1 du^2_1 + h^2_2 du^2_2 + h^2_3 du^2_3 \label{2.12} \tag{2.12} \end{equation} 體積元素與各個面的表面積則為(請參考圖1) \begin{eqnarray} dV&=&dx_1 dx_2 dx_3 = {\partial (x_1, x_2, x_3)\over \partial (u_1, u_2, u_3)} du_1 du_2 du_3 = {h_1 h_2 h_3} du_1 du_2 du_3 \label{2.13} \tag{2.13}\\ dS_1&=&h_2h_3du_2du_3,\qquad dS_2=h_3h_1du_3du_1,\qquad dS_3=h_1h_2du_1du_2 \label{2.14} \tag{2.14} \end{eqnarray} 讀者很容易可以驗證 $$ [dV]={\rm L}^3, \qquad [dS_i]={\rm L}^2, \quad i=1,2,3 $$

例題 2.1 {圓柱座標} 空間上任意一點 $(x,y,z)$ 可以表示為圓柱座標 $(r,\theta,z)$ \begin{equation} {\bf r}=(x,y,z)=(r\cos\theta,r\sin\theta,z) \label{2.15} \tag{2.15} \end{equation} \begin{equation} h_1=h_r=\bigg|\frac{\partial {\bf r}}{\partial \rho}\bigg|=1,\qquad h_2=h_\theta=\bigg|\frac{\partial {\bf r}}{\partial \theta}\bigg|=r,\qquad h_3=h_z=\bigg|\frac{\partial {\bf r}}{\partial z}\bigg|=1 \label{2.16} \tag{2.16} \end{equation} 所以 \begin{eqnarray} \displaystyle{\bf e}_r &=&\frac1{h_1}\frac{\partial {\bf r}}{\partial r} =(\cos\theta,\sin\theta,0)=\cos\theta{\bf i}+\sin\theta{\bf j}\nonumber\\ \displaystyle{\bf e}_\theta &=&\frac1{h_2}\frac{\partial {\bf r}}{\partial\theta} =(-\sin\theta,\cos\theta,0)=-\sin\theta{\bf i}+\cos\theta{\bf j}\label{2.17} \tag{2.17}\\ \displaystyle{\bf e}_z &=&\frac1{h_3}\frac{\partial {\bf r}}{\partial z}=(0,0,1) ={\bf k}\nonumber \end{eqnarray} 容易驗證 $\{{\bf e}_r,{\bf e}_\theta,{\bf e}_z\}$ 是彼此互相垂直, 而且滿足右手法則 \begin{equation} \begin{array}{ccc} \displaystyle{\bf e}_r\perp {\bf e}_\theta,& {\bf e}_\theta\perp {\bf e}_z,& {\bf e}_z\perp {\bf e}_r\\[3mm] \displaystyle{\bf e}_r\times {\bf e}_\theta={\bf e}_z,\quad~ & {\bf e}_\theta\times {\bf e}_z={\bf e}_r,\quad~& {\bf e}_z\times {\bf e}_r={\bf e}_\theta \end{array} \label{2.18} \tag{2.18} \end{equation} 為何其右手法則之順序是 $r\to \theta\to z$ 呢$?$ 可以這麼看 $:$ 從原點沿著半徑到 $r$, 而後逆時針轉 $\theta$, 最後再垂直爬升到 $z$ (參考圖2 圖形路徑 $1\to 2\to 3$), 是不是正好就是右手法則 $!$

圖2.圓柱座標

位置向量之全微分為 \begin{equation} d{\bf r}= dr {\bf e}_r + rd\theta {\bf e}_\theta + dz {\bf e}_z, \qquad [dr]=[rd\theta]=[dz]={\rm L} \label{2.19} \tag{2.19} \end{equation} 而圓柱座標之弧長元素、體積元素與各個面的表面積則為(請參考圖3) \begin{equation} \begin{array}{c} \displaystyle ds^2 = dr^2 + r^2 d\theta^2 + dz^2, \qquad dV = r dr d\theta dz \\[3mm] dS_r = r d\theta dz,\qquad dS_\theta= dr dz,\qquad dS_z= rdr d\theta \end{array} \label{2.20} \tag{2.20} \end{equation} 這裡沒有必要透過全微分與複雜的連鎖律來得到弧長元素 $ds^2$, 因為 $\{{\bf e}_r,{\bf e}_\theta,{\bf e}_z\}$ 是一正交系統, 無窮小方塊(蛋糕形狀)的長寬高為 $dr, rd\theta, dz$, 所以由畢氏定理容易看出來弧長元素 \eqref{2.20}$_1$, 一旦長寬高決定了則體積元素與各個面的表面積可由圖形直觀地看出來。

圖3.圓柱座標各個面

例題 2.2 {球座標} 空間上任意一點 $(x,y,z)$ 可以表示為球座標 $(\rho,\varphi,\theta)$ \begin{equation} {\bf r} =(x,y,z)=(\rho\sin\varphi\cos\theta, \rho\sin\varphi\sin\theta,\rho\cos\varphi) \label{2.21} \tag{2.21} \end{equation} \begin{equation} h_1=h_\rho =\bigg|\frac{\partial {\bf r}}{\partial \rho}\bigg|=1,\qquad h_2=h_\varphi=\bigg|\frac{\partial {\bf r}}{\partial \varphi}\bigg|=\rho,\qquad h_3=h_\theta =\bigg|\frac{\partial {\bf r}}{\partial \theta}\bigg|=\rho \sin\varphi \label{2.22} \tag{2.22} \end{equation} 因此 \begin{equation} \begin{array}{l} \displaystyle {\bf e}_\rho =\frac1{h_1}\frac{\partial {\bf r}}{\partial \rho} =(\sin\varphi\cos\theta,\sin\varphi\sin\theta,\cos\varphi) \\[12pt] \displaystyle {\bf e}_\varphi =\frac1{h_2}\frac{\partial {\bf r}}{\partial\varphi} =(\cos\varphi\cos\theta,\cos\varphi\sin\theta,-\sin\varphi) \\[12pt] \displaystyle {\bf e}_\theta =\frac1{h_3}\frac{\partial {\bf r}}{\partial\theta} =(-\sin\theta,\cos\theta,0) \end{array} \label{2.23} \tag{2.23} \end{equation} 實際上$\{{\bf e}_\rho,{\bf e}_\varphi,{\bf e}_\theta\}$ 形成一右手正交座標系統。 為何其右手法則之順序是 $\rho\to \varphi\to\theta$ 呢 $?$ 可以這麼看 $:$ 從原點順著半徑到 $\rho$, 而後垂直往下轉高低角 $\varphi$ $($緯度$!)$, 最後再水平逆時針轉水平角 $\theta$ $($經度$!)$, $($請參考圖4 圖形路徑 $1\to 2\to 3)$, 是不是正好就是右手法則! 位置向量之全微分為

圖4.球座標之右手法則
\begin{equation} d{\bf r} = d\rho {\bf e}_\rho + \rho d\varphi {\bf e}_\varphi + \rho\sin\varphi d\theta {\bf e}_\theta \label{2.24} \tag{2.24} \end{equation} 因為角度與三角函數不具有量綱, $[d\varphi]=[\sin \varphi]=[d\theta]=1$, 所以 $$ [d\rho]=[\rho d\varphi]=[\rho \sin \varphi d\theta]={\rm L} $$ 要特別提醒讀者的是第三個邊長是 $\rho\sin\varphi d\theta$ 並不是 $\rho d\theta$, 因為球面上在相同高低角$($相同緯度$)$兩點之弧長投影至赤道這個平面時其對應的圓的半徑等於 $\rho\sin\varphi$, 所以這兩點的距離為
$\rho\sin\varphi d\theta$, $\sin\varphi$ 這個量顯示的是投影的效應。對於球座標我們仍然先掌握無窮小方塊$($像西瓜切片但要切除尖的部分之形狀$)$ 的三邊長等於 $d\rho$, $\rho d\varphi$, $\rho\sin\varphi d\theta$, 所以球座標之弧長元素、 體積元素與各個面的表面積為 \begin{equation} \begin{array}{c} \displaystyle ds^2 = d\rho^2 + \rho^2 d\varphi^2 + \rho^2 \sin^2 \varphi d\theta^2, \qquad dV = \rho^2\sin \varphi d\rho d\varphi d\theta \\[10pt] \displaystyle dS_\rho = \rho^2\sin \varphi d\varphi d\theta,\qquad dS_\varphi= \rho\sin \varphi d\rho d\theta,\qquad dS_\theta = \rho d\rho d\varphi \end{array} \label{2.25} \tag{2.25} \end{equation}
圖5.球座標

三、方向導數 --- 梯度

正如本節之標題: 梯度本質上是方向導數。 給定方向 $\nu$ 則純量函數 $f$ 之方向導數定義為 \begin{equation} {\partial f\over \partial \nu} =\nabla f \cdot \nu \label{3.1} \tag{3.1} \end{equation} 雖然這是一個等式但左右兩邊卻有本質上的差異。 因為方向導數 ${\partial f\over \partial \nu}$ 與座標之選取無關, 所以要將梯度推廣到曲線座標系統必須從方向導數著手。

圖6.$(u_1, u_2, u_3)$ 到 $(u_1 +\triangle u_1, u_2 +\triangle u_2,u_3 +\triangle u_3)$ 之變化

由 \eqref{3.1} 直觀而言, 梯度 $\nabla f$ 就是方向導數 ${\partial f\over \partial \nu}$ 除以 $\nu$! 考慮 $f$ 從 $(u_1, u_2, u_3)$ 到
$(u_1 +\triangle u_1, u_2 +\triangle u_2,u_3 +\triangle u_3)$ 之變化並利用連接這兩點之弦的向量 $ \triangle {\bf r} = {h_1 \triangle u_1} {\bf e}_1 + {h_2 \triangle u_2} {\bf e}_2 + {h_3 \triangle u_3} {\bf e}_3 $ 得(請參考圖 6) \begin{eqnarray*} \triangle f &=& f(u_1 +\triangle u_1, u_2 +\triangle u_2,u_3 +\triangle u_3) -f(u_1, u_2, u_3)\\ &=& {\partial f \over \partial u_1} \triangle u_1 +{\partial f \over \partial u_2} \triangle u_2 +{\partial f \over \partial u_3} \triangle u_3 +\cdots\\ &=& {1\over {h_1}}{\partial f \over \partial u_1} {h_1}\triangle u_1 +{1\over {h_2}}{\partial f \over \partial u_2} {h_2}\triangle u_2 +{1\over {h_3}}{\partial f \over \partial u_3} {h_3}\triangle u_3 +\cdots\\ &=& \left({1\over h_1}{\partial f \over \partial u_1}{\bf e}_1 +{1\over h_2}{\partial f \over \partial u_2}{\bf e}_2 +{1\over h_3}{\partial f \over \partial u_3}{\bf e}_3\right) \cdot \triangle {\bf r} +\cdots \end{eqnarray*} 再除以長度 $\triangle s=|\triangle {\bf r}|=\sqrt{h_1^2 (\triangle u_1)^2 + h_2^2 (\triangle u_2)^2+ h_3^2 (\triangle u_3)^2}$ $$ {\triangle f \over \triangle s} =\left({1\over h_1}{\partial f \over \partial u_1}{\bf e}_1 +{1\over h_2}{\partial f \over \partial u_2}{\bf e}_2 +{1\over h_3}{\partial f \over \partial u_3}{\bf e}_3\right) \cdot {\triangle {\bf r}\over \triangle s} +\cdots $$ 取極限 $(\triangle s\to 0$) 將弦變為單位切向量 $\nu={d{\bf r} \over ds}$ 得方向導數 \begin{equation} \frac{\partial f}{\partial \nu}= \left({1\over h_1}{\partial f \over \partial u_1}{\bf e}_1 +{1\over h_2}{\partial f \over \partial u_2}{\bf e}_2 +{1\over h_3}{\partial f \over \partial u_3}{\bf e}_3\right)\cdot \nu \label{3.2} \tag{3.2} \end{equation} 比較 \eqref{3.1}、 \eqref{3.2} 兩式得 \begin{equation} \nabla f= {1\over h_1}{\partial f \over \partial u_1}{\bf e}_1 +{1\over h_2}{\partial f \over \partial u_2}{\bf e}_2 +{1\over h_3}{\partial f \over \partial u_3}{\bf e}_3 \label{3.3} \tag{3.3} \end{equation} 藉由類比, 實際上 \eqref{3.3} 這個公式是很容易猜得出來。 只需要將 $\{\bf i, {\bf j}, {\bf k}\}$ 換為 $\{{\bf e}_1, {\bf e}_2, {\bf e}_3\}$、 $\{x, y, z\}=\{x_1, x_2,x_3\}$ 換為 $\{u_1, u_2, u_3\}$, 唯一要注意的是 $\partial u_i$ 一定要伴隨著 $h_i$, 所以每一項都有 ${1\over h_i}$ 是自然的。 而背後真正的理由則是: 量綱平衡。

例題 3.1 圓柱與球座標之梯度 $($請參考圖 $7, 8)$ \begin{eqnarray*} \nabla f &=& {\partial f \over \partial x} \bf i +{\partial f \over \partial y} {\bf j} +{\partial f \over \partial z} {\bf k} \\ &=& {\partial f \over \partial r} {\bf e}_r +{1\over r}{\partial f \over \partial \theta}{\bf e}_\theta +{\partial f \over \partial z} {\bf e}_z \\ &=& {\partial f \over \partial \rho} {\bf e}_\rho +{1\over \rho}{\partial f \over \partial \varphi} {\bf e}_\varphi + {1\over \rho \sin\varphi}{\partial f \over \partial \theta} {\bf e}_\theta \end{eqnarray*} 簡單的量綱分析可以看出上式的合理性, 我們留給讀者去 看出來。 讀者當然可以由 \eqref{3.3} 直接代入驗證上式, 但是我們並不鼓勵這種被動的學習方式, 反而是從方向導數著手, 類似前面的推導方式自己走過一遍, 因為數學是非常的存在主義: 沒有人可以代替你相信!

圖7.$(r,\theta, z)$ 到 $(r +\triangle r, \theta +\triangle \theta, z +\triangle z)$ 之變化
圖8.$(\rho,\varphi,\theta)$到$(\rho +\triangle \rho, \varphi+\triangle \varphi, \theta +\triangle \theta)$ 之變化

四、 Gauss 散度定理 --- 散度

如果從散度的原始定義 \eqref{2.1}$_1$, 它只告訴我們微分而已甚麼東西都看不到。 要瞭解散度就必須從 Gauss 散度定理著手 \begin{equation} \iiint_V {\rm div} {\bf A} dV = \iint_S {\bf A} \cdot {\bf n} dS, \qquad S=\partial V \label{4.1} \tag{4.1} \end{equation} 兩邊同時除以體積 $|V|$, 然後再令 $|V|\to 0$ 得 \begin{equation} {\rm div} {\bf A} = \lim_{|V|\to 0}{1\over |V|}\iint_S {\bf A} \cdot {\bf n} dS \label{4.2} \tag{4.2} \end{equation} 所以散度 ${\rm div} {\bf A}$ 就是單位體積向量場 ${\bf A}$ 之通量 (flux)。 這個定義除了物理意義清楚之外, 也告訴我們散度與所選取之座標無關, 換言之 高斯散度定理提供我們絕佳的散度之定義。 $$ {散度}\quad\Longleftrightarrow\quad {通量(曲面積分)}\quad\Longleftrightarrow\quad {Gauss散度定理} $$ 簡單的量綱分析得 $$ [{\rm div} {\bf A}] ={1\over [|V|]}\iint_S [{\bf A} \cdot {\bf n} dS] ={1\over {\rm L}^3}[{\bf A}] {\rm L}^2= {[{\bf A}]\over {\rm L}} $$ 假設 ${\bf A}=A_1{\bf e}_1 +A_2 {\bf e}_2 + A_3 {\bf e}_3$, 我們先考慮 ${\bf e}_1$-方向之通量。 此時向量場由 $u_1$ 到 $u_1+\triangle u_1$, 利用 $\triangle S_1=(h_2\triangle u_2)(h_3 \triangle u_3)$ 得 \begin{eqnarray*} \iint_{S_1} {\bf A} \cdot {\bf n} dS &=& (A_1\triangle S_1)(u_1+\triangle u_1) - (A_1 \triangle S_1)(u_1)\\ &\approx& \left[(A_1h_2h_3)(u_1+\triangle u_1) -(A_1h_2h_3)(u_1)\right]\triangle u_2 \triangle u_3 \\ &=& {\partial \over \partial u_1}(A_1 h_2 h_3)\triangle u_1\triangle u_2 \triangle u_3\qquad {(均值定理)} \\ &=&{1\over h_1 h_2 h_3} {\partial \over \partial u_1}(A_1 h_2 h_3)h_1\triangle u_1 h_2\triangle u_2 h_3\triangle u_3 \end{eqnarray*} 這裡 $A_1, h_2, h_3$ 是 $u_1$-函數(因為 $u_2, u_3$ 固定)。 ${\bf e}_1$-方向有兩個面:前面與後面, 負號代表 後面。 同理可以考慮 ${\bf e}_2, {\bf e}_3$-方向之通量, 由於 $\triangle S_2=(h_3\triangle u_3)(h_1 \triangle u_1)$、 $\triangle S_3=(h_1\triangle u_1)(h_2 \triangle u_2)$ 得 \begin{eqnarray*} \iint_{S_2} {\bf A} \cdot {\bf n} dS &=& (A_2\triangle S_2)(u_2+\triangle u_2) - (A_2 \triangle S_2)(u_2)\\ &\approx& \left[(A_2h_3h_1)(u_2+\triangle u_2) -(A_2h_3h_1)(u_2)\right]\triangle u_3 \triangle u_1 \\ &=& {\partial \over \partial u_2}(A_2 h_3 h_1)\triangle u_2\triangle u_3 \triangle u_1 \\ &=&{1\over h_1 h_2 h_3} {\partial \over \partial u_2}(A_2 h_3 h_1)h_2\triangle u_2 h_3\triangle u_3 h_1\triangle u_1 \end{eqnarray*} 這裡 $A_2, h_3, h_1$ 視為 $u_2$-函數 (因為 $u_3, u_1$ 固定) \begin{eqnarray*} \iint_{S_3} {\bf A} \cdot {\bf n} dS &=& (A_3\triangle S_3)(u_3+\triangle u_3) - (A_3 \triangle S_3)(u_3)\\ &\approx& \left[(A_3h_1h_2)(u_3+\triangle u_3) -(A_3h_1h_2)(u_3)\right]\triangle u_1 \triangle u_2 \\ &=& {\partial \over \partial u_3}(A_3 h_1 h_2)\triangle u_3\triangle u_1\triangle u_2 \\ &=&{1\over h_1 h_2 h_3} {\partial \over \partial u_3}(A_3 h_1 h_2) h_3\triangle u_3 h_1\triangle u_1 h_2\triangle u_2 \end{eqnarray*} 這裡 $A_3, h_1, h_2$ 視為 $u_3$-函數(因為 $u_1, u_2$ 固定)。 將這三部分通量相加之後再除以體積
$h_1\triangle u_1 h_2\triangle u_2 h_3\triangle u_3$ 得 \begin{equation} {\rm div} {\bf A} = {1\over h_1 h_2 h_3}\left[ {\partial \over \partial u_1}(A_1 h_2 h_3) +{\partial \over \partial u_2}(A_2 h_3 h_1) +{\partial \over \partial u_3}(A_3 h_1 h_2)\right] \label{4.3} \tag{4.3} \end{equation} 只要記得方程式會講話則 \eqref{4.3} 實際上是非常直觀且容易記憶的。 以第一項 (${\bf e}_1$-方向)而言 $$ {1\over h_1 h_2 h_3}{\partial \over \partial u_1} (A_1 h_2 h_3) $$ 把 $h_1 h_2 h_3$ 視為體積、 $h_2 h_3$ 為底面積則 $A_1 h_2 h_3$ 是通量, 但有兩個面其 朝外法向量正好一正一負兩個通量必須相減, 也就是差分, 在無窮小的時候 ($u_1\to 0$) 差分變微分, 所以出現 ${\partial \over \partial u_1}$。 同理第二、三項也是如此分析。 表面上 \eqref{4.3} 雖然只有三項但本質上是六項, 正是無窮小方塊的六個面。 所以 \eqref{4.3} 告訴我們的正是: 單位體積之通量(三個方向)。

圖9.${\bf e}_i, i=1,2,3,$之通量

例題 4.1 圓柱座標之散度: ${\bf A}=A_r{\bf e}_r + A_\theta{\bf e}_\theta +A_z{\bf e}_z$ \begin{equation} {\rm div} {\bf A}= {1\over r}\left[{\partial \over \partial r} (A_r r)+{\partial \over \partial \theta} (A_\theta) +{\partial \over \partial z} (A_z r)\right] \label{4.4} \tag{4.4} \end{equation}

圖10.圓柱座標${\bf e}_r, {\bf e}_\theta,{\bf e}_z,$之通量

我們利用圖解$($請參考圖 $10)$ 來看出三個方向之通量: \begin{eqnarray*} \displaystyle {\bf e}_r:&& A_r \cdot r\cdot 1 ({r: 大}) - A_r \cdot r\cdot 1 ({r: 小}) \approx {\partial \over \partial r} (A_r r) \\ \displaystyle {\bf e}_\theta:&& A_\theta \cdot 1\cdot 1 ({\theta: 大}) - A_\theta \cdot 1\cdot 1 ({\theta: 小}) \approx {\partial \over \partial \theta} (A_\theta) \\ \displaystyle {\bf e}_z:&& A_z \cdot r\cdot 1 ({z: 大}) - A_z \cdot r\cdot 1 ({z: 小}) \approx {\partial \over \partial z} (A_z r) \end{eqnarray*} 再除以體積 $h_1h_2h_3=h_rh_\theta h_z=1\cdot r\cdot 1=r$ 就是 \eqref{4.4}。

例題 4.2 球座標之散度: $ {\bf A}=A_\rho{\bf e}_\rho + A_\varphi{\bf e}_\varphi +A_\theta{\bf e}_\theta$ \begin{equation} {\rm div} {\bf A}= {1\over \rho^2\sin\varphi}\left[ {\partial \over \partial \rho} (A_\rho \rho^2 \sin\varphi)+ {\partial \over \partial \varphi} (A_\varphi \rho\sin\varphi)+{\partial \over \partial \theta} (A_\theta \rho) \right] \label{4.5} \tag{4.5} \end{equation} 此時體積為 $h_1h_2h_3= h_\rho h_\varphi h_\theta=\rho^2\sin\varphi$, 而三個方向之通量為 $($請參考圖 $11, 12)$ \begin{eqnarray*} \displaystyle {\bf e}_\rho:&& A_\rho \cdot \rho\cdot \rho\sin \varphi ({\rho: 大}) - A_\rho \cdot \rho\cdot \rho\sin\varphi ({\rho: 小}) \approx {\partial \over \partial \rho} (A_\rho \rho^2 \sin\varphi) \\ \displaystyle {\bf e}_\varphi:&& A_\varphi \cdot \rho\sin\varphi\cdot 1 ({\varphi: 大}) - A_\varphi \cdot \rho\sin\varphi\cdot 1({\varphi: 小}) \approx {\partial \over \partial \varphi} (A_\varphi \rho\sin\varphi) \\ \displaystyle {\bf e}_\theta:&& A_\theta \cdot \rho\cdot 1 ({\theta: 大}) - A_\theta \cdot \rho\cdot 1 ({\theta: 小}) \approx {\partial \over \partial \theta} (A_\theta \rho)\hskip2.6cm \end{eqnarray*}

圖11.球座標 ${\bf e}_\rho$ 之通量
圖12.球座標 ${\bf e}_\varphi,{\bf e}_\theta$ 之通量

五、 Stokes 定理 --- 旋度

同樣從旋度的原始定義 \eqref{2.1}$_2$, 甚麼東西都看不到。 旋度直接聯想就是環流 (circulation) 或者是功 (work) 也就是線積分, 所以要瞭解旋度就必須從 Stokes 定理開始 \begin{equation} \oint_C {\bf A}\cdot d{\bf r} = \iint_S {curl } {\bf A} \cdot {\bf n} dS, \qquad C=\partial S \label{5.1} \tag{5.1} \end{equation} 要提醒讀者的是圍道$C$的方向與朝外法向量 ${\bf n}$ 必須滿足右手法則 (如果將 ${\bf n}$ 當作大拇指的方向則圍道 $C$ 是以逆時針方向為正)。 \eqref{5.1} 兩邊同時除以曲面面積 $|S|$, 然後再令 $|S|\to 0$ 得 \begin{equation} {curl } {\bf A} \cdot {\bf n} = \lim_{|S|\to 0}{1\over |S|}\oint_C {\bf A} \cdot d {\bf r} \label{5.2} \tag{5.2} \end{equation} 因為 $|{curl } {\bf A} \cdot {\bf n}|\le |{curl } {\bf A}|$, 所以旋度 ${curl } {\bf A}$ 可以詮釋為單位面積向量場 ${\bf A}$ 之最大環流。 由此可知旋度 ${curl {\bf A}}$ 之量綱 (或因次; dimension) 等於 ${\bf A}$ 之量綱除以長度 $$ [{curl {\bf A}}]= {[{\bf A}\cdot {\bf t}] {\rm L} \over {\rm L}^2} = {[{\bf A}]\over {\rm L}} $$ 因為線積分與座標無關, 所以旋度與所選取之座標無關, 換言之 Stokes 定理提供 我們絕佳的旋度之定義 $$ {旋度}\qquad \Longleftrightarrow \qquad {環流(線積分)}\qquad \Longleftrightarrow \qquad {Stokes定理} $$

假設 ${\bf A}=A_1{\bf e}_1 +A_2 {\bf e}_2 + A_3 {\bf e}_3$。 我們先考慮 ${\bf e}_1$-方向, 此時 $S_1$ 是前面 $$ {\bf e}_1\perp S_1, \qquad C_1=\partial S_1,\qquad |S_1|=(h_2\triangle u_2)( h_3\triangle u_3) $$ 按照右手法則 $C_1$ 有四個邊, 所以線積分有四項 \begin{eqnarray*} \oint_{C_1} {\bf A}\cdot d{\bf r} &\approx& (A_2h_2)(u_3)\triangle u_2 +(A_3h_3)(u_2+\triangle u_2)\triangle u_3\\ &&-(A_2h_2)(u_3+\triangle u_3)\triangle u_2 -(A_3h_3)(u_2)\triangle u_3 \\ &=&\left({\partial (A_3h_3)\over\partial u_2} - {\partial (A_2h_2)\over\partial u_3}\right)\triangle u_2 \triangle u_3 \\ &=&{1\over h_2 h_3}\left({\partial (A_3h_3)\over\partial u_2} - {\partial (A_2h_2)\over\partial u_3}\right)(h_2\triangle u_2)(h_3\triangle u_3) \end{eqnarray*} 這裡 $A_3h_3, A_2h_2$ 分別視為 $u_2$ 與 $u_3$ 的函數。 這個線積分再除以 $S_1$ 的面積 $|S_1|$ 得 $\nabla\times{\bf A}$ 在 ${\bf e}_1$-方向的分量 \begin{equation} \nabla\times{\bf A}|_{{\bf e}_1}={curl } {\bf A}|_{{\bf e}_1}= {1\over h_2 h_3}\left({\partial (A_3h_3)\over\partial u_2} - {\partial (A_2h_2)\over\partial u_3}\right) \label{5.3} \tag{5.3} \end{equation} 同理可以得 ${\bf e}_2, {\bf e}_3$-方向的分量。 令 $S_1$ 是 右面, 則 $$ {\bf e}_2\perp S_2, \qquad C_2=\partial S_2,\qquad |S_2|=(h_3\triangle u_3)( h_1\triangle u_1) $$ \begin{eqnarray*} \oint_{C_2} {\bf A}\cdot d{\bf r} &\approx& (A_3h_3)(u_1)\triangle u_3 +(A_1h_1)(u_3+\triangle u_3)\triangle u_1\\ && -(A_3h_3)(u_1+\triangle u_1)\triangle u_3 -(A_1h_1)(u_3)\triangle u_1 \\ &=&\left({\partial (A_1h_1)\over\partial u_3} - {\partial (A_3h_3)\over\partial u_1}\right)\triangle u_3 \triangle u_1 \\ &=&{1\over h_3 h_1}\left({\partial (A_1h_1)\over\partial u_3} - {\partial (A_3h_3)\over\partial u_1}\right)(h_3\triangle u_3)(h_1\triangle u_1) \end{eqnarray*} 這裡 $A_1h_1, A_3h_3$ 分別視為 $u_3$ 與 $u_1$ 的函數。 所以 $\nabla\times{\bf A}$ 在 ${\bf e}_2$-方向( 右面)的分量為 \begin{equation} \nabla\times{\bf A}|_{{\bf e}_2}={curl } {\bf A}|_{{\bf e}_2}= {1\over h_3 h_1}\left({\partial (A_1h_1)\over\partial u_3} - {\partial (A_3h_3)\over\partial u_1}\right) \label{5.4} \tag{5.4} \end{equation} 令 $S_3$ 是 上面, 則 $$ {\bf e}_3\perp S_3, \qquad C_3=\partial S_3,\qquad |S_3|=(h_1\triangle u_1)( h_2\triangle u_2) $$ \begin{eqnarray*} \oint_{C_3} {\bf A}\cdot d{\bf r} &\approx& (A_1h_1)(u_2)\triangle u_1 +(A_2h_2)(u_1+\triangle u_1)\triangle u_2\\ &&-(A_1h_1)(u_2+\triangle u_2)\triangle u_1 -(A_2h_2)(u_1)\triangle u_2 \\ &=&\left({\partial (A_2h_2)\over\partial u_1} - {\partial (A_1h_1)\over\partial u_2}\right)\triangle u_1 \triangle u_2 \\ &=&{1\over h_1 h_2}\left({\partial (A_2h_2)\over\partial u_1} - {\partial (A_1h_1)\over\partial u_2}\right)(h_1\triangle u_1)(h_2\triangle u_2) \end{eqnarray*} 這裡 $A_2h_2, A_1h_1$ 分別視為 $u_1$ 與 $u_2$ 的函數。 所以 $\nabla\times{\bf A}$ 在 ${\bf e}_3$-方向( 右面)的分量為 \begin{equation} \nabla\times{\bf A}|_{{\bf e}_3}={curl } {\bf A}|_{{\bf e}_3}= {1\over h_1 h_2}\left({\partial (A_2h_2)\over\partial u_1} - {\partial (A_1h_1)\over\partial u_2}\right) \label{5.5} \tag{5.5} \end{equation} 將 \eqref{5.3}--\eqref{5.5} 合併得 \begin{eqnarray} \displaystyle {curl {\bf A}} &=&{1\over h_2 h_3} \Big( {\partial \over \partial u_2}(A_3h_3) -{\partial \over \partial u_3}(A_2h_2)\Big) {\bf e}_1\hskip2cm \nonumber\\ &&\displaystyle +{1\over h_3 h_1} \Big( {\partial \over \partial u_3}(A_1h_1) -{\partial \over \partial u_1}(A_3h_3)\Big) {\bf e}_2 \nonumber\\ &&\displaystyle +{1\over h_1 h_2} \Big( {\partial \over \partial u_1}(A_2h_2) -{\partial \over \partial u_2}(A_1h_1)\Big) {\bf e}_3\label{5.6} \tag{5.6} \end{eqnarray} 再由 Laplace 降階法結論 \begin{equation} {curl {\bf A}}=\nabla\times{\bf A} ={1\over h_1 h_2 h_3}\left|\begin{array}{ccc} h_1{\bf e}_1 & h_2{\bf e}_2 & h_3{\bf e}_3\\ \frac{\partial}{\partial u_1} & \frac{\partial}{\partial u_2} & \frac{\partial}{\partial u_3}\\ A_1h_1 & A_2h_2 & A_3h_3 \end{array}\right| \label{5.7} \tag{5.7} \end{equation}

註解:

(i) 在計算過程中雖然我們只考慮三個面, 但一共有十二個邊, 因此這個立方體所有面之邊界的線積分都已經考慮到了, 所以最後之結果確實是全部之環流。

(ii) 因為$h_i$必須配合$\partial u_i$, 為了保持量綱平衡${\bf A}$的旋度, 例如第一項(${\bf e}_1$方向) $$ {1\over h_2 h_3} \Big( {\partial \over \partial u_2}(A_3h_3) -{\partial \over \partial u_3}(A_2h_2)\Big){\bf e}_1 $$ 分子分母的$h_3$、$h_2$分別彼此抵消而$\frac1{h_2}$、$\frac1{h_3}$則正好配合 $\frac{\partial }{\partial u_2}$、$\frac{\partial }{\partial u_3}$, 而且可以驗證其量綱(因次)等於$[A_1]/{\rm L}$ 與 $[{\rm curl }{\bf A}]=[{\bf A}]/{\rm L}$ 相同, 至於正負符號則由右手法則 來決定 \begin{eqnarray*} {\partial \over \partial u_2}(A_3h_3){\bf e}_1&\Longrightarrow& (2\to 3\to 1) \quad (+)\\ {\partial \over \partial u_3}(A_2h_2){\bf e}_1&\Longrightarrow& (3\to 2\to 1) \quad (-) \end{eqnarray*}

(iii) 由Stokes定理知道旋度是單位面積之環流, 事實上(5.6) 告訴我們的正是如此, 看第一項(${\bf e}_1$-方向, ${\bf e}_1$為法向量的面) $$ {1\over h_2 h_3} \Big( {\partial \over \partial u_2}(A_3h_3) -{\partial \over \partial u_3}(A_2h_2)\Big) $$ 此時 ${\bf e}_3$-方向之線積分(環流)為 $A_3h_3(u_2{大})$ 與 $-A_3h_3(u_2{小})$, 當$u_2\to 0$時這兩項合併得 ${\partial \over \partial u_2}(A_3h_3)$, ( 差分變微分!)。 同理 ${\bf e}_2$-方向之線積分(環流)為 $-A_2h_2(u_3{大})$ 與 $A_2h_2(u_2{小})$, 當 $u_3\to 0$時這兩項合併得 $-{\partial \over \partial u_3}(A_2h_2)$, 因此偏導數是說兩個線積分之差(正負兩個方向), $h_2 h_3$視為這個長方形之面積, 整個合併起來就是${\bf e}_1$-方向的單位面積之環流。同理可以討論${\bf e}_2$-方向與${\bf e}_3$-方向。

我們利用圖解(請參考圖 13)來看出三個方向之環流: \begin{eqnarray*} \oint_{C_1} {\bf A}\cdot d{\bf r} &\approx& A_2h_2 (u_3 {小})+A_3 h_3(u_2 {大})-A_2h_2 (u_3 {大})-A_3 h_3(u_2 {小})\\ &\approx& {\partial \over \partial u_2}(A_3h_3)-{\partial \over \partial u_3}(A_2h_2)\\ \oint_{C_2} {\bf A}\cdot d{\bf r} &\approx& A_3h_3 (u_1 {小})+A_1 h_1(u_3 {大})-A_3h_3 (u_1 {大})-A_1 h_1(u_3 {小}) \\ &\approx& {\partial \over \partial u_3}(A_1h_1)-{\partial \over \partial u_1}(A_3h_3) \\ \oint_{C_3} {\bf A}\cdot d{\bf r} &\approx& A_1h_1 (u_2 {小})+A_2 h_2(u_1 {大})-A_1h_1 (u_2 {大})-A_2 h_2(u_1 {小}) \\ &\approx& {\partial \over \partial u_1}(A_1h_1)-{\partial \over \partial u_2}(A_1h_1) \end{eqnarray*} 再分別除以 $h_2h_3, h_3 h_1, h_1 h_2$ 就是 \eqref{5.3}--\eqref{5.5}。

圖13.${\bf e}_i,i=1,2,3$ 三個方向之環流

例題 5.1 圓柱座標之旋度: ${\bf A}=A_r{\bf e}_r + A_\theta{\bf e}_\theta +A_z{\bf e}_z$ \begin{equation} {{\rm curl} {\bf A}}=\frac1r\left[\frac{\partial A_z}{\partial\theta} -r\frac{\partial A_\theta}{\partial z}\right]{\bf e}_r +\left[\frac{\partial A_r}{\partial z} -\frac{\partial A_z}{\partial r}\right]{\bf e}_\theta +\frac1r\left[\frac{\partial}{\partial r}(rA_\theta) -\frac{\partial A_r}{\partial\theta}\right]{\bf e}_z \label{5.8} \tag{5.8} \end{equation}

我們利用圖解 (請參考圖14) 推得三個方向之環流: \begin{eqnarray*} \oint_{C_r} {\bf A}\cdot d{\bf r} &\approx& A_\theta r (z {小})+A_z\cdot 1(\theta {大})-A_\theta r (z {大})-A_z (\theta {小}) \\ &\approx& {\partial \over \partial \theta}(A_z)-{\partial \over \partial z}(A_\theta r) \\ \oint_{C_\theta} {\bf A}\cdot d{\bf r} &\approx& A_z (r {小})+A_r (z {大})-A_z (r {大})-A_r (z {小}) \\ &\approx& {\partial \over \partial z}(A_r)-{\partial \over \partial r}(A_z) \\ \oint_{C_z} {\bf A}\cdot d{\bf r} &\approx& A_r (\theta {小})+A_\theta r(r {大})-A_r (\theta {大})-A_\theta r (r {小}) \\ &\approx& {\partial \over \partial r}(A_\theta r )-{\partial \over \partial \theta}(A_r) \end{eqnarray*} 再分別除以 $r, 1, r$ $($視為 $S_r, S_\theta, S_z$ 之面積$)$ 就是 \eqref{5.8}。

圖14.圓柱座標三個方向之環流

例題 5.2 球座標之旋度: ${\bf A}=A_\rho{\bf e}_\rho + A_\varphi{\bf e}_\varphi +A_\theta{\bf e}_\theta$ \begin{eqnarray} \displaystyle {{\rm curl} {\bf A}} &=&\frac1{\rho^2\sin\varphi}\left[\rho\frac{\partial}{\partial\varphi}(\sin\varphi A_\theta) -\rho\frac{\partial A_\varphi}{\partial\theta}\right]{\bf e}_\rho\nonumber\\ &&\displaystyle +\frac1{\rho\sin\varphi}\left[\frac{\partial A_\rho}{\partial\theta} -\sin\varphi\frac{\partial}{\partial \rho}(\rho A_\theta)\right]{\bf e}_\varphi +\frac1\rho\left[\frac{\partial}{\partial \rho}(\rho A_\varphi) -\frac{\partial A_\rho}{\partial\varphi}\right]{\bf e}_\theta\label{5.9} \tag{5.9} \end{eqnarray}

我們利用圖解 (請參考圖15) 推得三個方向之環流: \begin{eqnarray*} \oint_{C_\rho} {\bf A}\cdot d{\bf r} &\approx& A_\varphi \rho (\theta {小})+A_\theta \rho \sin \varphi(\varphi{大})-A_\varphi \rho (\theta {大})-A_\theta \rho\sin \varphi (\varphi {小}) \\ &\approx& {\partial \over \partial \varphi}(A_\theta \rho \sin \varphi)-{\partial \over \partial \theta}(A_\varphi \rho) \\ \oint_{C_\varphi} {\bf A}\cdot d{\bf r} &\approx& A_\theta \rho \sin \varphi (\rho {小})+A_\rho (\theta {大})-A_\theta \rho \sin \varphi (\rho {大})-A_\rho (\theta {小}) \\ &\approx& {\partial \over \partial \theta}(A_\rho)-{\partial \over \partial \rho}(A_\theta \rho \sin\varphi) \\ \oint_{C_\theta} {\bf A}\cdot d{\bf r} &\approx& A_\rho (\varphi {小})+A_\varphi \rho(\rho {大})-A_\rho (\varphi {大})-A_\varphi \rho (\rho {小}) \\ &\approx& {\partial \over \partial \rho}(A_\varphi \rho)-{\partial \over \partial \varphi}(A_\rho) \end{eqnarray*} 再分別除以 $\rho^2\sin\varphi, \rho\sin\varphi, \rho$ $($視為 $S_\rho, S_\varphi, S_\theta$ 之面積$)$ 就是 \eqref{5.9}。

圖15.球座標三個方向之環流

六、 Laplace 算子

根據 Laplace 算子的定義 \begin{equation} \Delta= {\rm div} \nabla=\nabla\cdot\nabla \label{6.1} \tag{6.1} \end{equation} 所以由梯度與散度的公式 \eqref{3.3} 與 \eqref{4.3} 得 \begin{eqnarray} \displaystyle \Delta f&=&{\rm div} \bigg({1\over h_1}{\partial f \over \partial u_1}{\bf e}_1 +{1\over h_2}{\partial f \over \partial u_2}{\bf e}_2 +{1\over h_3}{\partial f \over \partial u_3}{\bf e}_3\bigg)\hskip4.3cm \nonumber\\ \displaystyle &=& {1\over h_1h_2h_3}\bigg[{\partial \over \partial u_1}\Big({h_2h_3\over h_1} {\partial f\over \partial u_1}\Big) +{\partial \over \partial u_2}\Big({h_3h_1\over h_2} {\partial f\over \partial u_2}\Big) +{\partial \over \partial u_3}\Big({h_1h_2\over h_3} {\partial f\over \partial u_3}\Big)\bigg] \label{6.2} \tag{6.2} \end{eqnarray} 雖然這樣可以得到答案, 但對於瞭解 Laplace 算子並沒有幫助, 頂多是散度與梯度的組合。 回到 Gauss 散度定理 \begin{equation} \iiint_V \Delta f dV =\iiint_V {\rm div} (\nabla f) dV \\ =\iint_S \nabla f\cdot {\bf n} dS =\iint_S {\partial f\over \partial {\bf n}} dS \label{6.3} \tag{6.3} \end{equation} 兩邊同時除以體積 $|V|$, 然後再令 $|V|\to 0$ 得 \begin{equation} \Delta f= \lim_{|V|\to 0}{1\over|V|}\iint_S {\partial f\over \partial {\bf n}} dS \label{6.4} \tag{6.4} \end{equation} 所以 $\Delta f$ 就是單位體積 $f$ 之 法向量的方向導數 ${\partial f\over \partial {\bf n}}$ 之通量。 量綱分析得 $$ [\Delta f]={[f]\over {\rm L}^2} = {1\over[|V|]} \left[{\partial f\over \partial {\bf n}}\right] [dS] ={1\over {\rm L}^3} {[f]\over {\rm L}} {\rm L}^2={[f]\over {\rm L}^2} $$ 我們利用圖解來推得 \eqref{6.2}。 ${\partial f\over \partial {\bf n}}$ 在三個方向之通量為 $$ {{\bf e}_1-方向(S_1曲面):}\quad {1\over h_1} {\partial f\over \partial u_1} h_2h_3 (u_1{大})- {1\over h_1} {\partial f\over \partial u_1} h_2h_3 (u_1{小}) \approx {\partial \over \partial u_1}\Big({h_2h_3\over h_1} {\partial f\over \partial u_1} \Big) $$ $$ {{\bf e}_2-方向(S_2曲面):}\quad {1\over h_2} {\partial f\over \partial u_2} h_3h_1 (u_2{大})- {1\over h_2} {\partial f\over \partial u_2} h_3h_1 (u_2{小}) \approx {\partial \over \partial u_2}\Big({h_3h_1\over h_2} {\partial f\over \partial u_2} \Big) $$ $$ {{\bf e}_3-方向(S_3曲面):}\quad {1\over h_3} {\partial f\over \partial u_3} h_1h_2 (u_3{大})- {1\over h_3} {\partial f\over \partial u_3} h_1h_2 (u_3{小}) \approx {\partial \over \partial u_3}\Big({h_1h_2\over h_3} {\partial f\over \partial u_3} \Big) $$ 再除以 $h_1 h_2 h_3$ (視為無窮小方塊之體積) 就是 \eqref{6.2}。 類似於散度 \eqref{6.2} 可以詮釋如下: \begin{eqnarray*} {1\over h_1 h_2 h_3}\quad&\Longrightarrow&\quad{視為單位體積}\\ h_2 h_3, h_3h_1, h_1h_2\quad&\Longrightarrow&\quad{視為S_1, S_2, S_3-之面積}\\ {1\over h_1} {\partial f\over \partial u_1}, {1\over h_2} {\partial f\over \partial u_2}, {1\over h_3} {\partial f\over \partial u_3} \quad&\Longrightarrow&\quad{{\bf e}_1,{\bf e}_2,{\bf e}_3-方向之梯度}\\ {h_2h_3\over h_1} {\partial f\over \partial u_1}, {h_3h_1\over h_2} {\partial f\over \partial u_2}, {h_1h_2\over h_3} {\partial f\over \partial u_3} \quad&\Longrightarrow&\quad{{\bf e}_1,{\bf e}_2,{\bf e}_3-方向之通量}\\ {\partial \over \partial u_i},\quad i=1,2,3\quad&\Longrightarrow&\quad{差分變微分} \end{eqnarray*}

例題 6.1 極座標之 ${\rm Laplace}$ 算子 \begin{equation} \Delta f= {\partial^2 f \over \partial x^2}+{\partial^2 f \over \partial y^2} ={\partial^2 f \over \partial r^2} +{1\over r} {\partial f \over \partial r}+{1\over r^2}{\partial^2 f \over \partial \theta^2} \label{6.5} \tag{6.5} \end{equation} 我們仍然利用圖解 $($請參考圖16 推 \eqref{6.5} $$ {\partial f\over \partial r} r (r{大}) -{\partial f\over \partial r} r (r{小}) \approx {\partial \over \partial r}\Big(r {\partial f \over \partial r}\Big) $$ $$ {1\over r}{\partial f\over \partial \theta}(\theta {大}) -{1\over r}{\partial f\over \partial \theta}(\theta {小}) \approx {1\over r} {\partial^2 f\over \partial \theta^2} $$ 這兩式相加之後除以 $r$ 得 $$ \Delta f = {1\over r}{\partial \over \partial r}\Big(r {\partial f \over \partial r}\Big)+{1\over r^2}{\partial^2 f\over \partial \theta^2} ={\partial^2 f \over \partial r^2} +{1\over r} {\partial f \over \partial r}+{1\over r^2}{\partial^2 f \over \partial \theta^2} $$

圖16.極座標之 Laplace 算子

例題 6.2 圓柱座標之 ${\rm Laplace}$ 算子 \begin{equation} \Delta f= {\partial^2 f \over \partial x^2}+{\partial^2 f \over \partial y^2} +{\partial^2 f \over \partial z^2} ={\partial^2 f \over \partial r^2} +{1\over r} {\partial f \over \partial r}+{1\over r^2}{\partial^2 f \over \partial \theta^2}+{\partial^2 f \over \partial z^2} \label{6.6} \tag{6.6} \end{equation}

圓柱座標幾乎就是極座標我們留給讀者練習。

例題 6.3 球座標之 ${\rm Laplace}$ 算子 \begin{eqnarray} \displaystyle \Delta f&=&{1\over \rho^2 \sin\varphi}\bigg[{\partial \over \partial \rho}\bigg({\partial f\over \partial \rho}\rho^2 \sin \varphi\bigg)+{\partial \over \partial \varphi}\bigg({\partial f\over \partial \varphi}\sin \varphi\bigg)+{\partial \over \partial \theta}\bigg({1\over \sin \varphi}{\partial f\over \partial \theta}\bigg)\bigg] \nonumber\\ \displaystyle &=& {1\over \rho^2}{\partial \over \partial \rho}\bigg(\rho^2{\partial f\over \partial \rho}\bigg)+{1\over \rho^2 \sin\varphi}{\partial \over \partial \varphi}\bigg(\sin \varphi{\partial f\over \partial \varphi}\bigg)+ {1\over \rho^2\sin^2\varphi} {\partial^2 f\over \partial \theta^2}\quad\label{6.7} \tag{6.7} \end{eqnarray}

讀者可以驗證一下量綱關係 \begin{eqnarray*} {1\over \rho^2}{\partial \over \partial \rho}\bigg(\rho^2{\partial f\over \partial \rho}\bigg)&\approx& {1\over {\rm L}^2} {1\over {\rm L}}\bigg({\rm L}^2 {[f]\over {\rm L}}\bigg) ={[f]\over {\rm L}^2}\\ {1\over \rho^2 \sin\varphi}{\partial \over \partial \varphi}\bigg(\sin \varphi{\partial f\over \partial \varphi}\bigg) &\approx&{1\over {\rm L}^2 \cdot 1} {1\over 1}\bigg( 1 {[f]\over 1}\bigg) ={[f]\over {\rm L}^2} \\ {1\over \rho^2\sin^2\varphi} {\partial^2 f\over \partial \theta^2} &\approx& {1\over {\rm L}^2 1^2} {[f]\over 1^2}= {[f]\over {\rm L}^2} \end{eqnarray*} 首先由例題 3.1 得 $$ \nabla f= {\partial f\over \partial \rho}{\bf e}_\rho + {1\over \rho}{\partial f\over \partial \varphi}{\bf e}_\varphi +{1\over \rho \sin\varphi}{\partial f\over \partial \theta}{\bf e}_\theta $$ 所以 \eqref{6.7} 圖解如下: $$ {\partial f\over \partial \rho}\cdot\rho \cdot\rho\sin \varphi (\rho{大})- {\partial f\over \partial \rho}\cdot\rho \cdot\rho\sin \varphi (\rho{小}) \approx {\partial \over \partial \rho}\bigg({\partial f\over \partial \rho}\rho^2 \sin \varphi\bigg) $$ $$ {1\over \rho}{\partial f\over \partial \varphi} \cdot\rho\sin \varphi\cdot 1 (\varphi{大})- {1\over \rho}{\partial f\over \partial \varphi} \cdot\rho\sin \varphi\cdot 1 (\varphi{小}) \approx {\partial \over \partial \varphi}\bigg({\partial f\over \partial \varphi}\sin \varphi\bigg) $$ $$ {1\over \rho \sin\varphi}{\partial f\over \partial \theta}\cdot 1 \cdot\rho (\theta{大})- {1\over \rho \sin\varphi}{\partial f\over \partial \theta}\cdot 1 \cdot\rho (\theta{小}) \approx {\partial \over \partial \theta}\bigg({1\over \sin \varphi}{\partial f\over \partial \theta}\bigg) $$ 這三式相加再除以 $\rho^2 \sin \varphi$ 就是 \eqref{6.7}。

圖17.$n$ 維球座標之 Laplace 算子

例題 6.4 $n$ 維球座標之 ${\rm Laplace}$ 算子 \begin{equation} \Delta f = {1\over \omega_n r^{n-1}}\left[{\partial \over \partial r}\bigg({\partial f \over \partial r}\omega_n r^{n-1}\bigg)+\cdots\right] ={\partial^2 f \over \partial r^2} +{n-1\over r} {\partial f\over \partial r} + {1\over r^2}\Delta^* f \label{6.8} \tag{6.8} \end{equation} 這裡 $\omega_n$ 是 $n$ 維空間的單位球面之表面積 \begin{equation} \omega_n=|S^{n-1}(1)|={2\pi^{n\over 2} \over \Gamma({n\over 2})} \label{6.9} \tag{6.9} \end{equation} 而 $\Delta^*$ 則是單位球面 $S^{n-1}(1)$ 上的 ${\rm Laplace}$ 算子。

$n$ 維球座標其實不難想像。 與極座標的扇形方塊類比我們可以想像此時是一個圓錐形的冰淇淋半徑方向的長度還是 $dr$, 弧長 $rd\theta$ 則取代為半徑等於 $r$ 之球面的一部分其表面積可以設為 $k\omega_n r^{n-1}$ ($k$ 與頂角有關), 因此半徑方向之法向方向導數 ${\partial f\over \partial r}$ 的通量(請參考圖 17) $$ {\partial f \over \partial r}k\omega_n r^{n-1} (r{大})- {\partial f \over \partial r}k\omega_n r^{n-1} (r{小}) \approx k{\partial \over \partial r}\bigg({\partial f \over \partial r}\omega_n r^{n-1}\bigg) $$ 再除以體積 $1\cdot k\omega_n r^{n-1}$ 就是 \eqref{6.8} 第一個等式右邊的第一項。

註解:

(i) 從計算過程中可以看得出來在 \eqref{6.8} 出現的 $n-1$ 是半徑為 $r$ 的球表面積對 $r$ 微分的結果, 其實就是球面的維數 (dimension)。

(ii) 當 $n=2,3$ 時由 \eqref{6.5}, \eqref{6.7} 得 \begin{equation} \Delta^* f = {\partial^2 f \over \partial \theta^2},\qquad \Delta^* f = {\partial^2 f \over \partial \varphi^2} +{\cos \varphi\over \sin \varphi} {\partial f \over \partial \varphi}+{1\over \sin^2\varphi} {\partial^2 f \over \partial \theta^2} \label{6.10} \tag{6.10} \end{equation}

誌謝

本文是作者於2009年暑假在中研院數學所 《數學名題及其故事》 一系列演講之一。 後來更在中研院數學所星期一的分析討論班中詳細地以圖解的方式講述了一次, 在此特別謝謝李志豪教授的鼓勵與安排。

參考文獻

Michael J. Crowe, A History of Vector Analysis, The evolution of the Idea of a Vectorial System, Dover Publications, Inc. 1985.

學任何ㄧ門學問正如人生一樣, 偶而要問一下人生三大問題: 「我從哪裡來?我要往哪裡去?我現在活著要做什麼?」這樣才不致於迷失。 ㄧ門學問的歷史與其本身一樣重要, 從歷史中去揣摩前人想要解決的問題, 為何這問題重要; 還有他們的解決方法, 但最重要的還是去體會他們的心路歷程。 這本書對於向量分析的歷史介紹幾乎是百科全書式的, 裡面有Maxwell, Gibbs, $\ldots$ 等人的引言, 讀起來樂趣無窮, 我個人是極力推薦。

H. M. Schey, Div, Grad, Curl and all that, An Informal Text on Vector Calculus, Fourth Edition, W. W. Norton & Company Inc. 2005.

這本奇妙的小書具有武功祕笈的特色, 作者是學電機的, 說起來實在是讓我們這些數學科班出身的覺得慚愧, 書中以圖形直觀的介紹梯度、散度與旋度, 是整本書最引人入勝之處, 這也是我個人真正對梯度、散度與旋度有深入了解的開始。 這本書原先有中譯本「誰怕向量微積分」但因版權問題已不再出版。

Freeman Dyson, Missed Opportunities, Bulletin of Americal Mathematical Society, Vol. 78, 636--652, 1972.

這是著名物理學家Freeman Dyson 在 Gibbs 講座的演講紀錄, Dyson提了不少數學界所曾失去的大好機會, 最後他也提了一些值得數學家投入的好問題。讀數學的人應該與外界的人交流, 你會發現人家對於數學的要求與期待是何等地殷切。如果整天躲在數學館, 只會自艾自憐不知數學要何去何從。

林琦焜, 向量分析, 滄海書局, 2007。 林琦焜, Green定理與應用, 數學傳播(中央研究院數學所), Vol.~84, 25-41 (1997)。 林琦焜, 從量綱看世界, 數學傳播(中央研究院數學所), Vol.~131, p. 13--27 (2009)。 邵錦昌, 淺談 Stokes' 定理與電磁學, 數學傳播(中央研究院數學所), Vol.~72, 1-13 (1994)。 蔡聰明, 從醉月湖的面積談起---向量微積分簡介, 數學傳播(中央研究院數學所), Vol.~82, 3-16 (1997)。

我極力建議學生在大學部時就培養閱讀期刊的習慣, 上面所列是數學傳播裡面與向量分析有關聯的幾篇文章, 有興趣的學生可以上中央研究院數學所的數學傳播網站搜尋相關文獻, 另外有一份英文期刊, American Mathematical Monthly, 也是一份相當適合大學部學生的期刊。

---本文作者為交通大學應用數學系退休教授---