
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

39227661
71338
10.1038/s41598-024-71338-x
Article
Dynamic analysis of a novel hyperchaotic system based on STM32 and application in image encryption
Cheng XueFeng 1
Zhu Hongmei 1
Liu Lian 1
Mao Kunpeng 1
Liu Juan liujuan@cque.edu.cn

2345
1 School of Big Data and Information Industry, Chongqing City Management College, Chongqing, 401331 China
2 https://ror.org/02d06s578 grid.495238.1 0000 0000 8543 8239 School of Artificial Intelligence, Chongqing University of Education, Chongqing, 400065 China
3 https://ror.org/023rhb549 grid.190737.b 0000 0001 0154 0904 Chongqing University, Chongqing, 400044 China
4 grid.9227.e 0000000119573309 Chongqing Institute of Green and Intelligent Technology, Chinese Academy of Sciences, Chongqing, 400714 China
5 https://ror.org/05qbk4x57 grid.410726.6 0000 0004 1797 8419 Chongqing College, University of Chinese Academy of Sciences, Chongqing, 400722 China
3 9 2024
3 9 2024
2024
14 2045210 5 2024
27 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
This paper presents a novel 4D hyperchaotic system derived from a modified 3D Lorenz chaotic system. A key aspect of this system is the presence of a single equilibrium point, and its stability is carefully analyzed. The dynamic properties, including the Lyapunov exponent spectrum, bifurcation diagram, and chaotic attractors, are investigated using MATLAB simulations. The results reveal that the system displays hyperchaotic behavior across a wide range of the parameter d, with its chaotic attractor transitioning through four states: hyperchaotic, chaotic, periodic, and quasi-periodic, showcasing the system's complex dynamics. Experimental validation using STM32 embedded hardware successfully reproduces these four types of chaotic attractors, confirming the theoretical predictions. The proposed hyperchaotic system is then applied to image encryption, introducing a novel encryption method. The hyperchaotic key sequence generated by this system meets 15 tests of the NIST SP800-22 standard, and further experimental validation with STM32 hardware demonstrates the algorithm's effectiveness, simplicity, non-linearity, and high security. The encrypted images and sequences are rigorously tested key space analysis, histogram similarity analysis, information entropy analysis, statistical attack analysis, differential attack analysis, key sensitivity analysis, and correlation analysis, highlighting the robustness and reliability of the proposed system. This method is versatile and can be extended to various fields, including audio and video encryption, text encryption, IoT security, financial transaction security, and medical data protection.

Subject terms

Nonlinear phenomena
Electrical and electronic engineering
Science and Technology Project of Chongqing Municipal Education CommissionKJZD-K202303301 KJQN202203308 Cheng XueFeng Zhu Hongmei issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Since Rŏssler discovered the first hyperchaotic system1 in 1979, the evolution of hyperchaotic systems from discovery to application has encompassed multiple stages, including theoretical research2–5, numerical simulation6–9, experimental validation10–14, and application exploration15–18. In recent years, the rapid proliferation of digital images containing sensitive information has heightened the need for robust encryption techniques to protect against unauthorized access and illegal copying19–23. Traditional image encryption methods often struggle with issues such as limited parameter ranges, slow encryption speeds, and vulnerability to various attacks. To address these challenges, researchers have explored the use of hyperchaotic systems, which offer complex dynamic properties and high sensitivity to initial conditions, making them suitable for secure image encryption applications.

Several novel image encryption frameworks and algorithms have been proposed, leveraging different forms of hyperchaotic systems to enhance security and efficiency. A common approach involves multi-stage encryption processes that combine permutation and diffusion techniques. For instance, one framework24 utilizes a new hyperchaotic system based on STM32 and single neuron models (SNM), involving three stages that apply a substitution box (S-box) followed by an XOR operation with encryption keys derived from the numerical solutions of hyperchaotic maps and SNM. This method significantly enhances security with a vast key space and improves encryption efficiency through optimized parallel processing.

Another approach25 employs a 2D hyperchaotic map for designing an efficient pseudo-random number generator (PRNG) and a color image encryption algorithm. The hyperchaotic system generates highly random sequences validated through various analyses, with the PRNG passing all NIST tests. The encryption algorithm uses cross-channel permutation and diffusion to enhance security and speed, though the framework's complexity may require significant computational resources.

Further advancements include the use of a 6D hyperchaotic system combined with random signal insertion. This scheme26 leverages the sum of all plaintext pixels to generate initial values, ensuring high key sensitivity. The encryption process involves splitting each pixel, scrambling, cycle shifting, and diffusion, resulting in robust security against various attacks. Similarly, a fast and secure image encryption algorithm27 utilizes a time-delayed combinatorial hyperchaos map to enhance security and efficiency through simultaneous shuffling and diffusion processes.

Color image encryption algorithms28 have also been developed, incorporating hyperchaotic maps with DNA mutation operations. These methods enhance dynamic properties and introduce additional randomness and complexity, validated through extensive simulations demonstrating robustness against various attacks. A novel image encryption algorithm29 combining a seven-dimensional hyperchaotic map with stochastic signal injection further increases security, employing SHA-512 for initial value generation to ensure strong correlation with plaintext images.

Additionally, The logistic Feigenbaum nonlinear cross-coupled hyperchaotic map30 with dynamic correlation of plaintext pixels offers improved chaotic properties and reduced computational complexity, validated through extensive security analyses. Reservoir computing31 approaches for digital image encryption and hyperchaotic finance model forecasting provide faster and more accurate long-term predictions with reduced execution times, enhancing security and efficiency.

Optical image encryption32 algorithms combining 4D memristive hyperchaotic systems with compressed sensing reduce image size and transmission burden while maintaining high security, validated through extensive simulations. Deep learning-based compression33 combined with two-dimensional Sin-Linear-Cos hyperchaotic maps and matrix encoding ensures high-quality cipher images, fast processing speeds, and robustness against attacks.

Despite the advancements in hyperchaotic-based encryption algorithms, challenges remain in terms of implementation complexity and dependency on numerical precision. The development and application of novel four-dimensional hyperchaotic systems for image encryption demonstrate increased complexity and unpredictability, validated through extensive statistical tests and dynamic analyses16,34,35.

Critically, the cryptanalysis of existing hyperchaotic-based encryption algorithms36 reveals vulnerabilities to chosen plaintext attacks, highlighting the need for continuous improvement in algorithmic security structures. Lastly, the introduction of reservoir computing approaches and non-degenerate hyperchaotic systems for secure DNA-coding image optical communication underscores the potential for enhanced performance and security in optical access networks37.

The reviewed literature illustrates significant progress in developing robust, hyperchaotic-based image encryption algorithms, each offering unique enhancements in security, efficiency, and resistance to various attacks. However, the construction of new hyperchaotic systems and the comprehensive analysis of their dynamic behavior, as well as the engineering implementation of image encryption, remain key considerations for future research and practical applications. Additionally, there are few reported studies on utilizing the high randomness and complexity of hyperchaotic systems combined with embedded hardware STM32 for image encryption.

Based on the above analysis, the research content of this paper includes the following. Firstly, a novel 4D hyperchaotic system based on a modified 3D Lorenz chaotic system is devised and analyzed for its dynamic properties, including equilibrium point stability, Lyapunov exponents’ spectrum, bifurcation diagram, and chaotic attractors. MATLAB numerical simulations confirm its hyperchaotic behavior, and experimental validation using embedded hardware STM32 demonstrates consistency between the hardware results and MATLAB simulations through oscilloscope displays of the attractor phase portrait. Secondly, leveraging the hyperchaotic key sequence generated by this system, color image pixel encryption is achieved via two bit-XOR operations, with experimental findings revealing high randomness in the hyperchaotic key sequence and the encrypted image making it impossible to discern any information about the original image, meeting 15 tests of the NIST SP800-22 standard. Thirdly, analysis against statistical attacks, differential attacks, key sensitivity, correlation, and key space demonstrates that the encrypted data sequence resembles random noise, underscoring the high complexity and nonlinearity of the hyperchaotic key sequence. Lastly, the conclusions are drawn.

Dynamical analysis of a new hyperchaotic system

An 4D hyperchaotic system38 is described by (1):1 x˙=40y-x+wy˙= 10x+25y-xzz˙=xy- 3zw˙=dx

where x, y, z and w are driving variables, and d is a positive real parameter. Varying the parameter d can cause system (1) to exhibit hyperchaotic, chaotic, periodic, or quasi-periodic motion. Notably, the hyperchaotic system (1) has only one equilibrium point O0,0,0,0, which is unstable and acts as a saddle point.

Linearizing the system (1) and calculating the Jacobian matrix at the equilibrium point O0,0,0,0 is shown as follows:2 J0=-404001100d25000-30000

Let detJ0-λI=0, the eigenvalues of system (1) in the equilibrium point O0,0,0,0 is:3 λ1=-40,λ2=25,λ3=-3,λ4=0.

Herein, λ2=25 is positive, λ1=-40 and λ3=-3 are negative. So, the equilibrium point O0,0,0,0 is saddle point of system (1), and it is unstable. It indicates the new system (1) is dissipative, and can be obtained:4 ∇V=∂x˙∂x+∂y˙∂y+∂z˙∂z+∂w˙∂w=-40+25-3=-18.

That is, each volume containing the system orbit shrinks to zero as t→∞ at an exponential rate -18. Consequently, all orbits of system (1) will ultimately be confined to a specific subset of zero volume.

Numerical simulation results

MATLAB excels in numerical simulation with advanced mathematical libraries, matrix operation advantages, a flexible programming language, and support for parallel computing. Its extensive toolbox spans multiple domains, and powerful plotting and visualization tools aid in intuitively showcasing simulation results. Therefore, this paper employs MATLAB simulation software for the numerical analysis of the hyperchaotic system.

According to chaos theory, the adjacent trajectories of a chaotic attractor in the state space separate from each other exponentially. The spectrum of Lyapunov exponents quantitatively shows the contraction or expansion of these trajectories.

Figure 1 shows the Lyapunov exponents’ spectrum of system (1) with the changing of parameter d. It can be observed that as the parameter d changes, the four Lyapunov exponents of the constructed new system (1) also vary. The largest Lyapunov exponent remains positive over a significant range of parameter d, and the second-largest Lyapunov exponent exhibits positive as parameter d changes. Therefore, the system (1) demonstrates hyperchaotic characteristics.Fig. 1 Lyapunov exponents’ spectrum of system (1).

Figure 2 shows the bifurcation diagram of state variable x of system (1). It is evident that as the parameter d gradually changes, the bifurcation diagram reflects the instability and complexity of the system's solutions x, providing information about the transition from order to chaos. When the bifurcation diagram exhibits more and increasingly complex branches, it indicates that the dynamical behavior of the new system (1) becomes more intricate, potentially entering a chaotic state.Fig. 2 Bifurcation diagram of x.

It is found that both of the Lyapunov exponents’ spectrum and the bifurcation diagram match, indicating that the numerical simulation results are correct. The variation of system (1) with parameter d is depicted in Table 1. For 1<d≤23.5, system (1) exhibits hyperchaotic behavior.Table 1 Stability of system (1) under different parameter d.

Parameter d	Lyapunov exponents	Stability	
0<d≤1,

23.5<d≤78.6

	LE1>0,LE2=0,LE3<0,LE4<0	Chaotic	
1<d≤23.5	LE1>0,LE2>0,LE3=0,LE4<0	Hyperchaotic	
78.6<d≤113,

177.5<d≤200

	LE1=0,LE2<0,LE3<0,LE4<0	Periodic	
113<d≤177.5	LE1=0,LE2=0,LE3<0,LE4<0	Quasi-periodic	

Let d=15, the Lyapunov exponents of system (1) is 1.8451,0.1832,0,-19.9925, the hyperchaotic attractor is shown in Fig. 3a. Let d=70, the Lyapunov exponents of system (1) is 1.2823,0,-0.5963,-18.7196, the chaotic attractor is shown in Fig. 3b. Let d=100, the Lyapunov exponents of system (1) is 0,-1.0039,-0.9961,-16.0151, the periodic orbit attractor is shown in Fig. 3c. Let d=120, the Lyapunov exponents of system (1) is 0,0,-2.1678,-15.8563, the quasi-periodic orbit attractor is shown in Fig. 3d.Fig. 3 The attractors of system (1) ((a) hyperchaotic attractor (d=15); (b) chaotic attractor (d=70); (c) periodic attractor (d=100); (d) quasi-periodic attractor (d=120)).

Embedded hardware implementation results

In embedded hardware, common components include DSP and FPGA, with DSP widely applied in various embedded implementations due to its flexibility. The STM32 microcontroller, representing a typical embedded hardware, integrates numerous peripherals, prioritizes power optimization, and offers flexible configuration options. Extensively utilized in industrial control, automotive electronics, medical devices, and consumer electronics, it meets the requirements of diverse application scenarios.

In this experiment, an STM32F103ZET6 microcontroller is used to simulate and control the phase diagram of a hyperchaotic system. The results of the embedded system are displayed using a Tektronix TBS1072C digital oscilloscope. The implementation of hyperchaotic systems on STM32 has significant purposes and implications:

Hyperchaotic systems are complex nonlinear dynamical systems that exhibit chaotic behavior. Implementing them on STM32 can provide valuable insights into the characteristics and behavior of hyperchaotic systems, as well as challenges related to modeling, simulation, and control of nonlinear dynamical systems.

Hyperchaotic systems have wide-ranging applications in communication, encryption, control, and other fields. Implementing them on STM32 can provide a solid foundation for basic research and experimental validation in these fields.5 h1_x=40×yτ-xτ+wτ×hh1_y=10×xτ+25×yτ-xτ×zτ×hh1_z=xτ×yτ-3×zτ×hh1_w=d×xτ×htemp_x=xτ+h1_w/2temp_y=yτ+h1_y/2temp_z=zτ+h1_z/2temp_w=wτ+h1_w/2h2_x=40×temp_y-temp_x+temp_w×hh2_y=10×temp_x+25×temp_y-temp_x×temp_z×hh2_z=temp_x×temp_y-3×temp_z×hh2_w=d×temp_x×htemp_x=xτ+h2_x/2temp_y=yτ+h2_y/2temp_z=zτ+h2_z/2temp_w=wτ+h2_w/2h3_x=40×temp_y-temp_x+temp_w×hh3_y=10×temp_x+25×temp_y-temp_x×temp_z×hh3_z=temp_x×temp_y-3×temp_z×hh3_w=d×temp_x×htemp_x=xτ+h3_xtemp_y=yτ+h3_ytemp_z=zτ+h3_ztemp_w=wτ+h3_wh4_x=40×temp_y-temp_x+temp_w×hh4_y=10×temp_x+25×temp_y-temp_x×temp_z×hh4_z=temp_x×temp_y-3×temp_z×hh4_w=d×temp_x×hxτ+1=xτ+h1_x+2×h2_x+2×h3_x+h4_x/6yτ+1=yτ+h1_y+2×h2_y+2×h3_y+h4_y/6zτ+1=zτ+h1_z+2×h2_z+2×h3_z+h4_z/6wτ+1=wτ+h1_w+2×h2_w+2×h3_w+h4_w/6

Implementing hyperchaotic systems on STM32 requires overcoming technical challenges such as limited hardware resources and complex algorithm implementation, which can promote the development of embedded systems and algorithm optimization.

In experiment, the STM32 is used to verify the 4th-order Runge–Kutta method used for the discretization of the new hyperchaotic system (1), the iterative process is shown in Eq. (5). Where parameter d=15, iteration step h=0.02. τ represents the current moment, xτ, yτ, zτ and wτ are the values of the current moment, temp_x, temp_y, temp_z and temp_w are the intermediate variables, h1_x, h1_y, h1_z and h1_w are the slopes at the beginning of the time period, h2_x, h2_y, h2_z and h2_w are the slopes of the midpoint of the time period, the slopes h1_x, h1_y, h1_z and h1_w are used by Euler's method to determine the values of xτ, yτ, zτ and wτ at the point τ+h/2, h3_x, h3_y, h3_z and h4_w are the slopes of the midpoint, but this time the slopes h2_x, h2_y, h2_z and h2_w are used to determine the xτ, yτ, zτ and wτ values, h4_x, h4_y, h4_z and h4_w are the slopes of the end of the time period, and their xτ, yτ, zτ and wτ values are determined by h3_x, h3_y, h3_z and h3_w. Finally, the values xτ+1, yτ+1, zτ+1 and wτ+1 of the τ+1 moment are calculated.

Then, the digital signals are converted into analog signals through SPI, which is controlled by the DAC converter, and then displayed on the oscilloscope. The embedded system STM32 is depicted in Fig. 4.Fig. 4 Experimental platform for STM32 implementation.

When d=15, Fig. 5a shows the phase transition trajectory of hyperchaotic signal by oscilloscope. When d=70, Fig. 5b shows the phase transition trajectory of chaotic signal by oscilloscope. When d=100, Fig. 5c shows the phase transition trajectory of periodic orbit by oscilloscope. When d=120, Fig. 5d shows the phase transition trajectory of quasi-periodic orbit by oscilloscope.Fig. 5 Experimental results of STM32 implementation ((a) hyperchaotic orbit (d=15); (b) chaotic orbit (d=70); (c) periodic orbit (d=100); (d) quasi-periodic orbit (d=120)).

The alignment of simulation outcomes from MATLAB with the results derived from the hyperchaotic system implemented on the STM32, as illustrated in Figs. 3 and 5, establishes the coherence between the two datasets. This alignment not only affirms the precision and dependability of the simulation model but also underscores the efficacy of deploying the hyperchaotic system on the STM32 platform.

Image encryption

In today's digital society, image encryption, as a crucial branch of information security, plays a key role in safeguarding the privacy of image data and preventing unauthorized access. With the widespread use of images in communication, storage, and sharing, effectively encrypting them has become paramount. Image encryption aims to transform images into unintelligible forms to ensure the security of personal privacy, business secrets, and sensitive information through the utilization of complex mathematical algorithms and keys. The advancement of encryption technology provides essential tools for preventing image information leakage, safeguarding national security, and ensuring confidentiality in fields like medical imaging.

The image encryption technology based on hyperchaotic system-generated hyperchaotic keys exhibits significant necessity and superiority. The hyperchaotic system, with its high dimensionality, complexity, and nonlinear characteristics, forms the basis for generating keys with extremely high randomness, thereby enhancing the security of encryption algorithms. This technology effectively addresses modern challenges in image security by generating complex and unpredictable key sequences, protecting images from unauthorized access and information leakage. Compared to traditional encryption methods, image encryption based on hyperchaotic systems is more resilient to attacks. Supported by embedded hardware like STM32, it achieves efficient encryption and decryption processes, providing an advanced and practical solution for image encryption in the field of information security.

Encryption method

To highlight the simplicity and non-linear advantages of implementing hyperchaotic keys, as well as fully leveraging the benefits of ciphertext interleaving diffusion technology in image encryption to improve its ability to resist illegal attacks, a method of image encryption bit-XOR operation based on hyperchaotic key sequence is proposed. Characterized by suitability for image encryption, non-linear ciphertext, and easy implementation, this approach enhances the speed of ciphertext diffusion. The image encryption process based on the generation of a hyperchaotic key sequence from the hyperchaotic system is shown in Fig. 6, involves four steps:Fig. 6 The image encryption process using the hyperchaotic key sequence.

Firstly, construct a hyperchaotic system and analyze its dynamic characteristics to demonstrate its complexity, randomness, and nonlinearity, as described in section “Dynamical analysis of a new hyperchaotic system”.

Secondly, generate hyperchaotic key sequences that meet the requirements based on the constructed hyperchaotic system, with each hyperchaotic key represented by 16 bits, as detailed in subsection “Hyperchaotic key sequence”. Next, decompose the color image to be encrypted (with each pixel represented by 8 bits) into three matrices representing the RGB channels. Convert each matrix into a column vector and concatenate them to form a complete column vector, serving as the data sequence for the image, as shown in Fig. 7.Fig. 7 The process diagram of converting the original image into a data sequence.

Then, perform bitwise bit-XOR operations between each value of the image sequence and the corresponding value of the hyperchaotic sequence, obtaining an 8-bit encrypted sequence according to the following formula (6).6 R=S⊕KH⊕KL

where R sequence is the encrypted image sequence using bit-XOR operation, S sequence is the original image sequence, and K sequence is the hyperchaotic key sequence, KH and KL represent the high 8-bit and low 8-bit sub-sequences of the key sequence K.

Finally, output the encrypted image sequence through the serial port of the embedded hardware STM32. The decryption of the image follows the reverse process of encryption, with the decryption bit-XOR operation being consistent with the encryption process.

Hyperchaotic key sequence

Due to the inherent high randomness and complexity of the hyperchaotic system, the encryption algorithm's security is significantly strengthened. The generation of the hyperchaotic key sequence involves transforming the hyperchaotic sequence into the required key sequence through a specific quantization algorithm, outlined in the following steps:

Step 1: Pre-iterate the hyperchaotic system N1 times to eliminate transient effects as it enters a hyperchaotic state.

Step 2: Iterate the hyperchaotic system N2 times to obtain a new set of state values, A=Ax,Ay,Az,Aw, where Ax={x1,x2,x3,x4}T, Ay={y1,y2,y3,y4}T, Az={z1,z2,z3,z4}T, and Aw={w1,w2,w3,w4}T, 0<k≤N2. Scale down the state values A to 1/p and take q decimal places to generate new state values B=Bx,By,Bz,Bw, where Bx={bx1,bx2,bx3,bx4}T, By={by1,by2,by3,by4}T, Bz={bz1,bz2,bz3,bz4}T, and Bw={bw1,bw2,bw3,bw4}T, 0<k≤N2. The conversion formula is given by:7 bxi=xi∗10p-q-xi∗10p-q,ifxi≥0,i=1,2,⋯,kbxi=xi∗10p-q-xi∗10p-q,ifxi<0,i=1,2,⋯,kbyi=yi∗10p-q-yi∗10p-q,ifyi≥0,i=1,2,⋯,kbyi=yi∗10p-q-yi∗10p-q,ifyi<0,i=1,2,⋯,kbzi=zi∗10p-q-zi∗10p-q,ifzi≥0,i=1,2,⋯,kbzi=zi∗10p-q-zi∗10p-q,ifzi<0,i=1,2,⋯,kbwi=wi∗10p-q-wi∗10p-q,ifwi≥0,i=1,2,⋯,kbwi=wi∗10p-q-wi∗10p-q,ifwi<0,i=1,2,⋯,k

where · denotes the floor function, p and q are adjustable parameters, representing scaling down the value by 1/p and taking q decimal places as the new state value.

Step 3: Calculate the key sequence. Adjust the order of the state values B based on positive and negative changes in the state values A to generate the key sequence =K1,K2,K3,K4, 0<k≤N2, following the mapping relationship outlined in Table 2.Table 2 The mapping relationship illustrating the sorting rule of state values A and B in correlation with the generation of the key sequence K.

State values Ak	State values Bk	Key sequence K	
xk	yk	zk	wk	
+	+	+	+	xk,yk,zk,wk	bx1,by1,bz1,bw1,bx2,by2,bz2,bw2,bx3,by3,bz3,bw3,⋯⋯bxk,byk,bzk,bwk,

0<k≤N2

	
+	+	+	–	xk,yk,wk,zk	
+	+	–	+	xk,zk,yk,wk	
+	+	–	–	xk,zk,wk,yk	
+	–	+	+	yk,xk,zk,wk	
+	–	+	–	yk,xk,wk,zk	
+	–	–	+	yk,zk,xk,wk	
+	–	–	–	yk,zk,wk,xk	
–	+	+	+	zk,xk,yk,wk	
–	+	+	–	zk,xk,wk,yk	
–	+	–	+	zk,yk,xk,wk	
–	+	–	–	zk,yk,wk,xk	
–	–	+	+	wk,xk,yk,zk	
–	–	+	–	wk,xk,zk,yk	
–	–	–	+	wk,yk,xk,zk	
–	–	–	–	wk,yk,zk,xk	

The randomness of the hyperchaotic key sequence can be evaluated using the NIST SP800-22 standard. This standard, developed by the National Institute of Standards and Technology (NIST) in the United States, is one of the authoritative standards for pseudo randomness testing. It includes a total of 15 test metrics, examining the deviation of the tested sequence from ideal random sequences from different perspectives in terms of statistical characteristics. With a significance level of 0.01, a test group number of 100, a group length of 1024000 bits, and a confidence interval of 0.93,1, the test results are presented in Table 3. For items with a test frequency not less than 2, Table 3 provides P-values and the minimum values of pass rates. It can be observed that the hyperchaotic key sequence generated by the algorithm has successfully passed all 15-test metrics, further confirming its excellent randomness.Table 3 The NIST SP 800-22 test result for hyperchaotic key sequence.

No	Project name	Number of tests	P-values	Pass rates	Test results	
1	Frequency	1	0.000017	98%	Yes	
2	BlockFrequency	1	0.000014	98%	Yes	
3	CumulativeSums	2	0.000058	96%	Yes	
4	Runs	1	0.000033	97%	Yes	
5	LongestRun	1	0.514124	99%	Yes	
6	Rank	1	0.000051	99%	Yes	
7	FFT	1	0.071177	95%	Yes	
8	NonOverlappingTemplate	148	0.030806	98%	Yes	
9	OverlappingTemplate	1	0.021999	100%	Yes	
10	Universal	1	0.595549	100%	Yes	
11	ApproximateEntropy	1	0.389596	95%	Yes	
12	RandomExcursions	8	0.468595	97%	Yes	
13	RandomExcursionsVariant	18	0.350485	98%	Yes	
14	Serial	2	0.474986	97%	Yes	
15	LinearComplexity	1	0.319084	93%	Yes	

Experimental results

Since the STM32 achieves synchronized control of the encryption module and the interface module through the master clock, eliminating the need for separate control of the hyperchaotic system's encryption and decryption processes. This approach avoids the synchronization issues between the encryption and decryption modules that are commonly encountered in traditional encryption systems.

Specifically, in this experiment, image encryption based on the hyperchaotic key sequence was implemented using the embedded hardware STM32F103ZET6. Due to the absence of an integrated image acquisition and display module in the STM32F103ZET6, data input and output were facilitated in the form of a list. The resulting encrypted and decrypted image sequences were then transmitted to the PC via the serial port of the embedded hardware STM32F103ZET6. Subsequently, the data was visualized and printed using a serial port debugging tool, as shown in Fig. 8.Fig. 8 The image encrypted and decrypted sequences printed by the serial port debugging tool.

Although the embedded hardware STM32F103ZET6 can complete the experimental process, it lacks data analysis tools, making it difficult to visually compare the information before and after image encryption. Therefore, it is necessary to transmit the data sequence generated in the embedded hardware STM32F103ZET6 to the PC through the serial port and use simulation software MATLAB for analysis. Figure 9a and d show the color image before encryption, Fig. 9b and e display the color image after encryption, Fig. 9c and f depict the color image after decryption.Fig. 9 Image and histogram of encryption/decryption.

Comparing Fig. 9a, b, c, and d, e, f, the original images are the 256*256 color picture. After encryption with the hyperchaotic key sequences, the images are completely flooded with noise, making them impossible to discern the real information of the original image. After decryption, the images show no differences in the time domain compared to the original images. The serial port data printed by the embedded hardware STM32F103ZET6 also confirms that the image sequences before and after encryption are identical.

The real information of the original image is completely concealed after encryption, and the encrypted image sequence exhibits a Gaussian random white noise distribution, validating the high security of the image encryption method based on the hyperchaotic key sequence. The decrypted image is identical to the original image, demonstrating that image encryption based on the hyperchaotic key sequence does not have a negative impact on the original data, and using the STM32F103ZET6 embedded hardware system as the encryption device does not introduce noise. Therefore, the image encryption method based on the hyperchaotic key sequence generated by the hyperchaotic system is both secure and effective.

Performance evaluation

To evaluate the performance of this scheme, five key aspects were analyzed: key space analysis, histogram similarity analysis, information entropy analysis, statistical attack analysis, differential attack analysis, key sensitivity analysis, and correlation analysis. The performance analysis was conducted on a Windows operating system using Python software, with a PC configuration of an Intel i5-12450H 2.0 GHz CPU, 16 GB RAM, and 512 GB ROM. MATLAB 2019a was used for the analysis software.

Key space analysis

The key space size is a crucial factor in determining whether an encryption algorithm can resist brute-force attacks. Generally, an encryption algorithm is considered secure against such attacks when its key space exceeds 2100. For the hyperchaotic mapping, the state values x,y,z,w range between -100,100, and six decimal places are used as the key. With double-precision floating-point data, accurate to 15 decimal places, the key space can reach approximately 2×1015×2×1015×2×1015×2×1015≈2203, which is equivalent to a 203-bit key length. When including system parameters d and pre-iteration times N0, the key space expands even further, providing strong resistance against exhaustive attacks.

Table 4 presents the key space test results of different schemes, showing that the reference schemes have larger key spaces than the proposed method. However, those systems are more complex, with dimensions exceeding six, which affects encryption efficiency and storage requirements. Although the key space in this scheme is relatively smaller, it still far exceeds the minimum key space requirements for encryption systems.Table 4 Key space of different schemes.

Encryption scheme	Proposed	Ref. 28	Ref. 39	Ref. 40	Ref. 41	Ref. 42	
Key space	>2203	2398	2395	2246	2325	2356	

Histogram similarity analysis

To comprehensively demonstrate the effectiveness of the proposed hyperchaotic system for encryption, 71 images were selected from the MATLAB 2019a toolbox for encryption validation. The effectiveness of the encryption scheme was verified by analyzing the histogram similarity before and after encryption. The test results are shown in Table 5.Table 5 Histogram similarity before and after image encryption.

No	Test images	Image sizes	Histogram similarity	No	Test images	Image sizes	Histogram similarity	
1	baby	2240*3584*3	0.04420524	37	llama	1280*832*3	0.15345570	
2	bag	128*192*3	0.04760596	38	lowlight_1	704*960*3	0.18525304	
3	blobs	320*256*3	0.11474805	39	lowlight_2	1984*1280*3	0.13422752	
4	car_1	640*448*3	0.04177567	40	micromarket	3456*2304*3	0.04981033	
5	car_2	640*448*3	0.18803619	41	office_1	896*576*3	0.10119631	
6	car_3	640*448*3	0.1295514	42	office_2	896*576*3	0.16850714	
7	car_4	640*448*3	0.05781802	43	office_3	896*576*3	0.2954075	
8	car1	3456*2304*3	0.05509924	44	office_4	896*576*3	0.05546688	
9	car2	3456*2304*3	0.05189874	45	office_5	896*576*3	0.14683144	
10	circles	256*256*3	0.0089738	46	office_6	896*576*3	0.13143341	
11	circlesBrightDark	512*512*3	0.09751573	47	onion	256*256*3	0.05183778	
12	coins	256*192*3	0.07121657	48	parkavenue	2048*1536*3	0.14046591	
13	coloredChips	512*384*3	0.2908124	49	peacock	1024*768*3	0.13711488	
14	concordaerial	3008*1984*3	0.14349766	50	pears	704*448*3	0.14576219	
15	concordorthophoto	2944*2176*3	0.13939115	51	pillsetc	512*384*3	0.18243526	
16	DistortedImage	448*640*3	0.13410182	52	printedtext	1600*896*3	0.14413334	
17	eSFRTestImage	3072*1792*3	0.23765657	53	saturn	1152*1472*3	0.29165825	
18	fabric	640*448*3	0.14163200	54	sevilla	1600*896*3	0.05219821	
19	flamingos	1280*960*3	0.04767290	55	sherlock	960*640*3	0.14022593	
20	flower	256*256*3	0.27108517	56	snowflakes	320*64*3	0.13637836	
21	foggyroad	2048*1536*3	0.27971324	57	strawberries	1024*640*3	0.14117495	
22	foggysf1	3456*2304*3	0.05197918	58	tape	512*384*3	0.13746690	
23	foggysf2	3456*2304*3	0.04182551	59	testpat1	256*256*3	0.15047596	
24	foosball	3456*2304*3	0.14954731	60	text	256*256*3	0.00921987	
25	football	320*256*3	0.13200099	61	threads	448*448*3	0.2376722	
26	gantrycrane	384*256*3	0.23377732	62	tissue	768*448*3	0.14848739	
27	glass	256*128*3	0.14447569	63	toyobjects	320*320*3	0.08204885	
28	greens	448*256*3	0.13989351	64	toysflash	896*640*3	0.23470400	
29	hallway	3456*2304*3	0.2446188	65	toysnoflash	896*640*3	0.13711072	
30	hands1	192*320*3	0.2966372	66	trailer	1024*640*3	0.14140032	
31	hands1-mask	320*192*3	0.08397624	67	wagon	768*1024*3	0.2033131	
32	hands2	192*320*3	0.2891822	68	westconcordaerial	320*384*3	0.13978115	
33	hestain	256*192*3	0.23509654	69	westconcordorthophoto	320*320*3	0.2870106	
34	indiancorn	1024*768*3	0.23153046	70	yellowlily	1216*1600*3	0.27546623	
35	kobi	1600*1216*3	0.05287075	71	yellowlily-segmented	1216*1600*3	0.07700115	
36	lighthouse	448*640*3	0.13331444					

Histogram similarity is an effective metric for evaluating the performance of an image encryption scheme. A lower histogram similarity indicates a significant difference between the histograms of the original and encrypted images, suggesting that the encryption algorithm has successfully altered the pixel distribution, making it difficult to find any correlation between the two. In practical applications, the histogram similarity after image encryption should generally be less than 0.3. As shown in Table 5, among the 71 images tested, the highest histogram similarity was 0.2966372, and the lowest was 0.0089738, which meets the requirements for evaluating the effectiveness of image encryption based on histogram similarity.

Information entropy analysis

Information entropy is a crucial metric for measuring the amount of information or uncertainty within an image, reflecting the complexity of the image by evaluating the randomness of pixel values. In image encryption, high entropy indicates that the encrypted image possesses a high degree of randomness and uncertainty, making it extremely difficult to reverse-engineer the original image through statistical methods. Therefore, entropy is commonly used to assess the effectiveness of encryption algorithms. Ideally, the entropy of an encrypted image should approach the maximum value (typically close to 8 for an 8-bit image), indicating a uniform pixel distribution and that the image resembles random noise, effectively safeguarding the original image information.

Table 6 presents the entropy values before and after encryption for 71 images, showing that the entropy of the encrypted images is close to 8, demonstrating the high randomness and uncertainty achieved by the proposed encryption scheme.Table 6 Information entropy before and after image encryption.

No	Test images	Information entropy	No	Test images	Information entropy	
Original image	Encrypted image	Original image	Encrypted image	
1	baby	7.694291	7.996322	37	llama	7.364563	7.999322	
2	bag	7.489169	7.998253	38	lowlight_1	6.410127	7.999316	
3	blobs	1.303949	7.998321	39	lowlight_2	6.247875	7.996329	
4	car_1	7.658223	7.996329	40	micromarket	7.816366	7.997327	
5	car_2	6.473951	7.997325	41	office_1	4.870706	7.997320	
6	car_3	5.555290	7.996325	42	office_2	6.319726	7.998307	
7	car_4	1.109169	7.997319	43	office_3	7.205910	7.996309	
8	car1	7.524531	7.996324	44	office_4	7.561897	7.998328	
9	car2	7.742344	7.998326	45	office_5	7.667470	7.999335	
10	circles	0.797671	7.996289	46	office_6	6.883114	7.999316	
11	circlesBrightDark	1.288123	7.999320	47	onion	7.436057	7.997305	
12	coins	6.346112	7.998320	48	parkavenue	7.157566	7.998332	
13	coloredChips	7.204383	7.996335	49	peacock	7.447847	7.996323	
14	concordaerial	7.175397	7.997325	50	pears	7.264007	7.998329	
15	concordorthophoto	7.355761	7.999322	51	pillsetc	6.007417	7.998332	
16	DistortedImage	7.854813	7.996333	52	printedtext	7.679555	7.999326	
17	eSFRTestImage	6.739018	7.996324	53	saturn	4.985464	7.999323	
18	fabric	7.348980	7.998306	54	sevilla	7.839214	7.996324	
19	flamingos	7.359811	7.996319	55	sherlock	7.366483	7.999317	
20	flower	7.508741	7.996326	56	snowflakes	4.872856	7.998238	
21	foggyroad	6.612072	7.999319	57	strawberries	7.642410	7.998328	
22	foggysf1	7.450419	7.999330	58	tape	6.741722	7.998323	
23	foggysf2	7.180980	7.998328	59	testpat1	2.611444	7.999310	
24	foosball	7.547328	7.996327	60	text	0.454175	7.997289	
25	football	6.713418	7.999297	61	threads	5.186503	7.996315	
26	gantrycrane	6.683417	7.996299	62	tissue	7.534090	7.998314	
27	glass	7.340917	7.996271	63	toyobjects	1.597250	7.999326	
28	greens	7.422649	7.998303	64	toysflash	7.274874	7.998305	
29	hallway	7.286544	7.996329	65	toysnoflash	7.404918	7.996324	
30	hands1	6.848025	7.997260	66	trailer	7.362810	7.999330	
31	hands1-mask	1.052358	7.997301	67	wagon	7.707832	7.996315	
32	hands2	6.891528	7.999282	68	westconcordaerial	7.166493	7.999324	
33	hestain	7.277335	7.998270	69	westconcordorthophoto	7.739572	7.997343	
34	indiancorn	7.737863	7.999322	70	yellowlily	7.039695	7.999315	
35	kobi	7.420971	7.997329	71	yellowlily-segmented	1.599543	7.998321	
36	lighthouse	7.830373	7.997321					

Table 7 compares the entropy values of images after encryption using the proposed scheme with other schemes, further illustrating that the proposed method effectively confuses pixel values and enhances the security of the encrypted images.Table 7 Information entropy of different schemes.

No	Encryption scheme	Test images	Information entropy
(Encrypted image)	
1	Proposed	Average	7.997826	
2	Ref. 28	Average	7.998300	
3	Ref. 43	Baboon	7.977383	
4	Ref. 44	Peppers	7.996367	
5	Ref. 45	Baboon	7.998233	

Statistical attack analysis

The distribution histograms of the original image and the encrypted image are shown in Fig. 10a and c, b and d respectively.Fig. 10 Histograms distribution of images before and after encryption.

From Fig. 10, it can be observed that the distribution of the original image histogram is uneven, with a higher probability density in the low pixel segment and a lower probability density in the high pixel segment. In addition, the distribution of the encrypted image is flat and uniform, with similar probability densities in the low and high pixel segments. This indicates that the probability of pixel values in the encrypted image tends to be equal. Therefore, this encryption algorithm can effectively resist statistical analysis attacks.

Differential attack analysis

According to the principles of cryptography, the stronger the sensitivity of the key to plaintext, the stronger the resistance against differential attacks. The sensitivity of plaintext can be measured using the number of signals change rate (NSCR) and the unified average changing intensity (UACI). They respectively represent the proportion and degree of change in pixel values of the encrypted image when a pixel value of the original image is randomly altered. Assuming that the values of two ciphertext signals at point i are Ci and C′i, if Ci=C′i, then Di=0; if Ci≠C′i, then Di=1. Then the definitions of NSCR and UACI are as follows Eq. (8):8 NSCR=∑i=1NDiN×100\%UACI=1N∑i=1NC′i-Ci255×100\%

where N is the number of image pixels. In this experiment, 256 groups were selected for encryption, with each group consisting of 2 image sequences: one original image sequence and one encrypted image sequence.

The NSCR and UACI value curves for the 256 sets of ciphertext signals obtained from the Flower image are shown in Fig. 11a and b, respectively, with average values of 99.59% for NSCR and 35.22% for UACI. Similarly, the NSCR and UACI value curves for the 256 sets of ciphertext signals obtained from the Onion image are displayed in Fig. 11c and d, with average values of 99.60% for NSCR and 33.84% for UACI. These results indicate that the bit-XOR operation causes slight changes in the lower bits of the plaintext to induce significant changes in the higher bits of the ciphertext, while changes in the higher bits of the plaintext can simultaneously affect both the higher and lower bits of the ciphertext. This enhances the overall variation magnitude of the ciphertext. Additionally, the differential attack test results from other schemes28,43–47 are summarized in Table 8, demonstrating that the proposed algorithm not only resists differential attacks but also performs stably.Fig. 11 Plaintext sensitivity test.

Table 8 NPCR and UACI of different schemes.

No	Encryption scheme	Test images	NPCR	UACI	
1	Proposed	Flower	99.59%	35.22%	
Onion	99.60%	33.84%	
2	Ref. 28	Average	99.61%	33.47%	
3	Ref. 43	Lena	99.60%	33.88%	
4	Ref. 44	Average	99.62%	33.40%	
5	Ref. 45	Baboon	99.63%	33.51%	
6	Ref. 46	Lena	99.61%	33.44%	
7	Ref. 47	Lena	99.57%	33.44%	

Key sensitivity analysis

A good cryptographic algorithm must be highly sensitive to the key, meaning that even slight differences in encryption keys should result in significant changes in the ciphertext sequence for the same plaintext. Similarly, for the same ciphertext, slight differences in the decryption key should lead to vastly different decryption results.

In this experiment, the initial values x0,y0,z0,w0 of the system are used as the key. Initially, the encryption is performed using the key 1,1,1,1, followed by decryption using five sets of decryption keys with minor differences, including 1+10-15,1,1,1, 1,1+10-15,1,1, 1,1,1+10-15,1, 1,1,1,1+10-15, and 1+10-15,1+10-15,1+10-15,1+10-15.

Figure 12a is the original image, Fig. 12b is the decrypted image when the error in the initial decryption key x0 is 10-15, Fig. 12c is the decrypted image when the error in the initial decryption key y0 is 10-15, Fig. 12d is the decrypted image when the error in the initial decryption key z0 is 10-15, Fig. 12e is the decrypted image when the error in the initial decryption key w0 is 10-15, Fig. 12f is the decrypted image when the errors in the initial decryption keys x0, y0, z0, and w0 are all 10-15. Where it can be observed that the correctly decrypted signal matches the original image, indicating successful data recovery, while the remaining five sets of erroneous decryption signals resemble noise signals, demonstrating the extreme sensitivity of the algorithm to the key.9 MSE=1N∑i=1NP′i-Pi2

Fig. 12 Correctly decrypted image and incorrectly decrypted images.

Let the original image be denoted as P and the decrypted image as P′. Their mean square error is calculated using Eq. (9), and the mean square error between the erroneous decryption signal and the original image is presented in Fig. 13. It can be observed that even slight errors in the decryption key result in vastly different decryption results, thus numerically proving the sensitivity of the algorithm to the key.Fig. 13 The mean square error of error decryption to original image.

Correlation analysis

Pixel values in images are not independent, with small differences in adjacent pixel amplitudes leading to high correlation. Therefore, one of the objectives of encryption algorithms is to reduce the correlation between adjacent pixel values. Lower correlation implies better confusion effects and higher security. Correlation can be measured using correlation coefficients, calculated between adjacent amplitude values of plaintext and ciphertext, as fallowing Eq. (10).10 Corr=∑i=1Mxi-x¯yi-y¯∑i=1Mxi-x¯2∑i=1Myi-y¯2

Let xi and yi represent the i-th pair of adjacent pixel values randomly selected from the image pixel value matrix, and let x¯ and y¯ be their respective averages. M represents the total number of pairs, with M set to 8, 64, 128, and 256. The calculated results are presented in Table 9.Table 9 Correlation of adjacent flower image pixels.

Number of plaintext-ciphertext pairs M	Original image correlation	Encrypted image correlation	
Red	Green	Blue	Average	Red	Green	Blue	Average	
8	0.9964	0.9985	0.9983	0.9977	0.0240	0.4134	− 0.3087	0.0429	
64	0.9905	0.9796	0.9909	0.9870	− 0.0819	0.0687	0.1274	0.0381	
128	0.9889	0.9881	0.9947	0.9905	0.0626	0.0299	− 0.1489	− 0.0188	
256	0.9896	0.9797	0.9865	0.9853	0.0350	− 0.0140	− 0.0784	− 0.0191	

A comparison of the data reveals that adjacent original images have a high correlation Corr→1, while adjacent encrypted signals are almost uncorrelated Corr→0, indicating good confusion effects. As the correlation coefficient of plaintext increases, the correlation coefficient of ciphertext decreases. This approach effectively spreads the influence of the current ciphertext as the encryption initial parameters onto the next set, continuously enhancing the confusion performance. The correlation diagrams of adjacent pixels in the original image and the encrypted signal are shown in Fig. 14.Fig. 14 The distributions of adjacent pixels in the original image and encrypted image.

Figure 14 illustrates the distribution of neighboring pixels in both the original and encrypted images. In the original image, pixels are clustered along a line, indicating a strong correlation between neighboring pixels. In contrast, the encrypted image displays a uniform pixel distribution across the entire range. This uniformity results from the disruption and alteration of the original image’s structure during encryption, which effectively conceals patterns and correlations. As a result, statistical attacks are unable to extract meaningful information from the encrypted image, thereby enhancing the security and robustness of the proposed encryption scheme.

The purpose of image encryption is to disrupt the strong correlation between neighboring pixels in the original image. In this study, 256×256 randomly selected pixels were analyzed to compute the image correlations. Table 10 presents the calculated correlation coefficients before and after encryption. As shown in the table, the correlation between neighboring pixels is greatly diminished post-encryption, resulting in images with minimal detectable patterns or regularities.Table 10 Correlation coefficient values of different schemes.

No	Encryption scheme	Test images	Correlation coefficient values	
Original image	Encrypted image	
1	Proposed	Flower	0.9853	− 0.0191	
Onion	0.9959	− 0.0004	
2	Ref. 28	Peppers	0.9581	− 0.0186	
3	Ref. 43	Lena	0.9226	0.0350	
4	Ref. 44	Baboon	0.8172	0.1093	
5	Ref. 45	Baboon	0.8172	− 0.0163	
6	Ref. 46	Lena	0.9226	− 0.0250	
7	Ref. 47	Lena	0.9226	0.0830	

This section evaluates the performance of an image encryption scheme by analyzing key aspects such as key space, histogram similarity, information entropy, resistance to statistical and differential attacks, key sensitivity, and pixel correlation. The results demonstrate that the scheme effectively disrupts pixel distribution and reduces the correlation between neighboring pixels, exhibiting strong randomness and resistance to attacks, thereby ensuring the security and robustness of the encrypted images.

Conclusions

This paper introduces a novel 4D hyperchaotic system and comprehensively analyzes its fundamental dynamic behaviors, encompassing equilibrium point stability, Lyapunov exponent spectrum, bifurcation diagram, and chaotic attractors. The study reveals that the new hyperchaotic system features a single equilibrium point and can maintain two positive Lyapunov exponents across a wide parameter range, with dynamic behavior ranging hyperchaotic, chaotic, periodic, and quasi-periodic, by adjusting the parameter d. Experimental implementation on embedded hardware STM32 validates these dynamics, with oscilloscope displays illustrating phase trajectories, affirming the system's effectiveness.

Moreover, the paper integrates the designed hyperchaotic system with embedded hardware STM32 for image encryption. This encryption algorithm employs the hyperchaotic key sequence generated by the system, with each image pixel undergoing a bit-XOR operation with corresponding bits of the key sequence. Despite its simple implementation, the method demonstrates high algorithmic complexity. Experimental results show that all hyperchaotic key sequences pass the 15 randomness tests of the NIST SP800-22 standard, confirming their complete randomness. Key space analysis, histogram similarity analysis, information entropy analysis, statistical attack analysis, differential attack analysis, key sensitivity analysis, and correlation analysis further validates the effectiveness of the hyperchaotic key-based encryption algorithm for image encryption.

In future research, the aim is to leverage the highly random and complex nature of hyperchaotic key sequences, in conjunction with multiple embedded hardware STM32 units, to achieve secure communication for data encryption/decryption across various domains, such as audio encryption, video encryption, text encryption, IoT security, financial transaction security, and medical data protection.

Acknowledgements

This work was supported by Science and Technology Project of Chongqing Municipal Education Commission (Grant No. KJZD-K202303301, KJQN202203309 and KJQN202203308).

Author contributions

Conceptualization, X.F.C.; methodology, X.F.C. and J.L.; software, X.F.C.; validation, H.M.Z., L.L. and J.L.; formal analysis, X.F.C.; investigation, K.P.M. and J.L.; resources, X.F.C. and J.L.; data curation, H.M.Z. and L.L.; writing—original draft preparation, X.F.C.; writing—review and editing, J.L.; supervision, X.F.C.; project administration, X.F.C.; funding acquisition, X.F.C. and H.M.Z. All authors reviewed the manuscript.

Data availability

The data used to support the findings of this study are available from the corresponding author upon request.

Competing interests

The authors declare no competing interests.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

These authors contributed equally: XueFeng Cheng and Juan Liu.
==== Refs
References

1. Rössler OE An equation for continuous chaos Chaos Interdiscip. J. Nonlinear Sci. 1976 57 397 398 10.1016/0375-9601(76)90101-8
Rössler, O. E. An equation for continuous chaos. Phys. Lett. A 57 397–398. 10.1016/0375-9601(76)90101-8 (1976).10.1016/0375-9601(76)90101-8
2. Parker JP Ashtari O Schneider TM Predicting chaotic statistics with unsfig invariant tori Chaos Interdiscip. J. Nonlinear Sci. 2023 33 083111 10.1063/5.0143689
Parker, J. P., Ashtari, O. & Schneider, T. M. Predicting chaotic statistics with unsfig invariant tori. Chaos Interdiscip. J. Nonlinear Sci. 33, 083111. 10.1063/5.0143689 (2023).10.1063/5.0143689
3. Yu F Dynamic analysis and FPGA implementation of a new, simple 5D memristive hyperchaotic Sprott-C system Mathematics 2023 11 701 10.3390/math11030701
Yu, F. et al. Dynamic analysis and FPGA implementation of a new, simple 5D memristive hyperchaotic Sprott-C system. Mathematics 11, 701. 10.3390/math11030701 (2023).10.3390/math11030701
4. Liu Y Zhou Y Guo B Hopf bifurcation, periodic solutions, and control of a new 4D hyperchaotic system Mathematics 2023 11 2699 10.3390/math11122699
Liu, Y., Zhou, Y. & Guo, B. Hopf bifurcation, periodic solutions, and control of a new 4D hyperchaotic system. Mathematics 11, 2699. 10.3390/math11122699 (2023).10.3390/math11122699
5. Cui N Li J A new 4D hyperchaotic system and its control AIMS Math. 2023 8 905 923 10.3934/math.2023044
Cui, N. & Li, J. A new 4D hyperchaotic system and its control. AIMS Math. 8, 905–923. 10.3934/math.2023044 (2023).10.3934/math.2023044
6. Li J Cui N Dynamical behavior and control of a new hyperchaotic Hamiltonian system AIMS Math. 2022 7 5117 5132 10.3934/math.2022285
Li, J. & Cui, N. Dynamical behavior and control of a new hyperchaotic Hamiltonian system. AIMS Math. 7, 5117–5132. 10.3934/math.2022285 (2022).10.3934/math.2022285
7. Lin L Zhuang Y Xu Z Yang D Wu D Encryption algorithm based on fractional order chaotic system combined with adaptive predefined time synchronization Front. Phys. 2023 11 1202871 10.3389/fphy.2023.1202871
Lin, L., Zhuang, Y., Xu, Z., Yang, D. & Wu, D. Encryption algorithm based on fractional order chaotic system combined with adaptive predefined time synchronization. Front. Phys. 11, 1202871. 10.3389/fphy.2023.1202871 (2023).10.3389/fphy.2023.1202871
8. Karawia A Cryptographic algorithm using newton-raphson method and general bischi-naimzadah duopoly system Entropy 2021 23 57 10.3390/e23010057
Karawia, A. Cryptographic algorithm using newton-raphson method and general bischi-naimzadah duopoly system. Entropy 23, 57. 10.3390/e23010057 (2021).10.3390/e23010057
9. Chen TH Yang CH Region of interest encryption based on novel 2D hyperchaotic signal and bagua coding algorithm IEEE Access 2022 10 82751 82765 10.1109/ACCESS.2022.3190851
Chen, T. H. & Yang, C. H. Region of interest encryption based on novel 2D hyperchaotic signal and bagua coding algorithm. IEEE Access 10, 82751–82765. 10.1109/ACCESS.2022.3190851 (2022).10.1109/ACCESS.2022.3190851
10. Fu SM Cheng XF Liu J Dynamics, circuit design, feedback control of a new hyperchaotic system and its application in audio encryption Sci. Rep. 2023 13 19385 10.1038/s41598-023-46161-5 37938603
Fu, S. M., Cheng, X. F. & Liu, J. Dynamics, circuit design, feedback control of a new hyperchaotic system and its application in audio encryption. Sci. Rep. 13, 19385. 10.1038/s41598-023-46161-5 (2023).37938603 10.1038/s41598-023-46161-5
11. Cao H Chu R Cui Y Complex dynamical characteristics of the fractional-order cellular neural network and its DSP implementation Fractal Fract. 2023 7 633 10.3390/fractalfract7080633
Cao, H., Chu, R. & Cui, Y. Complex dynamical characteristics of the fractional-order cellular neural network and its DSP implementation. Fractal Fract. 7, 633. 10.3390/fractalfract7080633 (2023).10.3390/fractalfract7080633
12. Li X Mou J Banerjee S Wang Z Cao Y Design and DSP implementation of a fractional-order detuned laser hyperchaotic circuit with applications in image encryption Chaos Solitons Fractals 2022 159 11233 10.1016/j.chaos.2022.112133
Li, X., Mou, J., Banerjee, S., Wang, Z. & Cao, Y. Design and DSP implementation of a fractional-order detuned laser hyperchaotic circuit with applications in image encryption. Chaos Solitons Fractals 159 , 112133. 10.1016/j.chaos.2022.112133 (2022).10.1016/j.chaos.2022.112133
13. Jia SH Li YX Shi QY Huang X Design and FPGA implementation of a memristor-based multi-scroll hyperchaotic system Chin. Phys. B 2022 31 070505 10.1088/1674-1056/ac4a71
Jia, S. H., Li, Y. X., Shi, Q. Y. & Huang, X. Design and FPGA implementation of a memristor-based multi-scroll hyperchaotic system. Chin. Phys. B 31, 070505. 10.1088/1674-1056/ac4a71 (2022).10.1088/1674-1056/ac4a71
14. Wang Y FPGA-based implementation and synchronization design of a new five-dimensional hyperchaotic system Entropy 2022 24 1179 10.3390/e24091179 36141065
Wang, Y. et al. FPGA-based implementation and synchronization design of a new five-dimensional hyperchaotic system. Entropy 24, 1179. 10.3390/e24091179 (2022).36141065 10.3390/e24091179
15. Babu NR Kalpana M Balasubramaniam P A novel audio encryption approach via finite-time synchronization of fractional order hyperchaotic system Multimed. Tools Appl. 2021 80 18043 18067 10.1007/s11042-020-10288-8
Babu, N. R., Kalpana, M. & Balasubramaniam, P. A novel audio encryption approach via finite-time synchronization of fractional order hyperchaotic system. Multimed. Tools Appl. 80, 18043–18067. 10.1007/s11042-020-10288-8 (2021).10.1007/s11042-020-10288-8
16. Vaidyanathan S A new 4-D multi-stable hyperchaotic system with no balance point: Bifurcation analysis, circuit simulation, FPGA realization and image cryptosystem IEEE Access 2021 9 144555 144573 10.1109/ACCESS.2021.3121428
Vaidyanathan, S. et al. A new 4-D multi-stable hyperchaotic system with no balance point: Bifurcation analysis, circuit simulation, FPGA realization and image cryptosystem. IEEE Access 9, 144555–144573. 10.1109/ACCESS.2021.3121428 (2021).10.1109/ACCESS.2021.3121428
17. Hou W Li S He J Ma Y A novel image-encryption scheme based on a non-linear cross-coupled hyperchaotic system with the dynamic correlation of plaintext pixels Entropy 2020 22 779 10.3390/e22070779 33286550
Hou, W., Li, S., He, J. & Ma, Y. A novel image-encryption scheme based on a non-linear cross-coupled hyperchaotic system with the dynamic correlation of plaintext pixels. Entropy 22, 779. 10.3390/e22070779 (2020).33286550 10.3390/e22070779
18. Wang L Chen Z Sun X He C Color image ROI encryption algorithm based on a novel 4D hyperchaotic system Phys. Scr. 2024 99 015229 10.1088/1402-4896/ad14d1
Wang, L., Chen, Z., Sun, X. & He, C. Color image ROI encryption algorithm based on a novel 4D hyperchaotic system. Phys. Scr. 99, 015229. 10.1088/1402-4896/ad14d1 (2024).10.1088/1402-4896/ad14d1
19. Nguyen QD Pham QD Thanh NT Giap VN An optimal homogenous stability-based disturbance observer and sliding mode control for secure communication system IEEE Access 2023 11 27317 27329 10.1109/ACCESS.2023.3257854
Nguyen, Q. D., Pham, Q. D., Thanh, N. T. & Giap, V. N. An optimal homogenous stability-based disturbance observer and sliding mode control for secure communication system. IEEE Access 11, 27317–27329. 10.1109/ACCESS.2023.3257854 (2023).10.1109/ACCESS.2023.3257854
20. Rybin V Prototyping the symmetry-based chaotic communication system using microcontroller unit Appl. Sci. 2023 13 936 10.3390/app13020936
Rybin, V. et al. Prototyping the symmetry-based chaotic communication system using microcontroller unit. Appl. Sci. 13, 936. 10.3390/app13020936 (2023).10.3390/app13020936
21. Wang P Secure transmission for IoT wireless energy-carrying communication systems PLOS ONE 2023 18 e0289251 10.1371/journal.pone.0289251 37535589
Wang, P. et al. Secure transmission for IoT wireless energy-carrying communication systems. PLOS ONE 18, e0289251. 10.1371/journal.pone.0289251 (2023).37535589 10.1371/journal.pone.0289251
22. Wang M Niu Y Gao B Zou Q Hyperchaotic impulsive synchronization and digital secure communication J. Appl. Math. Phys. 2022 10 3485 3495 10.4236/jamp.2022.1012230
Wang, M., Niu, Y., Gao, B. & Zou, Q. Hyperchaotic impulsive synchronization and digital secure communication. J. Appl. Math. Phys. 10, 3485–3495. 10.4236/jamp.2022.1012230 (2022).10.4236/jamp.2022.1012230
23. He J Qiu W Cai J Synchronization of hyperchaotic systems based on intermittent control and its application in secure communication J. Adv. Comput. Intell. Intell. Inform. 2023 27 292 303 10.20965/jaciii.2023.p0292
He, J., Qiu, W. & Cai, J. Synchronization of hyperchaotic systems based on intermittent control and its application in secure communication. J. Adv. Comput. Intell. Intell. Inform. 27, 292–303. 10.20965/jaciii.2023.p0292 (2023).10.20965/jaciii.2023.p0292
24. Alexan W Chen YL Por LY Gabr M Hyperchaotic maps and the single neuron model: A novel framework for chaos-based image encryption Symmetry 2023 15 1081 10.3390/sym15051081
Alexan, W., Chen, Y. L., Por, L. Y. & Gabr, M. Hyperchaotic maps and the single neuron model: A novel framework for chaos-based image encryption. Symmetry 15, 1081. 10.3390/sym15051081 (2023).10.3390/sym15051081
25. Zhu S Deng X Zhang W Zhu C Construction of a new 2D hyperchaotic map with application in efficient pseudo-random number generator design and color image encryption Mathematics 2023 11 3171 10.3390/math11143171
Zhu, S., Deng, X., Zhang, W. & Zhu, C. Construction of a new 2D hyperchaotic map with application in efficient pseudo-random number generator design and color image encryption. Mathematics 11, 3171. 10.3390/math11143171 (2023).10.3390/math11143171
26. Sun S A new image encryption scheme based on 6D hyperchaotic system and random signal insertion IEEE Access 2023 11 66009 66016 10.1109/ACCESS.2023.3290915
Sun, S. A new image encryption scheme based on 6D hyperchaotic system and random signal insertion. IEEE Access 11 , 66009–66016. 10.1109/ACCESS.2023.3290915 (2023).10.1109/ACCESS.2023.3290915
27. Shen Y Fast and secure image encryption algorithm with simultaneous shuffling and diffusion based on a time-delayed combinatorial hyperchaos map Entropy 2023 25 753 10.3390/e25050753 37238508
Shen, Y. et al. Fast and secure image encryption algorithm with simultaneous shuffling and diffusion based on a time-delayed combinatorial hyperchaos map. Entropy 25, 753. 10.3390/e25050753 (2023).37238508 10.3390/e25050753
28. Gao X Sun B Cao Y Banerjee S Mou J A color image encryption algorithm based on hyperchaotic map and DNA mutation Chin. Phys. B 2021 32 030501 10.1088/1674-1056/ac8cdf
Gao, X., Sun, B., Cao, Y., Banerjee, S. & Mou, J. A color image encryption algorithm based on hyperchaotic map and DNA mutation. Chin. Phys. B 32 , 030501. 10.1088/1674-1056/ac8cdf (2023).10.1088/1674-1056/ac8cdf
29. Sun S Guo Y A new hyperchaotic image encryption algorithm based on stochastic signals IEEE Access 2021 9 144035 144045 10.1109/ACCESS.2021.3121588
Sun, S. & Guo, Y. A new hyperchaotic image encryption algorithm based on stochastic signals. IEEE Access 9, 144035–144045. 10.1109/ACCESS.2021.3121588 (2021).10.1109/ACCESS.2021.3121588
30. Hou W Li S He J Ma Y A new hyperchaotic image encryption algorithm based on stochastic signals Entropy 2020 22 779 10.3390/e22070779 33286550
Hou, W., Li, S., He, J. & Ma, Y. A novel image-encryption scheme based on a non-linear cross-coupled hyperchaotic system with the dynamic correlation of plaintext pixels. Entropy 22, 779. 10.3390/e22070779 (2020).33286550 10.3390/e22070779
31. Elsonbaty A Elsadany AA Adel W On reservoir computing approach for digital image encryption and forecasting of hyperchaotic finance model Fractal Fract. 2023 7 282 10.3390/fractalfract7040282
Elsonbaty, A., Elsadany, A. A. & Adel, W. On reservoir computing approach for digital image encryption and forecasting of hyperchaotic finance model. Fractal Fract. 7, 282. 10.3390/fractalfract7040282 (2023).10.3390/fractalfract7040282
32. Du Y Long G Jiang D Chai X Han J Optical image encryption algorithm based on a new four-dimensional memristive hyperchaotic system and compressed sensing Chin. Phys. B 2023 32 114203 10.1088/1674-1056/acef08
Du, Y., Long, G., Jiang, D., Chai, X. & Han, J. Optical image encryption algorithm based on a new four-dimensional memristive hyperchaotic system and compressed sensing. Chin. Phys. B10.1088/1674-1056/acef08 (2023).10.1088/1674-1056/acef08
33. Chen W Wang Y Xiao Y Hei X Explore the potential of deep learning and hyperchaotic map in the meaningful visual image encryption scheme IET Image Process. 2023 17 3235 3257 10.1049/ipr2.12858
Chen, W., Wang, Y., Xiao, Y. & Hei, X. Explore the potential of deep learning and hyperchaotic map in the meaningful visual image encryption scheme. IET Image Process. 17, 3235–3257. 10.1049/ipr2.12858 (2023).10.1049/ipr2.12858
34. Xu H Wang J New 4D hyperchaotic system’s application in image encryption J. Opt. 2024 26 065503 10.1088/2040-8986/ad3e0d
Xu, H. & Wang, J. New 4D hyperchaotic system’s application in image encryption. J. Opt. 26, 065503. 10.1088/2040-8986/ad3e0d (2024).10.1088/2040-8986/ad3e0d
35. Ding L Ding Q The establishment and dynamic properties of a new 4D hyperchaotic system with its application and statistical tests in gray images Entropy 2020 22 310 10.3390/e22030310 33286085
Ding, L. & Ding, Q. The establishment and dynamic properties of a new 4D hyperchaotic system with its application and statistical tests in gray images. Entropy 22, 310. 10.3390/e22030310 (2020).33286085 10.3390/e22030310
36. Jiang Q Yu S Wang Q Cryptanalysis of an image encryption algorithm based on two-dimensional hyperchaotic map Entropy 2023 25 395 10.3390/e25030395 36981284
Jiang, Q., Yu, S. & Wang, Q. Cryptanalysis of an image encryption algorithm based on two-dimensional hyperchaotic map. Entropy 25, 395. 10.3390/e25030395 (2023).36981284 10.3390/e25030395
37. Wen H Secure DNA-coding image optical communication using non-degenerate hyperchaos and dynamic secret-key Mathematics 2022 10 3180 10.3390/math10173180
Wen, H. et al. Secure DNA-coding image optical communication using non-degenerate hyperchaos and dynamic secret-key. Mathematics 10, 3180. 10.3390/math10173180 (2022).10.3390/math10173180
38. Liu J Cheng X Zhou P Circuit implementation synchronization between two modified fractional-order Lorenz Chaotic systems via a linear resistor and fractional-order capacitor in parallel coupling Math. Probl. Eng. 2021 2021 1 8 10.1155/2021/6771261
Liu, J., Cheng, X. & Zhou, P. Circuit implementation synchronization between two modified fractional-order Lorenz Chaotic systems via a linear resistor and fractional-order capacitor in parallel coupling. Math. Probl. Eng. 2021, 1–8. 10.1155/2021/6771261 (2021).10.1155/2021/6771261
39. Chen G Mao Y Chui CK A symmetric image encryption scheme based on 3D chaotic cat maps Chaos Solitons Fractals 2004 21 749 761 10.1016/j.chaos.2003.12.022
Chen, G., Mao, Y. & Chui, C. K. A symmetric image encryption scheme based on 3D chaotic cat maps. Chaos Solitons Fractals 21, 749–761. 10.1016/j.chaos.2003.12.022 (2004).10.1016/j.chaos.2003.12.022
40. Luo Y Du M Liu J A symmetrical image encryption scheme in wavelet and time domain Commun. Nonlinear Sci. Numer. Simul. 2015 20 447 460 10.1016/j.cnsns.2014.05.022
Luo, Y., Du, M. & Liu, J. A symmetrical image encryption scheme in wavelet and time domain. Commun. Nonlinear Sci. Numer. Simul. 20, 447–460. 10.1016/j.cnsns.2014.05.022 (2015).10.1016/j.cnsns.2014.05.022
41. Wang X-Y Zhang Y-Q Zhao Y-Y A novel image encryption scheme based on 2-D logistic map and DNA sequence operations Nonlinear Dyn. 2015 82 1269 1280 10.1007/s11071-015-2234-7
Wang, X.-Y., Zhang, Y.-Q. & Zhao, Y.-Y. A novel image encryption scheme based on 2-D logistic map and DNA sequence operations. Nonlinear Dyn. 82, 1269–1280. 10.1007/s11071-015-2234-7 (2015).10.1007/s11071-015-2234-7
42. Liu W Sun K He Y Yu M Color image encryption using three-dimensional sine ICMIC modulation map and DNA sequence operations Int. J. Bifurc. Chaos 2017 27 1750171 10.1142/S0218127417501711
Liu, W., Sun, K., He, Y. & Yu, M. Color image encryption using three-dimensional sine ICMIC modulation map and DNA sequence operations. Int. J. Bifurc. Chaos 27, 1750171. 10.1142/S0218127417501711 (2017).10.1142/S0218127417501711
43. Kaur G Agarwal R Patidar V Color image encryption scheme based on fractional Hartley transform and chaotic substitution–permutation Vis. Comput. 2022 38 1027 1050 10.1007/s00371-021-02066-w
Kaur, G., Agarwal, R. & Patidar, V. Color image encryption scheme based on fractional Hartley transform and chaotic substitution–permutation. Vis. Comput. 38, 1027–1050. 10.1007/s00371-021-02066-w (2022).10.1007/s00371-021-02066-w
44. Kaur G Agarwal R Patidar V Color image encryption system using combination of robust chaos and chaotic order fractional Hartley transformation J. King Saud Univ. Comput. Inf. Sci. 2022 34 5883 5897 10.1016/j.jksuci.2021.03.007
Kaur, G., Agarwal, R. & Patidar, V. Color image encryption system using combination of robust chaos and chaotic order fractional Hartley transformation. J. King Saud Univ. Comput. Inf. Sci. 34, 5883–5897. 10.1016/j.jksuci.2021.03.007 (2022).10.1016/j.jksuci.2021.03.007
45. Zhang D Chen L Li T Hyper-chaotic color image encryption based on transformed zigzag diffusion and RNA operation Entropy 2021 23 361 10.3390/e23030361 33802901
Zhang, D., Chen, L. & Li, T. Hyper-chaotic color image encryption based on transformed zigzag diffusion and RNA operation. Entropy 23, 361. 10.3390/e23030361 (2021).33802901 10.3390/e23030361
46. Chen C Sun K Xu Q A color image encryption algorithm based on 2D-CIMM chaotic map China Commun. 2020 17 12 20 10.23919/JCC.2020.05.002
Chen, C., Sun, K. & Xu, Q. A color image encryption algorithm based on 2D-CIMM chaotic map. China Commun. 17, 12–20. 10.23919/JCC.2020.05.002 (2020).10.23919/JCC.2020.05.002
47. Farah MAB Guesmi R Kachouri A Samet M A novel chaos based optical image encryption using fractional Fourier transform and DNA sequence operation Opt. Laser Technol. 2020 121 105777 10.1016/j.optlastec.2019.105777
Farah, M. A. B., Guesmi, R., Kachouri, A. & Samet, M. A novel chaos based optical image encryption using fractional Fourier transform and DNA sequence operation. Opt. Laser Technol. 121, 105777. 10.1016/j.optlastec.2019.105777 (2020).10.1016/j.optlastec.2019.105777
