藏在一條公式裡的三項特質
到現在你已能用三種不同方式解熱方程 u_t = k u_xx:在區間上用分離變數法展成一個衰減的傅立葉級數、在整條直線上與高斯熱核卷積、以及把它的穩態認作一個調和極限。這份指南做的事不一樣。它不再多產出一個解,而是退一步問:這些解全部有什麼共同點?擴散的性格是什麼?是哪一組無論你怎麼求出解、每個解都會穿戴的特質?
共有三項,而且它們不是各自獨立的裝飾——它們是同一套機制的三種觀看角度。熱方程會平滑:餵給它你想要的最鋸齒的資料,下一瞬間它就無限可微。它以無限速度傳播:任何地方的擾動會瞬間在各處被感受到。而它是不可逆的:你一般無法把它倒著跑來還原過去。平滑、無限速度與不可逆性,是每個拋物型方程的招牌,而一旦你看出它們如何從同一個事實掉出來,熱方程就不再是一條公式,而成為一個故事。
平滑:粗糙在一瞬間被消滅
想像一個帶尖角、甚至帶跳躍的溫度分布——一個帶兩道斷崖邊緣的方波熱脈衝。打開熱方程,這些角點不會在幾分鐘內漸漸變柔;就在時間變正的那一瞬間,分布就完美平滑,具有各階導數。這就是平滑效應,是擴散免費為你做的最令人驚訝的事。
你的三種解法中有兩種解釋了它,而且它們彼此一致。透過熱核:解 u(x,t) = integral of G(x - y, t) f(y) dy 是你的資料 f 與高斯 G 的卷積,而 G 無限平滑(事實上是實解析的)。卷積會繼承較平滑那一因子的平滑度,所以即使 f 只是帶跳躍,對每個 t > 0,u 在 x 上都是 C-無窮。透過傅立葉:粗糙的資料之所以粗糙,正是因為它由大量高頻模態組成,而因子 e^(-k n^2 t) 壓垮頻率 n 的模態。對任何 t > 0,指數裡的 n^2 幾乎完全殲滅高諧波,剩下的——一個沒有快速波折的函數——是平滑的。同一個特質,兩扇觀看它的窗。
這是 PDE 各型之間最深刻的分界線之一,所以把它對著它的反面來看。波動方程 u_tt = c^2 u_xx 不會平滑:撥動的弦初始形狀裡的一個尖角始終是尖角,並沿一條特徵線永遠行進——奇異點被傳輸,而非被抹去。奇異點如何表現的這個差別,正是拋物型對雙曲型的分裂。擴散瞬間原諒粗糙;波則永遠記得它。
無限速度:每一點都立刻感受到它
在一根很長很冷的金屬棒一端點燃火柴。常識說遠端會冷上一陣子。熱方程卻說了驚人的事:遠端會立刻升溫——升一個極其微小、無法察覺的量,但嚴格大於零,就在火柴點燃後的那一瞬間。這就是無限傳播速度,它源自熱核的一個性質。
看高斯 G(x,t) = (1 / sqrt(4 pi k t)) e^(-x^2 / (4 k t))。對任何 t > 0(無論多小)和任何 x(無論多遠),G(x, t) 都嚴格為正——一個實數的指數從不為零。所以在 u(x,t) = integral of G(x - y, t) f(y) dy 中,若初始熱 f 在任何地方為正,積分在之後每一刻、各處都為正。沒有波前,沒有「熱目前到達的最後一點、之外恰為零」這回事。熱方程沒有速度上限。
把它擺在波動方程旁邊,那是有限速度的範本。u_tt = c^2 u_xx 裡的訊號以確定速度 c 行進;它尚未到達的區域確確實實、可證明地未被擾動,由方程的特徵線圍起。波動方程有銳利的依賴域;熱方程沒有。拋物型方程瞬間散開;雙曲型方程服從一個有限速度上限——這是方程的「型」最乾淨的實務後果之一。
誠實面對這在物理上的意思:自然界沒有東西真的無限快,所以這是熱方程唯一說了個小謊的地方。這是傅立葉定律假設通量瞬間回應梯度的代價。模型之所以存活,是因為衝在前面的熱量呈指數級微小——高斯尾巴 e^(-x^2/(4 k t)) 在大 x 處極其微小,遠低於你能測到的任何水準,直到經過合理的擴散時間 x^2/k。若你真的需要有限速度,就必須改模型(雙曲型的電報方程加入一個鬆弛時間並把它恢復)。
不可逆性:時間之箭
把一滴奶油攪進咖啡。它混成均勻的淺褐色,而且不管你等多久,你永遠看不到淺褐色的咖啡自發地「反混合」回一滴銳利的奶油。擴散只往一個方向跑。擴散的不可逆性是這個日常確定性的數學面孔,而它不過就是從時間的另一端來看的平滑效應。
這就是兩者為何是同一個事實。向前時間時,每個模態以 e^(-k n^2 t) 衰減,所以細節穩定地被摧毀,許多不同的粗糙過去都流向相似的平滑未來。從過去到未來的映射丟失資訊——而一個丟失資訊的映射無法被反轉。要往後走,你得把每個模態乘以 e^(+k n^2 t),以一個爆炸式的因子復活丟失的高頻:現在 10^(-10) 量級的一個變化,會對應到推得過去的劇烈、無界的變化。向前平滑,正是讓過去無法復原的原因。
這是整門學科的警世故事。逆向熱方程 u_t = -k u_xx 看起來像個無辜的請求——找出模糊成這張圖之前的較早圖像——但它在哈達瑪的精確意義下是不適定的:資料任意微小的變化會產生答案任意大的變化,所以解不連續依賴於資料。留意這審慎的措辭:不適定不代表不存在解;對理想無雜訊資料可能存在唯一解。它的意思是,對任何真實、不完美的測量,答案都派不上用場,因為你試圖放大的丟失細節被淹沒在雜訊裡。
這也是為什麼許多誘人的真實問題在其原始形式下悄然無望。為照片去模糊、熱回溯以找出一個物體一小時前的溫度、從之後的感測器資料辨識熱源——這些全都是逆向擴散。沒有一個能靠字面上把方程倒著跑來解決;它們只能被正則化馴服,也就是你用一個鄰近的適定問題取代精確的逆向問題,加入額外假設(平滑度界、懲罰項)來壓制爆炸的高頻。正則化並沒有把時間之箭倒轉;它回答的是一個刻意修改過的問題。
一張把結繫起來的圖
一次看見這三項特質最乾淨的方式,是觀看一團粗糙的熱演化,並在每個階段問熱核對它做了什麼。那個會平滑的高斯卷積,正是那個會到達各處的嚴格為正的高斯,也正是那個無法復原的、摧毀資訊的平均。把底下這張圖讀作一個方波熱脈衝落在一條冷直線上。
initial f(x): ___|-----|___ (two sharp jumps: a square pulse)
t > 0 (any t): u(x,t) = convolve f with the Gaussian
G(x,t) = e^(-x^2/(4 k t)) / sqrt(4 pi k t)
SMOOTHING corners gone instantly -> C-infinity hump, for every t>0
INFINITE SPEED G(x,t) > 0 for ALL x -> u(x,t) > 0 everywhere at once
IRREVERSIBLE the average loses detail -> a backward step needs e^(+k n^2 t): blows up
later pictures: .-=#####=-. then .-===-. then .-_-. (spreading, fading bell)
width ~ sqrt(4 k t); area = total heat = conserved留意底行隨之而來的兩個安靜的守恆事實。寬度隨 sqrt(4 k t) 成長——你一路看到的 x ~ sqrt(k t) 擴散尺度——所以把散布加倍要花四倍時間;而鐘形下的面積保持固定,正是讓熱核成為「被重新分配但從未丟失的一單位熱量」的總熱量守恆。那個會平滑、會到達各處的同一條高斯,恰好也是每團局部的熱最終鬆弛而成的自相似分布,它已忘記自己起始於何處的細節。
這通往何方
這三項特質是定性的圖像;下一份指南會把其中之一變成一條滴水不漏的定理。無限速度與正性暗示熱方程無法憑空製造出一個新的熱點——熱只能從既有的尖峰往下流。說精確了,那就是拋物型最大值原理:解在一個時空區域上的最大值,總是在最一開始或在邊界上取得,絕不會在內部新鮮地冒出來。它是你一直憑信任使用的兩件事——唯一性與連續依賴——背後的嚴謹引擎,也因此是向前熱問題適定的真正原因。
那份最終指南也會回到源項。這裡的一切都假設棒子內沒有加熱器;打開一個,你就有了非齊次方程 u_t = k u_xx + f。杜阿梅原理處理它的方法,是把源視為一連串微小的瞬間脈衝,每一個接著就如同一個全新的初始條件那樣,恰好地平滑、散開並守恆。你剛認識的這三項特質不是一個觀光景點——它們是整個拋物型世界的工作文法,而你會在往後遇到的每個擴散方程上不斷讀出它們。