傅里叶变换的两种直观理解

本文主要记录了两种对傅里叶变换的直观理解,并比较其相似之处。

在网上看了一些有关傅里叶变换的视频,其中以3b1b 的视频以及 Zach 的视频最为直观且易于理解。他们两人是从不同的角度出发来解释傅里叶变换的意义的,其不完全相同但又有共通之处,因此在这里记录下来。

经过两人的讲解,我发现傅里叶变换的实际意义:可以将其想象成一种频率分离器。它提取频率的方法在于,其对不同的频率响应不同。对实际信号中存在的频率,其会产生一个幅值很大的响应。对信号中不存在的频率,其响应幅度就会很小。将不同频率的响应幅值作图,就是我们常见的傅里叶变换幅值图。

准确地说,傅里叶变换后的结果是一个复数,也就是,每一个频率 ω\omega 对应一个复数,我们可以根据此复数计算出幅值和相位。

3b1b 的解释

先写上傅里叶变换的公式

傅里叶变换公式

从中可以看出,原信号 f(t)f(t) 乘了一个函数 ejωte^{-j\omega t}. 根据欧拉公式,可以将 ejωte^{-j\omega t} 看作以频率 ω\omega 进行顺时针旋转。因而 f(t)ejωtf(t)e^{-j\omega t} 可以被看作原信号 f(t)f(t) 被以频率 ω\omega 缠绕在单位圆上。这正是 3b1b 想要解释的。

图1 信号缠绕

图2 信号以0.5圈每秒的频率缠绕

图3 信号以0.2圈每秒的频率缠绕

从图2和图3可以看出,同样频率的信号,可以以不同的频率缠绕在一个圆上。这个缠绕频率与信号本身的频率无关。当缠绕频率较小时,一个圆上可以缠绕更多周期的信号,因此信号看起来比较“瘦”。与之相反,缠绕频率较大时,信号看起来更“胖”。

尝试将原信号看作一个有质量的铁丝,并追踪缠绕后图像的质心。当缠绕频率 ω\omega 变化时,缠绕图像变化,其质心也不可避免地变换。

图4 每秒1.5圈

图5 每秒3圈

对比 图4 和 图5 可以发现,当缠绕的频率等于信号本身的频率时,其质心的横坐标会明显偏右,因此在缠绕频率和质心的图像上出现一个明显的尖峰。而根据这些尖峰,我们就能找到原信号里的频率。

在实际的傅里叶变换中,F(ω)\textit{F}(\omega) 是一个复数,可以看作质心在平面上的坐标。上面我们只跟踪了质心的横坐标,但是这里,我们跟踪它的横纵坐标。

在这里,质心的计算是通过积分来体现的。

这样,整个傅里叶变换公式都好理解了,f(t)ejωtf(t) e^{-j\omega t} 表示信号以频率 ω\omega 缠绕在圆上,而积分则表示计算缠绕后图像的质心。质心是分布在平面上的。因此 F(ω)F(\omega) 的输入是频率 ω\omega, 输出是质心在平面上的位置,用一个复数来表示。

这里,我们计算质心的公式与实际的傅里叶变换公式有一点小区别,就是实际的计算公式在进行积分后就结束了,而实际的计算质心的公式多了个除以积分区间的部分,
质心计算公式和傅里叶变换公式

可以看出,质心计算公式比傅里叶变换公式多了个除以信号持续时间的部分。也就是说,傅里叶变换不仅追踪了质心,还将其根据持续时间进行了缩放。

Zach Star 的解释

Zach 首先将傅里叶变换公式中的 ejωte^{-j\omega t} 分解成 cos(ωt)isin(ωt)\cos(\omega t)- \textit{i}\sin(\omega t),并将其带入原公式,进行适当的变形。这样,公式被拆分成两部分,可以将积分看作函数 f(t)cos(ωt)\textit{f(t)}\cos(\omega t)f(t)sin(ωt)\textit{f(t)}\sin(\omega t) 与横坐标围成的面积。这样,对某个具体的 ω\omega 值,我们就能得到一个复数。
公式拆分成两部分

假设原信号 f(t)\textit{f(t)} 是一个方波函数,则其乘以 cos(t)\cos(t)sin(t)\sin(t) 后的结果如下图所示
 为方波,w 为 3 时

接下来,扫面每个可能的 ω\omega 值,并记录其对应的幅值。如下所示,
不同的 w

考虑另一种更复杂的函数, etcos(5t)\textit{e}^{-\textit{t}}\cos(5\textit{t}), 其幅值如下图所示
复杂函数的傅里叶变换的幅值

由于傅里叶变换有幅值和相位这两种指标,其各自代表的意义不太一样。其中,幅值可以近似看作两边的面积组合有多大,如下图所示
幅值的意义
相位图告诉我们当对而言,哪边的面积更大,以及其中任意一个是否为负数。

接下来,考虑两个 coscos 函数相乘。这里将 cos(x)\cos(x) 这种面积随着 xx 在 0 两侧波动的情况记为面积为 0. 因此,cos(πx)cos(3x)\cos(\pi x)\cos(3x) 这种两个 cos\cos 函数相乘的情况,其面积也是横轴上下互相抵消,记为 0 的。
cos(pi*x)cos(3x) 函数的面积

cos(pi*x)cos(3x) 函数的面积

只有在两个函数周期相同的时候,其与横轴围成的面积变成 \infty.
频率相同时,面积为无穷

因此,当我们求 cos(πx)\cos(\pi x) 函数的傅里叶变换 F(ω)\textit{F}(\omega) 时,其分为两部分,一部分为 cos(πx)cos(ω)\cos(\pi x)\cos(\omega), 另一部分为 icos(πx)sin(ω)\textit{i}\cos(\pi x)\sin(\omega)。可以发现,cos(πx)sin(ω)\cos(\pi x)\sin(\omega) 是关于原点对称的,因此其面积一直为 0. 这样,当 ω=π\omega=\pi 时,值会变成 \infty. 这样,我们就可以通过记录所有的 ω\omega 值以及其对应的面积,来分解出原信号中的频率。

这不仅对有限多个 cos\cos 函数组成的函数成立,其对更为复杂的函数同样成立。比如一个周期方波函数,
w 为信号中不存在的频率时,面积为 0

w 为 pi 时,面积为无穷(因为 pi 存在)

w 为 5pi 时,面积为无穷(因为 5pi 存在)

上面所讨论的输入都是周期函数。针对非周期函数,比如一个非周期方波函数,我们发现当 ω=3\omega=3 时,其面积为 0.66. 这里就与上面周期函数 0 或者 \infty 不同,(因为 cos(3π)cos(3π)\cos(3\pi)\cos(3\pi) 的面积是 \infty)其是一个有限的非零值。为了将面积做成有限的非零值,我们将 a1a_{1} 设为一个无限接近 0 的数,这样,a1cos(3t)cos(3t)a_{1}\cos(3t)\cos(3t) 的面积为 0.66. 这可以通过将 a1a_{1} 设为一个无限小的值而达到。同样的,对其他任意 ω\omega, 都会出现面积为一个有限非零值的情况,操作方式与 ω=3\omega=3 时一样。
面积为 0.66

拓展:拉普拉斯变换和傅里叶变换的关系

一句话概括,傅里叶变换就像一个频率扫描器,其扫描出信号中的组成频率。而拉普拉斯变换除了扫描频率,还扫描信号中的衰减指数。
分析的过程同样可以用面积来思考。
傅里叶变换和拉普拉斯变换


文章作者: taosean
文章链接: https://taosean.github.io/2021/12/19/Fourier-Transform/
版权声明: 本博客所有文章除特别声明外,均采用 CC BY-NC-SA 4.0 许可协议。转载请注明来自 taosean's 学习之旅