Yanan Xing,Wenjie Dong,Jian Zeng,Pengteng Guo,Jing Zhang,Qun Ding
Study of Generalized Chaotic Synchronization Method Incorporating Error-Feedback Coefficients
chaotic synchronization,generalized synchronization,chaotic hiding and anti-hiding,parameter control,transmission system
In this paper, taking the generalized synchronization problem of discrete chaotic systems as a starting point, a generalized synchronization method incorporating error-feedback coefficients into the controller based on the generalized chaos synchronization theory and stability theorem for nonlinear systems is proposed. Two discrete chaotic systems with different dimensions are constructed in this paper, the dynamics of the proposed systems are analyzed, and finally, the phase diagrams, Lyapunov exponent diagrams, and bifurcation diagrams of these are shown and described. The experimental results show that the design of the adaptive generalized synchronization system is achievable in cases in which the error-feedback coefficient satisfies certain conditions. Finally, a chaotic hiding image encryption transmission system based on a generalized synchronization approach is proposed, in which an error-feedback coefficient is introduced into the controller.
Study of Generalized Chaotic Synchronization Method Incorporating Error-Feedback Coefficients

In this paper, taking the generalized synchronization problem of discrete chaotic systems as a starting point, a generalized synchronization method incorporating error-feedback coefficients into the controller based on the generalized chaos synchronization theory and stability theorem for nonlinear systems is proposed. Two discrete chaotic systems with different dimensions are constructed in this paper, the dynamics of the proposed systems are analyzed, and finally, the phase diagrams, Lyapunov exponent diagrams, and bifurcation diagrams of these are shown and described. The experimental results show that the design of the adaptive generalized synchronization system is achievable in cases in which the error-feedback coefficient satisfies certain conditions. Finally, a chaotic hiding image encryption transmission system based on a generalized synchronization approach is proposed, in which an error-feedback coefficient is introduced into the controller.

Chaos is a unique nonlinear dynamical phenomenon with the properties of ergodicity, initial sensitivity, and the long-term unpredictability of motion trajectories [1,2,3,4]. In recent years, the study of chaos has become very popular, and it is widely used in the field of secure communication [5,6,7]. Chaos control and synchronization theory, which has great potential for application in the field of chaos research, has also become a hot spot in the high-tech competition between countries [8,9]. From the point of view of chaotic system interactions, studies related to chaotic synchronization can be divided into the following categories: generalized synchronization, phase synchronization, hysteresis synchronization, and so on [10,11,12,13]. In addition, during the process of research, researchers have proposed complete synchronization, projective synchronization, and adaptive synchronization [14,15,16,17]. In practical studies, the problem of parameter selection is inevitable regarding the structural differences between the drive and response systems. The generalized synchronization problem for chaotic or hyperchaotic systems would be a more relevant and worthwhile approach, given that the problems mentioned above can be easily solved for generalized chaotic synchronization systems. Meanwhile, the development of generalized synchronization theory has provided new tools for constructing more secure communication systems. Generalized synchronization is the gradual convergence of the trajectory curves of two chaotic systems to a time-independent transformation relationship over time; that is, a functional relationship is determined between the state of the driven system and the state of the responding system, and the synchronization of the driven and responding systems is achieved by this functional relationship, which can be deterministic or nondeterministic [18,19,20,21,22]. This paper proposes a generalized chaotic synchronization method incorporating error-feedback coefficients into the process of determining the function relationships, which is based on the principle of using the relationships between the functions in the designed controller to synchronize the drive and response systems. Most of the systems used in practical engineering are high-dimensional nonlinear systems. With the continuous research in applied mathematical theory and the rapid development of computer technology, lower-dimensional chaotic systems in practical applications are facing more challenges; hence, high-dimensional hyperchaos with more than two Lyapunov exponents is of significant interest. Based on the above, new 3D and 6D discrete chaotic systems are constructed and proposed in this paper. The constructed new high-dimensional chaotic systems are used as the driving systems, and the response system is constructed by the proposed generalized chaotic synchronization method that incorporates error-feedback coefficients. The effectiveness of the synchronization method was confirmed by experiment. The paper is organized as follows: The generalized synchronization theoretic of discrete chaotic systems and the stability principle of error systems are analyzed in Section 2. A new 3D discrete chaotic system is proposed in Section 3, in which the dynamic behavior of phase diagrams, Lyapunov exponent diagrams, and bifurcation diagrams are depicted and analyzed. Subsequently, a generalized chaotic synchronization method incorporating error-feedback coefficients is proposed, with a new 3D discrete chaotic system as the driving system. The effectiveness of the method was verified by experimental simulations. In Section 4, a new 6D discrete chaotic system is proposed, and the dynamic behavioral properties of its phase diagram, Lyapunov exponent diagram, and bifurcation diagram are analyzed. Then, the new 6D discrete chaotic system is applied as the driving system through the proposed generalized chaotic synchronization method incorporating error-feedback coefficients; the effectiveness of the method was further demonstrated by performing simulations. In Section 5, a digital image transmission system based on 6D chaotic synchronization and encryption is proposed, the encryption and decryption processes are analyzed in detail, and encryption and decryption simulations are given. Then, in Section 6, security analyses are carried out based on the previously proposed encrypted image transmission system. Finally, the conclusion is given in the last section.

In our study of chaotic control problems, it is more important to convert the problem of chaotic synchronization into the analysis of system errors. The main idea is to consider the difference in the state between the drive and response systems, that is, the synchronization error of the system. Once a reasonable controller has been designed by parameter changes to make the system error asymptotically stable at the origin point, then the two systems can be considered synchronized with each other. Firstly, the mathematical model of generalized chaotic synchronization is proposed in this paper and described, as follows. Definition 1. Consider two n-dimensional nonlinear dynamical systems, and describe them using the following equations:where , , and as well as are n-dimensional nonlinear functions, and is an n-dimensional input control function. If the selectable function is applied such that , and thus , then it can be translated into the study of the error system , for which , and therefore . In this case, the drive system and response system can reach a generalized synchronization. Theorem 1. Define an invertible transform ();consequently, there is an incorporated error-feedback coefficient .Where the feedback coefficient satisfies the condition , it can enable the progressive stability of the zero solution of Error Equation (3) of the system, which is represented as follows: Because the zero solution of Equation (3) is gradually stable, by introducing a reasonable feedback coefficient (), the drive and response systems can be synchronized in a universal way. Hence, according to Theorem 1, it can be concluded that to synchronize the drive system (1) and response system (2), a nonlinear error system () needs to be constructed, and the progressive stability of the error equation of the system () needs to be guaranteed. Based on the above, the next major concern is to determine that the system error equation is asymptotically stable at the original point; therefore, the following lemmas are given: Lemma 1 ([23]). Given a linear discrete system, which can be defined as follows: where A is a coefficient matrix, and , we can draw the following conclusions: (1) Chaotic system (4) is progressively stable if the modulus of all eigenvalues of matrix A is not more than 1; (2) In case there is a matrix (Q > 0), so that the Lyapunov equation (ATPA − P = −Q) has a unique positive solution (P), system (4) is asymptotically stable. Proof of Lemma 1(1). Set , and then the tiny variables of can be calculated as follows: Because all the eigenvalues of matrix A have a value of modulo less than 1, all the eigenvalues of matrix are integers, which are less than 1; therefore, , and system (4) is asymptotically stable. □ Proof of Lemma 1(2). Set , where P is a positive definite matrix, given that , and then the tiny variables of can be calculated as follows: Furthermore, ; hence, system (4) is asymptotically stable. □ Based on the proof processes for the stability of linear discrete systems as related in Lemma 1, the determination processes for the stability of nonlinear discrete systems can be given through Lemma 2, which is described as follows: Lemma 2. For a nonlinear discrete system (), let (i.e., ) be the equilibrium point of the proposed system. Provided that the scalar function concerning satisfies the following: (1) , (2) . then is progressively stable. Proof of Lemma 2. For condition (1), let , and in the case of , , the first condition is proven. For condition (2), let . We can prove that P is a positive definite matrix from Lemma 1. Moreover, , and ; subsequently, the small changes () can be described (), which are calculated as follows: □ Thus, the proof of condition (2) is complete. Based on the above, it is concluded that nonlinear system (4) is asymptotically stable at the origin point. Thus, having proved Lemma 1 and Lemma 2, the proof of Theorem 1 can be obtained, which is as follows: Proof of Theorem 1. According to Equations (1) and (2), Equation (3) can be calculated as follows: where , and is a control function, which be represented as follows: Then, Equation (9) can be simplified as the following equation: Denote the scaled function of the nonlinear error system () represented as , and then the is calculated as follows: Therefore, when the parameter satisfies the condition , ; hence, and . □ Furthermore, . Thus, according to Lemma 2, the nonlinear error system () is asymptotically stable when , and, in turn, the drive and response systems are asymptotically synchronized.

A new 3D discrete chaotic system (12) is proposed in this paper, which is described as follows: where the , , and are iterative variables, and a is a parameter variable. The bifurcation diagram with the a of system (12) is represented by Figure 1. It is clear from the bifurcation diagram that the system is chaotic when a = 6.53, and after 1000 iterations, the Lyapunov exponents of system (12) are 0.7296, 0.1650, and 0.6226, which are all positive; thus, system (12) is a hyperchaotic system. Based on 64-bit Matlab software and double-floating-point representation, the initial values of the state variables of the chaotic system are set differently. The output chaotic sequences and their autocorrelation are shown in Figure 2, from which we can clearly see in Figure 2a,b,d,e that the chaotic sequences and have no periodicity, and Figure 2a,d demonstrate the proposed chaotic sequences with an initial value sensitivity. The initial values of the sequences are set as and , respectively, and the corresponding chaotic iteration diagrams are shown in Figure 2c,f, respectively, which show that system (12) is not in a chaotic state in this case. Therefore, the initial values of the chaotic system are what affect the output states of system (12). In addition, the chaotic attractor phase diagrams of the proposed 3D hyperchaotic mapping are shown in Figure 3 as , , and .

Let us assume that the system of responding systems of the drive system (12) is as follows: where , , and are iteration variables, and based on Theorem 1, the system errors are calculated as follows: Then, the error system equation can be expressed as follows: Let the control function be represented by Equation (16): Thus, system (14) can be simplified to expression (17): Let the Lyapunov exponent function of system (17) be represented by the following expressions: It is obvious that when the parameter is , then , and the error system is progressively stable at , according to Lemma 1. The zero solution of error system (13) is asymptotically stable so that generalized chaotic synchronization can be achieved. In summary, the incorporation of the error-system-feedback coefficient () into the design of the controller for generalized synchronization in this paper makes the design of the controller more flexible. Let the invertible transformation equation be , and matrix A is as follows: In the initial condition for the iterative variables , , and , with 1000 iterations, the dynamical curves of the status variables , , , , , and are displayed in Figure 4a,c,e, respectively, whereas the dynamical curves of , , and are represented in Figure 4b,d,f, respectively, from which it can be seen that the difference in the initial values does not affect the synchronization time.

High-dimensional chaotic systems have more complex dynamics than low-dimensional chaotic systems; thus, they are better able to resist the degradation of dynamics caused by the limited accuracy of computers. In this paper, a new 6D discrete chaotic system (21) is constructed by expanding on system (12), which can be presented as follows: where a = 3.5, b = −2.34, c = 6.53, and d = 3.46, after 1000 iterations, and the Lyapunov exponents of the proposed system are 1.0227, 0.3837, −0.2378, −0.2999, 0.4955, and 1.7245, respectively; thus, system (21) is a hyperchaotic system because the four Lyapunov exponents are positive. Based on 64-bit Matlab software and double-floating-point representation, for different initial values of the state variables, the output of chaotic sequences and their autocorrelations are shown in Figure 5, in which Figure 5a,b,d,e prove that the chaotic sequences and have no periodicity, and Figure 5a,d demonstrate the proposed chaotic sequences with an initial value sensitivity. The initial values of the sequences were set at , and , respectively, and the corresponding chaotic iteration diagrams are shown in Figure 5c,f, respectively, which show that system (21) is not in a chaotic state in this case. Therefore, the initial values of the chaotic system are what affect the output states of system (21). In addition, the chaotic attractor phase diagrams of the proposed 6D hyperchaotic mapping are shown in Figure 6 as , , , , , and .

Let the corresponding system of driving system (21) be considered as follows: Equation (23) is obtained from Lemma 2: The equation of the error system can then be expressed as follows: Let the control functions be calculated as follows: Then, system (24) can be expressed as follows: Consider the Lyapunov exponent functions of system (25) as follows: Obviously, according to Lemma 2, when the parameter and, hence, the error system is progressively stable at , the zero solution of the error system (25) is progressively stable, and as a result, the generalized chaotic synchronization can be implemented. Let the invertible transformation function be , and let the coefficient matrix A be expressed as follows: Moreover, set the initial conditions of the proposed driving system (21) as , , , , , and . After 1000 iterations, the dynamical diagrams of the change in the status of the iterations as k changes for systems (21), (22), and (24) are shown in Figure 7. The Figure 7a–d describe the dynamical diagrams of system (21), and the dynamical diagrams of system (22) are shown in Figure 7e–h, and Figure 7i–l display the dynamical diagrams of system (23).

The framework diagram of the proposed encryption and decryption transmission system constructed in this paper is shown in Figure 8. The proposed system is designed for the encrypted transmission of digital images with pixel matrix values of , , and for each component of the color digital image, where the number of pixels in the image is . Then, matrix components are converted into a sequence of integer values in row order, and the pixel values are selected in the range , converting each pixel value into an 8-bit binary number. Based on the above operations, the binary sequences , , and , which are based on the color image, can be obtained, where . The output sequences () of the proposed 6D generalized discrete chaotic system are quantized by the region to generate binary sequences for encryption. Therefore, the quantification process can be represented by the following equation: where is an arbitrary integer greater than 0, and are denoted as consecutive equal intervals on the interval of a range of real-valued sequences. If the output value of a chaotic sequence is in the odd interval, then it outputs 0, and if it is in the even interval, then it outputs 1. The workflow of the whole system is as follows, and the process of encryption on the transmitter side of the proposed system consists of the following parts: Step 1: The chaotic sequences generated by the proposed 6D discrete chaotic system are quantized and denoted as , , , and . Furthermore, the three color components of the original image are encrypted with chaotic sequences, and the calculation formula is Equation (31):Step 2: The encrypted sequences of the three color components are combined into using Equation (32): Step 3: Chaotic hiding of with the chaotic sequence . The resulting mixed signal () is transmitted in the common channel and is calculated as the following equation: On the receiving side of the proposed system, the response system will be in general synchronization with the driver system. Furthermore, the receiver will be able to decode all the state variables of the sender. Similarly, there are several parts to the decryption processes for the receiver of the proposed system. Step 1: Reconstructing the chaotic signal (), the sequences generated after quantization are denoted as , , , and , as we can see from Figure 8. The decryption process is the inverse of the encryption process; thus, it is important to perform the anti-hiding operation on the signal to obtain , which is calculated as follows: Step 2: Decompose into three color components. The formulation is calculated as Equation (35), and therefore the encrypted image is decoded using Equation (36): The standard Lena (256 × 256) image was used as an example for the encryption and decryption processes, and the results of the operation are shown in Figure 9.

As can be seen from the encryption process in Section 5.1, the entire encryption algorithm consists of simple operations, such as addition, subtraction, and iso-or. However, in the process of encrypting a three-dimensional matrix of the size 8 × m × n, it takes about operations to complete the bit-level operation of the encryption of the xor operation. Hence, the time complexity of the proposed algorithm in this paper is .

Security analyses of the transmission system for the proposed digital image encryption and decryption based on generalized chaotic synchronization are performed in this paper, which consist of image encryption histogram analysis, key space analysis, key sensitivity analysis, and correlation analysis.

A color histogram is a presentation of the statistical characteristics and distribution of the image pixels and is analyzed in terms of three colors: R, G, and B. As can be seen from the simulation results in Figure 10a–c, the histogram of the original image is unevenly distributed, whereas the histogram in Figure 10d–f, of the encrypted digital Lena image in the encryption transmission system with generalized chaotic synchronization, is uniformly distributed.

Cryptographic system security and resistance to exhaustive attacks are affected by the size of the key space. The sequences of the 6D discrete hyperchaotic system (21) are applied at the encryption stage, as the production of keys depends mainly on the parameters and initial conditions of the system. Hence, all the iterative variables of the proposed chaotic system () can be used as keys. Calculation is performed by means of the double-floating-point number with 64-bit precision; therefore, the key space is , which is much large than 2100. Consequently, the security of the image encryption system is improved and the resistance to exhaustive attacks is increased [24,25].

Key sensitivity is an important measure to evaluate the security of an encrypted image. Security is highly reliable if any minimal change to the key results in a large modification; thus, the more sensitive the key, the more secure the encryption system. To evaluate the sensitivity of the key of the proposed scheme, the image was decrypted with a keystream with slight differences, the initial values were set to , and the result of making minor changes to the initial values was . The results of the experiment are illustrated in Figure 11, in which it can be seen that the image could not be decrypted, even with minor changes to the key.

To prevent data from being attacked or lost, it is necessary to perform cutting attack analysis. The same encrypted Lena as in Section 5.1 was cropped with black pixels, which is shown in Figure 12a. Then, we decrypted the image with keys, and the decrypted image can be well identified, as shown in Figure 12b. Thus the encryption system proposed in this paper is resistant to attacks and data loss as well.

The analysis of image pixel correlation is one of the important indicators of the encryption effect of encrypted images. The color images were analyzed for correlations from R, G, and B, which are represented as red, green, and blue colors, respectively. In addition, the original image has a strong correlation between adjacent pixels, and an effective encryption system can reduce the image pixel correlation considerably. The correlation coefficient results for the Lena test images and their corresponding encrypted images are shown in Table 1. As can be seen from Table 1, the correlation coefficient between adjacent pixels in the original image is close to 1, with a strong correlation, whereas in the encrypted image, neighboring pixels are not correlated, as the correlation coefficient is close to 0. To demonstrate the correlation visually, we plotted scatter plots of all sampled pixel pairs. The correlations between the test image Lena and its corresponding R, G, and B color pixels are shown separately. Figure 13a–c,h–j,n–p show that the numerical points of the original images are clustered around the diagonal of the images; hence, there is a strong correlation between neighboring pixels about the original image, whereas the values of the points in Figure 13d–f,k–m and Figure 13q–s are evenly spread throughout the entire plane of the images, which indicates that there is virtually no correlation between neighboring pixels in the encrypted image.

In this paper, two discrete chaotic systems of different dimensions are constructed. Additionally, the dynamics of the new systems are analyzed, and the phase diagram, Lyapunov exponent diagram, and bifurcation diagram of the systems are presented and analyzed simultaneously. The proposed 3D and 6D discrete chaotic systems were constructed as drive systems, and the response systems were constructed by employing the new generalized synchronization method incorporating error-feedback coefficients. The experimental results show that the design of adaptive generalized synchronous systems can be realized provided that the feedback coefficient () of the error system satisfies certain conditions for the design of adaptive generalized synchronous systems. Further, the generalized synchronization method incorporating the error-feedback coefficient, and the incorporation of it into the controller, enables simpler and more flexible control of the generalized synchronization. Finally, a chaotic synchronization and encryption–decryption system for secure digital image transmission was constructed by applying the method of generalized synchronous chaotic systems incorporating the error-feedback coefficients devised in this paper. Due to the limited accuracy of the computer, the system proposed in this paper is more resistant to dynamic degradation and, hence, these features of high-dimensional chaotic systems play an active role and have very good theoretical value in image encryption as well as chaotic synchronization.