
PD03 = 'PD03';

%Start axis is min value of the BH file (BH is put on first)
%startx = min(PD03_BHAccXYZTime);
%End axis is the max value of the BH file
%endx = max(PD03_BHAccXYZTime);

%use the bioharness timing to set the size of figure
settime = PD03_BHAccXYZTime;

%axis(xlim,[startx endx]);

figure('Color','w')


subplot(5,1,1)
plot(PD03_BHAccXYZTime,PD03_BHAccXYZ_Xonly,'r')
hold on
plot(PD03_BHAccXYZTime,PD03_BHAccXYZ_Yonly,'b')
hold on
plot(PD03_BHAccXYZTime,PD03_BHAccXYZ_Zonly,'g')
title('BioHarness Acceleration (chest)')
xlim(settime([1 end]))
datetick('x','HH:MM:SSPM','keepticks')

subplot(5,1,2)
plot(PD03_BodyDynTime,PD03_BodyDyn_accXonly,'r')
hold on
plot(PD03_BodyDynTime,PD03_BodyDyn_accYonly,'b')
hold on
plot(PD03_BodyDynTime,PD03_BodyDyn_accZonly,'g')
title('BodyDyn Acceleration (ankle)')
xlim(settime([1 end]))
datetick('x','HH:MM:SSPM','keepticks')

subplot(5,1,3)
plot(PD03_BodyDynTime,PD03_BodyDyn_gyroXonly,'r')
hold on
plot(PD03_BodyDynTime,PD03_BodyDyn_gyroYonly,'b')
hold on
plot(PD03_BodyDynTime,PD03_BodyDyn_gyroZonly,'g')
title('BodyDyn Gyroscope (ankle)')
xlim(settime([1 end]))
datetick('x','HH:MM:SSPM','keepticks')


subplot(5,1,4)
plot(PD03_QSensorDataTime,PD03_QSensorData_Xonly,'r')
hold on
plot(PD03_QSensorDataTime,PD03_QSensorData_Yonly,'b')
hold on
plot(PD03_QSensorDataTime,PD03_QSensorData_Zonly,'g')
title('QSensor Acceleration (wrist)')
xlim(settime([1 end]))
datetick('x','HH:MM:SSPM','keepticks')


subplot(5,1,5)
plot(PD03_BHPeakAccTime,PD03_BHPeakAcc)
title('BioHarness Peak Acceleration (chest)')
xlim(settime([1 end]))
datetick('x','HH:MM:SSPM')
set(gca,'XTick',[settime(1):.025:settime(end)])
xlim(settime([1 end]))
datetick('x','HH:MM:SSPM')

% may have to run this loop outside of the script - for unknown reasons

for i = 1:5;
subplot(5,1,i)
xlim(settime([1 end]))
set(gca,'XTick',[settime(1):.025:settime(end)])

datetick('x','HH:MM:SSPM')
grid on
grid minor
end

hold on
