% Load female straight walking dataset
data = load(['A', filesep, '1.1.mat']);

% Find the index of the right knee joint in the cell array of trajectory
% labels where the individual letters in the label 'rKJZ_R' have the
% following meaning:
% r: revolute joint
% KJ: knee joint
% Z: rotation about local z-axis
% R: right side
indexKJ_R = find(strcmp(data.motion.trajectoryLabels, 'rKJZ_R'));

% Plot the right knee joint angle
timeKJ_R = 0:(1 / data.motion.frameRate):((data.motion.frames - 1) / data.motion.frameRate);
plot(timeKJ_R, data.motion.trajectory.q(indexKJ_R, :));
xlabel('Time in s');
ylabel('Right knee angle in rad');

% Load female straight running dataset
data = load(['A', filesep, '2.1.mat']);

% Find the index of the left ankle joint in the cell array of trajectory
% labels where the individual letters in the label 'rAJZ_L' have the
% following meaning:
% r: revolute joint
% AJ: ankle joint
% Z: rotation about local z-axis
% L: left side
indexAJ_L = find(strcmp(data.motion.trajectoryLabels, 'rAJZ_L'));

% Find the index of the left soleus muscle in the cell array of muscle
% labels where the individual letters in the label 'SOL_L' have the
% following meaning:
% SOL: soleus muscle
% L: left side
indexSOL_L = find(strcmp(data.muscle.muscleLabels, 'SOL_L'));

% Compute start and end indices for the 10th step
stepStart = data.events.eventStart_L(10);
stepEnd = data.events.eventStart_L(11);
startIndexAJ_L = round(stepStart * data.motion.frameRate) + 1;
endIndexAJ_L = round(stepEnd * data.motion.frameRate) + 1;
startIndexSOL_L = round(stepStart * data.muscle.frameRate) + 1;
endIndexSOL_L = round(stepEnd * data.muscle.frameRate) + 1;
startIndexGRFY_L = round(stepStart * data.force.frameRate) + 1;
endIndexGRFY_L = round(stepEnd * data.force.frameRate) + 1;

% Plot the left ankle joint acceleration, filtered activity of the left
% soleus muscle and left vertical ground reaction force of the 10th step
subplot(1, 3, 1);
timeAJ_L = stepStart:(1 / data.motion.frameRate):stepEnd;
plot(timeAJ_L, data.motion.trajectory.ddqddt(indexAJ_L, startIndexAJ_L:endIndexAJ_L));
xlabel('Time in s');
ylabel('Left ankle acceleration in rad/s^2');
subplot(1, 3, 2);
timeSOL_L = stepStart:(1 / data.muscle.frameRate):stepEnd;
plot(timeSOL_L, 1e3 * data.muscle.activities.filtered(indexSOL_L, startIndexSOL_L:endIndexSOL_L));
xlabel('Time in s');
ylabel('Activity of the left soleus muscle in mV');
subplot(1, 3, 3);
timeGRFY_L = stepStart:(1 / data.force.frameRate):stepEnd;
plot(timeGRFY_L, data.force.grfY_L(startIndexGRFY_L:endIndexGRFY_L));
xlabel('Time in s');
ylabel('Left ground reaction force in N');