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;