1.1 Fourier 與熱傳導 1807#
本章導言#
1807 年 12 月 21 日,在巴黎科學院(Institut de France)的會議廳裡,一位四十歲的前神學院學生、現任伊澤爾省(Isère)行政長官,宣讀了一篇關於「熱的傳播」的論文。這位作者是 Joseph Fourier(傅立葉,1768–1830),他論文的核心主張震驚了在場的兩位大數學家:任何函數——哪怕它在某些點是折線、在某些點是間斷的——都可以表示為無窮多個正弦與餘弦函數之和。這個主張後來被稱為傅立葉級數的起源,它不僅解決了熱傳導問題,更在此後兩百年間徹底改變了分析學、工程學與人類文明處理「訊號」的方式。本節將還原這段歷史現場,並說明熱傳導方程為何是傅立葉方法的原點。
Fourier 的時代與其人#
Joseph Fourier 的出身平凡而坎坷:1768 年生於法國 Auxerre,八歲成為孤兒,由教會撫養長大。他少年時原本準備成為本篤會(Benedictine)修士,卻在大革命爆發(1789)後投身政治,一度因立場溫和而被羅伯斯比(Robespierre)政權逮捕,差點上了斷頭台——1794 年 Robespierre 倒台前夕,Fourier 已被判死刑,因 Robespierre 失勢才倖免於難。
命運的轉折發生在 1798 年:拿破崙(Napoléon Bonaparte,1769–1821) 率軍遠征埃及,Fourier 被任命為隨行的科學團成員(他曾在師築學院受教於 Lagrange 與 Monge)。在埃及,Fourier 負責行政與考古調查,主持撰寫了巨著《埃及記述》(Description de l'Égypte,1809–1828)。1801 年返回法國後,他被任命為伊澤爾省省長,駐節格勒諾布爾(Grenoble)。正是在公務之餘的夜晚,這位省長完成了他改變數學史的研究——熱在固體中如何傳播。
這個研究動機本身值得注意:十八世紀末的物理學家對「熱」的本體尚無共識。主流學說是熱質說(caloric theory):熱是一種自相排斥、無重量的流體。Fourier 的論文並未正面回答熱的本體問題,而是採取了一個務實立場:不論熱是什麼,它的流動遵循可測量的定量規律。這種「先寫方程、後談本體」的態度,是數學物理方法的里程碑。
熱傳導方程的誕生#
Fourier 的物理模型如下:考慮一根均勻細桿(或一個均勻物體),記溫度為 \(u(x,t)\) 。實驗(Newton 的冷卻實驗、以及他自己的環形金屬板實驗)提示三個經驗規律:
- 熱量從高溫流向低溫,流速與溫度落差成正比(今日的 Fourier 熱傳導定律: \(\varphi = -k\,\dfrac{\partial u}{\partial x}\) , \(\varphi\) 為熱流密度);
- 微元體的溫度變化等於流入與流出熱量之差;
- 物質的熱容量 \(c\) 與熱導率 \(k\) 是材料常數。
綜合這三點,對一維桿(忽略側面散熱),對區間 \([x-h,\,x+h]\) 做熱量收支平衡,令 \(h \to 0\) ,即得著名的熱傳導方程:
這條方程的深刻性在於:它是一個時間不可逆的方程——熱會自動抹平溫度差,但倒放影片會出現「熱自發集中」的荒謬畫面。熱方程因此成為 19 世紀科學界討論時間之箭、以及後來熱力學第二定律的數學原點之一。
邊界條件與分離變數#
方程只有配上邊界條件(boundary condition)與初始條件(initial condition)才能唯一確定解。以兩端保持零溫的桿為例:
Fourier 的關鍵一招是分離變數(separation of variables):假設 \(u(x,t) = X(x)\,T(t)\) ,代入方程得
左邊只依賴 \(t\) 、右邊只依賴 \(x\) ,兩者相等就必須都是常數 \(-\lambda\) 。於是得到兩個常微分方程:
第二條方程配合邊界條件,是一個特徵值問題(eigenvalue problem):並非所有 \(\lambda\) 都允許非零解,只有
才行。這就是離散譜的第一次登場:邊界條件像琴弦的兩個固定端,只讓特定「音調」存在。每個對應的解 \(\sin\frac{n\pi x}{L}\,e^{-\alpha(n\pi/L)^2 t}\) 都滿足方程與邊界條件。
疊加與「任意函數」的爆炸性主張#
熱方程是線性的:解的任意倍數與相加仍是解。因此
仍滿足方程與邊界條件。剩下的問題只有一個:能否選擇係數 \(b_n\) ,使得 \(t=0\) 時 \(u(x,0) = f(x)\) ? 亦即:
Fourier 給出了求係數的方法(將在 1.3 節系統化),並且斷言:這個表示對任何初始溫度分布 \(f(x)\) 都有效——包括折線函數、分段常數函數、甚至有跳躍間斷的函數。
這個「任意函數可以展開為三角級數」的主張,是整個故事的引爆點。因為在當時的觀念裡,「函數」與「級數」是兩個階層分明的概念:函數是解析表達式(如 \(x^2\) 、 \(\sin x\) ),級數是函數的展開物。說一段折線、甚至一個間斷函數「等於」無窮多個光滑正弦波之和,等於宣告函數概念必須擴張——函數不必是解析式,它可以是一段任意的曲線、一張表格、一條折線。這一概念革命,被後世史家視為現代「函數論」的起點。
1807 年評審現場:大師們的反應#
Fourier 於 1807 年 12 月 21 日向科學院提交論文《熱在固體中的傳播》(Mémoire sur la propagation de la chaleur dans les corps solides)。科學院指定四位評審:Lagrange(拉格朗日,1736–1813)、Laplace(拉普拉斯,1749–1827)、Monge(蒙日,1746–1818)、S. F. Lacroix。他們的反應各自不同:
| 評審 | 態度 | 異議焦點 |
|---|---|---|
| Lagrange | 反對 | 三角級數能否表示任意(含間斷)函數;他數十年來堅持只有連續的解析函數才可展開 |
| Laplace | 保留 | 對熱傳導方程的推導(微分層次)有方法疑慮,後來在分子層面另尋模型 |
| Monge | 友善 | 曾是 Fourier 的師長,政治上也同屬拿破崙陣營 |
| Lacroix | 溫和 | 無重大反對意見 |
Lagrange 的反對尤為關鍵:他早在 1759 年研究振動弦時就與 Euler、d'Alembert 打過著名的「弦振動論戰」,立場正是「三角級數不可能表示含奇點的函數」。這場延續半個世紀的爭議,將在本書 1.2 節與 9.2 節詳述。
評審的結論是:論文的方法與計算「值得關注」,但其中部分內容的嚴密性不足,論文並未獲得出版。1808 年,Lagrange 與 Laplace 的擁護者 Poisson 與 Biot 還分別發表了針鋒相對的批評。Fourier 的回應是繼續工作:他在格勒諾布爾又寫了十幾年,於 1811 年將修訂版投給科學院的熱學大獎競賽,獲得首獎(1812 年頒發),但評審(含 Lagrange)仍以「嚴密性」為由不發表。最終,Fourier 把畢生成果熔鑄成《熱的解析理論》(Théorie analytique de la chaleur,1822 年出版)——這本書成為數學物理的經典,也是傅立葉分析的正式誕生證書。
一個具體例題#
例題 1.1.1:長為 \(L=\pi\) 的桿,兩端溫度恆為零,初始溫度分布為
(左半段加熱、右半段冷卻——一段有跳躍間斷的初始條件。)求溫度 \(u(x,t)\) 的級數解,並說明 \(t>0\) 時的行為。
解:由分離變數,特徵函數為 \(\sin(nx)\) ( \(L=\pi\) ),一般解為
其中係數需滿足 \(f(x) = \sum_{n=1}^{\infty} b_n \sin(nx)\) 。用將在 1.3 節證明的正交性公式,得
逐一列出: \(n=1\) : \(b_1 = \frac{2}{\pi}(1-0) = \frac{2}{\pi}\) ; \(n=2\) : \(b_2 = \frac{2}{2\pi}(1 - \cos\pi) = \frac{2}{\pi}\) ; \(n=3\) : \(b_3 = \frac{2}{3\pi}(1-0) = \frac{2}{3\pi}\) ; \(n=4\) : \(b_4 = \frac{2}{4\pi}(1 - \cos 2\pi) = 0\) 。因此
注意三個現象:(1)每項都隨時間指數衰減,且頻率越高的項衰減越快( \(e^{-\alpha n^2 t}\) )——因此 \(t>0\) 之後級數收斂得極快,間斷點被迅速「抹平」,這正是熱方程的平滑效應;(2)在 \(x = \pi/2\) 這個初始間斷點上,級數收斂到 \(\frac{1}{2}\) ,恰好是左右極限的平均——這個「中點值」現象將在 2.2 節 Dirichlet 定理中獲得嚴格證明;(3)在 \(t=0\) ,級數在間斷點附近收斂得很慢( \(b_n\) 只按 \(1/n\) 衰減),正弦波要「費力」疊出一個直角跳躍——這是 Gibbs 現象(2.3 節)的溫床。
這個例子完整展示了 Fourier 的野心:連一段「直角」的溫度分布,都能用平滑的琴弦音調疊出來。
歷史啟示#
- 應用問題孕育抽象理論:函數概念、積分概念、特徵值概念、無窮級數的收斂概念,都因熱傳導這個「不體面」的工程問題被逼上絕路而重生。數學史上「純出於內部邏輯」的重大概念其實罕見,Fourier 是最好的反例證明。
- 評審的保守與作者的堅持:Lagrange 並非愚蠢,他守護的「解析函數」觀念在他那個時代是自洽的;但 Fourier 用物理實驗與大量計算頂住了權威。科學的進步有時不靠說服,而靠熬時間——Fourier 等了十五年才讓成果問世。
- 公務員與數學家可以合一:Fourier 白天治理一省、推動沼澤排水與疫苗接種,晚上研究熱方程。數學研究未必需要職位名分。
與其他章節及領域的關聯#
- 1.2 節:Lagrange 為何反對?「弦振動論戰」的歷史前因與數學實質。
- 1.3 節:正交函數系與係數公式的系統化——本節例題中 \(b_n\) 公式的來源。
- 第 2 章:收斂問題的正式解決(Dirichlet 定理、Gibbs 現象)。
- 9.1 節:熱方程譜解法的完整現代敘述。
- 9.2 節:波方程與弦振動論戰——Fourier 方法的前身。
- 微積分教材 5.1–5.3 節:無窮級數與一致收斂的基礎概念在本書全面動用。
- 物理(熱學):Fourier 定律 \(\varphi = -k\nabla u\) 至今是熱傳工程的基本公式。
- 資訊(數位訊號處理):本節的「任意波形 = 音調疊加」是所有頻域分析的哲學起點。
練習#
- 用分離變數法求解:桿長 \(L = \pi\) 、兩端零溫、初始溫度 \(f(x) = \sin 2x\) 。為何此例不需要求係數積分?
- 對初始溫度 \(f(x) = x\) ( \(0 < x < \pi\) )計算係數 \(b_n = \frac{2}{\pi}\int_0^\pi x\sin(nx)\,dx\) (提示:分部積分),寫出級數解的前四項。
- 說明為何熱方程的解項 \(e^{-\alpha(n\pi/L)^2 t}\) 顯示「高頻成分先消失」,並用此解釋:為何一張模糊的照片可以理解為高頻資訊的丟失。
- 查閱並摘述 Fourier 1811 年得獎論文被「不發表」的經過,討論 Lagrange 反對的合理性與局限。
- 思考題:熱方程是時間不可逆的,但力學(Newton 方程)是時間可逆的。這個差異在數學上源於方程對 \(t\) 的幾階導數?對 \(x\) 呢?
進階研讀:熱之哲學的五個深化場景#
場景一:非齊次邊界條件。若兩端溫度是固定常數 \(u(0,t)=T_0\) 、 \(u(L,t)=T_1\) ,先找一個不隨時間變的「穩態解」 \(u_s(x) = T_0 + (T_1-T_0)x/L\) 滿足邊界條件,再令 \(v = u - u_s\) 為零邊界問題——這個「穩態 + 瞬態」分解是所有 PDE 求解的標準流程,讀者應在 9.1 節見到其一般形式。
場景二:絕緣端點。若桿端絕緣(熱流為零),邊界條件變成 \(X'(0) = X'(L) = 0\) ,特徵函數改為 \(\cos\frac{n\pi x}{L}\) (含常數 \(n=0\) 項)。同一條方程、不同邊界、不同譜——譜由邊界條件決定,這是 Sturm–Liouville 理論(1.3 節)的核心訊息。
場景三:環形桿(週期邊界)。若把桿彎成環, \(u(0,t) = u(L,t)\) 且 \(u'(0,t) = u'(L,t)\) ,特徵函數變成 \(\sin\frac{2n\pi x}{L}\) 與 \(\cos\frac{2n\pi x}{L}\) ——這正是本書第 2 章採用 \(2\pi\) 週期標準形式的物理來源:圓環是最自然的週期邊界。
場景四:指數衰減與長時間行為。 \(t\to\infty\) 時,級數由衰減最慢的第一項主導: \(u(x,t) \approx b_1 \sin\frac{\pi x}{L} e^{-\alpha(\pi/L)^2 t}\) 。這種「長時間由最低模態主導」的思想,在量子力學(基態)、結構工程(主振型)、化學動力學(慢反應通道)中反覆出現,值得讀者把 \(e^{-\alpha n^2 t}\) 的衰減表( \(n=1..5\) ,取 \(\alpha t = 0.1\) )動手算一遍: \(0.730,\ 0.207,\ 0.043,\ 0.0067,\ 0.0008\) ——五個模態後幾乎歸零。
場景五:反問題與時間倒推。熱方程正向解存在且唯一,但「反向」——由 \(t=T\) 的溫度推回 \(t=0\) ——是不適定(ill-posed)的:高頻成分要放大 \(e^{+\alpha n^2 T}\) ,測量誤差會爆炸。Hadamard 的不適定概念(1902)正是由反熱傳問題催生;今日 CT 影像重建、去模糊演算法都是與之搏鬥的現代戰場。
補充:熱質說的平反#
Fourier 使用熱質說語言(「熱流體在物體內流動」)建構方程,而熱質說在 1840 年代被 Joule 的熱功當量實驗推翻——但方程存活了下來。這展示數學物理的一個普遍規律:方程的正確性高於其導出時所依據的物理圖像。類似案例:Maxwell 從以太模型導出電磁方程、Schrödinger 從經典哈密頓類比導出波動方程——模型可棄,方程長存。讀者學習本書時也應採取這種態度:緊抓方程的結構性質(線性、疊加、對稱性),放寬對來源圖像的執著。
補充:Fourier 晚年與小插曲#
Fourier 在 1822 年當選科學院終身秘書,1830 年因心臟病(一說是他早年埃及生活的健康後遺症——他特別畏寒,常穿多層大衣、待在高溫房間)去世。他在晚年關照了年輕的 Sophie Germain(1776–1831) 與 Peter Gustav Lejeune Dirichlet(1805–1859)——後者於 1820 年代從德國來巴黎,隨身帶著一本被翻爛的《熱的解析理論》,最終用嚴格方法奠定了傅立葉級數收斂理論(2.2 節)。一門數學的血脈,就這樣從省長的書桌傳到普魯士的講堂。
名詞速查表(本節出現的核心術語)#
- 熱傳導方程: \(\dfrac{\partial u}{\partial t} = \alpha \dfrac{\partial^2 u}{\partial x^2}\) ,描述溫度隨時間的擴散。
- 熱擴散率: \(\alpha = k/(c\rho)\) ,決定熱傳播快慢的材料常數。
- 分離變數:設 \(u = X(x)T(t)\) ,將 PDE 化為 ODE 的方法。
- 特徵值問題:微分方程 + 邊界條件,僅在離散的參數值(特徵值)下有非零解。
- 離散譜:允許的 \(\lambda_n\) 形成可數集合,如 \((n\pi/L)^2\) 。
- 疊加原理:線性方程的解可加權相加。
- 穩態解:不隨時間變化的解,滿足 \(u_{xx} = 0\) 。
- 邊界條件/初始條件:在空間端點與時間起點對解施加的約束。
- 不適定問題:解不唯一、不連續依賴於資料或不存在的問題(反熱傳為範例)。
年表速查(本節涉及的關鍵年份)#
- 1768:Fourier 生於 Auxerre。
- 1789:法國大革命爆發,Fourier 投身政治。
- 1794:Fourier 差點上斷頭台,後入師築學院受教於 Lagrange。
- 1798–1801:拿破崙埃及遠征,Fourier 隨行主持科學與行政。
- 1802:任伊澤爾省省長,駐格勒諾布爾;開始熱學研究。
- 1807:12 月 21 日向巴黎科學院宣讀熱傳播論文,Lagrange 反對。
- 1811:以修訂論文獲科學院熱學大獎,仍未發表。
- 1822:《熱的解析理論》出版,傅立葉分析正式誕生。
- 1824–1826:發表「行星溫度」論文,成為溫室效應理論的先聲。
- 1829:Dirichlet 發表傅立葉級數收斂定理(見 2.2 節)。
- 1830:Fourier 逝世於巴黎。
深化問答三則#
問:Fourier 的「任意函數」主張錯了嗎? 答:以當時的語言,他確實陳述過度(例如對非週期函數的展開範圍、對「任意」的邊界都含糊)。但他的方向正確:對合理的 \(f\) (分段連續、有界變差等),三角級數表示確實成立。數學史的常態是先有大方向的正確斷言,再花半世紀劃定精確邊界——Dirichlet(1829)、Riemann(1854)、Lebesgue(1902)都在這條邊界劃定線上。
問:為什麼正弦函數特別「好」,值得讓一切函數向它們臣服? 答:三個原因。其一,正弦是常係數線性 ODE \(X'' + \lambda X = 0\) 的自然解——振動的物理;其二,它們互相正交(1.3 節),使求係數成為「投影」而非解方程組;其三,它們是平移本徵函數—— \(\sin(n(x+a))\) 仍是 \(\sin nx\) 與 \(\cos nx\) 的組合,這使捲積與濾波(4.3 節)在頻域下變成簡單乘法。三個原因合起來:正弦基把「微分」變成「乘法」,把「疊加」變成「加法」。
問:1807 年評審若全是年輕數學家,歷史會加快嗎? 答:可能加快出版,但不會加快嚴格化。Fourier 論文的真正障礙不在「能不能表示」而在「收斂是什麼意思」——這需要 Cauchy(1821)的極限語言與 Dirichlet(1829)的收斂定理才可能回答。一個概念的接受需要兩件事:動人的例子與可用的語言;例子 1807 年就有了,語言還要等二十年。
(本節完)