1 为什么要进行傅里叶分析
信号分析方法主流方法有:
(1)时域分析:以冲激信号为基本信号,任意输入信号可分解为一系列冲激信号;
(2)频域分析:以正弦信号和虚指数信号为基本信号,将任意输入信号分解为一系列不同频率的正弦信号或者虚指数信号之和。独立变量是频率,分析是在频率空间进行的,故称为频率域分析,或者频域分析。
傅里叶分析就是将时域信号转换成频域分析,其中傅里叶分析理论包括傅里叶级数、傅里叶积分和傅里叶变换。
2 傅里叶分析之傅里叶级数
从为什么信号都能够用正弦信号表示?一文我们可知信号可以使用完备正交函数集来分解。
2.1 三角函数形式的傅里叶级数
其中三角函数集{1, c o s ( n Ω t ) , s i n ( n Ω t ) , n = 1 , 2 , ⋯ 1 cos(n\Omega t),sin(n\Omega t),n=1,2,\cdots1 cos(nΩt),sin(nΩt),n=1,2,⋯1}
设周期函数
f
(
t
)
f(t)
f(t),其周期为T,角频率为
Ω
=
2
π
/
T
\Omega=2\pi/T
Ω=2π/T,满足狄里赫利条件时,可展开为三角形式傅里叶级数。
即
f
(
t
)
=
a
0
2
+
∑
n
=
1
∞
a
n
c
o
s
(
n
Ω
t
)
+
∑
n
=
1
∞
b
n
s
i
n
(
n
Ω
t
)
(
1
)
f(t)=\frac{a_0}{2}+\sum_{n=1}^{\infty}a_ncos(n\Omega t)+\sum_{n=1}^{\infty}b_nsin(n\Omega t)(1)
f(t)=2a0+n=1∑∞ancos(nΩt)+n=1∑∞bnsin(nΩt)(1)
直流分量
:
a
0
2
=
1
T
∫
−
T
2
T
2
f
(
t
)
d
t
直流分量: \frac{a_0}{2}=\frac{1}{T}\int_{-\frac{T}{2}}^{\frac{T}{2}}f(t)dt
直流分量:2a0=T1∫−2T2Tf(t)dt
余弦分量系数
:
a
n
=
2
T
∫
−
T
2
T
2
f
(
t
)
c
o
s
(
n
Ω
t
)
d
t
余弦分量系数:a_n=\frac{2}{T}\int_{-\frac{T}{2}}^{\frac{T}{2}}f(t)cos(n\Omega t)dt
余弦分量系数:an=T2∫−2T2Tf(t)cos(nΩt)dt
余弦分量系数
:
b
n
=
2
T
∫
−
T
2
T
2
f
(
t
)
s
i
n
(
n
Ω
t
)
d
t
余弦分量系数:b_n=\frac{2}{T}\int_{-\frac{T}{2}}^{\frac{T}{2}}f(t)sin(n\Omega t)dt
余弦分量系数:bn=T2∫−2T2Tf(t)sin(nΩt)dt
▶
\blacktriangleright
▶直流写成
a
0
2
?
?
?
\frac{a_0}{2}???
2a0???
实际上就相当于余弦分量系数a_n取n=0时候
直流分量
:
a
0
=
2
T
∫
−
T
2
T
2
f
(
t
)
d
t
直流分量: a_0=\frac{2} {T}\int_{-\frac{T}{2}}^{\frac{T}{2}}f(t)dt
直流分量:a0=T2∫−2T2Tf(t)dt
▶
\blacktriangleright
▶为什么傅里叶级数的系数是这样子的??
具体可参考保姆级推导傅里叶级数
这个时候大家就觉得麻烦,说有两个系数自己记不住,于是就想个方法统一一下,就有余弦形式的傅里叶级数产生,周期信号可以分解为直流和许多余弦分量。
主要使用和差化积将(1)式进行
s
i
n
和
c
o
s
sin和cos
sin和cos合并
f
(
t
)
=
A
0
2
+
∑
n
=
1
∞
A
n
c
o
s
(
n
Ω
t
+
φ
n
)
(
2
)
f(t)=\frac{A_0}{2}+\sum_{n=1}^{\infty}A_ncos(n\Omega t+ φ_n) (2)
f(t)=2A0+n=1∑∞Ancos(nΩt+φn)(2)
A
n
=
a
n
2
+
b
n
2
,
a
n
和
b
n
来自(
1
)式
A_n=\sqrt{a_n^2+b_n ^2},a_n 和b_n来自(1)式
An=an2+bn2,an和bn来自(1)式
φ
n
=
−
a
r
c
t
a
n
b
n
a
n
\varphi_n=-arctan\frac{b_n}{a_n}
φn=−arctananbn
其中
A
0
/
2
为直流分量,
A
1
c
o
s
(
Ω
t
+
φ
1
)
为基波或一次谐波
,
⋯
,
A
n
c
o
s
(
n
Ω
t
+
φ
n
)
为
n
次谐波
A_0/2为直流分量,A_1cos(\Omega t+ \varphi_1)为基波或一次谐波,\cdots,A_ncos(n\Omega t+ \varphi_n)为n次谐波
A0/2为直流分量,A1cos(Ωt+φ1)为基波或一次谐波,⋯,Ancos(nΩt+φn)为n次谐波
2.2指数形式的傅里叶级数
f
(
t
)
=
∑
n
=
−
∞
∞
F
n
e
j
n
Ω
t
f(t)=\sum_{n=-\infty}^{\infty}F_ne^{jn\Omega t}
f(t)=n=−∞∑∞FnejnΩt
F
n
=
1
T
∫
−
T
2
T
2
f
(
t
)
e
−
j
n
Ω
t
d
t
F_n=\frac{1}{T}\int_{-\frac{T}{2}}^{\frac{T}{2}}f(t)e^{-jn\Omega t}dt
Fn=T1∫−2T2Tf(t)e−jnΩtdt
F
n
F_n
Fn是频率为
n
Ω
n\Omega
nΩ的分量的系数,其中
F
0
=
A
0
/
2
F_0=A_0/2
F0=A0/2为直流分量
推导过程:
已知
f
(
t
)
=
A
0
2
+
∑
n
=
1
∞
A
n
c
o
s
(
n
Ω
t
+
φ
n
)
(
2
)
已知f(t)=\frac{A_0}{2}+\sum_{n=1}^{\infty}A_ncos(n\Omega t+ φ_n) (2)
已知f(t)=2A0+n=1∑∞Ancos(nΩt+φn)(2)
对其使用欧拉公式展开,将右边-n转换n:
f
(
t
)
=
A
0
2
+
∑
n
=
1
∞
A
n
2
[
e
j
(
n
Ω
t
)
+
φ
n
+
e
−
j
(
n
Ω
t
)
+
φ
n
]
=
A
0
2
+
1
2
∑
n
=
1
∞
A
n
e
j
(
n
Ω
t
)
e
φ
n
−
1
2
∑
n
=
1
∞
e
j
(
n
Ω
t
)
e
φ
n
=
1
2
∑
n
=
−
∞
∞
A
n
e
j
φ
n
e
j
n
Ω
t
\begin{aligned} f(t) & =\frac{A_0}{2}+\sum_{n=1}^{\infty}\frac{A_n}{2}[e^{j(n\Omega t)+\varphi_n}+e^{-j(n\Omega t)+\varphi_n}]\\ & =\frac{A_0}{2}+\frac{1}{2}\sum_{n=1}^{\infty}A_ne^{j(n\Omega t)}e^{\varphi_n}-\frac{1}{2}\sum_{n=1}^{\infty}e^{j(n\Omega t)}e^{\varphi_n}\\ & = \frac{1}{2}\sum_{n=-\infty}^\infty A_ne^{j\varphi_n}e^{jn\Omega t} \end{aligned}
f(t)=2A0+n=1∑∞2An[ej(nΩt)+φn+e−j(nΩt)+φn]=2A0+21n=1∑∞Anej(nΩt)eφn−21n=1∑∞ej(nΩt)eφn=21n=−∞∑∞AnejφnejnΩt
令复数
1
2
A
n
e
j
φ
n
=
∣
F
n
∣
e
j
φ
n
=
F
n
\frac{1}{2}A_ne^{j\varphi_n}=\left| F_n\right|e^{j\varphi_n}=F_n
21Anejφn=∣Fn∣ejφn=Fn
F n = 1 2 A n e j φ n = 1 2 ( A n c o s φ n + j A n s i n φ n ) = 1 2 ( a n − j b n ) = 1 T ∫ − T 2 T 2 f ( t ) c o s ( n Ω t ) d t − j 1 T ∫ − T 2 T 2 f ( t ) s i n ( n Ω t ) d t = 1 T ∫ − T 2 T 2 f ( t ) e − j n Ω t d t \begin{aligned} F_n&=\frac{1}{2}A_ne^{j\varphi_n}\\ &=\frac{1}{2}(A_ncos\varphi_n+jA_nsin\varphi_n)\\ &=\frac{1}{2}(a_n-jb_n)\\ &=\frac{1}{T}\int_{-\frac{T}{2}}^{\frac{T}{2}}f(t)cos(n\Omega t)dt-j\frac{1}{T}\int_{-\frac{T}{2}}^{\frac{T}{2}}f(t)sin(n\Omega t)dt\\ &=\frac{1}{T}\int_{-\frac{T}{2}}^{\frac{T}{2}}f(t)e^{-jn\Omega t}dt \end{aligned} Fn=21Anejφn=21(Ancosφn+jAnsinφn)=21(an−jbn)=T1∫−2T2Tf(t)cos(nΩt)dt−jT1∫−2T2Tf(t)sin(nΩt)dt=T1∫−2T2Tf(t)e−jnΩtdt
3 参考网上在线资料
[1]保姆级推导傅里叶级数
[2]通俗易懂的理解傅里叶变换(一)[收藏]
[3]傅里叶级数推导与理解
[4]傅里叶系列(一)傅里叶级数的推导
[5]西安电子科技大学——信号与系统(郭宝龙)