Compare commits
2 Commits
115e08f04d
...
81686bbf0d
Author | SHA1 | Date | |
---|---|---|---|
81686bbf0d | |||
ea3a7f95e8 |
@ -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 tat data is already reorderd to 1 = V; 2= ML, 3 = AP in an earlier stage;
|
else % we asume that data for CONTROLS is already detrended and in order 1 = AP, 2 = ML, 3 = VT in an earlier stage;
|
||||||
[B,A] = butter(2,20/(FS/2),'low');
|
[B,A] = butter(2,20/(FS/2),'low');
|
||||||
dataAcc = inputData;
|
dataAcc = inputData(:, [3,2,1]);
|
||||||
dataAcc_filt = filtfilt(B,A,dataAcc);
|
dataAcc_filt = filtfilt(B,A,dataAcc);
|
||||||
end
|
end
|
||||||
|
|
||||||
@ -101,6 +101,7 @@ 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.
@ -0,0 +1,2 @@
|
|||||||
|
<?xml version='1.0' encoding='UTF-8'?>
|
||||||
|
<Info />
|
Loading…
Reference in New Issue
Block a user