% bactreai.m

clear; clc;
[t X Y] = textread('bacteria_position.txt');
x = X;
y = Y;

for i = 2:length(t)-1
    vx(i) = (x(i+1)-x(i-1));
    vy(i) = (y(i+1)-y(i-1));
     V(i) = sqrt(vx(i)^2+vy(i)^2);
end

for i = 2:length(V)-1
    ax = (vx(i+1)-vx(i-1));
    ay = (vy(i+1)-vy(i-1));
    A(i) = sqrt(ax^2+ay^2);
end

subplot(3,2,1)
  plot(x,y,'*-')
  xlabel('X POSITION (Pixel)');ylabel('Y POSITION (Pixel)')
  grid on
subplot(3,2,3)
  plot(t,x,'*-')
  xlabel('TIME (s)');ylabel('X POSITION (Pixel)')
  grid on
subplot(3,2,4)
  plot(t,y,'*-')
  xlabel('TIME (s)');ylabel('Y POSITION (Pixel)')
  grid on
subplot(3,2,5)
  plot(t(2:end-1),V(2:end),'*-')
  xlabel('TIME (s)');ylabel('SPEED (Pixel/s)')
  grid on
subplot(3,2,6)
  plot(t(3:end-1),A(2:end),'*-')
  xlabel('TIME (s)');ylabel('ACCELERATION (Pixel/s^2)')
  grid on