傅里叶变换
一、傅里叶级数
傅里叶级数是通过三角函数和常数项来叠加逼近周期为T的函数f(x)
1、定义
一个以T为周期的函数
f
T
(
x
)
f_T(x)
fT(x),如果在
[
−
T
2
,
T
2
]
[\frac{-T}{2},\frac{T}{2}]
[2−T,2T]上满足狄利克雷条件,即在
[
−
T
2
,
T
2
]
[\frac{-T}{2},\frac{T}{2}]
[2−T,2T]上满足:1.
f
T
(
x
)
f_T(x)
fT(x)连续或只有有限个第一类间断点;2.只有有限个极值点。那么在
[
−
T
2
,
T
2
]
[\frac{-T}{2},\frac{T}{2}]
[2−T,2T]上可写成傅里叶级数:
f
T
(
t
)
=
a
0
2
+
∑
n
=
1
∞
(
a
n
c
o
s
n
ω
0
t
+
b
n
s
i
n
n
ω
0
t
)
(1.1)
f_T(t)=\frac{a_0}{2}+\sum_{n=1}^{\infty}(a_ncosn\omega_0t+b_nsinn\omega_0t)\tag{1.1}
fT(t)=2a0+n=1∑∞(ancosnω0t+bnsinnω0t)(1.1)
其中
ω
0
=
2
π
T
\omega_0=\frac{2\pi}{T}
ω0=T2π,
a
0
=
2
T
∫
−
T
2
T
2
f
T
(
t
)
d
t
a_0=\frac{2}{T}\int_{-\frac{T}{2}}^\frac{T}{2}f_T(t)dt
a0=T2∫−2T2TfT(t)dt
a n = 2 T ∫ − T 2 T 2 f T ( t ) c o s n ω 0 t d t , ( n = 1 , 2 , 3 , ⋯ ) a_n=\frac{2}{T}\int_{-\frac{T}{2}}^\frac{T}{2}f_T(t)cosn\omega_0tdt, (n=1,2,3,\cdots) an=T2∫−2T2TfT(t)cosnω0tdt,(n=1,2,3,⋯)
b n = 2 T ∫ − T 2 T 2 f T ( t ) s i n n ω 0 t d t , ( n = 1 , 2 , 3 , ⋯ ) b_n=\frac{2}{T}\int_{-\frac{T}{2}}^\frac{T}{2}f_T(t)sinn\omega_0tdt, (n=1,2,3,\cdots) bn=T2∫−2T2TfT(t)sinnω0tdt,(n=1,2,3,⋯)
第一类间断点:①左右极限相等,但不等于该点的函数值 f ( x 0 ) f(x_0) f(x0)或者该点无意义。此时称为可去间断点。②左右极限在该点不相等,此时称为跳跃间断点。
2、变形
利用Euler公式: c o s t = 1 2 ( e j t + e − j t ) , s i n t = 1 2 ( e j t − e − j t ) cost=\frac{1}{2}(e^{jt}+e^{-jt}), sint=\frac{1}{2}(e^{jt}-e^{-jt}) cost=21(ejt+e−jt),sint=21(ejt−e−jt)
则(1.1)可以变形为:
f
T
(
t
)
=
a
0
2
+
∑
n
=
1
∞
(
a
n
−
j
b
n
2
e
j
n
ω
0
t
+
a
n
+
j
b
n
2
e
−
j
n
ω
0
t
)
(1.2)
\begin{equation} f_T(t)=\frac{a_0}{2}+\sum_{n=1}^{\infty}(\frac{a_n-jb_n}{2}e^{jn\omega_0t}+\frac{a_n+jb_n}{2}e^{-jn\omega_0t}) \end{equation}\tag{1.2}
fT(t)=2a0+n=1∑∞(2an−jbnejnω0t+2an+jbne−jnω0t)(1.2)
令
C
0
=
a
0
2
=
∫
−
T
2
T
2
f
T
(
t
)
d
t
C
n
=
a
n
−
j
b
n
2
=
1
T
[
∫
−
T
2
T
2
f
T
(
t
)
c
o
s
n
ω
0
t
d
t
−
j
∫
−
T
2
T
2
f
T
(
t
)
s
i
n
n
ω
0
t
d
t
]
=
1
T
∫
−
T
2
T
2
f
T
(
t
)
e
−
j
n
ω
0
t
d
t
(
n
=
1
,
2
,
3
,
⋯
)
C
−
n
=
a
n
+
j
b
n
2
=
1
T
∫
−
T
2
T
2
f
T
(
t
)
e
j
n
ω
0
t
d
t
(
n
=
1
,
2
,
3
,
⋯
)
C_0=\frac{a_0}{2}=\int_{-\frac{T}{2}}^\frac{T}{2}f_T(t)dt\\ C_n=\frac{a_n-jb_n}{2}=\frac{1}{T}[\int_{-\frac{T}{2}}^\frac{T}{2}f_T(t)cosn\omega_0tdt-j\int_{-\frac{T}{2}}^\frac{T}{2}f_T(t)sinn\omega_0tdt]=\frac{1}{T}\int_{-\frac{T}{2}}^\frac{T}{2}f_T(t)e^{-jn\omega_0t}dt (n=1,2,3,\cdots)\\ C_{-n}=\frac{a_n+jb_n}{2}=\frac{1}{T}\int_{-\frac{T}{2}}^\frac{T}{2}f_T(t)e^{jn\omega_0t}dt (n=1,2,3,\cdots)
C0=2a0=∫−2T2TfT(t)dtCn=2an−jbn=T1[∫−2T2TfT(t)cosnω0tdt−j∫−2T2TfT(t)sinnω0tdt]=T1∫−2T2TfT(t)e−jnω0tdt(n=1,2,3,⋯)C−n=2an+jbn=T1∫−2T2TfT(t)ejnω0tdt(n=1,2,3,⋯)
则
C
n
=
1
T
∫
−
T
2
T
2
f
T
(
t
)
e
−
j
n
ω
0
t
d
t
(
n
=
±
1
,
±
2
,
±
3
,
⋯
)
C_n=\frac{1}{T}\int_{-\frac{T}{2}}^\frac{T}{2}f_T(t)e^{-jn\omega_0t}dt (n=\pm1,\pm2,\pm3,\cdots)
Cn=T1∫−2T2TfT(t)e−jnω0tdt(n=±1,±2,±3,⋯)
则(1.2)可写为:
f
T
(
t
)
=
C
0
+
∑
n
=
1
∞
(
C
n
e
j
n
ω
0
t
+
C
−
n
e
−
j
n
ω
0
t
)
=
∑
n
=
−
∞
+
∞
C
n
e
j
n
ω
0
t
=
1
T
∑
n
=
−
∞
+
∞
[
∫
−
T
2
T
2
f
T
(
τ
)
e
−
j
n
ω
0
τ
d
τ
]
e
j
n
ω
0
t
(1.3)
f_T(t)=C_0+\sum_{n=1}^{\infty}(C_ne^{jn\omega_0t}+C_{-n}e^{-jn\omega_0t})=\sum_{n=-\infty}^{+\infty}C_ne^{jn\omega_0t}=\frac{1}{T}\sum_{n=-\infty}^{+\infty}[\int_{-\frac{T}{2}}^\frac{T}{2}f_T(\tau)e^{-jn\omega_0\tau}d\tau]e^{jn\omega_0t}\tag{1.3}
fT(t)=C0+n=1∑∞(Cnejnω0t+C−ne−jnω0t)=n=−∞∑+∞Cnejnω0t=T1n=−∞∑+∞[∫−2T2TfT(τ)e−jnω0τdτ]ejnω0t(1.3)
由于
C
n
C_n
Cn可以看出仅与
n
ω
0
n\omega_0
nω0有关,因此可以将
C
n
C_n
Cn记为
C
n
=
F
(
n
ω
0
)
C_n=F(n\omega_0)
Cn=F(nω0)
那么
f
T
(
t
)
=
∑
n
=
−
∞
+
∞
C
n
e
j
n
ω
0
t
=
∑
n
=
−
∞
+
∞
F
(
n
ω
0
)
j
n
ω
0
t
(1.4)
\begin{equation} f_T(t)=\sum_{n=-\infty}^{+\infty}C_ne^{jn\omega_0t}=\sum_{n=-\infty}^{+\infty}F(n\omega_0)^{jn\omega_0t} \end{equation}\tag{1.4}
fT(t)=n=−∞∑+∞Cnejnω0t=n=−∞∑+∞F(nω0)jnω0t(1.4)
(1.4)式说明以T为周期的振动
f
T
(
t
)
f_T(t)
fT(t)是简谐振动
F
(
n
ω
0
t
)
e
j
n
ω
0
t
(
n
=
0
,
±
1
,
±
2
,
⋯
)
F(n\omega_0t)e^{jn\omega_0t}(n=0,\pm1,\pm2,\cdots)
F(nω0t)ejnω0t(n=0,±1,±2,⋯)叠加产生的运动,其中
ω
0
\omega_0
ω0为基频,
F
(
n
ω
0
)
F(n\omega_0)
F(nω0)为
f
T
(
t
)
f_T(t)
fT(t)的离散频谱,
∣
F
(
n
ω
0
)
∣
\vert{F(n\omega_0)}\vert
∣F(nω0)∣为离散幅谱。
二、从傅里叶级数到傅里叶变换
1、傅里叶积分公式
傅里叶级数针对的对象是周期函数,然而对于非周期函数,我们可以这样看待:任意一个非周期函数 f ( t ) f(t) f(t)可以看成当周期函数 f T ( t ) f_T(t) fT(t)的周期T趋向于无穷时的情况,即 lim T → + ∞ f T ( t ) = f ( t ) \lim_{T\rightarrow{+\infty}}f_T(t)=f(t) limT→+∞fT(t)=f(t)
通过推导(在此省略)可以有
f
(
t
)
f(t)
f(t)的傅里叶积分公式:
f
(
t
)
=
1
2
π
∫
−
∞
+
∞
[
∫
−
∞
+
∞
f
(
τ
)
e
−
j
ω
τ
d
τ
]
e
j
ω
t
d
ω
(1.5)
\begin{equation} f(t)=\frac{1}{2\pi}\int_{-\infty}^{+\infty}[\int_{-\infty}^{+\infty}f(\tau)e^{-j\omega\tau}d\tau]e^{j\omega t}d\omega \end{equation}\tag{1.5}
f(t)=2π1∫−∞+∞[∫−∞+∞f(τ)e−jωτdτ]ejωtdω(1.5)
2、傅里叶积分定理
若 f ( t ) f(t) f(t)在 ( − ∞ , + ∞ ) (-\infty,+\infty) (−∞,+∞)上满足如下条件:
(1) f ( t ) f(t) f(t)在任一有限区间上满足狄利克雷条件;
(2) f ( t ) f(t) f(t)在无限区间 ( − ∞ , + ∞ ) (-\infty,+\infty) (−∞,+∞)上绝对可积(即积分 ∫ − ∞ + ∞ ∣ f ( t ) ∣ d t \int_{-\infty}^{+\infty}|{f(t)}|dt ∫−∞+∞∣f(t)∣dt收敛),则有(1.5)式成立,而在 f ( t ) f(t) f(t)的间断点处,(1.5)式左端应为 f ( t + 0 ) + f ( t − 0 ) 2 \frac{f(t+0)+f(t-0)}{2} 2f(t+0)+f(t−0)。
3、傅里叶变换的概念
令 F ( ω ) = ∫ − ∞ + ∞ f ( t ) e − j ω t d t F(\omega)=\int_{-\infty}^{+\infty}f(t)e^{-j\omega t}dt F(ω)=∫−∞+∞f(t)e−jωtdt; f ( t ) = 1 2 π ∫ − ∞ + ∞ F ( ω ) e j ω t d ω f(t)=\frac{1}{2\pi}\int_{-\infty}^{+\infty}F(\omega)e^{j\omega t}d\omega f(t)=2π1∫−∞+∞F(ω)ejωtdω
则(1.5)式为傅里叶变换, F ( ω ) F(\omega) F(ω)为 f ( t ) f(t) f(t)的像函数,记为 F ( ω ) = Γ [ f ( t ) ] ; F(\omega)=\Gamma[f(t)]; F(ω)=Γ[f(t)];
f ( t ) f(t) f(t)为 F ( ω ) F(\omega) F(ω)的像原函数,记为 f ( t ) = Γ − 1 [ F ( ω ) ] f(t)=\Gamma^{-1}[F(\omega)] f(t)=Γ−1[F(ω)]
4、傅里叶变换的性质
(1)线性性质
Γ [ α f 1 ( t ) + β f 2 ( t ) ] = α Γ [ f 1 ( t ) ] + β Γ [ f 2 ( t ) ] = α F 1 ( ω ) + β F 2 ( ω ) \Gamma[\alpha f_1(t)+\beta f_2(t)]=\alpha \Gamma[f_1(t)]+\beta\Gamma[f_2(t)]=\alpha F_1(\omega)+\beta F_2(\omega) Γ[αf1(t)+βf2(t)]=αΓ[f1(t)]+βΓ[f2(t)]=αF1(ω)+βF2(ω)
Γ − 1 [ α F 1 ( ω ) + β F 2 ( ω ) ] = α f 1 ( t ) + β f 2 ( t ) \Gamma^{-1}[\alpha F_1(\omega)+\beta F_2(\omega)]=\alpha f_1(t)+\beta f_2(t) Γ−1[αF1(ω)+βF2(ω)]=αf1(t)+βf2(t)
(2)位移性质
Γ [ f ( t ± t 0 ) ] = e ± j ω t 0 Γ [ f ( t ) ] \Gamma[f(t\pm t_0)]=e^{\pm j\omega t_0}\Gamma[f(t)] Γ[f(t±t0)]=e±jωt0Γ[f(t)]
Γ − 1 [ F ( ω ∓ ω 0 ) ] = e ± j ω 0 t f ( t ) \Gamma^{-1}[F(\omega\mp\omega_0)]=e^{\pm j\omega_0 t}f(t) Γ−1[F(ω∓ω0)]=e±jω0tf(t)
说明:①时间函数
f
(
t
)
f(t)
f(t)沿
t
t
t轴向左或向右位移
t
0
t_0
t0后进行傅里叶变换等于
f
(
t
)
f(t)
f(t)的傅里叶变换乘以因子
e
ω
t
0
j
e^{\omega t_0j}
eωt0j或
e
−
ω
t
0
j
e^{-\omega t_0j}
e−ωt0j
②物理意义:当一个信号沿时间轴移动后,它的各频率成分的大小不变,但相位发生变化
(3)相似性质
Γ [ f ( k t ) ] = 1 ∣ k ∣ F ( ω k ) \Gamma[f(kt)]=\frac{1}{|{k}|}F(\frac{\omega}{k}) Γ[f(kt)]=∣k∣1F(kω), Γ − 1 [ F ( k ω ) ] = 1 ∣ k ∣ f ( t k ) \Gamma^{-1}[F(k\omega)]=\frac{1}{|{k}|}f(\frac{t}{k}) Γ−1[F(kω)]=∣k∣1f(kt) 其中k为非零实数
说明:对时间轴进行拉伸时,频谱也会相应进行拉伸,表现为相似性质
(4)微分性质
如果当 ∣ t ∣ → ∞ |{t}|\rightarrow\infty ∣t∣→∞时, f ( t ) → 0 f(t)\rightarrow0 f(t)→0,则 Γ [ f ′ ( t ) ] = j ω Γ [ f ( t ) ] \Gamma[f^{\prime}(t)]=j\omega\Gamma[f(t)] Γ[f′(t)]=jωΓ[f(t)]
推论:若 lim ∣ t ∣ → ∞ f ( k ) = 0 , k = 0 , 1 , 2 , ⋯ , n − 1 \lim_{|{t}|\rightarrow\infty}f^{(k)}=0,k=0,1,2,\cdots,n-1 lim∣t∣→∞f(k)=0,k=0,1,2,⋯,n−1
则有: Γ [ f n ( t ) ] = ( j ω ) n Γ [ f ( t ) ] \Gamma[f^{{n}}(t)]=(j\omega)^n\Gamma[f(t)] Γ[fn(t)]=(jω)nΓ[f(t)]
d n d ω n F ( ω ) = ( − j ) n Γ [ t n f ( t ) ] \frac{d^n}{d\omega^n}F(\omega)=(-j)^n\Gamma[t^nf(t)] dωndnF(ω)=(−j)nΓ[tnf(t)]
(5)积分性质
t → + ∞ t\rightarrow +\infty t→+∞时, ∫ − ∞ t f ( t ) d t → 0 \int_{-\infty}^tf(t)dt\rightarrow0 ∫−∞tf(t)dt→0,则 Γ [ ∫ − ∞ t f ( t ) d t ] = 1 j ω Γ [ f ( t ) ] \Gamma[\int_{-\infty}^tf(t)dt]=\frac{1}{j\omega}\Gamma[f(t)] Γ[∫−∞tf(t)dt]=jω1Γ[f(t)]
(6)乘积定理
若 F 1 ( ω ) = Γ [ f 1 ( t ) ] , F 2 ( ω ) = Γ [ f 2 ( t ) ] F_1(\omega)=\Gamma[f_1(t)],F_2(\omega)=\Gamma[f_2(t)] F1(ω)=Γ[f1(t)],F2(ω)=Γ[f2(t)],则
∫ − ∞ + ∞ f 1 ( t ) f 2 ( t ) d t = 1 2 π ∫ − ∞ + ∞ F 1 ( ω ) F 2 ( ω ) ‾ d ω = 1 2 π ∫ − ∞ + ∞ F 2 ( ω ) F 1 ( ω ) ‾ d ω \int_{-\infty}^{+\infty}f_1(t)f_2(t)dt=\frac{1}{2\pi}\int_{-\infty}^{+\infty}F_1(\omega)\overline{F_2(\omega)}d\omega=\frac{1}{2\pi}\int_{-\infty}^{+\infty}F_2(\omega)\overline{F_1(\omega)}d\omega ∫−∞+∞f1(t)f2(t)dt=2π1∫−∞+∞F1(ω)F2(ω)dω=2π1∫−∞+∞F2(ω)F1(ω)dω
其中 f 1 ( t ) , f 2 ( t ) f_1(t),f_2(t) f1(t),f2(t)为 t t t的实函数,而 F 1 ( ω ) ‾ , F 2 ( ω ) ‾ \overline{F_1(\omega)},\overline{F_2(\omega)} F1(ω),F2(ω)分别为 F 1 ( ω ) , F 2 ( ω ) F_1(\omega),F_2(\omega) F1(ω),F2(ω)的共轭函数。
(7)卷积与相关函数
设函数 f 1 ( t ) f_1(t) f1(t)与 f 2 ( t ) f_2(t) f2(t)在 ( − ∞ , + ∞ ) (-\infty,+\infty) (−∞,+∞)内有定义,若广义积分 ∫ − ∞ + ∞ f 1 ( τ ) f 2 ( t − τ ) d τ \int_{-\infty}^{+\infty}f_1(\tau)f_2(t-\tau)d\tau ∫−∞+∞f1(τ)f2(t−τ)dτ对任何实数 t t t收敛,则它在 ( − ∞ , + ∞ ) (-\infty,+\infty) (−∞,+∞)上定义了一个自变量为 t t t的函数,记为: f 1 ( t ) ∗ f 2 ( t ) = ∫ − ∞ + ∞ f 1 ( τ ) f 2 ( t − τ ) d τ f_1(t)*f_2(t)=\int_{-\infty}^{+\infty}f_1(\tau)f_2(t-\tau)d\tau f1(t)∗f2(t)=∫−∞+∞f1(τ)f2(t−τ)dτ,且有如下成立:
① ∣ f 1 ( t ) ∗ f 2 ( t ) ∣ ⩽ ∥ f 1 ( t ) ∣ ∗ ∣ f 2 ( t ) ∣ |{f_1(t)*f_2(t)}|\leqslant\|f_1(t)|*|{f_2(t)}| ∣f1(t)∗f2(t)∣⩽∥f1(t)∣∗∣f2(t)∣
② f 1 ( t ) ∗ f 2 ( t ) = f 2 ( t ) ∗ f 1 ( t ) f_1(t)*f_2(t)=f_2(t)*f_1(t) f1(t)∗f2(t)=f2(t)∗f1(t)
③ f 1 ( t ) ∗ [ f 2 ( t ) ∗ f 3 ( t ) ] = [ f 1 ( t ) ∗ f 2 ( t ) ] ∗ f 3 ( t ) f_1(t)*[f_2(t)*f_3(t)]=[f_1(t)*f_2(t)]*f_3(t) f1(t)∗[f2(t)∗f3(t)]=[f1(t)∗f2(t)]∗f3(t)
④ f 1 ( t ) ∗ [ f 2 ( t ) + f 3 ( t ) ] = f 1 ( t ) ∗ f 2 ( t ) + f 1 ( t ) ∗ f 3 ( t ) f_1(t)*[f_2(t)+f_3(t)]=f_1(t)*f_2(t)+f_1(t)*f_3(t) f1(t)∗[f2(t)+f3(t)]=f1(t)∗f2(t)+f1(t)∗f3(t)
(8)卷积定理
假定 f 1 ( t ) f_1(t) f1(t)和 f 2 ( t ) f_2(t) f2(t)都满足傅里叶积分定理中的条件,且 Γ [ f 1 ( t ) ] = F 1 ( ω ) , Γ [ f 2 ( t ) ] = F 2 ( ω ) \Gamma[f_1(t)]=F_1(\omega),\Gamma[f_2(t)]=F_2(\omega) Γ[f1(t)]=F1(ω),Γ[f2(t)]=F2(ω),则有
Γ [ f 1 ( t ) ∗ f 2 ( t ) ] = F 1 ( ω ) ⋅ F 2 ( ω ) \Gamma[f_1(t)*f_2(t)]=F_1(\omega)\cdot F_2(\omega) Γ[f1(t)∗f2(t)]=F1(ω)⋅F2(ω),或 Γ − 1 [ F 1 ( ω ) ⋅ F 2 ( ω ) ] = f 1 ( t ) ∗ f 2 ( t ) \Gamma^{-1}[F_1(\omega)\cdot F_2(\omega)]=f_1(t)*f_2(t) Γ−1[F1(ω)⋅F2(ω)]=f1(t)∗f2(t)
同理可得 Γ [ f 1 ( t ) ⋅ f 2 ( t ) ] = 1 2 π F 1 ( ω ) ∗ F 2 ( ω ) \Gamma[f_1(t)\cdot f_2(t)]=\frac{1}{2\pi}F_1(\omega)*F_2(\omega) Γ[f1(t)⋅f2(t)]=2π1F1(ω)∗F2(ω)