我搜索以显示具有原始离散时间信号的时频信号(采样步长 = 0.001 秒)。我使用 Python 和库Scipy.signal。我使用该函数使用复morlet 小波(或 gabor 小波)cwt(data, wavelet, widths)
进行连续小波变换。
第一步:获取比例转换信号。有疑问,我直接将数组“宽度”与可能不同比例的数组联系起来。因为,如果不是比例,我不明白什么是参数宽度。也许,你会告诉我“这是你当前小波的宽度”!但是,即使是这样,我也不知道宽度与比例是如何联系起来的……
我的第二个问题是找到并显示等价的频率。在文献中,我找到了这个公式:Fa = Fc / (s*delta),其中 Fa 是最终频率,Fc 是小波的中心频率,以 Hz 为单位,s 是尺度,delta 是采样周期。因此,对于比例(如果我找到与宽度的链接)和增量(= 0.001 秒)来说没问题,但是小波的中心频率更复杂。在 scipy 文档中,我发现:“这个小波 [morlet 小波] 的基频 (Hz) 由 f = 2*s*w*r / M 给出,其中 r 是采样率 [s 在这里是缩放因子,窗口化从 -s*2*pi 到 +s*2*pi。默认为 1;w 宽度;和 M 小波的长度]。我认为是中心频率,是吗?也许我的第一个问题的解决方案也在这里(比例因子和比例..?)..?
谢谢