從一條方程到一族算子
在第 1 份指南裡,你學會把熱問題完全不看作一個關於 x 和 t 的 PDE,而看作單一個點 u(t) 在一個函數空間——像 L^2 這樣的希爾伯特空間——中描出一條曲線,同時一個無界算子 A(帶邊界條件的負拉普拉斯算子)推著它前進。整件事塌縮成那條看似簡單的抽象柯西問題 u_t + A u = 0,且 u(0) = u_0。這份指南回答那條思路所招來的問題:如果這看起來恰恰像純量常微分方程 y' = -a y,我們能不能就用同樣的方式寫下它的解?
對純量常微分方程,答案是微積分裡最有名的公式:y(t) = e^(-a t) y(0)。數 e^(-a t) 是一個乘子,把時間 0 的狀態帶到時間 t 的狀態。這份指南的大膽主張是:同一條公式對 PDE 也成立,只要把純量 a 升格成算子 A——解就是 u(t) = e^(-tA) u_0。難處——也正是整門學科所在——在於 e^(-tA) 不再是一個你能用泰勒級數算出的數,因為 A 是無界的:它微分兩次,而你無法界定那會把一個函數拉伸多少。所以 e^(-tA) 必須先被定義出來,才能被使用。
下面這一步讓一切都站得住腳。忘掉指數的冪級數;只保留一個流必須具備的性質,並把它們當作公理來要求。對每個時間 t >= 0,令 S(t) 是把初始資料映到時間 t 的解的算子:S(t) u_0 = u(t)。於是「解出來再等」這件事的意義,逼出三個事實。等零秒什麼都不改變:S(0) = I,恆等算子。先等 s 秒再等 t 秒,與直接等 s + t 秒相同:S(t) S(s) = S(t + s)。而解應該連續依賴於你等待的時間:S(t) u_0 隨 t 變化而連續地移動。一族滿足這些的算子,就是一個 C0 半群。
在 t = 0 處讀出生成元
我們純粹從流的性質造出了半群 S(t),從未提到 A。那麼方程藏在哪裡?它藏在流如何起步之中。回到純量圖像:e^(-a t) 在 t = 0 等於 1,而它在 t = 0 的導數恰好是 -a。數 -a 是流在時間原點的瞬時變化率,而它足以決定整個指數。半群也一樣,只是這個變化率是一個算子。
在流開始的那一刻對它微分。半群的無窮小生成元,就是當 h 縮向 0 時,由 (S(h) u_0 - u_0) / h 的極限得到的算子——曲線 u(t) 離開 u_0 時的速度。對熱半群,那個極限恰好是 -A u_0:解開始移動的速率,是(負拉普拉斯)算子作用在資料上再取負號。生成元是 -A,而把它復原出來,正是你確認一個給定半群真的解的是你的方程、而非別的方程的方式。
現在整幅圖咬合到位,而且你可以雙向閱讀。生成元 -A 是流在 t = 0 處的導數;半群 e^(-tA) 是你把那個變化率往前積分而長出的流。知道生成元就是知道方程(生成元帶著算子 A 以及它的定義域——包含烤進「A 在何處有定義」之中的邊界條件);知道半群就是一次知道所有解。「方程 <-> 生成元 <-> 半群」這部字典,是整個演化方程觀點的脊梁。
指數為何存在:預解算子與希爾-吉田
我們有了一個流,也有了它的生成元 -A。深層的問題是反過來的那個:給定一個無界算子 A,它在什麼時候是某個半群的生成元——也就是說,e^(-tA) 何時真的存在、成為一族表現良好的算子?你無法靠把泰勒級數加起來來回答,因為那些項含有 A^n,而無界算子沒有可以控制這些冪的範數。古典的技巧是從不直接反轉 A,而是去反轉那個有界的算子 (lambda I + A)。
那個逆算子就是預解算子 R(lambda) = (lambda I + A)^(-1)。A 會無界地拉伸函數,它的預解算子卻溫和地收縮它們——它是有界的,而對擴散算子,它隨 lambda 增大而變小,對每個模態就像 1 / (lambda + a)。精妙之處在於,預解算子是半群的一種拉普拉斯變換:把 e^(-lambda t) S(t) 對 t >= 0 積分,會重現 R(lambda)。所以只要你能控制預解算子,你就控制了半群。希爾-吉田定理把這件事說精確。
希爾-吉田用乾淨的一句話說:-A 生成一個 C0 收縮半群,恰好當 A 是閉的、稠密定義的,且它的預解算子對所有 lambda > 0 都被 1 / lambda 所界。證明本身就是構造:把 A 換成有界的逼近 A_lambda = lambda A R(lambda)(吉田逼近),用收斂的冪級數安全地取它們的指數,再取極限。所以抽象解的存在不是憑信任斷言出來的——它是一個模態一個模態地,從你能在實際橢圓算子上驗證的預解算子界裡建造出來的。
對對稱的擴散算子,有一條你已經能想像的捷徑。因為 -A 是帶(譬如)狄利克雷條件的負拉普拉斯算子,它有一組正交歸一的特徵函數 phi_n,對應特徵值 a_n > 0——正是分離變數法裡那些空間特徵函數。在那組基底上,e^(-tA) 不過就是把第 n 個座標乘以 e^(-a_n t)。被去神秘化的抽象指數,無非就是你寫過上百次的那個衰減傅立葉級數——如今被認出來,是一個悄悄阻尼每個模態的算子。
解析半群:平滑住在哪裡
普通的 C0 半群一般到足以涵蓋波那樣、什麼都不平滑的運動。擴散半群是特殊的——它是解析的——而解析性正是上一份指南裡你親手證明的一切的算子理論名稱。一個解析半群,是那種你能把時間 t 延拓進複平面的一個楔形區域、並愛微分幾次就微分幾次的流;算子 e^(-tA) 對每個 t > 0,都把整個空間映進 A 的定義域、再進 A^2 的定義域,依此類推。
那一句話就是喬裝過的瞬間平滑。回想你逐模態看到的平滑:因子 e^(-a_n t) 把高模態壓得如此之狠,使得結果對任何 t > 0 都具有各階導數。用算子語言,平滑率是那條界 ||A e^(-tA)|| <= C / t——多作用一次 A(多一次導數)只花一個 1 / t 的因子,它在 t -> 0 時爆炸,但對每個正時間都有限。把它迭代,你就瞬間得到每一階導數。半群的解析性,與熱方程的 C-無窮平滑,是同一條定理的兩套戲服。
是什麼讓一個生成元變成解析的?答案又是預解算子,但帶一條更強的界。算子 A 被稱為扇形的,若它的預解算子不只在正實軸上、而是在複平面繞著它的一整個扇形區域內都受控——它的譜落在一個楔形裡,而 (lambda I + A)^(-1) 在整個區域都像 1 / |lambda| 那樣被界。扇形算子恰好就是解析半群的生成元。一致橢圓算子是扇形的;這就是為什麼每個誠實的拋物型方程都會平滑的抽象理由,遠遠超出常係數的拉普拉斯算子。
整幅圖,以及它為你換來什麼
把整個對應關係與你從大一微積分起就熟知的純量常微分方程並排攤開來看。算子故事的每一行,都是數的故事某一行的忠實升級——同樣的形狀,只是 A 取代了 a,外加一條關於定義域的審慎條款,以及一個禁止倒著跑的、缺失的負號。
scalar ODE y' = -a y | abstract PDE u_t + A u = 0 --------------------------------+---------------------------------------- solution y(t) = e^(-a t) y0 | u(t) = e^(-tA) u0 = S(t) u0 flow law e^(-a s) e^(-a t) | S(s) S(t) = S(s+t), S(0) = I rate at 0 d/dt|_0 = -a | d/dt|_0 S(t) u0 = -A u0 (generator -A) exists when (always, a a number)| Hille-Yosida: ||(lam I + A)^-1|| <= 1/lam resolvent 1/(lam + a) | R(lam) = (lam I + A)^(-1) spectrum eigenvalues a_n > 0 | e^(-tA) phi_n = e^(-a_n t) phi_n smoothing (none) | analytic: ||A e^(-tA)|| <= C/t, C-infinity
這層抽象除了優雅,究竟為你換來什麼?三件具體的事。第一,一致性:一條定理一次涵蓋熱方程、變係數擴散、反應-擴散與方程組,因為它們共享一個扇形生成元——你不必為每個新方程重證一次存在性。第二,給非齊次問題 u_t + A u = f 一個乾淨的歸宿,在那裡半群把常數變易公式變成對 e^(-(t-s)A) 的一個精確積分——那正是下一份指南所立足的溫和解公式。第三,免費得到長時間行為:衰減最慢的特徵值,即譜隙,控制著每個解鬆弛到平衡的速度,因為 e^(-tA) 以速率 a_n 阻尼第 n 個模態。
在你把這套帶往下一步之前,給你一個誠實的提醒。那幅乾淨的特徵函數圖像——e^(-tA) 只是把模態 n 乘以 e^(-a_n t)——倚賴 A 是自伴的,這對對稱擴散成立,但一旦你加進像對流這樣的輸運項、使 A 變成非正規算子,它就失效了。那時半群依然存在、依然平滑(解析性真正需要的是扇形性,而非對稱性),但特徵函數未必正交,而且即使每個特徵值都說「衰減」,短時間行為仍可能暫時增長。扇形算子,而非自伴算子,才是這套理論真正的邊界。