傅立葉分析 · 傅立葉級數

傅立葉級數

傅立葉級數將週期函數寫成正弦與餘弦的和:a₀/2 + Σ (aₖ cos kx + bₖ sin kx)。每個係數都能單獨量出:將函數乘上一個波,再求有向面積。

從圖中看傅立葉級數

上方圖形以墨色畫出週期為2π2\pi的函數ff,以黃色畫出部分和SNS_N:常數項加上頻率不超過NN的全部餘弦與正弦。中間的圖將ff乘上一個探針波,並以陰影表示乘積曲線與橫軸之間的有向面積。最下方排列這樣量出的係數,餘弦在上,正弦在下。

每個係數都是一個面積
ak=1π∫−ππf(x)cos⁡kx dx,bk=1π∫−ππf(x)sin⁡kx dxa_k = \frac{1}{\pi}\int_{-\pi}^{\pi} f(x)\cos kx\,dx, \qquad b_k = \frac{1}{\pi}\int_{-\pi}^{\pi} f(x)\sin kx\,dx

ff同時包含所有波,為什麼一個面積能挑出一個係數?在完整週期上,1,cos⁡x,sin⁡x,cos⁡2x,sin⁡2x,…1, \cos x, \sin x, \cos 2x, \sin 2x, \dots中任意兩個不同的波相乘,有向面積都是0;正弦或餘弦乘上自己,面積則是π\pi。因此,將這些波的和乘上sin⁡kx\sin kx再積分,除了bksin⁡kxb_k \sin kx外,其餘各項都消去,只留下πbk\pi b_k。在兩個正弦組成的精選視角中,中間的圖將乘積拆成兩部分,另一個頻率的部分會正負抵消。這些波就像互相垂直的座標軸,各係數是對應軸上的投影,因此加進新的一項不會改變已量出的係數。

項數足夠多時,部分和會在跳躍旁過衝。隨著NN增加,峰會變窄並靠近跳躍,但過衝不會消失:方波超出頂端的高度,趨近跳躍幅度的約8.949%。每個SNS_N都連續,因此必須在一小段距離內跨過跳躍;在任何包含跳躍的區間上,最大誤差都不可能趨近0。單靠連續性無法決定過衝大小:8.949%來自部分和的計算,分段平滑函數的每個跳躍最後都有相同的比例。

越平滑的函數需要的項越少。這裡的ff指無限重複的函數,所以±π\pm\pi兩端接不起來,也算一個跳躍。分部積分將微分從波移到ff上,並帶出因子1/k1/k。遇到跳躍時,第一次分部積分就留下含因子1/k1/k的邊界項:方波在kk為奇數時有bk=4/(kπ)b_k = 4/(k\pi),在kk為偶數時則為0。若ff連續,只有斜率跳躍,例如三角波的尖角,第一次分部積分不留下該邊界項,第二次才留下,因此係數按1/k21/k^2衰減。切換對數刻度,比較長條與1/k1/k、1/k21/k^2兩條虛線。

能量讀數將同樣的垂直關係用於平方。f2f^2的積分可分解為係數的平方和,沒有交叉項:1π∫−ππf(x)2 dx=a022+∑k≥1(ak2+bk2).\frac{1}{\pi}\int_{-\pi}^{\pi} f(x)^2\,dx = \frac{a_0^2}{2} + \sum_{k \ge 1} \left(a_k^2 + b_k^2\right).讀數顯示SNS_N已占總量的多少;剩下的就是以相同方式量測的誤差平方。對方波,S1S_1約占81%,S9S_9約占96%,過衝卻仍存在:很窄的峰只增加少許面積。部分和可以在平均意義上很接近ff,卻在跳躍附近仍有明顯偏差。

值得留意

  • 常數項寫成a0/2a_0/2,是為了讓餘弦公式也適用於k=0k = 0。常數探針直接量出a0/2a_0/2,也就是ff的平均值。
  • 方波這類奇函數只有正弦項,∣sin⁡x∣|\sin x|這類偶函數只有餘弦項,因此頻譜總有一排是空的。

操作方式

  • 按一下頻譜中的長條,即可將探針移到該處。
  • 選擇「手繪」,在上方圖形中拖曳畫過一個週期。畫出的點以直線段連接,並向兩側重複。

相關內容