^{1}

^{*}

^{2}

The study of the parameter space of chaotic systems is complicated by its high dimensionality (multi-parametricability). Two approaches to the study of chaotic systems are presented: multi-parameter analysis and optimal suppression of chaotic dynamics. For non-autonomous chaotic systems, this is the way to compare the effectiveness of various correction parameters that provide optimal removal of irregular dynamics. For the class of autonomous chaotic systems, this is the way to investigate the optimal conditions of super-stable behavior for the chaotic system.

A rather wide class of dynamic systems that demonstrate chaotic behavior is mathematically described by the systems of non-linear autonomous and non-autonomous differential equations [1-3]. As a rule, the dynamics of such systems is defined by the values of several parameters simultaneously.

A clear insight into the mechanisms of chaotic dynamics appearance (disappearance) is possible only through multi-parameter analysis of the system. It is shown in [

In this paper, we show two approaches to multi-parameter analysis of non-autonomous (dissipative nonlinear oscillators) and autonomous (Lorenz-like) chaotic systems. The paper is divided into two parts. In Part I we show, that Melnikov method [1,6,7] can be applied for two-parameter analysis of chaos suppression conditions. The connections between the key parameters of the system, revealed in the process, let us make an advance in the investigation of optimal “chaos-regular dynamics” transitions. The solution for the problem of multi-parameter optimization which allows considering the possibility of achieving the state of super-stability by chaotic system is presented in Section 2.

In this part, we present the results of a two-parametrical analysis of optimal chaotic dynamics suppression in dissipative nonlinear oscillators. As a reference model of this class of chaotic systems we chose a double-well Duffing oscillator

Among other non-autonomous chaos models this oscillator is knows as a paradigmatic one, as it clearly demonstrates principal characteristics of a wide class of nonlinear dissipative oscillators.

The parametric space of the oscillator (1) is five-dimensional and unites two groups of parameters that play different roles in dynamics. It contains internal parameters and parameters that characterize the environment of the system. The multidimentionality of parametric space results in multi-parametricability of Melnikov function. Thus the function calculation for (1) gives

where

.

If is a set of possible parameter values for oscillator (1) then, using the condition of sign change for Melnikov function, we can allocate in P: 1) regular dynamics region, that corresponds to the values of parameters for which A < C; 2) complex dynamics region, where ; 3) the boundary (extreme values of parameters) between the regions of regular and complex dynamics (RC boundary), which is defined by the condition. As a starting point, in this paper we use the configuration of parameters,

when oscillator (1) demonstrates typical chaotic behavior.

Convenient way for visual study of the structure of oscillator parametric space (1) is two-parameter analysis of the boundary of regular and complex dynamics regions.

Let’s modify the equality first as

and then as

These dependencies give us two points of view for the study of the system. For fixed values of parameters from the group function (2) defines the RCboundary on the plane of internal parameters of the oscillator. Function (3) allows studying the RCboundary at fixed values and on the plane of the parameters of chaotizing perturbation.

The RC-boundary built with the use of the functions (2) and (3) is shown in _{0}) on the plane of parameters. These versions correspond to the changes of frequency values, (

The results shown in

and

correspondingly, while the region of complex dynamics is specified by the conditions and. In this case variations of parameter values lead to various shifts of RC-boundary and change the position of minimum.

Consider the situation when initial values of parameters and are corrected. Having applied corrective amendments and we will write down the oscillator (1) equation as

It seems interesting to study the combined influence of amendments and that can provide chaos suppression. As minimal corrective amendments, , we can use the values that lie on the RC-boundary. With them the solution of Equation (4) will be stable instead of chaotic.

Let’s define the set of accessible values and as

where

the Melnikov function of oscillator (4). Then, by formalizing the demand of small corrective amendments, necessary for the optimal chaos suppression, we get the optimal correction problem:

To solve the problem (5) we use statement (2) instead of repeating the process of calculating Melnikov function M_{1}. By changing, we transform (2) as follows

As a dependence expression (6) demonstrates the relationship between the values of corrective amendments and, that lie on the RC-boundary.

Actually, each point on the curve can be a solution (5), but only one of them is the best. For example, if we correct parameters and individually, we get the pairs: and, correspondingly. It means that if parameter is being corrected, the system achieves stability with smaller amendment for this criterion (5) than for. To find optimal pair we should find the minimum for the function while moving along the RC-boundary.

The changes of quality criterion along the RC-boundary are also shown in

. Thus the solution of optimization problem (5) gives us interrelated values of amendments and which provide optimal transi-

tion from chaotic to stable region of parametrical space for the given configuration of parameters.

A rather wide class of chaotic Lorenz-like systems is described as follows:

where—the matrix of system parameters,—nonlinear part that characterizes nonlinear perturbation.

While suppressing chaotic dynamics, usually, we set a goal and try to find out under what conditions can the system (7) achieve stable behavior. However there exists a narrow class of systems which demonstrate a special type of regular behavior. They are super-stable systems (for essential information on super-stability see [

It is essential that, unlike stability, the super-stability is preserved at nonlinear perturbations. It grounds the possibility of analysis of super-stability achievement conditions with regard to system (7).

System (7) can have several states of equilibrium

, where—the number of states of equilibrium. Then Jacobian matrix for k-state of equilibrium of the system (7) will be written as

where—is the matrix with the elements. The generalization of super-stability condition for the systems with multiple states of equilibrium (7) leads to the condition

where

.

The value is called the super-stability degree. For every k-state of equilibrium the conditions are linear restrictions on J_{k} matrix coefficients that show that negative diagonal dominance should necessarily be present.

Let initial values of parameters of the system (7) lie in the area of parametric space that corresponds to chaotic dynamics. Then according to the conditions of dynamics chaotization the Jacobian matrix J_{k} of the system (7) will be unstable for every k and condition (8) will not be performed.

Consider the system stabilization problem in the following form: for the given unstable matrix should be found the closest super-stable matrix. Let’s define the closeness as the least value, where—

some matrix norm (from here on).

Then the task will be reduced to finding the optimal corrective matrix that provides super-stability of the matrix with condition

where C—the set of corrective matrixes, which allow performing conditions (8) for the coefficients of the corrected matrix.

The problem formulated above can be regarded as a generalization of the problem of static output stabilization of state space MIMO systems. The peculiarity of the problem is that we must provide not just stability, but super-stability. At the same time optimization criterion (9) is built with regard to the possibility of the system having several states of equilibrium. Hence the solution of the problem is the pair: the matrix and state of equilibrium correspondent to it. Together they define the corrected matrix which is the super-stable matrix closest to.

It may seem that in problem (9) the correction rule of the matrix should be presented as transformation of its coefficients

where—corrective amendments which define matrixes H_{k} and. However, as super-stability conditions (8) are strict there arises a question. Is there always exist the corrective matrix H_{k} for the system (7), which allows fulfilling super-stability conditions?

The work [_{k} coefficients allows finding the best way to change the parameters and modify a chaotic system into the super-stable one. The structural analysis of the matrix substantionally restricts the correction rule: 1) the matrix A is often sparse (that is some are equal to zero) and its zero coefficients should not be influenced by; 2) the matrix G_{k} can not be corrected at all. These restriction are caused by our attempts to apply multiparametrical correction to the system while preserving its properties. If these restrictions are not fulfilled the system undergo structural changes (new terms appear in the equation). With regard to these restrictions the correction rule will be written in the form

where H_{k}—corrective matrix structurally equivalent to A.

The rule (10) defines only the allowed method of correction, but does not provide us with any information about super-stability accessibility. Let’s write super-stability conditions of the corrected system:

Then super-stability accessibility means that there exists matrix H_{k} (according to (10)) with coefficients that satisfy these equations. It is easy to name the situations when super-stability is initially impossible. It happens when one of diagonal elements of the matrix J_{k} is equal to zero. Then the existing condition (11) results in an unavoidable contradiction.

If super-stability is accessible (i.e. set C in (9) is not empty) the problem (9) can be solved in two steps. At step I we should find the matrixes that is to solve s quadratic programming problem

with restrictions (11). At step II we find with matrix correspondent to it. To do this it is necessary to solve the problem. Once these steps are performed we get the optimal method of parametrical correction which provides super-stable dynamic behavior of the system.

Thus super-stability is a unique property peculiar to few chaotic systems. It can serve as a foundation for an additional classification (super-stabilizable systems). For example the optimal multiparametrical correction problem (9) can easily be solved for the subclass of chaotic systems investigated in [

where. If, and we get Lorenz, Chen and Lsystems, correspondingly. The existence of pairs, for these systems means that the character of their dynamics can be transformed from chaotic to the stable one through parametric correction. So the systems are super-stabilizable.

Two results presented in the paper illustrate the effectiveness of implementation of multi-parametrical analysis for the peculiarities of chaos dynamics suppression.

First of all it is shown that Melnikov criterion can be used for stating and solving important optimization problems which deepen our understanding of chaos suppression process. While searching the optimal method of turning the system into the non-chaotic one we reveal the connection between corrective amendments to the parameters and find their minimal values.

The second result was achieved during the study of two opposite states of autonomous dynamic systems— chaos and super-stability. While controlling the chaotic system it is normal to study the conditions of bringing the system to stable behavior. Still the question whether it is possible for the system to achieve super-stability remains unanswered. The conditions of achievement of the super-stable dynamics were studied for a class of Lorenz-like chaotic systems. Based on super-stability criteria the optimal multi-parametrical (matrix) correction problem was formulated and an effective method of solution was presented. The offered method allows finding super-stabilizable chaotic systems and transforming the chaotic dynamics to the super-stable one using optimal parametric correction.

Thus the named methods of multi-parametrical analysis are applied to a wide range of chaotic systems. They help compare the efficiency of different forms of parametrical perturbations and choose the ones which provide optimal suppression of chaotic dynamics.