Please share the algorithm of the particular data set. I want to find out average of the attached data of every 5 seconds. Thank You

1 次查看（过去 30 天）

显示更早的评论

SATYA PAL 2023-6-20

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1985794-please-share-the-algorithm-of-the-particular-data-set-i-want-to-find-out-average-of-the-attached-da

回答： Steven Lord 2024-2-15

采纳的回答： Mathieu NOE

Datas ain Excel are attached

7 个评论
显示 5更早的评论隐藏 5更早的评论

SATYA PAL 2023-6-22

Dear Sir

I want the averages of discrete blocks over 5-second intervals.

SATYA PAL 2023-6-22

but the number of data in each interval of 5 seconds are not same

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Mathieu NOE 2023-6-22

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1985794-please-share-the-algorithm-of-the-particular-data-set-i-want-to-find-out-average-of-the-attached-da#answer_1260499

编辑：Mathieu NOE 2023-6-22

在 MATLAB Online 中打开

hello

try this

data = readmatrix('Data.csv'); % Force (N),Position (mm),Time (sec)
Time = data(:,3);
t = linspace(Time(1),Time(end),size(data,1)); % recompute the time vector because the original data is rounded 
dt = mean(diff(t));
%% home made solution (you choose the amount of overlap)
buffer_size = round(5/dt); % how many samples for 5 seconds buffer ? 
overlap = 0; % overlap expressed in samples
%%%% main loop %%%%
[new_time,data_out] = my_movmean(t,data(:,1:2),buffer_size,overlap);
figure(2),
plot(t,data(:,1),new_time,data_out(:,1),'*-r');
title('Force');
legend('raw data','5s mean');
xlabel('Time(s)');
ylabel('(N)');
figure(3),
plot(t,data(:,2),new_time,data_out(:,2),'*-r');
title('Position');
legend('raw data','5s mean');
xlabel('Time(s)');
ylabel('(mm)');
%%%%%%%%%% my functions %%%%%%%%%%%%%%
function [new_time,data_out] = my_movmean(t,data_in,buffer_size,overlap)
% NB : buffer size and overlap are integer numbers (samples)
% data (in , out) are 1D arrays (vectors)
shift = buffer_size-overlap;    % nb of samples between 2 contiguous buffers  
[samples,~] = size(data_in);
nb_of_loops = fix((samples-buffer_size)/shift +1);
    for k=1:nb_of_loops
        start_index = 1+(k-1)*shift;
        stop_index = min(start_index+ buffer_size-1,samples);
        x_index(k) = round((start_index+stop_index)/2);  
        data_out(k,:) = mean(data_in(start_index:stop_index,:),1,'omitnan');  % 
        
    end
new_time = t(x_index); % time values are computed at the center of the buffer
end