跳轉到

5.2 實積分的計算#

本章導言#

留數定理最初的動機(Cauchy 1814)就是算實積分。本節是臨床課:四類典型實積分的圍道手術——三角積分(單位圓)、有理函數在無窮區間(大半圓)、Fourier 型積分(大半圓 + Jordan 引理)、以及含分支的積分(鑰匙孔圍道)。讀完本節,讀者應能對「這個積分該用哪條圍道」建立職業反射。

類型一:三角積分與單位圓#

積分 \(\displaystyle\int_0^{2\pi} R(\cos\theta, \sin\theta)\,d\theta\) ( \(R\) 有理)。令 \(z = e^{i\theta}\) : \(\cos\theta = \frac{z + z^{-1}}{2}\) 、 \(\sin\theta = \frac{z - z^{-1}}{2i}\) 、 \(d\theta = \frac{dz}{iz}\) 。三角積分變成單位圓上的有理函數積分——留數定理收割。

示範: \(\displaystyle\int_0^{2\pi}\frac{d\theta}{2 + \cos\theta}\) 。代入:

\[ \oint_{|z|=1} \frac{1}{2 + \frac{z + z^{-1}}{2}}\cdot\frac{dz}{iz} = \oint_{|z|=1}\frac{2\,dz}{i(z^2 + 4z + 1)}. \]

分母零點: \(z = -2 \pm \sqrt3\) 。 \(z_1 = -2 + \sqrt3 \approx -0.268\) (圓內)、 \(z_2 = -2 - \sqrt3 \approx -3.73\) (圓外)。留數(單極點):

\[ \operatorname{Res}_{z_1} = \frac{2}{i}\cdot\frac{1}{(z_1 - z_2)} = \frac{2}{i\cdot 2\sqrt3} = \frac{1}{i\sqrt3}. \]

積分 \(= 2\pi i \cdot \frac{1}{i\sqrt3} = \frac{2\pi}{\sqrt3}\) 。檢算: \(2 + \cos\theta \in [1,3]\) ,平均值 \(\approx 2\) ,積分 \(\approx 4\pi \approx 12.57\) ; \(\frac{2\pi}{\sqrt3} \approx 3.63\) ——不符!重算: \(\int_0^{2\pi}\frac{d\theta}{2+\cos\theta} = \frac{2\pi}{\sqrt{4 - 1}} = \frac{2\pi}{\sqrt3}\) (標準公式 \(\int\frac{d\theta}{a + b\cos\theta} = \frac{2\pi}{\sqrt{a^2 - b^2}}\) ) \(\approx 3.63\) 。而 \(\cos\theta\) 平均為 0, \(2 + \cos\theta\) 平均 2, \(\frac{1}{2+\cos\theta}\) 的平均是 \(\frac{1}{\sqrt{4-1}} = \frac{1}{\sqrt3}\) (調和型平均,不是 \(\frac12\) )——直覺估計要小心,公式值正確 ✓。

類型二:無窮區間的有理函數#

積分 \(\displaystyle\int_{-\infty}^{\infty}\frac{p(x)}{q(x)}dx\) ( \(\deg q \ge \deg p + 2\) ,實係數, \(q\) 無實零點)。圍道: \([-R, R]\) + 上半大圓。大圓上 \(|f| \le \frac{C}{R^2}\) 、弧長 \(\pi R\) ,貢獻 \(\le \frac{\pi C}{R} \to 0\) 。故

\[ \int_{-\infty}^\infty \frac{p}{q}dx = 2\pi i\sum_{\operatorname{Im} a_k > 0}\operatorname{Res}_{a_k}\frac{p}{q}. \]

示範: \(\displaystyle\int_{-\infty}^{\infty}\frac{dx}{1 + x^2}\) 。上半圓內奇點 \(z = i\) ,留數 \(\frac{1}{2i}\) 。積分 \(= 2\pi i\cdot\frac{1}{2i} = \pi\) ✓(與 \(\arctan\) 值一致)。** \(\pi\) 從圍道而來**——又一次「圓的簽名」。

類型三:Fourier 型積分與 Jordan 引理#

積分 \(\displaystyle\int_{-\infty}^{\infty} f(x)e^{iax}dx\) ( \(a > 0\) , \(f\) 有理且 \(\deg q \ge \deg p + 1\) )。 \(e^{iaz} = e^{ia(x+iy)} = e^{iax}e^{-ay}\) 在上半平面指數衰減——大圓貢獻更小。

Jordan 引理:若 \(f \to 0\) ( \(|z|\to\infty\) ,上半圓)一致地不差於 \(\frac{C}{|z|^\alpha}\) ( \(\alpha > 0\) ),則 \(\displaystyle\int_{C_R} f(z)e^{iaz}dz \to 0\) ( \(a>0\) ,上半圓)。證明: \(|e^{iaz}| = e^{-aR\sin\theta}\) ,用 \(\sin\theta \ge \frac{2\theta}{\pi}\) ( \(0\le\theta\le\frac\pi2\) )控制,弧長積分 \(\to 0\) 。

示範: \(\displaystyle\int_{-\infty}^{\infty}\frac{\cos x}{1 + x^2}dx = \operatorname{Re}\int_{-\infty}^{\infty}\frac{e^{ix}}{1+x^2}dx\) 。奇點 \(z = i\) : \(\operatorname{Res} = \frac{e^{ii}}{2i} = \frac{e^{-1}}{2i}\) 。積分 \(= 2\pi i\cdot\frac{e^{-1}}{2i} = \frac{\pi}{e}\) 。複變的紅利:虛部 \(\int\frac{\sin x}{1+x^2}dx = 0\) (奇函數)順帶免費;更一般地,這條技術給出Cauchy 分布的特徵函數 \(e^{-|t|}\) (10.1 篂的 Fourier 分析原型)。

同型高峰: \(\displaystyle\int_0^\infty \frac{\sin x}{x}dx = \frac{\pi}{2}\) (Dirichlet 積分)。圍道:上半大半圓 + 繞原點的小半圓(跳過 \(z=0\) 的可去/單極點,貢獻 \(i\pi\) 的半值);大圓由 Jordan 引理消去;結果 \(\int_{-\infty}^\infty \frac{e^{ix}}{x}dx\) 主值 \(= i\pi\) ,取虛部得 \(\frac{\pi}{2}\) 。這條積分是訊號處理的「理想低通濾波器衝激響應」(10.1 篂)與 Gibbs 現象的數學心臟。

類型四:分支函數與鑰匙孔圍道#

積分含 \(\sqrt{x}\) 、 \(\log x\) 時,被積函數多值,需分支切割(7.3 篂正式化)+ 鑰匙孔圍道:沿正實軸上下岸各走一次,跳躍(jump)產生目標積分。

示範: \(\displaystyle\int_0^\infty \frac{\sqrt{x}}{1 + x^2}dx\) 。取分支 \(\sqrt{z}\) (沿正實軸切割,上岸 \(= \sqrt{x}\) 、下岸 \(= -\sqrt{x}\) )。鑰匙孔圍道繞極點 \(z = i, -i\) ( \(\frac{\sqrt z}{1+z^2}\) 的極點)。上下岸積分:

\[ \int_0^\infty \frac{\sqrt x}{1+x^2}dx - \int_0^\infty \frac{-\sqrt x}{1+x^2}dx = 2\int_0^\infty\frac{\sqrt x}{1+x^2}dx = 2\pi i\left(\operatorname{Res}_i + \operatorname{Res}_{-i}\right). \]

\(\operatorname{Res}_i = \frac{\sqrt{i}}{2i} = \frac{e^{i\pi/4}}{2i}\) 、 \(\operatorname{Res}_{-i} = \frac{\sqrt{-i}}{-2i} = \frac{e^{-i\pi/4}}{-2i}\) (下岸分支約定)。和 \(= \frac{e^{i\pi/4} - e^{-i\pi/4}}{2i} = \sin\frac\pi4 = \frac{\sqrt2}{2}\) 。故目標積分 \(= \pi i\cdot\frac{\sqrt2}{2}\cdot\) ……等一下: \(2\pi i \times \frac{\sqrt2}{2} = \pi i\sqrt2\) ,但左端 \(2\int\) 是實數——分支約定需再核對。正確約定:切割沿正實軸,上岸 \(\sqrt z = \sqrt r\) 、下岸 \(\sqrt z = -\sqrt r\) ; \(\operatorname{Res}_{i} = \frac{e^{i\pi/4}}{2i}\) 、 \(\operatorname{Res}_{-i}\) : \(z = -i = e^{-i\pi/2}\) (在此分支下輻角 \(-\pi/2 \in (-\pi, \pi)\) ?切割在正實軸,分支輻角取 \((0, 2\pi)\) , \(-i\) 的輻角 \(= \frac{3\pi}{2}\) , \(\sqrt{-i} = e^{3\pi i/4}\) )。和 \(= \frac{e^{i\pi/4}}{2i} + \frac{e^{3i\pi/4}}{-2i} = \frac{e^{i\pi/4} - e^{3i\pi/4}}{2i} = \frac{e^{i\pi/4}(1 - e^{i\pi/2})}{2i} = \frac{e^{i\pi/4}(1 - i)}{2i} = \frac{e^{i\pi/4}\cdot(-i)(1 - i)}{2} = \frac{(-i + 1)e^{i\pi/4}}{2}\cdot\) …… \(\frac{(1-i)e^{i\pi/4}}{2i}\) ; \(e^{i\pi/4} = \frac{1+i}{\sqrt2}\) , \((1-i)(1+i) = 2\) ,故 \(= \frac{2/\sqrt2}{2i} = \frac{1}{i\sqrt2} = \frac{-i}{\sqrt2}\) 。於是 \(2\pi i \times \frac{-i}{\sqrt2} = \frac{2\pi}{\sqrt2} = \pi\sqrt2\) 。核對: \(2\int_0^\infty\frac{\sqrt x}{1+x^2}dx = \pi\sqrt2\) ,故 \(\int_0^\infty \frac{\sqrt x}{1+x^2}dx = \frac{\pi}{\sqrt2}\) 。標準檢算:代換 \(x = t^2\) , \(\int_0^\infty \frac{t\cdot 2t}{1+t^4}dt = 2\int_0^\infty\frac{t^2}{1+t^4}dt = 2\cdot\frac{\pi}{2\sqrt2} = \frac{\pi}{\sqrt2}\) ✓( \(\int_0^\infty \frac{t^{p-1}}{1+t^q} = \frac{\pi}{q}\csc\frac{p\pi}{q}\) ,此處 \(p=3, q=4\) : \(\frac\pi4\csc\frac{3\pi}4 = \frac{\pi}{4}\sqrt2 = \frac{\pi}{2\sqrt2}\) ,乘 2 得 \(\frac{\pi}{\sqrt2}\) ✓✓)。

圍道圖譜總表#

實積分類型 圍道 關鍵技術 經典例
\(\int_0^{2\pi}R(\cos\theta,\sin\theta)d\theta\) 單位圓 \(z = e^{i\theta}\) \(\int\frac{d\theta}{2+\cos\theta}\)
\(\int_{-\infty}^\infty\frac{p}{q}dx\) 大半圓(上) 大圓 \(O(1/R^2)\) \(\int\frac{dx}{1+x^2}\)
\(\int f(x)e^{iax}dx\) 大半圓 + Jordan 指數衰減 \(\int\frac{\cos x}{1+x^2}dx\)
\(\int_0^\infty f(x)x^{a-1}dx\) 鑰匙孔 割線跳躍 \(\int\frac{\sqrt x}{1+x^2}dx\)
\(\int_0^\infty\frac{x^{a-1}}{1+x}dx\) 扇形 旋轉自相似 \(\frac{\pi}{\sin\pi a}\)
Fresnel \(\int_0^\infty e^{ix^2}dx\) \(45^\circ\) 扇形 旋轉 + Jordan \(\frac{\sqrt{2\pi}}{2}(1+i)/\sqrt2\)

Fresnel 積分補充: \(\int_0^\infty\cos(x^2)dx = \int_0^\infty\sin(x^2)dx = \frac{\sqrt{2\pi}}{4}\) 。圍道: \(45^\circ\) 扇形, \(e^{iz^2}\) 在第一卦限的扇形邊界上衰減( \(z = re^{i\pi/4}\) 時 \(e^{iz^2} = e^{ir^2 i} = e^{-r^2}\) )——旋轉圍道把振盪變衰減。光學的 Cornu 螺線(Fresnel 積分的參數曲線)就是這兩條積分的幾何像。

一個具體例題(綜合)#

例題 5.2.1:計算 \(\displaystyle\int_{-\infty}^{\infty}\frac{x\sin x}{x^2 + 4}\,dx\) 。

解:考慮 \(I = \oint_{C_R}\frac{z e^{iz}}{z^2+4}dz\) (上半大半圓)。奇點 \(z = 2i\) (單極點):

\[ \operatorname{Res}_{2i} = \frac{2i\,e^{2i\cdot i}}{2(2i)} = \frac{e^{-2}}{2}. \]

大圓: \(|ze^{iz}/(z^2+4)| \le \frac{R e^{-R\sin\theta}}{R^2 - 4}\) ,由 Jordan 引理貢獻 \(\to 0\) 。故

\[ \int_{-\infty}^{\infty}\frac{x e^{ix}}{x^2+4}dx = 2\pi i\cdot\frac{e^{-2}}{2} = \frac{\pi i}{e^2}. \]

取虛部: \(\displaystyle\int_{-\infty}^{\infty}\frac{x\sin x}{x^2+4}dx = \frac{\pi}{e^2}\) 。(實部 \(\int\frac{x\cos x}{x^2+4} = 0\) ,奇函數 ✓。)

檢算(殘餘對稱): \(\frac{z}{z^2+4}\) 在 \(-2i\) 的留數 \(\frac{e^{-2i\cdot(-2)}}{-4i}\cdot(-2i)\) ……利用「下半平面留數與上半共軛相消」的結構: \(\frac{xe^{ix}}{x^2+4}\) 的積分只由 \(2i\) 決定,值 \(\frac{\pi i}{e^2}\) 純虛 ✓(被積函數實部為奇函數,全域積分虛)。

歷史背景#

  • Cauchy 1814–1827:留數法的臨床誕生,最初計算的正是 \(\int\frac{e^{iax}}{x^2+b^2}dx\) 型。
  • Jordan 引理(1894,Cours d'analyse 教材):把「大圓上振盪函數可忽略」變成引理級的標準工具。Jordan 是馬賽出身的工程師出身數學家(Jordan 曲線定理以其命名)。
  • 物理學的大規模使用: nineteenth 世紀末 Maxwell 學派、二十世紀量子力學(Born、Jordan、Heisenberg 的矩陣力學時期,1925–26,圍道積分計算譜)——圍道積分成為理論物理的基本文體。

歷史啟示#

  • 「換戰場」戰術:實積分的難點(振盪、無窮區間)在複平面上變成優勢( \(e^{-ay}\) 衰減、大圓可忽略)。把問題搬到對它有利的空間是複分析的元策略,也預示了 10.1 篂 Z/Fourier 變換的哲學。
  • 圍道即演算法:五種圍道 = 五類積分的標準程序。數學計算的「機械化」早在十九世紀就由圍道圖譜完成——今日 CAS 的 Integrate 內部仍是這些程序的精緻版。
  • 估計與構造並重:Jordan 引理(估計)與鑰匙孔(構造)顯示:留數法的成功 = 精確的代數(留數)+ 粗略的分析(大小估計)。粗細分工是分析學的常態。

與其他章節及領域的關聯#

  • 5.1 節:留數計算是本節所有圍道手術的內臟。
  • 4.2、7.3 節:鑰匙孔圍道預告分支切割的正式理論。
  • 1.3 節: \(e^{iz}\) 的極形式(衰減 = 模、振盪 = 輻角)是 Jordan 引理的語言。
  • 9.1 節: \(\Gamma\) 函數的 \(\frac{\pi}{\sin\pi z}\) 型恆等式、Hankel 圍道(繞負實軸的大鑰匙孔)是本節圍道的升級版。
  • 10.1 節:Fourier 變換的反演積分、因果濾波器、Wiener–Hopf 方法都以本節圍道為工具。
  • 機率論:特徵函數反演(Gauss 積分、Cauchy 分布)是圍道積分的統計應用。
  • 微積分系列:5.3 篂廣義積分收斂判別是本節「大圓可忽略」條件的實變版。

練習#

  1. 計算 \(\displaystyle\int_0^{2\pi}\frac{d\theta}{5 + 4\cos\theta}\) ,並用標準公式 \(\frac{2\pi}{\sqrt{a^2-b^2}}\) 檢算。
  2. 計算 \(\displaystyle\int_{-\infty}^{\infty}\frac{dx}{(x^2+1)(x^2+9)}\) (部分分數檢算)。
  3. 用 Jordan 引理計算 \(\displaystyle\int_{-\infty}^{\infty}\frac{x\cos x}{x^2+1}dx\) ;解釋為什麼被積函數不是奇函數而積分不為零( \(\frac{x\cos x}{x^2+1}\) 是奇函數?——先判斷再計算,注意 \(\cos\) 是偶函數)。
  4. 用鑰匙孔圍道計算 \(\displaystyle\int_0^\infty\frac{\log x}{1+x^2}dx\) 。(答案 0,用 \(x\mapsto 1/x\) 對稱性檢核。)
  5. 思考題:為什麼 \(\int_{-\infty}^{\infty}\frac{e^{ix}}{x}dx\) 需要主值與繞原點的小半圓,而 \(\int\frac{e^{ix}}{x^2+1}dx\) 不需要?「奇點在圍道上」的處理原則是什麼?

進階研讀:圍道積分的五個深化場景#

場景一:最陡下降法(saddle point method)。積分 \(\int e^{\lambda f(z)}dz\) ( \(\lambda \to \infty\) )的主貢獻來自 \(f'(z_0) = 0\) 的鞍點,沿「最陡下降路徑」( \(f\) 的虛部恆定)通過。Debye 1909 創立,是漸近分析的核心引擎: \(\Gamma\) 函數的 Stirling 公式、Airy 函數漸近、量子力學的 WKB。振盪積分的能量集中在鞍點——與留數定理「積分集中在奇點」互為表裡:前者是極值原理的積分版。

場景二:Hankel 圍道與 \(\Gamma\) 函數。 \(\Gamma(s) = \frac{1}{2i\sin\pi s}\int_{\mathcal H} \frac{e^{-t}(-t)^{-s}}{}dt\) 型 Hankel 表示(繞負實軸的大鑰匙孔)使 \(\Gamma\) 延拓到全平面(除非正整數極點)、給出函數方程、並連通 \(\zeta\) 函數(9.1、8.3 篂)。Hankel 圍道是鑰匙孔的「終極形態」:割線跳躍 \(\frac{1}{e^{2\pi i s} - 1}\) 直接編碼 \(\Gamma\) 的極點結構。

場景三:Wiener–Hopf 方法。半直線上的積分方程(如輻射傳輸、翼型繞流)化為「上半平面解析 + 下半平面解析」的因子分解問題——用切割線把函數拆成兩半。1931 年 Wiener–Hopf 創立,至今是波動散射(繞射理論)的標準武器。圍道的角色從「計算」升級為「分解結構」。

場景四:機率與特征函數。特徵函數 \(\varphi(t) = E[e^{itX}]\) 的反演(Lévy 反演公式)是 Fourier 積分;穩定分布(Lévy \(\alpha\) -stable)的密度只能以圍道積分表示(除 Cauchy、Gauss 外無封閉式)。大偏差理論的 Lugannani–Rice 公式用鞍點法計算尾機率——金融風險管理(Tail risk)的計算引擎是複圍道。

場景五:數值圍道積分與 trapezoidal 奇蹟。週期解析函數的複合梯形公式指數收斂(3.1 場景四);非週期積分 \(\int_0^\infty\) 先經指數/雙曲變換變週期,再梯形——tanh-sinh 求積(Takahasi–Mori 1974)正是「把實積分搬到複平面使被積函數解析化」的工業級實現,現代高精度計算(Mathematica NIntegrate 預設之一)靠它。圍道思想在數值分析中的第二次生命。

名詞速查表(本節出現的核心術語)#

  • 單位圓法: \(z = e^{i\theta}\) ,三角積分 → 圓上留數。
  • 大半圓圍道:無窮區間有理積分的標準路線。
  • Jordan 引理:上半圓上 \(f e^{iaz}\) 的大圓貢獻為零。
  • 鑰匙孔圍道:繞割線的圍道,跳躍 = \(2\times\) 目標積分。
  • Fresnel 積分: \(\int_0^\infty e^{ix^2}dx\) ;旋轉圍道技術。
  • 最陡下降法:鞍點主導的漸近計算。
  • Hankel 圍道: \(\Gamma\) 函數的繞負軸圍道。
  • 主值積分:p.v.,奇點在圍道上的對稱取極限。

年表速查(本節涉及的關鍵年份)#

  • 1814:Cauchy 首用複圍道算實積分(雙重積分論文)。
  • 1826:留數定理,圍道法系統化。
  • 1836:Dirichlet 積分 \(\int\frac{\sin x}{x} = \frac\pi2\) (實變證明)。
  • 1878:Fresnel 積分的複變處理成熟(實變 1818 Fresnel 波動論)。
  • 1894:Jordan 引理入教材。
  • 1909:Debye 最陡下降法。
  • 1931:Wiener–Hopf 因子分解。
  • 1974:tanh-sinh 數值求積。

深化問答三則#

問:為什麼三角積分「變數替換」後必須檢查極點在圓內外? 答: \(z = e^{i\theta}\) 把 \(\theta \in [0, 2\pi]\) 映為單位圓;實分母 \(a + b\cos\theta\) 在 \(\theta\) 上無零(若 \(|a|>|b|\) ),但對應的 \(z\) 多項式有兩個根——一內一外。實函數的「無奇點」不等於複化後的「無奇點」:複化使分母因式分解,根的位置由 \(a^2 - b^2\) 的符號決定(正 → 實圓內外各一;負 → 根在單位圓上,積分發散——這正是 \(\int\frac{d\theta}{1+\cos\theta}\) 發散的原因)。

問:Jordan 引理為什麼要 \(a > 0\) 與上半圓配套? 答: \(e^{iaz} = e^{iax - ay}\) :衰減因子 \(e^{-ay}\) 只在上半平面( \(y > 0\) )有效。 \(a < 0\) 時要用下半圓(極點留數取下半 + 積分號反轉)。指數的「優勢半平面」由 \(a\) 的符號決定——這就是 10.1 篂「因果性 = 解析半平面」的萌芽: \(e^{izt}\) 的解析半平面(下)對應 \(t > 0\) 的時間方向。

問:鑰匙孔圍道的「跳躍」為什麼恰是 \(2\times\) 目標積分? 答:割線兩岸的分支值差(如 \(\sqrt z\) : \(\sqrt x - (-\sqrt x) = 2\sqrt x\) )。圍道沿兩岸各走一次、方向相反、積分相加 = 差值 × 目標積分。跳躍 = 分支的多值性殘餘——「多值性」在此轉化為「可提取的資訊」,與 5.1 篂「奇點轉化為留數」同一哲學:把障礙變成訊號。

(本節完)