Identification of Structural Model for Chaotic Systems

This article is talking about the study constructive method of structural identification systems with chaotic dynamics. It is shown that the reconstructed attractors are a source of information not only about the dynamics but also on the basis of the attractors which can be identified and the mere sight of models. It is known that the knowledge of the symmetry group allows you to specify the form of a minimal system. Forming a group transformation can be found in the recon-structed attractor. The affine system as the basic model is selected. Type of a nonlinear system is the subject of calcula-tions. A theoretical analysis is performed and proof of the possibility of constructing models in the central invariant manifold reduced. This developed algorithm for determining the observed symmetry in the attractor. The results of identification used in real systems are an application.


Introduction
The first work on the reconstruction of the strange attractor from the time series has been publishing the results on hydrodynamics [1]. The article shows that you can get a satisfactory picture of the strange attractor of the geometric dimensions of a small, if the variables x, appearing in the equations of the dynamical system   d d x t F x  , use the m-dimensional vectors, derived from the elements of time series of the same principle, which in the problems of autoregression. That same year, F. Takens reported on his theorem, which was published a year later [2]. That it is the basis of all algorithms for time series analysis methods of nonlinear dynamics. The problem of determining the form of a dynamical system from its one-dimensional realization belongs to a class of incorrect problems. Unlike the problem of analyzing this issue is ambiguous, since there are infinitely many dynamical systems of various kinds which can play the existing signal with a given degree of accuracy.
The method of the global reconstruction of a dynamic system of equations for its one-dimensional realization was proposed in [3,4]. The algorithm is as follows. One-dimensional realization of the process is in a system, which is considered a "black box" recovered phase portrait on the Takens theorem, topologically equivalent to the attractor of the original system. A priori given equation is the method of least squares which is a set of un-known coefficients. Now there are considerable publications, developing and constantly improving the proposed the method [5][6][7][8][9][10][11]. For example, in R. Brown and others [12] to reconstruction dynamic equations on the experimental time series with a broadband continuous spectrum use additional information about the dynamic and statistical properties of the original system contained in the implementation. In obtaining equation takes into account the values of Lyapunov exponents and the probability density, calculated from the original time series. However, the resulting evolution equations have a very cumbersome, inconvenient to use. In [13] used the hidden variables to write a model equation. In [12] describes a method for synchronizing the model with the original data. In a number of O.L. Anosov proposed reconstruction algorithm scalar differential equations for systems with delay.
However, the feature of many studies is that the proposed methods are illustrated with examples of simple low-dimensional model systems when we know in advance what should be the result of global reconstruction. It does not show substantial benefits given by the modifications of the basic method [3]. Described in the publications of the algorithms are tested on a number of known model systems with small dimension and a simple form right sides. The efficiency of the method is demonstrated by the example of time series generated by the real "black boxes".
The challenge is the need to work with noisy data in the processing of experimental time series. On the one hand, more desirable is the use of sequential differentiation to restore the phase trajectory, because it can get a model that contains, in general, approximately n times smaller than the coefficients of the various non-linear than when using the method of delays. But differentiation will inevitably lead to increased noise components of high order. Without pre-filtering the time dependence of the second derivative can be noise-like process. In addition, methods of attachment are obvious flaws in the analysis significantly heterogeneous implementations, i.e., signals in which areas with fast motion alternating with areas of slow motions.
Arbitrary choice of nonlinearities, as a rule, does not allow for a successful reconstruction of the dynamic equations for real systems. In particular, in [14] indicated the presence of three typical cases: 1) Reconstruction of locally describe the phase trajectory of the original system. In this case, the reconstructed model is unstable in the sense that the solution of these equations reproduces the signal under investigation only in a short period of time.
2) There is poor local predictability of the phase trajectory, but there is visual similarity of the phase portraits. Reconstruction of solution is stable in the sense of Poisson. In this case, the attractor of the reconstructed model has metric characteristics similar to those of the original attractor.
3) There is a good local predictability of the phase trajectory of any point in time values exceeding the characteristic correlation time. The phase portrait reconstructed model is identical to the original, and the system is Poisson stable.

Group Properties for the Construction of Equivalent Mapping
Consider the research questions of the global structure of orbits of dynamical systems, which do not depend on the choice of coordinate system, i.e., admit a symmetry transformation. From a global point of view, the change of coordinates is a diffeomorphism (in the case of a smooth structure) or a homeomorphism (in the topological situation) phase spaces. Thus, we can introduce a natural equivalence relation between the dynamical systems associated with different classes of coordinate changes, and interpret the problem of describing the structure of orbits as a classification problem of dynamical systems up to these equivalence relations. Given the equivalence of the trajectories on the torus [15], we can show that every periodic point of the map, whose spectrum does not contain a unit, defines certain -modules. Since these periodic orbits are separated from each other, their spectra can be outraged, regardless, at least, for any finite set of points. The modules from a variety of periodic orbits a independent. 1


On the other hand, at least in some cases, the local spectrum is a complete invariant with respect to a smooth conjugacy. Different approaches to the problem of local smooth conjugacy shown in [16].
For the global structure of orbits there are ways to build structures of independent modules associated with periodic orbits. In the case of infinitely many periodic orbits, like expanding mapping and hyperbolic automorphism of the torus, there are infinitely many invariants of the local -equivalence. For these two cases, which represent the simplest examples of hyperbolic systems, the spectra of periodic orbits form a complete system of invariants for -and which can be extended to the whole real line by the formula Thus, the two diffeomorphism  1 and  2 , for which for some s in [0,1]. Constructed above the map transition depends on the factors that determine the linearization  in the vicinity of the endpoints, as well as the choice of the reference point а as well. It is obvious that the change does not alter the linearization map transition if a amended accordingly. Changing the reference point leads to the substitution mapping the transition to the equiva-lent.
The notion of smooth equivalence in accordance with [16] can be easily moved in case of continuous time. Flow equivalence is a contingency flows as a differentiable actions of a group of real numbers. The structure of the orbits of the flow, unlike the case of discrete-time systems, has two distinct aspects: 1) The relative behavior of points on different orbits and 2) evolution of the initial condition along the orbit over time.
There is a natural way to change the thread, keeping the first aspect of the structure of its orbits, and it does not change its orbit.


For the two streams t  and t  will define the time change by using the vector fields of infinitesimal generator of the symmetry group: From the uniqueness of solutions of differential equations, it follows that the zeros of the vector fields are fixed points of the corresponding flow. Thus, we con- is not a fixed point, then the tangent vector to the curve does not vanish and have the same direction.
It is natural to try to describe all the changes of time given flow modulo trivial substitutions. This problem is essentially equivalent to the problem of describing the space of all sufficiently smooth positive functions up to adding features that are derivatives of other smooth functions in the direction of flow.
It makes sense when determining equivalence to consider the conservation of structural stability for flows. One way of determining is based on the equivalence of all the perturbations [16]. Not being fully degenerate, this requirement is rarely performed, for example, in the presence of periodic orbits of periods are modules in this sense. We assume that the local topological equivalence should preserve the structural stability, thus transforming the homeomorphism can be quite close to the identity for small perturbations.
For all the above concepts in this section of the compact phase spaces is immaterial. In addition, a natural way to modify these definitions for the cases where for some points of the dynamical system is defined only for a finite time interval. Such a generalization leads to the concepts of local and semilocal structural stability.
For any contraction phase space can not have a smooth structure, so that these concepts are not directly applicable. However, the contraction mapping of a small disc in the Euclidean space, structural stability, as well as a hyperbolic linear map in the vicinity of the fixed point.
Symmetry rotation is not structurally stable. Since the topological conjugacy preserves the periodic orbits, and the transformation of rotation by an irrational angle can not be associated with the transformation of rotation, for which the corresponding number is rational. But since the rational numbers, and the set of irrational numbers are dense, then some arbitrarily small perturbations of the turn to rational angle you can find the rotation by an irrational angle, and vice versa.

Geometric Method for the Reduction of Systems
According to [15] on centrally-stable manifold, the system around an equilibrium point O takes the form ,  together with its first derivatives vanish at the origin. In this case, the right side of the system may depend on the control u either continuously (in this case, smooth manifolds, discussed below, depend continuously on u), or smooth. In the latter case, u is included in the number of "central" variables x and, thus, further investigated the variety and the foliation will have the smoothness of u, equal to the smoothness of x. Theorem 1. In a small neighborhood of the equilibrium state there exists a (m + k)-dimensional center-stable invariant manifold. Central-stable manifold is not uniquely defined, but for any two manifolds 1 Proof. Without loss of generality, we consider the second-order system. Let the system defined on a manifold gives rise to a one-parameter transformation group generated by the operator We need to show that the operator system defined by 2 sC W an equivalent ( -group option).
We write the transformation of the group as a Li series (the operator exponential): A similar formula holds for the second coordinate, which together form a differential equation: Solution of the resulting Cauchy problem can be obtained by expanding the operator  , in a Taylor series in powers of : thus, we have: Finally, a number will look like in the Hausdorff form: By Theorem from [15] function 1 sC  and 2 sC  determine the same expansion in Taylor series, which suggests field commutations: Thus, the theorem is proved. At the reference time t  -t of A, B and C go into respectively A, B and C. Thus, part of the spectrum of characteristic exponents corresponding to the variable z, is now left of the imaginary axis, and part of the spectrum, corresponding to the variable y-the right of it. To the system obtained from system (1) by time reversal, we can again apply the theorem about the Central stable manifold and get the following theorem on center-unstable manifold. Straightening of the central stable and center-unstable manifolds, as well as rectification of the strong stable and strong unstable invariant foliations on these manifolds leads to the good fact.
Theorem 3 [15]. With a -smooth transformation of the system (1) can locally be reduced to Here, the local center-unstable manifold given by the equation {y = 0}, the local center-stable manifold-{z = 0}, and a local center manifold-{y = 0, z = 0}. Strong stable foliation consists of surfaces {x = const, z = 0}, and the layers are strongly unstable foliation of the form {x = const, y = 0}.
A similar theory is constructed for discrete systems. For systems that admit the group of symmetries and construct models of systems, reduced to a central invariant manifold results can be summarized as follows [17].
In the local region of the qualitative dynamic behavior of systems typologically equivalent system, reduced to the central manifold where whose structure is determined by the symmetry transformation. For a discrete system, the local center manifold is determined by the system:  that are close to periodic, and the construction of finite-transformations transform one area to another. That is, the construction of the symmetry group of phase trajectories, which is characterized by the transformation of graphs: The resulting transformations determine the structure of the desired evolution equations.

Algorithms and Applications
The result of reconstruction algorithms [1] is a set of points belonging to the attractor of the system. The task is an automatic search for symmetries of the local sections of the phase trajectories. We seek the symmetry of translation, rotation, stretching and compression.
You must select the initial sequence collection sites so that when reducing them to a single scale, position and angle of rotation, they would be most similar to each other, as well as numerical indicators of transformations taking one piece to another without breaking of symmetry, to give a numerical estimate degree of symmetry breaking. The Figure 1 shows a schematic illustration of the problem.
The task of allocating a set of similar fragments in the circuit is defined as follows: Input: sequence of n points 1 2 , , , n v v v  , the function   1 2 , Dist V V , which gives an assessment of similarity between any two fragments of the original sequence 1 , , , Output: a set of disjoint fragments of the original sequence Number of different fragments, which can be identified in the original sequence: Number of arbitrary sets of fragments, determined by the size of the set of all subsets of the set 1 2 , , , m V V V  : Marking. In the simplest case for the selection of fragments of trajectories satisfying the decision to consider the fragments located between all possible pairs of contour points. In most problems, this approach is not possible because of limited computer resources. So it makes sense to consider as a start / end fragment only some points of the contour, for example, extremes of functions of each of the coordinates and evenly distributed along the length of the markers on the slowly varying parts (Figure 2). This approach provides a sufficient rational partition of the circuit, although it suffers from some redundancy.
Partitioning algorithm can be anything important to exclude monotonous plots and identify the markers characteristic sections of the circuit. Extracting fragments. To isolate fragments of contour to study the simplest case, you can use brute force pairs of markers. In this case, the number of fragments needed to be addressed:  where n -number of markers on the contour. It makes sense to limit the length of a fragment in the original contour points as the top and bottom as well as fragments, the length of the approaching to the original contour, are not of interest to identify patterns of distribution symmetry, and too small fragments do not reflect the specific behavior of the system. In addition, the introduction of restrictions on the length of the aboveincrease the number of fragments for consideration is linear, as opposed to quadratic, with the full brute force.
Selected fragments for the possibility of a comparison between a need to interpolate in a way that they all described the same fixed number of points (Figure 3).
Normalization. After reducing to one the number of points each of the fragments should be subjected to the normalization procedure. The aim of the procedure is to convert the fragment into a descriptor-the image that is invariant with respect to the transport, rotation, and scaling the original fragment, and obtaining the numerical parameters of this transformation (Figure 4).
The proposed normalization procedure is a sequence of the following: Isolation in a fragment of the so-called "axis"-pairs of points, the distance between the maximum in the current space projection. Rotation of the fragment so that its axis were parallel to the coordinate axes.
The construction restricting fragment of the n-dimensional parallelepiped and the shift of the fragment so that the center of the parallelepiped coincides with the center coordinates.
Scaling the fragment so that the length of its main axis is equal to 1.
Fragment of an n-dimensional the contour F consisting of m real points, represented as a matrix : , , , For a compact description of the transformation fragment as a single matrix, we give it to homogeneous coordinates using matrix transformations, obtain: where rt g -angle, a significant elements of the matrix are located in cells (i, i), (i, j), (j, i) и (j, j).

Normalization Algorithm Is as Follows
3) 13) Transfer fragment so that the center of the bounding parallelepiped it coincided with the origin. Conversion is performed by multiplying with the translation matrix sh M form (6) Figure 9.
The presence of the symmetry of the system determined the structure of the form (3). Parameter identification system using the m llowing result:     sion can be close enough to the identity for small perturbations.

Conclusion
The use of symmetric properties of evolution equations can be successfully used for the reconstruction of the system. Developing a constructive method for structural identification systems with chaotic dynamics based on geometric analysis of the experimentally obtained phase portrait of the system.