format long

gnuplot_binary ("gnuplot")

data = load ('/test/a.dat');

time = data(:,1);

signal = data(:,2);

T = time(2) - time(1);

Fs = 1/T;

L = length(signal);

NFFT = 2^nextpow2(L);

Y = fft(signal,NFFT)/L;

f= Fs/2*linspace(0,1,NFFT/2+1);

A(:,1) = f;

A(:,2) = (2*abs(Y(1:NFFT/2+1)));

B = sortrows(A,2);

figure(1, "visible", "off");

plot(A(:,1),A(:,2))

axis ([0.005 1])

print -deps "/test/FA.eps"

for i = 1:60

freq(i) = B(length(B)-(i-1),1);

Amp(i) = B(length(B)-(i-1),2);

Period(i) =1/freq(i);

file_id = fopen(/test/FrequencyA.dat', 'at');

fprintf(file_id, ' \n \t %u Max Amplitude:\t %f \t Dominant frequency %u: \t %f \t Period: \t %f \t\n',(i), Amp(i),(i), freq(i), Period(i));

fclose(file_id);

end

exit;