XX=accread2('DATA13.DAT');
%XX=accread2('DATA1.DAT',54e3,58e3);

xmin=1785;xmax=2288; g2x=xmax-xmin; e_X=2*9.8/g2x;
ymin=1789;ymax=2296; g2y=ymax-ymin; e_Y=2*9.8/g2y;


freq=1041
d=length(XX)
%d=1000000;
i=1:1:d;

time=i/freq;
z=XX(:,2); z=e_Y*(z(1:d)-z(1));
x=XX(:,1);x=e_X*(x(1:d)-x(1));
xz=sqrt(x.^2+z.^2);
figure; plot(time,(x),time,(z), time,xz)

%figure; plot(i,(x),i,(z), i,xz)

dx=1000;

Count=[];Count_X=[];Count_Z=[];time_count=[];
for i=1:dx:floor(d/dx)*dx
   
    Count_Z=[Count_Z ; sum(abs(z(i:(i+dx))))];
    Count_X=[Count_X ; sum(abs(x(i:(i+dx))))];
    Count=[Count ; sum(abs(xz(i:(i+dx))))];
    
    time_count=[time_count;(time(i)+time(floor(dx/2)))];
end

Total=sum(xz);Total_X=sum(abs(x));Total_Z=sum(abs(z));
Norm_Count=Count/max(Count)*Total;Norm_Count_X=Count_X/max(Count_X)*Total_X; Norm_Count_Z=Count_Z/max(Count_Z)*Total_Z;
figure;
bar(time_count,Norm_Count)
hold on
plot(time, cumsum(xz),'Linewidth',3,'Color','r')
plot(time, cumsum(abs(z)),'Linewidth',2,'Color','b');
plot(time, cumsum(abs(x)),'Linewidth',2,'Color','g');
title('ABSOLTNA VREDNOST X IN Z')

hold off;
figure; 
bar(time_count,Norm_Count_X)
hold on
plot(time, cumsum(abs(x)),'Linewidth',3,'Color','r')
title('X OS')

hold off;
figure; 
bar(time_count,Norm_Count_Z)
hold on
plot(time, cumsum(abs(z)),'Linewidth',3,'Color','r');
title('Z OS')