從封閉的圓到開放的線
上一篇指南裡,你藉著把實變數繞上單位圓,馴服了在整個週期上的積分:代入 z = e^(i theta),把從 0 到 2 pi 的三角積分變成一條封閉迴路上的圍道積分,正好可餵給留數定理。那行得通,是因為路徑本來就封閉——角度 theta 繞回了起點。這篇指南要對付更棘手的形狀:在整條實軸上、從負無窮到正無窮的積分。實軸是一條開放而無限的道路,不是迴路,所以沒有現成的封閉圍道可以餵給留數定理。整篇指南的核心想法,就是如何把它封閉起來。
畫面是這樣的。取一個很大的半徑 R,沿著實軸從 -R 走到 +R。接著別停下,改沿一條半徑為 R、向上鼓進上半平面的巨大半圓弧繞回起點。那段直線加上那條圓弧合起來構成一條封閉迴路——一個躺在平底上的大寫字母 D。如今留數定理便能用在這條封閉迴路上:它的值等於 2 pi i 乘以被困在內部的那些極點上的留數總和,也就是落在上半平面的那些被積函數的極點。
但請注意,我們真正想要的只有那段直線部分——沿實軸的積分。圓弧是我們黏上去封閉迴路的額外鷹架。於是這個策略天生背著一筆債:我們會用留數算出整條迴路的值,然後必須證明當 R 增大到無窮時,圓弧的貢獻會縮到零。若圓弧消失,封閉迴路的值就塌縮成實軸上的積分,我們便能讀出答案。若圓弧不消失,這套方法就失效,我們什麼也沒學到。這份誠實的帳目——留數減去一條消失的圓弧——正是每一次瑕積分計算的心跳。
有理函數:圓弧自己死去之時
最乾淨的情形,是有理函數 P(x)/Q(x) 的積分,其中 P 與 Q 為多項式,且 Q 在實軸上沒有零點。我們要的是這個有理函數在整條線上的瑕積分。如上構造 D 形圍道。圓弧的貢獻由 ML 不等式掌控:圓弧上的積分至多等於它的弧長 pi R,乘以被積函數在圓弧上的最大值。在半徑為 R 的圓上,|P(z)/Q(z)| 的大小約莫像 R^(deg P) / R^(deg Q)。所以圓弧的貢獻大致被 pi R 乘以 R^(deg P - deg Q) 所限制住。
細讀這個估計。這個界的數量級是 R^(1 + deg P - deg Q),而它在 R 趨於無窮時死去,恰好就在 1 + deg P - deg Q 為負數之時——也就是 deg Q 至少比 deg P 大 2 之時。分母的次數必須至少高過分子兩級。這正是讓那個實積分一開始就收斂的同一道差距:若 Q 只比 P 高一級,被積函數的衰減就像 1/x,其積分以對數方式發散,於是根本沒有有限的答案可找。「高兩級」這條規則不是技術細節,而是收斂性與圓弧消失用同一個聲音說出的同一件事。
- 檢查收斂:確認 deg Q >= deg P + 2,且 Q 在實軸上沒有零點,使被積函數在整條軸上都有限,積分收斂。
- 找出 P(z)/Q(z) 在上半平面的極點——即虛部為正的 Q 的零點。只有這些會被 D 形圍道圍住。
- 在每個這樣的極點上計算留數。對單極點 z_0,單極點公式給出 Res = P(z_0)/Q'(z_0),也就是商的捷徑。
- 把這些留數相加,再乘以 2 pi i。由於圓弧消失,這個值就直接等於整條實軸上的積分。
一個你能記在腦中的範例
取 1/(x^2 + 1) 在整條實軸上的積分。這裡 deg Q = 2、deg P = 0,所以 deg Q - deg P = 2,圓弧將會消失。被積函數 1/(z^2 + 1) = 1/((z - i)(z + i)) 在 z = i 與 z = -i 處有單極點。只有 z = i 落在上半平面,所以只有它的留數要算。用商的捷徑,z = i 處的留數是 1/(2z) 在 z = i 的值,也就是 1/(2i)。
f(z) = 1 / (z^2 + 1) = 1 / ((z - i)(z + i)) upper-half-plane pole: z = i (simple) Res(f; i) = 1/(2z) |_(z=i) = 1/(2i) integral over R = 2 pi i * (1/(2i)) = pi check: integral of 1/(x^2+1) dx = arctan(x), pi/2 - (-pi/2) = pi OK
相乘後,整條線上的積分是 2 pi i 乘以 1/(2i) = pi。而這一個你可以徒手驗算:1/(x^2 + 1) 的原函數是 arctan(x),它從 -pi/2 跑到 pi/2,得到 pi。圍道的機器與初等微積分一致,這正是你在信賴它去算那些完全沒有初等原函數的積分之前,所需要的信心基礎。這才是真正的獎賞:同一個 D 形圍道能算出像 1/(x^4 + 1) 的積分這樣的怪物,那裡用部分分式與 arctan 會是場噩夢,但留數卻只是一段短短的計算。
振盪的被積函數,以及為何天真的界會失效
現在把難度提高。像 cos(x)/(x^2 + 1) 或 sin(x)/x 在實軸上的積分,到處出現於物理與訊號處理——它們是傅立葉積分。自然的做法是把 cos(x) 或 sin(x) 換成 e^(i x),積分後再取實部或虛部,因為 cos x = Re(e^(i x))、sin x = Im(e^(i x))。在圍道上我們於是面對 e^(i z) = e^(i x - y),其模為 e^(-y)。在上半平面 y 為正,所以 e^(i z) 在那裡是指數地小——這對圓弧是天大的好消息。
但這裡有個陷阱,值得你敏銳地感受一下。假設那個有理因子只像 1/z 那樣衰減——例如 sin(x)/x,其有理部分為 1/z。天真的 ML 界把弧長 pi R 乘以被積函數的最大值。指數 e^(-y) 在圓弧頂端附近極小,但在 y = 0 的兩端恰等於 1。所以 |e^(i z)/z| 在圓弧上的最大值約為 1/R,而 pi R 乘以 1/R 等於 pi——一個不趨於零的常數。粗糙的 ML 估計在此無能為力。我們需要一個更細緻的論證,把指數的衰減在整條圓弧上認真對待,而不只是在它最糟的那單獨一點。
約當引理:精準的救援
在振盪情形裡賺到圓弧消失的工具,是約當引理。它說:若 g(z) 在上半圓上隨 R 增大而一致地趨於零——哪怕只像 1/R 那樣——那麼對任何正常數 a,g(z) e^(i a z) 在上半圓弧上的積分都會在 R 趨於無窮時趨於零。換句話說,約當引理讓指數 e^(i a z) 去做有理衰減本身做不到的工作:即使 g 只有微弱的 1/R 衰減,乘上那個又振盪又衰減的指數,也足以殺死圓弧。
證明內部的祕密,是關於正弦的一條誠實不等式。在圓弧 z = R e^(i theta) 上,指數的大小是 e^(-a R sin theta),而整條圓弧的貢獻牽涉到 e^(-a R sin theta) 對 theta 從 0 到 pi 的積分。關鍵估計是:在第一個四分之一上 sin theta >= 2 theta / pi——弦落在正弦曲線之下——這馴服了那個積分,逼它像 1/(a R) 那樣下降。於是指數的衰減,誠實地在整條圓弧上加總、而非在它最糟一點上瞎猜,便勝過了弧長那緩慢的 1/R 增長。那單獨一條正弦不等式就是整部引擎,其餘一切都只是記帳。
於是振盪情形的食譜,就是有理情形的食譜加上這道保險。對 cos(x)/(x^2 + 1) 的積分:把它寫成 e^(i x)/(x^2 + 1) 積分的實部;往上封閉;唯一的上半平面極點是 z = i;e^(i z)/(z^2 + 1) 在 z = i 的留數是 e^(-1)/(2i);乘以 2 pi i 得到 pi/e;而約當引理保證圓弧已掉落。取實部(此處即全部,因答案本就是實的)得到 cos(x)/(x^2 + 1) 的積分 = pi/e。那個被積函數沒有初等原函數——然而留數加上約當,交出了一個乾淨的封閉形式。
這留下了什麼尚未處理
還剩下兩道裂縫,而本階梯的最後一篇指南正是為填補它們而設。第一,我們假設了 Q 在實軸上沒有零點——被積函數在軸上有限。但著名的 sin(x)/x 積分帶有因子 1/x,有個極點正正坐在圍道上的 x = 0。你無法讓圍道穿過極點,於是你用一個小半圓繞過它、讀出一個主值,並從這道繞行中拾取半個留數。第二,被積函數若帶有 x 的非整數次冪或對數,便是多值的,需要在一條尊重分支切割的圍道上做留數計算。
貫穿這一切的主線,是一個有紀律的習慣:挑一條圍道,它的直線部分正是你要的積分,它的彎曲部分你能證明會消失(靠 ML 或靠約當),而它圍住的極點你能用本階梯開頭的留數公式算出。留數定理供給數值;你對圓弧的誠實供給使用它的權利。精通這組搭配,整座定積分的動物園——有理的、振盪的、分支切割的、主值的,乃至用留數加總的無窮級數——便都向一個單一、可重複的方法敞開。