這些函數從何而來
上一篇裡,伽瑪函數為階乘在每一個實數與複數自變量上都安了家。這正是本級階梯的統一套路:物理或統計中的某個問題反覆寫出同一個定積分,而這個積分沒有初等閉形式,於是我們給它起個名字,一勞永逸地研究它。本篇的四個函數正是如此——它們是一批你會反覆遇到的積分的、有名字、可製表、可畫圖的答案。
貝塔函數:伽瑪函數的雙臂兄弟
貝塔函數由一個在單位區間上的積分定義:B(x, y) = 從 0 到 1 對 t^{x-1} (1-t)^{y-1} dt 的積分,在 x > 0 且 y > 0 時有效。可以把它想成兩個冪之間的一次加權握手——一個因子偏向左端點 t = 0,另一個偏向右端點 t = 1。它是對稱的,B(x, y) = B(y, x),原因很簡單:把 t 換成 1 - t 就把兩個端點對調了。
最醒目的事實是:貝塔其實根本不是一個新函數:B(x, y) = Gamma(x) Gamma(y) / Gamma(x+y)。於是一張伽瑪函數表就為每一個貝塔值定了價。這個恆等式正是貝塔無處不在的原因——只要出現形如 t^{a}(1-t)^{b} 的積分(在機率論和組合學裡它到處都是),你就能把答案讀作一組伽瑪值之比。做個驗算:B(1, 1) = 從 0 到 1 對 1 dt 的積分 = 1,而 Gamma(1)Gamma(1)/Gamma(2) = 1·1/1 = 1,吻合。
誤差函數:馴服鐘形曲線
在第一卷你已認識高斯積分:從負無窮到正無窮對 e^{-x^2} dx 的積分 = sqrt(pi),從 0 到無窮對 e^{-x^2} dx 的積分 = sqrt(pi)/2。這個乾淨的總值藏著一樁煩惱:部分面積,即從 0 到 z 的積分,沒有初等原函數。對於『鐘形曲線在 0 與 z 之間佔了多少面積』,並不存在用指數和冪寫出的公式。於是我們給它起名。誤差函數為 erf(z) = (2/sqrt(pi)) · 從 0 到 z 對 e^{-t^2} dt 的積分,常數 2/sqrt(pi) 是這樣選的:使 erf(無窮) = 1。
想象它的圖象:erf 是一條奇對稱的 S 形曲線,erf(0) = 0,先迅速上升,再逐漸壓平、貼向直線 y = 1(左側則貼向 y = -1)。它的搭檔,互補誤差函數 erfc(z) = 1 - erf(z),度量 erf 留下的那條尾巴——z 之外的面積,也就是一個正態樣本落在遠離均值處的概率。統計學家就住在這條尾巴裡,所以 erfc 也配得上自己的名字與自己的按鈕。
我們究竟怎樣取出數值?在零附近,把被積函數展成泰勒級數——e^{-t^2} = 1 - t^2 + t^4/2 - ...——再逐項積分,它處處收斂,給出 erf(z) = (2/sqrt(pi))(z - z^3/3 + z^5/10 - ...)。對大的 z 這條級數沒用(需要巨量相消),此時改用尾部的漸近展開:erfc(z) ~ e^{-z^2}/(z sqrt(pi)) · (1 - 1/(2z^2) + ...)。這條漸近級數若把所有項都保留其實是發散的——可一旦在大 z 處只取前幾項截斷,它精確得令人咋舌。這種『小 z / 大 z』的二分法,是特殊函數計算反覆出現的節奏。
指數積分、正弦積分與餘弦積分
還有三個『弄丟了閉形式』的積分函數補全這一家族,每當一個看似無害的被積函數在分母上帶個 1/t 或 t,你就會遇到它們。指數積分為 Ei(x) = 從負無窮到 x 對 e^{t}/t dt 的積分(在 t = 0 處的極點取主值),以及與之密切相關的 E_1(x) = 從 x 到無窮對 e^{-t}/t dt 的積分。正弦積分與餘弦積分為 Si(x) = 從 0 到 x 對 sin(t)/t dt 的積分,以及 Ci(x) = 負的、從 x 到無窮對 cos(t)/t dt 的積分。
Si 最容易畫出來。它的被積函數 sin(t)/t 就是那著名的鼓包——在 t = 0 處等於 1,再以越來越小的波紋盪下去。當 x 向外延伸,Si(x) 上升、過衝,隨後以衰減振盪安頓到極限 Si(無窮) = pi/2——這個數恰好就是狄利克雷積分,從 0 到無窮對 sin(t)/t dt 的積分 = pi/2。Si 的過衝不是數值故障;它正是本卷稍後會成為傅立葉級數跳躍處吉布斯現象的那種振鈴。相比之下,Ci 在 x 趨於 0 時像 ln(x) 一樣跌向負無窮,再在大 x 處向 0 振盪。
它們都在哪兒出現?Ei 和 E_1 支配輻射輸運、線源散熱,以及那著名的、估計某數以下有多少素數的對數積分。Si 和 Ci 刻畫狹縫邊緣的繞射圖樣,以及濾波器和天線的響應。它們一點也不冷僻——它們正是物理與工程裡那種日常的『本該是閉形式卻偏偏不是』,這恰恰就是它們被命名的原因。
它們如何相連——以及如何使用
這些函數並非五件互不相干的奇珍;它們是同一棵樹上的葉子。貝塔是若干伽瑪函數的閉形式組合。誤差函數則是上限可移動的伽瑪——事實上 erf 正是同一個高斯積分的不完全版本,而那個積分一旦補全,就給出了上一篇的 Gamma(1/2) = sqrt(pi)。而 erf、Ei、Si、Ci 用的全是同一招:取一個原函數逃出初等函數範圍的被積函數,把流動的面積定義成一個新函數,再去研究它。認出這一招,才是本級真正要練的本領。
- 認出標誌性被積函數:[0,1] 上的 t^{a}(1-t)^{b} 意味著貝塔;e^{-t^2} 意味著誤差函數;e^{t}/t 意味著指數積分;sin(t)/t 或 cos(t)/t 意味著正弦或餘弦積分。
- 把你的積分揉成標準形:用換元做縮放與平移,讓積分限和被積函數與定義嚴絲合縫——大部分功夫不過是常數的記帳。
- 讀出數值:對貝塔,轉成伽瑪值;其餘則查表或呼叫 erf、erfc、Ei、Si、Ci,在 0 附近用小自變量級數、在大自變量處用漸近尾。
書架上緊鄰的是橢圓積分——下一篇的主題——它由橢圓的弧長和真實單擺的擺動逼出來。教訓一字不變地延續:當真實問題變硬時,答案往往不是『積不出來』的失敗,而是『還沒人告訴你它名字』的一個函數的登場。