Abstract

In this paper, a spectral entropy algorithm has been used successfully for complexity measure of a discrete permanent-magnet synchronous motor (PMSM) system. Firstly, the discrete PMSM system is achieved using the forward Euler scheme. Secondly, by adopting the bifurcation diagram, phase portraits, 0-1 test, and largest Lyapunov exponent, the chaotic dynamics of the discrete PMSM system are analyzed. The complexity of the discrete PMSM system is discussed by employing the spectral entropy algorithm. It shows that the spectral entropy complexity analysis is an efficient tool to study chaotic dynamics. Finally, we illustrate this result through numerical experiments.

1. Introduction

Due to its potential applications in the fields of secure communication and information encryption, chaos as a hot topic has been widely studied for nearly three decades, especially the chaotic dynamics and complexity analysis of chaotic sequences.

At present, the complexity measure of chaotic sequences can be analyzed by behavioral complexity and structure complexity. As we have known, there are several methods of measuring complexity of chaotic sequences, including the Lempl-Ziv algorithm [1, 2], ApEn algorithm [35], FuzzyEn algorithm [6, 7], LMC algorithm [810], SymEn algorithm [11], sample entropy algorithm [5], permutation entropy (PE) algorithm [12], entropy algorithm [1316], SCM algorithm [1618], and C0 algorithm [15, 1922], which are widely used to analyze chaotic sequences. We should point out that the spectral entropy algorithm is based on the power spectrum of the chaotic signal and describes the irregularity of the signal spectrum, and the computation cost of this method is low [16]. Therefore, the entropy algorithm is an efficient tool to deal with complexity of chaotic sequences.

In recent years, a discrete chaotic system, which is obtained by the forward Euler scheme for its continuous system, is widely studied. A natural question to follow is whether these discrete chaotic systems are more complicated than its continuous systems or not? How to select the system parameters and step size to get more complex dynamics? Therefore, it is worthwhile to take the fluctuation of step size into full account for the dynamic problem of a chaotic system. To our knowledge, the relevant researches of these questions keep open, and these topics deserve further investigation.

Based on the above questions, the major task of this paper is to deal with complexity of chaotic sequences by means of the spectral entropy algorithm. As a result, a discrete permanent-magnet synchronous motor (PMSM) system, which is obtained by a continuous PMSM system, is considered. The PMSM system is a kind of high-efficient and high-powered motor, which has been widely used in the industry [23]. The chaotic dynamics of the continuous PMSM system becomes one of most active research areas. Many important research works on the PMSM system are found in [2427] and the references cited therein.

This paper is organized as follows. In Section 2, the chaotic dynamics of the discrete PMSM system, which is obtained by using the forward Euler scheme, is investigated by means of the bifurcation diagram, phase portraits, 0-1 test, and largest Lyapunov exponent. In Section 3, the spectral entropy algorithm is described, and the complexity of chaotic sequences in the discrete PMSM system is calculated. Finally, conclusions and future direction are presented in Section 4.

2. Dynamics of the Three-Dimensional Discrete PMSM System

2.1. The Three-Dimensional Discrete PMSM System

This section illustrates the application of the forward Euler method to a continuous PMSM system.

Recently, the following three-dimensional continuous PMSM system [2327] is proposed. where and are positive parameters. When and , the PMSM system is chaotic. Its typical butterfly chaotic attractor is shown in Figure 1. When the forward Euler is employed, we obtain the following three-dimensional discrete PMSM system: where and are positive parameters and is the step size.

2.2. Fixed Points and Their Stability

The fixed points , , and of system (2) are obtained by solving the equations:

In order to study the stability of the fixed point , we calculate the Jacobian matrix of system (2) at as and its characteristic polynomial equation is where , , and . According to the Schur-Cohn criterion [28], we derive the stability conditions of as follows: (a) and (b) . As a result, by taking and , the stability condition of is . Here, we only consider stability of the fixed point , and the stability of the fixed points and can be studied in a similar manner as for .

2.3. Chaos in the Three-Dimensional Discrete PMSM System

In this section, by means of the phase portraits, bifurcation diagram, 0-1 test and largest Lyapunov exponent, the dynamics of system (2) is researched.

In particular, variations of the parameter are considered by keeping and . The largest Lyapunov exponent diagram of system (2) is obtained using the Jacobi matrix method as shown in Figure 2. Letting and fixing , the bifurcation diagram of system (2) in Figure 3 is obtained for in the range . From Figures 2 and 3, of particular interest is that it has a chaotic attractor with one positive Lyapunov exponent in a large range of parameters. As a result, starting with , a stable solution is obtained, and the corresponding phase diagram is given in Figure 4(a). By increasing the value of , a two-scroll chaotic attractor is obtained for as shown in Figure 4(b). By adopting 0-1 test [29], the corresponding trajectories in the plane are shown in Figures 5(a) and 5(b). Here, and are two real valued functions, respectively, where , , and is a random number in the range . Also note that bounded trajectories in the plane correspond to regular dynamics, whereas Brownianlike (unbounded) trajectories correspond to chaotic dynamics [29]. See Gottwald and Melbourne [29] for details.

From Figures 24, it is clear that the stable, periodic solutions and ensuing transition to chaos are observed.

3. Complexity Analysis

In this section, the complexity of system (2) is discussed by employing the spectral entropy algorithm.

The spectral entropy algorithm is a method that deals with the complexity of a time series. The spectral entropy reflects the disorder in the Fourier transformation domain, if a flatter spectrum has a larger value of spectral entropy, which shows a higher complexity of the time series [16]. Consider a set of discrete data with a length of , representing a one-dimensional observable data set. According to discrete Fourier transformation approach, the spectral entropy algorithm [1316] is briefly described as follows.

Step 1. Data processing. In order to better reflect the energy information of the data , transform the data as follows: where .

Step 2. Discrete Fourier transformation. Taking the discrete Fourier transform for the data and its corresponding discrete Fourier transform is defined as where and is the unit imaginary.

Step 3. Calculating relative power spectrum. If the power of a discrete power spectrum with the th frequency is , then the “probability” of this frequency is defined as

Step 4. Spectral entropy is calculated. Applying the , the spectral entropy is defined by

Based on the above four steps, the spectral entropy of a given time series with length is calculated, without choosing any other parameters [15].

Here, we choose parameter planes , , and to calculate its complexity, for example, to see how the complexity changes as the parameters , , and increase (or decrease). As a result, the length of data for spectral entropy is 4 × 104 after removing the first 1 × 104 points of data [15]. By taking , the spectral entropy complexity of series of system (2) is calculated and illustrated in Figure 6.

From Figure 6, it can be seen that the parameter has a larger range of choices when parameter is smaller. For example, the parameter is selected in the range ; the system (2) maintains a high complexity for parameter . When the parameters and are selected in the two light-colored areas in the upper right, the system (2) maintains the largest complexity. Hence, for this chaotic system (2), the complexity depends on parameters and sensitively.

Next, we calculate the complexity of system (2), when the step size is selected in the range . By taking , the spectral entropy complexity in Figure 7 is obtained for parameters and in the range and , respectively.

From Figure 7, it can be seen that the largest value of spectral entropy complexity of system (2) fluctuates within a small range for parameters and . Similarly, by taking , the spectral entropy complexity in Figure 8 is obtained for parameters and in the range and , respectively. According to Figure 8, it can be seen that the largest value of spectral entropy complexity of system (2) is 0.7. For parameters and , they are roughly in the range of and , respectively.

Based on the discussion above, we found that spectral entropy complexity analysis is a convenient method to choose parameters for a discrete system.

In the following discussion, two of the parameters are fixed and the spectral entropy complexity is calculated. By taking and , the spectral entropy complexity in Figure 9(a) is obtained for parameter in the range . Similarly, fixing and , the spectral entropy complexity is calculated for parameter in the range as shown in Figure 9(b). Combined with Figure 9, we find that the spectral entropy complexity calculation is similar to that of the Lyapunov exponent, which reflects the complexity of chaos. In Figure 9(a), letting , , and , the 0-1 test is obtained as shown in Figure 10. In Figure 9(b), letting , , and , the 0-1 test is obtained as shown in Figure 11.

In addition, it is worth noting that the spectral entropy complexity of system (2) increases with the increasing step size . In particular, variations of the step size is considered by keeping and . The spectral entropy complexity is shown in Figure 12.

4. Conclusions

In this paper, by means of the largest Lyapunov exponent, phase diagram, bifurcation diagram, 0-1 test, and spectral entropy algorithm, we investigated the chaotic dynamics and complexity of the discrete PMSM system obtained by the forward Euler scheme. The conclusions are drawn as follows.

The discrete PMSM system contains rich dynamical behaviors. The two-scroll stronger chaotic attractor is observed with the increasing parameter . Moreover, the discrete PMSM system is chaotic in large range for parameter . Furthermore, by calculating the chaos complexity of discrete PMSM system, we found that it has a high complexity. Therefore, the discrete PMSM system has higher application value in the real-word applications. For example, this system can be considered for chaotic encryption.

We also found an interesting phenomenon that the spectral entropy complexity of discrete PMSM system increases with increasing step size . Thus, the spectral entropy complexity analysis is a convenient method to choose step size. In fact, it also provides a reference for Euler method to study the chaotic dynamics of continuous dynamical system.

It is expected that the results enrich our knowledge of the dynamics of discrete systems, which is obtained by the forward Euler scheme. Furthermore, it suggests that the results are beneficial to researchers who explore other dynamical behaviors.

Data Availability

No additional data are available for the paper.

Conflicts of Interest

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work is supported by the Guizhou Science and Technology Cooperation Plan (Qian Ke He LH zi [2016]7032 and Qian Ke He LH zi [2016]7031) and Youth Science and Technology Talents Growth Projects of Department of Education of Guizhou Province (no. 12847, Qian Jiao He KY Zi [2018]315, [2018]313, [2017]256, and [2017]251 and Qian Ke He LH Zi [2017]7072).