You can refer to the the code below:
% Define the real and theoretical data tables
real_data = [126.75, 980;
240.50, 1020;
241, 300];
theoretical_data = [126.8, 0, 1;
126.8, 0, 2;
240.5, 1, 1;
243, 1, 2];
% Initialize the combined_data matrix and a flag to mark used theoretical data
combined_data = zeros(0, size(real_data, 2) + size(theoretical_data, 2) - 1);
n_theoretical_pts = size(theoretical_data, 1);
theory_used = zeros(n_theoretical_pts, 1);
% Loop through each row of the real data
for i = 1:size(real_data, 1)
% Compute the difference between real and theoretical data's first column
diff = theoretical_data(:, 1) - real_data(i, 1);
% Find the index of the minimum difference
[min_diff, idx] = min(abs(diff));
% Check if the corresponding theoretical data point has already been used
while theory_used(idx) % while true
% Replace the used point's value in the difference vector with Inf to ignore it
diff(idx) = Inf;
% Find the new minimum difference and its index
[min_diff, idx] = min(abs(diff));
end
% Mark the theoretical data point as used
theory_used(idx) = 1;
% Extract the matched data
real_entry = real_data(i, :);
theory_entry = theoretical_data(idx, 2:end);
% Append the matched data to the combined_data matrix
combined_entry = [real_entry, theory_entry];
combined_data = [combined_data; combined_entry];
end
% Display the resulting combined_data matrix
disp(combined_data);
