Geoscience Reference
In-Depth Information
Fig. 5.24
h e recurrence plot for the synthetic data derived from the distance matrix, as
shown in Fig. 5.23, at er applying a threshold of สต=1.
trajectory with a fully vectorized algorithm supported by MATLAB. For this
we need to transfer the trajectory vector into two test vectors whose element-
wise test will provide the pair-wise test of the trajectory vector:
x1 = repmat(xe,N2,1);
x2 = reshape(repmat(xe(:),1,N2)',N2*N2,m);
From these vectors we calculate the recurrence plot using the Euclidean norm
without any
FOR
loop (see Section 9.4 for details on Euclidean distances).
S = sqrt(sum((x1 - x2).^ 2,2 ));
S = reshape(S,N2,N2);
imagesc(t(1:N2),t(1:N2),S<10)
colormap([1 1 1;0 0 0])
xlabel('Time'), ylabel('Time')
axis xy
h is recurrence plot reveals many short diagonal lines (Fig. 5.25). h ese
