частоты в реальном БПФ Джулии - PullRequest
Новая
       10

частоты в реальном БПФ Джулии

0 голосов
/ 16 мая 2018

Я использую реализацию БПФ Джулии для выполнения двухмерного реального БПФ на нескольких массивах, но я не могу быть уверен в порядке частот на выходе.Рассмотрим 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).

У Джулии нет функции rfftfreq, как у Numpy, и тот факт, что выходные данные Джулии отличаются от выходных данных Numpy fft.rfftn, заставляет меня не хотеть использовать настройки Numpy по умолчанию.Я прочитал документацию, но не ясно, как организованы частоты, просто прочитав это.

Есть ли где-нибудь, что говорит нам порядок частот?

1 Ответ

0 голосов
/ 16 мая 2018

Я не уверен, что именно вы ищете, но если вы используете 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.
"""
...