# Analysis, Stabilization, and DSP-Based Implementation of a Chaotic System with Nonhyperbolic Equilibrium

Authors: Xuan-Bing Yang; Yi-Gang He; Chun-Lai Li; Chang-Qing Liu
Journal: Complexity
 (2020)

---

## Abstract

This paper reports an autonomous dynamical system, and it finds that one nonhyperbolic zero equilibrium and two hyperbolic nonzero equilibria coexist in this system. Thus, it is difficult to demonstrate the existence of chaos by Šil’nikov theorem. Consequently, the topological horseshoe theory is adopted to rigorously prove the chaotic behaviors of the system in the phase space of Poincaré map. Then, a single control scheme is designed to stabilize the dynamical system to its zero-equilibrium point. Besides, to verify the theoretical analyses physically, the attractor and stabilization scheme are further realized via DSP-based technique.


---

## Body

## 1. Introduction

Chaotic behaviors exist widely in biology, engineering, economy, and many other scientific disciplines [1–4]. Chaotic systems have the properties of unpredictability, topological mixing, ergodicity, and sensitivity to their initial values and control parameters. Owing to the noise-like spectrum and broad-band, chaotic signals are potentially applicable in engineering such as random sequence generation [5], secure communication [6], image encryption [7, 8], signal detection [9], radar and sonar systems [10], and so on. Therefore, it is significant to design and analyze new chaotic systems. In 1963, Lorenz proposed the first three-dimensional chaotic system when studying earth’s atmospheric convection [11]. Then, many low-dimensional mathematical models with positive Lyapunov exponents have been introduced, along with the analysis of a rich class of dynamical behaviors. New examples continue to be reported for publication in nonlinear dynamics journals for the physical applications. Nevertheless, it is necessary to develop chaotic systems with simple structure and rich dynamical behaviors from the perspective of application.As a striking chaos theory with symbolic dynamics obtained by Kennedy in continuous map [12], the topological horseshoe can provide an impactful tool for proving chaotic dynamics in hyperbolic or nonhyperbolic system. Up to present, many noteworthy theoretical progresses have been extended in finding the existence of horseshoe. For example, Yang introduced the remarkable criteria for finding the topological horseshoe in noncontinuous map [13, 14], which has been successfully applied to some practical systems for verifying chaos [15, 16]. Li presented a new method with three steps for finding horseshoes in dynamical systems by using several simple results on topological horseshoes [17]. However, it is still a tough work to find a topological horseshoe in a practical chaotic system [18].It is impossible to chronically predict the future behavior of chaotic system, but one can stabilize the future behavior into a certain range by using control technology. The seminal attempt of controlling chaos is the well-known OGY method, which applies small perturbations to system parameter to keep the system close to the target periodic orbit [19]. However, the experimental implementation of this method is restricted by the level of noise in the experimental data for the sake of the discrete nature of the control signal. Then, a continuous delayed feedback scheme is proposed consequently [20]. The control signal is the perturbations of the difference between the states of the current system and one period of the target orbit in the past. Therefore, the intensity of the perturbations will vanish when the system evolves to the desired orbit. Since then, a wide variety of control approaches emerged for different applications, such as backstepping control [21], sliding mode control [22], sampled-data control [23], and multiswitching combination control [24], just to name a few. However, the present schemes focused on the control problem need several controllers. From the points of both practical application and theoretical research, it is significant to design simple yet executable control technique.In this paper, an autonomous chaotic system with a simple algebraic structure of six terms is proposed. Basic dynamical properties of the system, including equilibrium point, phase portrait, Poincaré map, parameter bifurcation, and Lyapunov exponent, are studied in theory and numerical simulation. It is found that this system exhibits fruitful dynamic behaviors of dense periodic windows and coexistence of nonhyperbolic and hyperbolic equilibrium points. And to rigorously verify the emergence of chaos of this system in theory, the topological horseshoe is investigated in the phase space of Poincaré map. Then, based on the Lyapunov stability criterion, a single control scheme is designed to stabilize the chaotic system to its zero-equilibrium point. The implementation scheme of attractors and control scheme are discussed in detail and realized via DSP-based technique, confirming the validity and enforceability of the theoretical scheme.


## 2. The Proposed Dynamical System

The autonomous dynamical system considered here is given by(1)x˙1=−ax1+bx2,x˙2=ex2−fx1x3,x˙3=−cx33+dx1x2.Differential equations of system (1) are simple with a cubic term and two quadratic cross-product terms. The system parameters a, b, c, d, e, and f are all positive constants.


### 2.1. Dissipativity and Existence of Attractor

We first consider the general condition of dissipativity to ensure the chaotic property:(2)∇V=∂x˙1x1+∂x˙2x2+∂x˙3x3=−a+e−3cx32.Therefore, system (1) would be dissipative and will converge to a subset of measure zero volume according to dV/dt=e−a+e−3cx32, when satisfying −a+e−3cx32<0. This means that the volume will become V0e−a+e−3cx32t at time t through the flow generated by the system for an initial volume V0. Therefore, there exists an attractor in system (1) with −a+e−3cx32<0.


### 2.2. Equilibrium Points and Stability

Considering the condition of equilibrium pointx˙1=0, x˙2=0, and x˙3=0, we obtain three equilibrium points of system (1), as follows:(3)P00,0,0,P1abce3df3,a2b2ce3df3,aebf,P2−abce3df3,−a2b2ce3df3,aebf.When selectinga = 1, b = 6, c = 5, d = 1, e = 3, and f = 2, the equilibrium points and the corresponding eigenvalues are shown in Table 1, including the type of equilibrium points.Table 1
Equilibrium points and eigenvalues of system (1).
Equilibrium point
Eigenvalues
Type of equilibrium point
P
0 (0, 0, 0)
λ
1 = −1, λ2 = 3, λ3 = 0
Nonhyperbolic
P
1 (0.6847, 0.1141, 0.25)
λ
1 = 1.1087 + 0.6279i,λ2 = 1.1087 – 0.6279i,λ3 = −1.1549
Hyperbolic
P
2 (–0.6847, -0.1141, 0.25)
λ1 = 1.1087 + 0.6279i,λ2 = 1.1087 – 0.6279i,λ3 = −1.1549
HyperbolicFrom Table1, it is known that the three equilibrium points are all unstable with stable manifold and unstable manifold. According to [25], since the characteristic value λ3 of equilibrium point P0 equals to zero, the equilibrium point is nonhyperbolic type. However, equilibrium points P1 and P2 are hyperbolic since all the real parts of the corresponding eigenvalues of these two equilibrium points are nonzero. Therefore, this is a chaotic system in which hyperbolic and nonhyperbolic equilibrium points coexist. And, this kind of chaotic system does not belong to Šil’nikov sense of the chaotic system, and it is difficult to prove the existence of chaos by Šil’nikov theorem. At present, chaotic systems with special features, such as chaotic systems with nonequilibrium, with multistability and with hyperbolic and nonhyperbolic equilibrium coexisting, have attracted extensive attention by researchers [26, 27].


### 2.3. Phase Portrait and Chaotic Properties

When selectinga = 1, b = 6, c = 5, d = 1, e = 3, f = 2, and initial values x (0) = (0.01, 0.01, 0.05), the Lyapunov exponents of system (1) are depicted in Figures 1(a) and 1(b) with the values as 0.123017, 0.000029, and −11.414797. According to the three Lyapunov exponents, the Kaplan–Yorke dimension is DKY = 2 + (0.123017 + 0.000029)/11.414797 = 2.0108. Therefore, the Kaplan–Yorke dimension is fractional. The corresponding chaotic phase diagrams further reveal that the proposed system displays complicated chaotic behaviors, as depicted in Figure 2.Figure 1
Lyapunov exponents of system (1) when a = 1, b = 6, c = 5, d = 1, e = 3, and f = 2.
(a)
(b)Figure 2
Chaotic phase diagrams of the proposed system.
(a)
(b)
(c)As an important analytical technique, Poincaré map can reflect the bifurcation and folding properties of chaos. Selecting the parameter valuesa = 1, b = 6, c = 5, d = 1, e = 3, and f = 2 and taking the crossing planes x3 = 0, x1 = 0, and x2 = 0, we obtain the corresponding Poincaré maps illustrated in Figure 3. We can see the attractor structure from the dense dots.Figure 3
Poincaré maps on the plane of (a)x3 = 0, (b) x1 = 0, and (c) x2 = 0.
(a)
(b)
(c)


### 2.4. Influence of Parameter Variation

It is found that system (1) exhibits complicated dynamical behaviors with the variation of system parameters. As explication, we only consider the variation of parameters c and f in this section.When fixing parameters {a = 1, b = 6, d = 1, e = 3, f = 2} and letting c vary from 2 to 8, the corresponding bifurcation diagram of state variable x3 and the maximum Lyapunov exponent of system (1) versus c of numerical calculation are shown in Figures 4(a) and 4(b), respectively. As we know that system (1) shows rich dynamical behaviors, ranging from stable equilibrium points, periodic orbits to chaotic oscillations, depending on the parameter values. Furthermore, there emerge many visible periodic windows in the chaotic region. Then, we fix parameters a = 1, b = 6, c = 5, d = 1, and e = 3, while letting f vary from 1 to 3. Figures 5(a) and 5(b) depict the bifurcation diagram of state variable x3 and the maximum Lyapunov exponent of system (1) versus parameter f, respectively. It is known that with the variation of parameter f, system (1) ranges from stable equilibrium points, periodic orbits to chaotic oscillations, and there emerges many densely distributed periodic windows in the chaotic region, also showing the rich dynamics. The variation of the properties of system (1) with system parameters c and f is of great importance in image encryption.Figure 4
(a) Bifurcation diagram and (b) Lyapunov exponents of system (1) versus parameter c.
(a)
(b)Figure 5
(a) Bifurcation diagram and (b) Lyapunov exponents of system (1) versus parameter f.
(a)
(b)


## 2.1. Dissipativity and Existence of Attractor

We first consider the general condition of dissipativity to ensure the chaotic property:(2)∇V=∂x˙1x1+∂x˙2x2+∂x˙3x3=−a+e−3cx32.Therefore, system (1) would be dissipative and will converge to a subset of measure zero volume according to dV/dt=e−a+e−3cx32, when satisfying −a+e−3cx32<0. This means that the volume will become V0e−a+e−3cx32t at time t through the flow generated by the system for an initial volume V0. Therefore, there exists an attractor in system (1) with −a+e−3cx32<0.


## 2.2. Equilibrium Points and Stability

Considering the condition of equilibrium pointx˙1=0, x˙2=0, and x˙3=0, we obtain three equilibrium points of system (1), as follows:(3)P00,0,0,P1abce3df3,a2b2ce3df3,aebf,P2−abce3df3,−a2b2ce3df3,aebf.When selectinga = 1, b = 6, c = 5, d = 1, e = 3, and f = 2, the equilibrium points and the corresponding eigenvalues are shown in Table 1, including the type of equilibrium points.Table 1
Equilibrium points and eigenvalues of system (1).
Equilibrium point
Eigenvalues
Type of equilibrium point
P
0 (0, 0, 0)
λ
1 = −1, λ2 = 3, λ3 = 0
Nonhyperbolic
P
1 (0.6847, 0.1141, 0.25)
λ
1 = 1.1087 + 0.6279i,λ2 = 1.1087 – 0.6279i,λ3 = −1.1549
Hyperbolic
P
2 (–0.6847, -0.1141, 0.25)
λ1 = 1.1087 + 0.6279i,λ2 = 1.1087 – 0.6279i,λ3 = −1.1549
HyperbolicFrom Table1, it is known that the three equilibrium points are all unstable with stable manifold and unstable manifold. According to [25], since the characteristic value λ3 of equilibrium point P0 equals to zero, the equilibrium point is nonhyperbolic type. However, equilibrium points P1 and P2 are hyperbolic since all the real parts of the corresponding eigenvalues of these two equilibrium points are nonzero. Therefore, this is a chaotic system in which hyperbolic and nonhyperbolic equilibrium points coexist. And, this kind of chaotic system does not belong to Šil’nikov sense of the chaotic system, and it is difficult to prove the existence of chaos by Šil’nikov theorem. At present, chaotic systems with special features, such as chaotic systems with nonequilibrium, with multistability and with hyperbolic and nonhyperbolic equilibrium coexisting, have attracted extensive attention by researchers [26, 27].


## 2.3. Phase Portrait and Chaotic Properties

When selectinga = 1, b = 6, c = 5, d = 1, e = 3, f = 2, and initial values x (0) = (0.01, 0.01, 0.05), the Lyapunov exponents of system (1) are depicted in Figures 1(a) and 1(b) with the values as 0.123017, 0.000029, and −11.414797. According to the three Lyapunov exponents, the Kaplan–Yorke dimension is DKY = 2 + (0.123017 + 0.000029)/11.414797 = 2.0108. Therefore, the Kaplan–Yorke dimension is fractional. The corresponding chaotic phase diagrams further reveal that the proposed system displays complicated chaotic behaviors, as depicted in Figure 2.Figure 1
Lyapunov exponents of system (1) when a = 1, b = 6, c = 5, d = 1, e = 3, and f = 2.
(a)
(b)Figure 2
Chaotic phase diagrams of the proposed system.
(a)
(b)
(c)As an important analytical technique, Poincaré map can reflect the bifurcation and folding properties of chaos. Selecting the parameter valuesa = 1, b = 6, c = 5, d = 1, e = 3, and f = 2 and taking the crossing planes x3 = 0, x1 = 0, and x2 = 0, we obtain the corresponding Poincaré maps illustrated in Figure 3. We can see the attractor structure from the dense dots.Figure 3
Poincaré maps on the plane of (a)x3 = 0, (b) x1 = 0, and (c) x2 = 0.
(a)
(b)
(c)


## 2.4. Influence of Parameter Variation

It is found that system (1) exhibits complicated dynamical behaviors with the variation of system parameters. As explication, we only consider the variation of parameters c and f in this section.When fixing parameters {a = 1, b = 6, d = 1, e = 3, f = 2} and letting c vary from 2 to 8, the corresponding bifurcation diagram of state variable x3 and the maximum Lyapunov exponent of system (1) versus c of numerical calculation are shown in Figures 4(a) and 4(b), respectively. As we know that system (1) shows rich dynamical behaviors, ranging from stable equilibrium points, periodic orbits to chaotic oscillations, depending on the parameter values. Furthermore, there emerge many visible periodic windows in the chaotic region. Then, we fix parameters a = 1, b = 6, c = 5, d = 1, and e = 3, while letting f vary from 1 to 3. Figures 5(a) and 5(b) depict the bifurcation diagram of state variable x3 and the maximum Lyapunov exponent of system (1) versus parameter f, respectively. It is known that with the variation of parameter f, system (1) ranges from stable equilibrium points, periodic orbits to chaotic oscillations, and there emerges many densely distributed periodic windows in the chaotic region, also showing the rich dynamics. The variation of the properties of system (1) with system parameters c and f is of great importance in image encryption.Figure 4
(a) Bifurcation diagram and (b) Lyapunov exponents of system (1) versus parameter c.
(a)
(b)Figure 5
(a) Bifurcation diagram and (b) Lyapunov exponents of system (1) versus parameter f.
(a)
(b)


## 3. Topological Horseshoe in the Dynamical System

### 3.1. Review of Topological Horseshoe Theorems

It is still a challenge to find topological horseshoes in a concrete system, especially to select a suitable quadrilateral in the cross section. Before studying the horseshoe embedded in the dynamical system, some theorems on topological horseshoe are reviewed below [12–14, 17].LetD be a compact subset of S, which is a metric space, and there exists m mutually disjoint compact subsets D1, D2, …, Dm of D.Definition 1.
LetDi1,Di2⊂Di be two fixed disjoint compact subsets with 1≤i≤m. If γ∩Di1 and γ∩Di2 are compact and nonempty, we say that the connected subset γ of Di connects Di1 and Di2, and denote this by Di1↔γDi2.Definition 2.
Letγ be a connected subset of Di, we say that fγ is suitably across Di with respect to Di1 and Di2, if there is a connected subset γi⊂γ satisfying fγi⊂Di, and fγi∩Di1 and fγi∩Di2 are nonempty. In this case, it is denoted by fγ⟼Di.Definition 3.
Suppose thatσ:Σm⟶Σm and f:S⟶S are both continuous functions with topological spaces Σm and S, respectively. If there exists a continuous surjection h:Σm⟶S confirming to f∘h=h∘σ, it is said that f is topologically semiconjugate to function σ.Lemma 1.
IffpDi⟼Di, then we have fmpDi⟼Di, where m is a positive integer.Lemma 2.
IffpD1⟼D1, fpD1⟼D2, fqD2⟼D1, and fqD2⟼D2, then there would exist a compact invariant set M⊂D, such that fp+qM is semiconjugate to 2-shift dynamics and the topological entropy of f will satisfy entf≥1/p+qlog2.


### 3.2. Finding Topological Horseshoe in the Dynamical System

According to the theory above, a horseshoe will be found in the dynamical system by three steps [17]. In this process, we set the parameters of system (1) as a = 1, b = 6, c = 5, d = 1, e = 3, f = 2, and initial condition x (0) = (0.01, 0.01, 0.05).Step 1.
As shown in Figure6, we first select four vertices of the Poincaré section P on the plane Θ=x1,x2,x3∈R3:x2=0 as (–8, 0, –0.5), (–8, 0, 1.5), (8, 0, 1.5), and (8, 0, –0.5).Figure 6
Poincaré cross section of the proposed system.Step 2.
Then, after many trial-and-error numerical simulations, we carefully pick a quadrilateralD1 of quadrangle P, with the four vertices being(4)−7.618610595,0,1.062250623−7.643703532,0,1.049345387−7.239428439,0,1.018765586−7.215729554,0,1.031951372.
Let us suppose thatD11 denotes the left side while D12 denotes the right side of quadrilateral D1. It is known from the numerical result that the third return map H3D11 lies on the left side of D1, but the third return map H3D12 lies on the right side of D1. Thereby, under this return map, the image H3x (x∈D1) lies wholly across the quadrangle D1 with respect to the sides D11 and D12, seen in Figure 7(a).Figure 7
(a) The subsetD1 and its image; (b) and (c) the subset D2 and its image.
(a)
(b)
(c)Step 3.
Finally, we will take another quadrilateralD2 of P such that HD1⟼D2, HD2⟼D2 and HD2⟼D1. By a great deal of attempts, the four vertexes of D2 are picked as follows:(5)−7.865357807,0,1.082029302−7.896026952,0,1.066318579−7.661826208,0,1.048363466−7.633945167,0,1.063232544.
Analogously,D21 and D22 indicate the left and right sides of quadrilateral D2, respectively. The numerical simulations of the third return Poincaré map H3D2 and the enlarged view are depicted in Figures 7(b) and 7(c), respectively. It is shown from the figures that the return map H3x (x∈D2) suitably across the quadrangles D1 and D2, with H3D21 lying on the right side of D1 and H3D22 lying on the left side of D2.
Therefore, we conclude by virtue of Lemma2 that there would exist a compact invariant set M⊂D, such that H6M is semiconjugate to 2-shift dynamics, and we obtain entH≥log2/6>0. Thus, system (1) is proved to be chaotic in theory, with the parameters of a = 1, b = 6, c = 5, d = 1, e = 3, f = 2, and the initial condition of x (0) = (0.01, 0.01, 0.05).


## 3.1. Review of Topological Horseshoe Theorems

It is still a challenge to find topological horseshoes in a concrete system, especially to select a suitable quadrilateral in the cross section. Before studying the horseshoe embedded in the dynamical system, some theorems on topological horseshoe are reviewed below [12–14, 17].LetD be a compact subset of S, which is a metric space, and there exists m mutually disjoint compact subsets D1, D2, …, Dm of D.Definition 1.
LetDi1,Di2⊂Di be two fixed disjoint compact subsets with 1≤i≤m. If γ∩Di1 and γ∩Di2 are compact and nonempty, we say that the connected subset γ of Di connects Di1 and Di2, and denote this by Di1↔γDi2.Definition 2.
Letγ be a connected subset of Di, we say that fγ is suitably across Di with respect to Di1 and Di2, if there is a connected subset γi⊂γ satisfying fγi⊂Di, and fγi∩Di1 and fγi∩Di2 are nonempty. In this case, it is denoted by fγ⟼Di.Definition 3.
Suppose thatσ:Σm⟶Σm and f:S⟶S are both continuous functions with topological spaces Σm and S, respectively. If there exists a continuous surjection h:Σm⟶S confirming to f∘h=h∘σ, it is said that f is topologically semiconjugate to function σ.Lemma 1.
IffpDi⟼Di, then we have fmpDi⟼Di, where m is a positive integer.Lemma 2.
IffpD1⟼D1, fpD1⟼D2, fqD2⟼D1, and fqD2⟼D2, then there would exist a compact invariant set M⊂D, such that fp+qM is semiconjugate to 2-shift dynamics and the topological entropy of f will satisfy entf≥1/p+qlog2.


## 3.2. Finding Topological Horseshoe in the Dynamical System

According to the theory above, a horseshoe will be found in the dynamical system by three steps [17]. In this process, we set the parameters of system (1) as a = 1, b = 6, c = 5, d = 1, e = 3, f = 2, and initial condition x (0) = (0.01, 0.01, 0.05).Step 1.
As shown in Figure6, we first select four vertices of the Poincaré section P on the plane Θ=x1,x2,x3∈R3:x2=0 as (–8, 0, –0.5), (–8, 0, 1.5), (8, 0, 1.5), and (8, 0, –0.5).Figure 6
Poincaré cross section of the proposed system.Step 2.
Then, after many trial-and-error numerical simulations, we carefully pick a quadrilateralD1 of quadrangle P, with the four vertices being(4)−7.618610595,0,1.062250623−7.643703532,0,1.049345387−7.239428439,0,1.018765586−7.215729554,0,1.031951372.
Let us suppose thatD11 denotes the left side while D12 denotes the right side of quadrilateral D1. It is known from the numerical result that the third return map H3D11 lies on the left side of D1, but the third return map H3D12 lies on the right side of D1. Thereby, under this return map, the image H3x (x∈D1) lies wholly across the quadrangle D1 with respect to the sides D11 and D12, seen in Figure 7(a).Figure 7
(a) The subsetD1 and its image; (b) and (c) the subset D2 and its image.
(a)
(b)
(c)Step 3.
Finally, we will take another quadrilateralD2 of P such that HD1⟼D2, HD2⟼D2 and HD2⟼D1. By a great deal of attempts, the four vertexes of D2 are picked as follows:(5)−7.865357807,0,1.082029302−7.896026952,0,1.066318579−7.661826208,0,1.048363466−7.633945167,0,1.063232544.
Analogously,D21 and D22 indicate the left and right sides of quadrilateral D2, respectively. The numerical simulations of the third return Poincaré map H3D2 and the enlarged view are depicted in Figures 7(b) and 7(c), respectively. It is shown from the figures that the return map H3x (x∈D2) suitably across the quadrangles D1 and D2, with H3D21 lying on the right side of D1 and H3D22 lying on the left side of D2.
Therefore, we conclude by virtue of Lemma2 that there would exist a compact invariant set M⊂D, such that H6M is semiconjugate to 2-shift dynamics, and we obtain entH≥log2/6>0. Thus, system (1) is proved to be chaotic in theory, with the parameters of a = 1, b = 6, c = 5, d = 1, e = 3, f = 2, and the initial condition of x (0) = (0.01, 0.01, 0.05).


## 4. Stabilization for the Dynamical System

### 4.1. Control Scheme

In order to stabilize the proposed dynamical system, we add the single controlleru on the second equation. Thus, the controlled dynamical system is depicted as follows:(6)x˙1=−ax1+bx2,x˙2=ex2−fx1x3+u,x˙3=−cx33+dx1x2.The purpose of our design is to propose suitable control schemeu such that all the output variables of system (6) converge to the zero equilibrium point asymptotically.Theorem 1.
For the controlled system (6), we design the single controller u as(7)u=Ksignx2.If the control gain satisfiesK≤−b2/4a+eB2 with B2≥x2, then the output variables of the controlled system (6) converges to the zero equilibrium point asymptotically.Proof.
We choose the candidate Lyapunov function as(8)V=0.5x12+x22+fx32d.
The corresponding time derivative ofVx is deduced by(9)V˙=x1x˙1+x2x˙2+fx3x˙3d=−ax12+bx1x2+ex22−fx1x2x3+Kx2signx2−cfx34d+fx1x2x3=−ax12+bx1x2+ex22+Kx2−cfx34d=−ax12−bx1x2+b2ax22+b2ax22+ex22+Kx2−cfx34d=−ax1−b2ax22+b24a+ex2+Kx2−cfx34d.
Thus, whenK≤−b2/4a+eB2, we find that(10)V˙≤−ax1−b2ax22+b24a+eB2+Kx2−cfx34d≤−ax1−b2ax22−cfx34d≤0.
Therefore, the output variables of the controlled system (6) will converge to the zero equilibrium point asymptotically.


### 4.2. Numerical Verification

In this section, numerical simulations are executed by adopting the ODE45 method to verify the availability of the proposed control scheme. For comparing conveniently, we choose the system parameters asa = 1, b = 6, c = 5, d = 1, e = 3, and f = 2, and the initial states are set as x (0) = (0.01, 0.01, 0.01). With the parameter set, system (6) is chaotic before the controller is put into effect.The controlleru is exerted at 60th second, and the control gain is taken as K = –2 and K = –16. Figures 8(a) and 8(b) depict the control results, from which one can see that we can obtain good control effect with only taking a short time for the system to be stabilized at the fixed point. And, as we can also know only small control energy is needed to reach our control purpose.Figure 8
Control results of state variables with (a)K = –2 and (b) K = –16.
(a)
(b)


## 4.1. Control Scheme

In order to stabilize the proposed dynamical system, we add the single controlleru on the second equation. Thus, the controlled dynamical system is depicted as follows:(6)x˙1=−ax1+bx2,x˙2=ex2−fx1x3+u,x˙3=−cx33+dx1x2.The purpose of our design is to propose suitable control schemeu such that all the output variables of system (6) converge to the zero equilibrium point asymptotically.Theorem 1.
For the controlled system (6), we design the single controller u as(7)u=Ksignx2.If the control gain satisfiesK≤−b2/4a+eB2 with B2≥x2, then the output variables of the controlled system (6) converges to the zero equilibrium point asymptotically.Proof.
We choose the candidate Lyapunov function as(8)V=0.5x12+x22+fx32d.
The corresponding time derivative ofVx is deduced by(9)V˙=x1x˙1+x2x˙2+fx3x˙3d=−ax12+bx1x2+ex22−fx1x2x3+Kx2signx2−cfx34d+fx1x2x3=−ax12+bx1x2+ex22+Kx2−cfx34d=−ax12−bx1x2+b2ax22+b2ax22+ex22+Kx2−cfx34d=−ax1−b2ax22+b24a+ex2+Kx2−cfx34d.
Thus, whenK≤−b2/4a+eB2, we find that(10)V˙≤−ax1−b2ax22+b24a+eB2+Kx2−cfx34d≤−ax1−b2ax22−cfx34d≤0.
Therefore, the output variables of the controlled system (6) will converge to the zero equilibrium point asymptotically.


## 4.2. Numerical Verification

In this section, numerical simulations are executed by adopting the ODE45 method to verify the availability of the proposed control scheme. For comparing conveniently, we choose the system parameters asa = 1, b = 6, c = 5, d = 1, e = 3, and f = 2, and the initial states are set as x (0) = (0.01, 0.01, 0.01). With the parameter set, system (6) is chaotic before the controller is put into effect.The controlleru is exerted at 60th second, and the control gain is taken as K = –2 and K = –16. Figures 8(a) and 8(b) depict the control results, from which one can see that we can obtain good control effect with only taking a short time for the system to be stabilized at the fixed point. And, as we can also know only small control energy is needed to reach our control purpose.Figure 8
Control results of state variables with (a)K = –2 and (b) K = –16.
(a)
(b)


## 5. DSP-Based Realization of the Attractors and Stabilization Scheme

The implementation of continuous chaotic systems with electronic circuit has been widely adopted. The parameter tolerance of the components of analog electronic circuit will cause the trajectory change of the chaotic system, which restricts its practical engineering application. Therefore, DSP-based implement of continuous chaotic system can overcome these problems effectively. In this section, the attractors of system (1) and stabilization scheme for system (6) are implemented based on the DSP platform. In our experiments, the Texas Instrument DSP TMS320F28335 is employed to calculate the state variable and the control variable, which can run at 150 MHz and interfaces with a 12-bit quad-channel digital-to-analog converter DAC7724 by parallel bus (PB) mode. Control signals required by DAC7724 are generated by DSP and CPLD chips. The block diagram of hardware platform is shown in Figure 9.Figure 9
Block diagram for DSP implementation of chaotic attractors and its stabilization.In practice, we discretize the continuous system by the classical fourth-order Runge–Kutta algorithm with the sampling periodΔT. From Figure 2, we know that the amplitude of variable x3 (x3<2) is too small to affect the computational precision. Therefore, we proposed a general method to rescale the system variables by scaling factors ki, where i=1,2,3.Settingui=kixi,i=1,2,3, system (1) can be transformed to(11)u˙i=fia,b,c,d,e,f,k1,k2,k3,u1,u2,u3,where i=1,2,3.When introducing the linear transformation ofui⟶xi, the corresponding differential equation of system (11) can be expressed as(12)x˙1=−ax1+bk1k2x2,x˙2=ex2−fk2k1k3x1x3,x˙3=−ck32x33+dk3k1k2x1x2.Discretizing system (12) using the classical fourth-order Runge–Kutta algorithm, the following difference equations can be obtained:(13)xin+1=xin+ΔTKi,1+2Ki,2+2Ki,3+Ki,46,where i=1,2,3 and(14)K1,1=−ax1n+bk1k2x2n,K2,1=ex2n−fk2k1k3x1nx3n,K3,1=−ck32x33n+dk3k1k2x1nx2n,(15)K1,j=−ax1n+mK1,j−1+bk1k2x2n+mK2,j−1,K2,j=ex2n+mK2,j−1−fk2k1k3x1n+mK1,j−1x3n+mK3,j−1,K3,j=−ck32x3n+mK3,j−13+dk3k1k2x1n+mK1,j−1x2n+mK2,j−1.Whenj=2,3, m=0.5, and when j=4, m=1.The digital sequencesx1n,x2n,x3n are produced by DSP according to (13)–(15) and converted into the corresponding analog sequences by the DAC7724 and then transmitted to oscilloscope. With the set a = 1, b = 6, c = 5, d = 1, e = 3, f = 2, k1=1,k2=1,k3=2, x (0) = (0.01, 0.01, 0.01), and the control gain K = –6, the results of DSP-based realization of the attractor and stabilization scheme are depicted in Figure 10. By comparing with the attractors in Figure 2 and stabilization process in Figure 7 simulated by Matlab, it can be concluded that they have a good qualitative agreement.Figure 10
DSP-based realization for (a)–(c) chaotic attractors; (d) stabilization process; (e) experiment platform.
(a)
(b)
(c)
(d)
(e)


## 6. Conclusion

In this paper, we presented a three-dimensional dynamical system with a simple algebraic structure. Basic dynamical properties of the system, including equilibrium point, phase portrait, Poincaré map, parameter bifurcation, and Lyapunov exponent, are studied through theoretical analysis and numerical simulation. And, the theory of topological horseshoe is adopted to rigorously prove the chaotic emergence of the system theoretically. Then, based on the Lyapunov stability criterion, we designed a single control scheme to stabilize the chaotic system to its zero-equilibrium point. The attractor and stabilization process are realized via DSP-based technology, which have a good qualitative agreement to the Matlab simulation; thus it well confirmed the validity and enforceability of the theoretical scheme.



---

*Source: 6786832-2020-01-08.xml*