function [ESig  ELabel] = segmentEMG(Sig,WTime,STime,MTime,SFreq,MLabel)
% To segment signal by a certain rule to a cell, which is 1*n that n is for
% the number of motions. Every unit is a A * B *C matrix that A stands for 
% the number of data in an analysis window, B stands for the number of
% the channel, C stands for the number of samples in one motion.
% Caution: this function did not compute the transition time between
% motions. So it fits to compute the training data.
% Input variables:
%                Sig is original signal.                              data by channel matrix 
%                WTime is the time of an analysis window.             seconds
%                STime is the time of a sliding window.               seconds
%                MTime is the lasting time of a single motion.        seconds
%                SFreq is the Sample Frequency.
%                MLabel is the label of the motions.   vector
%
% Written by Jiayuan He
% email: jiauyuan.he@uwaterloo.ca

ELabel = 0;
winLen = floor(WTime * SFreq);
sldLen = floor(STime * SFreq);
monLen = floor(MTime * SFreq);
% winLen = WTime * SFreq;
% sldLen = STime * SFreq;
% monLen = MTime * SFreq;
k1 = size(MLabel,1);   % the times of all the motions
k2 = floor((monLen - winLen) / sldLen);  % the times of window sliding in one motion time,not computing the transition of motions

% 剔除动作转换时的数据
% ----------------- segment signal ---------------------
for i = 1: k1,  % fix the motion    
    for j = 0: k2,  % fix the sliding window
          ESig(:,:,(i-1)*(k2+1)+j+1) = Sig(round(((i-1) * monLen + j * sldLen + 1)):round(((i-1) * monLen + j * sldLen + winLen)),:)';
    end
    temp = repmat(MLabel(i),k2+1,1);
    ELabel = [ELabel;temp];
end
ELabel = ELabel(2:size(ELabel,1),:);

%%
%Sample
% for i = 1:10
% [D{i},L{i}] = SegmentEMG(D7((i-1)*27000+1:i*27000,:),0.2,0.025,3,1000,La);
% end;
% tsD6 = cat(3,D{1:i});
% tsL6 = cat(1,L{1:i});