Computing Laminarity

Everything about quantification of recurrence plots.

Computing Laminarity

Postby holistic » Wed Feb 19, 2014 12:39


I have 10 timeseries and I want to compute the Lamniarity for all of those without using the gui-interface. Unfortunately I don't know how the commands are to save to compute the crp and the LAM and save it into an array. Can anybody help me out with which matlab code to use?
Posts: 16
Joined: Wed Feb 19, 2014 09:53

Re: Computing Laminarity

Postby Norbert » Wed Feb 19, 2014 21:51

There are different ways, depending whether the time series are in files or are already in your Matlab workspace. Let's start simple and assume your time series are the columns of an array:

x = [x1, x2, x3, … x10];
size(x) gives something like N, 10, where N is the length of your time series.

Now you can use
Code: Select all
for i = 1:10
   Rtemp = crqa(x(:,i),m,t,e,'silent');
   LAM(i) = Rtemp(6);
(where m and t are the embedding parameters and e the recurrence threshold which you have to chose)

In case you need the time resolved values, i.e., applying moving windows, simply use
Code: Select all
for i = 1:10
   Rtemp = crqa(x(:,i),m,t,e,w,ws,'silent');
   LAM(:,i) = Rtemp(:,6);
(where w and ws are window size and window step)
User avatar
Posts: 173
Joined: Wed Jan 4, 2006 11:03
Location: Potsdam Institute for Climate Impact Research, Germany

Return to Quantification

Who is online

Users browsing this forum: No registered users and 1 guest