hello
I was not able to understand the second part of your code
If the target was to compute the RPM value vs time , look here after
see the attached function used here
FYI, I had to "regenerate" the time values (assuming this is sampled at fixed sampling frequency) , because du to rounding (short format) the time precision was insufficient to compute the RPM with good accuracy
Time= tsout.time;
Enc= tsout.Data(:,1);
samples = length(Time);
dt = (Time(samples) - Time(1))/ (samples-1);
Fs = 1/dt;
Time2 = (0:samples-1)*dt;
n = 1500;
% find crossing points / positive slope / threshold = 2.5 V)
threshold = 2.5 ;
[~,t0_pos,s0_pos,~,~,~] = crossing_V7(Enc(1:n),Time2(1:n),threshold,'linear');
figure(1),
plot (Time2(1:n),Enc(1:n),'b',t0_pos,s0_pos,'r+');
delta_time = diff(t0_pos);
rotation_freq = 1./delta_time;
RPM = 60*rotation_freq;
figure(2),
plot (t0_pos(2:end),RPM,'b');
xlabel('time (s)');
ylabel('RPM');
