% Wczytanie danych z pliku do macierzy A
A = load ("dane.dat");

% Ile pozycji jest w danych
L = length(A);

% Wykorzystanie wbudowanej funkcji fft. 
% Szybka transformata fouriera wymaga, 
% żeby liczba danych była potęgą dwójki
% dlatego wykorzystujemy funkcje 
% 'nextpow2' (patrz wyjasnienia w instrukcji Matlaba/Octave'a).
FFT = fft(A(:,7), 2^nextpow2(L)) / L ;

plot (FFT)
% Czestotliwosc i amplituda
freq = 1/2 * linspace(0,1,2^nextpow2(L)/2+1)';
ampl = 2 * abs( FFT( 1:2^nextpow2(L)/2+1));

% Mielismy dane godzinne a wyniki chcemy w cpd
freq = freq * 24;

% wyrzućmy pierwsze 40 wartości
freq = freq(40:end);
ampl = ampl(40:end);

% findpeaks(ampl)
% narysuj wyniki
plot(freq,ampl)
pause
