Next Article in Journal
A Low-Computational Burden Closed-Form Approximated Expression for MSE Applicable for PTP with gfGn Environment
Next Article in Special Issue
Fractal-Fractional-Order Modeling of Liver Fibrosis Disease and Its Mathematical Results with Subinterval Transitions
Previous Article in Journal
Existence, Uniqueness, and Stability of a Nonlinear Tripled Fractional Order Differential System
Previous Article in Special Issue
Fractional Operators and Fractionally Integrated Random Fields on Zν
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Application of Fractional-Order Multi-Wing Chaotic System to Weak Signal Detection

Physics Department, Changchun University of Science and Technology, Changchun 130022, China
*
Author to whom correspondence should be addressed.
Fractal Fract. 2024, 8(7), 417; https://doi.org/10.3390/fractalfract8070417
Submission received: 27 May 2024 / Revised: 11 July 2024 / Accepted: 13 July 2024 / Published: 16 July 2024

Abstract

:
This work investigates a fractional-order multi-wing chaotic system for detecting weak signals. The influence of the order of fractional calculus on chaotic systems’ dynamical behavior is examined using phase diagrams, bifurcation diagrams, and SE complexity diagrams. Then, the principles and methods for determining the frequencies and amplitudes of weak signals are examined utilizing fractional-order multi-wing chaotic systems. The findings indicate that the lowest order at which this kind of fractional-order multi-wing chaotic system appears chaotic is 2.625 at a = 4 , b = 8 , and c = 1 , and that this value decreases as the driving force increases. The four-wing and double-wing change dynamics phenomenon will manifest in a fractional-order chaotic system when the order exceeds the lowest order. This phenomenon can be utilized to detect weak signal amplitudes and frequencies because the system parameters control it. A detection array is built to determine the amplitude using the noise-resistant properties of both four-wing and double-wing chaotic states. Deep learning images are then used to identify the change in the array’s wing count, which can be used to determine the test signal’s amplitude. When frequencies detection is required, the MUSIC method estimates the frequencies using chaotic synchronization to transform the weak signal’s frequencies to the synchronization error’s frequencies. This solution adds to the contact between fractional-order calculus and chaos theory. It offers suggestions for practically implementing the chaotic weak signal detection theory in conjunction with deep learning.

1. Introduction

As a nonlinear signal detection method, chaotic weak signal detection has been studied for over thirty years. The difference between the chaotic system method and other classical signal processing techniques is that chaotic systems can resist noise through iteration. This allows the noise to be much larger than weak signals and achieve a negative signal-to-noise ratio, which is already superior to the most common amplification methods. Among numerous weak signal detection methods that can handle noise, the detection performance of chaotic detection methods is related to the design of chaotic systems. By increasing the effective number of control parameters in chaotic systems, weaker, smaller signals can be detected, resulting in higher detection performance limits. Since Birx [1] proved in 1992 that the Duffing chaotic system can be used for weak signal detection, he is considered the pioneer of the theory of chaotic weak signal detection. Wang [2] continued his work on Birx by showing how to detect weak signals using the initial value sensitivity property of the Duffing system and gave a complete procedure for doing so. In recent years, researchers have kept a close eye on Duffing systems for weak signal detection, and several novel theories for gauging the signal detection capability have emerged [3,4,5,6,7]. Meanwhile, many practical problems have been resolved by applying chaotic weak signal detection approaches [8,9,10]. This shows that the subject of chaos-based weak signal identification has gained popularity in engineering applications and can potentially tackle problems that are not well-suited for linear signal detection.
Researchers have employed more chaotic systems to detect signals and have also found solutions to a few real-world issues to enhance the detection performance further [11,12,13,14,15,16,17,18]. Furthermore, the multi-wing chaotic system has a distinct attractor topology beyond these additional systems. As a result, this chaotic system likewise displays extremely complicated features [19,20,21,22,23,24], which naturally prompts researchers to think of using it for the detection of weak signals. Li [25] created a nonlinear feedback controller from chaotic synchronization, stabilized the chaotic system with many wings to the equilibrium point, and then used the chaotic system’s synchronized state to detect weak signals with multiple frequencies. In Yan’s [26] study of a multi-wing chaotic system with an infinite number of equilibrium points, he found that weak signals might be identified by comparing the fluctuations between four-wing and two-wing chaotic attractors. The viability and superiority of employing a multi-wing chaotic system for weak signal identification are amply illustrated by these works.
Due to advancements in chaos theory, fractional-order chaotic systems have grown in importance within the field of chaos. Additionally, fractional-order calculus has introduced new features to the field, making fractional-order chaotic systems highly relevant as a research hotspot with significant potential for engineering applications [27,28,29,30,31]. In all chaotic applications, fractional-order chaotic systems for the detection of weak signals have emerged as a state-of-the-art part of chaotic weak signal detection [32,33,34,35,36].
Although fractional-order multi-wing systems have been partially studied and applied [37,38,39,40], research on how to use such systems for weak signal detection is rare, especially in demonstrating the advantages of fractional calculus. Secondly, how to apply the existing theory of chaotic weak signal detection to practical device design is also a major obstacle to the development of this theory. Finally, deep learning image recognition is a very mature technology, and in the state criteria of chaotic systems, the attractor graph method is a classic method. However, the development of combining the attractor graph method with deep learning image recognition and introducing deep learning into chaotic weak signal detection is still blank.
Therefore, in order to address the above issues, in this paper, we researched the application of fractional-order multi-wing chaotic systems in weak signal detection. Firstly, the multi-wing chaotic system will be rewritten into the fractional-order form, because the fractional-order system is more consistent with the actual physical model. Then the influence of order parameters in fractional-order calculus on the characteristics of chaotic systems is analyzed. Lastly, weak signal detection was accomplished using this fractional-order chaotic system. Using the difference in the number of four-wing and two-wing chaotic attractor wings, a chaotic array represents the difference between weak signals while detecting their amplitude. In determining the frequency of the weak signal, a fractional-order multi-wing chaotic synchronization system is designed which is built using the drive–response method. If the weak signal is added to the driving part, the synchronization error of this system will vary with different frequencies of the weak signal. This system then uses a multi-signal classification (MUSIC) algorithm to estimate the weak signal frequencies by converting the weak signal frequencies into a chaotic synchronization error.
The remainder of the paper is summarized as follows after the introduction: Section 2 presents the system model of the fractional-order multi-wing chaotic system and employs phase diagrams, equilibrium points, bifurcation diagrams, complexity, and other techniques to assess the dynamics features. Section 3 covers the identification of weak signals using a fractional-order multi-wing chaotic system. It divides the procedure into two parts: amplitude detection and frequency detection. It also provides the detection principle and discusses using the chaotic array and MUSIC methods for processing data in real-world scenarios. Section 4 concludes the paper at the end.

2. Fractional-Order Multi-Wing Chaotic System Model and Dynamic Characteristics Analysis

2.1. Fractional-Order Multi-Wing Chaotic System Model

A non-autonomous, multi-wing chaotic system was proposed by Yan [26]. It has complicated dynamical features, including the generation of symmetric attractors, an infinite number of equilibrium points, and a driving amplitude that affects the number of wings. The following is the equation:
x ˙ = a x + b y z y ˙ = c y x z z ˙ = z + x y + G
where G is the external driving force signal, a, b and c are positive parameters, x, y, and z are state variables, and the amplitude r and frequency ω are often included in G. As a result, this system can be employed as a weak signal detection system. Nevertheless, there are still issues with this detection system. The main one is that the frequency and amplitude r do not regulate the appearance of chaos in the numerical region of all standard weak signals, which leads to the creation of detection blind zones and lowers the detection performance.
Therefore, using fractional-order calculus, Equation (1) is enhanced in this study to create a new 3D fractional-order system that is more suited for weak signal detection. This new system can be stated as:
D q x = a x + b y z D q y = c y x z D q z = z + x y + G
where G is the external driving force signal, x, y, and z are the state variables, a, b, and c are the positive parameters, and q is the derivative order. D is the fractional-order differential established by Caputo. A few fundamental characteristics of system Equation (2) are examined next.

2.2. Dissipativity and the Existence of an Attractor

The divergence can be used to calculate the dissipativity of the chaotic equations [41], and the exponential constraint rate of system Equation (2) is:
Δ V = x ˙ x + y ˙ y + z ˙ z
For the proposed system, Δ V is a + c 1 , and for the chosen set of parameters, it is equal to 4 . This result means that all system orbits will eventually be confined to a subset of zero volume, i.e., system Equation (2) is dissipative. This dissipative nature of the system guarantees the existence of an attractor for the system.
So the numerical approximate solution of the fractional-order multi-wing chaotic system can be obtained through calculation, which can be used to draw the attractor phase diagram. For the approximate solution, we use the Grünwald–Letnikov approximation formula based on the Riemann–Liouville definition. The general numerical solution expression for fractional differential equations using this method is:
D t q a y ( t ) = f ( y ( t ) , t ) y ( t k ) = f ( y ( t k ) , t k ) × h q j = ν k c j ( q ) y ( t k j )
where
c 0 ( q ) = 1 , c j ( q ) = ( 1 1 + q j ) × c j 1 ( q ) .
According to the existing studies on chaotic systems of integer order, let q = 0.95 , a = 4 , b = 8 , c = 1 , G = r × s i n ( ω × t ) ; take r = 1.126 , and ω = 1 . The initial values are taken as [ x , y , z ] = [ 0.1 ; 0.1 ; 0.1 ] . This chaotic system forms a four-wing chaotic attractor as shown in Figure 1.

2.3. Equilibrium Analysis

Drive force signal G governs system Equation (2), which is non-autonomous. As a result, the left side of system Equation (2) must be set to zero to compute the equilibrium point.
0 = a x + b y z 0 = c y x z 0 = z + x y + G
Then, we may determine that, at this moment, the equilibrium points of Equations (1) and (2) are the same and that infinite equilibrium points exist with infinite driving signal G values. Consequently, as indicated in Table 1, they are grouped to obtain five different types of equilibrium points.
To get the Jacobian matrix, linearize system Equation (2) at the equilibrium points ( x * , y * , z * ):
J = 4 8 z * 8 y * z * 1 x * y * x * 1 .
Let d e t ( λ E J ) = 0 ; obtain the corresponding characteristic polynomial as:
λ 3 + A 2 λ 2 + A 1 λ + A 0 = 0 .
where
A 2 = 4 , A 1 = 1 8 ( y * ) 2 + 8 ( x * ) 2 + ( z * ) 2 , A 0 = 4 + 16 x * y * z * + 8 ( y * ) 2 + 4 ( x * ) 2 + 8 ( z * ) 2 .
The eigenvalues of the equilibrium point fulfill R e [ λ ] < 0 if it is stable. According to the Routh–Hurwitz criterion, the equilibrium point is only stable when A 2 > 0 , A 0 > 0 and A 2 × A 1 A 0 > 0 . Given the conditions of this paper, we select G = 0 , 1 , and 5, respectively, because the size of G impacts each equilibrium point, and G also affects the system’s ability to identify weak signals. Table 2 displays the equilibrium point’s stability.
According to Table 2, the equilibrium point S 0 belongs to the first class of saddle points, whereas the saddle points S 1 , S 2 , S 3 , and S 4 belong to the second class. The first class of saddle points connects the rings of chaotic attractors that were created around the second class of saddle points.
M. S. Tavazoei [42] provides the condition for the system to appear as a locally asymptotically stable equilibrium point, which can be computed to determine the minimal fractional calculus order in which chaos can arise. The results are as follows:
a r g ( e i g ( J ) ) < q π 2 .
As shown in Table 2, the necessary condition for the existence of chaotic attractors in system Equation (2) at a = 4 , b = 8 , and c = 1 is q > 0.875 , which means the total order is 2.625 , and this value decreases with the increase of G.

2.4. Dynamic Characteristics under the Influence of Order

The integer-order chaotic system can be considered a specific fractional-order system by the theory of fractional-order calculus. Except for order, the remaining parameters of a fractional-order chaotic system then influence the chaotic system in a way comparable to the law of an integer-order system. As a result, we solely consider how the additional order parameter affects the fractional-order multi-wing chaotic system’s dynamic properties. A bifurcation diagram and SE complexity diagram are used for this. The result is shown in Figure 2:
The chaotic phenomenon is evident in the range of [ 0.85 , 1.05 ] in Figure 2, which is consistent with the results computed in the preceding section. Furthermore, the chaotic system ceases to be chaotic when the order exceeds 1.05 , demonstrating how differing fractional-order calculus orders alter the remaining parameters’ control range. Thus, the advantage of the fractional-order chaotic system is demonstrated by the fact that, by altering the order, it can gain a more extensive chaos interval or period window than the integer-order chaotic system’s fixed influence range of each parameter.
Next, taking orders of 0.9 and 0.95 as examples, we will compare the influence of G values on the state of chaotic systems at different orders. The results are shown in Figure 3:
From Figure 3, it can be seen that after changing the order value from 1, the length of the chaotic region becomes significantly longer, and the number of occurrences of the period window becomes less. This suggests fewer blind spots for detection and a more comprehensive detection range for weak signal detection with this technique.

3. Weak Signal Detection

The study that is now available indicates that the multi-wing chaotic system Equation (1) has the property of having switchable numbers of wings, which is achieved by varying the driving force G. Consequently, this property of the fractional-order chaotic system remains. Furthermore, this chaotic system possesses the properties of a single-parameter chaotic system with rapid state change, which is compatible with the chaotic system properties required in the field of chaotic weak signal detection. Consequently, it is possible to convert the information of weak signals into fractional-order multi-wing chaotic systems’ dynamical phenomena and thereafter identify the weak signals.

3.1. Amplitudes Detection of Weak Signals

According to the conventional theory of chaotic weak signal detection, adding a weak signal will cause the chaotic system to enter a periodic state. At this time, the chaotic system will undergo chaotic control to revert to its chaotic state, offsetting the added weak signals and the altered control parameter to obtain the signal for measurement. In this process, the chaotic state has a higher noise resistance than the periodic state. Therefore, to circumvent the periodic state’s low robustness to noise, we used the difference between the four-winged and the two-winged chaotic states in this study to detect the weak signals.
First, we considered the case of a single fractional-order chaotic system as a detection system when the equations of the system are:
D q x = a x + b y z D q y = c y x z D q z = z + x y + G + M ( t )
where M ( t ) = g + N o s i e is the addition of a mixed signal made up of noise and the signal that has to be measured. A portion of the signal to be measured must be known in advance; g is the same form of the weak signal to be calculated as the drive signal. Considering Gaussian white noise as the noise and a sinusoidal weak signal with a frequency of 1 rad/s, M ( t ) represents the weak signal. The question is how to determine the amplitude of the weak signal at this particular instant. With q = 0.9 , a = 4 , b = 8 , c = 1 , and G = r   × sin ( t ) , the attractor of Equation (10) is displayed as follows when r is between the critical values of 1.12939796 and 1.12939797 .
Figure 4 shows that the attractor of the chaotic system Equation (10) will vary if the control parameters are altered by 0.00000001 . Consequently, it is possible to determine the variation pattern of control parameters by tracking changes in attractors. The influence of the signal to be tested is similar to altering the control parameters because both the form of G and the signal M ( t ) to be tested are consistent. This allows the attractor’s change to determine the amplitude of the weak signal that needs to be measured. Numerous studies have examined the effects of various noise on chaotic systems [43,44], and fractional-order multi-wing chaotic systems have the inherent ability to withstand noise.
In terms of the performance of individual chaotic detection, there is no essential difference between fractional-order chaotic systems and integer-order chaotic systems. However, when establishing a chaotic oscillator array, multiple arrays with order differences can be established by changing the fractional order. That is why this paper discusses the use of chaotic oscillator arrays for weak signal detection. It should be noted that fractional-order chaotic systems can use distinct fractional-order systems to form arrays, unlike integer-order systems, which can only modify one parameter. Consequently, the array’s dimensionality has increased. Once the array is established, variations can be used to determine the amplitude information of weak signals. Several arrays can be built to increase precision because this paper’s fractional-order multi-wing chaotic system can display chaotic events in fractional orders between 0.875 and 1. Here, taking the order equal to 0.9 as an example, we explain how to use a chaotic oscillator array for detection.
Figure 5 shows that the interval over which the four-wing chaos phenomenon occurs has a few complex windows. When utilizing an array to detect signals instead of a single chaotic oscillator, these windows can prevent errors from occurring. The chaotic array alters due to the addition of various weak signals, as seen in Figure 6.
Figure 6 illustrates how the addition of weak signals altered the chaotic array’s properties, with the direction of the change varying according to whether the signal was positive or negative. The chaotic array’s wing count will decrease with the addition of a positive weak signal and rise with the addition of a negative weak signal. Thus, in real-world applications, the measured signal’s amplitude can be acquired as long as the chaotic array is located.
From the above analysis, it can be seen that in order to use chaotic arrays to detect weak signals, it is necessary to identify the state of the chaotic array and obtain information on the changes in the chaotic array. If the common Lyapunov exponent method is used, each oscillator needs to be calculated, and the 4 × 4 size in this article requires 16 calculations. If the chaotic phase diagram method is used for discrimination, it reduces the calculation time and number of calculations, but there is no universally recognized criterion for qualitative analysis. Therefore, we consider using deep learning image recognition to recognize and distinguish chaotic phase diagrams.
The manuscript introduces deep learning image recognition into the theory of chaotic weak signal detection, which is used to quickly identify information in chaotic arrays. This method has been rare in previous research, and it uses the multi-wing features of the attractor to simultaneously identify the entire array, reducing the judgment time. In addition, the main problem of detecting weak signals with chaos is the need to judge the chaotic state, so currently, commonly used methods make it difficult to design actual signal detection devices. By establishing a chaotic array to transform weak signals into image changes in a chaotic array and then using more mature deep learning image recognition for processing, this design makes the emergence of actual signal detection devices possible.
There are many mature models for deep learning image recognition, and the T-Rex2 model is used in this manuscript. It is open-source and can be used on online websites. T-Rex2 is an interactive object counting model designed to first detect and then count any objects. It formulates object counting as an open-set object detection task with the integration of visual prompts. Users can specify the objects of interest by marking points or boxes on a reference image, and T-Rex2 then detects all objects with a similar pattern. Guided by the visual feedback from T-Rex2, users can also interactively refine the counting results by identifying missing or falsely detected objects.
As a result, we may learn the attractor graph based on two wings by using Figure 4 as sample data. Equation (7) presents the findings.
From Figure 7, it can be seen that the deep learning image recognition method can successfully obtain the number of wings of the chaotic attractor based on the number of samples. So the number of wings of the chaotic arrays in Figure 5 and Figure 6 are 21, 23, and 20, respectively. It is possible to determine the precise amplitude of the extra-weak signal by comparing the wing numbers. The actual amplitudes of the weak signals can be determined to be 0.01 × 2 = 0.02 and 0.01 × 1 = 0.01 , respectively, based on the order of magnitude involved in building a chaotic oscillator array.
Next, we will consider the detection performance of fractional-order multi-wing chaotic systems under the conditions selected for this paper. The lowest detectable signal-to-noise ratio is used to represent the performance, with the specific condition that when the noise is too large and the attractor phase diagram changes to the point where deep learning cannot recognize it, it is considered undetectable. When the noise levels are 0.329 and 0.330 , respectively, the results of deep learning image recognition using chaotic attractors are shown in Figure 8.
So according to the formula, the minimum detectable signal-to-noise ratio of a fractional-order multi-wing chaotic system under the selected conditions is:
S N R = 10 l o g ( 0.01 / 0.329 ) = 34.93 d b .

3.2. Frequencies Detection of Weak Signals

Currently, frequencies detection can be considered using the sweeping method, but this method is time-consuming and consumes many hardware resources. Therefore, researchers have proposed using the method of chaotic synchronization to detect the frequencies of weak signals. This is because the current research on chaos synchronization has been very in-depth, and many different synchronization methods have been proposed [45,46]. The drive–response method is a very classic chaos synchronization method. The use of the drive–response chaotic system to detect the frequencies of weak signals is based on the principle of using an applied weak signal to change the synchronization error of the chaotic system, and this error is linearly correlated with the applied weak signal. Therefore, the frequencies of the weak signal can be obtained by calculating the frequencies of the error value. The designed driving and controlled response systems are shown in Equations (11) and (12).
D q x = a x 1 + b y 1 z 1 D q y = c y 1 x 1 z 1 D q z = z 1 + x 1 y 1 + G
D q x = a x 2 + b y 2 z 2 u x D q y = c y 2 x 2 z 2 u y D q z = z 2 + x 2 y 2 + G u z
The terms u x , u y , and u z in system Equation (12) are nonlinear controllers. Let the synchronization error between system Equation (11) and system Equation (12) be e x = x 2 x 1 , e y = y 2 y 1 , and e z = z 2 z 1 . Then the synchronization error is the difference between Equations (11) and (12). That is, the following synchronization error of fractional-order chaotic systems is obtained:
D q e x = a e x + b ( e y e z + e y z 1 + e z y 1 ) u x D q e y = c e y ( e x e z + e x z 1 + e z x 1 ) u y D q e z = e z + e x e y + e x y 1 + e y x 1 u z
The synchronization controller is designed as follows:
u x = b ( e y e z + e y z 1 + e z y 1 ) u y = c e y ( e x e z + e x z 1 + e z x 1 ) + k y e y u z = e x e y + e x y 1 + e y x 1
Bringing system Equation (14) into system Equation (13) yields the final synchronization error, as shown in system Equation (15).
D q e x = a e x D q e y = k y e y D q e z = e z
Since the fractional-order order 0 < q < 1 in this paper, the Lyapunov function E is constructed according to the sufficient condition for stability of fractional-order systems [47,48]:
E ( e x , e y , e z ) = 1 2 ( e x 2 + e y 2 + e z 2 ) .
The derivative is then obtained as:
E ˙ = e x ˙ e x + e y ˙ e y + e z ˙ e z = a e x 2 k e y 2 e z 2 .
It is obvious that the synchronization error is asymptotically stable; thus, we are able to show that the drive and response systems Equations (11) and (12) are synchronized, as shown in Figure 9.
Therefore, synchronization errors can be used to observe changes in the drive–response system when detecting the signal’s frequencies. For instance, the equation for driving the system becomes as follows when a mixed signal, S ( t ) , is added to the third component in system Equation (11):
D q x = a x + b y z D q y = c y x z D q z = z + x y + G + S ( t )
Thus, S ( t ) will enter the chaotic system’s iteration in system (16), and S ( t ) will not change the chaotic system’s dynamic properties. On the other hand, if S ( t ) has the same form as G, then altering the system’s parameters equals modifying the chaotic system’s dynamic properties. As seen in Figure 10, this will modify the chaotic synchronization system’s synchronization error.
Figure 10 illustrates this point. While the chaotic system is stable, varying S ( t ) will result in varying synchronization errors, and the frequencies of fluctuations in synchronization errors is equal to that of S ( t ) . Thus, it is possible to determine the frequencies of S ( t ) by computing the frequencies of synchronization mistakes. Additionally, a weak signal frequencies was found. The frequencies inside the synchronization error can be obtained in various ways. Schmidt [49] proposed the MUSIC algorithm, which is used in this paper. This algorithm can assess the signal’s frequencies and successfully withstand noise. The following are the primary steps:
1. Use N data points stabilized by the synchronization error, then extract N M data points from them as the sample matrix.
2. Find the covariance matrix of the constructed sample. R = 1 N M n = 1 N M X ( n ) X H ( n ) .
3. Perform eigenvalue decomposition of the covariance matrix R. The eigenvalues are arranged from small to large, and the first m corresponding eigenvectors are taken to construct the signal eigenmatrix U S = [ v 1 , v 2 , v m ] . The rest are used to construct the noise characteristic matrix U N = [ v m + 1 , v m + 2 , v M ] .
4. Using the frequencies estimation formula, P ( f ) = 1 G ( f ) U N U N G ( f ) , where G ( f ) = [ 1 , e 2 π f j , , e 2 π ( M 1 ) f j ] , the maximum value of P ( f ) calculated is the frequencies of the synchronization error.
Using the above steps, calculate the synchronization error of Figure 10 and obtain Figure 11.
From Figure 11, it can be seen that the three peaks in the resulting graph are around 1, 5, and 10 rad/s, which is consistent with the frequencies of the weak signal we input. This result indicates that using a synchronous fractional-order multi-wing chaotic system to detect the frequencies of weak signals is feasible.

4. Conclusions

This paper addresses the potential application of fractional-order multi-wing chaotic systems for weak signal detection. First, the impact of order in fractional calculus on the dynamic properties of chaotic systems is examined. Based on the findings, a 0.9 fractional-order multi-wing chaotic system is chosen as the weak signal detection system. In the case of amplitude detection, the various features of chaotic array detection and single detection oscillators are examined independently. Suppose that frequencies detection of weak signals is required. In that case, a fractional-order multi-wing chaotic system with drive–response synchronization is devised to detect the frequencies of the measured signal through changes in the synchronization error. When processing data, the MUSIC technique measures the synchronization error and subsequently estimates the frequencies of weak signals. Deep learning image recognition is used to process the number of wings of chaotic arrays.
The novelty of the work in this paper consists of the following: First, the creation of chaotic oscillator detection arrays highlights the advantages of fractional-order chaotic systems. Secondly, using the difference between four-winged chaotic states and double-winged chaotic states for detection makes both states resistant to noise. Finally, the states of the chaotic array are identified using deep learning images: a method that significantly reduces the judgment time compared to the usual Lyapunov method.
This paper still requires some prior knowledge to identify weak signals. Furthermore, this paper’s correctness could still be enhanced. Therefore, the main problem of chaotic weak signal detection theory remains how to detect the frequencies and amplitude of weak signals simultaneously. The primary issue is that high accuracy increases the approximate calculation error of fractional calculus and complicates the construction of chaotic arrays and synchronization. Consequently, the engineering implementation of the suggested signal detection technique will be the main focus of future research.

Author Contributions

Conceptualization, H.M. and Z.Y.; Data curation, X.W.; Formal analysis, H.M.; Funding acquisition, C.G.; Investigation, H.M.; Methodology, H.M.; Project administration, Z.Y.; Resources, Y.F.; Software, H.M.; Supervision, Y.F., C.G. and Z.Y.; Validation, H.M., Y.F. and Z.Y.; Visualization, Z.Y.; Writing—original draft, H.M.; Writing—review and editing, X.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by grants from the Key Science and Technology Research Projects in Jilin Province (grant No. 20160204020GX) and the Jilin Province Science and Technology Development Plan Project (grant No. 20190201135JC).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors on request.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

The following abbreviations are used in this manuscript:
MUSICMulti-signal classification

References

  1. Birx, D.L.; Pipenberg, S.J. Chaotic oscillators and complex mapping feed forward networks(CMFFNs) for signal detection in noisy environments. In Proceedings of the [Proceedings 1992] IJCNN International Joint Conference on Neural Networks, Baltimore, MD, USA, 7–11 June 1992. [Google Scholar]
  2. Wang, G.Y.; He, S. Quantitative study on detection and estimation of weak signals by using chaotic duffifing oscillators. IEEE Trans. Circuits Syst. I 2003, 50, 945–953. [Google Scholar] [CrossRef]
  3. Wang, Q.B.; Yang, Y.J.; Zhang, X. Weak signal detection based on Mathieu-Duffing oscillator with time-delay feedback and multiplicative noise. Chaos Solitons Fractals 2020, 137, 109832. [Google Scholar] [CrossRef]
  4. Jiao, S.B.; Jiang, W.; Lei, S.; Huang, W.C.; Zhang, Q. Research on detection method of multi-frequency weak signal based on stochastic resonance and chaos characteristics of Duffing system. Chin. J. Phys. 2020, 64, 333–347. [Google Scholar]
  5. Zhao, Z.H.; Yang, S. Application of van der Pol-Duffing oscillator in weak signal detection. Comput. Electr. Eng. 2015, 41, 1–8. [Google Scholar]
  6. Li, Q.Y.; Shi, S. Research on weak signal detection method based on duffing oscillator in narrowband noise. In Proceedings of the Artificial Intelligence for Communications and Networks: Second EAI International Conference, Virtual Event, 19–20 December 2020. [Google Scholar]
  7. Akilli, M.; Yilmaz, N.; Gediz Akdeniz, K. Automated system for weak periodic signal detection based on Duffing oscillator. IET Signal Process. 2020, 14, 710–716. [Google Scholar] [CrossRef]
  8. Li, C.S.; Qu, L.S. Applications of chaotic oscillator in machinery fault diagnosis. Mech. Syst. Signal Process. 2007, 21, 257–269. [Google Scholar] [CrossRef]
  9. Chen, H.Y.; Lv, J.T.; Zhang, S.Q.; Zhang, L.G.; Li, J. Chaos weak signal detecting algorithm and its application in the ultrasonic Doppler bloodstream speed measuring. J. Phys. Conf. Ser. 2005, 13, 320. [Google Scholar] [CrossRef]
  10. Hu, G.; Wang, K.J.; Liu, L.L. Detection Line Spectrum of Ship Radiated Noise Based on a New 3D Chaotic System. Sensors 2021, 21, 1610. [Google Scholar] [CrossRef] [PubMed]
  11. Xiong, L.; Qi, L.W.; Teng, S.F.; Wang, Q.S.; Wang, L.; Zhang, X.G. A simplest Lorenz-like chaotic circuit and its applications in secure communication and weak signal detection. EPJ-Spec. Top. 2021, 230, 1933–1944. [Google Scholar] [CrossRef]
  12. Shi, M.; Jin, C.L. Applying Improved Chaos System to Weak Signal Detection. In Proceedings of the 2015 International Industrial Informatics and Computer Engineering Conference, Xi’an, China, 10–11 January 2015. [Google Scholar]
  13. Li, G.Z.; Tan, N.L.; Su, S.Q.; Zhang, C. Unknown frequency weak signal detection based on Lorenz chaotic synchronization system. J. Vib. Shock 2019, 38, 155–161. [Google Scholar]
  14. Li, G.Z.; Zhang, B. A Novel Weak Signal Detection Method via Chaotic Synchronization Using Chua’s Circuit. IEEE Trans. Ind. Electron. 2017, 64, 2255–2265. [Google Scholar] [CrossRef]
  15. Li, Y.; Li, F.G.; Lyu, S.X.; Xu, M.; Wang, S.Y. Blind extraction of ECG signals based on similarity in the phase space. Chaos Solitons Fractals 2021, 147, 110950. [Google Scholar] [CrossRef]
  16. Li, S.Y.; Lin, Y.C.; Tam, L.M. A smart detection technology for personal ECG monitoring via chaos-based data mapping strategy. Multimed. Tools Appl. 2021, 80, 6397–6412. [Google Scholar] [CrossRef]
  17. Yin, C.; Jiang, S.B.; Zhang, Y.N. A Photoacoustic Spectrum Detection System Based on Chaos Detection of Weak Signal. In Proceedings of the 2021 IEEE 5th Advanced Information Technology, Electronic and Automation Control Conference, Chongqing, China, 12–14 March 2021. [Google Scholar]
  18. Soong, C.Y.; Huang, W.T.; Lin, F.P.; Tzeng, P.Y. Controlling chaos with weak periodic signals optimized by a genetic algorithm. Phys. Rev. E 2004, 70, 016211. [Google Scholar] [CrossRef] [PubMed]
  19. Yang, Y.; Huang, L.L.; Xiang, J.H.; Bao, H.; Li, H.Z. Generating multi-wing hidden attractors with only stable node-foci via non-autonomous approach. Phys. Scr. 2021, 96, 125220. [Google Scholar] [CrossRef]
  20. Wu, Q.J.; Hong, Q.H.; Liu, X.Y.; Wang, X.P.; Zeng, Z.G. A novel amplitude control method for constructing nested hidden multi-butterfly and multiscroll chaotic attractors. Chaos Solitons Fractals 2020, 134, 109727. [Google Scholar] [CrossRef]
  21. Sahoo, S.; Roy, B.K. Design of multi-wing chaotic systems with higher largest Lyapunov exponent. Chaos Solitons Fractals 2022, 157, 111926. [Google Scholar] [CrossRef]
  22. Peng, X.N.; Zeng, Y.C. A simple method for generating mirror symmetry composite multiscroll chaotic attractors. Int. J. Bifurc. Chaos 2020, 30, 2050220. [Google Scholar] [CrossRef]
  23. Yan, D.W.; Ji’e, M.; Wang, L.D.; Duan, S.K.; Du, X.Y. Generating novel multi-scroll chaotic attractors via fractal transformation. Nonlinear Dyn. 2022, 107, 3919–3944. [Google Scholar] [CrossRef]
  24. Hong, Q.H.; Wu, Q.J.; Wang, X.P.; Zeng, Z.G. Novel Nonlinear Function Shift Method for Generating Multiscroll Attractors Using Memristor-Based Control Circuit. IEEE Trans. Very Large Scale Integr. (VLSI) Syst. 2019, 27, 1174–1185. [Google Scholar] [CrossRef]
  25. Li, W.J.; Li, P.; Jia, M.M. Chaos control and chaos synchronization of a multi-wing chaotic system and its application in multi-frequency weak signal detection. AIP Adv. 2021, 11, 095003. [Google Scholar] [CrossRef]
  26. Yan, S.H.; Sun, X.; Wang, Q.Y.; Ren, Y.; Shi, W.L.; Wang, E. A novel double-wing chaotic system with infinite equilibria and coexisting rotating attractors: Application to weak signal detection. Phys. Scr. 2021, 96, 125216. [Google Scholar] [CrossRef]
  27. Hammouch, Z.; Yavuz, M.; Ozdemir, N. Numerical Solutions and Synchronization of a Variable-Order Fractional Chaotic System. Math. Model. Numer. Simul. Appl. 2021, 1, 11–23. [Google Scholar] [CrossRef]
  28. Chen, L.P.; Yin, H.; Huang, T.W.; Yuan, L.G.; Zheng, S.; Yin, L.S. Chaos in fractional-order discrete neural networks with application to image encryption. Neural Netw. 2020, 125, 174–184. [Google Scholar] [CrossRef] [PubMed]
  29. Baleanu, D.; Zibaei, S.; Namjoo, M.; Jajarmi, A. A nonstandard finite difference scheme for the modeling and nonidentical synchronization of a novel fractional chaotic system. Adv. Differ. Equ. 2021, 2021, 308. [Google Scholar] [CrossRef]
  30. Wang, S.J.; He, S.B.; Yousefpour, A.; Jahanshahi, H.; Repnik, R.; Perc, M. Chaos and complexity in a fractional-order financial system with time delays. Chaos Solitons Fractals 2020, 131, 109521. [Google Scholar] [CrossRef]
  31. Sayed, S.; Amir, A.; Emile Franc Doungmo Goufo. Investigation of complex behaviour of fractal fractional chaotic attractor with mittag-leffler Kernel. Chaos Solitons Fractals 2021, 152, 111332. [Google Scholar]
  32. He, Y.Z.; Fu, Y.X.; Qiao, Z.J.; Kang, Y.M. Chaotic resonance in a fractional-order oscillator system with application to mechanical fault diagnosis. Chaos Solitons Fractals 2021, 142, 110536. [Google Scholar] [CrossRef]
  33. Huang, P.F.; Chai, Y.; Chen, X.L. Multiple dynamics analysis of Lorenz-family systems and the application in signal detection. Chaos Solitons Fractals 2022, 156, 111797. [Google Scholar] [CrossRef]
  34. Li, G.H.; Xie, R.T.; Yang, H. Detection method of ship-radiated noise based on fractional-order dual coupling oscillator. Nonlinear Dyn. 2024, 112, 2091–2118. [Google Scholar] [CrossRef]
  35. Qiao, Z.J.; He, Y.B.; Liao, C.R.; Zhu, R.H. Noise-boosted weak signal detection in fractional nonlinear systems enhanced by increasing potential-well width and its application to mechanical fault diagnosis. Chaos Solitons Fractals 2023, 175, 113960. [Google Scholar] [CrossRef]
  36. Dong, K.F.; Xu, K.; Zhou, Y.Y.; Zuo, C.; Wang, L.M.; Zhang, C.; Jin, F.; Song, J.; Mo, W.; Hui, Y. A memristor-based chaotic oscillator for weak signal detection and its circuitry realization. Nonlinear Dyn. 2022, 109, 2129–2141. [Google Scholar] [CrossRef]
  37. Zhang, S.; Zeng, Y.C.; Li, Z.J. One to four-wing chaotic attractors coined from a novel 3D fractional-order chaotic system with complex dynamics. Chin. J. Phys. 2018, 56, 793–806. [Google Scholar] [CrossRef]
  38. Cai, H.; Sun, J.Y.; Gao, Z.B.; Zhang, H. A novel multi-wing chaotic system with FPGA implementation and application in image encryption. J.-Real-Time Image Process. 2022, 19, 775–790. [Google Scholar] [CrossRef]
  39. Wang, C.H. Dynamic Behavior Analysis and Robust Synchronization of a Novel Fractional-Order Chaotic System with Multiwing Attractors. J. Mathematics. 2021, 2021, 6684906. [Google Scholar] [CrossRef]
  40. Zhang, S.; Zeng, Y.C.; Li, J.Z. A Novel Four-Dimensional No-Equilibrium Hyper-Chaotic System With Grid Multiwing Hyper-Chaotic Hidden Attractors. J. Comput. Nonlinear Dynam. 2018, 13, 090908. [Google Scholar] [CrossRef]
  41. Sahoo, S.; Roy, B.K. A new multi-wing chaotic attractor with unusual variation in the number of wings. Chaos Solitons Fractals 2022, 164, 112598. [Google Scholar] [CrossRef]
  42. Tavazoei, M.S.; Haeri, M. A necessary condition for double scroll attractor existence in fractional-order systems. Phys. Lett. A 2007, 367, 102–113. [Google Scholar] [CrossRef]
  43. Nie, C.Y.; Wang, Z.W. Application of Chaos in Weak Signal Detection. In Proceedings of the 2011 Third International Conference on Measuring Technology and Mechatronics Automation, Shanghai, China, 6–7 January 2011. [Google Scholar]
  44. Faber, J.; Bozovic, D. Noise-induced chaos and signal detection by the nonisochronous Hopf oscillator. Chaos 2019, 29, 043132. [Google Scholar] [CrossRef] [PubMed]
  45. Wang, S.F. Finite-time synchronization of fractional multi-wing chaotic system. Phys. Scr. 2023, 98, 115224. [Google Scholar] [CrossRef]
  46. Kumar, V.; Heiland, J.; Benner, P. Projective lag quasi-synchronization of coupled systems with mixed delays and parameter mismatch: A unified theory. Neural Comput. Applic. 2023, 35, 23649–23665. [Google Scholar] [CrossRef]
  47. Matignon, D. Stability Results For Fractional Differential Equations With Applications To Control Processing. Comput. Eng. Syst. Appl. 1997, 2, 963–968. [Google Scholar]
  48. Hu, J.B.; Han, Y.; Zhao, L.D. A novel stablility theorem for fractional systems and its applying in synchronizing fractional chaotic system based on back-stepping approach. Acta Phys. Sin. 2009, 58, 2235–2239. [Google Scholar]
  49. Schmidt, R. Multiple emitter location and signal parameter estimation. IEEE Trans. Antennas Propag. 1986, 34, 276–280. [Google Scholar] [CrossRef]
Figure 1. If q = 0.95 , a = 4 , b = 8 , c = 1 , G = r × sin ( ω × t ) , r = 1.126 , ω = 1 , and the starting value of ( x , y , z ) = [ 0.1 ; 0.1 ; 0.1 ] , then a fractional-order multi-wing chaotic system has four wings in its attractor.
Figure 1. If q = 0.95 , a = 4 , b = 8 , c = 1 , G = r × sin ( ω × t ) , r = 1.126 , ω = 1 , and the starting value of ( x , y , z ) = [ 0.1 ; 0.1 ; 0.1 ] , then a fractional-order multi-wing chaotic system has four wings in its attractor.
Fractalfract 08 00417 g001
Figure 2. Bifurcation diagram and SE complexity diagram of the fractional-order multi-wing chaotic system Equation (2) varying with parameter q, where a = 4 , b = 8 , c = 1 , ω = 1 , G = r   × sin ( ω × t ) , r = 1.126 , and ω = 1 . The value range of q is from 0.85 to 1.15.
Figure 2. Bifurcation diagram and SE complexity diagram of the fractional-order multi-wing chaotic system Equation (2) varying with parameter q, where a = 4 , b = 8 , c = 1 , ω = 1 , G = r   × sin ( ω × t ) , r = 1.126 , and ω = 1 . The value range of q is from 0.85 to 1.15.
Fractalfract 08 00417 g002
Figure 3. Bifurcation diagram of the fractional-order multi-wing chaotic system Equation (2) varying with parameter r, where G = r   × sin ( ω × t ) , a = 4 , b = 8 , c = 1 , ω = 1 , and q = 0.9 in the left panel, q = 0.95 in the right panel, and the value range of r is from 0 to 1.4.
Figure 3. Bifurcation diagram of the fractional-order multi-wing chaotic system Equation (2) varying with parameter r, where G = r   × sin ( ω × t ) , a = 4 , b = 8 , c = 1 , ω = 1 , and q = 0.9 in the left panel, q = 0.95 in the right panel, and the value range of r is from 0 to 1.4.
Fractalfract 08 00417 g003
Figure 4. Using the G-algorithm for the numerical approximation of Riemann–Liouville, the attractor of Equation (10) is projected in the x-y plane, where the parameters are a = 4 ,   b = 8 ,   c = 1 ,   w = 1 , and the time interval is 0.01; r = 1.12939796 on the left, and r = 1.12939797 on the right.
Figure 4. Using the G-algorithm for the numerical approximation of Riemann–Liouville, the attractor of Equation (10) is projected in the x-y plane, where the parameters are a = 4 ,   b = 8 ,   c = 1 ,   w = 1 , and the time interval is 0.01; r = 1.12939796 on the left, and r = 1.12939797 on the right.
Fractalfract 08 00417 g004
Figure 5. Chaos oscillator detection array, where the value of r ranges from 1.122 to 1.137 , simulation time is 500, and calculation step size is 0.005 .
Figure 5. Chaos oscillator detection array, where the value of r ranges from 1.122 to 1.137 , simulation time is 500, and calculation step size is 0.005 .
Fractalfract 08 00417 g005
Figure 6. Chaos oscillator detection array, where the weak signal added on the left is −0.02sin ( t ) , and on the right it is 0.01 sin ( t ) .
Figure 6. Chaos oscillator detection array, where the weak signal added on the left is −0.02sin ( t ) , and on the right it is 0.01 sin ( t ) .
Fractalfract 08 00417 g006
Figure 7. The yellow solid box indicates that users prompt once on one image. This article selects the double-wing part of the four-wing chaotic attractor. Other different color boxes are automatically annotated by T-Rex2 and display other images with similar object patterns to the prompt image. Using the left half of the four-wing attractor as an example, deep learning image recognition results were able to correctly identify both biplane and four-wing attractors.
Figure 7. The yellow solid box indicates that users prompt once on one image. This article selects the double-wing part of the four-wing chaotic attractor. Other different color boxes are automatically annotated by T-Rex2 and display other images with similar object patterns to the prompt image. Using the left half of the four-wing attractor as an example, deep learning image recognition results were able to correctly identify both biplane and four-wing attractors.
Fractalfract 08 00417 g007
Figure 8. The results of chaotic attractors obtained by adding different levels of noise in deep learning image recognition.
Figure 8. The results of chaotic attractors obtained by adding different levels of noise in deep learning image recognition.
Fractalfract 08 00417 g008
Figure 9. The synchronization error of the state variables x, y, and z when the system Equations (11) and (12) are synchronized.
Figure 9. The synchronization error of the state variables x, y, and z when the system Equations (11) and (12) are synchronized.
Fractalfract 08 00417 g009
Figure 10. After adding S ( t ) of different frequencies to system Equation (16), the error between the selected state variables z in the figure show the synchronization error of the chaotic synchronization system, where S ( t ) = 0.01 sin ( ω t ) , where ω is 1, 5, and 10 rad/s.
Figure 10. After adding S ( t ) of different frequencies to system Equation (16), the error between the selected state variables z in the figure show the synchronization error of the chaotic synchronization system, where S ( t ) = 0.01 sin ( ω t ) , where ω is 1, 5, and 10 rad/s.
Fractalfract 08 00417 g010
Figure 11. The result graph of frequencies estimation using the MUSIC algorithm for the stable synchronization error shows three peaks around 1, 5, and 10 rad/s.
Figure 11. The result graph of frequencies estimation using the MUSIC algorithm for the stable synchronization error shows three peaks around 1, 5, and 10 rad/s.
Fractalfract 08 00417 g011
Table 1. Equilibrium points of system Equation (2).
Table 1. Equilibrium points of system Equation (2).
Equilibrium PointState Variable xState Variable yState Variable z
S 0 000
S 1 1 2 G 1 2 G 2 1 2
S 2 1 2 G 1 2 G 2 1 2
S 3 1 + 2 G 1 + 2 G 2 1 2
S 4 1 + 2 G 1 + 2 G 2 1 2
Table 2. Stability of equilibrium points of system Equation (2).
Table 2. Stability of equilibrium points of system Equation (2).
Value of Driving Signal GEigenvalues at Equilibrium Point S 0 Eigenvalues at Equilibrium Point S 1   S 2 Eigenvalues at Equilibrium Point S 3   S 4
0−4, −1, 1 Unstable Saddle Point−4.7186, 0.3593 ± 1.8060i Unstable Saddle Point−4.7186, 0.3593 ± 1.8060i Unstable Saddle Point
1−1, −1.5 ± 1.3229i stable Focus0.8128, −2.4064 ± 1.5374i Unstable Saddle Point−5.8526, 0.9263 ± 2.3962i Unstable Saddle Point
5−1, −1.5 ± 13.9194i stable Focus2.5574, −3.2787 ± 5.2185i Unstable Saddle Point−8.3744, 2.1872 ± 3.2614i Unstable Saddle Point
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Mao, H.; Feng, Y.; Wang, X.; Gao, C.; Yao, Z. Application of Fractional-Order Multi-Wing Chaotic System to Weak Signal Detection. Fractal Fract. 2024, 8, 417. https://doi.org/10.3390/fractalfract8070417

AMA Style

Mao H, Feng Y, Wang X, Gao C, Yao Z. Application of Fractional-Order Multi-Wing Chaotic System to Weak Signal Detection. Fractal and Fractional. 2024; 8(7):417. https://doi.org/10.3390/fractalfract8070417

Chicago/Turabian Style

Mao, Hongcun, Yuling Feng, Xiaoqian Wang, Chao Gao, and Zhihai Yao. 2024. "Application of Fractional-Order Multi-Wing Chaotic System to Weak Signal Detection" Fractal and Fractional 8, no. 7: 417. https://doi.org/10.3390/fractalfract8070417

APA Style

Mao, H., Feng, Y., Wang, X., Gao, C., & Yao, Z. (2024). Application of Fractional-Order Multi-Wing Chaotic System to Weak Signal Detection. Fractal and Fractional, 8(7), 417. https://doi.org/10.3390/fractalfract8070417

Article Metrics

Back to TopTop