概述
傅里叶变换
一、傅里叶级数
傅里叶级数是通过三角函数和常数项来叠加逼近周期为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_ncosnomega_0t+b_nsinnomega_0t)tag{1.1}
fT(t)=2a0+n=1∑∞(ancosnω0t+bnsinnω0t)(1.1)
其中
ω
0
=
2
π
T
omega_0=frac{2pi}{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)cosnomega_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)sinnomega_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^{jnomega_0t}+frac{a_n+jb_n}{2}e^{-jnomega_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)cosnomega_0tdt-jint_{-frac{T}{2}}^frac{T}{2}f_T(t)sinnomega_0tdt]=frac{1}{T}int_{-frac{T}{2}}^frac{T}{2}f_T(t)e^{-jnomega_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^{jnomega_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^{-jnomega_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^{jnomega_0t}+C_{-n}e^{-jnomega_0t})=sum_{n=-infty}^{+infty}C_ne^{jnomega_0t}=frac{1}{T}sum_{n=-infty}^{+infty}[int_{-frac{T}{2}}^frac{T}{2}f_T(tau)e^{-jnomega_0tau}dtau]e^{jnomega_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
nomega_0
nω0有关,因此可以将
C
n
C_n
Cn记为
C
n
=
F
(
n
ω
0
)
C_n=F(nomega_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^{jnomega_0t}=sum_{n=-infty}^{+infty}F(nomega_0)^{jnomega_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(nomega_0t)e^{jnomega_0t}(n=0,pm1,pm2,cdots)
F(nω0t)ejnω0t(n=0,±1,±2,⋯)叠加产生的运动,其中
ω
0
omega_0
ω0为基频,
F
(
n
ω
0
)
F(nomega_0)
F(nω0)为
f
T
(
t
)
f_T(t)
fT(t)的离散频谱,
∣
F
(
n
ω
0
)
∣
vert{F(nomega_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_{Trightarrow{+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}{2pi}int_{-infty}^{+infty}[int_{-infty}^{+infty}f(tau)e^{-jomegatau}dtau]e^{jomega t}domega 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^{-jomega t}dt F(ω)=∫−∞+∞f(t)e−jωtdt; f ( t ) = 1 2 π ∫ − ∞ + ∞ F ( ω ) e j ω t d ω f(t)=frac{1}{2pi}int_{-infty}^{+infty}F(omega)e^{jomega t}domega 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)]+betaGamma[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(tpm t_0)]=e^{pm jomega 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(omegampomega_0)]=e^{pm jomega_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(komega)]=frac{1}{|{k}|}f(frac{t}{k}) Γ−1[F(kω)]=∣k∣1f(kt) 其中k为非零实数
说明:对时间轴进行拉伸时,频谱也会相应进行拉伸,表现为相似性质
(4)微分性质
如果当 ∣ t ∣ → ∞ |{t}|rightarrowinfty ∣t∣→∞时, f ( t ) → 0 f(t)rightarrow0 f(t)→0,则 Γ [ f ′ ( t ) ] = j ω Γ [ f ( t ) ] Gamma[f^{prime}(t)]=jomegaGamma[f(t)] Γ[f′(t)]=jωΓ[f(t)]
推论:若 lim ∣ t ∣ → ∞ f ( k ) = 0 , k = 0 , 1 , 2 , ⋯ , n − 1 lim_{|{t}|rightarrowinfty}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)]=(jomega)^nGamma[f(t)] Γ[fn(t)]=(jω)nΓ[f(t)]
d n d ω n F ( ω ) = ( − j ) n Γ [ t n f ( t ) ] frac{d^n}{domega^n}F(omega)=(-j)^nGamma[t^nf(t)] dωndnF(ω)=(−j)nΓ[tnf(t)]
(5)积分性质
t → + ∞ trightarrow +infty t→+∞时, ∫ − ∞ t f ( t ) d t → 0 int_{-infty}^tf(t)dtrightarrow0 ∫−∞tf(t)dt→0,则 Γ [ ∫ − ∞ t f ( t ) d t ] = 1 j ω Γ [ f ( t ) ] Gamma[int_{-infty}^tf(t)dt]=frac{1}{jomega}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}{2pi}int_{-infty}^{+infty}F_1(omega)overline{F_2(omega)}domega=frac{1}{2pi}int_{-infty}^{+infty}F_2(omega)overline{F_1(omega)}domega ∫−∞+∞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)dtau ∫−∞+∞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)dtau 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}{2pi}F_1(omega)*F_2(omega) Γ[f1(t)⋅f2(t)]=2π1F1(ω)∗F2(ω)
最后
以上就是哭泣纸鹤为你收集整理的傅里叶变换傅里叶变换二、从傅里叶级数到傅里叶变换的全部内容,希望文章能够帮你解决傅里叶变换傅里叶变换二、从傅里叶级数到傅里叶变换所遇到的程序开发问题。
如果觉得靠谱客网站的内容还不错,欢迎将靠谱客网站推荐给程序员好友。
发表评论 取消回复