%
%
g=9.81;
load  EXP17_bounce1.xyz
t=EXP17_bounce1(:,1);
N=length(t);
tN=t(N);
i=1:N;
x=EXP17_bounce1(:,2);
y=EXP17_bounce1(:,3);
z=EXP17_bounce1(:,4);
cal=sqrt(x.*x+y.*y+z.*z)/g;
%
%   plot each channel
%
figure(1)
subplot(4,1,1);plot(t,x);
axis([0,tN,-20.,00.])
ylabel('acceleration m/s**2')
title('CH1')
subplot(4,1,2);plot(t,y);
axis([0,tN,-12.,12.])
ylabel('acceleration m/s**2')
title('CH2')
subplot(4,1,3);plot(t,z);
axis([0,tN,-12.,12.])
ylabel('acceleration m/s**2')
title('CH3')
subplot(4,1,4);plot(t,cal);
axis([0,tN,0.5,1.5])
ylabel('|g|/g0')
xlabel('time (s) ')
%
figure(2)
subplot(2,1,1);plot(t,cal);
title('Oscillator, EXP 1')
axis([0,tN,0.0,2.0])
ylabel('|g|/g0')
xlabel('time (s) ')
