A Denoising Method of Ship Radiated Noise Signal Based on Modiﬁed CEEMDAN, Dispersion Entropy, and Interval Thresholding

: Due to the non-linear and non-stationary characteristics of ship radiated noise (SR-N) signal, the traditional linear and frequency-domain denoising methods cannot be used for such signals. In this paper, an SR-N signal denoising method based on modiﬁed complete ensemble empirical mode decomposition (EMD) with adaptive noise (CEEMDAN), dispersion entropy (DE), and interval thresholding is proposed. The proposed denoising method has the following advantages: (1) as an improved version of CEEMDAN, modiﬁed CEEMDAN (MCEEMDAN) combines the advantages of EMD and CEEMDAN, and it is more reliable than CEEMDAN and has less consuming time; (2) as a fast complexity measurement technology, DE can e ﬀ ectively identify the type of intrinsic mode function (IMF); and (3) interval thresholding is used for SR-N signal denoising, which avoids loss of amplitude information compared with traditional denoising methods. Firstly, the original signal is decomposed into a series of IMFs using MCEEMDAN. According to the DE value of IMF, the modes are divided into three types: noise IMF, noise-dominated IMF and pure IMF. After noise IMFs are removed, the noise-dominated IMFs are denoised using interval thresholding. Finally, the pure IMF and the processed noise-dominated IMFs are reconstructed to obtain the ﬁnal denoised signal. The denoising experiments with the Chen’s chaotic system show that the proposed method has a higher signal-to-noise ratio (SNR) than the other three methods. Applying the proposed method to denoise the real SR-N signal, the topological structure of chaotic attractor can be recovered clearly. It is proved that the proposed method can e ﬀ ectively suppress the high-frequency noise of SR-N signal.


Introduction
Ship radiated noise (SR-N) signals are caused by various vibrations and sound sources from a ship, including mechanical noise, propeller noise, and hydrodynamic noise [1]. SR-N signals contain many types of information about ship characteristics, such as the type of ship, speed, tonnage, etc. [2,3]. The navy collects SR-N signals using underwater acoustic detection equipment, torpedoes, etc., and uses the SR-N signal as a source of target information [4]. At the same time, SR-N signals are the basis for the development, testing and use of these naval equipment. Due to the marine environmental noise and the time-varying characteristics of the underwater acoustic channel, the measured SR-N signal is a nonlinear and non-stationary chaotic signal [5,6]. Noise reduction of SR-N signal is beneficial to further feature extraction, classification and detection [7]. However, traditional linear and frequency-domain denoising methods cannot be directly applied to SR-N signal. It is necessary to find a denoising method suitable for SR-N signal. At present, a data-driven adaptive signal decomposition method provides a new idea for underwater acoustic signal processing.

MCEEMDAN Method
In the decomposition process, NPE is first used to determine whether the signal contains abnormal components, and then to decide whether to add noise for assistance. Therefore, a NPE-based MCEEMDAN algorithm is proposed in this paper.

NPE
The algorithm steps of NPE are as follows [21]: For a time series x(n) = x(1), x(2), . . . , x(N) , with given embedding dimension m and time delay d, the embedding vectors are represented as: where j = 1, 2, . . . , N − (m − 1)d. Then each of the N − (m − 1)d vectors is arranged in ascending order: x( j + (k 1 − 1)d) ≤ x( j + (k 2 − 1)d) ≤ · · · ≤ x( j + (k m − 1)d) (2) If x( j + (k l1 − 1)d) = x( j + (k l2 − 1)d), then x( j + (k l1 − 1)d) ≤ x( j + (k l2 − 1)d), (l1 ≤ l2). Next, an ordinal pattern, which is one of m! possibilities, can be assigned to each X( j) and be denoted as: where 1 ≤ l ≤ m!. Compute the probability distribution p l of each ordinal pattern, then the NPE of time series x can be defined as: where h ≤ m!. Compare with PE, −p l ln(p l ) is replaced by p 2 l , and a constant − 1 m! is added. When p l = 1 m! , NPE reaches the smallest value 0, this means that the distance to white noise is 0. Li et al. [21] have proved that NPE is more stable than PE, and NPE has strong noise recognition ability.
a. Yes, go to step 5, and save 1st to k−1th IMFs as noise-assisted IMFs; b. No, do k = k + 1, go back to step 3 for next IMF; Step 5. Decompose r k−1 by EMD to obtain the remaining IMFs.
Finally, the combination of noise-assisted IMFs and the remaining IMFs decomposed by EMD is the decomposition result of MCEEMDAN. Compared with CEEMDAN, MCEEMDAN further reduces the amount of calculation. Through multiple tests, we set the threshold of NPE to 0.18.

DE
Given a univariate signal x = {x 1 , x 2 , . . . , x N }, the DE algorithm is as follows [30]: Step 1. Mapping of the N elements in the input signal x into c classes. The first part of the mapping applies the normal cumulative distribution function (NCDF) to x in order to obtain a normalized vector y (with values from 0 to 1). Then, y is linearly mapped to a vector z with integers values from 1 to c using the formula round(c · y + 0.5), where the rounding can be an increasing or a decreasing to the next digit; Step 2. Creation of n_emb = N − (m − 1)d embedding vectors (namely z_emb), where m is the number of embedding dimensions and d the time delay. Each element of the embedding vectors is defined choosing m elements from z at distance d starting from position i (that is, z_emb ij = z(i + ( j − 1)d), with i and j ranging from 1 to n emb and from 1 to m respectively). Each embedding vector is mapped to a dispersion pattern Step 3. Calculation of the probabilities p i of each of the c m dispersion patterns (i.e., the number of occurrences of each dispersion pattern out of the number of embedding vectors); Step 4. Calculation of the normalized dispersion entropy There are three parameters for DE, including the embedding dimension m, the time delay d, and number of classes c. This paper uses the same parameters as the literature [25], that is, m = 2, d = 1, and c = 6.

Interval Thresholding
In the EMD-based SR-N signal denoising method, noise evaluation criteria are often used to divide the modes into noise IMFs and pure IMFs, and then the pure IMFs are reconstructed to obtain denoised signal. However, there may be useful information in the noise IMF, so it is necessary to reduce the loss of useful information. In this paper, the EMD-IT method is introduced into the SR-N signal denoising, and the interval thresholding is used to denoise the noise-dominated IMFs. The main principle of EMD-IT is as follows [26,31]: the first step is to find the interval i m = [z i m , z i m+1 ] of each IMF, where z i m and z i m+1 are the mth and m+1th zeros of the ith IMF, and there is only one extreme value h m (r i m ) in this interval. The second step is to compare the magnitude of the extreme value and the threshold. If the extreme value is greater than the threshold, all amplitudes within the interval are retained, otherwise all amplitudes within the interval are set to zero. The interval thresholding can be expressed as: where m = 1, 2, . . . , N i , and N i is the total number of extreme points of the ith IMF. T i is the threshold value of the ith IMF, which can be calculated by: where N is the number of sampling points of the original signal. C represents the weight, generally taking 0.6 to 0.8 [32]. E i represents the noise energy in the ith IMF, which can be estimated by: where β and ρ are parameters that have be estimated by a large number of independent noise realizations and their IMFs. Flandrin et al. [33] proposed setting the values of these two parameters to 0.719 and 2.01, respectively. E 1 is the noise energy of the first IMF, which can be estimated by the noise standard deviation of IMF 1 , that is:

The Proposed Method
We propose a denoising method named MCEEMDAN-DE-IT for ship radiated noise signal. The flow chart of the proposed denoising method is shown in Figure 1.
Step 1. The SR-N signal S(t) is decomposed into a series of IMFs by MCEEMDAN.
where IMF i (t) is the ith IMF and K is the number of IMFs (the residue from the decomposition is treated as the last IMF); Step 2. Calculate the DE of each IMF; Step 3. According to the DE value, IMFs are divided into three types: noise IMFs, noise-dominated IMFs, and pure IMFs. If DE is greater than 0.88, the corresponding IMF is determined to be a noise IMF. Then look for the last IMF with DE greater than 0.65, and determine the IMF and the previous IMFs as noise-dominated IMFs. In addition, the remaining IMFs are pure IMFs; Step 4. Interval threshold denoising is performed on all noise-dominated IMFs; Step 5. Finally, the noise IMFs are removed, and the pure IMFs and the processed noise-dominated IMFs are reconstructed to obtain the final denoised signal. Therefore, the denoised signal S(t) can be expressed as: where L 1 and L 2 represent the order of the noise-dominated IMFs. I MF i (t) is the thresholded version of the noise-dominated IMF i (t).

Evaluation Criteria for SR-N Signal Noise Reduction
To quantitatively evaluate the denoising effect of the method on the SR-N signal, it is necessary to calculate some feature parameters before and after noise reduction. These feature parameters include noise intensity (NI) [5], spatial-dependence recurrence sample entropy (sdrSampEn) [34] and DE.

SdrSampEn
The traditional sample entropy (SampEn) is a tool to measure the complexity of time series and is widely used in chaotic signal analysis. However, the SampEn does not take sequential information, which is inherently useful, into its calculation of sample similarity. So, Pham et al. proposed a tool that is based on the mathematical principle of the SampEn and enables the capture of sequential information of a time series in the context of spatial dependence provided by the binary-level co-occurrence matrix (BLCM) of a recurrence plot (RP). The sdrSampEn algorithm can be described as follows [33]: Step 1. Given a time series and model parameters embedding dimension m, and spatial offset δ;

Evaluation Criteria for SR-N Signal Noise Reduction
To quantitatively evaluate the denoising effect of the method on the SR-N signal, it is necessary to calculate some feature parameters before and after noise reduction. These feature parameters include noise intensity (NI) [5], spatial-dependence recurrence sample entropy (sdrSampEn) [34] and DE.

NI
Assuming a time series of N samples in length, that is x(n) = x(1), x(2), . . . , x(N) . Its NI can be approximated by the standard deviation, i.e.,: x(n) represents the mean of the time series.

SdrSampEn
The traditional sample entropy (SampEn) is a tool to measure the complexity of time series and is widely used in chaotic signal analysis. However, the SampEn does not take sequential information, which is inherently useful, into its calculation of sample similarity. So, Pham et al. proposed a tool that is based on the mathematical principle of the SampEn and enables the capture of sequential information of a time series in the context of spatial dependence provided by the binary-level co-occurrence matrix (BLCM) of a recurrence plot (RP). The sdrSampEn algorithm can be described as follows [33]: Step 1. Given a time series and model parameters embedding dimension m, and spatial offset δ; Step 2. Use time series, m, to compute the RP by using the Chebyshev distance; Step 3. Use offset δ to construct the BLCM of the RP to obtain S m (δ); Step 4. Set m = m + 1, repeat steps 2-3 to obtain S m+1 (δ); Step 5. Finally, sdrSampEn(m, δ) = − log

Performance Analysis of MCEEMDAN
In this section, we will use a simulated signal to test the performance of MCEEMDAN and compare the proposed method with EMD, EEMD, CEEMD, and CEEMDAN.
where s1, s2, and s3 represent the three components of S. The time-domain waveform of simulated signal is shown in Figure 2. For EEMD, CEEMD, CEEMDAN and its modified version, the noise standard deviation and the ensemble size are set to 0.2 and 60, respectively. Figure 3 shows the decomposed results of simulation signals. Since the decomposition results of EEMD and CEEMD are similar, only CEEMD is shown in Figure 3b. Step 2. Use time series, m, to compute the RP by using the Chebyshev distance; Step 3. Use offset δ to construct the BLCM of the RP to obtain ( ) m S δ ; Step 4. Set 1 m m = + , repeat steps 2-3 to obtain 1 ( ) m S δ + ; Step 5. Finally, Pham et al.
[34] used SdrSampEn to measure the uncertainty of Lorenz signal. At the same time, it is pointed out that SdrSampEn have a better discriminative ability in measuring the irregularity of chaotic time series. Detailed calculation of SdrSampEn can be seen in [34].

Performance Analysis of MCEEMDAN
In this section, we will use a simulated signal to test the performance of MCEEMDAN and compare the proposed method with EMD, EEMD, CEEMD, and CEEMDAN.
s , 2 s , and 3 s represent the three components of S . The time-domain waveform of simulated signal is shown in Figure 2. For EEMD, CEEMD, CEEMDAN and its modified version, the noise standard deviation and the ensemble size are set to 0.2 and 60, respectively. Figure 3 shows the decomposed results of simulation signals. Since the decomposition results of EEMD and CEEMD are similar, only CEEMD is shown in Figure 3b.  Figure 3 shows obvious mode mixing of the IMF1 simulated with EMD. Although CEEMD improved the mode mixing, it divided 2 s into two IMFs with large amplitude changes. The IMF4 and IMF5 of CEEMDAN method have the same frequency components, and there are some pseudo components that have no physical meaning. Figure 3d shows that the modified CEEMDAN can restore the time-domain waveform of the simulated signal better.
In order to quantitatively evaluate the effect of these decomposition algorithms, index of orthogonality (IO), consuming time (CT), and reconstruction error order of magnitude (REOM) are calculated [18]. Where CT is the running time statistics of each decomposition algorithm. IO can be used to measure the accuracy of decomposition algorithms. If IO is smaller, then the decomposition accuracy is higher. Otherwise, the decomposition accuracy is lower. IO can be defined as [18]:  Figure 3 shows obvious mode mixing of the IMF1 simulated with EMD. Although CEEMD improved the mode mixing, it divided s2 into two IMFs with large amplitude changes. The IMF4 and IMF5 of CEEMDAN method have the same frequency components, and there are some pseudo components that have no physical meaning. Figure 3d shows that the modified CEEMDAN can restore the time-domain waveform of the simulated signal better.
In order to quantitatively evaluate the effect of these decomposition algorithms, index of orthogonality (IO), consuming time (CT), and reconstruction error order of magnitude (REOM) are calculated [18]. Where CT is the running time statistics of each decomposition algorithm. IO can be used to measure the accuracy of decomposition algorithms. If IO is smaller, then the decomposition accuracy is higher. Otherwise, the decomposition accuracy is lower. IO can be defined as [18]: where x(l) represents the original time series and its sample length is L. The original time series is decomposed into K IMFs. IMF i (l) and IMF j (l) represent the ith and jth IMF, respectively. The decomposition performances of these decomposition algorithms are shown in Table 1.
where ( ) x l represents the original time series and its sample length is L. The original time series is decomposed into K IMFs. ( ) i IMF l and ( ) j IMF l represent the ith and jth IMF, respectively. The decomposition performances of these decomposition algorithms are shown in Table 1.   Table 1 shows that the CEEMDAN method overcomes the mode mixing problem to a certain extent, but because of the multiple addition of white Gaussian noise, the CT is the largest. The above results show that the MCEEMDAN method has a small IO, and its computational efficiency is higher than that of CEEMDAN.

Comparison of SE, PE, and DE
To compare the performance of SE, PE and DE, give an example. We added white Gaussian noise with different signal-to-noise ratios (SNRs) to the sinusoidal signal, and then calculated the SE, PE, and DE of the simulation signal separately. The simulation signals ( ) S t can be determined by: where ( ) n t represents the white Gaussian noise, and the sampling frequency is 100 Hz. ( ) S t is added with −10 dB to 35 dB white Gaussian noise with an interval of 5 dB. The time-domain waveform of the noisy simulation signal is shown in Figure 4. The values of the three types of entropy are shown in Figure 5, and the CT of each entropy algorithm is also listed in Table 2.    Table 1 shows that the CEEMDAN method overcomes the mode mixing problem to a certain extent, but because of the multiple addition of white Gaussian noise, the CT is the largest. The above results show that the MCEEMDAN method has a small IO, and its computational efficiency is higher than that of CEEMDAN.

Comparison of SE, PE, and DE
To compare the performance of SE, PE and DE, give an example. We added white Gaussian noise with different signal-to-noise ratios (SNRs) to the sinusoidal signal, and then calculated the SE, PE, and DE of the simulation signal separately. The simulation signals S(t) can be determined by: where n(t) represents the white Gaussian noise, and the sampling frequency is 100 Hz. S(t) is added with −10 dB to 35 dB white Gaussian noise with an interval of 5 dB. The time-domain waveform of the noisy simulation signal is shown in Figure 4. The values of the three types of entropy are shown in Figure 5, and the CT of each entropy algorithm is also listed in Table 2.      Figure 4 shows that when the SNR is less than 0 dB, the sinusoidal signal is almost overwhelmed by noise. When the SNR is greater than 5 dB, the influence of noise on the sinusoidal signal is small. Figure 5 shows that when the SNR is less than 0 dB, the DE is greater than 0.88. With the increase of SNR, the DE of simulation signal gradually decreases, and finally approaches 0.6. When SNR is less than 15 dB, the PE of the simulation signal is basically greater than 0.95. Compared to PE, SE and DE are more robust to noise, so they can better describe the complexity of noise sinusoidal signals. However, the value of SE is between 0 and 3.5, which limits its use because it is difficult to normalize the SE. In addition, from the CT of the three algorithms, SE has the lowest computational efficiency. In summary, DE can not only accurately measure the complexity of time series, but also has the highest computational efficiency, so we introduce DE into the denoising of SR-N signals. In the following, if the DE of the IMF is greater than 0.88, it is determined to be a noise IMF.    Figure 4 shows that when the SNR is less than 0 dB, the sinusoidal signal is almost overwhelmed by noise. When the SNR is greater than 5 dB, the influence of noise on the sinusoidal signal is small. Figure 5 shows that when the SNR is less than 0 dB, the DE is greater than 0.88. With the increase of SNR, the DE of simulation signal gradually decreases, and finally approaches 0.6. When SNR is less than 15 dB, the PE of the simulation signal is basically greater than 0.95. Compared to PE, SE and DE are more robust to noise, so they can better describe the complexity of noise sinusoidal signals. However, the value of SE is between 0 and 3.5, which limits its use because it is difficult to normalize the SE. In addition, from the CT of the three algorithms, SE has the lowest computational efficiency. In summary, DE can not only accurately measure the complexity of time series, but also has the highest computational efficiency, so we introduce DE into the denoising of SR-N signals. In the following, if the DE of the IMF is greater than 0.88, it is determined to be a noise IMF.

Denoising of Chen's Chaotic Signal
Chen's system can be expressed as: when the parameters a = 35, b = 3 and c = 28, the system is in typical chaotic state. The equation is integrated by using a fourth-order Runge-Kutta method with a fixed step size of 0.01, and the initial values of the equation are x(0) = −1, y(0) = 0, z(0) = 1. The y component signal with a length of 2048 points is selected as a chaotic signal, and the signal are added the white Gaussian noise with SNR = 11 dB. The time-domain waveforms of the noisy signal and the noiseless signal are shown in Figure 6.  First, the noisy Chen's signal is decomposed by MCEEMDAN. Figure 7 shows that the Chen's signal is decomposed into 12 IMFs. It can be seen that the frequency decreases with the increase of IMF order. The DE of each IMF is calculated, and DE values are listed in Table 3. First, the noisy Chen's signal is decomposed by MCEEMDAN. Figure 7 shows that the Chen's signal is decomposed into 12 IMFs. It can be seen that the frequency decreases with the increase of IMF order. The DE of each IMF is calculated, and DE values are listed in Table 3. First, the noisy Chen's signal is decomposed by MCEEMDAN. Figure 7 shows that the Chen's signal is decomposed into 12 IMFs. It can be seen that the frequency decreases with the increase of IMF order. The DE of each IMF is calculated, and DE values are listed in Table 3.   IMF1   IMF2   IMF3   IMF4   IMF5   IMF6   IMF7   IMF8   IMF9   IMF10   IMF11   As can be seen from Table 3, since DEs of IMF1 to IMF3 are both greater than 0.88, they are judged as noise IMFs. The last IMF with DE greater than 0.65 is IMF5, so IMF4 and IMF5 are determined to be noise-dominated IMFs. IMF6-IMF12 are pure IMFs. The noise-dominated IMFs are denoised by the interval thresholding method. The time-domain waveforms of IMF4 and IMF5 before and after noise reduction are shown in Figure 8. It can be seen that in IMF4, the lower amplitude intervals are filtered out as noise, while the other intervals are retained. For IMF5, the extreme values in most of the intervals are greater than the threshold, so there is no significant difference before and after noise reduction.  As can be seen from Table 3, since DEs of IMF1 to IMF3 are both greater than 0.88, they are judged as noise IMFs. The last IMF with DE greater than 0.65 is IMF5, so IMF4 and IMF5 are determined to be noise-dominated IMFs. IMF6-IMF12 are pure IMFs. The noise-dominated IMFs are denoised by the interval thresholding method. The time-domain waveforms of IMF4 and IMF5 before and after noise reduction are shown in Figure 8. It can be seen that in IMF4, the lower amplitude intervals are filtered out as noise, while the other intervals are retained. For IMF5, the extreme values in most of the intervals are greater than the threshold, so there is no significant difference before and after noise reduction. The pure IMFs and the processed noise-dominated IMFs are restructured to obtain the denoised Chen's signal. In order to prove the superiority of MCEEMDAN-DE-IT method, this method is compared with EMD-DE-IT, EEMD-DE-IT, and CEEMD-DE-IT. Among them, EMD-DE-IT uses the EMD-IT code provided in [26], and noise evaluation criteria is the same as this paper. For the sake of evaluating the performance of these methods, the SNR and root mean square error (RMSE) are adopted and defined as follows [28]: where x and x are the noiseless signal and the denoised signal; N represents the length of the signal. The time-domain waveform before and after noise reduction for noisy Chen's signal when SNR = 11 dB are shown in Figure 9.  The pure IMFs and the processed noise-dominated IMFs are restructured to obtain the denoised Chen's signal. In order to prove the superiority of MCEEMDAN-DE-IT method, this method is compared with EMD-DE-IT, EEMD-DE-IT, and CEEMD-DE-IT. Among them, EMD-DE-IT uses the EMD-IT code provided in [26], and noise evaluation criteria is the same as this paper. For the sake of evaluating the performance of these methods, the SNR and root mean square error (RMSE) are adopted and defined as follows [28]: where x andx are the noiseless signal and the denoised signal; N represents the length of the signal. The time-domain waveform before and after noise reduction for noisy Chen's signal when SNR = 11 dB are shown in Figure 9. The pure IMFs and the processed noise-dominated IMFs are restructured to obtain the denoised Chen's signal. In order to prove the superiority of MCEEMDAN-DE-IT method, this method is compared with EMD-DE-IT, EEMD-DE-IT, and CEEMD-DE-IT. Among them, EMD-DE-IT uses the EMD-IT code provided in [26], and noise evaluation criteria is the same as this paper. For the sake of evaluating the performance of these methods, the SNR and root mean square error (RMSE) are adopted and defined as follows [28]: where x and x are the noiseless signal and the denoised signal; N represents the length of the signal. The time-domain waveform before and after noise reduction for noisy Chen's signal when SNR = 11 dB are shown in Figure 9.    Figure 9 shows that the EMD-DE-IT method loses some important information of the noiseless signal, and its noise reduction performance is the worst. Both EEMD-DE-IT and CEEMD-DE-IT have good noise reduction effects, and CEEMD-DE-IT is superior to EEMD-DE-IT, but there is still some residual noise in some details. In addition, Table 4 shows that the proposed method has a better noise reduction effect for the Chen's signal with SNR = 11 dB. As shown in Figure 10, in order to compare noise reduction effects more intuitively, we also display the phase diagrams before and after noise reduction. Table 4. Denoising result of the Chen's signal when SNR is 11 dB. By comparing the phase diagram in Figure 10, it can be seen that the attractor of the noisy Chen's signal is annihilated by the noise, and the fractal features are hardly observed. After the noise reduction by the EMD-DE-IT method, the shape of the attractor is quite different from the noiseless signal, indicating that some important information is actually lost. Both EEMD-DE-IT and CEEMD-DE-IT methods roughly restore the phase diagram of the pure attractor, but there are still deformed orbits. Since the intensity of the noise determines the damage degree of the attractor's self-similar structure, it indicates that the denoised signal still has some noise. Using the proposed method, the geometry of the attractor structure is more regular and the degree of self-similarity is higher. The phase diagram after noise reduction is closer to that of the noiseless signal, which indicates that the proposed method can effectively filter out the noise of the noisy signal. Therefore, Figure 10 objectively shows that the denoising effect of the proposed method is better than the other three methods. Figure 9 shows that the EMD-DE-IT method loses some important information of the noiseless signal, and its noise reduction performance is the worst. Both EEMD-DE-IT and CEEMD-DE-IT have good noise reduction effects, and CEEMD-DE-IT is superior to EEMD-DE-IT, but there is still some residual noise in some details. In addition, Table 4 shows that the proposed method has a better noise reduction effect for the Chen's signal with SNR = 11 dB. As shown in Figure 10, in order to compare noise reduction effects more intuitively, we also display the phase diagrams before and after noise reduction. By comparing the phase diagram in Figure 10, it can be seen that the attractor of the noisy Chen's signal is annihilated by the noise, and the fractal features are hardly observed. After the To further compare the noise reduction effect, we added −5 dB, 0 dB, 5 dB, 10 dB, and 15 dB white Gaussian noise in the Chen's signal, respectively. Figure 11 presents the comparison of the SNR and the RMSE of the above four methods. The results are averaged from 10 random experiments. It can be seen that the performance of EEMD-DE-IT is similar to that of CEEMD-DE-IT, and the performance of EMD-IT method is poor at high SNR. Experiments show that the proposed method has better noise suppression ability for the Chen's signal with low SNR. From the RMSE point of view, the reconstruction accuracy of this method is better than other methods at low SNR. noise reduction by the EMD-DE-IT method, the shape of the attractor is quite different from the noiseless signal, indicating that some important information is actually lost. Both EEMD-DE-IT and CEEMD-DE-IT methods roughly restore the phase diagram of the pure attractor, but there are still deformed orbits. Since the intensity of the noise determines the damage degree of the attractor's self-similar structure, it indicates that the denoised signal still has some noise. Using the proposed method, the geometry of the attractor structure is more regular and the degree of self-similarity is higher. The phase diagram after noise reduction is closer to that of the noiseless signal, which indicates that the proposed method can effectively filter out the noise of the noisy signal. Therefore, Figure 10 objectively shows that the denoising effect of the proposed method is better than the other three methods.

EMD-DE-IT EEMD-DE-IT CEEMD-DE-IT
To further compare the noise reduction effect, we added −5 dB, 0 dB, 5 dB, 10 dB, and 15 dB white Gaussian noise in the Chen's signal, respectively. Figure 11 presents the comparison of the SNR and the RMSE of the above four methods. The results are averaged from 10 random experiments. It can be seen that the performance of EEMD-DE-IT is similar to that of CEEMD-DE-IT, and the performance of EMD-IT method is poor at high SNR. Experiments show that the proposed method has better noise suppression ability for the Chen's signal with low SNR. From the RMSE point of view, the reconstruction accuracy of this method is better than other methods at low SNR.

Denoising of Real Ship Radiated Noise Signal
In this paper, three different types of SR-N signals are selected as sample data, namely the Ship-I, the Ship-II, and the Ship-III. The SR-N signals were measured in the South China Sea. Each type of SR-N signals has 100 sample data. Each sample length is 2048 points and sampling interval is 0.05 ms. Before the noise reduction experiment, the sample data have been filtered, normalized, and sampled. SR-N signals and denoised SR-N signals and their attractor trajectories are shown in Figures 12-14.

Denoising of Real Ship Radiated Noise Signal
In this paper, three different types of SR-N signals are selected as sample data, namely the Ship-I, the Ship-II, and the Ship-III. The SR-N signals were measured in the South China Sea. Each type of SR-N signals has 100 sample data. Each sample length is 2048 points and sampling interval is 0.05 ms. Before the noise reduction experiment, the sample data have been filtered, normalized, and sampled. SR-N signals and denoised SR-N signals and their attractor trajectories are shown in Figures 12-14.

Denoising of Real Ship Radiated Noise Signal
In this paper, three different types of SR-N signals are selected as sample data, namely the Ship-I, the Ship-II, and the Ship-III. The SR-N signals were measured in the South China Sea. Each type of SR-N signals has 100 sample data. Each sample length is 2048 points and sampling interval is 0.05 ms. Before the noise reduction experiment, the sample data have been filtered, normalized, and sampled. SR-N signals and denoised SR-N signals and their attractor trajectories are shown in  x(n+τ)   It can be seen from the time-domain waveforms in Figures 12-14 that the original SR-N signals contain ocean background noise, and the high-frequency noise is filtered out after noise reduction. The phase diagram shows that the attractor of the original signal appears as a disorderly pseudo-random property. After noise reduction, the orbits of attractors become smooth, clear and regular. The proposed method does reduce the noise, enhances the inherent deterministic component of the SR-N signal, and reveals the inherent dynamic characteristics of chaotic system.
We calculated the feature parameters of the SR-N signal before and after noise reduction. Table  5 shows that the feature parameters after noise reduction have been improved, indicating that the complexity of the SR-N signal has been reduced. In summary, the MCEEMDAN-DE-IT method can be used for noise reduction of real SR-N signals.

Conclusions
An SR-N signal denoising method based on modified CEEMDAN, dispersion entropy, and interval thresholding is proposed. The original signal is decomposed into a series of IMFs by MCEEMDAN. According to DE, the IMFs are divided into three parts: noise IMFs, noise-dominated IMFs and pure IMFs. The noise-dominated IMFs are processed using interval thresholding. Finally, It can be seen from the time-domain waveforms in Figures 12-14 that the original SR-N signals contain ocean background noise, and the high-frequency noise is filtered out after noise reduction. The phase diagram shows that the attractor of the original signal appears as a disorderly pseudo-random property. After noise reduction, the orbits of attractors become smooth, clear and regular. The proposed method does reduce the noise, enhances the inherent deterministic component of the SR-N signal, and reveals the inherent dynamic characteristics of chaotic system.
We calculated the feature parameters of the SR-N signal before and after noise reduction. Table 5 shows that the feature parameters after noise reduction have been improved, indicating that the complexity of the SR-N signal has been reduced. In summary, the MCEEMDAN-DE-IT method can be used for noise reduction of real SR-N signals.

Conclusions
An SR-N signal denoising method based on modified CEEMDAN, dispersion entropy, and interval thresholding is proposed. The original signal is decomposed into a series of IMFs by MCEEMDAN. According to DE, the IMFs are divided into three parts: noise IMFs, noise-dominated IMFs and pure IMFs. The noise-dominated IMFs are processed using interval thresholding. Finally, the noise-dominated IMFs and the pure IMFs are reconstructed to obtain the denoised SR-N signal. The experimental results of the simulated signal and the real SR-N signal show that the denoising method has the following advantages: