JOVANA
Explore Library Glossary Getting Started Three Levels Fields How it works Mission
Join the mission
All guides

傅立葉正弦與餘弦級數

分離變數法交給我們一組正弦或餘弦模態作為基底;現在我們要問它們能搭出什麼。傅立葉正弦或餘弦級數把區間上任何合理的形狀重寫成那些模態之和——而你該伸手拿哪一種,完全由邊界條件決定。

為什麼我們需要一個級數

上一篇指南跑完了分離變數法,並以一個懸念作結。代入乘積擬設 u(x,t) = X(x)T(t) 把偏微分方程拆成一個空間本徵值問題與一條時間常微分方程,由分離常數綁在一起。邊界條件接著把被允許的形狀篩成一個離散的家族——對冰水兩端而言,就是正弦模態 sin(n pi x / L)。但單一個模態幾乎永遠匹配不上你實際拿到的那個雜亂初始分布 f(x)。於是自然的問題變成:區間上一個任意的形狀,能不能被寫成這些模態之和?

令人驚訝地,可以——而這個肯定的答案正是傅立葉級數的全部要點。其主張是:[0, L] 上幾乎任何函數 f(x),都能藉由把正弦波(或在一個姊妹版本裡,餘弦波)以精心選定的振幅加總起來而重現。疊加原理保證:一旦 f 被這樣寫出來,把每個模態乘上它自己的時間因子再重新求和,就給出偏微分方程一個真正的解。這就是從積木通往真正答案的橋;搭這座橋正是本篇與接下來兩篇指南的工作。

正弦級數:被釘住的兩端

從你已經見過的情形開始。當兩端都被維持在零——齊次的狄利克雷條件 u(0) = u(L) = 0——本徵值問題逼使餘弦那一塊消亡,只留下正弦。於是我們把 f 展開成一個傅立葉正弦級數:f(x) = b_1 sin(pi x / L) + b_2 sin(2 pi x / L) + b_3 sin(3 pi x / L) + ...。每一項在 x = 0 與 x = L 處都自動歸零,因為清單裡每個正弦本身在那裡就是零。這個級數只能表示在兩端都歸零的函數——而那恰好就是邊界條件所允許的那個家族。

想像這些模態本身。第一個,sin(pi x / L),是一個從零升起又回到零的單一平滑隆起——半個波長橫跨整根桿子。第二個,sin(2 pi x / L),先上升、在中間穿過零、再下沉——一個整波。第 n 個模態有 n 個隆起與 n - 1 個內部過零點。以振幅 b_n 把它們相加,就像在一台圖示等化器上調整一排推桿:低階的 n 定下寬闊的輪廓,高階的 n 雕刻出細緻、皺褶的細節。有了夠多的推桿,你幾乎能塑造出任何從地板起、又回到地板的曲線。

餘弦級數:隔熱的兩端

現在改變兩端的物理。把它們隔熱起來,而不是浸冰:沒有熱量跨越流動,所以溫度的斜率在那裡歸零。這是齊次的諾伊曼條件,u_x(0) = u_x(L) = 0。在這些條件下重跑分離變數法,本徵值問題會保留餘弦而殺掉正弦——因為餘弦正是那個在端點上導數為零的函數。形狀變成 cos(n pi x / L),而展開式是一個傅立葉餘弦級數:f(x) = a_0/2 + a_1 cos(pi x / L) + a_2 cos(2 pi x / L) + ...。

餘弦級數的兩個特徵值得再看一眼。第一,它帶著一個正弦級數所沒有的額外項:常數 a_0/2,也就是 n = 0 的模態 cos(0) = 1。那個孤零零的常數是 f 在區間上的平均值,而它在熱問題裡永不衰減——隔熱的兩端把總熱量困住,所以桿子安頓到它的平均溫度,而不是零。第二,因為餘弦只在遠離兩端處才有非零斜率,這個級數自然地表示那些以平坦姿態抵達兩面牆的形狀,而這正是無通量邊界所要求的。

boundary at the two ends      ->  which modes survive  ->  series

  u = 0          (Dirichlet)   ->  sin(n pi x / L)      ->  sine series
  u_x = 0        (Neumann)     ->  cos(n pi x / L)      ->  cosine series
  u, u_x periodic              ->  sin & cos together   ->  full Fourier series
決定你該用哪一族模態來展開的,是邊界條件,而不是個人偏好。

一張秘密的圖:把函數延拓

為什麼邊界能這麼乾淨地在正弦與餘弦之間挑選?有一張單一的圖能同時解釋兩者。那些正弦對 x = 0 而言全是函數;那些餘弦全是函數。所以一個正弦級數暗地裡把 f 當作往負側做奇延拓那樣對待,而一個餘弦級數則把它做偶延拓。這就是偶延拓與奇延拓,也是這兩個級數底下隱藏的機械裝置。

把奇延拓的邏輯走到底。把 f 對 x = 0 做奇反射以填滿 [-L, 0],然後沿著軸把整塊 [-L, L] 反覆鋪排——這是一個週期為 2L 的週期延拓。奇反射逼使函數在 x = 0 處(以及經由鋪排在 x = L 處)穿過零,這恰好就是免費出現的狄利克雷條件。偶反射則以平坦——零斜率——的姿態抵達那同樣的兩點,免費把諾伊曼條件交給你。你施加的邊界,不過就是你所選那個延拓的對稱性。

把振幅讀出來

我們承諾了一個級數,卻還沒說怎麼找出振幅 b_n 或 a_n。那個訣竅——在緊接著的下一篇指南裡完整發展——就是正交性:不同的正弦模態在 [0, L] 上彼此積分,會抵銷成零,不同的餘弦模態也是如此。這意味著你能一次孤立出一個振幅。把 f 的級數乘上一個選定的模態,在區間上積分,除了相配的那一項之外,每一項都崩塌掉,留下那單一個係數的一條乾淨公式。

  1. 讀邊界。 兩端值為零意味著正弦(正弦級數);斜率為零意味著餘弦(餘弦級數)。
  2. 寫下展開式。 把相配的模態以未知振幅相加——b_n sin(n pi x / L),或 a_0/2 加上 a_n cos(n pi x / L)。
  3. 投影出每個振幅。 利用正交性,b_n = (2/L) 乘以 f(x) sin(n pi x / L) 在 [0,L] 上的積分,餘弦的 a_n 同理。
  4. 重新組裝。 把每個模態的時間因子接上去再重新求和,得到 u(x,t)——這就是匹配初始條件的動作。

在你盲目信任這個級數之前,有兩個誠實的提醒。第一,f 與它的級數相等並非處處逐點免費成立:在一個跳躍處,級數收斂到那道缺口的中點,而不是任一側,且跳躍附近的過衝永遠不會完全縮小。第二,正弦與餘弦並不是故事的終點——它們是某個特別簡單的幾何的本徵函數。在圓盤或球面上,同樣的邏輯產生的是貝索函數或斯圖姆-劉維爾本徵函數,而本階的最後一篇指南會展示那個更寬廣的世界。正弦與餘弦級數,是一個遠為一般的想法之中最乾淨的範例。