FFT of discrete time domain data

Question

shravankumar 2011-4-19

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/5823-fft-of-discrete-time-domain-data

评论： Imantha Ahangama 2018-10-21

Dear all,

I have vibration data from a setup in time domain. I need to convert into frequency domain,for which i have coded as below

n=length(data);
ts=t(2)-t(1);
p=abs(fft(data));
p=p(1:n/2).^2;
freq=[1:n/2]/ts;
figure(2)
plot(freq,p)

But the frequency components obtained in the plot, are not matching with actual ones,which i think the error in frequency axis of fft plot.Is any normalisation of axis needed . how to rectify this error .

Regards.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Arnaud Miege 2011-4-19

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/5823-fft-of-discrete-time-domain-data#answer_8122

在 MATLAB Online 中打开

I would suggest following the example in the fft documentation:

n = length(data); 
Ts = t(2)-t(1); % sample time
Fs = 1/Ts; % sampling frequency
NFFT = 2^nextpow2(n); % Next power of 2 from length of data
Y = fft(data,NFFT)/n;
f = Fs/2*linspace(0,1,NFFT/2+1);
% Plot single-sided amplitude spectrum.
plot(f,2*abs(Y(1:NFFT/2+1))) 
title('Single-Sided Amplitude Spectrum of y(t)')
xlabel('Frequency (Hz)')
ylabel('|Y(f)|')

HTH,

Arnaud

3 个评论
显示 1更早的评论隐藏 1更早的评论

Ella 2015-12-1

Hi Shravankumar do you get the same amplitude of vibration as your time domain ? I tried using the same code to get FFT of my signal but the amplitude drops. How can I solve this ? your suggestion is appreciated

Imantha Ahangama 2018-10-21