傅立叶变换通常产生复杂的光谱数据。在某些技术条件下,它们是双射。通过傅里叶变换,您可以唯一地恢复一个信号。但是,在查看频谱时,情况有所不同:信号可以具有相同的幅度谱和非常不同的相位,如下例所示:
作为双射,傅里叶变换不会丢失时间定位;然而,幅度谱不知何故。一个问题是,相位以一种通常不容易阅读或破译的方式对时间或空间定位进行编码。一种解释是,对于不太稳定的数据,相位可以非常快速地变化,并且在]−π,π],它看起来像一个被截断的信号,很难解开。这尤其适用于图像。下面是 Lenna 的一排,Cornouaille 船的一排。
第二列是幅度谱。摇摇欲坠,但可读。第三列是展开阶段。它似乎根本没有结构。然而,它编码了原始图像中的大部分变化。您可以通过交换模数和相位来测试它:这是第四列。如果你取Lenna光谱的模数和船的相位,做一个逆傅里叶变换,你会得到一艘幽灵船,几乎没有Lenna的痕迹。您可以在底行看到相反的情况。相位分量承载着很多信息,但如果不进行本地化,则不容易理解。
定位光谱的一种方法是使用窗口h,你很容易得到短时傅里叶变换:
Ss(τ,f;h)=∫s(t)h∗(t−τ)e−ı2πftdt
但如果你考虑h作为一个无限常数窗口,那么h∗(t−τ)是常数,你恢复标准傅里叶变换:
Ss(⋅,f;1)=∫s(t)e−ı2πftdt
但人们仍在试图理解复杂的相位模式,即使是短时傅里叶。窗口化并不能揭示一切。
[编辑] 在时间频率或时间尺度表示中找到最佳表示是一件困难的事情,原因有几个(并非详尽无遗):
- “最佳”通常与客观指标相关。传统的权衡(Weyl-Pauli-Heisenberg)被表述为L2范数,非常适合双时/频向量空间。然而,它并没有告诉你表示是如何简单、稀疏的,而且许多其他信息度量(熵、散度)在这方面提供了很多信息,但计算起来却非常困难。
- 所有数据的“最佳”意义不大,应该限制在更具体的数据模型和空间,并在这些类中推导出“最佳”:分段多项式、峰值总和......
- 数字信号需要对上述连续方案进行离散化,而离散化的方式会产生很大影响。
- 实际测量不服从标准假设:一切都不是线性的,实现不精确地遵循理论分布,短数据不具有渐近特性。
- 最终,理论上的最佳值并不是您可以轻松参数化以获得可重复结果的最佳值。
如果您的数据形态多样(如尖峰 + 振荡),您可能不得不尝试掌握不同的工具:框架、基数并集、数据转换等。