Initial Value Determination of Chua System with Hidden Attractors and Its DSP Implementation

In this paper, a method for determining the initial value of the hidden attractors in the Chua system is studied. *e initial value of the hidden attractors can be calculated quickly and accurately by the proposed method, and the hidden attractors can be found by numerical simulation. *en, the initial values of the hidden attractors are set accurately by digital signal processor (DSP), so as to the circuit realization of the chaotic system with hidden attractors is performed. *e results show that the numerical simulation results of Matlab are consistent with the experimental results of DSP.

Because the domain of attraction of the hidden attractor does not intersect with any small neighbourhood of the equilibrium point, there is no general method to predict the existence of the hidden attractor, so it is of great theoretical and practical significance to study the hidden attractor in the field of machinery and so on [24,25]. In 2011, Leonov et al. proposed a locating algorithm for hidden attractors [26] and used the algorithm to find hidden attractors of the Chua system. Since then, the research on hidden attractors has attracted extensive interest of scholars. In 2012, Leonov et al. used the algorithm to find the hidden attractor in the Chua system with hyperbolic tangent function as nonlinear function [27]. In 2014, Zhao et al. used the algorithm to find the hidden attractor in a generalized autonomous Pol-Duffing system [28]. In the same year, Li et al. found the twin hidden attractors in the Chua system [29]. In 2017, Zhao et al. found the hidden attractor in a modified Chua system [30]. In the same year, Kuznetsov et al. also found coexistence limit cycle and symmetric hidden attractors in the Chua system [31]. Stankevich et al. analysed the scenario of the birth of the hidden attractor from its attractor basin in the Chua system [32]. ese attractors are different from the classical Lorenz attractors, Chua attractors, and Chen attractors.
ey are not near the equilibrium point and cannot be calculated by traditional methods. Leonov et al. proposed an algorithm to determine the initial value of the hidden attractors and found the hidden attractors.
In 2016, Bao et al. designed the chaotic circuit of the chaotic system and found the hidden attractor of the system by PSIM simulation [33]. In [34], a three-dimensional autonomous chaotic circuit is designed, and the hidden attractor of the system is found by PSpice simulation. In [35], a 5-D memristor chaotic circuit is designed. e hidden multiscroll attractors and hidden multiwing attractors are found by PSpice simulation. In [36], a 5D extreme multistable chaotic circuit is designed, and the hidden attractor of the system is found by PSpice simulation. In [37], a fourdimensional chaotic circuit is designed, and the coexistence hidden attractor of the system is found by PSpice simulation. In [38], a new chaotic circuit is designed by introducing the cosh function into the system in [37], and the coexistence hidden attractor of the system is found by PSpice simulation. In [39], by introducing a memristor to improve Fitzhugh-Nagumo circuit, a new memristor chaotic circuit is designed, and the hidden attractor is found by PSIM simulation. However, the hidden attractor is studied by circuit simulation software, but the initial state of the experimental circuit is random, so the initial value of the hidden attractor cannot be set accurately.
In this paper, we study the method to determine the initial value of the hidden attractors in the Chua system. Its initial value of the hidden attractors can be set accurately by DSP, and the circuit realization of the chaotic system with hidden attractors is performed. e results show that the numerical simulation results of Matlab are consistent with the experimental results of DSP. e rest of this work is organized as follows. Section 2 describes initial value determining algorithm for the Chua system with hidden attractors. Section 3 calculates initial values of hidden attractors and finds its hidden attractors.
e Chua system with hidden attractors is implemented by DSP in Section 4. Finally, we conclude in Section 5.

Initial Value Determining Algorithm for Hidden Attractors
According to the initial value determining algorithm for the chaotic system with hidden attractors in [26], the Chua system with hidden attractors is x, y, and z are system variables, and a, b, c, m, and n are system constants. Now, system (1) is rewritten as lure system: where P � Let k be the coefficient of harmonic linearization, and ε be an infinitesimal number, and equation (2) can be rewritten as where Using nonsingular linear transformation x � Sy, equation (3) can be transformed as where e transfer function of equation (4) can be expressed as e transfer functions of system (3) can be expressed as where p is complex variables, ω 0 is the initial frequency, which can be calculated by ImWH(ω 0 )50, and k is the harmonic linearization coefficient, which can be calculated by k � -(ReW H iω 0 ) -1 . From the equivalence of the transfer functions of systems (3) and (4), it can be concluded: System (3) is transformed by nonsingular linear transformation, and it can be concluded:  (1 + a − d))/a).

Complexity
For the infinitesimal number ε, the initial value of (4) is From equation (9), the relationship between the initial values of equations (3) and (4) can be obtained: In this way, the initial value of system (1) is where the description function of a 0 can be calculated as and the description function satisfies Φ(g 0 ) � 0 and b 1 (dΦ(g)/dg)| g�g 0 ≠ 0.

Numerical Simulation
According to the algorithm in [26], the initial values are Based on the calculated initial value, the phase diagram is shown in Figure 1, and the attractor basin is shown in Figure 2. From Figure 1, it can be seen that the hidden attractor can be found according to the initial value calculated by the algorithm in [26]. From Figure 2, the blue centre region is stable equilibrium, the red region is period-1 limit cycle, and the cyan region is divergent.

DSP Implementation of Chua System with
Hidden Attractors e realization of the chaotic system by hardware circuit is the most common method to verify the design of new chaotic system, including analog circuit and digital circuit. Analog circuit mainly adopts discrete components [13] and integrates circuit (IC) [40][41][42] design method, while digital circuit mainly adopts FPGA [43] and DSP [44][45][46]. It is difficult to design and debug chaotic circuit with discrete components, and the circuit is bulky. When using IC to design the chaotic oscillator, the chip area is greatly reduced, but the design of IC requires high chip technology, and the number of wings or scrolls of the attractor is difficult to control. Because the analog circuit cannot accurately set the initial state of the system and cannot reach the calculated initial value of the hidden attractors. FPGA and DSP have high-speed data processing capability and can realize various processing algorithms through software programming [47,48], which can conveniently realize the nonlinear characteristics of the chaotic system with hidden attractors.

Implementation of the Chaotic System.
In this part, the chaotic system with hidden attractors is implemented on DSP platform. e block diagram of the working principle is shown in Figure 3. In the experiments, the Texas Instrument DSP device TMS320F28335 is employed. It is a 32 bit DSP running at 150 MHz with floating point operations. Such a high-speed clock rate is considered to be sufficient. In order to observe the phase diagrams of the attractor on the oscilloscope, the digital chaotic sequences generated on DSP are converted into analog signals. DAC8552, a 16 bit digitalto-analog converter with dual channels, is adopted. It connects DSP through SPI (serial peripheral interface). e flowchart of the programming is shown in Figure 4. In the program, in order to reduce the effect of finite computing precision in digital circuits, all data types are defined as long double. After initializing DSP, we set the initial conditions, including initial values of state variables, and the system parameters. Iterative computation is started according to the initial values [x 0 , y 0 , z 0 ]. To keep the iteration not being affected by data processing, it is necessary to push the results of each iteration into the stack. Data processing includes two steps. Firstly, an appropriate positive number is added to all data to make sure all data is greater than zero. Secondly, all data is rescaled and truncated to make the output adapting the 16 bit digital-to-analog converter.

Runge-Kutta4 (RK4) Algorithm.
According to the required iterative equation in Figure 4, we use the RK4 algorithm to realize the iterative equation of the Chua system with hidden attractors. RK4 is a derivative of Runge-Kutta basic model, which is used to solve ordinary differential equations with high accuracy, and mostly has proved itself superior to other solutions. RK4 algorithm is expressed as According to the values of h, x i , and y i and calculating the values of k 1 , k 2 , k 3 , and k 4 , we can attain the value of y i+1 . ree equations of system (1) are substituted into equation (13), and the three state variables (x, y, z) of system (1) are solved, respectively:        Figure 5 shows the hardware part of DSP. Phase diagrams of the system are captured randomly by the oscilloscope, as shown in Figure 6. When a � 8.4, b � 12.1, c � 0.005, m � -1.1, and n � 0.1, its phase diagrams is shown in Figure 7. It indicates that the Chua system with hidden attractors is realized successfully on the DSP platform. From Figures 5-7, it is observed that the DSP circuit can generate two hidden attractors.

Conclusions
In this paper, we calculate the initial values of the Chua system with hidden attractors, find its hidden attractors, and obtain its phase diagram and attractor basin. Since the analog circuit cannot accurately set its initial state and cannot achieve the calculated initial value of the hidden attractor, this paper uses DSP to realize the chaos system with hidden attractors. e results show that the numerical simulation is consistent with the experimental results of DSP, which provides a practical method for the circuit implementation of the chaotic system with hidden attractors. In the next work, we will study the hidden attractor applied to secure communication.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that they have no conflicts of interest.

Authors' Contributions
e work presented in this paper was a collaboration of all authors. Xianming WU contributed the idea and wrote the paper. Weijie Tan and Huihai Wang did the simulation analysis and reviewed the paper.