*Article* **Tracking Control of a Hyperchaotic Complex System and Its Fractional-Order Generalization**

**Feng Liang 1,2, Lu Lu 2,† , Zhengfeng Li 3 , Fangfang Zhang 3, \* and Shuaihu Zhang 3**


**Abstract:** Hyperchaotic complex behaviors often occur in nature. Some chaotic behaviors are harmful, while others are beneficial. As for harmful behaviors, we hope to transform them into expected behaviors. For beneficial behaviors, we want to enhance their chaotic characteristics. Aiming at the harmful hyperchaotic complex system, a tracking controller was designed to produce the hyperchaotic complex system track common expectation system. We selected sine function, constant, and complex Lorenz chaotic system as target systems and verified the effectiveness by mathematical proof and simulation experiments. Aiming at the beneficial hyperchaotic complex phenomenon, this paper extended the hyperchaotic complex system to the fractional order because the fractional order has more complex dynamic characteristics. The influences order change and parameter change on the evolution process of the system were analyzed and observed by MATLAB simulation.

**Keywords:** hyperchaotic complex system; tracking control; fractional order

**Citation:** Liang, F.; Lu, L.; Li, Z.; Zhang, F.; Zhang, S. Tracking Control of a Hyperchaotic Complex System and Its Fractional-Order Generalization. *Processes* **2022**, *10*, 1244. https://doi.org/10.3390/ pr10071244

Academic Editor: Zhiwei Gao

Received: 25 May 2022 Accepted: 16 June 2022 Published: 22 June 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Copyright:** © 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

#### **1. Introduction**

Chaos is a complex nonlinear phenomenon in nature. The chaotic system has extensive applications in the fields of secure communication [1], industrial process [2], ecosystem [3], and so on. With the application of chaos theory in increasingly more fields, there are increasingly more requirements for chaotic systems. For example, people expect higher dimensional and more complex chaotic systems to describe industrial processes. Therefore, following this, scholars put forward the hyperchaotic system, complex chaotic system, and hyperchaotic complex system [4–6]. The hyperchaotic complex system in particular has higher dimensions and more controllable parameters [7,8] that can more accurately describe some chaotic phenomena in the industrial process. The chaotic system will experience four states: stable point, period, chaos, and divergence. When the system is in the periodic or period doubling state, the phase plane will form a closed trajectory, that is, the limit cycle. Generally speaking, the chaotic system will evolve from a periodic state to a chaotic state, and the corresponding phase diagram will evolve from limit cycle to chaotic attractor [9].

However, some chaotic behaviors are harmful, while some are beneficial. As for harmful behaviors, we hope to transform them into expected behaviors. For example, permanent magnet motor under some parameters can produce chaotic behaviors and disturb the normal operation of the motor. Therefore, it is necessary to add a controller to make it track the desired motion trajectory [10,11]. This means the tracking control for chaos, which can be used to obtain the desired output and improve the performance of the system [12,13]. The research on the tracking control of hyperchaotic complex systems is of great value.

At present, most of the research focuses on the tracking control of real chaotic systems [14–17], and few scholars have studied the tracking control of hyperchaotic complex systems. Gao proposed a novel tracking control method for Lorenz systems by using single-state feedback [18]. Loria addressed the problem of controlled synchronization of a class of uncertain chaotic systems [19]. Zhang presented the tracking control method and the parameter identification procedure, aiming at CVCSs with complex parameters [20]. Chaudhary et al. investigated a hybrid projective combination–combination synchronization scheme (HPCCSS) in four different hyperchaotic (HC) systems via active control technique (ACT) [21]. Abbasi proposed a robust resilient design methodology for stabilization and tracking control for a class of chaotic dynamical systems [22]. Zhao realized tracking control and synchronization of the fractional hyper chaotic Lorenz system [23]. Nagy et al. investigated the combination synchronization phenomena of various fractionalorder systems using the scaling matrix [24]. In this paper, the tracking control of complex hyperchaotic system was realized.

As for beneficial chaotic behavior, we hope to enhance its chaotic characteristics. For example, more complex chaotic behaviors can make the stirring more sufficient in industrial process. Moreover, in secure communication, complex chaotic signal can increase the confidentiality of transmitted signal.

In order to obtain more complicated chaotic behaviors, we can extend the integer-order hyperchaotic complex systems to fractional order, because fractional order can increase the degree of freedom of the chaotic system and make its dynamic behaviors more complicated. As an extension of the integer-order complex hyperchaos system, the fractional- order complex hyperchaos system also has higher accuracy in describing processes in many fields and can more accurately describe various irregular physical phenomena. Therefore, many scholars have carried out a large amount of research on the fractional order. Ma investigated a new 4D incommensurate fractional-order chaotic system [25]. Jahanshahi investigated a multi-stable fractional-order chaotic system [26]. Rahman presented a new three-dimensional fractional-order complex chaotic system [27].

However, the above references fractional-order extension of real chaotic systems, but few scholars have studied the fractional-order hyperchaotic complex systems. In this paper, the complex hyperchaotic system was extended to fractional order.

On the basis of the above discussion, the main innovations of this paper are as follows:


The rest of this paper is structured as follows: In Section 2, we introduce the model of hyperchaotic complex system. In Section 3, we designed a controller to realize the tracking control of the hyperchaotic complex system and verified its feasibility from two aspects of a mathematical proof and simulation experiment. In Section 4, we extended the complex hyperchaotic system to fractional order and analyzed its initial value sensitivity and the system evolution process of the fractional complex hyperchaotic system with order and parameters. In the last section, we conclude the paper.

#### **2. The Model of the Hyperchaotic Complex System** = − + +

In 2021, Li et al. constructed a new hyperchaotic complex system by adding feedback term and introducing complex variables [7]. The mathematical model is as follows: = − − + = + −

$$\begin{cases}
\dot{\mathbf{x}}\_1 = a(\mathbf{x}\_2 - \mathbf{x}\_1) + \mathbf{x}\_2 \mathbf{x}\_3 + \mathbf{x}\_4 \\
\dot{\mathbf{x}}\_2 = c\mathbf{x}\_1 - \mathbf{x}\_2 - \mathbf{x}\_1 \mathbf{x}\_3 + \mathbf{x}\_4 \\
\dot{\mathbf{x}}\_3 = 1/2(\overline{\mathbf{x}}\_1 \mathbf{x}\_2 + \mathbf{x}\_1 \overline{\mathbf{x}}\_2) - b\mathbf{x}\_3 \\
\dot{\mathbf{x}}\_4 = 1/2(\overline{\mathbf{x}}\_1 \mathbf{x}\_2 + \mathbf{x}\_1 \overline{\mathbf{x}}\_2) - d\mathbf{x}\_4
\end{cases} \tag{1}$$

where . *x*<sup>1</sup> = *u*<sup>1</sup> + *ju*2, . *x*<sup>2</sup> = *u*<sup>3</sup> + *ju*<sup>4</sup> is a complex variable and . *x*<sup>3</sup> = *u*5, . *x*<sup>4</sup> = *u*<sup>6</sup> is a real variable. We separated the real part and imaginary part of the variable to obtain the equivalent mathematical model, as shown in system (2): = − + + = − + 

$$\begin{cases} \dot{u}\_1 = a(u\_3 - u\_1) + u\_3 u\_5 + u\_6 \\ \dot{u}\_2 = a(u\_4 - u\_2) + u\_4 u\_5 \\ \dot{u}\_3 = cu\_2 - u\_3 - u\_1 u\_5 + u\_6 \\ \dot{u}\_4 = cu\_2 - u\_4 - u\_2 u\_5 \\ \dot{u}\_5 = u\_1 u\_3 + u\_2 u\_4 - bu\_5 \\ \dot{u}\_6 = u\_1 u\_3 + u\_2 u\_4 - du\_6 \end{cases} \tag{2}$$

When *a* = 10, *b* = 8/3, *c* = 20, *d* = 15 and the initial value is (1, 2, 3, 4, 5, 6), the system shows obvious chaotic characteristics. The Lyapunov exponents of system (2) is as follows: *LE*<sup>1</sup> = 1.7555, *LE*<sup>2</sup> = 0.1175, *LE*<sup>3</sup> = 0, *LE*<sup>4</sup> = −11.6914, *LE*<sup>5</sup> = −14.8827, *LE*<sup>6</sup> = −22.3471, which is (+, +, 0, −, −, −), so the system is in a hyperchaotic state. The attractor phase diagram of system (2) is shown in Figure 1. <sup>1</sup> = 1.7555, <sup>2</sup> = 0.1175, <sup>3</sup> = 0, <sup>4</sup> = −11.6914, <sup>5</sup> = −14.8827, <sup>6</sup> = −22.3471 − − −

<sup>1</sup> − <sup>5</sup> <sup>3</sup> − 6 **Figure 1.** Attractor phase diagram of system (2). (**a**) *u*<sup>1</sup> − *u*<sup>5</sup> phase diagram of system (2); (**b**) *u*<sup>3</sup> − *u*<sup>6</sup> phase diagram of system (2).

#### **3. Tracking Control**

In this section, the controller is designed so that the three state variables of the system (2) can track the sine function, constant 5, and the fourth dimension of the complex Lorenz chaotic system individually.

The system model of complex Lorenz chaotic system is as follows:

$$\begin{cases}
\dot{\mathbf{x}}\_1 = a\_1(\mathbf{x}\_2 - \mathbf{x}\_1) \\
\dot{\mathbf{x}}\_2 = a\_2\mathbf{x}\_1 - \mathbf{x}\_1\mathbf{x}\_3 - a\_3\mathbf{x}\_2 \\
\dot{\mathbf{x}}\_3 = -a\_4\mathbf{x}\_3 + 1/2(\overline{\mathbf{x}}\_1\mathbf{x}\_2 - \mathbf{x}\_1\overline{\mathbf{x}}\_2)
\end{cases} \tag{3}$$

Separating the real part and imaginary part of the complex Lorenz chaotic system, the following equivalent mathematical model can be obtained:

$$\begin{cases} \dot{\mathbf{x}}\_1^r = a\_1(\mathbf{x}\_2^r - \mathbf{x}\_1^r) \\ \dot{\mathbf{x}}\_1^i = a\_1(\mathbf{x}\_2^i - \mathbf{x}\_1^i) \\ \dot{\mathbf{x}}\_2^r = a\_2\mathbf{x}\_1^r - \mathbf{x}\_1^r\mathbf{x}\_3 - a\_3\mathbf{x}\_2^r \\ \dot{\mathbf{x}}\_2^i = a\_2\mathbf{x}\_1^i - \mathbf{x}\_1^i\mathbf{x}\_3 - a\_3\mathbf{x}\_2^i \\ \dot{\mathbf{x}}\_3 = -a\_4\mathbf{x}\_3 + (\mathbf{x}\_1^r\mathbf{x}\_2^r + \mathbf{x}\_1^i\mathbf{x}\_2^i) \end{cases} \tag{4}$$

We add the controller to system (2) and obtain system (5)

$$\begin{cases} \dot{u}\_1 = a(u\_3 - u\_1) + u\_3 u\_5 + u\_6 + v\_1 \\ \dot{u}\_2 = a(u\_4 - u\_2) + u\_4 u\_5 + v\_2 \\ \dot{u}\_3 = cu\_2 - u\_3 - u\_1 u\_5 + u\_6 + v\_3 \\ \dot{u}\_4 = cu\_2 - u\_4 - u\_2 u\_5 + v\_4 \\ \dot{u}\_5 = u\_1 u\_3 + u\_2 u\_4 - bu\_5 + v\_5 \\ \dot{u}\_6 = u\_1 u\_3 + u\_2 u\_4 - du\_6 + v\_6 \end{cases} \tag{5}$$

where *V* = (*v*1, *v*2, *v*3, *v*4, *v*5, *v*6) *T* is the designed controller vector. We obtain the following Theorem 1.

**Theorem 1.** *As for system (5), if we design the controller as system (6),*

$$\begin{cases} v\_1 = \dot{r}\_1 + r\_1 + (a-1)u\_1 - au\_3 - u\_3u\_5 - u\_6 \\ v\_2 = \dot{r}\_1 + r\_1 + (a-1)u\_2 - au\_4 - u\_4u\_5 \\ v\_3 = \dot{r}\_2 + r\_2 - cu\_2 + u\_1u\_5 - u\_6 \\ v\_4 = \dot{r}\_3 + r\_3 - cu\_2 + u\_2u\_5 \\ v\_5 = \dot{r}\_2 + r\_2 - u\_1u\_3 - u\_2u\_4 + (b-1)u\_5 \\ v\_6 = \dot{r}\_2 + r\_2 - u\_1u\_3 - u\_2u\_4 + (d-1)u\_6 \end{cases} \tag{6}$$

*then the state variable u*<sup>1</sup> *can track sine function, the state variable u*<sup>2</sup> *can track constant 5, and the state variable can track x<sup>i</sup>* 2 *of the complex Lorenz chaotic system.*

**Proof:** Set the error as *e*1,*e*2,*e*3, then

.

$$e\_1 = \mu\_1 - r\_1, e\_2 = \mu\_3 - r\_2, e\_3 = \mu\_4 - r\_3 \tag{7}$$

The expected goals are expressed as follows,

$$\begin{cases} \begin{aligned} r\_1 &= \sin t \\ r\_2 &= \mathbf{5} \\ r\_3 &= \mathbf{x}\_2^i \end{aligned} \Rightarrow \begin{cases} \dot{r}\_1 = \cos t \\ \dot{r}\_2 = 0 \\ \dot{r}\_3 &= a\_2 \mathbf{x}\_1^i - a\_3 \mathbf{x}\_2^i - \mathbf{x}\_1^i \mathbf{x}\_3 \end{aligned} \tag{8}$$

Select Lyapunov function *V* = <sup>1</sup> 2 (*e*<sup>1</sup> + *e*<sup>2</sup> + *e*3) <sup>2</sup> > 0. Substituting (5)–(8) into Lyapunov function *V*, we can obtain

$$\begin{array}{l} \dot{V} &= e\_1 \dot{e}\_1 + e\_2 \dot{e}\_2 + e\_3 \dot{e}\_3 \\ &= (u\_1 - r\_1)(u\_1 - r\_1)' + (u\_3 - r\_2)(u\_3 - r\_2)' + (u\_4 - r\_3)(u\_4 - r\_3)' \\ &= (u\_1 - \sin t)(\sin t - u\_1) + (u\_3 - r\_2)(5 - u\_3) \\ &+ (u\_4 - \dot{x}\_2^i) \left(cu\_2 - u\_4 - u\_2 u\_5 + v\_4 - a\_2 \dot{x}\_1^i + a\_3 \dot{x}\_2^i + \dot{x}\_1^i \mathbf{x}\_3\right) \\ &= (u\_1 - \sin t)(\sin t - u\_1) + (u\_3 - \mathbf{5})(\mathbf{5} - u\_3) + \left(u\_4 - \dot{x}\_2^i\right) \left(\dot{x}\_2^i - u\_4\right) \\ &= -\left(u\_1 - \sin t\right)^2 - \left(u\_3 - \mathbf{5}\right)^2 - \left(u\_4 - \dot{x}\_2^i\right)^2 < 0 \end{array}$$

According to the Lyapunov stability theorem, the error of tracking control approaches 0, and the proof is completed.

= + +

= + +

= + +

( ) ( ) ( ) ( ) ( ) ( )

− − − − − −

( ) ( ) ( ) ( ) ( ) ( )

− − − − − −

( ) ( ) ( ) ( ) ( ) ( )

− − − − − −

( ) ( )

( ) ( )

( ) ( )

= − − + − − + − −

= − − + − − + − −

= − − + − − + − −

+ − − − + − + +

+ − − − + − + +

+ − − − + − + +

( ) ( ) ( ) ( ) ( ) ( )

( ) ( ) ( ) ( ) ( ) ( )

( ) ( ) ( ) ( ) ( ) ( )

–

–

–

( ) ( ) ( ) ( )

( ) ( ) ( ) ( )

( ) ( ) ( ) ( )

= − − + − −

= − − + − −

= − − + − −

( ) ( ) ( )

□

□

□

= − − − − − −

( ) ( ) ( )

( ) ( ) ( )

= − − − − − −

= − − − − − −

**Figure 2.** Tracking sine function.

**Figure 3.** Tracking constant 5.

**Figure 4.** Tracking the fourth dimension of the complex Lorenz chaotic system.

The simulation results are consistent with the theoretical analysis, and the effectiveness of the controller is verified from two aspects of simulation and mathematical analysis.

#### **4. Fractional-Order Generalization**

In order to make the beneficial chaotic behavior in industrial process more complicated, this paper extended the hyperchaotic complex system to fractional order, which is a useful and simple method to enhance the beneficial chaotic behaviors.

#### *4.1. Mathematical Background*

Fractional order refers to any order of calculus. In a sense, fractional calculus is a generalized form of integer calculus. For fractional calculus, there are three main definitions: Grunwald Letnikov definition, Riemanu Liouville definition, and Caputo definition. Since Caputo definition includes initial conditions and initial values, Caputo calculus is considered in engineering calculations. In this paper, we chose the Caputo definition.

**Definition 1.** *Caputo fractional differential form is as follows*:

$$\, \_0^\mathbb{C} D\_t^a f(t) = \frac{d^m}{dt^m} I\_{m-a} = \begin{cases} \frac{1}{\Gamma(m-a)} \int\_0^t \frac{f(m)(\tau)}{(t-\tau)^{a-m+1}} d\tau, & m-1 < a < m\\ & \frac{d^m}{dt^m} f(t), & a = m \end{cases} \tag{9}$$

*where <sup>α</sup>* = [*m*]+1, [*m*] *is the integer part of m,* <sup>Γ</sup>(∗) *is the gamma function, and <sup>D</sup><sup>α</sup> t is the gamma function α order differential operator. In this paper, D<sup>α</sup>* ∗ *is used to represent C* <sup>0</sup> *D<sup>α</sup>* ∗ *, and we mainly consider the case of* 0 < *α* < 1.

#### *4.2. Fractional-Order System Model*

In this section, the hyperchaotic complex system is extended to fractional order, and the following fractional-order new hyperchaotic complex system is constructed.

$$\begin{cases} D\_{\ast}^{\theta\_1} \mathbf{x}\_1 = a(\mathbf{x}\_2 - \mathbf{x}\_1) + \mathbf{x}\_2 \mathbf{x}\_3 + \mathbf{x}\_4 \\\ D\_{\ast}^{\theta\_2} \mathbf{x}\_2 = c\mathbf{x}\_1 - \mathbf{x}\_2 - \mathbf{x}\_1 \mathbf{x}\_3 + \mathbf{x}\_4 \\\ D\_{\ast}^{\theta\_3} \mathbf{x}\_3 = 1/2(\overline{\mathbf{x}}\_1 \mathbf{x}\_2 + \mathbf{x}\_1 \overline{\mathbf{x}}\_2) - b\mathbf{x}\_3 \\\ D\_{\ast}^{\theta\_4} \mathbf{x}\_4 = 1/2(\overline{\mathbf{x}}\_1 \mathbf{x}\_2 + \mathbf{x}\_1 \overline{\mathbf{x}}\_2) - d\mathbf{x}\_4 \end{cases} \tag{10}$$

where *D αl* <sup>∗</sup> is the Caputo operation of order *α<sup>l</sup>* , and *α<sup>l</sup>* is the order of relevant variables of *xl*(*l* = 1, 2, 3, 4). According to fractional linear operation, we can obtain *D α*1 <sup>∗</sup> *x*<sup>1</sup> = *D α*1 <sup>∗</sup> (*u*<sup>1</sup> + *ju*2) = *D α*1 <sup>∗</sup> *u*<sup>1</sup> + *jD α*1 <sup>∗</sup> *u*2, *D α*2 <sup>∗</sup> *x*<sup>2</sup> = *D α*2 <sup>∗</sup> (*u*<sup>3</sup> + *ju*4) = *D α*2 <sup>∗</sup> *u*<sup>3</sup> + *jD α*2 <sup>∗</sup> *u*<sup>4</sup> *D α*3 <sup>∗</sup> *x*<sup>3</sup> = *D α*3 <sup>∗</sup> *u*5, *D α*4 <sup>∗</sup> *x*<sup>4</sup> = *D α*4 <sup>∗</sup> *u*6.

The above system can be transformed into the following forms:

$$\begin{cases} D\_{\*}^{\underline{\alpha}}u\_1 = a(u\_3 - u\_1) + u\_3u\_5 + u\_6 \\ D\_{\*}^{\underline{\alpha}}u\_2 = a(u\_4 - u\_2) + u\_4u\_5 \\ D\_{\*}^{\underline{\alpha}}u\_3 = cu\_1 - u\_3 - u\_1u\_5 + u\_6 \\ D\_{\*}^{\underline{\alpha}}u\_4 = cu\_2 - u\_4 - u\_2u\_5 \\ D\_{\*}^{\underline{\alpha}}u\_5 = u\_1u\_3 + u\_2u\_4 - bu\_5 \\ D\_{\*}^{\underline{\alpha}}u\_6 = u\_1u\_3 + u\_2u\_4 - du\_6 \end{cases} \tag{11}$$

= 10, = 8/3, =

#### *4.3. Fractional-Order Attractor*

We select the initial value of system (11) as (1, 1, 1, 1, 1, 1), *a* = 10, *b* = 8/3, *c* = 30, *d* = 12, and fractional-order *α<sup>l</sup>* = 0.95(*l* = 1, 2, 3, 4). On the basis of the definition of Caputo, the system (11) is simulated by MATLAB, and the attractor phase diagram of system (11) is obtained as shown in Figure 5. It can be seen that the attractor of the system presents obvious chaotic characteristics. Comparing Figure 5 with Figure 1, the interval and shape of the attractor are found to be different. 30, = 12 = 0.95( = 1,2,3,4)

<sup>1</sup> − <sup>5</sup> <sup>3</sup> − 6 **Figure 5.** Attractor phase diagram of system (11). (**a**) *u*<sup>1</sup> − *u*<sup>5</sup> phase diagram of system (11); (**b**) *u*<sup>3</sup> − *u*<sup>6</sup> phase diagram of system (11).

"0 1 test"

"0 1"test.

#### *4.4. 0-1 Test*

6

Gottwald and Melbourne proposed a reliable and effective binary test method to test whether the system is chaotic, which is called the "0-1 test" [28]. The basic idea is to establish a stochastic dynamic process for data and then study how the scale of the stochastic process changes with time. Next, we used this method to test and analyze the chaotic characteristics of system (11), as shown in Figure 6.

"0 1 test"

= 10, = 8/3, =

<sup>1</sup> − <sup>5</sup> <sup>3</sup> −

30, = 12 = 0.95( = 1,2,3,4)

"0 1"test. **Figure 6.** The "0-1"test.

It can be seen from Figure 6 that the new hyperchaotic complex system shows obvious unbounded motion, similar to Brownian motion. Therefore, it is chaotic.

#### *4.5. Order α Impact on System Status*

In this section, we all took *u*<sup>1</sup> − *u*<sup>5</sup> to observe the evolution of the attractor of the system (11).

When *α* < 0.82, the system is in a divergent state.

When *α* ∈ (0.82, 0.94), the system converges to a stable point.

When *α* = 0.95, the system changes from stable point to chaotic state, and at this time, it presents obvious butterfly attractor shape.

When *α* = 0.99, the system is still in chaos, but under the same number of cycles, the shape of the attractor is fuller.

When *α* = 1.01, the order of the system is a fractional order greater than 1. At this time, the system is still in a chaotic state, but the attractor forms are different and sparse.

When *α* = 1.03, the system is in a chaotic state, and the attractor shape becomes fuller with the increase in order.

When *α* = 1.05, the system is in period doubling limit cycle state.

When *α* > 1.07, system divergence occurs.

The detailed evolution process is shown in Figure 7. Through observation, it is found that the fractional complex hyperchaotic system shows a more complex system evolution process with the change of order, such as the position and shape of the attractor having changed to some extent.

*Order α Impact on* 

 α < 0.82 α ∈ (0.82,0.94) α = 0.95 α = 0.99 α = 1.01 α = 1.03 α = 1.05 α > 1.07

<sup>1</sup> − <sup>5</sup>

l(l=1,2,3,4) = 0.82 l(l=1,2,3,4) = 0.94 l(l=1,2,3,4) = 0.95 l(l=1,2,3,4) = 0.99 l(l=1,2,3,4) = 1.01 l(l=1,2,3,4) = 1.03 l(l=1,2,3,4) = 1.05 l(l=1,2,3,4) = 1.07 **Figure 7.** The evolution process of system (11) with order. (**a**) *α*l(l=1,2,3,4) = 0.82; (**b**) *α*l(l=1,2,3,4) = 0.94; (**c**) *α*l(l=1,2,3,4) = 0.95; (**d**) *α*l(l=1,2,3,4) = 0.99; (**e**) *α*l(l=1,2,3,4) = 1.01; (**f**) *α*l(l=1,2,3,4) = 1.03; (**g**) *α*l(l=1,2,3,4) = 1.05; (**h**) *α*l(l=1,2,3,4) = 1.07.

= 0.95

#### *4.6. Influence of Parameter Change on System Attractor*

 = 0.95 <sup>1</sup> − <sup>5</sup> It can be seen from Section 4.4 that when order *α* = 0.95, the system attractor presents an obvious chaotic attractor form. In this section, we selected order *α* = 0.95 and changed the values of system parameters *a*, *b*, *c*, and *d* individually to observe the influence of parameters on the system state and system attractor. In this section, we all took the phase diagram of *u*<sup>1</sup> − *u*<sup>5</sup> to observe the evolution of the attractor of the system.

#### 4.6.1. Parameter *a* Change

We kept parameters *b*, *c*, and *d* unchanged; changed the value of parameter *a*; and selected the same initial value as the integer order to observe the evolution process of the system with parameter *a*, as shown in Figure 8. It can be seen that the system entered the chaotic state from the limit cycle state and continued to evolve in the chaotic state. The attractor changed from sparse to full and then to sparse, and then returned to the limit cycle state and finally diverged.

**Figure 8.** The evolution process of system (11) with parameter *a*. (**a**) *a* = 6, limit cycle; (**b**) *a* = 8, chaos; (**c**) *a* = 10, chaos; (**d**) *a* = 20, chaos; (**e**) *a* = 23, chaos; (**f**) *a* = 28, limit cycle.

#### 4.6.2. Parameter *b* Change

We kept the parameters *a*, *c,* and *d* unchanged; changed the value of parameter *b*; and selected the same initial value as the integer order to observe the evolution process of the system with the change of parameter *b*, as shown in Figure 9. It can be seen that the system entered the chaotic state from the limit cycle state and continued to evolve in the chaotic state. The attractor changed from sparse to full, then returned to the limit cycle state and finally diverged.

**Figure 9.** The evolution process of system (11) with parameter *b*. (**a**) *b* = 0.2, limit cycle; (**b**) *b* = 0.5, chaos; (**c**) *b* = 1.5, chaos; (**d**) *b* = 8/3, chaos; (**e**) *b* = 3, chaos; (**f**) *b* = 4, limit cycle.

#### 4.6.3. Parameter *c* Change

We kept parameters *a*, *b,* and *d* unchanged; changed the value of parameter *c*; and selected the same initial value as the integer order to observe the evolution process of the system with the change of parameter *c*, as shown in Figure 10. It can be seen that the system entered the chaotic state from the limit cycle state, evolved continuously in the chaotic state, and finally returned to the limit cycle state and diverged.

**Figure 10.** The evolution process of system (11) with parameter *c*. (**a**) *c* = 22, limit cycle; (**b**) *c* = 23, chaos; (**c**) *c* = 30, chaos; (**d**) *c* = 40, chaos; (**e**) *c* = 44, chaos; (**f**) *c* = 44.8, limit cycle.

#### 4.6.4. Parameter *d* Change

We kept parameters *a*, *b*, and *c* unchanged; changed the value of parameter *d*; and selected the same initial value as the integer order to observe the evolution process of the system with parameter *d*, as shown in Figure 11. It can be seen that the system entered the chaotic state from the limit cycle state and evolved continuously in the chaotic state. It can be seen that when the system was in the chaotic state, compared with parameters *a*, *b*, *c*, parameter *d* had the largest value range, and finally the system returned to the limit cycle state and finally diverged.

**Figure 11.** The evolution process of system (11) with parameter *d*. (**a**) *d* = 6, limit cycle; (**b**) *d* = 9, chaos; (**c**) *d* = 10, chaos; (**d**) *d* = 60, chaos; (**e**) *d* = 100, chaos; (**f**) *d* = 151.6, limit cycle.

It was found that the fractional hyperchaotic complex system showed a more complex system evolution process with the change of parameters, such as the position and shape of the attractor having changed to some extent.

#### **5. Conclusions**

In this paper, aiming at the harmful hyperchaotic complex behavior in industrial process, a tracking controller was designed for the hyperchaotic complex system so that the three state variables of the hyperchaotic complex system can track the controller of sine function, constant 5, and complex Lorenz chaotic system individually, and its stability was proven to realize the tracking control of the system. In order to make the intentional chaotic behavior in the industrial process more complex, we extended the hyperchaotic complex system to fractional order. The effects of initial value, order, and parameter changes on the fractional hyperchaotic complex system were discussed and studied. When the order and parameters changed, the detailed evolution process of the system state was given. It was found that there were no coexistence attractors and parameter attractors in the system.

In this paper, the application of the hyperchaotic complex system was studied. For the harmful chaotic system, the controller was designed to convert it into the desired system; for the beneficial chaotic system, this paper extended it to fractional order, which made its chaotic behavior more complex.

There are several prospects for the study of chaos theory: (1) research on the physical background of the chaotic system, or it can show more abundant dynamic behavior; (2) chaos theory can be applied to some complex systems, such as weather forecasting and industrial processes; (3) on the basis of chaotic systems, new chaotic cryptographic algorithms or chaotic neural networks can be formed and applied in various fields.

**Author Contributions:** Conceptualization, F.L., Z.L. and F.Z.; methodology, Z.L. and L.L.; software, F.L., Z.L. and S.Z.; validation, F.L., L.L. and Z.L.; formal analysis, F.L. and L.L.; investigation, F.L. and F.Z.; resources, F.L. and F.Z.; data curation, F.L. and Z.L.; writing—original draft preparation, F.L., L.L. and Z.L.; writing—review and editing, F.L., Z.L. and F.Z.; visualization, F.L., L.L., Z.L. and S.Z.; supervision, F.Z.; project administration, F.Z.; funding acquisition, F.L., Z.L. and F.Z. All authors have read and agreed to the published version of the manuscript.

**Funding:** This work is in part supported by the Shandong Provincial Natural Science Fund (2022HWY Q-081), the National Science Foundation of China Key Project (U21A20488), Academic Promotion Project of Shandong First Medical University, the International Collaborative Research Project of Qilu University of Technology (QLUTGJHZ2018020), and major scientific and technological innovation projects of Shandong Province under grant number (2019JZZY010731 and 2020CXGC010901).

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Not applicable.

**Conflicts of Interest:** The authors declare no conflict of interest.

#### **References**

