Eşit ve tek FFT uzunluğu için frekans ekseni nasıl yapılır?

Negatif frekanstan pozitif frekansa (Hertz cinsinden) frekans ekseninin nasıl yapılacağı konusunda yardım alabilir miyim, bu bir FFT sonucundaki x ekseni olacaktır, ancak çift uzunluklu FFT veya tek uzunluklu FFT verildi. MATLAB'da bazı sıkıntılar yaşıyorum. (Örnekleme frekansını f_s bildiğinizi varsayalım).

Bir yaklaşım basitçe, kaydırılmamış DFT çıkışı için frekans vektörünü hesaplamaktır (yani fft(), bir an yapmadan doğrudan MATLAB işlevinden alacağınız şey fftshift()), ardından eksenin negatif tarafındaki konumlara karşılık gelen frekansları yeniden ifade etmektir. Misal:

% assume input signal "x", sampling frequency "fs"
% calculate FFT
X = fft(x,Nfft);
% calculate frequency spacing
df = fs / Nfft;
% calculate unshifted frequency vector
f = (0:(Nfft-1))*df;
% move all frequencies that are greater than fs/2 to the negative side of the axis
f(f >= fs/2) = f(f >= fs/2) - fs;
% now, X and f are aligned with one another; if you want frequencies in strictly
% increasing order, fftshift() them
X_normal_order = fftshift(X);
f_normal_order = fftshift(f);

Learnvst tarafından verilen cevap da işe yaramalıdır; bu, çift / tek DFT boyutları için herhangi bir özel kasa gerektirmeyen başka bir düşünme şeklidir.

Olumlu bir frekans spektrumu oldukça basit bir şekilde yapabilirsiniz ( fsörnekleme oranı ve NFFTfft kutularının sayısı nerede ). FFT algoritmasının Matlab uygulamasında, ilk eleman her zaman DC bileşenidir, bu nedenle dizi neden sıfırdan başlar. Bu NFFT'nin tek ve çift değerleri için geçerlidir.

%//Calculate frequency axis
df = fs/NFFT;
fAxis = 0:df:(fs-df);

Frekans spektrumunu sarmanız gerekiyorsa, tek numaralı NFFT'niz olup olmadığını dikkate almanız gerekir. Her zaman bir DC bileşeni olması gerekir . .

df = fs/NFFT;
fAxis = (0:df:(fs-df)) - (fs-mod(NFFT,2)*df)/2;

Pozitif frekans ekseninin hesaplamasının yukarıyla nasıl aynı olduğuna dikkat edin, ancak FFT-kaydırılan terim, çift veya tek FFT uzunluklarını karşılamak için değişir.

Bu kod parçacıkları burada bulunan SO'da (ilginç bulabileceğiniz) uzun bir cevaptan alınmıştır: /programming/9694297/matlab-fft-xaxis-limits-messing-up-and-fftshift/ 9699983 # 9699983