傅里葉級數

將周期函數分解為更簡單的正弦形式的和

數學中,傅里葉級數(英語:Fourier series/ˈfʊri, -iər/)是把類似的函數表示成簡單諧波的方式。更正式地說,對於滿足狄利克雷定理周期函數,其傅里葉級數是由一組正弦餘弦函數加權和表示的方法。傅里葉級數與用來找出無周期函數的頻率信息的傅里葉變換有密切的關係。

方波的傅里葉級數前四項的部份總和。隨着增加更多的諧波,這個部份總和收斂成方波。

傅里葉級數是傅里葉分析的一個研究分支,也是採樣定理原始證明的核心。傅里葉級數在數論組合數學信號處理、概率論統計學密碼學聲學光學等領域都有着廣泛的應用。

歷史

傅里葉級數得名於法國數學家約瑟夫·傅里葉(1768年–1830年),他提出任何函數都可以展開三角級數。此前數學家歐拉達朗貝爾克萊羅,已發現在認定一個函數有三角級數展開後,通過積分方法計算其係數的公式,而拉格朗日等人已經找到了一些非周期函數的三角級數展開。將周期函數分解為簡單振盪函數的總和的最早想法,可以追溯至公元前3世紀古代天文學家的均輪和本輪學說。

傅里葉的工作得到了丹尼爾·伯努利的贊助[1],傅里葉介入三角級數用來解熱傳導方程,其最初論文雖經西爾維斯特·拉克魯瓦英語Sylvestre François Lacroix加斯帕爾·蒙日同意[2],但在1807年經拉格朗日拉普拉斯勒讓德評審後被拒絕出版,他的現在被稱為傅里葉逆轉定理英語Fourier inversion theorem的理論後來發表於1822年出版的《熱的解析理論》[3]

定義

 
上圖展示了藍色無周期函數 ,它只定義在紅色區間 之上。可以認為傅里葉級數,是對下圖中最初函數的「周期延拓」的分析,傅里葉級數總是周期函數,即使最初函數 不是周期函數。

傅里葉級數可以用不同的形式來表達,下面將周期 的一個周期函數 表達為不同形式的傅里葉級數。

正弦-餘弦形式

人們常用  三角級數來表示 ,就是將所有 諧波  ,乘以其各自在 中的權重,求得它們的總和;這些 階諧波的權重稱爲傅立葉級數係數,它們可以藉由如下積分來獲得:

傅立葉級數係數
  Eq.1

符號 表示在選定區間上的積分,典型的選擇為 或者 。注意 是函數 平均值英語Mean of a function[A],這個性質擴展到了類似的變換比如傅里葉變換

通過這些係數定義傅里葉級數為:

傅立葉級數,正弦-餘弦形式
  Eq.2

這裡使用符號 ,表示傅里葉級數的求和不一定總是等於 。普遍來說 是理論上趨近於無限大的,但是就算趨近於無限大,對所有的 (例如在某一點上不連續),傅立葉級數也不一定收斂到  。儘管不收斂的可能性始終存在,在科學和工程領域中經常將Eq. 2中的 直接替代為 

在傅里葉級數係數中的整數索引 ,是級數中相應的  ,在這個函數的周期 中,形成的圓周(cycle)的數目。因此對應於  的項有着:

  • 波長等於 ,並且有着同於 的單位。
  • 頻率等於 ,並且有着 的倒數單位。

指數形式

下面藉由歐拉公式 ,將傅里葉級數係數簡化成複數指數形式。

根據定義,我們可以得到:

複數傅里葉級數係數
  Eq. 3

通過將等式Eq. 1代入Eq. 3,可以證實[4]

複數傅里葉級數係數

 

給定複數傅里葉級數係數,可以用公式復原出  

複數傅里葉級數係數

 

通過這些定義,傅里葉級數可以寫為:

傅立葉級數,指數形式
  Eq. 4

這是可推廣到複數值域函數的慣用形式。 的負數值對應於負頻率

複數值函數

人們習慣將 值域普遍化到複數上,設 是一個複數值函數,它的實部和虛部,都是實數值函數:

 

定義 則:

 
 
 

對於這個複數值函數,它的傅里葉級數的實部,是它的實部的傅里葉級數;它的傅里葉級數的虛部,是它的虛部的傅里葉級數:

 

振幅-相位形式

還可以利用三角恆等式 ,把正弦-餘弦形式中後面的正弦函數跟餘弦函數合併起來:

 

然後定義振幅 相位 ,這裡的  對應正弦-餘弦形式中    平均值英語Mean of a function 

傅立葉級數,振幅-相位形式
  Eq. 5

部份求和算子

在描述傅里葉級數行為的時候,經常會為一個函數 介入部份求和算子 [5]

  Eq. 6

這裡的  的傅里葉係數。不同於微積分中的級數,傅里葉級數的部份求和必須採用對稱形式,否則收斂結果可能不成立。

假設  是在 上的可積函數,   捲積 為:

 

周期為 的函數 的傅立葉級數的部份求和,可以經由 狄利克雷核 摺積來表示:

 

收斂性概要

  近似了 ,該近似程度會隨着 逐漸改善。這個無窮和 叫做  的傅里葉級數表示。傅里葉級數的收斂性取決於函數有限數量的極大值和極小值,這就是通常稱為傅里葉級數的狄利克雷條件。參見傅里葉級數的收斂性英語Convergence of Fourier series之一。對於廣義函數或分布也可以用範數或弱收斂英語Weak convergence (Hilbert space)定義傅里葉係數。在 不可導點上,如果我們只取無窮級數中的有限項求和,那麼在這些點上會有幅度不隨 增大而持續變小的起伏,這叫做吉布斯現象,一個簡單的例子是方波信號

工程應用中,一般假定傅里葉級數除了在不連續點以外處處收斂,原因是工程上遇到的函數比數學家提供的這個假定的反例表現更加良好。特別地,傅里葉級數絕對收斂一致收斂 ,只要在 的導數(或許不會處處存在)是平方可積的[6]。如果一個函數在區間 上是平方可積的,那麼此傅里葉級數在幾乎處處的點都收斂於該函數。

其他常用表示法

符號 在討論多個不同函數的傅里葉係數時是不夠用的。因此習慣上將其替代為函數(這裡是函數 )的某種修改形式,即採用函數式符號比如  ,來替代下標式符號:

 常用的數學符號
 常用的工程符號

在工程上,特別是在變量 表示時間的時候,係數序列叫做頻域表示。經常使用方括號來強調這個函數的定義域是頻率的離散集合。

另一個常用頻域表示,使用傅里葉級數係數,調製像梳子一樣的狄拉克採樣函數英語Dirac comb

 

這裡的 表示連續頻域。在變量 以秒為單位的時候, 赫茲為單位。採樣的間隔為基本頻率  倍(即為諧波)。  可以通過逆傅里葉變換英語Fourier inversion theorem從這種表示恢復出來:

 

構造出的函數 ,因而通常稱為「傅里葉變換」,即使一個周期函數的傅里葉積分在這個諧波頻率上不收斂[B]

常用的傅里葉級數

下表列出常用的周期函數及其傅里葉級數係數。

  •  指示周期 的周期函數。
  •    指示周期函數 的傅里葉級數係數(正弦-餘弦形式)。
時域
 
繪圖 頻域(正弦-餘弦形式)
 
注釋 引用
 
 
  全波整流正弦 [8]:p. 193
 
 
  半波整流正弦 [8]:p. 193
 
 
   
 
 
  鋸齒函數 [8]:p. 192
 
 
  反鋸齒函數 [8]:p. 192
 
 
  反全波整流 [8]:p. 193

基本性質

下表展示在時域中的一些數學運算及其對應的在傅里葉級數係數上的效果。

  • 複數共軛指示為上標星號 
  •   指示周期為 的函數或只定義在 中的函數。
  •   指示  的傅里葉級數係數(指數形式)。
性質 時域 頻域(指數形式) 注釋 引用
線性      
時間反轉/頻率反轉     [9]:p. 610
時間共軛     [9]:p. 610
時間反轉且共軛    
時間實部    
時間虛部    
頻率實部    
頻率虛部    
時間移位/頻率調製       [9]:p. 610
頻率移位/時間調製       [9]:p. 610

對稱性質

所有的函數都可以分解成唯一性的偶部和奇部英語Even–odd decomposition ,這裡的  。實數參數的複數值函數 ,對於所有 ,如果 則稱其為「偶對稱」,如果 則稱其為「奇對稱」,這裡 的上頂橫線指示複數共軛

一個複數值函數的實部和虛部,分解成各自的偶部和奇部,就有了四個分量,分別用下標標明為RE、RO、IE和IO。一個複數值時間參數函數的四個分量,與它的複數頻率變換的四個分量之間,有着一一映射[10]

 

由此可見,各種關係是顯而易見的,例如:

  • 實數值函數sRE + sRO的變換,是偶對稱函數SRE + i SIO。反過來說,偶對稱變換蘊含了實數值時域。
  • 虛數值函數i sIE + i sIO的變換,是奇對稱函數SRO + i SIE,反過來說也成立。
  • 偶對稱函數sRE + i sIO的變換,是實數值函數SRE + SRO,反過來說也成立。
  • 奇對稱函數sRO + i sIE的變換,是虛數值函數i SIE + i SIO,反過來說也成立。

範例

一個簡單的傅里葉級數

 
鋸齒波周期函數的圖
 
前五個部分傅里葉級數的動態圖

我們現在用上面的公式給出一個簡單函數的傅里葉級數展開式。考慮一個鋸齒波:

 
 

在這種情況下,傅里葉級數為:

 

可以證明,當 可微時,傅立葉級數在每個點 都收斂於 ,於是:

  Eq.1

 時,傅里葉級數收斂於 ,為在  的左極限和右極限之和的一半。這是傅里葉級數的狄利克雷定理的特例。

這個例子為我們引出了巴塞爾問題的一種解法。

傅里葉誘導

 
金屬板內的熱分布,使用傅里葉方法求解

在上例中我們的函數的傅里葉級數展開式看起來不比 簡單,因此人們需要傅里葉級數的原因也就不會立即顯現出來。但還有很多應用,我們舉用傅里葉誘導解熱方程的例子。考慮邊長為 米的方形金屬版,坐標為 。如果板內沒有熱源,並且四個邊中三個都保持在 攝氏度,而第四條邊 ,對於 ,保持在溫度梯度 攝氏度。在這種情況下,穩態(或者說很長時間過後的)熱分布函數 不能得出解析解,但卻可以證明:

 

這裡的 雙曲正弦函數。熱方程的這個解是通過將 的傅里葉級數的每一項乘以 得到的。儘管示例的函數 的傅里葉級數似乎很複雜,用傅里葉的方法卻可以求解這個熱分布問題。

其他例子

我們也可以應用傅立葉級數去證明等周不等式,或是構造處處連續處處不可微的函數。

收斂性

至今還沒有判斷傅里葉級數的收斂性充分必要條件,但是對於實際問題中出現的函數,有很多種判別條件可用於判斷收斂性。比如 的可微性或級數的一致收斂性。在閉區間上滿足狄利克雷條件的函數表示成的傅里葉級數都收斂。狄利克雷條件如下:

  1. 在定義區間上, 絕對可積
  2. 在任一有限區間中, 只能取有限個極值點;
  3. 在任何有限區間上, 只能有有限個第一類間斷點

滿足以上條件的 傅里葉級數都收斂,且:

1.當  的連續點時,級數收斂於 
2.當  的間斷點時,級數收斂於 

1966年,里納特·卡爾松證明了勒貝格二次可積函數的傅立葉級數一定是幾乎處處收斂的,即級數在除了一個勒貝格零測集外均收斂。

傅立葉級數收斂證明

假設一個函數在  上是平方可積,則會有:

  

證明的第一步:

考慮一系列正交基底, ,其中 ,且有

 

然後有 

特別的有, 的傅立葉級數的部分和 

然後根據  以及畢氏定理,可以有:

  替換一下後有  

如果右邊第一項收斂到0,再根據正交的性質,可以看出上述式子中的右手邊第二項:

 ,這就證明了帕塞瓦爾定理

證明的第二步:

回到證明右邊第一項,因為函數 可積,找到一個連續函數 ,然後根據最佳逼近引理,可以找到一個三角多項式p(x),使得

 

故當 ,函數  的差為0。

其他性質

傅立葉級數的唯一性

如果有一個定義在 的函數  ,其中函數  的傅立葉係數 還有 相同,且傅立葉級數都收斂到函數本身,那麼可以證明此傅立葉級數具有唯一性,也就是 。換句話說,如果函數  上可積,傅立葉係數 為0,對所有的 ,那麼函數 

卷積定理

給定周期為 的函數  ,它們具有傅里葉級數係數  ,這裡的 

  • 逐點乘積 ,也是周期為 ,並且它的傅里葉級數係數是序列  離散卷積 
  • 周期卷積 ,也是周期為 ,它具有傅里葉級數係數: 
  •  中的雙無限序列 ,是在 中的傅里葉係數的序列,當且僅當它是在 中的兩個序列的卷積[11]

微分性質

我們說 屬於在    如果 是一個在實數上以 為週期的函數,且 次可微而且 階連續。

  • 如果 屬於在 ,那麼 傅立葉係數 可以被用 傅立葉係數 的表示,藉由公式 
  • 如果 屬於在  。特別的,當固定 ,我們有 趨近於0當 ,且有 

黎曼-勒貝格定理

如果 可積函數,則   

帕塞瓦爾定理

如果函數 屬於在 之中,那麼便有 

普朗歇爾定理

如果 是係數,並且 ,則有一個唯一的函數 使得對於所有 有着 

延伸

希爾伯特空間的解讀

 
正弦和餘弦形成了正交集合。正弦、餘弦及其乘積的積分,當mn不同或二函數不同時是0(綠色和紅色區域相等抵消),僅當mn相等並且函數相同時為π。

所謂的兩個不同向量正交是指它們的內積為0,這也就意味着這兩個向量之間沒有任何相關性,例如,在三維歐氏空間中,互相垂直的向量之間是正交的。事實上,正交是垂直在數學上的一種抽象化和一般化。一組n個互相正交的向量必然是線性無關的,所以必然可以張成一個n維空間,也就是說,空間中的任何一個向量可以用它們來線性表出。

希爾伯特空間釋義下,函數的集合{en = einx; nZ}是[−π, π]平方可積函數L2([−π, π])的正交基。這個空間實際上是一個希爾伯特空間,有著針對任何兩個的元素fg的如下內積:

 

三角函數族的正交性用公式表示出來就是:

 
 

(這裡的δmn克羅內克函數),而

 

參閱

注釋

  1. ^ 一些作者定義了與此不同的 ,使得可以用相同的積分定義  。這改變了Eq. 2使得第一項需要除以 
  2. ^ 因為周期函數的傅里葉變換的積分定義不是收斂的,需要將周期函數和它的變換視為分布。在這種意義上, 是一個狄拉克δ函數,它是分布的是例子。

引用

  1. ^ 詳見莫里斯·克萊因《古今數學思想》,第20章無窮級數,第5節三角級數;第28章十九世紀的偏微分方程,第5節熱方程與傅里葉級數。
    see here, pg.s 209 & 210, 頁面存檔備份,存於網際網路檔案館
  2. ^ 李狗嗨. 如何给文科生解释傅里叶变换?. 知乎專欄. 2019-07-25 [2020-02-07]. (原始內容存檔於2020-10-24) (中文). 
  3. ^ Théorie analytique de la chaleur. Paris: Firmin Didot Père et Fils. 1822. OCLC 2688081 (法語). 
  4. ^ Pinkus, Allan; Zafrany, Samy. Fourier Series and Integral Transforms 1st. Cambridge, UK: Cambridge University Press. 1997: 42–44. ISBN 0-521-59771-4. 
  5. ^ Katznelson, Yitzhak. An introduction to Harmonic Analysis 2nd corrected. New York, NY: Dover Publications, Inc. 1976: 46. ISBN 0-486-63331-4. 
  6. ^ Georgi P. Tolstov. Fourier Series. Courier-Dover. 1976. ISBN 0-486-63317-9. 
  7. ^ 這裡頁面存檔備份,存於網際網路檔案館
  8. ^ 8.0 8.1 8.2 8.3 8.4 Papula, Lothar. Mathematische Formelsammlung: für Ingenieure und Naturwissenschaftler [Mathematical Functions for Engineers and Physicists]. Vieweg+Teubner Verlag. 2009. ISBN 978-3834807571 (德語). 
  9. ^ 9.0 9.1 9.2 9.3 Shmaliy, Y.S. Continuous-Time Signals. Springer. 2007. ISBN 978-1402062711. 
  10. ^ Proakis, John G.; Manolakis, Dimitris G. Digital Signal Processing: Principles, Algorithms, and Applications  3rd. Prentice Hall. 1996: 291. ISBN 978-0-13-373762-2. 
  11. ^ Characterizations of a linear subspace associated with Fourier series. MathOverflow. 2010-11-19 [2014-08-08]. 

延伸閱讀

外部連結

本條目含有來自PlanetMathexample of Fourier series》的內容,版權遵守知識共享協議:署名-相同方式共享協議