en>Selinger: /* References */ clarified Benabou citation

2008-12-22T21:12:34Z

References: clarified Benabou citation

New page

In descriptive [[statistics]] and [[chaos theory]], a '''recurrence plot''' ('''RP''') is a plot showing, for a given moment in time, the times at which a [[phase space]] trajectory visits roughly the same area in the [[phase space]]. In other words, it is a graph of

:<math>\vec{x}(i)\approx \vec{x}(j),\,</math>

showing <math>i</math> on a horizontal axis and <math>j</math> on a vertical axis, where <math>\vec{x}</math> is a [[phase space]] trajectory.

==Background==

Natural processes can have a distinct recurrent behaviour, e.g. periodicities (as [[Seasonality|seasonal]] or [[Milankovich cycle]]s), but also irregular cyclicities (as [[El Niño]] Southern Oscillation). Moreover, the recurrence of states, in the meaning that states are again arbitrarily close after some time of [[divergence]], is a fundamental property of [[deterministic]] [[dynamical systems]] and is typical for [[nonlinear]] or [[chaotic system]]s (cf. [[Poincaré recurrence theorem]]). The recurrence of states in nature has been known for a long time and has also been discussed in early work (e.g. [[Henri Poincaré]] 1890).

==Detailed description==

Eckmann et al. (1987) introduced recurrence plots, which provide a way to visualize the periodic nature of a trajectory through a [[phase space]]. Often, the phase space does not have a low enough dimension (two or three) to be pictured, since higher-dimensional phase spaces can only be visualized by projection into the two or three-dimensional sub-spaces. However, making a recurrence plot enables us to investigate certain aspects of the ''m''-dimensional phase space trajectory through a two-dimensional representation.

A '''recurrence''' is a time the trajectory returns to a location it has visited before. The recurrence plot depicts the collection of pairs of times at which the trajectory is at the same place, i.e. the set of <math>(i,j)</math> with <math>\vec{x}(i) = \vec{x}(j)</math>. This can show many things, for instance, if the trajectory is strictly periodic with period <math>T</math>, then all such pairs of times will be separated by a multiple of <math>T</math> and visible as diagonal lines. To make the plot, continuous time and continuous phase space are discretized, taking e.g. <math>\vec{x}(i)</math> as the location of the trajectory at time <math>i \tau</math> and counting as a recurrence any time the trajectory gets sufficiently close (say, within ε) to a point it has been previously. Concretely then, recurrence/non-recurrence can be recorded by the binary function

:<math>R(i,j) = \begin{cases} 1 \quad &\text{if} \quad \| \vec{x}(i) - \vec{x}(j)\| \le \varepsilon \\ 0 \quad & \text{otherwise}, \end{cases}</math>

and the recurrence plot puts a (black) point at coordinates <math>(i,j)</math> if <math>R(i,j)=1</math>.

Caused by characteristic behaviour of the phase space trajectory, a recurrence plot contains typical small-scale structures, as single dots, diagonal lines and vertical/horizontal lines (or a mixture of the latter, which combines to extended clusters). The large-scale structure, also called ''texture'', can be visually characterised by ''homogenous'', ''periodic'', ''drift'' or ''disrupted''. The visual appearance of an RP gives hints about the dynamics of the system.

[[Image:Rp examples740.gif|thumb|center|740px|Typical examples of recurrence plots (top row: [[time series]] (plotted over time); bottom row: corresponding recurrence plots). From left to right: uncorrelated stochastic data ([[white noise]]), [[harmonic oscillation]] with two frequencies, chaotic data with linear trend ([[logistic map]]) and data from an [[autoregressive process|auto-regressive process]].]]

The small-scale structures in RPs are used by the [[recurrence quantification analysis]] (Zbilut & Webber 1992; Marwan et al. 2002). This quantification allows to describe the RPs in a quantitative way, and to study transitions or nonlinear parameters of the system. In contrast to the heuristic approach of the recurrence quantification analysis, which depends on the choice of the embedding parameters, some [[dynamical invariant]]s as [[correlation dimension]], [[K2 entropy]] or [[mutual information]], which are independent on the embedding, can also be derived from recurrence plots. The base for these dynamical invariants are the recurrence rate and the distribution of the lengths of the diagonal lines.

Close returns plots are similar to recurrence plots. The difference is that the relative time between recurrences is used for the <math>y</math>-axis (instead of absolute time).

The main advantage of recurrence plots is that they provide useful information even for short and non-stationary data, where other methods fail.

==Extensions==

Multivariate extensions of recurrence plots were developed as '''cross recurrence plots''' and '''joint recurrence plots'''.

Cross recurrence plots consider the phase space trajectories of two different systems in the same phase space (Marwan & Kurths 2002):

:<math>\mathbf{CR}(i,j) = \Theta(\varepsilon - \| \vec{x}(i) - \vec{y}(j)\|), \quad \vec{x}(i),\, \vec{y}(i) \in \Bbb{R}^m, \quad i=1, \dots, N_x, \ j=1, \dots, N_y.</math>

The dimension of both systems must be the same, but the number of considered states (i.e. data length) can be different. Cross recurrence plots compare the simultaneous occurrences of ''similar states'' of two systems. They can be used in order to analyse the similarity of the dynamical evolution between two different systems, to look for similar matching patterns in two systems, or to study the time-relationship of two similar systems, whose time-scale differ (Marwan & Kurths 2005).

Joint recurrence plots are the [[Matrix product#Hadamard product|Hadamard product]] of the recurrence plots of the considered sub-systems (Romano et al. 2004), e.g. for two systems <math>\vec{x}</math> and <math>\vec{y}</math> the joint recurrence plot is

:<math>\mathbf{JR}(i,j) = \Theta(\varepsilon_x - \| \vec{x}(i) - \vec{x}(j)\|) \cdot \Theta(\varepsilon_y - \| \vec{y}(i) - \vec{y}(j)\|), \quad \vec{x}(i) \in \Bbb{R}^m, \quad \vec{y}(i) \in \Bbb{R}^n,\quad i,j=1, \dots, N_{x,y}.</math>

In contrast to cross recurrence plots, joint recurrence plots compare the simultaneous occurrence of ''recurrences'' in two (or more) systems. Moreover, the dimension of the considered phase spaces can be different, but the number of the considered states has to be the same for all the sub-systems. Joint recurrence plots can be used in order to detect [[phase synchronization|phase synchronisation]].

==Example==

[[Image:Rp soi.gif|thumb|center|377px|Recurrence plot of the [[Southern Oscillation]] index.]]

== See also==
* [[Recurrence quantification analysis]], a heuristic approach to quantify recurrence plots.
* [[Recurrence period density entropy]], an information-theoretic method for summarising the recurrence properties of both deterministic and stochastic dynamical systems.
* [[Self-similarity matrix]]
* [[Poincaré plot]]

{{More footnotes|date=December 2009}}

==References==
* {{cite journal
| author=J. P. Eckmann, S. O. Kamphorst, [[David Ruelle|D. Ruelle]]
| title=Recurrence Plots of Dynamical Systems
| journal=Europhysics Letters
| volume=5
| issue=9
| pages=973–977
| year=1987
| url=http://iopscience.iop.org/0295-5075/4/9/004/
| doi=10.1209/0295-5075/4/9/004
|bibcode = 1987EL......4..973E }}
* {{cite journal
| author=N. Marwan, M. C. Romano, M. Thiel, J. Kurths
| title=Recurrence Plots for the Analysis of Complex Systems
| journal=Physics Reports
| volume=438
| issue=5-6
| year=2007
| doi=10.1016/j.physrep.2006.11.001
| pages=237
|bibcode = 2007PhR...438..237M }}
* {{cite journal
| author=N. Marwan
| title=A historical review of recurrence plots
| journal=European Physical Journal ST
| volume=164
| issue=1
| year=2008
| pages=3–12
| url=http://www.springerlink.com/content/5412256211633127/
| doi=10.1140/epjst/e2008-00829-1
|bibcode = 2008EPJST.164....3M }}

== External links==
* http://www.recurrence-plot.tk/
* http://www.scitopics.com/Recurrence_plot.html

[[Category:Plots (graphics)]]
[[Category:Signal processing]]
[[Category:Dynamical systems]]
[[Category:Visualization (graphic)]]
[[Category:Chaos theory]]
[[Category:Scaling symmetries]]

Bicategory - Revision history

en>Selinger: /* References */ clarified Benabou citation