Kind of long code
clc, clear
noise = 0.1;
n = 10; % number of areas
jump = 1;
cur_jump = 0;
t = linspace(0,30,500);
signal = round(t./10) + rand(1,length(t))*noise;
ds = diff(signal);
t = t(abs(ds)<noise);
s = signal(abs(ds)<noise);
ar = zeros(length(signal),n);
time = zeros(length(signal),n);
ns = 0;
while cur_jump < max(s)
nj = 0;
ns = ns + 1;
for i = 1:length(s)
if cur_jump-jump/2 < s(i) && s(i) < cur_jump+jump/2
nj = nj + 1;
ar(nj,ns) = s(i);
time(nj,ns) = t(i);
end
end
cur_jump = cur_jump + jump;
end
subplot(3,2,1)
plot(t,s)
hold on
for i = 1:ns
subplot(3,2,i+1)
plot(time(:,i), ar(:,i),'r')
end
hold off
