Calculating with different date times

2 ビュー (過去 30 日間)
Christopher
Christopher 2024 年 5 月 27 日
コメント済み: Christopher 2024 年 7 月 6 日
Hi,
I have one array that contains my speed during running and the corresponding time it was recorded at. Measurements are about every one second.
Then I have a second array that contains my heart rate measurements and the corresponding time they were recorded at. These measurements contain data from about every 6-10 seconds.
This is shown by the exemplary screenshot.
What I am trying to do is divide my speed by heart rate at the corresponding times to get something like 'performance'. So I need to average the speed of all measurements that were done for one measurement of heart rate. Something like this:
mean(Speed(1:6))/Heartrate(1)
mean(Speed(7:16))/Heartrate(2)
I cannot do this by hand for the amount of data but have problems coming up with code that selects the correct time range to average the speed (the time range that corresponds to every one measurement of heart rate) and then divides these values.
I appreciate your ideas/ help!
PS Sorry for the bad question title, I couldn't think of a better one :/
  2 件のコメント
Ayush Modi
Ayush Modi 2024 年 5 月 27 日
編集済み: Ayush Modi 2024 年 5 月 27 日
Hi Christopher,
Try this ->
% I have assumed the names of arrays as below
speed_time = duration(15, 15, 24:45);
speed_values = [3, 4, 5, 3, 4, 5, 3, 4, 5, 3, 3, 4, 5, 6, 7, 5, 3, 4, 3, 3, 2, 2]; % I extended the array from the image to show the result
hr_time = duration(15, 15, [24, 30, 39]);
hr_values = [143, 145, 134];
performance = [];
for i = 1:length(hr_time)-1
indices = find(speed_time >= hr_time(i) & speed_time < hr_time(i+1));
mean_speed = mean(speed_values(indices));
performance = [performance, mean_speed / hr_values(i)];
end
% Handling the last heart rate measurement
indices = find(speed_time >= hr_time(end));
mean_speed = mean(speed_values(indices));
performance = [performance, mean_speed / hr_values(end)];
disp(performance);
0.0280 0.0307 0.0235
Christopher
Christopher 2024 年 7 月 6 日
This worked for me. Thanks a lot!

サインインしてコメントする。

回答 (2 件)

Mario Malic
Mario Malic 2024 年 5 月 27 日
編集済み: Mario Malic 2024 年 5 月 27 日
Hello Christopher,
Here's the code that does it
speedTime = datetime("now");
speedTime = dateshift(speedTime, "start", "second", 0:15);
speed = 1:numel(speedTime);
speedTable = table(speedTime', speed');
speedTable.Properties.VariableNames = {'Time', 'Speed'};
heartTime = datetime("now");
heartTime = dateshift(heartTime, "start", "second", 0:5:15);
bpm = linspace(130, 150, numel(heartTime));
heartTable = table(heartTime', bpm');
heartTable.Properties.VariableNames = {'Time', 'BPM'};
head(speedTable)
Time Speed ____________________ _____ 27-May-2024 18:33:05 1 27-May-2024 18:33:06 2 27-May-2024 18:33:07 3 27-May-2024 18:33:08 4 27-May-2024 18:33:09 5 27-May-2024 18:33:10 6 27-May-2024 18:33:11 7 27-May-2024 18:33:12 8
head(heartTable)
Time BPM ____________________ ______ 27-May-2024 18:33:05 130 27-May-2024 18:33:10 136.67 27-May-2024 18:33:15 143.33 27-May-2024 18:33:20 150
for i = 1 : height(heartTable) - 1
idxFrame = isbetween(speedTable.Time, heartTable.Time(i), heartTable.Time(i + 1)); % gets indices between i-th and i+1-th time
if any(idxFrame) % There is overlap in time between two sensors
performance = mean(speedTable.Speed(idxFrame)) / heartTable.BPM(i)
else
error("measurement error")
end
end
performance = 0.0269
performance = 0.0622
performance = 0.0942

Peter Perkins
Peter Perkins 2024 年 5 月 28 日
編集済み: Peter Perkins 2024 年 5 月 28 日
I would think you would want to compute performance at each speed measurement by interpolating heart rate, rather than computing at each heart rate measurement by averaging speed. In any case, either becomes easy with a timetables and retime or synchronize.
Time = datetime(2024,5,28,16,30,0) + seconds(0:30)';
Speed = (1:numel(Time))';
speedTable = timetable(Time,Speed)
speedTable = 31x1 timetable
Time Speed ____________________ _____ 28-May-2024 16:30:00 1 28-May-2024 16:30:01 2 28-May-2024 16:30:02 3 28-May-2024 16:30:03 4 28-May-2024 16:30:04 5 28-May-2024 16:30:05 6 28-May-2024 16:30:06 7 28-May-2024 16:30:07 8 28-May-2024 16:30:08 9 28-May-2024 16:30:09 10 28-May-2024 16:30:10 11 28-May-2024 16:30:11 12 28-May-2024 16:30:12 13 28-May-2024 16:30:13 14 28-May-2024 16:30:14 15 28-May-2024 16:30:15 16
Time = datetime(2024,5,28,16,30,0) + seconds(0:5:30)';
BPM = linspace(130, 150, numel(Time))';
heartTable = timetable(Time,BPM)
heartTable = 7x1 timetable
Time BPM ____________________ ______ 28-May-2024 16:30:00 130 28-May-2024 16:30:05 133.33 28-May-2024 16:30:10 136.67 28-May-2024 16:30:15 140 28-May-2024 16:30:20 143.33 28-May-2024 16:30:25 146.67 28-May-2024 16:30:30 150
Version 1: average the speed every 5 seconds
data = synchronize(speedTable,heartTable,"last","mean");
data.Performance = data.Speed ./ data.BPM
data = 7x3 timetable
Time Speed BPM Performance ____________________ _____ ______ ___________ 28-May-2024 16:30:00 3 130 0.023077 28-May-2024 16:30:05 8 133.33 0.06 28-May-2024 16:30:10 13 136.67 0.095122 28-May-2024 16:30:15 18 140 0.12857 28-May-2024 16:30:20 23 143.33 0.16047 28-May-2024 16:30:25 28 146.67 0.19091 28-May-2024 16:30:30 31 150 0.20667
Version 2: Interpolate the BPM at each second
data = synchronize(speedTable,heartTable,"first","spline");
data.Performance = data.Speed ./ data.BPM
data = 31x3 timetable
Time Speed BPM Performance ____________________ _____ ______ ___________ 28-May-2024 16:30:00 1 130 0.0076923 28-May-2024 16:30:01 2 130.67 0.015306 28-May-2024 16:30:02 3 131.33 0.022843 28-May-2024 16:30:03 4 132 0.030303 28-May-2024 16:30:04 5 132.67 0.037688 28-May-2024 16:30:05 6 133.33 0.045 28-May-2024 16:30:06 7 134 0.052239 28-May-2024 16:30:07 8 134.67 0.059406 28-May-2024 16:30:08 9 135.33 0.066502 28-May-2024 16:30:09 10 136 0.073529 28-May-2024 16:30:10 11 136.67 0.080488 28-May-2024 16:30:11 12 137.33 0.087379 28-May-2024 16:30:12 13 138 0.094203 28-May-2024 16:30:13 14 138.67 0.10096 28-May-2024 16:30:14 15 139.33 0.10766 28-May-2024 16:30:15 16 140 0.11429
You may choose to bin things differently at the edges, but that's the general idea.

カテゴリ

Help Center および File ExchangeMatrix Indexing についてさらに検索

タグ

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by