Julia 真实 FFT 中的频率答案

【问题标题】：frequencies in Julia's real FFTJulia 真实 FFT 中的频率
【发布时间】：2018-10-25 22:28:07
【问题描述】：

我正在使用 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 的默认值。我阅读了文档，但仅通过阅读就不清楚频率是如何组织的。

有什么地方可以告诉我们频率的顺序吗？

【问题讨论】：

标签： julia fft fftw

【解决方案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.
"""

【讨论】：