兩個不肯就範的古老問題
先從一個你以為早已掌握的問題入手:橢圓的周長。半徑為 r 的圓,周長是 2 pi r——乾淨、初等,在大一微積分裡就解決了。於是你搬出弧長公式,代入橢圓 x = a cos theta、y = b sin theta,自信地寫下這個定積分。結果是從 0 到 2 pi 的 square root of (a^2 sin^2 theta + b^2 cos^2 theta) d theta。然後一切就卡住了。那個根號沒有初等原函數。橢圓——僅次於圓的最基本曲線——它的周長,竟無法用多項式、根式、指數、對數和三角函數的任何有限組合來表示。
鐘擺講的是同一個故事。你最初遇到小角度擺,是把它當作簡諧振子:theta'' = -(g/L) sin theta,而課本悄悄地把 sin theta 換成 theta,得到那個乾淨的週期 2 pi square root of (L/g),且與擺幅無關。但對任何真實的擺動,這都是個謊言。保留真正的 sin theta,去問一個從振幅 theta_0 釋放的擺的精確週期,你又會落到形如 1 over square root of (某物 - 某物 times sin^2) 的積分上。誠實的鐘擺和誠實的橢圓,交給你的是同一個形狀——而這個形狀有自己的名字。
為答案命名:勒讓德的四種標準形
當微積分拒絕給出公式時,成熟的做法是為我們所需的函數命名——就像當年對數被發明為 1/x 的原函數,又像上一篇裡 erf 為鐘形曲線而命名一樣。這裡這一族棘手的積分範圍很廣:任何關於 theta 與一個三次或四次式的平方根的有理函數的積分,都叫作橢圓積分——參見橢圓積分(非初等)。勒讓德證明,其中每一個都可以透過標準代換化歸為屈指可數的幾個基本積木。我們只需兩種基本形狀。
這兩種形狀就是第一類和第二類橢圓積分,寫法裡帶一個自由的上限 phi(稱為幅角)和一個介於 0 和 1 之間的數 k(稱為模,它決定根號的離心程度)。讓上限 phi 保持自由,就得到不完全積分;把上限一路推到 pi/2——整整一個四分之一週期——就得到完全積分 K(k) 和 E(k)。這樣總共有四個命名對象:不完全的 F 與 E,完全的 K 與 E。記號略有重疊——第二類在完全和不完全兩種情形都用 E——所以一定要看清是否帶有幅角參數。
Incomplete (upper limit phi is free):
F(phi, k) = integral_0^phi d theta / sqrt(1 - k^2 sin^2 theta) (first kind)
E(phi, k) = integral_0^phi sqrt(1 - k^2 sin^2 theta) d theta (second kind)
Complete (set phi = pi/2, a full quarter turn):
K(k) = F(pi/2, k) = integral_0^{pi/2} d theta / sqrt(1 - k^2 sin^2 theta)
E(k) = E(pi/2, k) = integral_0^{pi/2} sqrt(1 - k^2 sin^2 theta) d theta
Endpoints: K(0) = E(0) = pi/2 ; E(1) = 1 ; K(k) -> +infinity as k -> 1一個把它繫回大一微積分的合理性檢驗:令 k = 0。根號塌縮為 1,於是 F(phi, 0) = integral from 0 to phi of d theta = phi,而 E(phi, 0) 也等於 phi。所以在模為零時,橢圓積分不過是恆等函數——它們本就是初等情形,只是換了身打扮。模 k 恰好度量你離開初等世界有多遠。k 很小意味著近乎圓形的橢圓或極小的擺幅,你可以把熟悉的答案作為關於 k 的級數的首項重新找回來。
讓它們上崗:鐘擺與橢圓
鐘擺的回報來了,用一幅圖來推。用能量守恆把角 theta 處的角速度表示出來,分離變量,再對四分之一擺程積分。經過那個乾淨的代換 sin(theta/2) = sin(theta_0/2) sin(phi)——這一招與三角換元同屬一族——那個雜亂的四分之一擺程積分恰好折疊成 F(pi/2, k),也就是 K(k),其中模 k = sin(theta_0/2)。
於是精確週期為 T = 4 square root of (L/g) times K(sin(theta_0/2))。讀懂它在說什麼。當振幅 theta_0 很小,k 接近 0,K(0) = pi/2,T 便塌回那個熟悉的 2 pi square root of (L/g)——小角度答案作為極限情形自然落出,而非偶然。但當你讓鐘擺擺得更寬,k 增大,K(k) 增大,真實週期就比課本值更長:從 90 度釋放的擺,每個週期比小角度公式預測的約長 18%。當 theta_0 趨向 180 度(在頂端達到平衡)時,k 趨向 1,K(k) 對數式發散,週期奔向無窮——這恰恰正確,因為一個在頂端豎直平衡的擺,永遠不會回來。
橢圓則是第二類的故事。它的周長是 4 a times E(e),其中 a 是半長軸,e 是離心率,E 是完全第二類積分——即完全橢圓積分。圓是 e = 0 的情形,此時 E(0) = pi/2,周長為 4 a times pi/2 = 2 pi a,正是你早已熟悉的周長。而從橢圓頂端到某一般點的弧長,則是一個不完全第二類積分 E(phi, e)——參見不完全橢圓積分——因為此時上限只走到一半就停了。「完全」對「不完全」,無非就是「整整一個四分之一」對「你走到哪兒了」。
如何計算它們,以及別怕「非初等」這個詞
人們很容易覺得非初等積分低人一等——覺得只要答案不是個整潔公式,自己就失敗了。把這種感覺扔掉。K(k) 是用數值數學中最美的算法之一來計算的:算術-幾何平均。取兩個數,反覆把它們替換成它們的普通平均和幾何平均(兩數乘積的平方根),這兩個數列便奔向一個共同極限,正確位數大致每步翻倍。三到四次迭代就已給出機器精度。
- 認出形狀。每當一個積分裡含有一個三次或四次多項式的平方根——或一個你能整理成 1 over square root of (1 minus k^2 sin^2 theta) 的表達式——就要懷疑這是橢圓積分,而不是你技巧的失敗。
- 化歸為勒讓德標準形。透過標準代換,把積分推入 F、E 標準形之一,並讀出幅角 phi 和模 k。這才是真正的本領——做到這一步,活兒就幹完了。
- 交給現成的程式。呼叫內建的 K、E、F,或算術-幾何平均,收下十五位有效數字。一個命名的非初等函數,求值起來並不比 sin 或 log 更難。
反轉積分:橢圓函數登場
現在是那個深刻的轉折。回想一個你也許沒這樣表述過的 Volume I 事實:普通正弦是一個積分的反函數。若你寫下 u = integral from 0 to x of dt over square root of (1 minus t^2),這個積分不過是 arcsin x,於是 x = sin u。正弦正是反轉那個最簡單的弧長積分所得到的。那麼自然要問:若改為反轉那個橢圓弧長積分,又會得到什麼?
你得到的是一對全新的週期函數。令 u = F(phi, k) 再反轉:定義 sn(u, k) = sin phi、cn(u, k) = cos phi、dn(u, k) = square root of (1 minus k^2 sin^2 phi)。這就是雅可比橢圓函數——參見雅可比橢圓函數。它們是正弦餘弦更豐富的表親。它們滿足貌似熟悉的恆等式 sn^2 + cn^2 = 1 與 k^2 sn^2 + dn^2 = 1,以及乾淨的導數 d(sn)/du = cn times dn。當 k = 0 時它們恰好塌回 sn = sin、cn = cos、dn = 1;當 k = 1 時它們退化為雙曲函數 sn = tanh、cn = dn = sech。
標誌性的性質是雙週期性。普通正弦只有一個週期 2 pi,沿實軸行進。允許 sn 的自變量取複數,第二個獨立的週期就在虛方向上出現。於是該函數在鋪滿複平面的一整片平行四邊形格子上重複——它同時沿兩個方向週期,而這是任何非常數的初等函數都做不到的。這種格子結構最乾淨的化身是魏爾斯特拉斯 P 函數,它直接由格子構造,並滿足微分方程 (P')^2 = 4 P^3 minus g_2 P minus g_3——這正是一條橢圓曲線的方程,是從微積分這一隅通向數論與密碼學的橋樑。
於是迴路又閉合回鐘擺。非線性擺的精確角運動——不是小角度的贗品,而是真東西——可以寫成含 sn 的封閉形式 theta(t)。同樣這些函數描述自旋剛體、淺水的橢圓餘弦波與孤立波,以及非線性電路中的振盪。每當小角度線性化太過粗糙、你需要真正的非線性振盪時,sn、cn、dn 就是那個誠實的答案,是普通正弦餘弦根本給不了的。