scipy.fft.
hfft#
- scipy.fft.hfft(x, n=None, axis=-1, norm=None, overwrite_x=False, workers=None, *, plan=None)[源代码]#
计算具有埃尔米特对称性(即,实频谱)的信号的 FFT。
- 参数:
- xarray_like
输入数组。
- nint, 可选
输出的变换轴的长度。对于 n 个输出点,需要
n//2 + 1
个输入点。如果输入长度大于此值,则会裁剪。如果短于此值,则会用零填充。如果未给出 n,则将其取为2*(m-1)
,其中m
是输入沿 axis 指定的轴的长度。- axisint, 可选
计算 FFT 的轴。如果未给出,则使用最后一个轴。
- norm{"backward", "ortho", "forward"}, 可选
归一化模式(请参阅
fft
)。默认值为“backward”。- overwrite_xbool, 可选
如果为 True,则可以销毁 x 的内容;默认值为 False。有关更多详细信息,请参阅
fft
。- workersint, 可选
用于并行计算的最大工作线程数。如果为负数,则该值从
os.cpu_count()
环绕。有关更多详细信息,请参阅fft
。- planobject, 可选
此参数保留用于传递下游 FFT 供应商提供的预计算计划。它目前在 SciPy 中未使用。
在版本 1.5.0 中添加。
- 返回:
- outndarray
沿 axis 指示的轴变换的截断或零填充输入,如果未指定 axis,则沿最后一个轴变换。变换轴的长度为 n,或者,如果未给出 n,则为
2*m - 2
,其中m
是输入的变换轴的长度。要获得奇数个输出点,必须指定 n,例如,在典型情况下为2*m - 1
,
- 引发:
- IndexError
如果 axis 大于 a 的最后一个轴。
说明
hfft
/ihfft
是一对类似于rfft
/irfft
的组合,但适用于相反的情况:此处信号在时域中具有埃尔米特对称性,在频域中是实数。因此,这里是hfft
,如果结果为奇数,则必须提供结果的长度。* 偶数:ihfft(hfft(a, 2*len(a) - 2) == a
,在舍入误差范围内,* 奇数:ihfft(hfft(a, 2*len(a) - 1) == a
,在舍入误差范围内。示例
>>> from scipy.fft import fft, hfft >>> import numpy as np >>> a = 2 * np.pi * np.arange(10) / 10 >>> signal = np.cos(a) + 3j * np.sin(3 * a) >>> fft(signal).round(10) array([ -0.+0.j, 5.+0.j, -0.+0.j, 15.-0.j, 0.+0.j, 0.+0.j, -0.+0.j, -15.-0.j, 0.+0.j, 5.+0.j]) >>> hfft(signal[:6]).round(10) # Input first half of signal array([ 0., 5., 0., 15., -0., 0., 0., -15., -0., 5.]) >>> hfft(signal, 10) # Input entire signal and truncate array([ 0., 5., 0., 15., -0., 0., 0., -15., -0., 5.])