矩形函数:
Π ( t ) = { 1 , ∣ t ∣ < 1 2 0 , ∣ t ∣ ≥ 1 2 \Pi(t) =\begin{cases} 1, &\left\lvert t \right\rvert<\frac{1}{2}\\ 0, &\left\lvert t \right\rvert\geq\frac{1}{2} \end{cases} Π(t)={1,0,∣t∣<21∣t∣≥21
对其进行傅里叶变换,我们有(需要运用欧拉公式的变形公式,详见傅里叶级数推导过程(2)、(3)式):
F Π ( s ) = ∫ − ∞ ∞ e − 2 π i s t Π ( t ) d t = ∫ − 1 2 1 2 e − 2 π i s t d t = − 1 2 π i s e − 2 π i s t ∣ t = − 1 2 t = 1 2 = 1 π s ( e π i s − e − π i s ) 2 i = sin ( π s ) π s \begin{aligned} \mathscr{F}\Pi (s) &= \int_{-\infty}^\infty e^{-2\pi ist} \Pi (t)dt\\ &= \int_{-\frac{1}{2}}^{\frac{1}{2}} e^{-2\pi ist}dt\\ &= \left. -\frac{1}{2\pi is} e^{-2\pi ist}\right | _{t=-\frac{1}{2}}^{t=\frac{1}{2}}\\ &=\frac{1}{\pi s} \frac{(e^{\pi is} - e^{-\pi is})}{2i}\\ &=\frac{\sin(\pi s)}{\pi s} \end{aligned} FΠ(s)=∫−∞∞e−2πistΠ(t)dt=∫−2121e−2πistdt=−2πis1e−2πist∣∣∣∣t=−21t=21=πs12i(eπis−e−πis)=πssin(πs)
我们定义
s
i
n
c
(
x
)
=
sin
(
π
x
)
π
x
sinc(x) = \frac{\sin(\pi x)}{\pi x}
sinc(x)=πxsin(πx),我们称其为归一化
s
i
n
c
sinc
sinc 函数,在数字信号处理中常用;
而对于
s
i
n
c
(
x
)
=
sin
(
x
)
x
sinc(x) = \frac{\sin(x)}{x}
sinc(x)=xsin(x),我们称其为未归一化的
s
i
n
c
sinc
sinc 函数,在数学中常用。
二者十分类似,图像如下:
故
F Π = s i n c \huge\mathscr{F}\Pi = sinc FΠ=sinc
观察推导过程,不难发现对于傅里叶逆变换,我们可以得到相同的结果,即
F − 1 Π = s i n c \huge\mathscr{F}^{-1}\Pi = sinc F−1Π=sinc
更一般地,设
Π p ( t ) = { 1 , ∣ t ∣ < p 2 0 , ∣ t ∣ ≥ p 2 \Pi_p(t) = \begin{cases} 1, &\left\lvert t \right\rvert<\frac{p}{2}\\ 0, &\left\lvert t \right\rvert\geq\frac{p}{2} \end{cases} Πp(t)={1,0,∣t∣<2p∣t∣≥2p
则有(证明类似,略):
F Π p ( s ) = sin ( π p s ) π s = p ⋅ s i n c ( p s ) \huge\mathscr{F}\Pi_p (s) = \frac{\sin \left( \pi ps \right)}{\pi s} = p \cdot sinc \left( ps \right) FΠp(s)=πssin(πps)=p⋅sinc(ps)
F − 1 Π p ( t ) = sin ( π p t ) π t = p ⋅ s i n c ( p t ) \huge\mathscr{F}^{-1}\Pi_p (t) = \frac{\sin \left( \pi pt \right)}{\pi t} = p \cdot sinc \left( pt \right) F−1Πp(t)=πtsin(πpt)=p⋅sinc(pt)
三角函数
Λ ( t ) = { 1 − ∣ t ∣ , ∣ t ∣ ≤ 1 0 , ∣ t ∣ > 1 \Lambda (t) = \begin{cases} 1-\left\lvert t \right\rvert , &\left\lvert t \right\rvert \leq 1\\ 0 , &\left\lvert t \right\rvert > 1 \end{cases} Λ(t)={1−∣t∣,0,∣t∣≤1∣t∣>1
对其进行傅里叶变换,我们有(需要运用分部积分法,详见微积分相关内容):
F Λ ( s ) = ∫ − ∞ ∞ e − 2 π i s t Λ ( t ) d t = ∫ − 1 0 e − 2 π i s t ( 1 + t ) d t + ∫ 0 1 e − 2 π i s t ( 1 − t ) d t = ∫ − 1 0 ( 1 + t ) d ( e − 2 π i s t − 2 π i s ) + ∫ 0 1 ( 1 − t ) d ( e − 2 π i s t − 2 π i s ) = ( 1 + t ) ( e − 2 π i s t − 2 π i s ) ∣ − 1 0 − ∫ − 1 0 ( e − 2 π i s t − 2 π i s ) d ( 1 + t ) + ( 1 − t ) ( e − 2 π i s t − 2 π i s ) ∣ 0 1 − ∫ 0 1 ( e − 2 π i s t − 2 π i s ) d ( 1 − t ) = ( 1 + t ) ( e − 2 π i s t − 2 π i s ) ∣ − 1 0 − ∫ 0 1 ( e − 2 π i s ( k − 1 ) − 2 π i s ) d k + ( 1 − t ) ( e − 2 π i s t − 2 π i s ) ∣ 0 1 − ∫ 1 0 ( e − 2 π i s ( 1 − k ) − 2 π i s ) d k = ( 1 + t ) ( e − 2 π i s t − 2 π i s ) ∣ − 1 0 − ( e 2 π i s − 2 π i s e − 2 π i s k − 2 π i s ) ∣ 0 1 + ( 1 − t ) ( e − 2 π i s t − 2 π i s ) ∣ 0 1 − ( e − 2 π i s − 2 π i s e 2 π i s k 2 π i s ) ∣ 1 0 = 2 − e 2 π i s − e − 2 π i s 4 π 2 s 2 \begin{aligned} \mathscr{F}\Lambda (s) &= \int_{-\infty}^\infty e^{-2\pi ist} \Lambda (t)dt\\ &= \int_{-1}^{0} e^{-2\pi ist}(1+t)dt + \int_{0}^{1} e^{-2\pi ist}(1-t)dt\\ &= \int_{-1}^{0} (1+t)d(\frac{e^{-2\pi ist}}{-2\pi is}) + \int_{0}^{1} (1-t)d(\frac{e^{-2\pi ist}}{-2\pi is})\\ &= \left. (1+t)(\frac{e^{-2\pi ist}}{-2\pi is})\right |_{-1}^0 - \int_{-1}^{0} (\frac{e^{-2\pi ist}}{-2\pi is})d(1+t) + \left. (1-t)(\frac{e^{-2\pi ist}}{-2\pi is})\right |_0^1 - \int_0^1 (\frac{e^{-2\pi ist}}{-2\pi is})d(1-t)\\ &= \left. (1+t)(\frac{e^{-2\pi ist}}{-2\pi is})\right |_{-1}^0 - \int_0^1 (\frac{e^{-2\pi is(k-1)}}{-2\pi is})dk + \left. (1-t)(\frac{e^{-2\pi ist}}{-2\pi is})\right |_0^1 - \int_1^0 (\frac{e^{-2\pi is(1-k)}}{-2\pi is})dk\\ &= \left. (1+t)(\frac{e^{-2\pi ist}}{-2\pi is})\right |_{-1}^0 - \left.(\frac{e^{2\pi is}}{-2\pi is} \frac{e^{-2\pi isk}}{-2\pi is}) \right |_0^1 + \left. (1-t)(\frac{e^{-2\pi ist}}{-2\pi is})\right |_0^1 - \left.(\frac{e^{-2\pi is}}{-2\pi is} \frac{e^{2\pi isk}}{2\pi is}) \right |_1^0\\ &=\frac{2-e^{2\pi is}- e^{-2\pi is}}{4\pi ^2 s^2} \end{aligned} FΛ(s)=∫−∞∞e−2πistΛ(t)dt=∫−10e−2πist(1+t)dt+∫01e−2πist(1−t)dt=∫−10(1+t)d(−2πise−2πist)+∫01(1−t)d(−2πise−2πist)=(1+t)(−2πise−2πist)∣∣∣∣−10−∫−10(−2πise−2πist)d(1+t)+(1−t)(−2πise−2πist)∣∣∣∣01−∫01(−2πise−2πist)d(1−t)=(1+t)(−2πise−2πist)∣∣∣∣−10−∫01(−2πise−2πis(k−1))dk+(1−t)(−2πise−2πist)∣∣∣∣01−∫10(−2πise−2πis(1−k))dk=(1+t)(−2πise−2πist)∣∣∣∣−10−(−2πise2πis−2πise−2πisk)∣∣∣∣01+(1−t)(−2πise−2πist)∣∣∣∣01−(−2πise−2πis2πise2πisk)∣∣∣∣10=4π2s22−e2πis−e−2πis
利用欧拉公式,我们有:
F Λ ( s ) = 2 − ( cos 2 π s + i sin 2 π s ) − ( cos 2 π s − i sin 2 π s ) 4 π 2 s 2 = 2 − 2 cos 2 π s 4 π 2 s 2 = 2 − 2 ( 1 − 2 sin 2 π s ) 4 π 2 s 2 = sin 2 π s π 2 s 2 = s i n c 2 ( s ) \begin{aligned} \mathscr{F}\Lambda (s) &= \frac{2-(\cos2\pi s +i \sin2\pi s) -(\cos2\pi s - i \sin2\pi s)}{4\pi ^2 s^2} \\ &=\frac{2-2\cos2\pi s}{4\pi ^2 s^2} \\ &=\frac{2-2(1-2\sin^2 \pi s)}{4\pi ^2 s^2} \\ &=\frac{\sin^2 \pi s}{\pi ^2 s^2} \\ &=sinc^2 (s) \end{aligned} FΛ(s)=4π2s22−(cos2πs+isin2πs)−(cos2πs−isin2πs)=4π2s22−2cos2πs=4π2s22−2(1−2sin2πs)=π2s2sin2πs=sinc2(s)
即
F Λ = s i n c 2 \huge\mathscr{F}\Lambda = sinc^2 FΛ=sinc2
最后我们来看高斯函数的傅里叶变换
高斯分布(又称正态分布)公式如下:
f ( x ) = 1 2 π σ e − ( x − μ ) 2 2 σ 2 f(x) = \frac{1}{\sqrt{2\pi}\sigma} e^{-\frac{(x-\mu)^2}{2\sigma^2}} f(x)=2πσ1e−2σ2(x−μ)2
即
X
X
X ~
N
(
μ
,
σ
2
)
N(\mu, \sigma^2)
N(μ,σ2)
我们令
μ
=
0
\mu = 0
μ=0,
σ
2
=
1
2
π
\sigma^2 = \frac{1}{2\pi}
σ2=2π1,用
t
t
t 代替
x
x
x 成为自变量,则有以下高斯函数特例:
f ( t ) = e − π t 2 (1) f(t) = e^{-\pi t^2}\tag1 f(t)=e−πt2(1)
对 ( 1 ) (1) (1) 式进行傅里叶变换,为了方便,我们将 F f \mathscr{F}f Ff 写作 F F F,则有
F ( s ) = F f ( s ) = ∫ − ∞ + ∞ e − 2 π i s t f ( t ) d t (2) F(s) = \mathscr{F}f(s) = \int_{-\infty}^{+\infty} e^{-2\pi ist}f(t)dt \tag2 F(s)=Ff(s)=∫−∞+∞e−2πistf(t)dt(2)
两边同时对 s s s求导,
F ′ ( s ) = ∫ − ∞ + ∞ d d s ( e − 2 π i s t ) e − π t 2 d t = ∫ − ∞ + ∞ ( − 2 π i t ) e − 2 π i s t e − π t 2 d t = i ∫ − ∞ + ∞ e − 2 π i s t d ( e − π t 2 ) = i e − 2 π i s t e − π t 2 ∣ t = − ∞ t = + ∞ − i ∫ − ∞ ∞ e − π t 2 d ( e − 2 π i s t ) (3) \begin{aligned} F^\prime (s) &= \int_{-\infty}^{+\infty} \frac{d}{ds}(e^{-2\pi ist}) e^{-\pi t^2}dt\\ &= \int_{-\infty}^{+\infty} (-2\pi it) e^{-2\pi ist} e^{-\pi t^2}dt\\ &=i\int_{-\infty}^{+\infty} e^{-2\pi ist}d(e^{-\pi t^2})\\ &=i\left. e^{-2\pi ist}e^{-\pi t^2}\right|_{t=-\infty}^{t=+\infty}-i\int_{-\infty}^\infty e^{-\pi t^2}d(e^{-2\pi ist})\tag3 \end{aligned} F′(s)=∫−∞+∞dsd(e−2πist)e−πt2dt=∫−∞+∞(−2πit)e−2πiste−πt2dt=i∫−∞+∞e−2πistd(e−πt2)=ie−2πiste−πt2∣∣∣t=−∞t=+∞−i∫−∞∞e−πt2d(e−2πist)(3)
对于
(
2
)
(2)
(2) 式,我们将等式右边拆分成两个部分进行计算
对于第一部分,我们不难发现,
i
e
−
2
π
i
s
t
=
i
cos
(
2
π
s
t
)
+
sin
(
π
s
t
)
ie^{-2\pi ist}=i\cos(2\pi st) + \sin (\pi st)
ie−2πist=icos(2πst)+sin(πst),由于
sin
x
\sin x
sinx、
cos
x
\cos x
cosx 是有限的,即取值范围都在
[
−
1
,
1
]
[-1,1]
[−1,1],因此,
i
e
−
2
π
i
s
t
ie^{-2\pi ist}
ie−2πist 的取值也是有限的,即在复平面内模长为1的圆上取值;而当
t
→
∞
t\to \infty
t→∞ 时
e
−
π
t
2
→
0
e^{-\pi t^2}\to 0
e−πt2→0。一个有限的数与零相乘,其结果为0,故
i e − 2 π i s t e − π t 2 ∣ t = − ∞ t = + ∞ = 0 (4) i\left. e^{-2\pi ist}e^{-\pi t^2}\right|_{t=-\infty}^{t=+\infty} = 0\tag4 ie−2πiste−πt2∣∣∣t=−∞t=+∞=0(4)
第二部分:
i ∫ − ∞ + ∞ e − π t 2 d ( e − 2 π i s t ) = i ∫ − ∞ + ∞ − 2 π i s e − 2 π i s t e − π t 2 d t = 2 π s ∫ − ∞ + ∞ e − 2 π i s t e − π t 2 d t = 2 π s F ( s ) (5) \begin{aligned} i\int_{-\infty}^{+\infty} e^{-\pi t^2}d(e^{-2\pi ist}) &= i\int_{-\infty}^{+\infty} -2\pi is e^{-2\pi ist} e^{-\pi t^2}dt\\ &=2\pi s \int_{-\infty}^{+\infty} e^{-2\pi ist} e^{-\pi t^2}dt\\ &=2\pi s F(s) \tag5 \end{aligned} i∫−∞+∞e−πt2d(e−2πist)=i∫−∞+∞−2πise−2πiste−πt2dt=2πs∫−∞+∞e−2πiste−πt2dt=2πsF(s)(5)
由 ( 3 ) (3) (3)、 ( 4 ) (4) (4)、 ( 5 ) (5) (5)式,我们可以得到:
F ′ ( s ) = 0 − 2 π s F ( s ) = − 2 π s F ( s ) F^\prime (s) = 0 - 2\pi s F(s) = - 2\pi s F(s) F′(s)=0−2πsF(s)=−2πsF(s)
运用分离变量法解该常微分方程,我们可以得到:
F ( s ) = F ( 0 ) e − π s 2 (6) F(s) = F(0) e^{-\pi s^2}\tag6 F(s)=F(0)e−πs2(6)
根据 ( 2 ) (2) (2) 式,我们可以得到:
F ( 0 ) = ∫ − ∞ + ∞ e − 2 π i 0 t f ( t ) d t = ∫ − ∞ + ∞ e − π t 2 d t (7) F(0) = \int_{-\infty}^{+\infty} e^{-2\pi i 0 t}f(t)dt = \int_{-\infty}^{+\infty} e^{-\pi t^2}dt\tag7 F(0)=∫−∞+∞e−2πi0tf(t)dt=∫−∞+∞e−πt2dt(7)
下面我们先来计算定积分
∫
−
∞
∞
e
−
π
t
2
d
t
\int_{-\infty}^\infty e^{-\pi t^2}dt
∫−∞∞e−πt2dt
令
g
(
x
)
=
∫
−
∞
∞
e
−
π
x
2
d
x
g(x) = \int_{-\infty}^\infty e^{-\pi x^2}dx
g(x)=∫−∞∞e−πx2dx,则
g 2 ( x ) = ( ∫ − ∞ + ∞ e − π x 2 d x ) 2 = ∫ − ∞ + ∞ e − π x 2 d x ∫ − ∞ + ∞ e − π y 2 d y = ∫ − ∞ + ∞ ∫ − ∞ + ∞ e − π x 2 e − π y 2 d x d y = ∫ − ∞ + ∞ ∫ − ∞ + ∞ e − π ( x 2 + y 2 ) d x d y \begin{aligned} g^2(x) &= (\int_{-\infty}^{+\infty} e^{-\pi x^2}dx)^2\\ &=\int_{-\infty}^{+\infty} e^{-\pi x^2}dx \int_{-\infty}^{+\infty} e^{-\pi y^2}dy\\ &=\int_{-\infty}^{+\infty} \int_{-\infty}^{+\infty} e^{-\pi x^2}e^{-\pi y^2}dxdy\\ &=\int_{-\infty}^{+\infty} \int_{-\infty}^{+\infty} e^{-\pi (x^2+y^2)}dxdy \end{aligned} g2(x)=(∫−∞+∞e−πx2dx)2=∫−∞+∞e−πx2dx∫−∞+∞e−πy2dy=∫−∞+∞∫−∞+∞e−πx2e−πy2dxdy=∫−∞+∞∫−∞+∞e−π(x2+y2)dxdy
我们对 x x x、 y y y 进行极坐标变换,可得:
g 2 ( x ) = ∫ 0 2 π ∫ 0 + ∞ e − π ρ 2 ρ d ρ d θ = ∫ 0 + ∞ 2 π e − π ρ 2 ρ d ρ = π ∫ 0 + ∞ e − π ρ 2 d ( ρ 2 ) = π ( 1 − π e − π ρ 2 ∣ 0 + ∞ ) = − ( 0 − 1 ) = 1 (8) \begin{aligned} g^2(x) &= \int_0^{2\pi} \int_0^{+\infty} e^{-\pi \rho^2} \rho d\rho d\theta \\ &= \int_0^{+\infty} 2\pi e^{-\pi \rho^2} \rho d\rho \\ &= \pi \int_0^{+\infty} e^{-\pi \rho^2} d(\rho^2) \\ &= \pi (\left. \frac{1}{-\pi} e^{-\pi \rho^2} \right|_0^{+\infty})\\ &= -(0-1)\\ &=1\tag8 \end{aligned} g2(x)=∫02π∫0+∞e−πρ2ρdρdθ=∫0+∞2πe−πρ2ρdρ=π∫0+∞e−πρ2d(ρ2)=π(−π1e−πρ2∣∣∣∣0+∞)=−(0−1)=1(8)
由 ( 6 ) (6) (6)、 ( 7 ) (7) (7)、 ( 8 ) (8) (8)式,我们有
F ( s ) = F f ( s ) = e − π s 2 (9) F(s) = \mathscr{F}f(s) = e^{-\pi s^2}\tag9 F(s)=Ff(s)=e−πs2(9)
这个结果表示
高斯函数 f ( t ) = e − π t 2 f(t) = e^{-\pi t^2} f(t)=e−πt2的傅里叶变换结果为 F f ( s ) = e − π s 2 \mathscr{F}f(s) = e^{-\pi s^2} Ff(s)=e−πs2,该图像与原图像形状是一样的
(本文中部分图片来源于wikipedia)