一個改變函數構成的積分
你已經花了整整一級,學著去*求解*微分方程——猜測形式、建立特徵方程、變易參數。拉普拉斯變換給出的是一筆全然不同的交易。它不在微積分自己的語言裡去解微分方程,而是把整個問題馱過一座橋,送進一個新世界,在那裡方程根本不再是微分的了——它只是普通的代數。你在那邊解掉容易的代數,再把答案馱回來。這座橋就是一個積分,而學會走過它,便是本級的全部內容。
橋本身在這裡。給定一個對 t >= 0 有定義的函數 f(t)——把 t 想成時間,t = 0 就是你把系統打開的那一刻——它的拉普拉斯變換是 F(s) = 從 0 到無窮對 e^{-s t} f(t) dt 積分。這就是拉普拉斯變換的全部定義。你拿起你那個關於時間的函數,乘上衰減的指數 e^{-s t},再把乘積一路積分到無窮。出來的東西不再依賴 t——t 已被積掉。它只依賴那個新變量 s,於是我們把結果記作 F(s),並稱 s 為變換變量。
請留意這個運算的形狀。積分對 t 進行,給出一個關於 s 的答案——這恰恰是你在 伽馬函數把一個對 t 的積分變成它的參數的函數時見過的樣式。變換是一台機器,吃進一整個函數,吐出一整個函數:餵它 f(t),收回 F(s)。我們把它寫作 L{f(t)} = F(s),你應當想像兩個平行世界——函數活在時間裡的 t 域,以及它們的變換所棲居的 s 域。我們所做的一切,就是把一個問題從前者搬到後者,再搬回來。
親手算出最初幾個變換
定義是一個反常積分,那就讓我們當真算一個,感受它如何運作。取最簡單的函數,f(t) = 1 對一切 t >= 0。於是 F(s) = 從 0 到無窮對 e^{-s t} dt 積分。e^{-s t} 的原函數是 -(1/s) e^{-s t},而當 t 跑向無窮,e^{-s t} 衰減到零——*前提是 s 為正*,這樣指數才真的在衰減。上限處此項消失;在 t = 0 處它等於 -(1/s)。相減,F(s) = 0 - (-(1/s)) = 1/s。所以 L{1} = 1/s。一個時間上的常值函數,變成了 s 裡的一個簡單的倒數。
現在是最重要的那一對。取 f(t) = e^{a t},一個增長或衰減的指數。則 e^{-s t} f(t) = e^{-s t} e^{a t} = e^{-(s - a) t},這積分正是我們剛做過的那個,只是把 s 換成了 s - a。所以 L{e^{a t}} = 1/(s - a),這回當 s > a 即 s - a > 0 讓指數保持衰減時成立。盯著它看:指數函數,那個*就是*自己的導數(差一個常數)、在每一個常係數常微分方程裡都是主角的函數,竟變換成了死板簡單的代數對象 1/(s - a)。這一個事實——指數變成 s 裡的簡單極點——正是變換之所以能駕馭微分方程的那台秘密引擎。
L{e^{a t}} = integral_0^infinity e^{-s t} e^{a t} dt
= integral_0^infinity e^{-(s - a) t} dt
[ -1 ]^{t -> infinity}
= [ ------- e^{-(s-a)t} ]
[ s - a ]_{t = 0}
= 0 - ( -1/(s - a) ) ( needs s - a > 0 )
= 1/(s - a)
Set a = 0 -> L{1} = 1/s.
Replace a by i*omega and take parts -> L{cos(omega t)} = s/(s^2 + omega^2),
L{sin(omega t)} = omega/(s^2 + omega^2).積分被允許棲居之處:收斂性
你已經感到那個圈套了:每次計算都附帶一條小字條件——常數要 s > 0,指數要 s > a。這不是麻煩,而是定義中真實且必要的一部分。一個伸向無窮的反常積分,只有當被積函數縮得足夠快時才給出有限的數。因子 e^{-s t} 是我們的收縮劑,它能否取勝,取決於 f(t) 增長得有多快。使定義積分收斂的那個 s 的集合,稱為收斂域,在它之外,變換 F(s) 乾脆就不作為一個數而存在。
這個規律美妙地簡單。對一個指數階的函數——即在大 t 時增長不快於某個 C e^{a t} 的函數——收斂域是一個半平面,即一切實部 Re(s) > a 的 s。這個門檻 a 稱為收斂橫坐標;它恰好坐落在 f(t) 增長最快之處的右側。對 f(t) = e^{3 t},收斂域是 Re(s) > 3;對像 sin(t) 這樣的有界函數,是 Re(s) > 0;對多項式(增長慢於任何指數),又是 Re(s) > 0。想像 s 平面裡一條豎直的線向右掃過:變換棲居於這條「f 不再可馴服」之線右側的一切之上。
變換對表,以及它為何就足夠了
實踐中你幾乎從不去算那個定義積分。一旦知道了少數幾個基本變換,其餘的一切都由它們拼裝出來,因為變換直接從積分繼承了線性:L{c_1 f + c_2 g} = c_1 L{f} + c_2 L{g}。於是你備一張簡短的變換對表,兩個方向都去查它。核心條目正是我們上面掙來的那些及其近親:t 的冪、指數,以及把一個虛指數餵進 L{e^{a t}} 而得到的那對三角函數。
下面是那些主力,各自待在自己的半平面上。L{1} = 1/s,更一般地,對整數 n 有 L{t^n} = n!/s^{n+1}——而對非整數冪,階乘換成伽馬函數,L{t^p} = Gamma(p+1)/s^{p+1},正是本卷兩條線索相遇的一處妙地。再有 L{e^{a t}} = 1/(s - a);L{cos(omega t)} = s/(s^2 + omega^2);以及 L{sin(omega t)} = omega/(s^2 + omega^2)。一個增長或衰減的振盪,e^{a t} 乘一個正弦或餘弦,只不過把餘弦那對平移一下:L{e^{a t} cos(omega t)} = (s - a)/((s - a)^2 + omega^2)。最後這一步——在時間裡乘 e^{a t} 變成把 s 換成 s - a——就是第一平移定理,整張表裡用得最多的那個捷徑。
為什麼這把微積分變成了代數
現在是那個為整座橋正名的回報。把變換作用到一個導數上,分部積分一次:L{f'(t)} = s F(s) - f(0)。仔細看發生了什麼——*t 世界裡的求導,變成了 s 世界裡乘以 s*,而初值 f(0) 從邊界項裡掉出來,分文不取地直接烤了進去。這就是導數的變換,用兩次便給出 L{f''(t)} = s^2 F(s) - s f(0) - f'(0)。每一個你本得用微積分去追的導數,都變成了一個你能用代數去追的 s 的冪。
於是想像一個線性常微分方程,比如 y'' + 3 y' + 2 y = f(t),初值已知。把每一項都變換。每個導數都化成 s 的冪乘 Y(s) 再減去初始數據;方程變成 (s^2 + 3 s + 2) Y(s) = (由 f 和初值搭成的某個東西)。任何地方都不再有導數——只有 s 與 Y(s),被一個普通的代數方程拴在一起。你就像中學時解 x 那樣去解它:相除。Y(s) = (右端)/(s^2 + 3 s + 2)。微分方程溶解成了一個分式。
- 變換。把 L 作用到微分方程兩邊;每個導數都變成 s 的冪,初始條件經由導數法則自動進入。
- 解代數。結果是一個關於 Y(s)、再無導數的普通方程——像解任何未知數那樣,通過相除把 Y(s) 孤立出來。
- 求逆。用部分分式把 Y(s) 揉成與表匹配的若干塊,再把每一塊讀回成 t 的函數——那個 f(t) 就是你的答案,初始條件已經滿足。
這三步節奏——變換、解代數、求逆——是前路一切的脊梁。後面幾篇指南為每一步添上血肉,並拓寬它的覆蓋:如何變換開關與驟然的衝擊(海維賽德階躍與狄拉克 delta),好讓連不連續的驅動都乾淨地通過;卷積定理如何處置任意的輸入;以及求逆在複 s 平面裡究竟如何作為一個圍道積分運作。但方法的內核已經握在你手中:一個積分架起通往某個世界的橋,在那裡求導不過是乘以 s,而一張小表載著你來回跨過它。