java - 希尔伯特-黄变换: Instantaneous frequency

Question

我正在用 Java 实现“经验模式分解”。下一步是希尔伯特-黄变换，我需要找到如何获得“瞬时频率”。有人知道如何获得吗？

使用 EMD 方法，任何复杂的数据集都可以分解为有限且通常数量较少的分量，即固有模态函数 (IMF) 的集合。下一步是希尔伯特-黄变换。由此我计算信号的相位，在必须计算瞬时频率之后，我在 Matlab 代码中进行搜索，如下链接所示:http://read.pudn.com/downloads100/sourcecode/math/408870/emd/instfreq/archive/ifreq.m__.htm

但是我不明白。

谢谢

Answer 1

实例频率可以通过Python计算

def hilb(s, unwrap=False):
    from scipy.signal import hilbert
    H = hilbert(s)
    amp = np.abs(H)
    phase = np.arctan2(H.imag, H.real)
    if unwrap: phase = np.unwrap(phase)
    return amp, phase

inst_amp, phase = hilb(imf, unwrap=True)
inst_freq = np.diff(phase)