scipy.fft.hfft¶
- scipy.fft.hfft(x, n=None, axis=- 1, norm=None, overwrite_x=False, workers=None, *, plan=None)[源代码]¶
计算具有厄米对称性的信号的FFT,即实谱。
- 参数
- xarray_like
输入数组。
- n整型,可选
输出的转换轴的长度。为 n 输出点,
n//2 + 1
输入点是必需的。如果输入的长度超过此长度,则会将其裁剪。如果它比这个短,则用零填充。如果 n 不是给定的,它被认为是2*(m-1)
,在哪里m
是输入沿由指定的轴的长度 axis 。- axis整型,可选
要在其上计算FFT的轴。如果未指定,则使用最后一个轴。
- norm{“向后”,“正向”,“向前”},选填
标准化模式(请参见
fft
)。默认值为“向后”。- overwrite_x布尔值,可选
如果为True,则 x 可以销毁;默认值为false。看见
fft
了解更多详细信息。- workers整型,可选
用于并行计算的最大工作进程数。如果为负值,则值从
os.cpu_count()
。看见fft
了解更多详细信息。- plan对象,可选
此参数保留用于传递由下游FFT供应商提供的预计算计划。它目前没有在本科学计划中使用。
1.5.0 新版功能.
- 退货
- outndarray
沿由指示的轴转换的截断或填零的输入 axis ,或者最后一个,如果 axis 未指定。变换后的轴的长度为 n ,或者,如果 n 不是给出的,
2*m - 2
,在哪里m
是输入的转换轴的长度。为了获得奇数个输出点, n 必须指定,例如,指定为2*m - 1
在典型情况下,
- 加薪
- IndexError
如果 axis 大于的最后一个轴 a 。
注意事项
hfft
/ihfft`是类似于 `rfft /‘irfft’,但是对于相反的情况:这里的信号在时域中具有厄米对称性,并且在频域中是实数。所以,在这里,它是hfft
,如果结果是奇数,则必须提供结果的长度。 * even:ihfft(hfft(a, 2* LEN(A)-2)==a
,在舍入误差内, * odd:ihfft(hfft(a, 2* LEN(A)-1)==a
,在舍入误差范围内。示例
>>> from scipy.fft import fft, hfft >>> 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.])