Code: Select all

load('data')
load('time')

Fs = 2000;
L= numel(time);
Fn = Fs/2;
NFFT = 2^nextpow2(L);
FTdata = fft(data,NFFT)/L;
Fv = linspace(0, 1, NFFT/2+1)*Fn;
Iv = 1:numel(Fv);

figure
plot(Fv, abs(FTdata(Iv))*2)
grid
xlabel('Frequency (Hz)')
ylabel('Amplitude')