Получение и построение волновых векторов из fft

Question

У меня есть набор данных, который описывает мгновенную скорость среды в диапазоне от x=0 до x=0.4 (с 32 узлами пространственное приращение составляет 0,0125).

Я хотел бы применить FFT к этому, это я в следующем фрагменте кода, где data содержит скоростей:

data = [-0.0005 
    -0.0004 
    -0.0003 
    -0.0002 
    -0.0001 
    -0.0000 
    0.0001 
    0.0001 
    0.0001 
    0.0002 
    0.0002 
    0.0002 
    0.0002 
    0.0002 
    0.0002 
    0.0002 
    0.0002 
    0.0002 
    0.0003 
    0.0004 
    0.0005 
    0.0006 
    0.0007 
    0.0009 
    0.0010 
    0.0011 
    0.0011 
    0.0012 
    0.0011 
    0.0011 
    0.0011 
    0.0010 ]; 


N = 100; 
X=abs(fft(data,N)); 
X = fftshift(X); 
F = [-N/2:N/2-1]/N; 

figure(1) 
plot(F,X) 

Однако, я не уверен, о том, какие значения волнового вектора x -аксис представляет собой график. В принципе, я хочу волновые векторы, но есть ли способ получить это от БПФ?

Answer 1

Ваш F, работающий от -0,5 (включительно) до +0,5 (исключая), представляет частоту относительно базовой частоты дискретизации (частота Найквиста). Поскольку ваш исходный сигнал X отбирается при равномерном расстоянии 0,0125 м на образец или частота выборки из 80 отсчетов на метр, F/0.0125 представляет вашу реальную пространственную частоту в единицах циклов на метр. Имеет смысл?