在數學 中,傅立葉級數 (英語:Fourier series , )是把類似波 的函數表示成簡單諧波 的方式。更正式地說,對於滿足狄利克雷定理 的週期函數 ,其傅立葉級數 是由一組正弦 與餘弦函數 的加權 和表示的方法。傅立葉級數與用來找出無週期函數的頻率資訊的傅立葉轉換 有密切的關係。
方波 的傅立葉級數前四項的部份總和。隨著增加更多的諧波,這個部份總和收斂成方波。
傅立葉級數是傅立葉分析 的一個研究分支,也是取樣定理 原始證明的核心。傅立葉級數在數論 、組合數學 、訊號 處理、機率論 、統計學 、密碼學 、聲學 、光學 等領域都有著廣泛的應用。
傅立葉級數得名於法國數學家約瑟夫·傅立葉 (1768年–1830年),他提出任何 函數都可以展開 為三角級數 。此前數學家歐拉 、達朗貝爾 和克萊羅 ,已發現在認定一個函數有三角級數展開後,通過積分方法計算其係數的公式,而拉格朗日 等人已經找到了一些非週期函數的三角級數展開。將週期函數分解為簡單振盪 函數的總和的最早想法,可以追溯至公元前3世紀古代天文學家的均輪和本輪 學說。
傅立葉的工作得到了丹尼爾·伯努利 的贊助[ 1] ,傅立葉介入三角級數用來解熱傳導方程式 ,其最初論文雖經西爾維斯特·拉克魯瓦 、加斯帕爾·蒙日 同意[ 2] ,但在1807年經拉格朗日 、拉普拉斯 和勒讓德 評審後被拒絕出版,他的現在被稱為傅里葉逆轉定理 的理論後來發表於1822年出版的《熱的解析理論》[ 3] 。
下面藉由歐拉公式
e
i
x
=
cos
x
+
i
sin
x
{\displaystyle \ e^{ix}=\cos x+i\sin x\ }
,將傅立葉級數係數簡化成複數 指數 形式。
根據定義,我們可以得到:
複數傅立葉級數係數
c
0
=
A
0
c
n
=
(
A
n
−
i
B
n
)
/
2
for
n
>
0
c
n
=
(
A
−
n
+
i
B
−
n
)
/
2
for
n
<
0
{\displaystyle {\begin{aligned}c_{0}&=A_{0}&\\c_{n}&=(A_{n}-iB_{n})/2\qquad &{\text{for }}n>0\\c_{n}&=(A_{-n}+iB_{-n})/2\qquad &{\text{for }}n<0\end{aligned}}}
Eq. 3
通過將等式Eq. 1 代入Eq. 3 ,可以證實[ 4] :
複數傅立葉級數係數
c
n
=
1
P
∫
P
s
(
x
)
e
−
2
π
i
n
x
P
d
x
for
n
∈
Z
{\displaystyle c_{n}={\frac {1}{P}}\int _{P}s(x)e^{-{\frac {2\pi inx}{P}}}\,dx\qquad {\text{for}}\ n\in \mathbb {Z} }
給定複數傅立葉級數係數,可以用公式復原出
A
n
{\displaystyle A_{n}}
和
B
n
{\displaystyle B_{n}}
:
複數傅立葉級數係數
A
0
=
c
0
A
n
=
c
n
+
c
−
n
for
n
>
0
B
n
=
i
(
c
n
−
c
−
n
)
for
n
>
0
{\displaystyle {\begin{aligned}A_{0}&=c_{0}&\\A_{n}&=c_{n}+c_{-n}\qquad &{\textrm {for}}~n>0\\B_{n}&=i(c_{n}-c_{-n})\qquad &{\textrm {for}}~n>0\end{aligned}}}
通過這些定義,傅立葉級數可以寫為:
傅立葉級數,指數形式
s
(
x
)
∼
∑
n
=
−
∞
∞
c
n
⋅
e
2
π
i
n
x
P
{\displaystyle s(x)\sim \sum _{n=-\infty }^{\infty }c_{n}\cdot e^{\frac {2\pi inx}{P}}}
Eq. 4
這是可推廣到複數值域 函數的慣用形式。
n
{\displaystyle n}
的負數值對應於負頻率 。
鋸齒波 週期函數的圖
前五個部分傅立葉級數的動態圖
我們現在用上面的公式給出一個簡單函數的傅立葉級數展開式。考慮一個鋸齒波:
s
(
x
)
=
x
π
,
f
o
r
−
π
<
x
<
π
{\displaystyle s(x)={\frac {x}{\pi }},\quad \mathrm {for} -\pi <x<\pi }
s
(
x
+
2
π
k
)
=
s
(
x
)
,
f
o
r
−
∞
<
x
<
∞
and
k
∈
Z
{\displaystyle s(x+2\pi k)=s(x),\quad \mathrm {for} -\infty <x<\infty {\text{ and }}k\in \mathbb {Z} }
在這種情況下,傅立葉級數為:
A
n
=
1
π
∫
−
π
π
s
(
x
)
cos
(
n
x
)
d
x
=
0
,
n
≥
0
B
n
=
1
π
∫
−
π
π
s
(
x
)
sin
(
n
x
)
d
x
=
−
2
π
n
cos
(
n
π
)
+
2
π
2
n
2
sin
(
n
π
)
=
2
(
−
1
)
n
+
1
π
n
,
n
≥
1
{\displaystyle {\begin{aligned}A_{n}&={\frac {1}{\pi }}\int _{-\pi }^{\pi }s(x)\cos(nx)\,dx=0,\quad n\geq 0\\[4pt]B_{n}&={\frac {1}{\pi }}\int _{-\pi }^{\pi }s(x)\sin(nx)\,dx\\[4pt]&=-{\frac {2}{\pi n}}\cos(n\pi )+{\frac {2}{\pi ^{2}n^{2}}}\sin(n\pi )\\[4pt]&={\frac {2\,(-1)^{n+1}}{\pi n}},\quad n\geq 1\end{aligned}}}
可以證明,當
s
{\displaystyle s}
可微時,傅立葉級數在每個點
x
{\displaystyle x}
都收斂於
s
(
x
)
{\displaystyle s(x)}
,於是:
s
(
x
)
=
A
0
+
∑
n
=
1
∞
(
A
n
cos
(
n
x
)
+
B
n
sin
(
n
x
)
)
=
2
π
∑
n
=
1
∞
(
−
1
)
n
+
1
n
sin
(
n
x
)
,
f
o
r
x
−
π
∉
2
π
Z
{\displaystyle {\begin{aligned}s(x)&=A_{0}+\sum _{n=1}^{\infty }\left(A_{n}\cos \left(nx\right)+B_{n}\sin \left(nx\right)\right)\\[4pt]&={\frac {2}{\pi }}\sum _{n=1}^{\infty }{\frac {(-1)^{n+1}}{n}}\sin(nx),\quad \mathrm {for} \quad x-\pi \notin 2\pi \mathbb {Z} \end{aligned}}}
Eq.1
當
x
=
π
{\displaystyle x=\pi }
時,傅立葉級數收斂於
0
{\displaystyle 0}
,為在
x
=
π
{\displaystyle x=\pi }
處
s
{\displaystyle s}
的左極限和右極限之和的一半。這是傅立葉級數的狄利克雷定理 的特例。
這個例子為我們引出了巴塞爾問題 的一種解法。
金屬板內的熱分布,使用傅立葉方法求解
在上例中我們的函數的傅立葉級數展開式看起來不比
s
(
x
)
=
x
π
{\displaystyle s(x)={\tfrac {x}{\pi }}}
簡單,因此人們需要傅立葉級數的原因也就不會立即顯現出來。但還有很多應用,我們舉用傅立葉誘導解熱方程式 的例子。考慮邊長為
π
{\displaystyle \pi }
米的方形金屬版,坐標為
(
x
,
y
)
∈
[
0
,
π
]
×
[
0
,
π
]
{\displaystyle (x,y)\in [0,\pi ]\times [0,\pi ]}
。如果板內沒有熱源,並且四個邊中三個都保持在
0
{\displaystyle 0}
攝氏度,而第四條邊
y
=
π
{\displaystyle y=\pi }
,對於
x
∈
(
0
,
π
)
{\displaystyle x\in (0,\pi )}
,保持在溫度梯度
T
(
x
,
π
)
=
x
{\displaystyle T(x,\pi )=x}
攝氏度。在這種情況下,穩態(或者說很長時間過後的)熱分布函數
T
(
x
,
y
)
{\displaystyle T(x,y)}
不能得出解析解 ,但卻可以證明:
T
(
x
,
y
)
=
2
∑
n
=
1
∞
(
−
1
)
n
+
1
n
sin
(
n
x
)
sinh
(
n
y
)
sinh
(
n
π
)
{\displaystyle T(x,y)=2\sum _{n=1}^{\infty }{\frac {(-1)^{n+1}}{n}}\sin(nx){\sinh(ny) \over \sinh(n\pi )}}
這裡的
sinh
{\displaystyle \sinh }
是雙曲正弦 函數。熱方程式的這個解是通過將
π
s
(
x
)
{\displaystyle \pi s(x)}
的傅立葉級數的每一項乘以
sinh
(
n
y
)
sinh
(
n
π
)
{\displaystyle {\tfrac {\sinh(ny)}{\sinh(n\pi )}}}
得到的。儘管示例的函數
s
(
x
)
{\displaystyle s(x)}
的傅立葉級數似乎很複雜,用傅立葉的方法卻可以求解這個熱分布問題。
我們也可以應用傅立葉級數去證明等周不等式 ,或是構造處處連續處處不可微的函數。
假設一個函數在
f
(
x
)
{\displaystyle f(x)}
在
[
0
,
2
π
]
{\displaystyle [0,2\pi ]}
上是平方可積,則會有:
1
2
π
∫
0
2
π
|
f
(
x
)
−
S
N
(
f
)
(
x
)
|
2
d
x
→
0
{\displaystyle {\frac {1}{2\pi }}\int _{0}^{2\pi }|f(x)-S_{N}(f)(x)|^{2}dx\rightarrow 0}
當
N
→
∞
{\displaystyle N\rightarrow \infty }
證明的第一步:
考慮一系列正交基底,
{
e
n
}
n
∈
Z
{\displaystyle \{e_{n}\}_{n\in \mathbb {Z} }}
,其中
e
n
(
x
)
=
e
−
i
n
x
{\displaystyle e_{n}(x)=e^{-inx}}
,且有
(
e
n
,
e
m
)
=
{
1
,
if
n
=
m
0
,
if
n
≠
m
{\displaystyle (e_{n},e_{m})={\begin{cases}1,&{\text{if }}n=m\\0,&{\text{if }}n\neq m\end{cases}}}
然後有
(
f
,
e
n
)
=
1
2
π
∫
0
2
π
f
(
x
)
e
−
i
n
x
d
x
=
f
^
(
n
)
{\displaystyle (f,e_{n})={\frac {1}{2\pi }}\int _{0}^{2\pi }f(x)e^{-inx}dx={\hat {f}}(n)}
特別的有,
f
(
x
)
{\displaystyle f(x)}
的傅立葉級數的部分和
S
N
(
f
)
(
x
)
=
∑
|
n
|
≤
N
f
^
(
n
)
e
n
{\displaystyle S_{N}(f)(x)=\sum _{|n|\leq N}{\hat {f}}(n)e_{n}}
然後根據
f
=
f
−
∑
|
n
|
≤
N
f
^
(
n
)
e
n
+
∑
|
n
|
≤
N
f
^
(
n
)
e
n
{\displaystyle f=f-\sum _{|n|\leq N}{\hat {f}}(n)e_{n}+\sum _{|n|\leq N}{\hat {f}}(n)e_{n}}
以及畢氏定理,可以有:
|
|
f
|
|
2
=
|
|
f
−
∑
|
n
|
≤
N
f
^
(
n
)
e
n
|
|
2
+
|
|
∑
|
n
|
≤
N
f
^
(
n
)
e
n
|
|
2
{\displaystyle ||f||^{2}=||f-\sum _{|n|\leq N}{\hat {f}}(n)e_{n}||^{2}+||\sum _{|n|\leq N}{\hat {f}}(n)e_{n}||^{2}}
替換一下後有
|
|
f
|
|
2
=
|
|
f
−
S
N
(
f
)
(
x
)
|
|
2
+
|
|
∑
|
n
|
≤
N
f
^
(
n
)
e
n
|
|
2
{\displaystyle ||f||^{2}=||f-S_{N}(f)(x)||^{2}+||\sum _{|n|\leq N}{\hat {f}}(n)e_{n}||^{2}}
如果右邊第一項收斂到0,再根據正交的性質,可以看出上述式子中的右手邊第二項:
|
|
∑
|
n
|
≤
N
f
^
(
n
)
e
n
|
|
2
=
∑
|
n
|
≤
N
|
f
^
(
n
)
|
2
{\displaystyle ||\sum _{|n|\leq N}{\hat {f}}(n)e_{n}||^{2}=\sum _{|n|\leq N}|{\hat {f}}(n)|^{2}}
,這就證明了帕塞瓦爾定理 。
證明的第二步:
回到證明右邊第一項,因為函數
f
(
x
)
{\displaystyle f(x)}
可積,找到一個連續函數
g
(
x
)
{\displaystyle g(x)}
,然後根據最佳逼近引理,可以找到一個三角多項式p(x),使得
|
f
−
S
N
(
f
)
(
x
)
|
≤
|
f
(
x
)
−
g
(
x
)
|
+
|
g
(
x
)
−
S
N
(
f
)
(
x
)
|
{\displaystyle |f-S_{N}(f)(x)|\leq |f(x)-g(x)|+|g(x)-S_{N}(f)(x)|}
故當
N
→
∞
{\displaystyle N\rightarrow \infty }
,函數
f
(
x
)
{\displaystyle f(x)}
跟
S
N
(
f
)
(
x
)
{\displaystyle S_{N}(f)(x)}
的差為0。
正弦和餘弦形成了正交集合。正弦、餘弦及其乘積的積分,當m 與n 不同或二函數不同時是0(綠色和紅色區域相等抵消),唯若m 和n 相等並且函數相同時為π。
所謂的兩個不同向量 正交是指它們的內積 為0,這也就意味著這兩個向量之間沒有任何相關性,例如,在三維歐氏空間中,互相垂直的向量之間是正交的。事實上,正交是垂直在數學上的一種抽象化和一般化。一組n個互相正交的向量必然是線性獨立 的,所以必然可以張成一個n維空間,也就是說,空間中的任何一個向量可以用它們來線性表出。
在希爾伯特空間 釋義下,函數的集合{en = einx ; n ∈ Z }是[−π, π]平方可積函數L 2 ([−π, π])的正交基 。這個空間實際上是一個希爾伯特空間,有著針對任何兩個的元素f 和g 的如下內積:
⟨
f
,
g
⟩
=
d
e
f
1
2
π
∫
−
π
π
f
(
x
)
g
(
x
)
¯
d
x
.
{\displaystyle \langle f,\,g\rangle \;{\stackrel {\mathrm {def} }{=}}\;{\frac {1}{2\pi }}\int _{-\pi }^{\pi }f(x){\overline {g(x)}}\,dx.}
三角函數族的正交性用公式表示出來就是:
∫
−
π
π
cos
(
m
x
)
cos
(
n
x
)
d
x
=
π
δ
m
n
,
m
,
n
≥
1
,
{\displaystyle \int _{-\pi }^{\pi }\cos(mx)\,\cos(nx)\,dx=\pi \delta _{mn},\quad m,n\geq 1,\,}
∫
−
π
π
sin
(
m
x
)
sin
(
n
x
)
d
x
=
π
δ
m
n
,
m
,
n
≥
1
{\displaystyle \int _{-\pi }^{\pi }\sin(mx)\,\sin(nx)\,dx=\pi \delta _{mn},\quad m,n\geq 1}
(這裡的δmn 是克羅內克函數 ),而
∫
−
π
π
cos
(
m
x
)
sin
(
n
x
)
d
x
=
0
;
{\displaystyle \int _{-\pi }^{\pi }\cos(mx)\,\sin(nx)\,dx=0;\,}
一些作者定義了與此不同的
A
0
{\displaystyle A_{0}}
,使得可以用相同的積分定義
A
0
{\displaystyle A_{0}}
和
A
n
{\displaystyle A_{n}}
。這改變了Eq. 2 使得第一項需要除以
2
{\displaystyle 2}
。
因為週期函數的傅立葉轉換的積分定義不是收斂的,需要將週期函數和它的轉換視為分布 。在這種意義上,
F
{
e
i
2
π
n
x
P
}
{\displaystyle {\mathcal {F}}\{e^{i{\frac {2\pi nx}{P}}}\}}
是一個狄拉克δ函數 ,它是分布的是例子。
Papula, Lothar. Mathematische Formelsammlung: für Ingenieure und Naturwissenschaftler [Mathematical Functions for Engineers and Physicists]. Vieweg+Teubner Verlag. 2009. ISBN 978-3834807571 (德語) .
William E. Boyce; Richard C. DiPrima. Elementary Differential Equations and Boundary Value Problems 8th. New Jersey: John Wiley & Sons, Inc. 2005. ISBN 0-471-43338-1 .
Joseph Fourier, translated by Alexander Freeman. The Analytical Theory of Heat. Dover Publications. 2003. ISBN 0-486-49531-0 . 2003 unabridged republication of the 1878 English translation by Alexander Freeman of Fourier's work Théorie Analytique de la Chaleur , originally published in 1822.
Enrique A. Gonzalez-Velasco. Connections in Mathematical Analysis: The Case of Fourier Series . American Mathematical Monthly. 1992, 99 (5): 427–441. JSTOR 2325087 . doi:10.2307/2325087 .
Fetter, Alexander L.; Walecka, John Dirk. Theoretical Mechanics of Particles and Continua . Courier. 2003. ISBN 978-0-486-43261-8 .
Felix Klein , Development of mathematics in the 19th century . Mathsci Press Brookline, Mass, 1979. Translated by M. Ackerman from Vorlesungen über die Entwicklung der Mathematik im 19 Jahrhundert , Springer, Berlin, 1928.
Walter Rudin . Principles of mathematical analysis 3rd. New York: McGraw-Hill, Inc. 1976. ISBN 0-07-054235-X .
A. Zygmund . Trigonometric Series third. Cambridge: Cambridge University Press. 2002. ISBN 0-521-89053-5 . The first edition was published in 1935.
電機電子類科《工程數學》,ISBN 978-957-584-377-9 ,作者 陳錫冠、曾致煌,高立出版社。