文章首发于:https://zhuanlan.zhihu.com/p/342982724
主要参考自 DR_CAN 的B站教程 纯干货傅里叶变换,某些地方加了一点自己的理解。。
0x01 三角函数系的正交性与傅里叶级数
我们首先看一下这个定理: 组成三角级数的函数系
0
(
sin
0
x
)
,
1
(
cos
0
x
)
,
sin
x
,
cos
x
,
sin
2
x
,
cos
2
x
,
…
,
sin
n
x
,
cos
n
x
\begin{aligned} 0(\text{sin}0x),1(\text{cos}0x),\text{sin}x,\text{cos}x,\text{sin}2x,\text{cos}2x,\dots,\text{sin}nx,\text{cos}nx \end{aligned}
0(sin0x),1(cos0x),sinx,cosx,sin2x,cos2x,…,sinnx,cosnx
在
[
−
π
,
π
]
[-\pi,\pi]
[−π,π] 上正交,即其中任意两个不同的函数之积在
[
−
π
,
π
]
[-\pi,\pi]
[−π,π] 上的积分等于0,证明见附录1。
我们可以先记住这个结论,就是从上面的那组三角函数中随便选两个不同的三角函数,将他们的乘起来,然后在 [ − π , π ] [-\pi,\pi] [−π,π] 下做积分,结果一定是0。
并且在三角函数系中两个相同的函数的乘积在
[
−
π
,
π
]
[-\pi,\pi]
[−π,π] 上的积分不等于0,且有:
∫
−
π
π
cos
n
x
⋅
cos
m
x
dx
=
∫
−
π
π
cos
2
n
x
dx
=
∫
−
π
π
1
2
(
1
+
cos
2
n
x
)
dx
=
1
2
[
∫
−
π
π
1
dx
+
∫
−
π
π
cos
2
n
x
dx
]
=
π
\begin{aligned} \int_{-\pi}^{\pi} \text{cos}nx·\text{cos}mx\text{dx} &= \int_{-\pi}^{\pi} \text{cos}^2nx\text{dx} \\ &=\int_{-\pi}^{\pi} \frac{1}{2}\left(1 +\text{cos}2nx\right)\text{dx} \\ &=\frac{1}{2} \left[\int_{-\pi}^{\pi}1\text{dx} + \int_{-\pi}^{\pi}\text{cos}2nx\text{dx}\right] \\ &=\pi \\ \end{aligned}
∫−ππcosnx⋅cosmxdx=∫−ππcos2nxdx=∫−ππ21(1+cos2nx)dx=21[∫−ππ1dx+∫−ππcos2nxdx]=π
∫
−
π
π
sin
2
n
x
dx
=
π
∫
−
π
π
1
⋅
1
dx
=
2
π
\begin{aligned} \int_{-\pi}^{\pi} \text{sin}^2nx\text{dx} =\pi \\ \int_{-\pi}^{\pi} 1·1\text{dx}=2\pi \end{aligned}
∫−ππsin2nxdx=π∫−ππ1⋅1dx=2π
我们知道简单的周期运动可写作: y = A sin ( ω t + φ ) y=A\text{sin}(\omega t+\varphi) y=Asin(ωt+φ) ,其中 A A A 为振幅, ω \omega ω 为角频率, φ \varphi φ 为初相。下图为 sin ( x + π ) \sin\left(x+\pi\right) sin(x+π):
复杂的周期运动 f ( t ) f(t) f(t) 可以看作简单周期运动的叠加。
下图为: sin ( x + π ) + 1 3 sin ( 3 x + π ) + 1 5 sin ( 5 x + π ) + 1 7 sin ( 7 x + π ) + 1 9 sin ( 9 x + π ) \sin\left(x+\pi\right)+\frac{1}{3}\sin\left(3x+\pi\right)+\frac{1}{5}\sin\left(5x+\pi\right)+\frac{1}{7}\sin\left(7x+\pi\right)+\frac{1}{9}\sin\left(9x+\pi\right) sin(x+π)+31sin(3x+π)+51sin(5x+π)+71sin(7x+π)+91sin(9x+π)
当累加的足够多,也就形成了一个周期性的方波。
因此周期函数
f
(
t
)
f(t)
f(t) 的表达式可写作:
f
(
t
)
=
∑
n
=
0
∞
A
n
sin
(
n
ω
t
+
φ
n
)
=
A
0
sin
φ
0
+
∑
n
=
1
∞
A
n
sin
(
n
ω
t
+
φ
n
)
=
A
0
sin
φ
0
+
∑
n
=
1
∞
[
A
n
sin
n
ω
t
cos
φ
n
+
A
n
cos
n
ω
t
sin
φ
n
]
\begin{aligned} f(t)&=\sum_{n=0}^{\infty}A_n\text{sin}{(n\omega t+\varphi_n)} \\ &=A_0\text{sin}{\varphi_0}+\sum_{n=1}^{\infty}A_n\text{sin}{(n\omega t+\varphi_n)} \\ &=A_0\text{sin}{\varphi_0}+\sum_{n=1}^{\infty}\left[A_n\text{sin}{n\omega t}\text{cos}{\varphi_n} + A_n\text{cos}{n\omega t}\text{sin}{\varphi_n}\right] \end{aligned}
f(t)=n=0∑∞Ansin(nωt+φn)=A0sinφ0+n=1∑∞Ansin(nωt+φn)=A0sinφ0+n=1∑∞[Ansinnωtcosφn+Ancosnωtsinφn]
我们令
a
n
=
A
n
sin
φ
n
a_n=A_n\text{sin}{\varphi_n}
an=Ansinφn ,
b
n
=
A
n
cos
φ
n
b_n=A_n\text{cos}{\varphi_n}
bn=Ancosφn,
ω
t
=
x
\omega t=x
ωt=x(这里我们**假设周期为
2
π
2\pi
2π,**即
ω
=
1
\omega=1
ω=1),代换后的周期函数表达式如下:
f
(
x
)
=
A
0
sin
φ
0
+
∑
n
=
1
∞
A
n
sin
(
n
ω
t
+
φ
n
)
=
a
0
+
∑
n
=
1
∞
(
a
n
cos
n
x
+
b
n
sin
n
x
)
(1)
\begin{aligned} f(x)&=A_0\text{sin}{\varphi_0}+\sum_{n=1}^{\infty}A_n\text{sin}{(n\omega t+\varphi_n)} \\ &=a_0+\sum_{n=1}^{\infty}(a_n\text{cos}{nx}+b_n\text{sin}{nx}) \\ \end{aligned} \tag{1}
f(x)=A0sinφ0+n=1∑∞Ansin(nωt+φn)=a0+n=1∑∞(ancosnx+bnsinnx)(1)
实际上,(1)式只是一个猜想式,我们需要找到
a
0
a_0
a0、
a
n
a_n
an 和
b
n
b_n
bn 式该表达式成立。
首先求
a
n
a_n
an的表达式,可将(1)式两边同时乘以
c
o
s
m
x
cosmx
cosmx,然后两边同时在
[
−
π
,
π
]
[-\pi,\pi]
[−π,π] 做积分,得到:
∫
−
π
π
f
(
x
)
cos
m
x
dx
=
∫
−
π
π
a
0
cos
m
x
dx
+
∫
−
π
π
∑
n
=
1
∞
a
n
cos
n
x
cos
m
x
dx
+
∫
−
π
π
∑
n
=
1
∞
b
n
sin
n
x
cos
m
x
dx
=
a
0
∫
−
π
π
cos
m
x
dx
+
a
n
∫
−
π
π
∑
n
=
1
∞
cos
n
x
cos
m
x
dx
+
b
n
∫
−
π
π
∑
n
=
1
∞
sin
n
x
cos
m
x
dx
(2)
\begin{aligned} \int_{-\pi}^{\pi}f(x)\text{cos}mx\text{dx}&=\int_{-\pi}^{\pi}a_0\text{cos}mx\text{dx}+\int_{-\pi}^{\pi}\sum_{n=1}^{\infty}a_n\text{cos}{nx}\text{cos}mx\text{dx}+\int_{-\pi}^{\pi}\sum_{n=1}^{\infty}b_n\text{sin}{nx}\text{cos}mx\text{dx} \\ &=a_0\int_{-\pi}^{\pi}\text{cos}mx\text{dx}+a_n\int_{-\pi}^{\pi}\sum_{n=1}^{\infty}\text{cos}{nx}\text{cos}mx\text{dx}+b_n\int_{-\pi}^{\pi}\sum_{n=1}^{\infty}\text{sin}{nx}\text{cos}mx \text{dx} \end{aligned} \tag{2}
∫−ππf(x)cosmxdx=∫−ππa0cosmxdx+∫−ππn=1∑∞ancosnxcosmxdx+∫−ππn=1∑∞bnsinnxcosmxdx=a0∫−ππcosmxdx+an∫−ππn=1∑∞cosnxcosmxdx+bn∫−ππn=1∑∞sinnxcosmxdx(2)
由前面证明的三角函数正交性可得:
∫
−
π
π
cos
m
x
dx
=
∫
−
π
π
cos
m
x
sin
0
x
dx
=
0
∫
−
π
π
∑
n
=
1
∞
sin
n
x
cos
m
x
dx
=
0
∫
−
π
π
∑
n
=
1
∞
cos
n
x
cos
m
x
dx
=
0
(
n
≠
m
)
∫
−
π
π
∑
n
=
1
∞
cos
n
x
cos
m
x
dx
=
π
(
n
=
m
)
\begin{aligned} &\int_{-\pi}^{\pi}\text{cos}mx\text{dx}=\int_{-\pi}^{\pi}\text{cos}mx\text{sin}0x\text{dx}=0 \\ &\int_{-\pi}^{\pi}\sum_{n=1}^{\infty}\text{sin}{nx}\text{cos}mx \text{dx}=0 \\ &\int_{-\pi}^{\pi}\sum_{n=1}^{\infty}\text{cos}{nx}\text{cos}mx\text{dx}=0(n\neq m) \\ &\int_{-\pi}^{\pi}\sum_{n=1}^{\infty}\text{cos}{nx}\text{cos}mx\text{dx}=\pi(n=m) \end{aligned}
∫−ππcosmxdx=∫−ππcosmxsin0xdx=0∫−ππn=1∑∞sinnxcosmxdx=0∫−ππn=1∑∞cosnxcosmxdx=0(n=m)∫−ππn=1∑∞cosnxcosmxdx=π(n=m)
因此 (2) 式可以写作如下形式
(
n
=
m
)
(n=m)
(n=m):
∫
−
π
π
f
(
x
)
cos
m
x
dx
=
∫
−
π
π
f
(
x
)
cos
n
x
dx
=
a
n
∫
−
π
π
∑
n
=
1
∞
cos
n
x
cos
m
x
dx
=
a
n
π
(3)
\begin{aligned} \int_{-\pi}^{\pi}f(x)\text{cos}mx\text{dx}&=\int_{-\pi}^{\pi}f(x)\text{cos}nx\text{dx} \\ &=a_n\int_{-\pi}^{\pi}\sum_{n=1}^{\infty}\text{cos}{nx}\text{cos}mx\text{dx}\\ &=a_n\pi \\ \end{aligned} \tag{3}
∫−ππf(x)cosmxdx=∫−ππf(x)cosnxdx=an∫−ππn=1∑∞cosnxcosmxdx=anπ(3)
可以求得:
a
n
=
1
π
∫
−
π
π
f
(
x
)
cos
n
x
dx
(4)
\begin{aligned} a_n=\frac{1}{\pi}\int_{-\pi}^{\pi}f(x)\text{cos}nx\text{dx} \end{aligned} \tag{4}
an=π1∫−ππf(x)cosnxdx(4)
同理,我们将(1)式两边同时乘以
s
i
n
m
x
sinmx
sinmx,然后两边同时在
[
−
π
,
π
]
[-\pi,\pi]
[−π,π] 做积分,可以求得:
b
n
=
1
π
∫
−
π
π
f
(
x
)
sin
n
x
dx
(5)
\begin{aligned} b_n=\frac{1}{\pi}\int_{-\pi}^{\pi}f(x)\text{sin}nx\text{dx} \end{aligned} \tag{5}
bn=π1∫−ππf(x)sinnxdx(5)
同理,对(1)式两边同时在
[
−
π
,
π
]
[-\pi,\pi]
[−π,π] 做积分,可以求得:
a
0
=
1
2
π
∫
−
π
π
f
(
x
)
dx
(6)
\begin{aligned} a_0=\frac{1}{2\pi}\int_{-\pi}^{\pi}f(x)\text{dx} \end{aligned} \tag{6}
a0=2π1∫−ππf(x)dx(6)
对比(4)式和(6)式,不难看出当(4)式的
n
=
0
n=0
n=0 时,它和
a
0
a_0
a0 只差了一个
1
2
\frac{1}{2}
21,因此为了将它们统一在一个表达式中,通常将(1)式写作如下的形式:
f
(
x
)
=
a
0
2
+
∑
n
=
1
∞
(
a
n
cos
n
x
+
b
n
sin
n
x
)
(7)
\begin{aligned} f(x)&=\frac{a_0}{2}+\sum_{n=1}^{\infty}(a_n\text{cos}{nx}+b_n\text{sin}{nx}) \\ \end{aligned} \tag{7}
f(x)=2a0+n=1∑∞(ancosnx+bnsinnx)(7)
最终得到的 (7) 式,就是三角级数的表示形式。当
a
n
a_n
an 和
b
n
b_n
bn 由上面的 (4)(5) 式确定时,(7) 式即为
f
(
x
)
f(x)
f(x) 的傅里叶级数,称此时的
a
n
a_n
an 和
b
n
b_n
bn 为傅里叶系数。(这里
f
(
x
)
f(x)
f(x) 周期为
2
π
2\pi
2π)
0x02 傅里叶级数的三角形式
上面我们推导傅里叶级数时,假设了周期为
2
π
2\pi
2π ,得到了如下形式的傅里叶级数:
f
(
x
)
=
a
0
2
+
∑
n
=
1
∞
(
a
n
cos
n
x
+
b
n
sin
n
x
)
{
a
n
=
1
π
∫
−
π
π
f
(
x
)
cos
n
x
dx
b
n
=
1
π
∫
−
π
π
f
(
x
)
sin
n
x
dx
(8)
\begin{aligned} f(x)&=\frac{a_0}{2}+\sum_{n=1}^{\infty}(a_n\text{cos}{nx}+b_n\text{sin}{nx}) \\ &\left\{\begin{array}{ll} a_n=\frac{1}{\pi}\int_{-\pi}^{\pi}f(x)\text{cos}nx\text{dx} \\ b_n=\frac{1}{\pi}\int_{-\pi}^{\pi}f(x)\text{sin}nx\text{dx} \end{array}\right. \end{aligned} \tag {8}
f(x)=2a0+n=1∑∞(ancosnx+bnsinnx){an=π1∫−ππf(x)cosnxdxbn=π1∫−ππf(x)sinnxdx(8)
现在可以进一步推广至周期为
T
T
T,即
f
(
t
)
=
f
(
t
+
T
)
f(t)=f(t+T)
f(t)=f(t+T),我们进行换元
x
=
2
π
T
t
x=\frac{2\pi}{T}t
x=T2πt,
t
=
T
2
π
x
t=\frac{T}{2\pi}x
t=2πTx,然后令
f
(
t
)
=
f
(
2
L
π
x
)
:
=
g
(
x
)
f(t)=f(\frac{2L}{\pi}x):= g(x)
f(t)=f(π2Lx):=g(x),则有
g
(
x
)
=
g
(
x
+
2
π
)
g(x)=g(x+2\pi)
g(x)=g(x+2π)。
所以可以将
g
(
x
)
g(x)
g(x) 用 (8) 式傅里叶级数展开:
g
(
x
)
=
a
0
2
+
∑
n
=
1
∞
(
a
n
cos
n
x
+
b
n
sin
n
x
)
{
a
n
=
1
π
∫
−
π
π
g
(
x
)
cos
n
x
dx
b
n
=
1
π
∫
−
π
π
g
(
x
)
sin
n
x
dx
(9)
\begin{aligned} g(x)&=\frac{a_0}{2}+\sum_{n=1}^{\infty}(a_n\text{cos}{nx}+b_n\text{sin}{nx}) \\ &\left\{\begin{array}{ll} a_n=\frac{1}{\pi}\int_{-\pi}^{\pi}g(x)\text{cos}nx\text{dx} \\ b_n=\frac{1}{\pi}\int_{-\pi}^{\pi}g(x)\text{sin}nx\text{dx} \end{array}\right. \end{aligned} \tag {9}
g(x)=2a0+n=1∑∞(ancosnx+bnsinnx){an=π1∫−ππg(x)cosnxdxbn=π1∫−ππg(x)sinnxdx(9)
然后对 (9) 式再进行一次换元
x
=
2
π
T
t
x=\frac{2\pi}{T}t
x=T2πt,得到:
f
(
t
)
=
a
0
2
+
∑
n
=
1
∞
(
a
n
cos
n
2
π
L
t
+
b
n
sin
n
2
π
L
t
)
=
a
0
2
+
∑
n
=
1
∞
(
a
n
cos
n
ω
t
+
b
n
sin
n
ω
t
)
(10)
\begin{aligned} f(t)&=\frac{a_0}{2}+\sum_{n=1}^{\infty}(a_n\text{cos}{n\frac{2\pi}{L}t}+b_n\text{sin}{n\frac{2\pi}{L}t}) \\ &=\frac{a_0}{2}+\sum_{n=1}^{\infty}(a_n\text{cos}{n \omega t}+b_n\text{sin}{n \omega t}) \end{aligned} \tag{10}
f(t)=2a0+n=1∑∞(ancosnL2πt+bnsinnL2πt)=2a0+n=1∑∞(ancosnωt+bnsinnωt)(10)
同理对
a
n
a_n
an 和
b
n
b_n
bn 进行换元得到:
a
n
=
1
π
∫
−
π
π
g
(
x
)
cos
n
x
dx
=
1
π
∫
−
T
2
T
2
f
(
t
)
cos
n
2
π
T
t
d
2
π
T
t
=
1
π
2
π
T
∫
−
T
2
T
2
f
(
t
)
cos
n
2
π
T
t
dt
=
2
T
∫
0
T
f
(
t
)
cos
n
ω
t
dt
(11)
\begin{aligned} a_n&=\frac{1}{\pi}\int_{-\pi}^{\pi}g(x)\text{cos}nx\text{dx} \\ &= \frac{1}{\pi}\int_{-\frac{T}{2}}^{\frac{T}{2}}f(t)\text{cos}{n\frac{2\pi}{T}t}\text{d}\frac{2\pi}{T}t \\ &=\frac{1}{\pi}\frac{2\pi}{T}\int_{-\frac{T}{2}}^{\frac{T}{2}}f(t)\text{cos}{n\frac{2\pi}{T}t}\text{dt} \\ &=\frac{2}{T}\int_{0}^{T}f(t)\text{cos}{n \omega t}\text{dt} \end{aligned} \tag {11}
an=π1∫−ππg(x)cosnxdx=π1∫−2T2Tf(t)cosnT2πtdT2πt=π1T2π∫−2T2Tf(t)cosnT2πtdt=T2∫0Tf(t)cosnωtdt(11)
b n = 2 T ∫ 0 T f ( t ) sin n ω t dt (12) \begin{aligned} b_n&=\frac{2}{T}\int_{0}^{T}f(t)\text{sin}{n \omega t}\text{dt} \end{aligned} \tag {12} bn=T2∫0Tf(t)sinnωtdt(12)
综上,也就说明了任何周期函数都能用傅里叶级数来表示为三角函数的形式。对于周期函数
f
(
t
)
f(t)
f(t),周期为
T
T
T ,角频率为
ω
0
=
2
π
f
0
=
2
π
T
\omega_0=2\pi f_0=\frac{2\pi}{T}
ω0=2πf0=T2π,其三角形式的傅里叶级数如下:
f
(
t
)
=
a
0
2
+
∑
n
=
1
∞
(
a
n
cos
n
ω
0
t
+
b
n
sin
n
ω
0
t
)
{
a
n
=
2
T
∫
0
T
f
(
t
)
cos
n
ω
0
t
dt
b
n
=
2
T
∫
0
T
f
(
t
)
sin
n
ω
0
t
dt
(13)
\begin{aligned} f(t)&=\frac{a_0}{2}+\sum_{n=1}^{\infty}(a_n\text{cos}{n \omega_0 t}+b_n\text{sin}{n \omega_0 t}) \\ &\left\{\begin{array}{ll} a_n=\frac{2}{T}\int_{0}^{T}f(t)\text{cos}{n \omega_0 t}\text{dt} \\ b_n=\frac{2}{T}\int_{0}^{T}f(t)\text{sin}{n \omega_0 t}\text{dt} \end{array}\right. \end{aligned} \tag {13}
f(t)=2a0+n=1∑∞(ancosnω0t+bnsinnω0t){an=T2∫0Tf(t)cosnω0tdtbn=T2∫0Tf(t)sinnω0tdt(13)
0x03 傅里叶级数的复指数形式
上一节我们推导了三角形式的傅里叶级数,其物理含义明确;但是实际操作上,指数形式的傅里叶级会更加便于数学处理,而且很容易与后面介绍的傅里叶变换统一起来。因此这一节我们就用欧拉公式将 (13) 式的傅里叶级数表示为复指数的形式。(欧拉公式的证明可见附录2)
首先,由欧拉公式
e
i
θ
=
cos
θ
+
i
sin
θ
e^{i\theta}=\text{cos}{\theta}+i\text{sin}{\theta}
eiθ=cosθ+isinθ ,分别将
θ
\theta
θ 和
−
θ
-\theta
−θ 代入其中再消元,得到三角函数的复指数形式如下:
cos
θ
=
1
2
(
e
i
θ
+
e
−
i
θ
)
sin
θ
=
−
1
2
i
(
e
i
θ
−
e
−
i
θ
)
(14)
\begin{aligned} \text{cos}{\theta} = \frac{1}{2}\left(e^{i\theta}+e^{-i\theta}\right) \\ \text{sin}{\theta} = -\frac{1}{2}i\left(e^{i\theta}-e^{-i\theta}\right) \end{aligned} \tag {14}
cosθ=21(eiθ+e−iθ)sinθ=−21i(eiθ−e−iθ)(14)
将 (14) 式的代入 (13) 式的傅里叶级数中得到:
f ( t ) = a 0 2 + ∑ n = 1 ∞ [ a n 1 2 ( e i n ω t + e − i n ω t ) − b n 1 2 i ( e i n ω t − e − i n ω t ) ] = a 0 2 + ∑ n = 1 ∞ ( a n − i b n 2 e i n ω t + a n + i b n 2 e − i n ω t ) = a 0 2 + ∑ n = 1 ∞ a n − i b n 2 e i n ω t + ∑ n = 1 ∞ a n + i b n 2 e − i n ω t = ∑ n = 0 0 a 0 2 e 0 ⋅ i ω t + ∑ n = 1 ∞ a n − i b n 2 e i n ω t + ∑ n = − 1 − ∞ a − n + i b − n 2 e i n ω t = ∑ n = − ∞ ∞ C n e i n ω t w h e r e C n = { a 0 2 , n = 0 a n − i b n 2 , n = 1 , 2 , 3 , … a − n + i b − n 2 , n = − 1 , − 2 , − 3 , … (15) \begin{aligned} f(t)&=\frac{a_0}{2}+\sum_{n=1}^{\infty}\left[a_n\frac{1}{2}\left(e^{in \omega t}+e^{-in \omega t}\right)-b_n\frac{1}{2}i\left(e^{in \omega t}-e^{-in \omega t}\right)\right] \\ &=\frac{a_0}{2}+\sum_{n=1}^{\infty}\left(\frac{a_n-ib_n}{2}e^{in \omega t}+\frac{a_n+ib_n}{2}e^{-in \omega t}\right) \\ &=\frac{a_0}{2} +\sum_{n=1}^{\infty}\frac{a_n-ib_n}{2}e^{in \omega t} +\sum_{n=1}^{\infty}\frac{a_n+ib_n}{2}e^{-in \omega t} \\ &=\sum_{n=0}^{0}\frac{a_0}{2}e^{0·i \omega t}+\sum_{n=1}^{\infty}\frac{a_n-ib_n}{2}e^{in \omega t} +\sum_{n=-1}^{-\infty}\frac{a_{-n} +ib_{-n}}{2}e^{in \omega t} \\ &=\sum_{n=-\infty}^{\infty} C_n e^{in \omega t} \\ \\ where \space C_n=&\left\{\begin{array}{ll} \frac{a_0}{2},n=0 \\ \frac{a_n-ib_n}{2},n=1,2,3,\dots \\ \frac{a_{-n}+ib_{-n}}{2},n=-1,-2,-3,\dots \end{array}\right. \end{aligned} \tag {15} f(t)where Cn==2a0+n=1∑∞[an21(einωt+e−inωt)−bn21i(einωt−e−inωt)]=2a0+n=1∑∞(2an−ibneinωt+2an+ibne−inωt)=2a0+n=1∑∞2an−ibneinωt+n=1∑∞2an+ibne−inωt=n=0∑02a0e0⋅iωt+n=1∑∞2an−ibneinωt+n=−1∑−∞2a−n+ib−neinωt=n=−∞∑∞Cneinωt⎩⎨⎧2a0,n=02an−ibn,n=1,2,3,…2a−n+ib−n,n=−1,−2,−3,…(15)
上面第 4 个等号第三项将 n n n 换元为 − n -n −n,并将第一项写成单个累和的形式( e 0 ⋅ i ω t = 1 e^{0·i \omega t}=1 e0⋅iωt=1);最后一个等号,可以看到三项的累和已经从 − ∞ -\infty −∞ 加到了 + ∞ +\infty +∞,我们将 e i n ω t e^{in \omega t} einωt 提出来, 剩余项项用 C n C_n Cn 表示。
进一步计算 C n C_n Cn 的形式,将 (13) 式的傅里叶系数代入 C n C_n Cn 中得到:
(1)当
n
=
0
n=0
n=0 时:
C
n
=
a
0
2
=
1
2
2
T
∫
0
T
f
(
t
)
dt
=
1
T
∫
0
T
f
(
t
)
dt
=
1
T
∫
0
T
f
(
t
)
e
−
i
n
ω
t
dt
(16)
\begin{aligned} C_n=\frac{a_0}{2}&=\frac{1}{2}\frac{2}{T}\int_{0}^{T}f(t)\text{dt} \\ &=\frac{1}{T}\int_{0}^{T}f(t)\text{dt} \\ &=\frac{1}{T}\int_{0}^{T}f(t)e^{-in \omega t}\text{dt} \end{aligned} \tag{16}
Cn=2a0=21T2∫0Tf(t)dt=T1∫0Tf(t)dt=T1∫0Tf(t)e−inωtdt(16)
上面的最后一个等号,因为在此情况下
n
=
0
n=0
n=0,所以加上
e
−
i
n
ω
t
e^{-in \omega t}
e−inωt 不改变结果。
(2)当
n
=
1
,
2
,
3
,
…
n=1,2,3,\dots
n=1,2,3,… 时:
C
n
=
a
n
−
i
b
n
2
=
1
2
(
2
T
∫
0
T
f
(
t
)
cos
n
ω
t
dt
−
i
2
T
∫
0
T
f
(
t
)
sin
n
ω
t
dt
)
=
1
T
∫
0
T
f
(
t
)
(
cos
n
ω
t
−
i
sin
n
ω
t
)
dt
=
1
T
∫
0
T
f
(
t
)
e
−
i
n
ω
t
dt
(17)
\begin{aligned} C_n&=\frac{a_n-ib_n}{2}\\ &=\frac{1}{2}\left(\frac{2}{T}\int_{0}^{T}f(t)\text{cos}{n \omega t}\text{dt}-i\frac{2}{T}\int_{0}^{T}f(t)\text{sin}{n \omega t}\text{dt}\right) \\ &=\frac{1}{T}\int_{0}^{T}f(t) \left(\text{cos}{n \omega t}-i\text{sin}{n \omega t}\right)\text{dt} \\ &=\frac{1}{T}\int_{0}^{T}f(t)e^{-in \omega t}\text{dt} \\ \end{aligned} \tag{17}
Cn=2an−ibn=21(T2∫0Tf(t)cosnωtdt−iT2∫0Tf(t)sinnωtdt)=T1∫0Tf(t)(cosnωt−isinnωt)dt=T1∫0Tf(t)e−inωtdt(17)
上面第3个等式可由欧拉公式得到:
cos
n
ω
t
−
i
sin
n
ω
t
=
cos
(
−
n
ω
t
)
+
i
sin
(
−
n
ω
t
)
=
e
−
i
n
ω
t
\text{cos}{n \omega t}-i\text{sin}{n \omega t}=\text{cos}{(-n \omega t)}+i\text{sin}{(-n \omega t)}=e^{-in \omega t}
cosnωt−isinnωt=cos(−nωt)+isin(−nωt)=e−inωt
(3)当
n
=
−
1
,
−
2
,
−
3
,
…
n=-1,-2,-3,\dots
n=−1,−2,−3,… 时:
C
n
=
a
−
n
+
i
b
−
n
2
=
1
2
(
2
T
∫
0
T
f
(
t
)
cos
(
−
n
ω
t
)
dt
+
i
2
T
∫
0
T
f
(
t
)
sin
(
−
n
ω
t
)
dt
)
=
1
2
(
2
T
∫
0
T
f
(
t
)
cos
n
ω
t
dt
−
i
2
T
∫
0
T
f
(
t
)
sin
n
ω
t
dt
)
=
1
T
∫
0
T
f
(
t
)
(
cos
n
ω
t
−
i
sin
n
ω
t
)
dt
=
1
T
∫
0
T
f
(
t
)
e
−
i
n
ω
t
dt
(18)
\begin{aligned} C_n&=\frac{a_{-n}+ib_{-n}}{2}\\ &=\frac{1}{2}\left(\frac{2}{T}\int_{0}^{T}f(t)\text{cos}{(-n \omega t)}\text{dt}+i\frac{2}{T}\int_{0}^{T}f(t)\text{sin}{(-n \omega t)}\text{dt}\right) \\ &=\frac{1}{2}\left(\frac{2}{T}\int_{0}^{T}f(t)\text{cos}{n \omega t}\text{dt}-i\frac{2}{T}\int_{0}^{T}f(t)\text{sin}{n \omega t}\text{dt}\right) \\ &=\frac{1}{T}\int_{0}^{T}f(t) \left(\text{cos}{n \omega t}-i\text{sin}{n \omega t}\right)\text{dt} \\ &=\frac{1}{T}\int_{0}^{T}f(t)e^{-in \omega t}\text{dt} \\ \end{aligned} \tag{18}
Cn=2a−n+ib−n=21(T2∫0Tf(t)cos(−nωt)dt+iT2∫0Tf(t)sin(−nωt)dt)=21(T2∫0Tf(t)cosnωtdt−iT2∫0Tf(t)sinnωtdt)=T1∫0Tf(t)(cosnωt−isinnωt)dt=T1∫0Tf(t)e−inωtdt(18)
观察 (16)(17)(18) 式,可以看到三个情况的最终结果竟然统一了!也就是对于 $ \forall n=0,±1,±2,±3,\dots$,有:
C
n
=
1
T
∫
0
T
f
(
t
)
e
−
i
n
ω
t
dt
(19)
\begin{aligned} C_n=\frac{1}{T}\int_{0}^{T}f(t)e^{-in \omega t}\text{dt} \end{aligned} \tag{19}
Cn=T1∫0Tf(t)e−inωtdt(19)
综上,对于周期函数
f
(
t
)
f(t)
f(t),周期为
T
T
T,角频率为
w
0
=
2
π
f
0
=
2
π
T
w_0=2\pi f_0=\frac{2\pi}{T}
w0=2πf0=T2π,其复指数形式的傅里叶级数为:
f
(
t
)
=
∑
n
=
−
∞
∞
C
n
e
i
n
ω
0
t
,
C
n
=
1
T
∫
0
T
f
(
t
)
e
−
i
n
ω
0
t
dt
,
n
=
0
,
±
1
,
±
2
,
±
3
,
…
(20)
\begin{aligned} &f(t)=\sum_{n=-\infty}^{\infty} C_n e^{in \omega_0 t}, \\ &C_n=\frac{1}{T}\int_{0}^{T}f(t)e^{-in \omega_0 t}\text{dt},n=0,±1,±2,±3,\dots \end{aligned} \tag {20}
f(t)=n=−∞∑∞Cneinω0t,Cn=T1∫0Tf(t)e−inω0tdt,n=0,±1,±2,±3,…(20)
0x04 傅里叶变换
傅里叶变换实际上就是对傅里叶级数的进一步扩展,即不再要求 f ( t ) f(t) f(t) 为周期函数,换句话也就是设周期 T T T 趋于 ∞ \infty ∞。
现在假设有一个周期为
T
T
T 的函数
f
T
(
t
)
=
f
T
(
t
+
T
)
f_T(t)=f_T(t+T)
fT(t)=fT(t+T),我们用复指数形式的傅里叶级数将其展开得到:
f
T
(
t
)
=
∑
n
=
−
∞
∞
C
n
e
i
n
ω
0
t
,
C
n
=
1
T
∫
0
T
f
T
(
t
)
e
−
i
n
ω
0
t
dt
,
n
=
0
,
±
1
,
±
2
,
±
3
,
…
(21)
\begin{aligned} &f_T(t)=\sum_{n=-\infty}^{\infty} C_n e^{in \omega_0 t}, \\ &C_n=\frac{1}{T}\int_{0}^{T}f_T(t)e^{-in \omega_0 t}\text{dt},n=0,±1,±2,±3,\dots \end{aligned} \tag {21}
fT(t)=n=−∞∑∞Cneinω0t,Cn=T1∫0TfT(t)e−inω0tdt,n=0,±1,±2,±3,…(21)
其中
ω
0
=
2
π
T
\omega_0=\frac{2 \pi}{T}
ω0=T2π 称为基频率。可以看到,实际上是由
C
n
C_n
Cn 区分不同函数的,因为对于任意的傅里叶变换都是
C
n
C_n
Cn 乘对应的复指数,再累和。假设
f
T
(
t
)
f_T(t)
fT(t) 的时域图像如下:
我们画出它在复平面上的频域图像如下:
如果觉得上面这张三维复平面下的频域图不好理解,可以参考下面这张图(来自@李泽光一个相关回答里的图,图中
X
n
X_n
Xn 与这里的
C
n
C_n
Cn 含义类似):
所谓傅里叶级数就是将时域下的周期连续函数(上图),转换为频域上的非周期离散函数(下图),也就是说对于频域图像做运算: … , C − 2 e − 2 i ω 0 t + C − 1 e − i ω 0 t + C 0 e 0 + C 1 e i ω 0 t + C 2 e 2 i ω 0 t , … \dots,C_{-2}e^{-2i \omega_0t}+C_{-1}e^{-i \omega_0t}+C_{0}e^{0}+C_{1}e^{ i\omega_0t}+C_{2}e^{2i \omega_0t},\dots …,C−2e−2iω0t+C−1e−iω0t+C0e0+C1eiω0t+C2e2iω0t,… 也就能得到时域下的函数图像。
而傅里叶变换要做的就是将时域下的非周期连续函数转换到频域上的非周期连续函数。

前面说了非周期函数可以将其看作周期无限大( T → ∞ T→\infty T→∞),也就是 lim T → ∞ f T ( t ) = f ( t ) \lim_{T→\infty}f_T(t)=f(t) limT→∞fT(t)=f(t)
我们令频域图像上的间隔为 Δ ω = ( n + 1 ) ω 0 − n ω 0 = ω 0 = 2 π T \Delta \omega=(n+1)\omega_0-n\omega_0=\omega_0=\frac{2\pi}{T} Δω=(n+1)ω0−nω0=ω0=T2π,当 T → + ∞ T→+\infty T→+∞ 时, Δ → − ∞ \Delta→-\infty Δ→−∞。从图像上看也就意味着每个 C i C_i Ci 之间的距离越来越来小,极限情况频域图像每个 C i C_i Ci 就挨在一起了,也就变成了连续的形式;并且离散情况下的横轴 n ω 0 n\omega_0 nω0,在连续的情况下就可以用一个连续的变量 ω \omega ω 表示。
进一步推导,将 (19) 式中的
C
n
C_n
Cn 代入
f
T
(
t
)
f_T(t)
fT(t) ,且
T
=
2
π
Δ
ω
T=\frac{2\pi}{\Delta \omega}
T=Δω2π:
f
(
t
)
=
∑
n
=
−
∞
∞
1
T
∫
0
T
f
(
t
)
e
−
i
n
ω
0
t
dt
e
i
n
ω
0
t
=
∑
n
=
−
∞
∞
Δ
ω
2
π
∫
0
T
f
T
(
t
)
e
−
i
n
ω
0
t
dt
e
i
n
ω
0
t
(22)
\begin{aligned} f(t)&=\sum_{n=-\infty}^{\infty} \frac{1}{T}\int_{0}^{T}f(t)e^{-in \omega_0 t}\text{dt} e^{in \omega_0 t} \\ &=\sum_{n=-\infty}^{\infty} \frac{\Delta \omega}{2\pi}\int_{0}^{T}f_T(t)e^{-in \omega_0 t}\text{dt} e^{in \omega_0 t} \end{aligned} \tag {22}
f(t)=n=−∞∑∞T1∫0Tf(t)e−inω0tdteinω0t=n=−∞∑∞2πΔω∫0TfT(t)e−inω0tdteinω0t(22)
对于 (22) 式,当我们令
T
→
∞
T→\infty
T→∞ 时,相应的有
∫
0
T
d
t
→
∫
−
∞
+
∞
\int_{0}^{T}dt→\int_{-\infty}^{+\infty}
∫0Tdt→∫−∞+∞,
n
ω
0
→
ω
n\omega_0→\omega
nω0→ω,
∑
n
=
−
∞
∞
Δ
ω
→
∫
−
∞
+
∞
d
ω
\sum_{n=-\infty}^{\infty}\Delta\omega→\int_{-\infty}^{+\infty}d\omega
∑n=−∞∞Δω→∫−∞+∞dω,得到:
f
(
t
)
=
1
2
π
∫
−
∞
∞
∫
−
∞
∞
f
(
t
)
e
−
i
ω
t
dt
e
i
ω
t
d
ω
(23)
\begin{aligned} f(t)&=\frac{1}{2\pi}\int_{-\infty}^{\infty}\int_{-\infty}^{\infty}f(t)e^{-i \omega t}\text{dt} e^{i \omega t}\text{d}\omega \end{aligned} \tag {23}
f(t)=2π1∫−∞∞∫−∞∞f(t)e−iωtdteiωtdω(23)
最后,我们令上面 (23) 式的中间部分为为
F
(
ω
)
F(\omega)
F(ω) 如下:
F
(
w
)
=
∫
−
∞
∞
f
(
t
)
e
−
i
ω
t
dt
(24)
\begin{aligned} F(w)&=\int_{-\infty}^{\infty}f(t)e^{-i \omega t}\text{dt} \end{aligned} \tag {24}
F(w)=∫−∞∞f(t)e−iωtdt(24)
得到的 (24) 式即为傅里叶变换,而原来的 (23) 式即为傅里叶逆变换,可写成如下形式:
f
(
t
)
=
1
2
π
∫
−
∞
∞
F
(
ω
)
e
i
ω
t
d
ω
(25)
\begin{aligned} f(t)&=\frac{1}{2\pi}\int_{-\infty}^{\infty}F(\omega) e^{i \omega t}\text{d}\omega \end{aligned} \tag {25}
f(t)=2π1∫−∞∞F(ω)eiωtdω(25)
附录1:三角函数正交性证明
当
n
≠
m
n≠m
n=m时,由积化和差及牛顿莱布尼茨公式得:
∫
−
π
π
cos
n
x
⋅
cos
m
x
dx
=
∫
−
π
π
1
2
[
cos
(
n
x
−
m
x
)
+
cos
(
n
x
+
m
x
)
]
d
x
=
1
2
[
∫
−
π
π
cos
(
n
−
m
)
x
dx
+
∫
−
π
π
cos
(
n
+
m
)
x
dx
]
=
1
2
[
1
n
−
m
sin
(
n
−
m
)
x
∣
−
π
π
+
1
n
+
m
sin
(
n
+
m
)
x
∣
−
π
π
]
=
0
\begin{aligned} \int_{-\pi}^{\pi} \text{cos}nx·\text{cos}mx\text{dx} &= \int_{-\pi}^{\pi} \frac{1}{2}\left[\text{cos}(nx-mx) + \text{cos}(nx+mx)\right]dx \\ &=\frac{1}{2} \left[\int_{-\pi}^{\pi}\text{cos}(n-m)x\text{dx} + \int_{-\pi}^{\pi}\text{cos}(n+m)x\text{dx}\right] \\ &=\frac{1}{2} \left[\frac{1}{n-m}\text{sin}(n-m)x\bigg|_{-\pi}^{\pi}+ \frac{1}{n+m}\text{sin}(n+m)x\bigg|_{-\pi}^{\pi} \right] \\ &=0 \end{aligned}
∫−ππcosnx⋅cosmxdx=∫−ππ21[cos(nx−mx)+cos(nx+mx)]dx=21[∫−ππcos(n−m)xdx+∫−ππcos(n+m)xdx]=21[n−m1sin(n−m)x∣∣∣∣−ππ+n+m1sin(n+m)x∣∣∣∣−ππ]=0
∫ − π π cos n x ⋅ sin m x dx = ∫ − π π 1 2 [ sin ( n x + m x ) − sin ( n x − m x ) ] d x = 1 2 [ ∫ − π π sin ( n + m ) x dx + ∫ − π π sin ( n − m ) x dx ] = 0 \begin{aligned} \int_{-\pi}^{\pi} \text{cos}nx·\text{sin}mx\text{dx} &= \int_{-\pi}^{\pi} \frac{1}{2}\left[\text{sin}(nx+mx) - \text{sin}(nx-mx)\right]dx \\ &=\frac{1}{2} \left[\int_{-\pi}^{\pi}\text{sin}(n+m)x\text{dx} + \int_{-\pi}^{\pi}\text{sin}(n-m)x\text{dx}\right] \\ &=0 \end{aligned} ∫−ππcosnx⋅sinmxdx=∫−ππ21[sin(nx+mx)−sin(nx−mx)]dx=21[∫−ππsin(n+m)xdx+∫−ππsin(n−m)xdx]=0
同理:
∫
−
π
π
sin
n
x
⋅
cos
m
x
dx
=
0
∫
−
π
π
sin
n
x
⋅
sin
m
x
dx
=
0
\begin{aligned} \int_{-\pi}^{\pi} \text{sin}nx·\text{cos}mx\text{dx} &= 0 \\ \int_{-\pi}^{\pi} \text{sin}nx·\text{sin}mx\text{dx} &= 0 \end{aligned}
∫−ππsinnx⋅cosmxdx∫−ππsinnx⋅sinmxdx=0=0
当
n
=
m
n=m
n=m 时,
∫
−
π
π
cos
n
x
⋅
sin
m
x
dx
=
∫
−
π
π
sin
n
x
⋅
cos
m
x
dx
=
∫
−
π
π
sin
n
x
⋅
cos
n
x
dx
=
∫
−
π
π
1
2
[
sin
2
n
x
+
s
i
n
0
x
]
dx
=
1
2
∫
−
π
π
sin
2
n
x
dx
=
0
\begin{aligned} \int_{-\pi}^{\pi} \text{cos}nx·\text{sin}mx\text{dx} &= \int_{-\pi}^{\pi} \text{sin}nx·\text{cos}mx\text{dx} \\ &=\int_{-\pi}^{\pi} \text{sin}nx·\text{cos}nx\text{dx} \\ &=\int_{-\pi}^{\pi} \frac{1}{2}\left[\text{sin}2nx+sin0x\right]\text{dx} \\ &=\frac{1}{2}\int_{-\pi}^{\pi}\text{sin}2nx\text{dx} \\ &=0 \end{aligned}
∫−ππcosnx⋅sinmxdx=∫−ππsinnx⋅cosmxdx=∫−ππsinnx⋅cosnxdx=∫−ππ21[sin2nx+sin0x]dx=21∫−ππsin2nxdx=0
附录2:欧拉公式的证明
欧拉公式: e i θ = cos θ + i sin θ e^{i\theta}=\text{cos}{\theta}+i\text{sin}{\theta} eiθ=cosθ+isinθ
证明:
首先,令
f
(
θ
)
=
e
i
θ
cos
θ
+
i
sin
θ
\begin{aligned} f(\theta)=\frac{e^{i\theta}}{\text{cos}{\theta}+i\text{sin}{\theta}} \end{aligned}
f(θ)=cosθ+isinθeiθ
对上式进行求导得到:
f
′
(
θ
)
=
i
e
i
θ
(
cos
θ
+
i
sin
θ
)
−
e
i
θ
(
−
sin
θ
+
i
cos
θ
)
(
cos
θ
+
i
sin
θ
)
2
=
i
e
i
θ
cos
θ
−
e
i
θ
sin
θ
+
e
i
θ
sin
θ
−
i
e
i
θ
cos
θ
(
cos
θ
+
i
sin
θ
)
2
=
0
\begin{aligned} f'(\theta)&=\frac{ie^{i\theta}(\text{cos}{\theta}+i\text{sin}{\theta})-e^{i\theta}(-\text{sin}{\theta}+i\text{cos}{\theta})}{(\text{cos}{\theta}+i\text{sin}{\theta})^2} \\ &=\frac{ie^{i\theta}\text{cos}{\theta}-e^{i\theta}\text{sin}{\theta}+e^{i\theta}\text{sin}{\theta}-ie^{i\theta}\text{cos}{\theta}}{(\text{cos}{\theta}+i\text{sin}{\theta})^2} \\ &=0 \end{aligned}
f′(θ)=(cosθ+isinθ)2ieiθ(cosθ+isinθ)−eiθ(−sinθ+icosθ)=(cosθ+isinθ)2ieiθcosθ−eiθsinθ+eiθsinθ−ieiθcosθ=0
因为
f
′
(
θ
)
=
0
f'(\theta)=0
f′(θ)=0,说明
f
(
θ
)
f(\theta)
f(θ) 为一个常数,不妨设
θ
=
0
\theta=0
θ=0,则:
f
(
θ
)
=
f
(
0
)
=
e
i
0
cos
0
+
i
sin
0
=
1
\begin{aligned} f(\theta)=f(0)=\frac{e^{i0}}{\text{cos}{0}+i\text{sin}{0}}=1 \end{aligned}
f(θ)=f(0)=cos0+isin0ei0=1
所以得证:
f
(
θ
)
=
e
i
θ
cos
θ
+
i
sin
θ
=
1
e
i
θ
=
cos
θ
+
i
sin
θ
\begin{aligned} &f(\theta)=\frac{e^{i\theta}}{\text{cos}{\theta}+i\text{sin}{\theta}}=1 \\ &e^{i\theta}=\text{cos}{\theta}+i\text{sin}{\theta} \end{aligned}
f(θ)=cosθ+isinθeiθ=1eiθ=cosθ+isinθ