We don't know anything about your files, but I don't think you want to combine your j and za for loops the way you did. I think you want to keep these as two separte loops.
Also note that in your calculation of sumMR you are using t0HR instead of t0MR.
clear all
close all
startpath ='1EXAMENOPDRACHT\Assignment_V2\Data';
Datapath = [startpath '\Data'];
% rechtersubjecten
FirstSwingExcursionAngle_HR = zeros(5,5)
FirstSwingExcursionAngle_MR = zeros(5,5)
sumHR = 0;
sumMR = 0;
avgHR = zeros(1,5);
avgMR = zeros(1,5);
for i = 1:5
subject = [3 5 7 9 14];
dataFolder = [startpath '\TD' num2str(subject(i))];
dataFile_HR = fullfile(dataFolder, '\IK\HR*.mot'); %naam van de map een bestandsnaam geven
pathwayIK_HR = dir(dataFile_HR);
dataFile_MR = fullfile(dataFolder, '\IK\MR*.mot');
pathwayIK_MR = dir(dataFile_MR)
for k = 1:length(pathwayIK_HR)
IKHRdata = importdata(fullfile(pathwayIK_HR(k).folder,pathwayIK_HR(k).name));
DataHRlinks = IKHRdata.data(:,strcmp('knee_angle_l',IKHRdata.colheaders));
TimeHRlinks = IKHRdata.data(:,strcmp('time',IKHRdata.colheaders));
dataFile_HR = fullfile(dataFolder, '\Event\HR*.mat'); %naam van de map een bestandsnaam geven
pathwayEvent_HR = dir(dataFile_HR);
EventHRdata = importdata(fullfile(pathwayEvent_HR(k).folder,pathwayEvent_HR(k).name));
t0knielinks_HR = EventHRdata.t0;
for j = 1:length(TimeHRlinks) % vervangen door functie
if t0knielinks_HR == TimeHRlinks(j)
t0HR = j
break
end
end
FirstSwingExcursionAngle_HR(k,i)=DataHRlinks(t0HR) % plot
sumHR = sumHR + DataHRlinks(t0HR);
end
avgHR(i)=sumHR/length(pathwayIK_HR) % staaf diagram van het gemiddelde
sumHR = 0; % i gebruiken voor de juiste subplots
for k = 1:length(pathwayIK_MR)
IKMRdata = importdata(fullfile(pathwayIK_MR(k).folder,pathwayIK_MR(k).name));
DataMRlinks = IKMRdata.data(:,strcmp('knee_angle_l',IKMRdata.colheaders));
TimeMRlinks = IKMRdata.data(:,strcmp('time',IKMRdata.colheaders));
dataFile = fullfile(dataFolder, '\Event\MR*.mat'); %naam van de map een bestandsnaam geven
pathwayEvent_MR = dir(dataFile);
EventMRdata = importdata(fullfile(pathwayEvent_MR(k).folder,pathwayEvent_MR(k).name));
t0knielinks_MR = EventMRdata.t0;
for j = 1:length(TimeMRlinks) % vervangen door functie
if t0knielinks_MR == TimeMRlinks(j)
t0MR = j
break
end
end
FirstSwingExcursionAngle_MR(k,i)=DataMRlinks(t0MR) % plot
sum = sum + DataMRlinks(t0MR);
end
avgHR(i)=sum/length(pathwayIK_MR) % staaf diagram van het gemiddelde
sum = 0; % i gebruiken voor de juiste subplots
end
filename='AverageFirstSwingExcursionAngle_MR.xlsx'
writematrix(avgMR,filename)