Question

我正在使用Julia的FFT实现在几个阵列上执行2D实数FFT，但我不能确定输出中频率的顺序。考虑MWE

N=64
U = rand(Float64, N, N);

FFTW.set_num_threads(2)
prfor = plan_rfft(U, (1,2), flags=FFTW.MEASURE);

size(prfor*U)

输出是一个大小为(33, 64)的数组。

Julia没有像Numpy那样的rfftfreq函数，而且Julia的输出与Numpy的fft.rfftn默认输出不同，这让我不想在这里使用Numpy的默认值。我阅读了文档但不清楚如何通过阅读它来组织频率。

有没有告诉我们频率顺序的地方？

Answer 1

我不确定你究竟在寻找什么，但是如果你使用DSP.jl，它的util.jl文件可能有你需要的东西：

https://github.com/JuliaDSP/DSP.jl/blob/master/src/util.jl

"""
    rfftfreq(n, fs=1)
Return discrete fourier transform sample frequencies for use with
`rfft`. The returned Frequencies object is an AbstractVector
containing the frequency bin centers at every sample point. `fs`
is the sample rate of the input signal.
"""

朱莉娅真实FFT中的频率

1 个答案: