Dividing cyclical data in array

Hello,

I have acceleration data of diferent segments of a moving horse. I would like to split it so that I have a different array for each of the strides of the horse, so that I can later time normalize each stride and average them to obtain one averaged stride. I can identify each stride by locating the peaks of acceleration of the horse's foot.

Is there an elegant way to split the array in such way?

Thank you so much.

Answer 1

[Mohammad Sami]

Assuming you can get the locations of the peak, you can create an id variable.

% acc = ... m x 1 array
%locationidxofpeak = somefuntion....
strideid = zeros(length(acc),1);
strideid(locationidxofpeak) = 1;
strideid = cumsum(strideid);
% now stride id would be like [0 0 0 0 0 1 1 1 1 1 2 2 2 2 2 ....]
for i = min(strideid):max(strideid)
    strideacc = acc(strideid == i);
    % your code.
end

[Marc Elmeua]

wow! this worked great! shokran