Compare commits
No commits in common. "81686bbf0df1926724e259f96559c6b29d9eeb6e" and "115e08f04dc06fa3a6c3c1d82bfd16427eb5ddaf" have entirely different histories.
81686bbf0d
...
115e08f04d
@ -25,9 +25,9 @@ minPeakHeight2 = 1.2*std(magNoGfilt); % based
|
|||||||
[pks, locs] = findpeaks(magNoGfilt, 'MINPEAKHEIGHT', minPeakHeight2); % for step detection
|
[pks, locs] = findpeaks(magNoGfilt, 'MINPEAKHEIGHT', minPeakHeight2); % for step detection
|
||||||
numStepsOption2_filt = numel(pks); % counts number of steps;
|
numStepsOption2_filt = numel(pks); % counts number of steps;
|
||||||
|
|
||||||
diffLocs = diff(locs); % calculates difference in step location
|
diffLocs = diff(locs); % calculates difference in step location
|
||||||
avg_diffLocs = mean(diffLocs); % average distance between steps
|
avg_diffLocs = mean(diffLocs); % average distance between steps
|
||||||
std_diffLocs = std(diffLocs); % standard deviation of distance between steps
|
std_diffLocs = std(diffLocs); % standard deviation of distance between steps
|
||||||
|
|
||||||
figure;
|
figure;
|
||||||
findpeaks(diffLocs, 'MINPEAKHEIGHT', avg_diffLocs, 'MINPEAKDISTANCE',9); % these values have been chosen based on visual inspection of the signal
|
findpeaks(diffLocs, 'MINPEAKHEIGHT', avg_diffLocs, 'MINPEAKDISTANCE',9); % these values have been chosen based on visual inspection of the signal
|
||||||
|
@ -65,9 +65,9 @@ if ApplyRealignment % apply relignment as described in Rispens S, Pijnappels M,
|
|||||||
dataAcc = RealignedAcc;
|
dataAcc = RealignedAcc;
|
||||||
[B,A] = butter(2,20/(FS/2),'low');
|
[B,A] = butter(2,20/(FS/2),'low');
|
||||||
dataAcc_filt = filtfilt(B,A,dataAcc);
|
dataAcc_filt = filtfilt(B,A,dataAcc);
|
||||||
else % we asume that data for CONTROLS is already detrended and in order 1 = AP, 2 = ML, 3 = VT in an earlier stage;
|
else % we asume tat data is already reorderd to 1 = V; 2= ML, 3 = AP in an earlier stage;
|
||||||
[B,A] = butter(2,20/(FS/2),'low');
|
[B,A] = butter(2,20/(FS/2),'low');
|
||||||
dataAcc = inputData(:, [3,2,1]);
|
dataAcc = inputData;
|
||||||
dataAcc_filt = filtfilt(B,A,dataAcc);
|
dataAcc_filt = filtfilt(B,A,dataAcc);
|
||||||
end
|
end
|
||||||
|
|
||||||
@ -101,7 +101,6 @@ if ApplyRemoveSteps
|
|||||||
clear LocsSteps;
|
clear LocsSteps;
|
||||||
|
|
||||||
% Change window length if ApplyRemoveSteps (16-2-2021 LD)
|
% Change window length if ApplyRemoveSteps (16-2-2021 LD)
|
||||||
|
|
||||||
WindowLen = size(dataAccCut,1);
|
WindowLen = size(dataAccCut,1);
|
||||||
|
|
||||||
else;
|
else;
|
||||||
|
Binary file not shown.
@ -1,2 +0,0 @@
|
|||||||
<?xml version='1.0' encoding='UTF-8'?>
|
|
||||||
<Info />
|
|
Loading…
Reference in New Issue
Block a user