積分住在哪裡:一條垂直線
你早在學複分析之前就見過拉普拉斯變換:它把一個定義在 t >= 0 的時間函數 f(t),送成新變數 s 的函數 F(s) = integral from 0 to infinity of e^(-s t) f(t) dt。它的用處在於槓桿。對時間微分變成乘以 s;一道彆扭的、關於 f(t) 的微分方程坍縮成關於 F(s) 的普通代數,你一行就能解開。麻煩在回程。一旦握著 F(s),你要怎麼還原出它從何而來的 f(t)?這個反演步驟,正是複平面終於非請出不可之處,也是本篇的全部主題。
驚喜在這裡。變換 F(s) 不只是一條實數公式——餵給它一個複數 s,那積分收斂,並把 s 定義成一個貨真價實的 亞純函數,在積分有意義之處全純、在別處綴著一些孤立奇點。積分在半平面 Re s > c 收斂,其中 c 是收斂橫坐標,剛好坐落在 f 最劇烈增長的右側。所以 F(s) 是一個你可以在平面上搬動的解析物件,而反演公式正反映了這點:它沿著垂直線 Re s = gamma 把 F(s) e^(s t) 往上積,其中 gamma 是任何大於 c 的實數——一條安全地畫在所有奇點右側的線。
F(s) = integral from 0 to infinity of e^(-s t) f(t) dt (the transform)
1 gamma + i*infinity
f(t) = ------- integral F(s) e^(s t) ds (the inverse)
2 pi i gamma - i*infinity
the path is the vertical line Re s = gamma, with gamma > c布羅米奇圍道:把直線合成迴路
一條從 gamma - i*infinity 往上伸到 gamma + i*infinity 的垂直線並不是閉迴路,所以留數定理碰不到它——然而一條光禿禿的無窮直線積分,恰恰是那套機器最愛吞下的對象。補救之道,正是替整套方法命名的那一步。加上一段大圓弧,把直線合成迴路——對 t > 0,是一條向左鼓出、伸進半平面 Re s < gamma 的大半圓——使直線段與圓弧合起來圍住一塊區域。那條閉合路徑就是布羅米奇圍道,而反演問題終於變成一條你讀得出來的圍道積分。
為什麼向左鼓、又只在 t > 0 這麼做?看看你一路拖著的因子 e^(s t)。在圓弧上寫 s = gamma + R e^(i theta);那麼 e^(s t) 帶著大小 e^(t * Re s),而在左半邊(Re s < gamma),隨著 Re s 下降,這指數被往下壓。對正的 t,把圓弧往左推會使 e^(s t) 衰減,遠方圓弧的貢獻便可弄成消失。對 t < 0 符號翻轉:你反而必須向右閉合,那裡 F 解析,你收集到的每個留數都是零,答案是 f(t) = 0——這恰恰正確,因為 f 本來就只定義在 t >= 0。這幾何正在強制因果律。
讓圓弧消失:若爾當引理
人很想揮揮手說 1/s 在縮小、把圓弧打發掉,但光這樣太弱:圓弧長度如 pi R 般增長,把 1/R 的衰減和 R 大小的路徑粗略相配,得到一個不肯死去的常數。誠實的工具是若爾當引理,它利用 e^(s t) 裡的振盪,而不只是它的大小。在向左閉合的圓弧上,大小 e^(t * Re s) 不只是在某處小——它離實軸越遠就衰減得越陡,以致你一旦在整段圓弧上積分,答案就被真切地擠到零。引理只要求當 |s| -> infinity 沿圓弧時 F(s) -> 0 一致成立,這是個溫和條件,多數體面函數的變換都滿足。
有一條值得明說的附帶條款,因為它常絆倒人。這個形式的若爾當引理嚴格要求 t > 0。在 t = 0,e^(s t) 因子就只是 1,提供不了衰減,乾淨的論證崩潰;那裡的反演在跳躍處回傳 f 左右極限的平均,而非單一值——正是你或許記得的、傅立葉級數在不連續處的中點行為。又若 F 只如 1/s 般衰減、再無餘裕,你就正好站在引理允許的邊界上:你保住了 e^(s t) 的振盪,卻負擔不起再丟掉一個 1/s,所以這筆帳必須小心算,不能反射性地揮手。
從極點讀出答案
圓弧一旦消失,畫面美得簡單。直線積分等於 2 pi i 乘以留數之和,前頭的 1/(2 pi i) 恰好把那因子消掉,剩下的是 f(t) = F(s) e^(s t) 在直線左側每個奇點處的留數之和。在 F 是多項式之比的常見情形,那些奇點就只是極點,每一個都貢獻時間信號的一塊乾淨碎片。s = a 處的單極點貢獻一個正比於 e^(a t) 的項;a 的實部定下增長或衰減速率,虛部定下振盪頻率。F 的極點,字面上就是 f(t) 的建構積木。
- 定位 F(s) e^(s t) 的每個奇點——對有理的 F,這些就是分母的零點,即極點。
- 確認它們全落在直線 Re s = gamma 的左側,使向左的布羅米奇圍道把它們圍住。
- 確認在大圓弧上 F(s) -> 0,使若爾當引理在 t > 0 時消掉圓弧。
- 在每個被圍的極點處算 F(s) e^(s t) 的留數並相加;其和即 f(t)。
拿最溫和的例子看它運作。取 F(s) = 1/(s - a),是單一指數的變換。它唯一的奇點是 s = a 處的單極點,而對任何 gamma > Re a,布羅米奇圍道都把它圍住。F(s) e^(s t) = e^(s t)/(s - a) 在 s = a 的留數,依單極點法則,就只是 e^(a t)。把這一個留數加起來得 f(t) = e^(a t)——而 e^(a t) 的拉普拉斯變換確實是 1/(s - a)。機器原封不動地還原了我們出發的函數,極點位置 a 直接化身為指數速率。
當畫面不全是極點時
極點的整潔故事是這門學問友善的那一半。許多誠實的變換改帶一個分支點——譬如 F(s) = 1/sqrt(s),其反演是 1/sqrt(pi t);或源自熱方程、含 sqrt(s) 的變換。在分支點你不能光把留數相加,因為 F 在那裡甚至不是單值的;你必須從該點拉出一條分支切割、把布羅米奇圍道繞著它變形,將反演化成一條像鑰匙孔般纏住切割的分支切割積分。同樣的「閉合再讀」哲學依舊成立,只是這筆帳更豐富,答案往往以一條實積分、而非有限多個指數之和的形式現身。
在你反射性地伸手去抓留數之前,兩個誠實的警告。其一,若 F 有無窮多個極點——分佈式系統的轉移函數正是如此——你會收集到一條無窮的指數項級數,你便欠下一個收斂論證,而非一個一廂情願的和。其二,本性奇點確實超出簡單留數配方的範圍;那裡局部行為狂野,你必須親手展開。這些都不會破壞布羅米奇方法——它仍是閉合圍道、論證圓弧、加總局部貢獻——只是提醒你:留數定理是一件精密儀器,不是一句口號,它的每一條假設都得付出代價。