Chaotic Signal Denoising Based on Adaptive Smoothing Multiscale Morphological Filtering

Nonlinear time series denoising is the prerequisite for extracting effective information from observation sequence. An effective chaotic signal denoising method not only has a good signal-to-noise ratio (SNR) enhancement performance, but also can remain as a good unpredictable denoised signal. However, the inherent characteristics of chaos, such as extreme sensitivity to initial values and broadband spectrum, pose challenges for noise reduction of polluted chaotic signals. To address these issues, an adaptive smoothing multiscale morphological filtering (ASMMF) is proposed to reconstruct chaotic signals. In the process of noise reduction for contaminated chaotic signals, firstly, a multiscale morphological filter is constructed adaptively according to the multiscale permutation entropy (MPE) of morphological filter residuals, and the contaminated signals are filtered. Secondly, the weight coefficients of scale structural element are calculated by the residual sum of squares operation, and the chaotic signals are reconstructed. Finally, the resampled filter signals are smoothed by cubic B-spline interpolation operation. In the experiment, the Lorenz signal with white Gaussian noise, the measured sunspot, and the chaotic vibration signal are reconstructed by four comparison methods. 0e test results show that the proposed ASMMF method has obvious advantages in noise suppression and topological trajectory restoration.


Introduction
Chaos is a seemingly random irregular motion that occurs in a deterministic system [1]. It is vital to preprocess observation data to effectively extract desired chaotic information. In the process of analyzing the existing chaotic phenomena, the measured chaotic behavior is susceptible to noise interference. erefore, it is desirable to preprocess the noise of observation data without distorting the dynamic characteristics of underlying signals [2,3]. Owing to chaotic signals usually having a broadband spectrum that overlaps with the spectrum of noise [4], the conventional linear low-pass filtering method is ineffective for chaotic signal denoising. If these methods are applied in noise reduction of chaos signals, the clean chaotic signals may be oversmoothed severely [5]. erefore, it is of great significance to study the corresponding noise reduction methods according to the characteristics of chaotic signals.
Over the past decades, a number of approaches for chaotic signal noise reduction have been proposed [3,[6][7][8][9][10][11][12][13]. e singular spectrum analysis method [6] and local projection method [7,8] are denoising algorithms based on phase space reconstruction; due to their different dynamic characteristics after reconstructing, the chaotic signals are separated from noise. But it is difficult to accurately determine the noise boundary points of noise reduction by singular spectrum analysis, and the parameter selection method of local projection has limitations. erefore, the noise reduction ability will decrease if the noises of denoising chaotic signals are severe. e threshold denoising method based on time-frequency analysis, such as the wavelet threshold (WT) denoising method [3], and empirical mode decomposition (EMD) [9,10] are noise suppression algorithms which decompose the noisy signal in time-frequency domain and separate the signal from the noise by threshold processing. However, the denoising performance of the WT method is greatly impacted by the wavelet basis and the number of layers, which reduces the adaptability of the method. Although the EMD threshold denoising overcomes the issue that WT must select the appropriate wavelet basis according to the signal characteristics, it is still hard to determine the threshold properly. e local curve fitting method [11], which smooths the noisy signal by segments, uses the least square polynomial fitting to realize the signal denoising, and it is arduous to reconstruct precisely by local linear approximation due to the high nonlinearity of chaos. Cooperative filtering denoising [12] uses the self-similarity of chaotic signals in time domain to group similar blocks, and the noise suppression effect of the reconstructed signal after filtering is obvious. Nevertheless, for the reason that the selection of filtering parameters of this method is not adaptive, the application of this method is still limited. Wang et al. [13] further improved the cooperative filtering denoising method; the complexity of chaotic signals is analyzed and the optimal filtering parameters are adaptively selected according to the smallest permutation entropy (PE) [14]. Most of the existing noise suppression methods are proposed on the basis of low-pass filter [15,16]. Analyzing the above research results, the study trend of chaotic signal noise suppression is developing from global noise reduction to local noise reduction, and the noise reduction methods proposed for dynamic characteristics of chaotic signal have better effect.
Morphological filtering (MF) is a nonlinear filtering technique [17]. It can be employed to make translation matching and local correction of the original signal from front to back by constructing a certain structuring element (SE), and the morphological characteristics of the original signal are preserved while the noises are suppressed [18]. Currently, morphological filtering has been widely used in image processing [19,20], signal analysis [21], feature extraction [22], and other fields. In this paper, the MF is first introduced to reduce the noise of contaminated chaotic signal. Combined with the self-similar characteristics of chaotic signals, an adaptive smooth multiscale morphological filtering (ASMMF) method for chaotic signal denoising is proposed on the basis of existing algorithms. Within this scheme, the scale range of structural element is determined in accordance with the self-similarity of temporal waveform of chaotic signal. e multiscale permutation entropy (MPE) [23] of the residual signal is calculated to select the optimal combination of morphological filtering scales, and the residual variance is calculated to determine the weight coefficients of different scale filtering results. e chaotic signal is reconstructed by weighted optimal scale combination of morphological filtering results. To solve the problem of topping distortion in morphological filtering, the reconstructed signal is resampled and smoothed by cubic B-spline interpolation. In order to distinguish ASMMF from the adaptive multiscale morphological filtering without smoothing, the latter is called multiscale morphological filtering for short. e ASMMF is an improved algorithm on the basis of AMMF. e compendium of this paper is organized as follows. In Section 2, the process of ASMMF is described in detail. Section 3 analyses the parameters of the algorithm and proposes the corresponding optimization scheme. Moreover, a noisy Lorenz signal is taken to describe the optimizing process. Section 4 applies the proposed method and three comparison methods to denoise different noisy chaotic signals. e ultimate conclusions are drawn in Section 5.

Adaptive Smoothing Multiscale Morphological Filtering
In order to denoise the contaminated chaotic signals, a novel method called ASMMF is presented in this paper. And this algorithm is divided into three steps: multiscale morphological filtering, signal reconstruction, and signal smoothing. Figure 1 shows the noise reduction process.

Multiscale Morphological
Filtering. Set f(n) as a onedimensional (1D) signal, which is a discrete function over the domain F � (0, . . . , N − 1), set g(m) as the SE, which is a function over the domain G � (0, . . . , M − 1), and the stipulated value of N should be larger than that of M. e dilation operation and the erosion operation are two kinds of basic operations of mathematical morphology (MM) [24], which can be denoted as symbols ⊕ and Θ, respectively, and the signal process can be represented as follows [17]: (1) A number of basic operators are used in MF. e opening operation and closing operation of a set f(n) by SE g(m) are the most widely used; they can be defined as where symbols ∘ and • stand for the opening operation and closing operation, respectively. e opening operation flattens the positive impulses and matches the negative ones, and the closing operation can flatten the negative impulses and matches the positive ones [20]. On the basis of the sequential operation of opening and closing, the openclosing operator and the close-opening operator are defined as [25,26] e output signals of the open-closing operator and the close-opening operator have the problem of median bias; thus, the combined morphological filter (CMF) is used in signal processing. e combined morphological filter is able to smooth the signal by eliminating the positive and the negative impacts in the signal simultaneously. e combination morphological filter is further defined as CMF(f(n)) � (FOC(f(n)) + FCO(f(n))) 2 .
For complex noisy signals, morphological filtering of single-scale structural element has limited accuracy in 2 Complexity reconstructing. It is necessary to construct multiscale morphological filters to denoise target signals at different levels and scales. Let λg be the SE at scale λ, which is obtained by λ − 1 times dilation operation; λg is described as [27] λg � g ⊕ g ⊕ · · · ⊕ g √√√√√ √√√√√ e multiscale CMF can be redefined as CMF noise reduction results can effectively restore the original signal waveform. Also, the weighted average operation can restrain the statistical bias of denoising results. Indeed, compared with a single-scale CMF, the multiscale CMF is more effective in noise reduction.

Signal Reconstruction.
ese results of morphological filtering at different scales can be obtained by multiscale morphological analysis of noisy signal. Each filtering result at single scale reflects a specific component contained in the signal. e different components of chaotic signals should appear in multiple scales as a processed signal. e reconstructed signal is obtained by aggregating the filtering results of different scales as follows: where y j is a result of morphological filtering on scale j, ω j is the weight coefficient of y j , and Y is the reconstruction result. e weight coefficient directly affects the effect of multiscale morphological filtering. For the purpose of Complexity 3 obtaining a better denoising result, the square operation of residual amplitude difference (RAD) before and after filtering processing is used to determine the weight coefficients adaptively. e RAD can be defined as e weight coefficient can be expressed as

Signal Smoothing.
It is found that the waveform of the aggregated signal is not smooth enough, and the local waveform is slightly rough. Owing to the morphological filter with large-scale SE blurring the detail information of signal and bringing the top-weakening distortion in the process of reducing noise, the signal directly gathered by the results of morphological filtering at different scales is not smooth enough. However, the occurrence of the topweakening distortion is regular all the way. e distortion only made some local waveforms a little rough, but did not change the general shape of the whole waveform. In order to reduce the top-weakening distortion, a new method of resampling the filtered signal is proposed according to the regularity of the distortion. e cubic B-spline interpolation method is proposed to smooth the filtered signal after resampling.
Let d i (i � 0, 1, . . . , n) be the de Boor point, where i is the serial number and n is the number of the de Boor point. e k-th B-spline function is defined as follows [28,29]: where t 0 , t 1 , . . . , t n+k+1 are node vectors and B i,k (t) is a base function of the node vector. When the denominator in equation (10) is 0, the value is specified as 0. When k is 3, the cubic uniform B-spline function defined in [0, 1] is described as Equation (11) is a cubic polynomial divided into uniform segments; the signal processed by cubic B-spline interpolation has the characteristics of detail smoothness and linear independence.

Algorithmic Parameter Analysis and Adaptive Optimization
In this section, the Lorenz signal with superimposed Gaussian white noise is taken as an example to discuss the influence of algorithm parameters on noise reduction performance. e equation of Lorenz system is expressed as [ where a � 10, b � 8/3, and c � 28; the system is in chaos. e fourth-order Runge-Kutta method is used to solve the equation, and the sampling time is set to 0.01. Chaotic signals are generated by the state variable x. Figure 2(a) shows the waveform of state variables changing with time. e chaotic signal is superimposed with white Gaussian noise of different intensities. Figure 2(b) shows a Lorenz signal with a superimposed signal to noise ratio (SNR) of 15 dB noise. e SNR is defined as [31] SNR � 10 log 10 where s(t) is the original chaotic time series. When calculating the signal-to-noise ratio of the input signal (SNR in ), the x(t) is regarded as a noisy chaotic signal. Nevertheless, if the calculated object is the output signal to noise ratio (SNR out ), the x(t) is regarded as the denoised signal.
e ASMMF method makes full use of the self-similar structure characteristics of chaotic signals and has good chaotic attractor reduction performance. Just like all the MF based algorithms, ASMMF algorithm also needs to specify the shape and select the optimal scale combination of SE.

e Shape Selection of Structural Element.
e SE can be used as a filter window in the process of MF, and its geometric feature matches the signal much better, and then main waveform of this part of the signal can be retained. e shape and size of structural element play an important role in the application of mathematical morphology. e triangular SE, the sinusoidal SE, and the linear SE are widely used SEs for signal analysis [32][33][34][35][36].
Morphological operation is a numerical operation of adding or subtracting extreme values. Flat SE is defined as the linear SE with height of 0, which can retain the shape feature of the impulse as much as possible [34]. And it has the simplest algorithm and higher operation efficiency [35]. However, if the raw signal is processed by a scale far away from the theoretical central scale, the filtered result is often heavily polluted. Fortunately, multiscale morphological filtering with adaptive weighting can solve this problem [36]. erefore, the flat SE is therefore chosen in filtering chaotic signal.

e Selection of Optimal Scale Combination.
Set the initial SE to {0, 0, 0}. Scale of flat SE is used to represent the length of SE, and define the scale of initial SE to be λ � 1. e difference between flat structural element at different scales is the number of zeros. e selection of scale coefficients of structural element of morphological filters has a significant effect on the denoising results. If the scale of structural element is relatively small, as shown in Figure 3(a), the effect of morphological filter on noise suppression is not obvious. On the contrary, when the scale of structural element is too large, the effect of noise reduction is obvious as shown in Figure 3(b), but the waveform has serious distortion. e experimental analysis shows that the optimal scale combination is affected by the signal characteristics and noise level, and structural element directly affects the effect of morphological filtering. In order to achieve the best filtering result, it is necessary to select the appropriate scale combination of structural element. is optimal filtering parameter is constantly changing in practical application. How to realize the automatic optimization of filter parameters determines the adaptability of the algorithm. In order to increase the adaptability of the algorithm, an adaptive method for determining the optimal scale combination parameter is proposed in this paper. e flowchart to determine the optimal size combination is shown in Figure 4. Figure 4, the adaptive selection process of optimal scale combination is divided into two steps.

As shown in
Step 1.
e selection of the maximum scale structural element. An appropriate maximum scale SE is able to reduce the computational complexity and operating time. e maximum scale of structural element is denoted as λ max . e experiment shows that the optimal length of the SE is closely related to the sampling frequency of the signal and the number of sampling points in the fault period. λ max should be smaller than the sampling point in the fault period, which is corresponding to a period of Lorenz signal. e number of sampling points corresponding to a period is denoted by f s /f 0 , where f s is sampling frequency (f s � 100 Hz) and f 0 is the characteristic frequency of signal. According to the calculation method of Reference [37], the characteristic frequency is defined as  Complexity 5 where N(T) is the number of turns performed in T.
where symbol ⌊•⌋ represents round down the number. e characteristic frequency of Lorenz signal is f ≈ 7.04 Hz; therefore, f s /f ≈ 1.40 Hz Hz. λ max should be less than 15, so set λ max to 14 adaptively.
e multiscale permutation entropy of residual signal is used to select the optimal scale combination of SE adaptively. MPE is an average entropy parameter developed on the basis of PE [14] to describe the irregularity degree and self-similarity of time series at different scales. e calculation process of MPE is described as follows.
Set x(t), t � 1, 2, . . . , N { } to be a time series, and the scale τ represents the window length of coarsening. In the first place, the time series is coarsely granulated under the condition of scale factor τ. And the resampling of scale τ is carried out; the sequence after coarse granulation is denoted as where j � 1, 2, . . . , ⌊N/τ⌋, and ⌊•⌋ denotes downward rounding.
en, the phase space of time series z (τ) (j), j � 1, 2, . . . , N is reconstructed as where m indicates the embedded dimension and ε indicates the delay time. e m data points of z (τ) (j) are rearranged in ascending order as where j d , d � 1, 2, . . . , m is the column index of each element in the reconstructed component. In case there are m! symbol sequences at most, k ≤ m!.
e frequency of occurrence of each permutation is recorded as f (τ) k , k � 1, 2, . . . k, k f k � 1; afterwards, the probability of each sequence in ε scale is described as According to Reference [14,15], the MPE is ultimately defined as In order to facilitate the analysis, the MPE is normalized by the maximum value of ln m!, that is, e MPE vector can be described as P � [P (1) , P (2) , . . .]. e MPE consists of four parameters: embedding dimension m, time delay ε, length of time series N, and scale factor τ. e PE value highly depends on the selected embedding dimension m and the scale factor τ. e PE method works with the embedding dimension 3 ≤ m ≤ 7 [38], and the embedding dimension m and time series N should satisfy the condition 5m! ≤ N [39]; then, the embedding dimension is set to 4 in this paper. According to Reference [40], the delay time is set to 1. e selection criterion of scale factor is given as τ ≤ ⌊N/(m + 1)!⌋ in Reference [41], and the maximum value of scale factor is 10 in this paper. Figure 5 shows the MPE analysis of morphological filtering results at different scales.
According to the analysis, with the increase of structural element scale, we find that the MPE values decrease gradually. And the complexity of the sequence decreases with the increase of the scale. e increase in scale makes the filtered sequence more orderly; however, it does not mean that the signal reduction effect is better. As shown in Figure 3(b), the filtering results with too large scale of structural element will also have smaller PE values.
In this paper, the optimal filtering scale is determined by the residual multiscale permutation entropy of the filtering The maximum scale selection Step 1: Step 2: 6 Complexity result. If the residual MPE is the largest, it means that the randomness of residual is the strongest and the filtering effect is the best. e MPE of multiscale morphological filtering residuals is shown in Figure 6.
As shown, PE of the residuals corresponding to different τ changes with the increase of λ, and all go through the process of gradually increasing from small to a maximum and then gradually decreasing. λ corresponding to the maximum of τ scale can be used as the optimal filtering scale coefficient of chaotic signal at this scale λ bestτ . Λ � λ best1 , λ best2 , . . . , λ bestτ is initially determined to be the optimal combination of scale coefficients.
ose repeat scales should be removed to ensure that each scale appears only one time. Finally, the optimal combination of scale coefficients for noisy signals as shown in Figure 2(b) turns out to be Λ � 1, 4,6,7,8,9,13 { }. After determining the optimal combination of filtering scales, the weighting coefficients of filtering results are calculated according to equation (9), and the chaotic signal is aggregated according to equation (7). e aggregation result is the noise reduction result by AMMF, which is shown in Figure 7(a). Compared with the noisy chaotic signal and the original signal, the waveform of the aggregated signal is similar to the original signal as a whole, but the local details are slightly rough. However, after smoothing, as shown in Figure 7(b), the reconstructed signal waveform is basically the same as the original signal.

Experimental Analysis
In order to verify the effectiveness of proposed ASMMF for different chaotic signals, the existing EMD-CIT [9], WT [3], and AMMF methods are applied for comparison, and the noise suppression results of the Lorenz signals with different intensities of Gaussian white noises, the sunspot chaotic signals, and chaotic vibration signals are, respectively, shown in this section.

Lorenz Chaotic Signal Denoising.
e total length of the original chaotic Lorenz signal is set at 8000. SNR out and the root mean square error (RMSE) [42] indexes are both used to evaluate the noise reduction performance. RMSE can be defined as follows: Figure 8 shows the denoising results of different noisy signals processed by four different algorithms. In the results, when the signal-to-noise ratio of input signal is low, SNR out of the unsmoothed adaptive multiscale morphological filtering is slightly lower than that of the EMD threshold method and the wavelet threshold method. As the noise intensity increases gradually, its performance index gradually exceeds the wavelet threshold method, but the overall denoising performance is poor. e morphological filtering algorithm after smoothing is superior to the other three methods in both the aspects.
In order to intuitively reflect the denoising results of different algorithms, the phase diagrams based on the chaotic characteristic are drawn ulteriorly. Figure 9 shows the denoising output results by four methods for noisy chaotic signal with input signal-to-noise ratio of 10 dB. e two-dimensional phase diagrams of the clean signal and the noisy signal are plotted in Figures 9(a) and 9(b). e noise reduction results of WT and EMD-CIIT are shown in Figures 9(c) and 9(d). According the phase diagrams, both of them can restore the attractor shape of chaotic signals smoothly, but the self-similarity structure of the system cannot be clearly expressed due to the orbital deformation. Unsmoothed morphological filtering result is shown in Figure 9(e) e shape of the attractor restored is not smooth enough, and the effect is not good enough. Furthermore, the cancellation distortion of local waveform in the process of noise reduction is obvious. e ASMMF is an improvement Complexity algorithm of AMMF. After smoothing the rough waveform of morphological filtering, the shape and self-similar structure of chaotic attractors can be restored more clearly.   8 Complexity changes of the Earth's magnetic field and climate. Recently, chaos theory has become an important means to study the rule of sunspot activity. e observed sunspot signal is always chaotic and noisy [10,31], so it is necessary to effectively reduce the noise of the observed data. e sunspot sequence from January 1749 to September 2016 provided by NASA is selected as experimental data (http://solarscience. msfc.nasa.gov/green-wch.shtml). e Lyapunov exponent (LE) of the attractor represents the predictability and the sensitivity of dynamical system. erefore, the LE is an important feature for the study of chaotic behavior. For a chaotic system, the largest Lyapunov exponent (LLE) is positive [43]. Referring to the calculation method in Reference [36], the LLE is 0.7615.

Sunspot
To investigate the method for different kinds of noise, EMD-CIT, WT, AMMF, and ASMMF are used to denoise the sunspot data. Compared with the original sunspot observed data shown in Figure 10 noise reduction results is not smooth enough. And the waveform of ASMMF after noise reduction is the smoothest. e projection of the attractor reconstructed from the data sequence on the two-dimensional plane before and after noise reduction is shown in Figure 11. It is clear that the chaotic attractor phase diagram of the measured lunar sunspot signal is very disordered due to the noise interference, and the attractor trajectory in the phase diagram space is not smooth. No regular geometric structure is shown in Figure 11(a). Noise reduction results of four mentioned methods are shown in Figures 11(b)-11(e); all of them obviously reduce the noise and obtain the geometric structure of attractor with strong regularity. Moreover, the attractor orbits obtained by the ASMMF method are smoother and clearer than those of the other three methods. It shows that this method has certain advantages in reducing noise and enhancing the inherent deterministic components of chaotic sequences.
In order to intuitively verify the advantages and disadvantages of noise reduction methods, MPE is used to describe the complexity and regularity of time series on multiple time scales. Figure 12 shows the analysis results by MPE.
By direct comparison, the four methods all perform obvious effect in chaotic signal denoising. e result of noise reduction by these methods is that the MPE value decreases compared with the signal before noise reduction, and the signal after noise reduction is more regular. If the time scale is small, the PE of the unsmoothed adaptive multiscale morphological filtering is lower than that of the EMD threshold method and the wavelet threshold method. As the time scale increases gradually, its performance index gradually exceeds the wavelet threshold method, but the overall denoising performance is poor. e denoising result using ASMMF algorithm has the lowest complexity in different time scales. It indicates that the intrinsic deterministic components of chaotic sequences are enhanced after noise reduction. e simulation results show that compared with WT, EMD-CIIT, and AMMF, ASMMF has a better denoising effect on the observed sunspot sequence, with a smoother phase diagram, and higher self similarity.

Chaotic Vibration Signal
Denoising. Detection of faults of rotating machinery by the complex noisy vibration signals is very important and difficult. In recent years, features of chaotic behavior are very good for fault diagnosis and identification. But the features of the LLE, approximate entropy, and correlation dimension are easily affected by noise. Clear attractor orbits and signals that do not damage nonlinear dynamical structures are helpful to obtain chaotic behavior characteristics accurately. So, detecting chaotic vibration signals in noise environment is a significant task.
In the experiment, the Bently RK-4 rotor system test bed was used to simulate the fault; it is equipped with a signal preadapter, speed control device, and bearing and oil pump system for oil film instability, which can realize the oil film whirl and oil film oscillation test. Eddy current sensors were installed on both sides of the single-disk rotor to measure radial vibration displacement of rotating shaft. ZonicBook/ 618E produced by Iotech Company in the United States collects vibration signals of rotors during rotation with sampling frequency of 1280 Hz. When the rotation speed is maintained at 4560 r/min, the rotor system shows obvious resonance oil film oscillation phenomenon, and the timedomain waveform is shown in Figure 13.
As mentioned in Section 4.2, the LE is an important feature for the study of chaotic behavior. For a chaotic system, the LLE is positive. Figure 14(a) shows the largest Lyapunov exponents, and phase diagram for the vibration signal is positive. So, the collected vibration signal is chaotic. Figure 14(b) shows the phase diagram for the vibration signal. Due to the influence of noise, the chaotic attractors cannot be clearly displayed.
Due to the influence of noise, the chaotic attractors cannot be clearly displayed. It is necessary to denoise the collected vibration signal. In the experiment, the contaminated chaotic vibration signal is denoised by the 4 noise reduction methods mentioned. Figure 15 shows the denoising results. e phase diagram of the ASMMF method after noise reduction is smoother, clearer, and more regular than other methods. As shown in Figure 16, as the time scale increases  Complexity gradually, with the increase of time scale, the PE values of AMMF, WT, and EMD CIIT are not significantly lower than before noise reduction. e PE value of noise reduction result by WT algorithm on some time scales is even greater than that before noise reduction, and the overall denoising performance is poor. e MPE of the ASMMF method's denoising result is the minimum, which represents that the result of denoising by this algorithm is the most regular. e   experimental results show that this method can also achieve good noise reduction effect on the vibration signal, and the denoising performance is better than the proposed three other methods.

Conclusions
We have proposed an adaptive smoothing morphological filtering method for denoising contaminated chaotic signals. is algorithm solves the problem of top-weakening distortion in morphological filtering and can filter the signal in multilevel and multiscale adaptively. In this paper, AMMF, EMD-CIIT, and WT are compared with ASMMF. e polluted Lorenz signals, the sunspot chaotic signals, and chaotic vibration signals are denoised by 4 methods, respectively. e ASMMF is superior to the other three methods in SNR out , RMSE, and MPE for the noise suppression results of contaminated chaotic signals. Moreover, when the noise of signal is strong, the attractor trajectory of the EMD-CIIT and the WT denoising results have no regular geometric structure, but the result obtained by the proposed method has smoother and more regular attractor orbits.
e comparison results show that the topology structure of the chaotic attractor restored by ASMMF is closer to the real signal. So, the proposed method has a certain reference value for actual engineering application. In the future, denoising of different kinds of chaotic signals will still be a hot topic. We will apply the proposed algorithm to the noise reduction of hyperchaotic signals and fractional chaotic signals.

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

Conflicts of Interest
e authors declare no conflicts of interest.