Physics:Takens's theorem

From HandWiki
Short description: Conditions under which a chaotic system can be reconstructed by observation
Rössler attractor reconstructed by Taken's theorem, using different delay lengths. Orbits around the attractor have a period between 5.2 to 6.2.

In the study of dynamical systems, a delay embedding theorem gives the conditions under which a chaotic dynamical system can be reconstructed from a sequence of observations of the state of that system. The reconstruction preserves the properties of the dynamical system that do not change under smooth coordinate changes (i.e., diffeomorphisms), but it does not preserve the geometric shape of structures in phase space.

Takens' theorem is the 1981 delay embedding theorem of Floris Takens. It provides the conditions under which a smooth attractor can be reconstructed from the observations made with a generic function. Later results replaced the smooth attractor with a set of arbitrary box counting dimension and the class of generic functions with other classes of functions.

It is the most commonly used method for attractor reconstruction.[1]

Delay embedding theorems are simpler to state for discrete-time dynamical systems. The state space of the dynamical system is a ν-dimensional manifold M. The dynamics is given by a smooth map

[math]\displaystyle{ f: M \to M. }[/math]

Assume that the dynamics f has a strange attractor [math]\displaystyle{ A\sub M }[/math] with box counting dimension dA. Using ideas from Whitney's embedding theorem, A can be embedded in k-dimensional Euclidean space with

[math]\displaystyle{ k \gt 2 d_A. }[/math]

That is, there is a diffeomorphism φ that maps A into [math]\displaystyle{ \R^k }[/math] such that the derivative of φ has full rank.

A delay embedding theorem uses an observation function to construct the embedding function. An observation function [math]\displaystyle{ \alpha : M \to \R }[/math] must be twice-differentiable and associate a real number to any point of the attractor A. It must also be typical, so its derivative is of full rank and has no special symmetries in its components. The delay embedding theorem states that the function

[math]\displaystyle{ \varphi_T(x) = \bigl(\alpha(x), \, \alpha(f(x)), \, \dots, \, \alpha(f^{k-1}(x)) \, \bigr) }[/math]

is an embedding of the strange attractor A in [math]\displaystyle{ \R^k. }[/math]

Simplified version

Suppose the [math]\displaystyle{ d }[/math]-dimensional state vector [math]\displaystyle{ x_t }[/math] evolves according to an unknown but continuous and (crucially) deterministic dynamic. Suppose, too, that the one-dimensional observable [math]\displaystyle{ y }[/math] is a smooth function of [math]\displaystyle{ x }[/math], and “coupled” to all the components of [math]\displaystyle{ x }[/math]. Now at any time we can look not just at the present measurement [math]\displaystyle{ y(t) }[/math], but also at observations made at times removed from us by multiples of some lag [math]\displaystyle{ \tau: y_{t+\tau}, y_{t+2\tau} }[/math], etc. If we use [math]\displaystyle{ k }[/math] lags, we have a [math]\displaystyle{ k }[/math]-dimensional vector. One might expect that, as the number of lags is increased, the motion in the lagged space will become more and more predictable, and perhaps in the limit [math]\displaystyle{ k \to \infty }[/math] would become deterministic. In fact, the dynamics of the lagged vectors become deterministic at a finite dimension; not only that, but the deterministic dynamics are completely equivalent to those of the original state space (precisely, they are related by a smooth, invertible change of coordinates, or diffeomorphism). In fact, the theorem says that determinism appears once you reach dimension [math]\displaystyle{ 2d+1 }[/math], and the minimal embedding dimension is often less.[2][3]

Choice of delay

Takens' theorem is usually used to reconstruct strange attractors out of experimental data, for which there is contamination by noise. As such, the choice of delay time becomes important. Whereas for data without noise, any choice of delay is valid, for noisy data, the attractor would be destroyed by noise for delays chosen badly.

The optimal delay is typically around one-tenth to one-half the mean orbital period around the attractor.[4][5]

See also

References

  1. Sauer, Timothy D. (2006-10-24). "Attractor reconstruction" (in en). Scholarpedia 1 (10): 1727. doi:10.4249/scholarpedia.1727. ISSN 1941-6016. 
  2. Shalizi, Cosma R. (2006). "Methods and Techniques of Complex Systems Science: An Overview". in Deisboeck, ThomasS; Kresh, J.Yasha. Complex Systems Science in Biomedicine. Topics in Biomedical Engineering International Book Series. Springer US. pp. 33–114. doi:10.1007/978-0-387-33532-2_2. ISBN 978-0-387-30241-6. https://archive.org/details/complexsystemssc00kres. 
  3. Barański, Krzysztof; Gutman, Yonatan; Śpiewak, Adam (2020-09-01). "A probabilistic Takens theorem". Nonlinearity 33 (9): 4940–4966. doi:10.1088/1361-6544/ab8fb8. ISSN 0951-7715. https://iopscience.iop.org/article/10.1088/1361-6544/ab8fb8. 
  4. Strogatz, Steven (2015). "12.4 Chemical chaos and attractor reconstruction". Nonlinear dynamics and chaos: with applications to physics, biology, chemistry, and engineering (Second ed.). Boulder, CO. ISBN 978-0-8133-4910-7. OCLC 842877119. 
  5. Fraser, Andrew M.; Swinney, Harry L. (1986-02-01). "Independent coordinates for strange attractors from mutual information". Physical Review A 33 (2): 1134–1140. doi:10.1103/PhysRevA.33.1134. PMID 9896728. https://link.aps.org/doi/10.1103/PhysRevA.33.1134. 

Further reading

External links

  • [1] Scientio's ChaosKit product uses embedding to create analyses and predictions. Access is provided online via a web service and graphic interface.
  • [2] Empirical Dynamic Modelling tools pyEDM and rEDM use embedding for analyses, prediction, and causal inference.