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()