9.1 熱方程的譜解法#

本章導言#

本書以熱傳導開篇（1.1 節：Fourier 1807），現在以它的完整現代解法回歸：譜解法（spectral method）——把偏微分方程的解展開為空間特徵函數（正弦、餘弦或更一般的 Sturm–Liouville 函數）的級數，把偏微分方程化為一組獨立的常微分方程。熱方程 \(\frac{\partial u}{\partial t} = \alpha u_{xx}\) 在譜域化為 \(\dot{u}_n = -\alpha\lambda_n u_n\) ——每個模態獨立指數衰減，衰減率由特徵值決定。本節將完整解決熱方程（三種典型邊界條件）、證明「 \(t > 0\) 即時光滑」的奇蹟性質、引入 Green 函數與最大值原理，並以譜方法（數值版）收尾——傅立葉方法最深的用武之地是 PDE，本章與下一節展示其威力。

熱方程的完整譜解法#

問題： \(\frac{\partial u}{\partial t} = \alpha\frac{\partial^2 u}{\partial x^2}\) ， \(0 < x < L\) ， \(t > 0\) ；邊界 \(u(0,t) = u(L,t) = 0\) ；初始 \(u(x,0) = f(x)\) 。

第一步：空間特徵問題（1.1 節已推導）： \(X'' + \lambda X = 0\) 、 \(X(0) = X(L) = 0\) 給 \(\lambda_n = \left(\frac{n\pi}{L}\right)^2\) 、 \(X_n = \sin\frac{n\pi x}{L}\) 。正交性（1.3 節）： \(\int_0^L X_nX_m = \frac{L}{2}\delta_{nm}\) 。

第二步：時間模態：設 \(u = \sum_n T_n(t)X_n(x)\) ，代入方程並投影到 \(X_m\) （正交性使係數分離）：

\[ T_n'(t) = -\alpha\lambda_n T_n(t) \;\Longrightarrow\; T_n(t) = b_n\,e^{-\alpha\lambda_n t}, \qquad b_n = \frac{2}{L}\int_0^L f\,X_n\,dx. \]

解：

\[ u(x,t) = \sum_{n=1}^{\infty} b_n\,\sin\frac{n\pi x}{L}\;e^{-\alpha(n\pi/L)^2t}. \]

三種邊界條件的譜表：

邊界條件	特徵函數	特徵值	物理意義
Dirichlet（ \(u=0\) 兩端）	\(\sin\frac{n\pi x}{L}\)	\(\left(\frac{n\pi}{L}\right)^2\)	定溫端
Neumann（ \(u_x=0\) 兩端）	\(\cos\frac{n\pi x}{L}\) （含 \(n=0\) ）	\(\left(\frac{n\pi}{L}\right)^2\)	絕熱端
週期（ \(u, u_x\) 連續）	\(e^{2\pi i nx/L}\)	\(\left(\frac{2\pi n}{L}\right)^2\)	環形桿

注意 Neumann 的 \(n = 0\) 模態： \(\lambda_0 = 0\) 、 \(T_0 =\) 常數——總熱量守恆（絕熱系統的熱量不變），解長時間收斂到常數（平均溫度）而非零。守恆律在譜域現形為零特徵值——這條對應（守恆律 ↔ kernel）是現代 PDE 理論的基本直覺。

長時間行為與「即時光滑」奇蹟#

長時間主導： \(t \to \infty\) 時 \(u \approx b_1\sin\frac{\pi x}{L}e^{-\alpha\pi^2 t/L^2}\) ——最慢衰減的第一模態統治一切（1.1 節場景四）。物理：高頻（小尺度）溫度起伏以 \(e^{-\alpha n^2 t}\) 極速消失——擴散是高頻殺手。

即時光滑（instantaneous smoothing）：設 \(f\) 僅分段連續（粗糙！）。對 \(t > 0\) ：係數變為 \(b_n e^{-\alpha\lambda_n t}\) ——指數衰減壓制一切，級數可對 \(x\) 微分任意次（逐項微分合法：微分給 \(n^k\) 因子，仍被 \(e^{-\alpha n^2t}\) 鎮壓）。結論：對任何 \(t > 0\) ， \(u(\cdot, t)\) 是 \(C^\infty\) ——甚至解析。粗糙的初始資料，瞬間產生無窮光滑的解——這是拋物方程（熱型）的標誌性奇蹟，與雙曲方程（波型，9.2 節：奇性以波速傳播、永不消失）形成本書最鮮明的對比：

	熱方程（拋物）	波方程（雙曲）
特徵值	\(\lambda_n \sim n^2\)	\(\sqrt{\lambda_n}\sim n\) （頻率）
高頻抑制	\(e^{-\alpha n^2t}\) （超強）	無抑制（僅振盪）
初始奇性	\(t>0\) 即抹平	永遠傳播（特徵線）
時間箭頭	不可逆	可逆
有限傳播速度	無（瞬時傳遍全域）	有（波速 \(c\) ）

有限傳播速度的悖論：熱方程的解在 \(t > 0\) 時全域非零（級數每一項都觸及全區間）——「熱瞬時傳遍宇宙」（經典悖論）。數學真相：非零但指數小（小於任何儀器靈敏度）；物理修正：微觀模型（隨機行走、相對論熱傳）恢復有限速度。譜域看：拋物型的「瞬時全域」= 係數全域非零；雙曲型的「有限速度」= 係數的支撐結構（9.2 節）。

Green 函數：熱核#

把解改寫為單一積分（熱核）：

\[ u(x,t) = \int_0^L G(x, \xi; t)\,f(\xi)\,d\xi, \qquad G(x,\xi;t) = \frac{2}{L}\sum_{n=1}^{\infty}\sin\frac{n\pi x}{L}\sin\frac{n\pi\xi}{L}\,e^{-\alpha(n\pi/L)^2t}. \]

\(G\) = δ 輸入的響應（5.1 節物理意義： \(x = \xi\) 處的點熱源的擴散）。 \(G\) 的性質：（1） \(t \to 0^+\) 時 \(G \to \delta(x - \xi)\) （分配意義）；（2）對稱 \(G(x,\xi) = G(\xi,x)\) （自伴算子的報酬）；（3） \(G > 0\) （最大值原理的表現）；（4）半無限域的極限恰是 Gauss 核（4.1 節補充的 Fourier 原始計算）：

\[ G_\infty(x, t) = \frac{1}{\sqrt{4\pi\alpha t}}e^{-x^2/4\alpha t}. \]

熱核 = 高斯模糊 = 捲積算子（4.3 節）——三個名字、同一個對象：擴散、模糊、平滑化是同一件事在 PDE、影像、分析三界的名字。

一個具體例題#

例題 9.1.1：桿 \(L = \pi\) 、 \(\alpha = 1\) 、兩端零溫，初始溫度 \(f(x) = x(\pi - x)\) （拋物線分布，最大 \(\frac{\pi^2}{4}\) 於中點）。（1）求級數解；（2）估計 \(t = 0.1\) 時保留兩項的誤差；（3）討論最大溫度的時間演化與「熱量哪裡去了」。

解：

（1）係數： \(b_n = \frac{2}{\pi}\int_0^\pi x(\pi - x)\sin nx\,dx\) 。分部積分兩次（ \(f(0) = f(\pi) = 0\) 使邊界項乾淨）： \(b_n = \frac{8}{n^3\pi}\) （**僅奇數 \(n\) **，因 \(f\) 關於 \(x = \frac{\pi}{2}\) 對稱——1.3 節的對稱性技巧）。解：

\[ u(x,t) = \frac{8}{\pi}\sum_{k=0}^{\infty}\frac{1}{(2k+1)^3}\sin\big((2k+1)x\big)\,e^{-(2k+1)^2t}. \]

（2） \(t = 0.1\) ：項衰減因子 \(e^{-0.1}, e^{-0.9}, e^{-2.5}, \ldots = 0.905, 0.407, 0.082, \ldots\) 。係數 \(\frac{8}{\pi}, \frac{8}{27\pi}, \frac{8}{125\pi} = 2.55, 0.094, 0.020\) 。兩項後殘餘 \(\lesssim 0.020\times0.082\times3 \approx 0.005\) （最大溫度的 0.4%）——** \(t > 0\) 的譜級數收斂極快（與 \(t = 0\) 時係數 \(\frac{1}{n^3}\) 的多項式收斂對照：擴散把條件收斂變成指數收斂**）。

（3）最大溫度（中點 \(x = \frac{\pi}{2}\) ， \(\sin\) 值 \(= \pm1\) 交替）：

\[ u\left(\frac{\pi}{2}, t\right) = \frac{8}{\pi}\left(e^{-t} - \frac{e^{-9t}}{27} + \frac{e^{-25t}}{125} - \cdots\right). \]

\(t = 0.1\) ： \(2.31 - 0.015 + 0.0003 = 2.29\) ； \(t = 1\) ： \(2.94\,e^{-1}\approx 1.08\) ； \(t\to\infty\) ： \(\to 0\) 。熱量去向：全域溫度 \(\to 0\) （Dirichlet 端點吸熱——邊界是熱匯）。能量核算： \(\int_0^\pi u\,dx = \frac{16}{\pi}\sum_k\frac{e^{-(2k+1)^2t}}{(2k+1)^2}\) （每項 \(\int\sin nx = \frac{2}{n}\) 僅奇數）， \(t = 0\) 時 \(= \frac{16}{\pi}\cdot\frac{\pi^2}{8} = 2\pi = \int f\) ✓，隨 \(t\) 指數衰減——逐模態的能量收支清晰可讀：這是譜方法「會計透明」的典範。

歷史啟示#

譜解法 = 「對角化 + 指數衰減」的模板：任何自伴算子的發展方程 \(\dot{u} = Au\) 在譜域都是逐模態的標量 ODE。Schrödinger 方程（ \(\dot\psi = -iH\psi\) ）是同一模板的酉版（特徵值變相位而非衰減）——量子力學的形式語言直接繼承了傅立葉的熱學方法（Heisenberg–Schrödinger 1926 年皆自覺引用）。
光滑性是自由參數：拋物方程「任意初始、即時光滑」使它成為數值上最友善的方程（譜方法的理想試驗場），也是反問題最兇險的方程（逆向要放大 \(e^{+\alpha n^2 t}\) ——1.1 節場景五）。
邊界條件決定譜：同一方程，Dirichlet/Neumann/週期給出三套譜——「方程 + 邊界」才是完整的數學對象（Sturm–Liouville 理論的教訓，1.3 節）。

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

1.1 節：本節是開篇故事的完整技術版。
1.3 節：正交性與 Sturm–Liouville——譜解法的代數基礎。
2.2 節： \(t = 0\) 時的收斂（Dirichlet 條件）； \(t > 0\) 的指數收斂（本節）。
4.1、4.3 節：半無限域熱核 = Gauss 核 = 高斯模糊。
5.2 節：Green 函數 = δ 驅動的解；基本解理論。
6.3 節：譜方法的 FFT 實現（週期邊界）。
9.2 節：波方程——拋物/雙曲對照的另一極。
9.3 節：Laplace 變換解 ODE——時間方向的姊妹方法。
10.1 節：半群理論（ \(e^{tA}\) ）——本節解的公理化。
數學物理：Schrödinger、Laplace（穩態極限 \(\Delta u = 0\) ）、Fokker–Planck（隨機過程）全是本節模板的變奏。

練習#

解絕熱桿（Neumann）： \(L = \pi\) 、 \(f(x) = x\) 。求級數解，驗證 \(n = 0\) 項 = 平均溫度 \(\frac{\pi}{2}\) ，且 \(t\to\infty\) 解趨於該常數（總熱量守恆）。
半無限桿 \(x > 0\) 、端點保持零溫：用奇延拓（把 \(f\) 延拓為奇函數）把問題化為全直線，驗證解為 \(u(x,t) = \int_0^\infty [G_\infty(x - \xi) - G_\infty(x + \xi)]f(\xi)d\xi\) ——邊界條件 = 對稱化技巧（與 1.3 節餘弦級數同理）。
對例題 9.1.1 的解，計算 \(t = 1\) 時保留一項的誤差界（用 \(\frac{8}{125\pi}e^{-25}\cdot3\) 估計），體會「長時間下第一模態的暴政」。
證明最大值原理的一個版本：Dirichlet 熱方程的解滿足 \(0 \le u \le \max f\) （對 \(0 \le f\) ）。提示：設最大值在內點 \(t_0 > 0\) 達到，用 \(u_t = 0, u_{xx} \le 0\) 導出矛盾（或用譜表達 + \(G > 0\) ）。
思考題：非齊次方程 \(u_t = \alpha u_{xx} + g(x, t)\) （持續熱源）：用Duhamel 原理（常數變易法的 PDE 版）寫出解 \(u = (\text{齊次演化}) + \int_0^t(\text{熱核} * g(\cdot, s))\,ds\) ，並解讀「源的歷史逐層疊加」。這與 ODE 的 \(y = e^{At}y_0 + \int e^{A(t-s)}g(s)ds\) 完全平行——半群語言（10.1 節）的初等預演。

進階研讀：譜方法與熱方程的五個深化場景#

場景一：譜方法的數值實務。週期邊界 + FFT（6.3）：把 \(u\) 的 FFT 逐模態乘 \(e^{-\alpha n^2t}\Delta t\) （顯式）或用精確積分因子（ETD，exponential time differencing——對剛性問題（stiffness，高頻極快衰減）特別有效）。剛性： \(n_{\max}^2\) 量級的衰減率使顯式時間步受限 \(\Delta t \lesssim \frac{\Delta x^2}{2\alpha}\) （CFL 型條件）——譜方法空間精確、時間仍受剛性折磨，ETD/隱式格式是解藥。**空間誤差指數衰減（3.2 節譜精度）+ 時間誤差 \(O(\Delta t^p)\) ** 的分解是譜方法誤差分析的標準框架。

場景二：多維與分離的極限。二維矩形域： \(u_{mn}\) 、 \(\lambda_{mn} = \lambda_m^x + \lambda_n^y\) （可分離）；二維圓盤：Bessel 函數 \(J_m(k r)\) （不可分離但可分變數——Sturm–Liouville 的 Bessel 版本，9.2 節繼續）；一般域：譜幾何問題（「能否聽出鼓的形狀？」Kac 1966——特徵值決定域的幾何到什麼程度？1992 年 Gordon–Webb–Wolpert 證明不能完全——存在同譜異形域！）。譜是域的指紋，但不唯一——熱方程把幾何問題轉為數列問題的能力與限度。

場景三：隨機視角——Brown 運動與 Fokker–Planck。熱方程 = Brown 運動的機率密度方程： \(u(x,t) = E[f(x + \sqrt{2\alpha t}W)]\) （W = 維納過程）——譜解法的隨機對偶：每個特徵模態對應隨機行走的低頻份額。Feynman–Kac 公式把熱方程的解表為隨機泛函的期望——PDE 與機率過程在熱方程處是同一枚硬幣（Kolmogorov 1931、Feynman 1948、Kac 1949）。量子力學的路徑積分是它的虛時間版。

場景四：逆時間與不適定。反熱傳（由 \(u(\cdot, T)\) 推 \(f\) ）：放大 \(e^{+\alpha n^2 T}\) ——任何測量噪聲在 \(n\) 大處爆炸（Hadamard 不適定，1.1 節場景五）。對策：正則化（Tikhonov：截斷高頻 + 罰項；11 世紀的方法論：「用偏差換方差」）——地熱勘探、衛星熱成像的反演都是這場與指數的戰爭。抛物方程的正反不對稱是時間之箭的最純數學形式。

場景五：非線性熱方程與 Blow-up。 \(u_t = \alpha u_{xx} + u^2\) （反應項）：解可能有限時間爆破（ \(u \to \infty\) 於 \(T^* < \infty\) ，如 \(u(0,x) = \frac{1}{(T^* - t)}\) 型解）——擴散 vs. 反應的拔河：擴散拉平（ \(n^2\) 衰減）、反應放大（非線性乘法），孰勝取決於維數與初始條件（Fujita 臨界指數 1966： \(p > 1 + \frac{2}{d}\) 時任何非平凡初始都爆破）。譜方法 + 能量法是爆破分析的標準工具——非線性世界裡「模態分離」失效（模態互相饋能），譜只提供座標、物理靠能量估計——線性方法的邊界在非線性處清晰顯形。

補充：Fourier 的 1824 年「行星溫度」論文——溫室效應的數學黎明#

Fourier 在 1824 年（與 1827 年修訂）發表〈地球及其行星空間溫度的概略〉，用半無限介質的熱方程（本節場景二的邊值版）推理：地球接收太陽輻射、以紅外線向太空散熱，若大氣對兩種波長的透明度不同（「如同一個罩子」，他的比喻：玻璃罩、高山湖泊），地表溫度將高於無大氣的輻射平衡值。這是溫室效應理論的第一份定量數學論證（現代精確機制——吸收頻譜——由 Tyndall 1859、Arrhenius 1896 完成，但問題的形式——熱平衡 + 波長選擇性——是 Fourier 的）。傅立葉分析與氣候科學的緣分：從 1824 年的熱方程到今日的輻射傳輸方程與衛星資料的頻譜反演，兩百年同一條數學血脈。

名詞速查表（本節出現的核心術語）#

譜解法：特徵函數展開 + 逐模態 ODE。
Dirichlet / Neumann / 週期譜：三種邊界的三套特徵值。
零特徵值 ↔ 守恆律：Neumann 的總熱量守恆。
即時光滑：拋物方程 \(t > 0\) 的 \(C^\infty\) 奇蹟。
熱核 / Green 函數：δ 響應； \(G = \frac{2}{L}\sum\sin\sin e^{-\alpha\lambda t}\) 。
高斯核：半無限域極限；擴散 = 模糊。
最大值原理：解界於初始與邊界的最大值。
Duhamel 原理：非齊次源的疊加。
剛性（stiffness）：多尺度衰減率的數值災難；ETD 對策。
逆時間不適定： \(e^{+\alpha n^2t}\) 的放大災難；正則化對策。
半群 \(e^{tA}\) ：線性發展方程的解算子（10.1 節）。
Fujita 指數：非線性熱方程爆破的臨界維數指數。

年表速查（本節涉及的關鍵年份）#

1807 / 1822：Fourier 的熱傳導理論與《熱的解析理論》。
1824 / 1827：Fourier 行星溫度論文（溫室效應先聲）。
1836–38：Sturm–Liouville 理論。
1859 / 1896：Tyndall 吸收光譜；Arrhenius 溫室定量。
1905–10：Einstein、Smoluchowski 的 Brown 運動理論。
1931：Kolmogorov 擴散過程的 Fokker–Planck 框架。
1948–49：Feynman 路徑積分；Kac 公式。
1950s–60s：譜方法數值化（Lanczos、Fox）；剛性理論（Curtiss–Hirschfelder 1952）。
1966：Kac「聽鼓」問題；Fujita 爆破臨界。
1992：同譜異形反例（Gordon–Webb–Wolpert）。
2000s–：ETD、隱式譜方法在流體與氣候模式的日常化。

深化問答三則#

問：為什麼熱方程的解「無限快」傳播而波方程不會？ 答：譜域最透徹：熱方程每個模態 \(e^{-\alpha n^2t}\) 從 \(t = 0^+\) 就全域非零（正弦函數無緊支撐）；波方程解 \(u = \frac{\phi(x - ct) + \phi(x + ct)}{2}\) 的支撐以速度 \(c\) 擴張（d'Alembert 公式的平移結構）。拋物 = 無限速度 + 指數衰減；雙曲 = 有限速度 + 無衰減——兩種「傳播哲學」。微觀上，熱傳導是隨機行走（每步有限但無停頓）、波是確定彈道——隨機性沒有時速。

問：譜方法與有限差分該選哪個？ 答：光滑問題用譜（指數收斂，3.2 節）、奇性/複雜幾何用差分/有限元素（魯棒、局部自適應）。判據就是 2.1 節價碼表：解的光滑度 → 係數衰減 → 譜收斂率。實務混合（譜元素法、激波捕捉 + 譜基座）是現代計算流體的主流。方法的選擇是對「解的光滑度先驗」的押注。

問：本節與 1.1 節同一個例子，學生兩次學到什麼不同？ 答：1.1 節學「怎麼算」（分離變數、係數、級數）；本節學「為什麼對、為什麼好」：正交性的合法性（1.3）、收斂的意義（第 2 章）、即時光滑與 Green 函數（本節新）、數值與反問題（本節新）。同一個例子是本書的螺旋主軸——每次重訪，理論裝備升級一層。這也是 Fejér 的教學哲學（3.3 節補充）：從具體例子出發，讓理論長在例子上。

（本節完）