27.12.2012 Views

Oscillations, Waves, and Interactions - GWDG

Oscillations, Waves, and Interactions - GWDG

Oscillations, Waves, and Interactions - GWDG

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

attractor in the<br />

unknown<br />

state space M<br />

flow φ t<br />

x<br />

measurement<br />

h: M IR<br />

observable s=h(x)<br />

Complex dynamics of nonlinear systems 421<br />

h<br />

s<br />

F<br />

measured time series<br />

reconstruction<br />

of the attractor in IR<br />

d<br />

Figure 11. Delay reconstruction of states from scalar time series.<br />

t<br />

y<br />

delay coordinates<br />

y(t) = (s(t), s(t - τ), ... , s(t -( d-1)<br />

τ))<br />

d reconstruction dimension<br />

τ delay time<br />

about the reconstruction of states, in particular in the presence of noise, may be<br />

found in Refs. [46,47].<br />

If the data are measured with a high sampling rate Broomhead-King-coordinates [48,<br />

49] may be advantageous. With this method a very high-dimensional reconstruction<br />

is used <strong>and</strong> then a new coordinate system is introduced where the origin is shifted<br />

to the center of mass of the reconstructed states <strong>and</strong> the axes are given by the<br />

(dominant) principal components of the distribution of points (states). This new<br />

coordinate system is based on a Karhunen-Loève transformation 6 that may be computed<br />

by a singular-value decomposition. A discussion of the advantages (e. g., noise<br />

reduction) <strong>and</strong> disadvantages of this “post-processing” of the reconstructed states<br />

may, for example, be found in Ref. [50].<br />

For time series that consist of a sequence of sharp spikes (e. g. from firing neurons)<br />

delay embedding may lead to very inhomogeneous sets of points in the reconstructed<br />

state space that are difficult to analyse. As an alternative one may use in this<br />

case the time intervals between the spikes as components for (re-)constructed state<br />

vectors [51–53].<br />

4.2 Forecasting <strong>and</strong> Modelling<br />

After successful state space reconstruction one can approximate the dynamics in reconstruction<br />

space to forecast or control the underlying dynamical process. Very<br />

simple but efficient algorithms for nonlinear prediction are nearest neighbours methods<br />

(also called local models). Let’s assume that we want to forecast the future<br />

evolution of a given (reconstructed) state for some time horizon T . Available (i. e.,<br />

6 Also called Proper Orthogonal Decomposition (POD) or Principle Component Analysis<br />

(PCA).

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!