FFT of discrete time domain data

3 次查看（过去 30 天）
shravankumar
shravankumar 2011-4-19
评论： 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.

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

采纳的回答

Arnaud Miege
Arnaud Miege 2011-4-19
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更早的评论
Ella
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
Imantha Ahangama 2018-10-21
@Ella, Try
data = data - mean(data);
before typing the above code.

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Fourier Analysis and Filtering 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by