scipy.signal.windows.

kaiser_bessel_derived

scipy.signal.windows.kaiser_bessel_derived(M, beta, *, sym=True, xp=None, device=None)

返回 Kaiser-Bessel 派生窗口。

参数:

Mint

输出窗口中的点数。如果为零，则返回一个空数组。当为负数时，会抛出一个异常。请注意，此窗口仅针对偶数个点定义。

betafloat

Kaiser 窗口形状参数。

symbool, 可选

此参数的存在只是为了符合其他窗口函数提供的接口，并且可以被 get_window 调用。当为 True（默认）时，生成一个对称窗口，用于滤波器设计。

xparray_namespace, 可选

可选的数组命名空间。应与数组 API 标准兼容，或由 array-api-compat 支持。默认值：numpy

device: any

输出的可选设备规范。应与 xp 中支持的设备规范之一匹配。

返回值:

wndarray

窗口，已标准化以满足 Princen-Bradley 条件。

参见

kaiser

注释

它被设计为适用于修改后的离散余弦变换 (MDCT)，主要用于音频信号处理和音频编码。

在 1.9.0 版本中添加。

参考文献

[1]

Bosi, Marina, and Richard E. Goldberg. Introduction to Digital Audio Coding and Standards. Dordrecht: Kluwer, 2003.

[2]

Wikipedia，“Kaiser window”，https://en.wikipedia.org/wiki/Kaiser_window

示例

在你的浏览器中尝试一下！

根据维基百科参考资料 [2] 绘制 Kaiser-Bessel 派生窗口

>>> import numpy as np
>>> from scipy import signal
>>> import matplotlib.pyplot as plt
>>> fig, ax = plt.subplots()
>>> N = 50
>>> for alpha in [0.64, 2.55, 7.64, 31.83]:
...     ax.plot(signal.windows.kaiser_bessel_derived(2*N, np.pi*alpha),
...             label=f"{alpha=}")
>>> ax.grid(True)
>>> ax.set_title("Kaiser-Bessel derived window")
>>> ax.set_ylabel("Amplitude")
>>> ax.set_xlabel("Sample")
>>> ax.set_xticks([0, N, 2*N-1])
>>> ax.set_xticklabels(["0", "N", "2N+1"])  
>>> ax.set_yticks([0.0, 0.2, 0.4, 0.6, 0.707, 0.8, 1.0])
>>> fig.legend(loc="center")
>>> fig.tight_layout()
>>> fig.show()

返回在新标签页中打开