Averaging data
1 次查看(过去 30 天)
显示 更早的评论
Hello,
I'm having difficult time trying to average my data. I have data which look like this
time1 = [0.14 0.99 1.67 2.45 3.10 3.70 4.24 5.00]
res1 = [4 6 2 9 1 0 4 7]
time2 = [0.11 0.69 1.00 1.45 1.66 2.04 2.24 2.99 3.11 3.43 4.25
4.55 5.00]
res2 = [8 1 5 3 7 1 3 10 9 5 3 2 1]
time3 = [0.09 0.33 1.13 1.44 2.10 2.70 3.24 4.00 4.80 5.00]
res3 = [0 3 4 7 2 6 3 4 9 8]
The period of time is always 5 min but there is different number of responses given during this time. These responses showing only a change at particular time and the value between time intervals is constant. I'd like to plot average of the data using a stairs function. Any help highly appreciated.
0 个评论
采纳的回答
Oleg Komarov
2011-8-9
time1 = [0.14 0.99 1.67 2.45 3.10 3.70 4.24 5.00];
res1 = [4 6 2 9 1 0 4 7];
axis([0 5 0 10])
hold on
% Plot stairs
stairs(time1,res1)
% Plot average
avg = sum(diff([0 time1])/time1(end) .* res1);
line([0,time1(end)],[avg avg],'Color','r')
One clarification, is the first element of res1 - 4 - valid from 0 to 0.14 or as it is it's ok?
EDIT
% unique time vector
untime = unique([time1 time2 time3]);
% Idx
[trash,idx1] = histc(untime,[0 time1(1:end-1) inf]);
[trash,idx2] = histc(untime,[0 time2(1:end-1) inf]);
[trash,idx3] = histc(untime,[0 time3(1:end-1) inf]);
rep = [res1(idx1); res2(idx2); res3(idx3)];
avg = mean(rep);
hold on
% Plot all the lines
stairs(untime,rep.')
% Plot the average in black
stairs([0 untime], [9 avg],'Color','k','Linew',2)
8 个评论
更多回答(1 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Data Type Conversion 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!