Guo Shusen Yu Xianwen Long Fengyang Wang Jiafu

(School of Transportation, Southeast University, Nanjing 211189, China)

Abstract：A filter method that combines ensemble empirical modal decomposition(EEMD)and wavelet analysis methods was proposed to separate and correct the global navigation satellite system(GNSS)multipath error more effectively.In this method, the GNSS signal is first decomposed into several intrinsic mode functions(IMFs)and a residual through EEMD.Then, the IMFs and residual are classified into noise terms, mixed terms, and useful terms according to a combined classification criterion.Finally, the mixed term denoised by wavelet and the useful term are reconstructed to obtain the multipath error and thus enable an error correction model to be built.The measurement data provided by the Curtin GNSS Research Center were used for processing and analysis.Results show that the proposed method can separate multipath error from GNSS data to a great extent, thereby effectively addressing the defects of EEMD and wavelet methods on multipath error weakening.The error correction model established with the separated multipath error has a higher accuracy and provides a certain reference value for research on related signal processing.

Key words：ensemble experience modal decomposition(EEMD); wavelet analysis; multipath error; global navigation satellite system(GNSS)

G lobal navigation satellite system(GNSS)uses static measurement for deformation monitoring[1-2].Weakening the multipath effect through difference is difficult because of its low correlation between stations[3-4].The repetitive correction method has an efficient weakening effect on the multipath error in GNSS deformation monitoring because of the diurnal repeatability of the multipath effect[5-6], in which multipath error separation plays a key role in the whole method.Satirapod et al.[7]and Su et al.[8]used wavelet to separate the multipath error from GPS signal data, and their results verify the effectiveness of wavelet.However, the wavelet effect is limited by choice of wavelet basis, threshold, and decomposition order, which are usually determined by experience[6,9].Huang et al.[10]proposed empirical mode decomposition(EMD)to adaptively decompose a signal into a series of intrinsic mode functions(IMFs)and a residual, allowing useful information to be extracted quickly from the decomposition without prior information[11-12].Dai et al.[13]used EMD to separate GPS multipath errors, and they found that EMD can separate multipath errors more directly and effectively than wavelet.Yan et al.[14]proposed a combination of EMD and recursive least squares(RLS)filtering to separate the multipath error from GNSS data, maximizing the advantages of EMD and RLS.However, EMD is easily susceptible to the mode aliasing problem, resulting in degraded decomposition accuracy[15].Wu et al.[16]proposed ensemble empirical mode decomposition(EEMD)to overcome mode aliasing, but the effect of EEMD is limited by the amount of white noise added[17].Therefore, this paper proposed a novel method that combines wavelet and EEMD to weaken the multipath errors more effectively.In this paper, the principle of the proposed method is introduced first, and then its results, as evaluated by GNSS data, are presented.

1 Multipath Effect

Fig.1 shows that the multipath effect is caused by the interference between direct and reflected signals.With the assumption that the receiver receives a set of signals from a satellite, the superimposed signal formed by the reflected signal and the direct signal in this set of signals can be expressed as

S=AScos(ωt+φ)

(1)

whereASis the amplitude of the superimposed signal;ωis the angular frequency;φis the multipath error, and its formula is

(2)

whereαis the reflection factor of the reflector, andΔis the phase delay caused by the path difference between the reflected signal and the direct signal, which can be obtained according to Fig.1.

Fig.1 Principle of multipath effect

(3)

whereλis the wavelength of the signal;θis the reflection angle; andHis the receiver height.

The typical frequencyfof the multipath effect can be expressed as

(4)

As a result of signal attenuation, the reflected signal can be neglected at a reflected distance greater than 50 m, implying a low multipath effect in the frequency domain[3].The operation cycle of GNSS satellites is basically maintained at 12 sidereal times.The multipath effect for two consecutive days shows a strong correlation[4], i.e., diurnal repeatability, when the surrounding environment of the station receiver changes slightly.

2 Improved EEMD-Wavelet Combined Filter

2.1 Wavelet analysis

According to the principle of wavelet, the wavelet basis is stretched and shifted, and the wavelet basis is used to fit the signal to be processed at multiple scales.The wavelet transform for a signal seriesf(t)can be expressed as[7]

(5)

whereais the stretch factor;bis the shift factor;ψ(t)is the wavelet basis; andψ*(t)is the conjugate function ofψ(t).

f(t)can be reconstructed by using inverse wavelet transform, which is expressed as

(6)

whereCψis the admissibility condition of the wavelet basis.

The following are the specific steps for separating the multipath error from the GNSS signal using wavelet:

1)The wavelet basis and decomposition layers are determined and are used to decompose the GNSS signal.

2)The rule for estimating threshold is determined and is used for soft thresholding of the high-frequency wavelet coefficients of each layer; the low-frequency wavelet coefficients of each layer are not processed.

3)The high-frequency wavelet coefficients after soft thresholding and the low-frequency wavelet coefficients are combined to reconstruct the useful signal, i.e., multipath error.

Knowing the characteristics of the signal in advance and determining the optimal wavelet parameters according to experience are generally necessary because wavelet lacks adaptive processing ability.Typically, the GNSS signal exhibits multiscale characteristics because of the superposition of various noises, leading to a significantly increased workload of multipath error separation by wavelet.

2.2 Ensemble empirical mode decomposition

EEMD is an improved EMD method based on the characteristics of non-correlation, zero-mean, and uniform distribution of Gaussian white noise in time-frequency space.The EEMD effectively overcomes mode mixing with the addition of Gaussian white noise and multiple EMD processing[16].

The following are the specific steps for separating the multipath error from the GNSS signal by using EEMD:

1)Gaussian white noise is added to the GNSS signal to obtain the signal to be processed.

Ym(t)=Y(t)+nm(t)

(7)

wherenm(t)is the Gaussian white noise series, andY(t)is the GNSS signal series.

2)The signal series are decomposed into a series of IMFs and one residual by EMD[10].

(8)

wherenis the order of EMD decomposition, i.e., the number of IMFs.

3)Steps 1 and 2 are repeatedNtimes to obtainNsets of IMFs and residuals.

4)The average values of IMFs and residuals are calculated as the final decomposition result.

(9)

(10)

5)The useful signal(multipath error)is reconstructed by analyzing the frequency and amplitude characteristics of the decomposition result.

According to Chen et al.[17], EEMD works best when the amplitude factor of white noise is 0.01 to 0.5, and the number of EMD processing is 100 to 300.However, theoretical and data support for determining the above two parameters in multipath error separation is lacking, which most likely will lead to an unsatisfactory separation.

2.3 Improvement of EEMD-wavelet combined filter

In accordance with the characteristics and shortcomings of wavelets and EEMD, a combined filter was proposed for GNSS multipath error weakening, and its process is shown in Fig.2.

Fig.2 Process of the proposed method

The steps of the process are as follows:

1)The GNSS signal is decomposed into several IMFs and a residual using EEMD, and the decomposition is subdivided into noiseY(t)noise, mixedY(t)mixedand useful termsY(t)useful.

(11)

wheret=1,2,…,m,mis the epochs.

2)Y(t)noiseis abandoned, and wavelet is used to filterY(t)mixed.Then, theY(t)mixedafter filtering andY(t)usefulare reconstructed to obtain the useful signal, i.e., the multipath error.

(12)

In Eq.(11),k1andk2are the combined classification criterion indexes proposed in this paper, which are calculated as follows.

The reconstructed signal can be expressed as

(13)

wherek=1, 2, …,n.

The continuous mean square error(CMSE)criterion[11]is used to determine the first indexk1of the classification criteria.CMSE criterion is defined as follows:

(14)

wherek=1, 2, …,n-1;Nis the signal length.

Eq.(14)indicates that CMSE measures the squared Euclidean distance between two consecutive reconstructions of the signal, which is equivalent to the energy density of thek-th IMF.With this quantity, the IMF order where the first significant change in energy occurs can be determined because the energy of the mixed and useful terms is much higher than that of noise.kis taken as the starting point of the mixed term when the value of thek-th CMSE is the first local minimum,k1=k.

(15)

The second indexk2is determined by the energy coefficient, which is defined by the product of energy density and average period, where the calculation for thek-th IMF’s product of energy density and average period is as follows

(16)

(17)

(18)

whereNis the length of thek-th IMF;Okis the number of extreme points of thek-th IMF.

The energy coefficientCkis defined as[12]

(19)

wherek=2,3,…,n.

The product of the IMF dominated by white noise is a constant.Ckis greater than the given thresholdσwhenPkincreases geometrically compared with the previous one, which implies thek-th IMF is non-constant whilst the previous IMFs are mainly dominated by noise.Accordingly,kserves as the starting point of the useful term,k2=k, as shown in Fig.3.In this paper,σ=3[13].

Fig.3 Determination process of k2

3 GNSS Measured Data Analysis

3.1 Data sources and analysis

The analysis in this subsection employed the GNSS data from three stations(CUT0, CUTB and CUTC)from July 5 to 7, 2020(DOY187, DOY188, DOY189).They were collected at 30 s intervals with millimeter accuracy and can be downloaded for free from the Curtin GNSS Research Center.First, we calculated the baseline vectors of CUTB-CUT0 and CUTB-CUTC by using the processing strategy shown in Tab.1.Then, we obtained the residual series by subtracting the true values of the baseline vectors.

Tab.1 GNSS data processing strategy

The residual series can be considered to consist of only multipath error and noise because CUTB-CUT0 and CUTB-CUTC are ultra-short baselines with lengths of 4.27 and 6.15 m, respectively.In this work, we considered only the residual series with a length of 2 500 epochs in the vertical direction because the height accuracy is the most interesting and important part of the deformation monitoring, where each residual series has a length of 2 500 epochs.We also calculated the trend of the residual series by using the moving average method with a moving window of 50 epochs.The residual series and their trend lines are shown in Fig.4.

Fig.4 shows that the residual series contains high-frequency noise and low-frequency multipath error, and some abnormal jump phenomena occur randomly.We calculated the Pearson’s correlation coefficientRof the residual series in Fig.4 for correlation analysis, and the results are shown in Tab.2.

(a)

Tab.2 Correlation analysis results for residual series

According to Tab.2, theRvalues of the residual series for two adjacent days are maintained at around 0.8, and theRvalues of the residual series with an interval of 1 d are maintained at around 0.68.This condition indicates the existence of significant diurnal repeatability among the residual series.

3.2 Multipath error separation and analysis

In this subsection, wavelet, EEMD, and the proposed method were used to separate the multipath from the residual series, and their results were compared.Wavelet adopted the Sym6 wavelet basis with five decomposition layers and the soft threshold function based on the Heursure rule.EEMD combined the mean of the standardized accumulated mode criterion[6].The proposed method is consistent with the above in the wavelet part.We used the root-mean-square error(RMSE)to evaluate the separation effect.

(20)

whereNis the length of the residual series;y(t)is the residual series; andu(t)is the filtered noise series in this subsection.

The separated multipath error series is shown in Fig.5.The trend of the multipath error series obtained by the proposed method is closer to that obtained by wavelet; their multipath errors are smoother than those of EEMD.The multipath error series obtained by wavelet exhibits some significant jump phenomena, whereas the multipath error series obtained by the proposed method does not exhibit such jump phenomena.

(a)

The RMSE values of residual series before and after filtering are shown in Tab.3.The ratio of RMSE before and after filtering of the three methods remains at 0.6-0.9, indicating that the multipath effect dominates in the residual series.The RMSE values of wavelet and the proposed method are lower than those of EEMD, with some RMSE values of the proposed method being slightly greater than those of wavelet.

Tab.3 RMSE of residual series before and after filtering

The aforementioned analysis cannot fully determine which approach is the most effective in multipath error separation because the true value of the multipath error in the residual series used is unknown.Therefore, further analysis was performed on the separation effect of the three methods according to the low-frequency and diurnal repeatability of the multipath effect.

Spectrum analysis was conducted for the residual series and its multipath error series, and the results of DOY187 are shown in Fig.6.The spectrogram of the residual series shows the presence of obvious noise in the whole frequency domain, while the multipath error is mainly concentrated in the low-frequency part of the rectangular window and its amplitude value is significantly greater than that of noise.The spectrum analysis results of the three methods in the window suggests that the multipath error of wavelet is mainly distributed at 0-0.02 and 0.06-0.2 Hz, and the multipath error of EEMD is mainly distributed at 0-0.06Hz.The multipath error of the proposed method is the lowest among the three, being mainly distributed at 0-0.02Hz.

(a)

The correlation analysis of the multipath error series is shown in Tab.4.The multipath effect is highly reproducible for the two adjacent days.TheRvalues of the proposed method remain at 0.66-0.85, where theRvalue between the second and third days of CUTB-CUTC is slightly lower than that of wavelet.

Tab.4 Correlation analysis results for multipath errors

In summary,the proposed method has better performance than wavelet and EEMD in multipath error separation, as indicated by an analysis based on intuition, low-frequency, and diurnal repeatability.

3.3 Multipath error correction

The correction effect of the proposed method on multipath error was determined by first using the multipath error series of the first day for modeling to correct the residual series of the second and third days.Then, the multipath error series for the first two days was used for mean modeling to correct the third day.In this subsection, we used Eq.(20)to calculate the RMSE to evaluate the correction effect, whereu(t)is the separated multipath error series.

The RMSE values of the residual series after multipath error correction are shown in Tab.5.The correction effect for the second day is better than that for the third day, indicating that repetitive modeling based on the three methods can effectively correct the multipath error at a relatively short time interval of multipath error.The correction effect of the multipath error begins to decline as the time interval increases, while the mean modeling further corrects the multipath error compared with the former.Thus, repetitive modeling based on mean processing can effectively suppress the effect of time interval extension.The RMSE values after correction by the proposed method are lower than those of wavelet and EEMD, and we can calculate that the proposed method attains about 21.07% and 41.90% multipath error correction in the residual series of CUTB-CUT0 and CUTB-CUTC, respectively.Compared with wavelet and EEMD, the proposed method improves the correction of the residual series of CUTB-CUT0 and that of CUTB-CUTC by 5.50% and 4.45% on average, respectively.

Tab.5 RMSE of residual series before and after multipath correction

4 Conclusions

1)The shortcomings of EEMD and wavelet on multipath weakening were analyzed in this work.A filter method that combines EEMD and wavelet was proposed for multipath error weakening.

2)According to the GNSS data analysis, the evaluation indices show that the effect of the proposed method on multipath separation and correction is better than that of EEMD and wavelet.

3)The proposed method achieved good results on multipath error weakening but has an unclear improvement effect on multipath error correction.Its effectiveness in practical application is limited to a certain extent, and thus it needs further improvement in terms of theory and algorithm.