R1=0;              % Starting Row
C1=3;              % Starting Column
%% uvoz
filename='F0000CH1.csv';
dataCsv = csvread(filename,R1,C1);
dataCsv(:,end)=[]; % remove last column

%% vreme - signal
vreme = dataCsv(:,1);
Ts=vreme(end)-vreme(end-1);
fs=1/Ts;
signal = dataCsv(:,2);

%% odsecanje signala
vreme=vreme(1:2000);
signal=signal(1:2000);

%FFT sirovog signala
[Ns,rezolucija,ucesatnostXosa,amplituda,faza] = fft_dz( signal,fs);

% prozorska funkcija
signal_prozor=signal.*hamming(Ns);
[Ns,rezolucija,ucesatnostXosa,amplituda_prozor,faza_prozor] = fft_dz( signal_prozor,fs);
%% plots
%Plot in time
figure(1),plot(vreme,signal)
title(['Plot of ',filename])
xlabel('vreme [s]'),ylabel('amplituda')

% FFT
figure(2),plot(ucesatnostXosa,amplituda,'b.-',...
               ucesatnostXosa,amplituda_prozor,'r.-');
xlabel('ucestanost [Hz]'),ylabel('amplituda')
title(['Ceo spektar - ','rezolucijaXosa=[',num2str(0),':',num2str(rezolucija),':',num2str(fs/2),']'])
legend('bez prozora','sa prozorom')
figure(3)
plot(ucesatnostXosa,amplituda,'b.-',...
               ucesatnostXosa,amplituda_prozor,'r.-');
xlim([0, 1e3])
xlabel('ucestanost [Hz]'),ylabel('amplituda')
title(['Uvelicannje po x-osi - ','rezolucijaXosa=[',num2str(0),':',num2str(rezolucija),':',num2str(fs/2),']'])
legend('bez prozora','sa prozorom')