我正在尝试解决一个与本文中讨论的问题非常相似的问题
我有一个宽带信号,其中包含一个频率随时间变化的分量。我需要随着时间的推移监控这个组件的阶段。我能够通过频谱图中的峰值跟踪(一种有点蛮力的方法)来跟踪频移。我需要“清理”这个随时间变化的峰值附近的信号以提取希尔伯特相位(或者,我需要一种跟踪不涉及希尔伯特变换的相位的方法)。
总结上一篇文章:及时改变 FIR/IIR 滤波器的系数会导致坏事发生(它不仅会改变通带,还会以导致令人惊讶的瞬态的方式完全混淆滤波器状态)。但是,可能有一些方法可以及时调整滤波器系数(可能通过以某种智能方式联合修改滤波器系数和滤波器状态)。这超出了我的专业知识,但我愿意接受任何解决方案。
有两类解决方案似乎是合理的:一类是使用具有时变频率的谐振器滤波器(基本上是由信号驱动的阻尼谐波振荡器)。该模型非常简单,可以避免令人惊讶的滤波器瞬变。我会试试这个——但是谐振器在阻带中的衰减很差(如果它们甚至可以说有阻带?)。这让我很紧张,因为我不能 100% 确定共振滤波器的表现。
另一个建议是使用滤波器组并根据频率在各种带通滤波信号之间进行平滑插值。这种方法看起来很吸引人,但我怀疑它有一些隐藏的警告。我想线性混合两个带通滤波信号可能并不总是像你期望的那样,并且可能会导致奇怪的事情?但是,这不是我的专业领域,所以如果在过滤器组上混合被认为是一种安全的解决方案(之前已经分析和发布过的解决方案),我会使用它。
我想到了另一类潜在的解决方案,即仅从滑动短时傅里叶变换中的频率峰值获取相位(可以是加窗、多锥等)。如果有人知道这方面的任何先前文献,我会非常感兴趣。相关的,将是从感兴趣频带上的滑动复数 Morlet 小波变换中获取频率功率峰值处的相位。
所以,我想,基本上我想到了三类解决方案。1. 具有时变频率的谐振器滤波器。2. 使用过滤器组,可能与混合?3. 从 STFT 或 CWT 中提取相位(这些可以被认为是滤波器组方法的子集)
我的猜测是,在 (2,3) 中,相位会不时发生令人惊讶的事情,而在 (1) 中,我们可能无法拒绝尽可能多的噪音。我不清楚这个问题甚至有一个完美的解决方案(时频分辨率中的不确定性原理?)。
无论如何,如果有人以前解决过这个问题,而且......更好的是,如果有人知道任何听起来直接适用于此的论文,我将不胜感激。