diff --git a/usedMcode/waveformLengthAll.m b/usedMcode/waveformLengthAll.m
index 2effb0c..a8600fd 100644
--- a/usedMcode/waveformLengthAll.m
+++ b/usedMcode/waveformLengthAll.m
@@ -1,8 +1,14 @@
 function [EMG]=waveformLengthAll(sig,frequency,windowEMG,windowShift)
+
+    %Filter around 50Hz and below 2 Hz
+    [A,B]= butter(2,[48 52]/(bci_sf/2),'stop');
+    [C,D]= butter(2,1/(bci_sf/2),'high');
+    [E,F]= butter(2,[148 152]/(bci_sf/2),'stop');
+
     signalWindow=frequency*windowEMG;
     shiftProp=windowEMG/windowShift;
     EMG=zeros((floor(size(sig,1)/signalWindow)-1)*shiftProp+1,size(sig,2));
     parfor i=1:size(sig,2)
-        EMG(:,i)=waveformLength(sig(:,i),frequency,windowEMG,windowShift);
+        EMG(:,i)=waveformLength(filtfilt(double(E),double(F),filtfilt(double(C),double(D),filtfilt(double(A),double(B),double(sig(:,i))))),frequency,windowEMG,windowShift);
     end
 end
\ No newline at end of file