I need to calculate joint probability of reccurence, so first I have to calculate delayed JRP. According to Norbert's article I have found, I have to calculate delayed reccurence matrix R[i+t,j+t]. My matlab code:
Code: Select all
d = 15;
S= zeros(N-d);
for i= 1 : N-d
for j = 1 : N-d
S(i+d,j+d) = sqrt(sum((x(i,:)-x(j,:)).^2,2)); ?
end
end
figure
imagesc(t(1:N-d), t(1:N-d), S < e)
Code: Select all
d = 15;
x1d = circshift(x1,d);
x2d = circshift(x2,d);
Sd = sqrt(sum( (x1d - x2d) .^ 2, 2 ));
Sd = reshape(Sd, N2, N2);
figure
imagesc(t(1:N2), t(1:N2), Sd < 2)
axis square
colormap([1 1 1;0 0 0])
xlabel('Time (sec)'), ylabel('Time (sec)')