Professor Marwan,
I have some doubts about the use of command line reccurrence plots:
Let's see the use of rp like the official page:
dos('rp -i temp_data -s -m 3 -t 5 -e 0.7 -o rqa.dat');
Doubts:
1) Is this threshold parameter an absolute value of the radius or a percentage value of the maximum diameter of the phase space?
2) Is there a way to vary the parameter m and t? I mean, instead of using "-m 3 -t 7", Could I to use something like that (adapted from the official page):
Code: Select all
% import some data
x=load('soi3.dat');
% data length
N = length(x);
% window length and window step size
w = 100; s = 10;
for i = 1:s:N-w
m=i; tau=round(N/i);
% make sub-string of data
sub_x = x(i:i+w-1);
% save sub-string of data
save('temp_data','-ascii','sub_x')
% call extern rp programme
dos('./rp -i temp_data -s -m m(i) -t tau(i) -e 0.7 -o rqa.dat');
end
delete('temp_data')
Instead of use numbers, could I to use m(i) and tau(i)? I have tested like this and it did not work.
3) The similar use of command line rp in CRP toolbox is like that: (?)
Code: Select all
X=temp_data;
l_min=2; v_min=2;
crqa(X,3,5, 0.7, [], [], l_min, v_min, 1, ...
'euclidean', 'nonormalize', 'silent');
Thanks in advance!
Best regards,
Vinicius Veira.