python - 如何找到 DWT 信号变换的频段？

Question

我是信号处理的新手，我想找出 Daubechies 小波“db4”变换输出的每个电平的频率范围。转换是使用 PyWavelets 完成的。我正在使用 python，下面的代码输出 5 个细节级别和 1 个近似值，但是我不确定每个级别描述的频率范围。

data = pywt.wavedec(X_train_ch1_raw[0], 'db4', level=5)

Answer 1

你的问题比看起来更棘手。

简短的回答是：使用 pywt 的scale2freq内置函数以给定的比例返回与给定小波相关的频率。例如，下面的代码返回 Daubechies 4 小波的频率，比例为 5 (0.14285714285714285)：

import pywt pywt.scale2frequency('db4',5)

您可以通过计算 db4 小波的中心频率 (0.7142857142857143) 然后除以比例 (5) 来获得相同的结果

import pywt pywt.central_frequency('db4')/5

请注意，这不是信号的实际中心频率！这个量被称为伪频率，因为它独立于被分析的信号。

为了恢复信号的中心频率，需要将伪频率除以信号的采样率：

import pywt pywt.scale2frequency('db4',5)/dt

dt你的采样率在哪里。

我希望这有帮助！

PS：我建议绘制重建信号的频谱，以说服自己中心频率与上述解析公式输出的值相匹配。