scipy.cwt函数中的Morlet小波

Question

我使用的是scipy的连续小波变换。

在文档中，signal.ricker函数用于小波：

https://docs.scipy.org/doc/scipy-0.16.0/reference/generated/scipy.signal.cwt.html

代码使用signal.ricker函数，该函数使用墨西哥帽小波：

from scipy import signal
import matplotlib.pyplot as plt
import numpy as np
import pywt


sig  = data
widths = np.arange(1, 31)

cw = signal.cwt(sig, signal.ricker, widths)

现在我想用morlet小波代替signal.ricker和墨西哥帽小波。

到目前为止，我已经尝试了这些morlet函数，但没有任何成功：

morl1 = signal.morlet(125, w=5.0, s=1.0, complete=True)

和：

wavelet = pywt.ContinuousWavelet('morl')
zA, zD = wavelet.wavefun(level=8)
morl2 = zA

有没有什么方法可以将morlet小波用于连续小波变换？

cw = signal.cwt(sig, morlet, widths)

Answer 1

原始signal.morlet不适合与signal.cwt-function一起使用。

您可以改用scipy.signal.morlet2，这将创建一个专门为signal.cwt创建的复杂morlet小波。