clc;
clear all;
close all;
id = fopen('LPdas_1636715098861_1.Data');
data = fread(id,'float');
[b,a] = butter(4, [0.5/125,20/125],'bandpass');
fd = filtfilt(b,a,data);
cut = fd(1:4000);
figure,plot(cut);
figure,findpeaks(cut,'minpeakdistance',100);
[pks,loc] = findpeaks(cut,'minpeakdistance',100);
delta = diff(loc)