How can I get the FFT (Amplitude - Frequency) from Time (s) - Acceleartion(m/s^2) data?

3 次查看（过去 30 天）

Bob 2016-3-21

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/274678-how-can-i-get-the-fft-amplitude-frequency-from-time-s-acceleartion-m-s-2-data

评论： Star Strider 2016-3-21

clc;
clear all;
close all;
Time = xlsread('Data','A2:A513'); % Time (s)
Acceleration = xlsread('Data','B2:B513'); % Acceleration (m/s^2)
n = length(Acceleration); 
Ts = Time(2)-Time(1);            % Sample Time
Fs = 1/Ts;                       % Sampling Frequency
NFFT = 2^nextpow2(n);            % Next power of 2 from length of data
Y = fft(Acceleration,NFFT)/n;
f = Fs/2*linspace(0,1,NFFT/2+1);
plot(f,Y) 
xlabel('Frequency (Hz)')
ylabel('Amplitude (m)')

I am trying to get the FFT from Time Acceleration Data.

Error using plot
Vectors must be the same length.

P.S I have attached the Data.xlsx

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

请先登录，再进行评论。

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

采纳的回答

Star Strider 2016-3-21

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/274678-how-can-i-get-the-fft-amplitude-frequency-from-time-s-acceleartion-m-s-2-data#answer_214436

在 MATLAB Online 中打开

In your code, ‘f’ is half the length of ‘y’ so you need to create an ‘index vector’.

Try this:

... CODE ...
Y = fft(Acceleration,NFFT)/n;
f = Fs/2*linspace(0,1,NFFT/2+1);
Iv = 1:length(f);                   % Index Vector
plot(f,abs(Y(Iv))) 
xlabel('Frequency (Hz)')
ylabel('Amplitude (m)')