Python SciPy signal.morlet用法及代码示例

本文简要介绍 python 语言中 scipy.signal.morlet 的用法。

用法:

scipy.signal.morlet(M, w=5.0, s=1.0, complete=True)#

复 Morlet 小波。

参数 ：：

M： int

小波的长度。

w： 浮点数，可选

欧米茄0。默认值为 5

s： 浮点数，可选

缩放因子，窗口从 -s*2*pi 到 +s*2*pi 。默认值为 1。

complete： 布尔型，可选

是使用完整版还是标准版。

返回 ：：

morlet： (M,) ndarray

注意：

标准版：

pi**-0.25 * exp(1j*w*x) * exp(-0.5*(x**2))

这种常用的小波通常简称为 Morlet 小波。请注意，此简化版本可能会在 w 值较低时导致可接纳性问题。

完整版：

pi**-0.25 * (exp(1j*w*x) - exp(-0.5*(w**2))) * exp(-0.5*(x**2))

此版本有一个更正术语以提高可接受性。对于 w 大于 5，校正项可以忽略不计。

请注意，返回小波的能量没有根据以下公式进行归一化：s.

该小波的基频(以 Hz 为单位)由下式给出f = 2*s*w*r / M其中r是采样率。

注意：该函数是在 cwt 之前创建的，与其不兼容。

例子：

>>> from scipy import signal
>>> import matplotlib.pyplot as plt

>>> M = 100
>>> s = 4.0
>>> w = 2.0
>>> wavelet = signal.morlet(M, s, w)
>>> plt.plot(wavelet.real, label="real")
>>> plt.plot(wavelet.imag, label="imag")
>>> plt.legend()
>>> plt.show()