# Chaos in Essentially Singular 3D Dynamical Systems with Two Quadratic Nonlinearities

Authors: Vasiliy Belozyorov
Journal: Journal of Nonlinear Dynamics
 (2016)

---

## Abstract

A new class of 3D autonomous quadratic systems, the dynamics of which demonstrate a chaotic behavior, is found. This class is a generalization of the well-known class of Lorenz-like systems. The existence conditions of limit cycles in systems of the mentioned class are found. In addition, it is shown that, with the change of the appropriate parameters of systems of the indicated class, chaotic attractors different from the Lorenz attractor can be generated (these attractors are the result of the cascade of limit cycles bifurcations). Examples are given.


---

## Body

## 1. Introduction

Chaos is a very interesting nonlinear phenomenon, which has been intensively studied in the last four decades. Many potential applications have come true in secure communication, laser and biological systems, hydrodynamics and chemical technologies, and other areas (see, e.g., [1, 2] and many references cited therein).Notice that interest in the study of the chaotic systems increases only (see, e.g., [3–19]). We will specify a few scientific directions, in which application of methods of chaotic dynamics appeared most effective.In the scientific literature, chaotic economic systems have got much attention because their complex dynamic behaviors cannot be described by traditional models. Recently, a few researches on the usage of these systems in cryptographic algorithms have been conducted. A new image encryption algorithm based on a chaotic economic map is proposed. Implementation of the proposed algorithm on a plain image based on the chaotic map is performed [3].The paper [4] investigates the multipulse heteroclinic bifurcations and chaotic dynamics of a laminated composite piezoelectric rectangular plate by using an extended Melnikov method in the resonant case.A chaotic system arising from double-diffusive convection in a fluid layer is investigated in [5] based on the theory of dynamical systems. A five-dimensional model of chaotic systems is obtained using the Galerkin truncated approximation. The results showed that the transition from steady convection to chaos via the Hopf bifurcation produced a limit cycle, which may be associated with a homoclinic explosion at a slightly subcritical value of the Rayleigh number.The paper [6] investigates the problem of two kinds of function projective synchronization of financial chaotic system with definite integration scaling function, which are not fully considered in the existing researches. In this paper, the following two questions (unlike previous methods) are investigated: (1) two kinds of the definite integration scaling function projective synchronization are given and (2) the upper and lower limit of the definite integral scaling function are the bound dynamical systems.In this connection, the important problem is the establishment of theoretical results, where it is possible to extend the wide classes of the dynamical systems. One of such classes is the class of 3D quadratic systems.Vast scientific literature is devoted to researches of different types of systems of this class. For instance, 3-scroll and 4-scroll chaotic attractors were investigated in [7]. New existence conditions of heteroclinic and homoclinic orbits were found in [8–10]. New three-wing and four-wing chaotic attractors were opened in [11–14]. The article [15] summarizes investigations of quadratic dynamical systems, which are indicated in [7–14].Another interesting domain of nonlinear dynamics is the control by the chaotic systems. Novel controllers for stabilization and tracking of chaotic and hyperchaotic Lorenz systems using extended backstepping techniques were proposed in [16]. Based on the proposed approach, the generalized weighted controllers were designed to control the chaotic behavior as well as achieve synchronization in chaotic and hyperchaotic Lorenz systems [16].Above, we mentioned researches of the class of 3D quadratic dynamical systems. The well-known Lorenz system, which is described by the equations(1)
x
˙
t
=
σ
y
t
-
x
t
,
y
˙
t
=
r
x
t
-
y
t
-
x
t
z
t
,
z
˙
t
=
-
b
z
t
+
x
t
y
t
,is a basic representative of this class (here, σ
=
10, b
=
8
/
3, and r
>
24).System (1) is a three-dimensional autonomous system with only two quadratic terms in its nonlinearity, which is very simple in the algebraic structure and yet is fairly complex in dynamical behaviors. The observation of these two seemingly contradictory aspects of the Lorenz system thereby triggered a great deal of interest from the scientific community to seek closely related Lorenz-like systems, by different motivations and from various perspectives [15].It should be said (as it was noted in [15]) that almost all known chaotic attractors of 3D quadratic dynamical systems can be generated (at the appropriate values of parameters) by the following system:(2)
x
˙
t
=
a
11
x
t
+
a
12
y
t
,
y
˙
t
=
a
21
x
t
+
a
22
y
t
+
m
1
x
t
z
t
+
m
2
y
t
z
t
,
z
˙
t
=
a
33
z
t
+
m
3
x
t
y
t
+
m
4
x
2
t
+
m
5
y
2
t
+
m
6
z
2
t
+
c
,where a
11
,
…
,
a
33, m
1
,
…
,
m
6, and c are real parameters.The primary purpose, which was pursued by the research program that was begun in this work, is to set communication between algebraic and geometric properties of solutions of system (2).In spite of the vast researches, which are developed in connection with the dynamics of the Lorenz system, a lot of the unsolved problems yet remain. A strange attractor corresponds to the complex nonperiodic mode of the Lorentz system. The system also generates spatial limit cycles in three-dimensional space. A few papers are devoted to the research of the complex closed trajectories which are also connected to chaos, because a transition to the chaotic behavior is sometimes realized in consequence of bifurcations of periodic trajectories on saddle trajectories. Recently (see [17–19]), a row of results describing the behavior of the Lorenz-like systems appeared. The present article is one of such results.In conclusion, we explain the name of the article. At first sight, it seems that such systems are simpler than 3D systems in which all three equations are nonlinear. However, the existence of one linear equation in system (2) makes application of the theory of positive (or negative) definite quadratic forms impossible. This definiteness was widely used in [17, 18].In the present paper, some special cases of system (2) will be considered. Note that in system (1) two nonlinearities exist: there are x
·
y and x
·
z. In the system which will be studied below, we also will have two nonlinearities: x
·
(
d
y
y
+
d
z
z
) and y
·
(
d
y
y
+
d
z
z
), where d
y and d
z are real numbers such that d
y
2
+
d
z
2
≢
0. In our opinion, this is some generalization of system (1). In addition, in this work, the researches begun in [17, 18] will be continued.


## 2. Bounded Solutions of Quadratic Dynamical Systems

Consider the following 3D autonomous quadratic system:(3)
x
˙
t
=
H
x
+
f
x
,where x
=
(
x
,
y
,
z
)
T; H
=
{
h
i
j
}, i
,
j
=
1
,
…
,
3, is a real (
3
×
3
)-matrix; and (4)
f
x
=
f
1
x
,
y
,
z
,
f
2
x
,
y
,
z
,
f
3
x
,
y
,
z
T
∈
R
3
,
f
1
x
,
y
,
z
=
a
12
x
y
+
a
22
y
2
+
a
13
x
z
+
a
23
y
z
+
a
33
z
2
,
f
2
x
,
y
,
z
=
b
12
x
y
+
b
22
y
2
+
b
13
x
z
+
b
23
y
z
+
b
33
z
2
,
f
3
x
,
y
,
z
=
c
12
x
y
+
c
22
y
2
+
c
13
x
z
+
c
23
y
z
+
c
33
z
2are real quadratic polynomials.Suppose that the matrix(5)
a
12
a
13
b
12
b
13
c
12
c
13has rank 1 or 2. Then, by suitable linear transformations of variables x
→
x
1
+
α
1
y
1
+
α
2
z
1 (α
1
,
α
2
∈
R), y
→
y
1, and z
→
z
1, system (3) can be represented in the same form as (3), where H
→
H
¯
=
{
h
i
j
¯
}, i
,
j
=
1
,
…
,
3, and (6)
f
¯
1
x
1
,
y
1
,
z
1
=
a
¯
22
y
1
2
+
a
¯
23
y
1
z
1
+
a
¯
33
z
1
2
,
f
¯
2
x
1
,
y
1
,
z
1
=
b
¯
12
x
1
y
1
+
b
¯
22
y
1
2
+
b
¯
13
x
1
z
1
+
b
¯
23
y
1
z
1
+
b
¯
33
z
1
2
,
f
¯
3
x
1
,
y
1
,
z
1
=
c
¯
12
x
1
y
1
+
c
¯
22
y
1
2
+
c
¯
13
x
1
z
1
+
c
¯
23
y
1
z
1
+
c
¯
33
z
1
2
,and b
¯
12
≠
0 or b
¯
13
≠
0. (Note that if the rank of matrix (5) is equal to 1, then the case c
¯
12
=
c
¯
13
=
0 is not excepted.)Thus, it can be considered that, for system (3), the conditions(7)
a
12
=
a
13
=
0
,
b
12
∨
b
13
≠
0are fulfilled.Introduce into system (3) (taking into account (7)) new variables ρ and ϕ under the formulas y
=
ρ
cos
⁡
ϕ and z
=
ρ
s
i
n
ϕ, where ρ
>
0. Then, after replacement of variables and multiplication of the second and third equations of system (3) on the matrix (8)
cos
⁡
ϕ
t
s
i
n
ϕ
t
-
s
i
n
ϕ
t
ρ
t
cos
⁡
ϕ
t
ρ
t
, we get(9)
x
˙
t
=
h
11
x
t
+
h
12
cos
⁡
ϕ
t
+
h
13
s
i
n
ϕ
t
ρ
t
+
a
22
cos
2
⁡
ϕ
t
+
a
23
cos
⁡
ϕ
t
s
i
n
ϕ
t
+
a
33
s
i
n
2
ϕ
t
ρ
2
t
,
ρ
˙
t
=
h
21
c
o
s
ϕ
t
+
h
31
s
i
n
ϕ
t
x
t
+
h
22
c
o
s
2
ϕ
t
+
h
33
s
i
n
2
ϕ
t
+
h
32
+
h
23
c
o
s
ϕ
t
s
i
n
ϕ
t
ρ
t
+
b
12
c
o
s
2
ϕ
t
+
b
13
+
c
12
c
o
s
ϕ
t
s
i
n
ϕ
t
+
c
13
s
i
n
2
ϕ
t
x
t
ρ
t
+
b
22
c
o
s
3
ϕ
t
+
b
23
+
c
22
c
o
s
2
ϕ
t
s
i
n
ϕ
t
+
b
33
+
c
23
c
o
s
ϕ
t
s
i
n
2
ϕ
t
+
c
33
s
i
n
3
ϕ
t
ρ
2
t
,
ϕ
˙
t
=
-
h
21
s
i
n
ϕ
t
+
h
31
c
o
s
ϕ
t
x
t
ρ
t
+
h
32
c
o
s
2
ϕ
t
-
h
23
s
i
n
2
ϕ
t
+
h
33
-
h
22
c
o
s
ϕ
t
s
i
n
ϕ
t
-
b
13
s
i
n
2
ϕ
t
+
b
12
-
c
13
s
i
n
ϕ
t
c
o
s
ϕ
t
-
c
12
c
o
s
2
ϕ
t
x
t
-
-
c
22
c
o
s
3
ϕ
t
+
b
22
-
c
23
c
o
s
2
ϕ
t
s
i
n
ϕ
t
+
b
23
-
c
33
c
o
s
ϕ
t
s
i
n
2
ϕ
t
+
b
33
s
i
n
3
ϕ
t
ρ
t
.Consider the first and second equations of system (9):(10)
x
˙
t
=
h
11
x
+
f
cos
⁡
ϕ
,
s
i
n
ϕ
ρ
+
f
22
cos
⁡
ϕ
,
s
i
n
ϕ
ρ
2
,
ρ
˙
t
=
g
cos
⁡
ϕ
,
s
i
n
ϕ
x
+
h
cos
⁡
ϕ
,
s
i
n
ϕ
ρ
+
g
12
cos
⁡
ϕ
,
s
i
n
ϕ
x
ρ
+
g
22
cos
⁡
ϕ
,
s
i
n
ϕ
ρ
2
,where ϕ is a real parameter; and (11)
f
cos
⁡
ϕ
,
s
i
n
ϕ
=
h
12
cos
⁡
ϕ
+
h
13
s
i
n
ϕ
,
f
22
cos
⁡
ϕ
,
s
i
n
ϕ
=
a
22
cos
2
⁡
ϕ
+
a
23
cos
⁡
ϕ
s
i
n
ϕ
+
a
33
s
i
n
2
ϕ
,
g
cos
⁡
ϕ
,
s
i
n
ϕ
=
h
21
cos
⁡
ϕ
+
h
31
s
i
n
ϕ
,
h
cos
⁡
ϕ
,
s
i
n
ϕ
=
h
22
cos
2
⁡
ϕ
+
h
33
s
i
n
2
ϕ
+
h
23
+
h
32
cos
⁡
ϕ
s
i
n
ϕ
,
g
12
cos
⁡
ϕ
,
s
i
n
ϕ
=
b
12
cos
2
⁡
ϕ
+
b
13
+
c
12
cos
⁡
ϕ
s
i
n
ϕ
+
c
13
s
i
n
2
ϕ
,
g
22
cos
⁡
ϕ
,
s
i
n
ϕ
=
b
22
cos
3
⁡
ϕ
+
b
23
+
c
22
cos
2
⁡
ϕ
s
i
n
ϕ
+
b
33
+
c
23
cos
⁡
ϕ
s
i
n
2
ϕ
+
c
33
s
i
n
3
ϕ
.In [17, 18], conditions of appearance of chaos in system (3) were analyzed. As a result of this analysis, the following theorem (with small corrections by comparison to [16, 17]) was proved.LetΔ
1
(
cos
⁡
ϕ
,
s
i
n
ϕ
)
≡
h
11
·
h
(
cos
⁡
ϕ
,
s
i
n
ϕ
)
-
f
(
cos
⁡
ϕ
,
s
i
n
ϕ
)
·
g
(
cos
⁡
ϕ
,
s
i
n
ϕ
), Δ
2
(
cos
⁡
ϕ
,
s
i
n
ϕ
)
≡
f
22
(
cos
⁡
ϕ
,
s
i
n
ϕ
)
·
g
12
(
cos
⁡
ϕ
,
s
i
n
ϕ
), and Δ
3
(
cos
⁡
ϕ
,
s
i
n
ϕ
)
≡
g
22
2
(
cos
⁡
ϕ
,
s
i
n
ϕ
)
+
4
Δ
2
(
cos
⁡
ϕ
,
s
i
n
ϕ
) be the bounded functions.Theorem 1 (see [17, 18]).
Leth
11
<
0. Suppose also that ∀
ϕ
∈
R for system (10) the following conditions are fulfilled:(i)
Δ
2
(
u
,
v
) is a nonsingular form of the variables u
≡
cos
⁡
ϕ and v
≡
s
i
n
ϕ.
(ii)
EitherΔ
1
(
cos
⁡
ϕ
,
s
i
n
ϕ
) is a negative constant or Δ
1
(
cos
⁡
ϕ
,
s
i
n
ϕ
) is a periodic nonpositive function or Δ
1
(
cos
⁡
ϕ
,
s
i
n
ϕ
) is a periodic alternating in sign on (
-
∞
,
∞
) function.
(iii)
Δ
2
(
cos
⁡
ϕ
,
s
i
n
ϕ
) and Δ
3
(
cos
⁡
ϕ
,
s
i
n
ϕ
) are periodic nonpositive functions.
Suppose that the condition(12)
l
i
m
i
n
f
t
→
∞
ρ
t
=
0is also valid. (From this condition, it follows that ∀
ϵ
>
0 there exists a numerical sequence t
n
→
∞ as n
→
∞ such that, ∀
n
ρ
(
t
n
)
<
ϵ.)
Then, in system (9) (and system (3)), there are chaotic dynamics.


## 3. Definition of Essentially Singular 3D Dynamical Systems

The general definition of the singular system was given in [19]. In the present paper, we are restricted by a narrower definition.Consider the following 2D autonomous quadratic system:(13)
x
˙
t
=
a
11
x
t
+
a
12
y
t
+
b
11
x
2
t
+
2
b
12
x
t
y
t
+
b
22
y
2
t
,
y
˙
t
=
a
21
x
t
+
a
22
y
t
+
c
11
x
2
t
+
2
c
12
x
t
y
t
+
c
22
y
2
t
,where a
11
,
…
,
a
22, b
11
,
…
,
b
22, and c
11
,
…
,
c
22 are real numbers.Introduce the following real(
2
×
2
)-matrices:(14)
A
=
a
11
a
12
a
21
a
22
,
T
1
=
b
11
b
12
c
11
c
12
,
T
2
=
b
12
b
22
c
12
c
22
.Replace variablesx, y in system (13) by new variables x
1, y
1 under the formula (15)
x
y
⟶
S
·
x
1
y
1
,where S is a linear transformation from the group G
L
(
2
,
R
) of all linear inverse transformations of the space R
2 [19]. In this case, the triple of matrices (
A
,
T
1
,
T
2
) transforms into triple S
∘
(
A
,
T
1
,
T
2
)
=
(
S
-
1
A
S
,
(
S
-
1
T
1
,
S
-
1
T
2
)
·
(
S
⊗
S
)
).Remember that a scalar polynomialf
(
A
,
T
1
,
T
2
) is called an invariant of weight l of the group G
L
(
2
,
R
), if ∀
S
∈
G
L
(
2
,
R
) and ∀
A
,
T
1
,
T
2
;
we have
f
(
S
∘
(
A
,
T
1
,
T
2
)
)
=
(
d
e
t
S
)
l
×
f
(
A
,
T
1
,
T
2
)
, where l
≥
0 is some integer.With the help of matricesT
1, T
2, we construct the auxiliary not depending on A invariants of weight 2 [19]: (16)
I
1
=
d
e
t
t
r
T
1
,
t
r
T
2
·
T
1
t
r
T
1
,
t
r
T
2
·
T
2
,
J
2
=
d
e
t
T
1
T
2
-
T
2
T
1
,
K
3
=
d
e
t
t
r
T
1
,
t
r
T
2
t
r
T
1
,
t
r
T
2
·
T
1
T
2
-
T
2
T
1
,where t
r
Pis a trace of the square matrix P.Now, we can introduce the main invariant of the present paper(17)
D
=
I
1
+
27
J
2
-
5
K
3of weight 2 [19].In system (13), we replace the variable y by the variable ρ. Now, we can replace the coefficients of system (13) by appropriate coefficients of system (10). (In this case, the invariants J
2 and D will be the functions of the variable ϕ
:
D
(
cos
⁡
ϕ
,
s
i
n
ϕ
)
,
J
2
(
cos
⁡
ϕ
,
s
i
n
ϕ
).)Letu
≡
cos
⁡
ϕ and v
≡
s
i
n
ϕ.Definition 2 (see [19]).
System (3) is called essentially singular if for arbitrary real number ϕ the conditions D
(
u
,
v
)
≡
D
(
c
o
s
ϕ
,
s
i
n
ϕ
)
≤
0 and J
2
(
u
,
v
)
≡
J
2
(
c
o
s
ϕ
,
s
i
n
ϕ
)
≤
0 are fulfilled. In addition, there exist the numbers ϕ
i such that D
(
c
o
s
ϕ
i
,
s
i
n
ϕ
i
)
=
J
2
(
c
o
s
ϕ
i
,
s
i
n
ϕ
i
)
=
0, i
=
1,2
,
….We compute the invariantsJ
2
(
u
,
v
) and D
(
u
,
v
) for system (3) under conditions c
12
=
c
13
=
c
22
=
c
23
=
c
33
=
0. Then, we have (18)
J
2
u
,
v
=
u
2
b
12
u
+
b
13
v
2
×
-
a
12
u
+
a
13
v
b
22
u
3
+
b
23
u
2
v
+
b
33
u
v
2
+
u
b
12
u
+
b
13
v
a
22
u
2
+
a
23
u
v
+
a
33
v
2
,
D
u
,
v
=
u
2
b
12
u
+
b
13
v
2
×
a
12
u
+
a
13
v
-
b
22
u
3
-
b
23
u
2
v
-
b
33
u
v
2
2
+
4
u
b
12
u
+
b
13
v
a
22
u
2
+
a
23
u
v
+
a
33
v
2
=
u
2
b
12
u
+
b
13
v
2
a
12
u
+
a
13
v
+
b
22
u
3
+
b
23
u
2
v
+
b
33
u
v
2
2
+
4
J
2
u
,
v
.From here, it follows that if D
(
u
,
v
)
≤
0 (D
(
u
,
v
)
<
0), then J
2
(
u
,
v
)
≤
0 (J
2
(
u
,
v
)
<
0).


## 4. Lorenz-Like Systems

Consider the following system:(19)
x
˙
t
=
h
11
x
t
+
y
t
a
22
y
t
+
a
23
z
t
,
y
˙
t
=
h
21
x
t
+
h
22
y
t
+
h
23
z
t
-
x
t
k
a
22
y
t
+
k
a
23
z
t
,
z
˙
t
=
h
31
x
t
+
h
32
y
t
+
h
33
z
t
.Here, we will consider that h
11
<
0, k
>
0, h
32
≠
0, and the equilibrium (
0,0
,
0
)
T is a saddle point. In addition, in this section, we will consider that h
21
=
h
31
=
0
.The direct calculations of the invariantsD
(
cos
⁡
ϕ
,
s
i
n
ϕ
) and J
2
(
cos
⁡
ϕ
,
s
i
n
ϕ
) of system (19) show that J
2
(
cos
⁡
ϕ
,
s
i
n
ϕ
)
=
-
k
3
·
(
cos
⁡
ϕ
)
4
·
(
a
22
cos
⁡
ϕ
+
a
23
s
i
n
ϕ
)
4 and D
(
cos
⁡
ϕ
,
s
i
n
ϕ
)
=
4
J
2
(
cos
⁡
ϕ
,
s
i
n
ϕ
). Thus, we have D
(
cos
⁡
ϕ
i
,
s
i
n
ϕ
i
)
=
J
2
(
cos
⁡
ϕ
i
,
s
i
n
ϕ
i
)
=
0 at ϕ
i
=
π
/
2
+
2
π
·
i or ϕ
i
=
-
a
r
c
t
a
n
(
a
22
/
a
23
)
+
π
·
i, i
=
0,1
,
2
,
…, and it is clear that system (19) is essentially singular.System (19) is a special case of system (3). (If we put h
11
=
-
b
=
-
8
/
3, h
22
=
-
1
,
h
23
=
r
>
0
,
h
32
=
-
σ
=
-
10
,
h
33
=
σ
=
10, h
12
=
h
13
=
0, a
22
=
0, a
23
=
1, and k
=
1, then system (19) becomes the classic chaotic Lorenz system.) In coordinates y
=
ρ
cos
⁡
ϕ and z
=
ρ
s
i
n
ϕ, system (19) takes the following form (see (9)):(20)
x
˙
t
=
h
11
x
t
+
cos
⁡
ϕ
t
a
22
cos
⁡
ϕ
t
+
a
23
s
i
n
ϕ
t
ρ
2
,
ρ
˙
t
=
h
22
cos
2
⁡
ϕ
t
+
h
23
+
h
32
cos
⁡
ϕ
t
s
i
n
ϕ
t
+
h
33
s
i
n
2
ϕ
t
ρ
t
-
cos
⁡
ϕ
t
k
a
22
cos
⁡
ϕ
t
+
k
a
23
s
i
n
ϕ
t
x
t
ρ
t
,
ϕ
˙
t
=
h
32
cos
2
⁡
ϕ
t
+
h
33
-
h
22
cos
⁡
ϕ
t
s
i
n
ϕ
t
-
h
23
s
i
n
2
ϕ
t
+
s
i
n
ϕ
t
k
a
22
cos
⁡
ϕ
t
+
k
a
23
s
i
n
ϕ
t
x
t
.For this system, we haveΔ
2
(
cos
⁡
ϕ
,
s
i
n
ϕ
)
= -
-
J
2
(
cos
⁡
ϕ
,
s
i
n
ϕ
)
/
k
= -
k
(
cos
⁡
ϕ
·
(
a
22
cos
⁡
ϕ
+
a
23
s
i
n
ϕ
)
)
2
≤
0, Δ
3
(
cos
⁡
ϕ
,
s
i
n
ϕ
)
= 4
Δ
2
(
cos
⁡
ϕ
,
s
i
n
ϕ
)
= -
4
-
J
2
(
cos
⁡
ϕ
,
s
i
n
ϕ
)
/
k
= -
4
k
(
cos
⁡
ϕ
·
(
a
22
cos
⁡
ϕ
+
a
23
s
i
n
ϕ
)
)
2
≤
0. Thus, conditions (i)–(iii) of Theorem 1 are fulfilled.From the last equation of system (20), we have (21)
ϕ
t
=
ϕ
0
+
∫
0
∞
h
32
cos
2
⁡
ϕ
τ
+
h
33
-
h
22
cos
⁡
ϕ
τ
s
i
n
ϕ
τ
-
h
23
s
i
n
2
ϕ
τ
d
τ
+
∫
0
∞
s
i
n
ϕ
τ
k
a
22
cos
⁡
ϕ
τ
+
k
a
23
s
i
n
ϕ
τ
x
τ
d
τ
.Theorem 3.
Let the equilibrium point(
0,0
,
0
) of system (19) be saddle. Suppose also that for this system the following conditions are valid:(i)
h
11
<
0, h
21
=
h
31
=
0.
(ii)
k
·
(
a
22
2
+
a
23
2
)
>
0.
(iii)
There exists a finite limit(22)
l
i
m
t
→
∞
ϕ
t
t
=
p
≠
0
.
(iv)
(
3
a
22
2
+
a
23
2
)
1
+
(
p
/
h
11
)
2
>
(
a
22
2
+
a
23
2
+
2
a
22
)
2.
Then, all solutions of system (19) are bounded.Proof.
Construct the two following quadratic forms:(23)
f
u
,
v
=
h
22
u
2
+
h
23
+
h
32
u
v
+
h
33
v
2
,
g
u
,
v
=
h
32
u
2
+
h
33
-
h
22
u
v
-
h
23
v
2
.
Introduce the designations(24)
H
=
h
22
h
23
h
32
h
33
,
F
=
h
22
h
23
+
h
32
2
h
23
+
h
32
2
h
33
,
G
=
h
32
h
33
-
h
22
2
h
33
-
h
22
2
-
h
23
.(Here, the symmetric matrices F and G are matrices of the forms f
(
u
,
v
) and g
(
u
,
v
).)
Compute the eigenvalues of the matricesF and G. By λ
m
i
n
(
F
) and λ
m
a
x
(
H
) (λ
m
i
n
(
G
) and λ
m
a
x
(
G
)), denote minimal and maximal eigenvalues of the matrix F (the matrix G). Then, we have (25)
λ
m
a
x
,
m
i
n
F
o
r
λ
m
i
n
,
m
a
x
G
=
h
22
+
h
33
2
±
h
22
2
+
h
33
2
+
h
23
2
+
h
32
2
-
2
d
e
t
H
2
,
λ
m
a
x
,
m
i
n
G
o
r
λ
m
i
n
,
m
a
x
G
=
h
32
-
h
23
2
±
h
22
2
+
h
33
2
+
h
23
2
+
h
32
2
-
2
d
e
t
H
2
.
Since the point(
0,0
,
0
) is saddle, then d
e
t
H
≠
0. Besides, in order for condition (22) to be satisfied, it is necessary that the form g
(
u
,
v
) was negative definite. In this case, we can consider the following two possibilities:(26)
0
<
λ
m
i
n
F
<
λ
m
a
x
F
,
λ
m
i
n
G
<
λ
m
a
x
G
<
0
;
(27)
λ
m
i
n
F
<
0
<
λ
m
a
x
F
,
λ
m
i
n
G
<
λ
m
a
x
G
<
0
.
For research of system (20), we use the known in theory of differential equations Comparison Principle (see [17, 18]). Then, with the help of inequalities (26) and (27), we get instead of system (20) the following system:(28)
x
˙
t
=
h
11
x
t
+
cos
⁡
ϕ
t
a
22
cos
⁡
ϕ
t
+
a
23
s
i
n
ϕ
t
ρ
2
,
ρ
˙
t
=
λ
m
a
x
F
ρ
t
-
cos
⁡
ϕ
t
k
a
22
cos
⁡
ϕ
t
+
k
a
23
s
i
n
ϕ
t
x
t
ρ
t
,
ϕ
˙
t
=
λ
m
a
x
G
+
s
i
n
ϕ
t
k
a
22
cos
⁡
ϕ
t
+
k
a
23
s
i
n
ϕ
t
x
t
.
By virtue of condition (22) at t
→
∞, system (28) can be transformed to the following system with periodic coefficients:(29)
x
˙
t
=
h
11
x
t
+
cos
⁡
p
t
a
22
cos
⁡
p
t
+
a
23
s
i
n
p
t
ρ
2
t
,
ρ
˙
t
=
λ
m
a
x
F
ρ
t
-
cos
⁡
p
t
k
a
22
cos
⁡
p
t
+
k
a
23
s
i
n
p
t
x
t
ρ
t
.
Now, we will consider that for some values of parameters system (29) has a periodic solution. (This means that systems (19) and (20) have also the periodic solutions.)
Suppose that the numbert
m
≥
0 satisfies the conditions p
t
m
=
p
t
0
+
m
π and the number c
o
s
(
p
t
m
)
(
a
22
c
o
s
(
p
t
m
)
+
a
23
s
i
n
(
p
t
m
)
)
≠
0, where m
=
0,1
,
2
,
….
Further, we will consider the solutions of system (29) in a small neighborhood O
m of the point t
m: t
∈
O
m, m
=
0,1
,
2
,
…. With initial conditions x
m
0, ρ
m
0 for system (29), the solutions of system (28) in the point t
m are appointed.
Suppose that the timet
0 also satisfies the condition (30)
x
˙
t
0
=
h
11
x
t
0
+
cos
⁡
p
t
0
a
22
cos
⁡
p
t
0
+
a
23
s
i
n
p
t
0
ρ
2
t
0
=
0
. From here, it follows that (31)
x
˙
t
m
=
h
11
x
t
m
+
cos
⁡
p
t
m
a
22
cos
⁡
p
t
m
+
a
23
s
i
n
p
t
m
ρ
2
t
m
=
0
,
(32)
x
t
m
=
-
cos
⁡
p
t
0
a
22
cos
⁡
p
t
0
+
a
23
s
i
n
p
t
0
ρ
2
t
m
h
11
.
If we will take into consideration formula (32), the second equation of system (29) can be rewritten as(33)
ρ
˙
t
=
λ
m
a
x
F
ρ
t
+
k
cos
2
⁡
p
t
a
22
cos
⁡
p
t
+
a
23
s
i
n
p
t
2
h
11
ρ
3
t
,
t
∈
O
m
.
Let(34)
a
=
λ
m
a
x
F
>
0
,
b
t
=
-
k
cos
2
⁡
p
t
a
22
cos
⁡
p
t
+
a
23
s
i
n
p
t
2
h
11
≥
0
. Then, from (33), we have the equation(35)
ρ
˙
t
=
a
ρ
t
-
b
t
ρ
3
t
.The solution of this equation is(36)
ρ
t
=
ρ
0
e
x
p
-
2
a
t
+
2
ρ
0
2
e
x
p
-
2
a
t
∫
0
t
e
x
p
2
a
τ
b
τ
d
τ
.
In our case, we can introduce the function(37)
c
t
=
-
h
11
k
e
x
p
-
2
a
t
∫
0
t
e
x
p
2
a
τ
b
τ
d
τ
=
e
x
p
-
2
a
t
∫
0
t
e
x
p
2
a
τ
cos
2
⁡
p
τ
a
22
cos
⁡
p
τ
+
a
23
s
i
n
p
τ
2
d
τ
=
3
a
22
2
+
a
23
2
16
a
+
a
22
2
8
-
a
23
2
8
2
a
cos
⁡
4
p
t
+
4
p
s
i
n
4
p
t
4
a
2
+
16
p
2
+
a
22
2
2
2
a
cos
⁡
4
p
t
+
2
p
s
i
n
4
p
t
4
a
2
+
4
p
2
+
a
22
a
23
4
-
4
p
cos
⁡
4
p
t
+
2
a
s
i
n
4
p
t
4
a
2
+
16
p
2
+
a
22
a
23
2
-
2
p
cos
⁡
2
p
t
+
2
a
s
i
n
2
p
t
4
a
2
+
4
p
2
=
3
a
22
2
+
a
23
2
16
a
+
a
22
2
8
-
a
23
2
8
s
i
n
4
p
t
+
α
4
a
2
+
16
p
2
-
a
22
a
23
4
cos
⁡
4
p
t
+
α
4
a
2
+
16
p
2
+
a
22
2
2
s
i
n
2
p
t
+
β
4
a
2
+
4
p
2
-
a
22
a
23
2
cos
⁡
2
p
t
+
β
4
a
2
+
4
p
2
,where (38)
s
i
n
α
=
a
a
2
+
4
p
2
,
s
i
n
β
=
a
a
2
+
p
2
.
We take advantage of the known Cauchy inequality:∀
w
1
,
w
2
∈
R
3
,
〈
w
1
,
w
2
〉
≤
∥
w
1
∥
·
∥
w
2
∥; here, 〈
w
1
,
w
2
〉 is a scalar product of the vectors w
1 and w
2. Then, it is possible to get the following estimation of the function c
(
t
): (39)
c
t
>
3
a
22
2
+
a
23
2
16
a
-
a
22
2
+
a
23
2
16
a
2
+
4
p
2
-
4
a
22
a
22
2
+
a
23
2
16
a
2
+
p
2
≥
3
a
22
2
+
a
23
2
16
a
-
a
22
2
+
a
23
2
+
4
a
22
a
22
2
+
a
23
2
16
a
2
+
p
2
≥
3
a
22
2
+
a
23
2
16
a
-
a
22
2
+
a
23
2
+
2
a
22
2
16
a
2
+
p
2
>
0
.Thus, from here, we get the fourth condition of Theorem 3.
Further, formula (36) can be rewritten as (40)
ρ
t
=
ρ
0
e
x
p
-
2
a
t
-
2
k
ρ
0
2
c
t
/
h
11
.
Above, it was shown that if∀
t
≥
0 we have b
(
t
)
≥
0, then ∫
0
t
e
x
p
(
2
a
τ
)
b
(
τ
)
d
τ
>
0. Therefore, we can consider that there exist the positive constants c
m
i
n and c
m
a
x such that 0
<
c
m
i
n
≤
c
(
t
)
≤
c
m
a
x. Thus, from here and the condition t
→
∞, it follows that (41)
1
-
2
k
c
m
a
x
/
h
11
≤
ρ
t
≤
1
-
2
k
c
m
i
n
/
h
11and, therefore, the solution ρ
(
t
) of system (29) is bounded. Since h
11
<
0, then from the first equation of system (29) it follows that the solution x
(
t
) is also bounded. In turn, this means the boundedness of solutions of systems (19) and (20). The proof is finished.Note that ifa
22
=
0 and k
=
1, then system (19) is the classic Lorentz-like system and condition (iv) becomes trivial.Theorem 4.
Let the equilibrium point(
0,0
,
0
) of system (19) be a saddle-focus. Then, under the conditions of Theorem 3 in system (19), there exists a limit cycle.Proof.
Let us calculate Lyapunov’s exponentΛ for a real function f
(
t
) [17, 18]: (42)
Λ
f
=
lim
t
→
∞
-
⁡
1
t
ln
⁡
f
t
f
t
0
.
We take advantage of the following properties Lyapunov’s exponents:(i)
Λ
[
f
1
(
t
)
·
f
2
(
t
)
]
≤
Λ
[
f
1
(
t
)
]
+
Λ
[
f
2
(
t
)
]
.
(ii)
IfΛ
[
f
(
t
)
]
<
0, then Λ
[
∫
t
∞
f
(
τ
)
d
τ
]
≤
Λ
[
f
(
t
)
]
.
(iii)
Λ
[
f
1
(
t
)
+
f
2
(
t
)
]
≤
m
a
x
(
Λ
[
f
1
(
t
)
]
,
Λ
[
f
2
(
t
)
]
)
.
(iv)
Λ
[
d
·
f
(
t
)
]
=
Λ
[
f
(
t
)
]
(
d
≠
0
).
Taking into account the boundedness of solutionsx
(
t
) and ρ
(
t
) of system (28), we have Λ
[
x
]
≤
0, Λ
[
ρ
]
≤
0.
We write the first equation of system (28) in the integral form (43)
x
t
=
x
0
e
x
p
h
11
t
+
∫
0
t
e
x
p
h
11
t
-
τ
cos
⁡
ϕ
τ
a
22
cos
⁡
ϕ
τ
+
a
23
s
i
n
ϕ
τ
ρ
2
τ
d
τ
≡
x
0
e
x
p
h
11
t
+
∫
0
t
e
x
p
h
11
t
-
τ
ξ
τ
d
τ
;
ξ
t
=
cos
⁡
ϕ
t
a
22
cos
⁡
ϕ
t
+
a
23
s
i
n
ϕ
t
ρ
2
t
.
From here and (i), it follows thatΛ
(
e
x
p
(
h
11
t
)
ξ
(
t
)
)
≤
Λ
(
e
x
p
(
h
11
t
)
)
+
Λ
(
ξ
(
t
)
)
≤
h
11
+
0
=
h
11 or Λ
(
e
x
p
(
h
11
t
)
ξ
(
t
)
)
=
-
∞. Thus, from (ii) and (iii), it follows that (44)
Λ
∫
t
∞
exp
⁡
h
11
t
-
τ
ξ
τ
d
τ
≤
h
11
. Therefore, from boundedness of x
(
t
) and conditions (i), (iii), and (iv), we have Λ
[
x
(
t
)
]
≤
h
11
<
0.
Further, the functionρ
(
t
) is bounded. Consequently, we get that Λ
[
ρ
(
t
)
]
≤
0. In addition, from conditions (22), it follows that ϕ
(
t
)
→
p
t. Thus, we have Λ
[
ϕ
(
t
)
]
=
0. Since for the limit cycle another Lyapunov exponent is equal to zero, then from here it follows that if Λ
[
ρ
(
t
)
]
<
0, then this cycle is stable; if Λ
[
ρ
(
t
)
]
=
0, then we have a limit torus.Let again the point(
0,0
,
0
) be a saddle-focus. We note that if in system (19) we do the replacement of variables y
→
α
y
1
+
β
z
1
,
z
→
γ
y
1
+
δ
z
1 (α
δ
-
β
γ
≠
0) such that (45)
h
22
h
23
h
32
h
33
⟶
p
q
-
q
p
;
p
>
0
,
q
≠
0
,then the form g
(
u
,
v
) becomes positive definite. In this case, formula (21) can be represented in the form ϕ
(
t
)
=
ϕ
0
+
p
t
+
⋯ and condition (22) is always fulfilled. From here, it follows that, in system (28), which was used in the Comparison Principle, there exists a limit cycle.Let(46)
S
≔
x
,
y
,
z
∈
R
3
∣
x
+
s
2
Q
2
+
y
2
k
Q
2
+
z
2
Q
2
≤
1 be an ellipsoid with the center in the point (
-
s
,
0,0
), where (47)
s
=
-
k
h
32
+
h
23
k
a
23
,
Q
2
=
s
2
1
+
m
a
x
k
h
11
s
k
a
22
+
h
22
,
h
11
h
33
.Theorem 5.
Assume thath
21
=
h
31
=
0 and a
23
≠
0. Suppose also that h
11
<
0
,
s
k
a
22
+
h
22
<
0, and h
33
<
0. If the start point (
x
0
,
y
0
,
z
0
)
∈
S, then ∀
t
>
0 the trajectory (
x
(
t
,
x
0
)
,
y
(
t
,
y
0
)
,
z
(
t
,
z
0
)
) of system (19) is completely contained in S. Thus, the ellipsoid S is a basin of system (19).Proof.
LetV
(
t
)
=
(
x
(
t
)
+
s
)
2
+
y
2
(
t
)
/
k
+
z
2
(
t
)
≥
0 be a positive definite function. Then, from system (19), we have (48)
1
2
V
˙
t
=
x
+
s
x
˙
+
y
y
˙
k
+
z
z
˙
=
x
+
s
h
11
x
+
y
a
22
y
+
a
23
z
+
h
22
y
2
k
+
h
23
y
z
k
-
x
y
a
22
y
+
a
23
z
+
z
h
32
y
+
h
33
z
=
h
11
x
+
s
2
2
+
s
a
22
+
h
22
k
y
2
+
h
33
z
2
-
h
11
s
2
4
.
In order for the conditionV
˙
(
t
)
≤
0 to be fulfilled ∀
t
≥
0, it is necessary that(49)
-
h
11
x
+
s
2
2
-
s
k
a
22
+
h
22
k
y
2
-
h
33
z
2
≥
-
h
11
s
2
4
.
Let(50)
B
≔
x
,
y
,
z
∈
R
3
∣
x
+
s
/
2
2
P
2
+
y
2
k
P
2
+
z
2
P
2
<
1
,
w
h
e
r
e
P
2
=
Q
2
4
, be an open ellipsoid with the center in the point (
-
s
/
2,0
,
0
).
Let also(51)
B
1
≔
x
,
y
,
z
∈
R
3
∣
x
+
s
/
2
2
D
x
2
+
y
2
D
y
2
+
z
2
D
z
2
<
1 be the open ellipsoid with the center in the point (
-
s
/
2,0
,
0
) and semiaxes (52)
D
x
=
s
2
4
,
D
y
=
k
h
11
s
2
4
s
k
a
22
+
h
22
,
D
z
=
h
11
s
2
4
h
33
. (The ellipsoid B
1 is an interior part of domain (49); thus, domain (49) is the set R
3
-
B
1.)
It is clear thatB
1
⊂
B
⊂
S. This means that trajectories incoming in the ellipsoid S at all t
≥
0 remain in this ellipsoid. Consequently, these trajectories are bounded. The proof is finished.Notice that Theorem5 is weaker than Theorems 1, 3, and 4. However, verification of its conditions is constructive.


## 5. Examples

We will consider that for system (19) the conditions of either Theorem 3 or Theorem 4 are satisfied.At first, we compute equilibrium points of system (19). We have the equilibrium point (
0,0
,
0
) and two equilibrium points (
x
e
,
y
e
,
z
e
), where (53)
x
e
=
d
e
t
H
k
·
a
22
h
33
-
a
23
h
32
,
y
e
=
∓
h
33
h
32
-
h
11
h
33
h
32
2
d
e
t
H
k
·
a
22
h
33
-
a
23
h
32
2
,
z
e
=
±
-
h
11
h
33
h
32
2
d
e
t
H
k
·
a
22
h
33
-
a
23
h
32
2
.Thus, ifh
33
d
e
t
H
<
0, then there are three equilibrium points (including the point (
0,0
,
0
)); if h
33
d
e
t
H
>
0, then there is only one equilibrium point (
0,0
,
0
). It is easy to check that all equilibrium points will be either saddle-nodes or saddle-focuses or such points, Jacobi’s matrix of which has three positive eigenvalues.Further, we will study not only situations described in Theorems3 and 4, but also the cases h
21
≠
0 and h
31
≠
0. It is necessary with the changes of the parameters h
21 and h
31 to realize the cascade of the limit cycle bifurcations at transition from this cycle to the chaotic attractor.A special case arises, when the conditiona
22
h
33
-
a
23
h
32
=
0 takes place. In this case, there exists only one unstable equilibrium point (
0,0
,
0
).Consistently, consider all indicated situations:The case
a
22
h
33
-
a
23
h
32
≠
0.
The case
a
22
h
33
-
a
23
h
32
=
0. In this case, if the discriminant (
h
22
+
h
33
)
2
-
4
d
e
t
H
<
0, then system (19) demonstrates the existence of limit cycles. If (
h
22
+
h
33
)
2
-
4
d
e
t
H
>
0, then in system (19) no stable attractor was found.


## 6. Conclusion

In the present paper, an attempt to investigate the chaotic dynamics of system (19) was made. This system is a generalization of the known Lorenz system (1). Theorems 3 and 4, by which it is possible not only to build bounded solutions of system (19), but also to predict the existence of spatial limit cycles in this system, were also proved (see Figures 1(a), 2(a)–2(c), and 3(a)–3(c)). As a result, two new chaotic attractors (see Figures 3(d) and 4(f)) were found. Their chaotic behavior is confirmed by Theorem 1 (see Figure 4(d)). (At the same time, Figures 4(c) and 4(e) demonstrate the presence of the limit cycle in system (19) for the specific values of parameters; besides, Figures 4(a) and 4(b) confirm that condition (iv) of Theorem 3 is fulfilled.)Figure 1
The phase portrait of a quasi-periodic behavior (a) and the graph of the functionρ
(
t
) (b) of system (19) at h
11
=
-
0.1,  h
22
=
h
33
=
4,  h
23
=
-
h
32
=
40,  k
=
1,  a
22
=
-
1, and a
23
=
1.
(a)
(b)Figure 2
The phase portraits (a)–(d) of system (19) at h
11
=
-
2, h
22
=
-
0.5, h
23
=
30.5, h
32
=
9, h
33
=
-
2, k
=
1, and a
23
=
10 and different values of the parameter a
22.
(a)
a
22
=
1.0341
(b)
a
22
=
1.0342
(c)
a
22
=
1.0343
(d)
a
22
=
1.03432Figure 3
The phase portraits (a)–(d) of system (19) at h
11
=
-
2, h
22
=
-
0.5, h
23
=
30.5, h
32
=
5, h
33
=
-
2, k
=
1, a
22
=
0.1, and a
23
=
1 and different values of the parameter a
23 in the first equation.
(a)
a
23
=
1.0
(b)
a
23
=
5.0
(c)
a
23
=
9.0
(d)
a
23
=
10.0Figure 4
The graphs (a, b) of the functionϕ
(
t
); the graphs (c, d) of the function ρ
(
t
); and the phase portraits (e, f) of system (19) at h
11
=
-
1.2, h
22
=
h
33
=
1, h
23
=
-
h
32
=
4, k
=
1, a
22
=
-
1, and a
23
=
1 and different values of the parameters h
21
a
n
d
h
31.
(a)
h
21
=
h
31
=
0
(b)
h
21
=
-
3.0, h
31
=
6.64
(c)
h
21
=
h
31
=
0
(d)
h
21
=
-
3.0, h
31
=
6.64
(e)
h
21
=
-
3.0, h
31
=
6
(f)
h
21
=
-
3.0, h
31
=
6.64We note that Theorem1 was the basic result of [17, 18]. In the present work, for the systems of type (19), substantial strengthening of Theorem 1 was got. The point is that condition (12) is not constructible. Nevertheless, if, for the proof of the state of chaos for system (19), we will use new theorems (Theorems 3, 4, and 5), then the verification of condition (12) can be ignored.In addition, we did not begin to write the value of coefficients of system (19) (different from classic) at which in this system there is the Lorentz attractor. Researches of this attractor in detail are described in numerous articles and books.



---

*Source: 2410806-2016-09-20.xml*