scipy.signal.ShortTimeFFT.

from_dual#

classmethod ShortTimeFFT.from_dual(dual_win, hop, fs, *, fft_mode='onesided', mfft=None, scale_to=None, phase_shift=0)[源代码]#

仅提供双窗口，实例化一个 ShortTimeFFT。

如果 STFT 可逆，则可以从给定的双窗口 dual_win 计算窗口 win。所有其他参数与 ShortTimeFFT 的初始化器中的含义相同。

如 短时傅里叶变换 部分的 SciPy 用户指南 中所述，可逆 STFT 可以解释为时移和频调双窗口的级数展开。例如，级数系数 S[q,p] 属于该项，该项将 dual_win 平移 p * delta_t，并将其乘以 exp( 2 * j * pi * t * q * delta_f）。

另请参阅

from_window

通过包装 get_window 创建实例。

ShortTimeFFT

使用标准初始化器创建实例。

示例

以下示例讨论将信号分解为时移和频移高斯函数。将使用由 51 个样本组成的标准差为 1 的高斯函数

>>> import numpy as np
>>> import matplotlib.pyplot as plt
>>> from scipy.signal import ShortTimeFFT
>>> from scipy.signal.windows import gaussian
...
>>> T, N = 0.1, 51
>>> d_win = gaussian(N, std=1/T, sym=True)  # symmetric Gaussian window
>>> t = T * (np.arange(N) - N//2)
...
>>> fg1, ax1 = plt.subplots()
>>> ax1.set_title(r"Dual Window: Gaussian with $\sigma_t=1$")
>>> ax1.set(xlabel=f"Time $t$ in seconds ({N} samples, $T={T}$ s)",
...        xlim=(t[0], t[-1]), ylim=(0, 1.1*max(d_win)))
>>> ax1.plot(t, d_win, 'C0-')

以下图显示了重叠为 41、11 和 2 个样本时 hop 间隔如何影响窗口 win 的形状

>>> fig2, axx = plt.subplots(3, 1, sharex='all')
...
>>> axx[0].set_title(r"Windows for hop$\in\{10, 40, 49\}$")
>>> for c_, h_ in enumerate([10, 40, 49]):
...     SFT = ShortTimeFFT.from_dual(d_win, h_, 1/T)
...     axx[c_].plot(t + h_ * T, SFT.win, 'k--', alpha=.3, label=None)
...     axx[c_].plot(t - h_ * T, SFT.win, 'k:', alpha=.3, label=None)
...     axx[c_].plot(t, SFT.win, f'C{c_+1}',
...                     label=r"$\Delta t=%0.1f\,$s" % SFT.delta_t)
...     axx[c_].set_ylim(0, 1.1*max(SFT.win))
...     axx[c_].legend(loc='center')
>>> axx[-1].set(xlabel=f"Time $t$ in seconds ({N} samples, $T={T}$ s)",
...             xlim=(t[0], t[-1]))
>>> plt.show()
../../_images/scipy-signal-ShortTimeFFT-from_dual-1_00_00.png
../../_images/scipy-signal-ShortTimeFFT-from_dual-1_00_01.png

除了以 t = 0 为中心的窗口 win 之外，还描绘了前一个 (t = -delta_t) 和后续窗口 (t = delta_t)。可以看出，对于较小的 hop 间隔，窗口是紧凑且平滑的，在 STFT 中具有良好的时频集中度。对于 4.9 秒的较大 hop 间隔，窗口在 t = 0 附近的值很小，这些值没有被相邻窗口的重叠所覆盖，这可能导致数值不准确。此外，窗口开头和结尾的尖峰形状表明带宽更高，导致 STFT 的时频分辨率较差。因此，hop 间隔的选择将是时频分辨率和较小的 hop 大小所需的内存需求之间的折衷方案。

本页内容