在圓盤上求解：帕松積分公式

從抽象的解到顯式的解

到現在，你對調和函數已經懂得不少，卻還一個也沒真正寫出來過。你知道狄利克雷問題——在一個區域的邊界一圈規定溫度，求內部的穩態溫度——恰好有唯一一個解，也知道均值性質與最大值原理擔保了它唯一且行為良好。但存在性的證明與唯一性的證明，都是欠條：它們發誓答案就在某處，卻從不讓你見它一面。這篇導讀要把這張欠條兌現。在一個特別的形狀上——圓盤——我們會把解完全寫出來，化成一個你能徒手算的積分。

為什麼是圓盤，而不是正方形或一團不規則形狀？因為圓盤有一種對稱——旋轉——恰好與拉普拉斯方程自身的對稱相符，後者對一切方向一視同仁。正是這種相符，才讓分離變數法咬得住。在圓盤上，自然的座標是極座標 (r, theta)：從中心向外的半徑 r，與繞一圈的角度 theta。邊界於是再簡單不過——它就是圓 r = 1，由單一的角度 theta 參數化。規定邊界資料，就是交出一個關於邊緣角度的函數 f(theta)，而最終的獎賞，是一個把這個 f 直接化成內部溫度 u(r, theta) 的公式。

從傅立葉級數把答案建起來

引擎在此。在極座標下，拉普拉斯方程 u_xx + u_yy = 0 變成 u_rr + (1/r) u_r + (1/r^2) u_(theta theta) = 0。試一個乘積 u = R(r) T(theta)。分離把它拆成兩個常微分問題：T 必須以 2*pi 為週期（轉滿一整圈後你回到出發點），這迫使 T 是 cos(n*theta) 或 sin(n*theta)，其中 n 為整數；而 R 接著滿足一個簡單方程，其在中心保持有限的解是 R = r^n。因此圓盤上調和的基本積木是 r^n cos(n*theta) 與 r^n sin(n*theta)——而它們的任意和仍然是調和的。

現在做疊加。把內部的解寫成對所有這些積木的求和：u(r, theta) = a_0/2 + 對 n 求和 r^n (a_n cos(n*theta) + b_n sin(n*theta))。令 r = 1，它必須等於邊界資料 f(theta)。但這恰恰就是 f 的一個傅立葉級數：係數 a_n 與 b_n 不過就是這個邊緣函數的傅立葉係數。於是原則上配方已清楚——把 f 沿邊緣做傅立葉分析，再給第 n 個模態掛上一個因子 r^n，就完成了。因子 r^n 道盡了內部的一切：因為 0 <= r < 1，每個越高的模態都乘上越來越小的 r 的冪，所以邊界資料中那些抖動的高頻細節，你越往中心深入，被壓抑得越強。

把級數求和成一個核

傅立葉級數是個漂亮的配方，卻是個笨重的對象——你得算出無窮多個係數，再把無窮多項加起來。圓盤的魔法在於：這個無窮和可以一勞永逸地、以封閉形式、對所有資料同時完成。把 a_n 與 b_n 的積分公式代回級數，交換求和與積分的次序，那個看似幾何級數的 r^n 之和便塌縮了。掉出來的，就是 帕松積分公式：一個把邊界函數 f 拿進去、把內部值 u(r, theta) 交回來的單一積分，再沒有級數要加。

u(r, theta) = (1 / (2*pi)) * integral_0^(2*pi) P(r, theta - phi) f(phi) d phi

Poisson kernel (unit disk, radius 1):

                 1 - r^2
P(r, psi) = -----------------------------
            1 - 2 r cos(psi) + r^2

interior value  =  weighted average of the rim values f(phi),
the weight P leaning toward the nearest boundary point.

單位圓盤上的帕松積分公式。核 P(r, psi) 是權重；psi = theta - phi 是內部點的方向與邊緣點之間的夾角。

把這公式讀成一個加權平均，它便不再像代數。內部某點處的值 u(r, theta)，是所有邊緣值 f(phi) 的一個平均，但是有偏的：每個邊緣點以權重 P(r, theta - phi) 投票，這權重就是 帕松核。它永遠為正，且對離你位置最近的邊緣點最大、對正對面那個最小。所以你所站之處的溫度，主要是你肩頭那段邊界的回聲，遠側只有更微弱的貢獻——正是「最近的邊界最要緊」的直覺。分子裡的 1 - r^2 是那個旋鈕：在中心 r = 0 時它讓權重均勻，而當 r 攀向 1 時，權重越來越尖銳地堆到那個唯一最近的邊緣點上。

公式免費證明了什麼

帕松公式不只是個計算器；它一舉兌現了前幾篇導讀的每一張欠條。令 r = 0，即中心。在那裡核塌縮成 P(0, psi) = 1，一個平坦的權重，於是積分變成 f 沿整圈邊緣的純粹平均。這就是均值性質——中心值等於邊界平均——如今是被推導出來的，而非被假設的。核嚴格為正，意味著 u(r, theta) 是邊界值的一個道地的加權平均，所以它既不可能超過其中最大者、也不可能低於最小者：最大值原理就在此，直接從公式上讀出。

它也一併解決了圓盤狄利克雷問題的存在性，這是先前任何抽象論證都還沒做到的。交出任何合理的邊界函數 f——甚至一個不連續的，比方上半圓溫熱、下半圓冰冷——這個積分都會在內部產出一個道地的調和函數，而當你走向邊緣時它趨近於 f（在每個 f 連續之處）。存在性與唯一性合在一起，使這在圓盤上成為一個真正已解、適定的問題：解存在、是唯一的，且連續地依賴於資料，因為把 f 稍微推一下，它的加權平均也只被推一點點。

給定邊緣資料 f(theta)，大可完全略過係數——直奔積分。
選定你想求溫度的內部點 (r, theta)。
把帕松核 P(r, theta - phi) = (1 - r^2) / (1 - 2 r cos(theta - phi) + r^2) 寫成邊緣角度 phi 的函數。
用這個權重對 f 求平均：u = (1/(2*pi)) 對 phi 積分 P(r, theta - phi) f(phi)。那個數就是答案。

限制、誠實，與什麼能推廣

把這份禮物涵蓋與不涵蓋的，說清楚。第一，這公式解的是狄利克雷問題——規定邊界值——而非諾伊曼問題，後者規定的是邊界通量；那個問題需要它自己（密切相關）的核，以及你先前遇過的相容條件。第二，上面的封閉形式是圓盤的特權，用它的旋轉對稱換來。在正方形或不規則區域上並沒有這麼整齊的 P，你只能退回傅立葉型級數、數值格式，或格林函數——而帕松核其實正是格林函數的邊界跡，圓盤是那唯一一個其格林函數有顯式表達的形狀。

關於邊界，有一個微妙的誠實。內部是輝煌地平滑的，但這公式只在 f 連續的邊緣處復原 f。在一個跳躍處——溫熱與冰冷兩個半圓的接縫——內部的解平順地把差額對半分，當你趨近接縫時滑向兩個單側值的平均，從不真正等於那個跳躍本身。這不是缺陷；又是平滑化在告訴你：沒有調和函數能把一個真正的不連續帶進它的內部。而整個構造都仰賴資料住在區域的邊界上：若把資料規定在圓盤內部的一條曲線上、柯西式地給，那就會是本學習階梯適定性故事所警告的那種不適定設置。