我搜索用离散时间信号绘制时频信号(采样步长=0.001秒)。我使用Python和Scipy.signal库。我使用函数cwt(data,wavelet,widths),它返回一个矩阵,对复morlet小波(或gabor小波)进行连续小波变换。不幸的是,这方面的文献并不多。我发现最好的是:
-this对于Matlab(我试图找到相同的标度时间结果),但我自然无法访问相同的函数,
-以及this解释了什么是连续小波变换,没有详细的小波参数。
第一步:获取刻度转换信号。毫无疑问,我将数组“宽度”与可能的不同比例的数组直接关联起来。因为,如果参数宽度不是比例的话,我不知道它是什么。也许,你会告诉我“这是你当前小波的宽度”!但是,即使现在,我也不确定链接的宽度和比例如何…在Scipy的Morlet文档中,链接可能是:“s:缩放因子,从-s*2*pi加窗到+s*2*pi”,所以,我认为宽度=4*pi*scale(宽度=窗口的宽度)。但当我画小波时,尺度越大,小波的视觉宽度就越小。。。
我的第二个问题是找到并画出与频率相等的。在文献中,我发现这个公式:Fa=Fc/(s*delta),其中Fa是最终频率,Fc是以赫兹为单位的小波中心频率,s是尺度,delta是采样周期。所以,对于标度(如果我找到与宽度相关的链接)和delta(=0.001秒)来说是可以的,但是对于小波的中心频率来说更为复杂。在scipy文档中,我发现:“这个小波的基频[morlet小波]单位为赫兹,由f=2*s*w*r/M给出