基於頻域頻譜的FFT

Question

我有一個功能，我想根據頻域繪制FFT。

t=0:0.01:60;

y = sin(2*pi*0.33*t)+sin(2*pi*0.23*t)+sin(2*pi*0.1*t);

如您所見，峰值應該在0.1、0.2、0.33頻率上，但是當我嘗試繪制它時，峰值不在這些點上。

我嘗試了FFT，DFT，但沒有一個能如我所願。 另外，我希望它們的FFT峰值在這些頻率下等於1。

有什么想法嗎？

Answer 1

如果Fs是采樣率（聽到100 Hz），則頻率倉的長度為（t）/ Fs Hz。 由於您要使垃圾箱間隔100 Hz，因此需要更長的信號。 嘗試這個

Fs=100; %sample rate
t=0:1/Fs:(100-1/Fs);  %10,000 samples long
F=[0:length(t)-1]/length(t)*Fs; %where the bins are
y = sin(2*pi*0.33*t)+sin(2*pi*0.23*t)+sin(2*pi*0.1*t); 
Y=abs(fft(y))/(length(y)/2);  %normalize by the signal length
figure(1);
plot(F,Y)
xlim([0 1]); %zoom in to the part of interest