Next Article in Journal
Optimal Linear Error Correcting Delivery Schemes for Two Optimal Coded Caching Schemes
Previous Article in Journal
Thermodynamic Analysis of a Hybrid Trigenerative Compressed Air Energy Storage System with Solar Thermal Energy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Denoising Method for Fiber Optic Gyroscope Based on Variational Mode Decomposition and Beetle Swarm Antenna Search Algorithm

Collage of Automation, Harbin Engineering University, Harbin 150001, China
*
Author to whom correspondence should be addressed.
Entropy 2020, 22(7), 765; https://doi.org/10.3390/e22070765
Submission received: 13 June 2020 / Revised: 30 June 2020 / Accepted: 11 July 2020 / Published: 13 July 2020

Abstract

:
Fiber optic gyroscope (FOG) is one of the important components of Inertial Navigation Systems (INS). In order to improve the accuracy of the INS, it is necessary to suppress the random error of the FOG signal. In this paper, a variational mode decomposition (VMD) denoising method based on beetle swarm antenna search (BSAS) algorithm is proposed to reduce the noise in FOG signal. Firstly, the BSAS algorithm is introduced in detail. Then, the permutation entropy of the band-limited intrinsic mode functions (BLIMFs) is taken as the optimization index, and two key parameters of VMD algorithm, including decomposition mode number K and quadratic penalty factor α , are optimized by using the BSAS algorithm. Next, a new method based on Hausdorff distance (HD) between the probability density function (PDF) of all BLIMFs and that of the original signal is proposed in this paper to determine the relevant modes. Finally, the selected BLIMF components are reconstructed to get the denoised signal. In addition, the simulation results show that the proposed scheme is better than the existing schemes in terms of noise reduction performance. Two experiments further demonstrate the priority of the proposed scheme in the FOG noise reduction compared with other schemes.

1. Introduction

Inertial navigation is an autonomous navigation method. Because of its uniqueness, the inertial navigation system (INS) is widely used in many fields [1,2,3]. Fiber optic gyroscope (FOG) is one of the most common hardware for inertial system, it has the advantages of high measurement accuracy and large measurement range [4,5]. However, due to the influence of external environment and signal sampling conditions, the output signal of FOG often contains a lot of random noise [6,7]. It will greatly affect the accuracy of INS. Therefore, it is necessary to suppress the random noise of FOG signal.
In general, there are three different ways to deal with random noise of FOG signal. The first one is to establish autoregressive moving average (ARMA) model based on the random noise sequence, and then the model is optimized by filtering methods [8,9]. Although this method performs well under the premise of accurate noise model is established, it is difficult to get accurate random noise model in practical application, so the noise reduction effect will be greatly affected [10,11].
Using filtering algorithm to reduce noise directly is the second common method. Narasimhappa et al. proposed a robust adaptive Kalman filtering algorithm to filter the FOG signal directly [12]. In this method, a new covariance matrix of measurement noise is designed to improve the filtering effect of the adaptive Kalman filtering algorithm. He et al. proposed a hybrid algorithm which combining the advantages of adaptive chirp mode pursuit (ACMP) algorithm and adaptive multiscale Savitzky–Golay filter (AMSGF) algorithm [13]. In this method, the ACMP algorithm is uesd to denoise the FOG signal roughly, and then the AMSGF algorithm is utilized to filter it further. All in all, although the second common method does not need to establish the mathematical model of noise sequence, the effect of the denoising is not ideal because of the relatively simple denoising process.
The third way to deal with FOG signal is to decompose the original signal first, and then analyze the decomposed components. The common methods for decomposing signal include wavelet transform (WT) [14] and empirical mode decomposition (EMD) [15]. WT can decompose the signal into components of different frequency bands. Zhang et al. put forward an improved wavelet threshold transform to process the random noise of FOG signal [16]. The experimental results show that this method can get lower mean square error. Ma et al. applied the optimized WT filter to the magnetic levitation gyroscope (GAT) and successfully reduced the vibration noise in actual use [17]. The WT method can effectively reduce the random noise in FOG signal, but this method needs to preset wavelet basis function and decomposition level, so the adaptability is not good. EMD can decompose the original signal into a series of intrinsic mode functions (IMFs) by analyzing the local characteristic time scale of the signal. Liu et al. developed a denoising method combining EMD and interval threshold methods to suppress the random error of FOG signal successfully [18]. Wang et al. presented a new method which combined EMD and recursive least square (RLS) method to reduce the noise of FOG signal. Firstly, EMD is used to decompose the noisy signal, then RLS algorithm is used to filter the selected IMFs [19]. Different from the WT method, the EMD algorithm can decompose different data adaptively, but the disadvantages of mode mixing and end effect can not be ignored.
Variational mode decomposition (VMD) is a new signal decomposition method proposed in recent years [20]. It overcomes the problem of mode mixing in the EMD method and can deal with non recursive signals better. Wu et al. employed VMD to decompose the FOG signal, and then the decomposed signal was filtered based on the generalized morphological filter (CGMF). They proved the feasibility of the proposed scheme through experiments [21]. Zhang et al. first decomposed the gyroscope signal by utilizing VMD, then filtered the component signal by the adaptive Sigmoid function-based tracking differentiator (STD), and finally realized an efficient denoising [22]. Although VMD algorithm performs well in signal decomposition, its performance is affected by two parameters, namely decomposition mode number K and quadratic penalty factor α . A traditional trial and error method that depends on experience will lead to an unreliable result and will waste a lot of time. It will greatly limits the performance of the VMD algorithm. Besides, the selection of relevant modes is also a problem.
In order to solve the difficulties of the VMD algorithm in processing signal noise, a new method based on the beetle swarm antenna search (BSAS) algorithm and the VMD algorithm is proposed in this paper. The other five section of this paper are arranged as follows. In Section 2, the BSAS algorithm, which is a new intelligent optimization algorithm, is described in detail. In addition, VMD algorithm will also be explained in this section. The denoising procedures for the FOG signal based on the proposed method will be introduced in Section 3. In Section 4, through a simulation study, the performance of the proposed scheme is compared with the other three common schemes, and the superiority of the proposed scheme is confirmed by analyzing the simulation results. Next, in Section 5, static and dynamic rotation test experiments based on FOG further prove the feasibility and effectiveness of the proposed scheme. Finally, in the last section, the conclusion of this paper is presented, and the future work is prospected.

2. Theoretical Background

2.1. Beetle Swarm Antenna Search Algorithm

The beetle swarm antenna search (BSAS) algorithm is a new intelligent optimization algorithm which was proposed in 2018 [23,24]. It is a variant algorithm based on a beetle antennae search (BAS) algorithm which was proposed in 2017 [25]. The BAS algorithm optimizes the optimization index by imitating the predatory behavior of the beetle. When the beetle is preying, the two antennae on the left and right sides of its head will receive food odour from two directions respectively. According to the odour intensity in two directions, the beetle will choose to approach the direction with stronger odour intensity.
The BSAS algorithm is an improved algorithm of BAS. It expands the number of beetles from a single one in the BAS algorithm to a beetle group, which greatly improves the optimization ability and stability of the algorithm. Besides, compared with the BAS algorithm, the BSAS algorithm adds some additional parameters and improves the step size update strategy, which greatly improves the ability to solve complex problems.
The intensity of food odor can be expressed by the function to be optimized, (i.e., the fitness function), and then the biological behavior expressed by BSAS algorithm can be expressed by mathematical formula. The search direction of each beetle in the beetle group is randomly initialized at each search, as shown in Equation (1).
d i r i = r n d n , 1 i r n d n , 1 i 2 ,
where i is the currently selected beetle and n is the dimension of the variable to be optimized. r n d ( · )  denotes random function which can generate a random direction d i r . The two antennae used to receive odor intensity can be expressed as follows:
x r i k = x i k + d k d i r i x l i k = x i k d k d i r i
In Equation (2), k represents the current iteration times, and x l i and x r i represent the left and right antennae of the ith beetle. d k is sensing distance which denotes the length between the antennae and the head of beetle. d k can be expressed as Equation (3), where d 0 is the minimum resolution of sensing distance, it needs to be set in advance. c d is the decay coefficient of sensing distance.
d k = c d d k 1 + d 0 c d ( 0 , 1 ) .
Based on the concept of beetle group in BSAS algorithm, the position x i k of each beetle in the beetle group can be obtained as shown in Equation (4)
x i k = x i k 1 + s k · d i r i · s i g n f x r i k 1 f x l i k 1 ,
where f ( · ) represents the odor intensity (i.e., the fitness function). s i g n is symbolic function, which can choose a more ideal optimization direction for the next iteration. s k is search step size, as shown in Equation (5), s 0 is the minimum resolution of the search step size, and c s is the decay coefficient of the search step size.
s k = c s s k 1 + s 0 c s ( 0 , 1 ) .
In addition, compared with the BAS algorithm, two probability constants p s t and p p o s , which are between 0 and 1, are added to the BSAS algorithm. Where p s t is used to control the step size update. It means that when no better position is found, the probability that the step size will be updated in the next iteration is p s t , and the probability that it will not be updated is ( 1 p s t ) . By setting the p s t parameter, the search step size will no longer be unconditionally updated and reduced, thus it can improve the ability of global optimization. Moreover, in order to ensure the efficiency of the algorithm and avoid the bad results caused by improper p s t parameter setting, the BSAS algorithm proposes the n s t parameter. n s t represents the maximum number of invalid searches by using the same step size. When the cumulative number exceeds n s t , in the next iteration, the search step size will be forced to update.
Similar to p s t , p p o s is used to control beetle position update. It means that in each iteration, if M beetles in the beetle group find a better position, then in the next iteration, the probability of beetle being updated to the best position is p p o s . On the contrary, the probability that the beetle will be updated to a better position rather than the best position is ( 1 p p o s ) . By setting p p o s parameter, the disturbance of local optimal value can be effectively reduced, and the ability of global optimization can be improved.
In summary, combined with the mathematical expression and biological characteristics of BSAS algorithm, the pseudocode of BSAS algorithm is shown in Algorithm 1.
Algorithm 1 BSAS algorithm
Input: Define the fitness function f ( · ) . Set the optimization variable x to be optimized and determine its dimension n.
Output: x b e s t and f b e s t .
 1: Initialize:
   Initialize the position of beetle group x i 0 and probability constant p p o s ;
   Initialize the search step size s 0 and decay coefficient c s ;
   Initialize probability constant p s t and set parameter n s t ;
   Initialize the sensing distance d 0 and decay coefficient c d ;
   Initialize the maximum number of iterations k m a x ;
   Set initial optimization results for x b e s t and f b e s t .
 2: if k < = k m a x then
 3:    Generate random directions d i r i for each beetle by Equation (1).
 4:    Calculate the antennae position x l i k and x r i k of each beetle by Equation (2).
 5:    Calculate the position of each beetle x i k by Equation (4), and calculate the fitness function value f ( x i k ) .
 6:    Compare all f ( x i k ) and f b e s t in this iteration.
 7:    if f ( x i k ) < f b e s t then
 8:        if r a n d ( 1 ) = a < p p o s then
 9:            x b e s t = a r g m i n ( f ( x i k ) )
 10:            f b e s t = m i n ( f ( x i k ) )
 11:        else
 12:            x b e s t = s a m p l e ( x i k )
 13:            f b e s t = ( f ( x s a m p l e k ) )
 14:        end if
 15:    else
 16:        if r a n d ( 1 ) = b < p s t | | i > n s t then
 17:           Search step s k is updated by Equation (5).
 18:           Sensing distance d k is updated by Equation (3).
 19:        else
 20:            i = i + 1
 21:        end if
 22:    end if
 23:     k = k + 1
 24:end if
Compared with the BAS algorithm, the BSAS algorithm has greatly improved its optimization and stability ability by increasing the number of beetles. Furthermore, the BSAS algorithm introduces additional parameters in the step size strategy and position update strategy, which makes it more likely to find the global optimal value. Because of the advantages mentioned above, the BSAS algorithm is widely used in many fields [26,27,28].

2.2. Variational Mode Decomposition Algorithm

Dragomiretskiy et al. proposed a variational mode decomposition algorithm (VMD) in 2014, which can decompose the signal into K band-limited intrinsic mode functions (BLIMFs) [29]. The center frequency of each modal component is ω k .
Firstly, the estimation of the center frequency and bandwidth for BLIMF components is transformed into a variational problem with constraints by using VMD algorithm. Then, by introducing penalty factor and Lagrange function, the problem with constraint condition is transformed into the problem without constraint condition. Finally, the alternative multiplier algorithm is used to get the optimal solution. The VMD algorithm can be regarded as an improved algorithm of Wiener filter [30]. It effectively solves the problem of modal aliasing in the EMD algorithm [31]. The specific VMD algorithm can be expressed as follows:
To begin with, the signal x ( t ) is decomposed into K BLIMF components u k ( t ) , k = 1 , 2 , , K , as shown in Equations (6) and (7).
x ( t ) = i = 1 K u k ( t )
u k ( t ) = A k ( t ) cos φ k ( t ) .
The instantaneous frequency ω k ( t ) of u k ( t ) is presented as Equation (8):
ω k ( t ) = d φ k ( t ) d t .
After that, the Hilbert transform is applied to each mode component, and the frequency spectrum of each BLIMF component is modulated to the frequency band with the center frequency ω k , as shown in Equation (9):
δ ( t ) + j π t * u k ( t ) · e j ω k t ,
where δ ( t ) is the Dirac distribution, and ω k is the center frequency. * represents symbol of convolution operation. Next, the Gaussian smoothness of the transformed signal can be used to estimate the bandwidth of each component. The variational problem with constraints is given as Equation (10):
min u k , ω k k = 1 K t δ ( t ) + j π t * u k ( t ) · e j ω k t 2 2 s . t . k = 1 K u k = x ( t ) .
By introducing quadratic penalty factor α and Lagrange multiplication operator λ ( t ) , the constrained variational problem can be transformed into unconstrained variational problem. It should be noted that the augmented Lagrangian is expressed as Equation (11)
L u k , ω k , λ = α k = 1 K t δ ( t ) + j π t * u k ( t ) · e j ω k t 2 2 + x ( t ) k = 1 K u k ( t ) 2 2 + λ ( t ) , x ( t ) k = 1 K u k ( t ) ,
u k n + 1 ( t ) = argmin u k X α t δ ( t ) + j π t * u k ( t ) · e j ω k t 2 2 + x ( t ) i u i ( t ) + λ ( t ) 2 2 2 .
Alternate direction method of multipliers (ADMM) is used to solve the unconstrained variational problem, and the saddle point is obtained by updating u k n + 1 , ω k n + 1 , and λ n + 1 in frequency respectively. The specific formula is denoted as Equations (13)–(15). Where component u k n + 1 ( t ) can be expressed as Equation (12).
u ^ k n + 1 ( ω ) = x ^ ( ω ) i 1 u ^ i ( ω ) + λ ^ ( ω ) 2 · 1 1 + 2 α ω ω k 2
ω k n + 1 = 0 ω u ^ k ( ω ) 2 d ω 0 u ^ k ( ω ) 2 d ω
λ ^ n + 1 ( ω ) = λ ^ n ( ω ) + τ x ^ ( ω ) k u ^ k n + 1 ( ω ) ,
where represents the Fourier transform, n is the number of iterations, and τ denotes the parameter of noise tolerance. Finally, according to Equations (13)–(15), u k n + 1 , ω k n + 1 , and λ n + 1 are continuously updated until the iteration termination condition is met. The iteration termination condition is shown in Equation (16):
k u ^ k n + 1 u ^ k n 2 2 u ^ k n 2 2 < e .
In conclusion, the detailed steps of the VMD algorithm are shown in Algorithm 2.
Algorithm 2 VMD algorithm
1:
Initialize u ^ k , w k , λ ^ and n.
2:
Update u ^ k n + 1 ( ω ) by Equation (13).
3:
Update ω k n + 1 by Equation (14).
4:
Update λ ^ n + 1 ( ω ) by Equation (15).
5:
Repeat the above updating steps until the iteration termination condition Equation (16) is met.

3. Methodology

This section will first introduce the concept of VMD parameter optimization. Then the fitness function based on permutation entropy is explained. Next, a new method based on Hausdorff distance (HD) between the probability density function (PDF) of all BLIMFs and that of the original signal is proposed to determine the relevant modes. Finally, the detailed steps of the proposed method are given.

3.1. VMD Parameter Optimization

The VMD algorithm is not a recursive signal decomposition method, and it can adaptively decompose the signal in frequency domain and effectively get the BLIMFs of the signal. Compared with the EMD algorithm, the VMD algorithm has a clearer mathematical theory and better noise robustness [32]. Because of the advantages of the VMD algorithm, the VMD algorithm is widely used in many fields [33,34,35].
However, some studies have proved that the data processing results by using the VMD method are affected by quadratic penalty factors α and the number of decomposition mode K [36]. As for the selection of two parameters, the traditional selection method is depending on experience [29]. For the sake of better determining the combination of parameters, this paper proposes a parameter optimization method based on the BSAS algorithm and permutation entropy.

3.2. Fitness Function Based on Permutation Entropy

In order to optimize the two parameters K and α of VMD by using BSAS algorithm, the optimization index (i.e., the fitness function) should be determined first.
Permutation entropy (PE) is proposed by Bandt et al., and it is a method to judge the randomness and complexity of time series [37]. It has small calculation and strong anti-interference ability, and it is applied in many fields as fitness function or feature selection [38,39]. The specific introduction of PE is as follows:
The matrix shown in Equation (17) is the phase space reconstruction of time series { x ( i ) , i = 1 , 2 , , n } .
x ( 1 ) x ( 1 + τ ) x ( 1 + ( m 1 ) τ ) x ( 2 ) x ( 2 + τ ) x ( 2 + ( m 1 ) τ ) x ( j ) x ( j + τ ) x ( j + ( m 1 ) τ ) x ( K ) x ( K + τ ) x ( K + ( m 1 ) τ ) ,
where m and τ are embedding dimension and delay time respectively, and j = 1 , 2 , 3 , , K . Each row in the matrix represents a reconstruction component, so there are K reconstruction components in total, K = n ( m 1 ) τ .
The jth component ( x ( j ) , x ( j + τ ) , , x ( j + ( m 1 ) τ ) in the reconstruction matrix is rearranged in ascending order, as shown in Equation (18). j 1 , j 2 , , j m represents the index number of each element in the reconstruction component.
x i + j 1 1 τ x i + j 2 1 τ x i + j m 1 τ .
If there are two identical reconstruction components, for example:
x i j p 1 τ = x i j q 1 τ .
Then it should be sorted by the size of the values of j p and j q . Therefore, for any vector x ( i ) composed of time series data, a set of symbol sequences can be obtained:
S ( l ) = j 1 , j 2 , , j m .
It should be noted that m dimensional phase space has m ! different sequence of symbols, and S ( l ) is one of them. If the probability of each symbol sequence is calculated as P 1 , P 2 , , P k , then the permutation entropy of time sequence x ( i ) can be defined as Equation (21):
H P = j = 1 k P j ln P j .
The value of H P indicates the randomness of the time series x ( i ) . The smaller the value of H P is, the more orderly it is; otherwise, the more random it is. It means that after the original signal is decomposed by VMD algorithm, the less noise the BLIMF component contains, the less PE it has.
VMD algorithm is used to decompose the original signal in each iteration. The permutation entropy of each BLIMF component is calculated, where the minimum value is taken as the optimal result of the current iteration. In the optimization process, the fitness function is given as follows:
f ( K , α ) = min H p ( B L I M F s ) .

3.3. Selecting Relevant Modes

According to the principle of the VMD algorithm mentioned above, the original signal can be decomposed into several BLIMF components from low to high frequency by using the VMD algorithm. Generally speaking, the random noise is mainly concentrated in the high frequency component. In order to reduce the noise of the original signal better by using the VMD algorithm, it is necessary to find the demarcation point between the noise BLIMF component and the noiseless BLIMF component, which also means to determine the relevant modes.
In traditional VMD noise reduction method, the demarcation point is determined by calculating the continuous mean square error (CMSE) of each BLIMFs. Where the BLIMF with the first minimum value of CMSE is the demarcation point k th [40]. The calculation formula of CMSE and the definition of demarcation point are shown in Equations (23) and (24). Where u k ( t ) is the component obtained by using VMD algorithm to decompose the original signal, and x ^ k is the reconstructed signal.
CMSE x ^ k , x ^ k + 1 = 1 N i = 1 N x ^ k ( t i ) x ^ k + 1 ( t i ) 2 = 1 N i = 1 N u k ( t i ) 2 ,
k th = argmin 1 k N 1 CMSE x ^ k , x ^ k + 1 .
However, the effect of this method is often not ideal, and it greatly reduces the noise reduction effect of the VMD algorithm [13,20]. In this paper, the Hausdorff distance (HD) between the probability density function (PDF) of all BLIMFs and that of the signal to be decomposed is proposed to determine the demarcation point.
HD was proposed by Huttenlocher in 1993 [41]. It is often used to calculate the similarity of two different sample sets. In this paper, HD is utilized to calculate the similarity between each BLIMFs and the original signal.
According to the principle of HD, the HD of two sample sets A= a 1 , , a i and B= b 1 , , b j are as follows:
HD ( A , B ) = max { h d ( A , B ) , h d ( B , A ) } .
Equation (25) indicates bidirectional Hausdorff distance between A = a 1 , , a i and B = b 1 , , b j , and it is the basic formula of HD. Where h d ( A , B ) and h d ( B , A ) are called unidirectional Hausdorff distance of sets A to B and B to A, respectively, as shown below:
h d ( A , B ) = max a A min b B a b
h d ( B , A ) = max b B min a A b a .
In Equations (26) and (27), a and b represent elements in sets A and B respectively, and · refers to l 2 -norm between two elements.
The probability density function (PDF) can represent the distribution of data, and it can also be used as a similarity measure for different data [42]. Therefore it can be combined with HD to calculate the similarity for signals, as follows:
S ( i ) = HD P D F ( x ( t ) ) , P D F B L I M F i
As shown in Equation (28), where S ( i ) is HD based PDF between each BLIMFs and the original signal x ( t ) .
To sum up, the selection of relevant modes (i.e., components determined by the demarcation point), can be solved by Equation (29):
k th = arg max 1 i K { S ( i + 1 ) S ( i ) } ,
where k th is the demarcation point, and it also represents the maximum slope based on HD of the PDF for two adjacent BLIMF components [43]. i represents the ith BLIMF, and there are K in total.

3.4. Proposed Methodology

According to the above analysis and statement, the detailed steps of the proposed method can be described as follows:
Step 1:
The parameters of the BSAS algorithm are initialized. At the same time, the parameters in permutation entropy are initialized too.
Step 2:
Firstly, the VMD algorithm is used to decompose the original signal. After that, the fitness function value based on permutation entropy is calculated, and the parameters K and α are optimized by using the BSAS algorithm.
Step 3:
Determine whether the termination condition is met. If it is, optimal parameter combination ( K , α ) is saved, otherwise step 2 is repeated.
Step 4:
The original signal is decomposed by using the VMD algorithm based on optimal parameter combination.
Step 5:
The HD between the PDF of all BLIMF and that of the original signal is calculated to determine the demarcation point.
Step 6:
According to the demarcation point, the selected BLIMF components are retained as the relevant modes, and the unselected BLIMF components are removed.
Step 7:
The relevant modes are reconstructed, and finally the denoised signal is obtained.
Based on the above introduction, the flowchart of the proposed scheme is shown in Figure 1. It should be pointed out that the algorithm proposed in this paper mainly includes two parts, one is the parameter optimization of the VMD algorithm, the other is the decomposition and reconstruction of signals. Due to the large amount of calculation in the process of parameter optimization, this part usually needs to be processed offline. However, the proposed method is still superior to the traditional trial and error method in terms of result reliability and time consumption. Therefore, this paper will mainly study the effect of noise reduction, and not do too much analysis on the timeliness and delay analysis of the algorithm.

4. Simulation and Analysis

In this section, the performance of the proposed scheme is demonstrated by simulation. Firstly, the simulation environment is introduced, which including initialization of algorithm parameters and selection of simulation signals. After that, the simulation results are analyzed and compared with the results of several common denoising schemes.

4.1. Simulation Environment

Bumps signal is a kind of continuously changing signal, which has a wide range of frequency coverage. It is similar to the characteristics of the FOG signal. The original Bumps signal is shown in Figure 2.
In practice, the FOG signal is often accompanied by some small vibration or bumpy signal and a large number of random noise in actual use. So in order to better verify the performance of the proposed algorithm, the impulse signal, as shown in Figure 3, is added to the Bumps signal. In addition, Gaussian white noise is also added to the whole signal, where SNR = 10 dB and signal length N = 1024. The final simulation signal is as shown in Figure 4.
Besides, some parameters of the proposed algorithm need to be set before the simulation. To begin with, the parameters of the BSAS algorithm are initialized as follows—the number of beetles is 10 and the maximum number of iterations is 50; The initial value of step size and sensing distance are 5 and 4 respectively; The decay coefficients of step and sensing distance are both 0.95. The probability constants of position updating and step updating are both set to 0.85, and the maximum number of invalid search is 2. Then, according to experience, the embedding dimension of permutation entropy is set to 6, and the delay time is initialized to 1.

4.2. Simulation Results and Analysis

According to the steps of the proposed algorithm, firstly, two parameters K and α of VMD are optimized by the BSAS algorithm. The convergence curve based on the BSAS algorithm is shown in Figure 5. In the first 10 iterations, the speed of curve convergence is the fastest. After that, the convergence tends to be stable.
It can be seen from the optimization process that the minimum value of entropy permutation is 0.3338, and the corresponding K and α are 7 and 1530 respectively.
After that, the simulation signal is decomposed by VMD algorithm based on the optimal parameter combination which obtained by BSAS algorithm, and all BLIMFs are as shown in the Figure 6:
Then, according to step 5 of the proposed algorithm, the Hausdorff distance between the PDF of all BLIMFs and that of the original signal should be calculated to determine the demarcation point.
As shown in Figure 7a, in the light of the definition of demarcation point proposed by Equations (28) and (29), on the one hand, BLIMF1 and BLIMF2 are reserved as relevant modes for signal reconstruction. On the other hand, BLIMF3 to BLIMF7 are removed as unselected modes. Figure 7b shows the result of traditional scheme based on CMSE.
Finally, the selected components are reconstructed to get the denoised signal. In order to evaluate the performance of the proposed algorithm, two evaluation indexes, signal-to-noise ratio (SNR) and root mean square error (RMSE), are given as Equations (30) and (31).
SNR = 10 · log i = 1 N | x ( i ) | 2 i = 1 N | x ^ ( i ) x ( i ) | 2 ,
RMSE = i = 1 N | x ^ ( i ) x m o d e l ( i ) | 2 N ,
where N is the length of the signal. x is the noiseless signal, x ^ is the denoised signal, and x m o d e l represents ideal value.
Figure 8a shows the noise reduction results based on the proposed scheme. At the same time, three other common denoising methods, including the traditional VMD scheme, the traditional EMD scheme and the scheme based on wavelet transformation (WT), are also used to process the simulation signals. As shown in Figure 8, the noise reduction effect is worst by using the traditional EMD method, and it is obvious that a lot of noise exists in its results. The SNR and RMSE of this method are 15.2305 and 0.3116 respectively. The noise reduction effect based on WT method and traditional VMD method is satisfactory, but both of them have some small noise signals. The SNR and RMSE of the two methods are given in Table 1. The best dnoised effect is obtained by the proposed scheme, which has the smoothest and most stable noise reduction result curve. Among all methods, the proposed method has the largest SNR and the smallest RMSE, which are 18.3232 and 0.2183 respectively.
Based on the Figure 8 and two evaluation indexes shown in Table 1, it can be concluded that after processing the original simulation signal with the proposed scheme, the denoised signal has the highest accuracy and the best stability. Therefore, the proposed scheme is superior to the other three schemes in noise reduction performance.

5. Experimental Analysis

In this section, static test and dynamic rotation test experiments based on FOG are carried out respectively to prove the feasibility of the proposed method. The equipment which will be used in experiments is a high precision three-axis turntable, as shown in Figure 9b. In the experiment, the temperature of room was kept at 26 degrees centigrade.
It is worth noting that the FOG-based inertial measurement unit (IMU) used in two experiments is independently designed and produced by Harbin Engineering University. The parameters of the FOG are given in Table 2.

5.1. Static Test Experiment

In the static test experiment, the FOG is placed on a three-axis turntable firstly. Then, the signal of measurement axis, which is aligned with the upward direction, is collected. The sampling frequency is 100 Hz. In order to avoid the measurement error of the instrument in the start-up stage, 12,000 sampling points are selected in the middle stage of the whole measurement process, as shown in Figure 10.
The raw signals are denoised by the proposed scheme, and the other three schemes are also implemented separately. The results of the four methods are shown in Figure 11.
It can be seen from Figure 11a that all four methods can effectively reduce the noise of FOG static signal. Figure 11b is a local enlarged graphs, which shows part of the noise reduction result. The scheme proposed in this paper is represented by red curve, which has the best effect and stability in noise reduction performance. The noise reduction effect of the WT method and the traditional VMD method are the second and third respectively, and their results are similar. Although the traditional EMD method, represented by green curve, can also reduce the noise of the original signal, its effect is the worst compared with the other three methods.
In order to better compare the performance of different methods, in addition to the RMSE, the noise intensity (NI) [44] is proposed as an evaluation index for all experiments. The definition of NI is similar to the standard deviation, which can be expressed as follows:
NI = i = 1 N x ^ i x ¯ 2 N .
Table 3 records the noise intensity and root mean square error for four method. The minimum value of the root mean square error comes from the proposed method, which is 1.2939 × 10 4 /s. For the other three methods, the root mean square error is 2.2630 × 10 4 /s, 2.8105 × 10 4 /s and 1.6716 × 10 4 /s respectively. In terms of the noise intensity, the proposed scheme is 7.3872 × 10 5 /s, which is also the minimum in the four methods. The maximum noise intensity obtained by traditional EMD method is 2.6109 × 10 4 /s. The noise intensity of traditional VMD method and WT method are 1.9982 × 10 4 /s and 1.2901 × 10 4 /s respectively.
In summary, the static test experiment demonstrate that the proposed scheme can effectively reduce the noise of the FOG static signal. In terms of noise reduction effect, the proposed scheme is better than the other three schemes.

5.2. Dynamic Rotation Test Experiment

In order to further verify the feasibility of the proposed scheme, a dynamic rotation test experiment is designed in this paper. First of all, the FOG-based IMU is placed on the three-axis turntable, and its measuring axis is aligned with the upward direction. After that, the FOG is rotated around the measuring axis by using the three-axis turntable, and the rotation speed is 5 /s and 10 /s respectively. As the same as the static test, 12,000 sampling points are selected for processing.
Figure 12 and Figure 13 show the noise reduction results of the four methods respectively at rotation speeds of 5 /s and 10 /s. Two groups of experimental results show that the noise reduction effect of the proposed scheme is significantly better than the other three schemes at different rotation speeds. From the the partial enlarged graphs, Figure 12b and Figure 13b, it can be seen that the results of the proposed scheme are smoother, and has the least signal burr. More information about the performance evaluation of the four schemes is shown in Table 4.
Table 4 gives the noise reduction results of four methods at different rotation speeds 5 /s and 10 /s. It can be seen from the table that the NI of the proposed method is 1.2759 × 10 4 /s and 1.3437 × 10 4 /s respectively, and RMSE is 1.4448 × 10 4 /s and 1.8412 × 10 4 /s, which are all the minimum values compared with other methods. The results of traditional EMD method are the worst. NI and RMSE are 4.8796 × 10 4 /s and 4.9211 × 10 4 /s respectively in 5 /s; 4.9692 × 10 4 /s and 5.1202 × 10 4 /s respectively in 10 /s.
From the processing results of the original signal at different rotation speeds above, it can be concluded that in different test environments, the proposed scheme in this paper always has better noise reduction performance.

6. Conclusions

This paper presents a new method for the denoising of a FOG signal. This method is based on the BSAS algorithm and permutation entropy to optimize two parameters K and α of the VMD algorithm firstly. Where the fitness function is designed by calculating the permutation entropy of each BLIMF. After that, in order to better find the relevant modes, the Hausdorff distance between the probability density function of all BLIMFs and that of the original signal is proposed in this paper to determine the demarcation point. Finally, the selected modes are reconstructed into denoised signal.
For the sake of demonstrating the performance of proposed method, simulations and experiments are performed in this paper. The simulation results show that the proposed scheme is superior to the traditional VMD method, the traditional EMD method and the WT method in terms of noise reduction effect and stability. In addition, the feasibility of proposed scheme is further proved by static test experiment and dynamic rotation test experiment.
However, the method proposed in this paper is still in the stage of development and exploration. Due to the large amount of calculation in the process of optimization, it cannot be well adapted to online application. Therefore, the future work will focus on the improvement of the new algorithm to achieve online noise reduction.

Author Contributions

Conceptualization and methodology, P.W. and Y.G.; validation, P.W.; data curation, P.W. and F.Z.; formal analysis, F.Z. and M.W.; writing, review and editing, P.W. and C.Q.; funding acquisition, G.L. and Y.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research work is supported by National Natural Science Foundation of China (NSFC) (No. 61803118), the Science and Technology Research Program of Chongqing Municipal Education Commission (KJZD-K201804701).

Acknowledgments

The authors would like to thank everyone who has helped us with this article. It is your opinion that makes this article better.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Eisele, J.; Song, Z.; Nelson, K.; Mohseni, K. Visual-Inertial Guidance with a Plenoptic Camera for Autonomous Underwater Vehicles. IEEE Robot. Autom. Lett. 2019, 43, 2777–2784. [Google Scholar] [CrossRef]
  2. Carlone, L.; Karaman, S. Attention and Anticipation in Fast Visual-Inertial Navigation. IEEE Trans. Robot. 2019, 35, 1–20. [Google Scholar] [CrossRef] [Green Version]
  3. Gross, J.N.; Gu, Y.; Rhudy, M.B. Robust UAV relative navigation with DGPS, INS, and peer-to-peer radio ranging. IEEE Trans. Autom. Sci. Eng. 2015, 12, 935–944. [Google Scholar] [CrossRef]
  4. Chung, H.; Ojeda, L.; Borenstein, J. Accurate mobile robot dead-reckoning with a precision-calibrated fiber-optic gyroscope. IEEE Trans. Robot. Autom. 2001, 17, 80–84. [Google Scholar] [CrossRef]
  5. Ren, W.; Luo, Y.; He, Q.N.; Zhou, X.; Ren, G. Stabilization Control of Electro-Optical Tracking System With Fiber-Optic Gyroscope Based on Modified Smith Predictor Control Scheme. IEEE Sensors J. 2018, 18, 8172–8178. [Google Scholar] [CrossRef]
  6. Feng, S.; Chao, L.; Wei, G.; Qi, N. Research on modeling and compensation method of Fiber Optic Gyro’Random error. In Proceedings of the IEEE International Conference Mechatronics and Automation, Niagara Falls, ON, Canada, 29 July–1 August 2005. [Google Scholar]
  7. Wang, L.; Zhang, C.; Lin, T.; Li, X.; Wang, T. Characterization of a Fiber Optic Gyroscope in a Measurement While Drilling System with the Dynamic Allan Variance. Measurement 2015, 75, 263–272. [Google Scholar] [CrossRef]
  8. Narasimhappa, M.; Nayak, J.; Terra, M.H.; Sabat, S.L. ARMA model based adaptive unscented fading Kalman filter for reducing drift of fiber optic gyroscope. Sens. Actuator A-Phys. 2016, 251, 42–51. [Google Scholar] [CrossRef]
  9. Narasimhappa, M. Modeling of Inertial Rate Sensor Errors Using Autoregressive and Moving Average (ARMA) Models. In Gyroscopes-Principles and Applications; IntechOpen: London, UK, 2019. [Google Scholar]
  10. Nassar, S.; Schwarz, K.P.; El-Sheimy, N.; Noureldin, A. Modeling Inertial Sensor Errors Using Autoregressive (AR) Models. Navigation 2004, 51, 259–268. [Google Scholar] [CrossRef]
  11. Nassar, S.; Schwarz, K.; Elsheimy, N. INS and INS/GPS Accuracy Improvement Using Autoregressive (AR) Modeling of INS Sensor Errors. In Proceedings of the Ion National Technical Meeting, San Diego, CA, USA, 26–28 January 2004; pp. 936–944. [Google Scholar]
  12. Narasimhappa, M.; Sabat, S.L.; Nayak, J. Fiber-Optic Gyroscope Signal Denoising Using an Adaptive Robust Kalman Filter. IEEE Sens. J. 16, 3711–3718. [CrossRef]
  13. He, J.; Sun, C.; Wang, P. Noise Reduction for MEMS Gyroscope Signal: A Novel Method Combining ACMP with Adaptive Multiscale SG Filter Based on AMA. Sensors 2019, 19, 4382. [Google Scholar] [CrossRef] [Green Version]
  14. Daubechies, I. The wavelet transform, time-frequency localization and signal analysis. IEEE Trans. Inf. Theory 1990, 36, 961–1005. [Google Scholar] [CrossRef] [Green Version]
  15. Flandrin, P.; Rilling, G.; Goncalves, P. Empirical mode decomposition as a filter bank. IEEE Signal Process. Lett. 2004, 11, 112–114. [Google Scholar] [CrossRef] [Green Version]
  16. Zhang, Q.; Wang, L.; Gao, P.; Liu, Z. An Innovative Wavelet Threshold Denoising Method for Environmental Drift of Fiber Optic Gyro. Math. Probl. Eng. 2016, 2016, 1–8. [Google Scholar] [CrossRef] [Green Version]
  17. Ma, J.; Yang, Z.; Shi, Z.; Zhang, X.; Liu, C. Application and Optimization of Wavelet Transform Filter for North-Seeking Gyroscope Sensor Exposed to Vibration. Sensors 2019, 19, 3624. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Liu, C.; Yang, Z.; Shi, Z.; Ma, J.; Cao, J. A Gyroscope Signal Denoising Method Based on Empirical Mode Decomposition and Signal Reconstruction. Sensors 2019, 19, 5064. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Wang, D.; Xu, X.; Zhang, T.; Zhu, Y.; Tong, J. An EMD-MRLS de-noising method for fiber optic gyro Signal. Optik 2019, 183, 971–987. [Google Scholar] [CrossRef]
  20. Wang, Y.; Markert, R. Filter bank property of variational mode decomposition and its applications. Signal Process. Process. 2016, 120, 509–521. [Google Scholar] [CrossRef]
  21. Wu, Y.; Shen, C.; Cao, H.; Che, X. Improved Morphological Filter Based on Variational Mode Decomposition for MEMS Gyroscope De-Noising. Micromachines 2018, 9, 246. [Google Scholar] [CrossRef] [Green Version]
  22. Zhang, X.; Cao, H.; Wang, C.; Kou, Z.; Shao, X.; Li, J.; Liu, J.; Shen, C. Delay-Free Tracking Differentiator Design Based on Variational Mode Decomposition: Application on MEMS Gyroscope Denoising. J. Sensors 2019, 2019, 1–13. [Google Scholar] [CrossRef]
  23. Wang, J.; Chen, H.; Yuan, Y.; Huang, Y. A novel efficient optimization algorithm for parameter estimation of building thermal dynamic models. Build. Environ. 2019, 153, 233–240. [Google Scholar] [CrossRef]
  24. Wang, J.; Chen, H. BSAS: Beetle Swarm Antennae Search Algorithm for Optimization Problems. arXiv 2018, arXiv:1807.10470. [Google Scholar]
  25. Jiang, X.; Li, S. BAS: Beetle antennae search algorithm for optimization problems. arXiv 2017, arXiv:1710.10724. [Google Scholar] [CrossRef]
  26. Xie, S.; Garofano, V.; Chu, X.; Negenborn, R.R. Model predictive ship collision avoidance based on Q-learning beetle swarm antenna search and neural networks. Ocean. Eng. 2019, 193, 106609. [Google Scholar] [CrossRef]
  27. Lin, X.; Liu, Y.; Wang, Y. Design and Research of DC Motor Speed Control System Based on Improved BAS. In Proceedings of the 2018 Chinese Automation Congress (CAC), Xi’an, China, 30 November–2 December 2018; pp. 3701–3705. [Google Scholar]
  28. Wang, P.; Gao, Y.; Wu, M. In-Field Calibration of Triaxial Accelerometer Based on Beetle Swarm Antenna Search Algorithm. Sensors 2020, 20, 947. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Dragomiretskiy, K.; Zosso, D. Variational Mode Decomposition. IEEE Trans. Signal Process. 2014, 62, 531–544. [Google Scholar] [CrossRef]
  30. Norman, L. The Wiener (Root Mean Square) Error Criterion in Filter Design and Prediction. J. Math. Phys. 1946, 25, 261–278. [Google Scholar]
  31. Wu, Z.; Huang, N. Ensemble Empirical Mode Decomposition: A Noise-Assisted Data Analysis Method. Adv. Adapt. Data Anal. 2009, 1, 1–41. [Google Scholar] [CrossRef]
  32. Lahmiri, S. Comparative study of ECG signal denoising by wavelet thresholding in empirical and variational mode decomposition domains. Healthc. Technol. Lett. 2014, 1, 104–109. [Google Scholar] [CrossRef] [Green Version]
  33. Miao, Y.; Zhao, M.; Lin, J. Identification of mechanical compound-fault based on the improved parameter-adaptive variational mode decomposition. ISA Trans. 2019, 84, 82–95. [Google Scholar] [CrossRef]
  34. Yang, K.; Wang, G.; Dong, Y.; Zhang, Q.; Sang, L. Early chatter identification based on an optimized variational mode decomposition. Mech. Syst. Signal Process. 2019, 115, 238–254. [Google Scholar] [CrossRef]
  35. Chen, X.; Yang, Y.; Cui, Z.; Shen, J. Vibration fault diagnosis of wind turbines based on variational mode decomposition and energy entropy. Energy 2019, 174, 1100–1109. [Google Scholar] [CrossRef]
  36. Mohanty, S.; Gupta, K.K.; Raju, K.S. Bearing fault analysis using variational mode decomposition. In Proceedings of the 9th international conference on industrial and information, Gwalior, India, 15–17 December 2014. [Google Scholar]
  37. Bandt, C.; Pompe, B. Permutation Entropy: A Natural Complexity Measure for Time Series. Phys. Rev. Lett. 2002, 88, 174102. [Google Scholar] [CrossRef] [PubMed]
  38. Li, Y.; Li, Y.; Chen, X.; Yu, J. A novel feature extraction method for ship-radiated noise based on variational mode. decomposition and multi-scale permutation entropy. Entropy 2017, 19, 342. [Google Scholar]
  39. Zhou, J.; Guo, X.; Wang, Z. Research on Fault Extraction Method of Variational Mode Decomposition Based on Immunized Fruit Fly Optimization Algorithm. Entropy 2019, 21, 400. [Google Scholar] [CrossRef] [Green Version]
  40. Lu, Q.; Pang, L.; Huang, H.; Shen, C.; Cao, H.; Shi, Y.; Liu, J. High-G Calibration Denoising Method for High-G MEMS Accelerometer Based on EMD and Wavelet Threshold. Micromachines 2019, 10, 134. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Huttenlocher, D.; Klanderman, G.; Rucklidge, W. Comparing images using the Hausdorff distance. IEEE Trans. Pattern Anal. Mach. Intell. 1993, 15, 850–863. [Google Scholar] [CrossRef] [Green Version]
  42. Komaty, A.; Boudraa, A.; Augier, B.; Dareemzivat, D. EMD-Based Filtering Using Similarity Measure Between Probability Density Functions of IMFs. IEEE Trans. Instrum. Meas. 2014, 63, 27–34. [Google Scholar] [CrossRef] [Green Version]
  43. Chang, J.; Zhu, L.; Li, H.; Xu, F.; Liu, B.; Yang, Z. Noise reduction in Lidar signal using correlation-based EMD combined with soft thresholding and roughness penalty. Opt. Commun. 2018, 407, 290–295. [Google Scholar] [CrossRef]
  44. Li, G.; Yang, Z.; Yang, H. Noise reduction method of underwater acoustic signals based on uniform phase empirical mode decomposition, amplitude-aware permutation entropy, and Pearson correlation coefficie. Entropy 2018, 20, 918. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Flowchart of the proposed method.
Figure 1. Flowchart of the proposed method.
Entropy 22 00765 g001
Figure 2. Bumps original signal.
Figure 2. Bumps original signal.
Entropy 22 00765 g002
Figure 3. Impulse signal.
Figure 3. Impulse signal.
Entropy 22 00765 g003
Figure 4. Simulation signal.
Figure 4. Simulation signal.
Entropy 22 00765 g004
Figure 5. Convergence curve of the beetle swarm antenna search (BSAS) algorithm.
Figure 5. Convergence curve of the beetle swarm antenna search (BSAS) algorithm.
Entropy 22 00765 g005
Figure 6. Decomposition result of the simulation signal.
Figure 6. Decomposition result of the simulation signal.
Entropy 22 00765 g006
Figure 7. Determination of demarcation point.
Figure 7. Determination of demarcation point.
Entropy 22 00765 g007
Figure 8. Denoised results.
Figure 8. Denoised results.
Entropy 22 00765 g008
Figure 9. Experimental equipment. (a) FOG-based IMU. (b) Three-axis turntable.
Figure 9. Experimental equipment. (a) FOG-based IMU. (b) Three-axis turntable.
Entropy 22 00765 g009
Figure 10. Fiber optic gyroscope (FOG) static test raw signal.
Figure 10. Fiber optic gyroscope (FOG) static test raw signal.
Entropy 22 00765 g010
Figure 11. Denoising results of four different methods for static test experiment.
Figure 11. Denoising results of four different methods for static test experiment.
Entropy 22 00765 g011
Figure 12. Denoising results of four different methods at a rotational speed of 5 /s.
Figure 12. Denoising results of four different methods at a rotational speed of 5 /s.
Entropy 22 00765 g012
Figure 13. Denoising results of four different methods at a rotational speed of 10 /s.
Figure 13. Denoising results of four different methods at a rotational speed of 10 /s.
Entropy 22 00765 g013
Table 1. Noise reduction performance index of four methods.
Table 1. Noise reduction performance index of four methods.
MethodSignal-to-Noise Ratio (SNR/dB)Root Mean Square Error (RMSE)
Proposed method18.32320.2183
Traditional VMD17.28770.2459
Traditional EMD15.23050.3116
Wavelet transform17.36860.2436
Table 2. Parameters of the FOG.
Table 2. Parameters of the FOG.
Parameter ItemParameter Values
FOG dynamic range ( /s) 600
FOG bias stability ( /h)0.03
FOG random bias ( / h 1 / 2 )0.003
Table 3. Noise intensity (NI) and root mean square error (RMSE) comparison for the four methods.
Table 3. Noise intensity (NI) and root mean square error (RMSE) comparison for the four methods.
MethodNoise Intensity (NI)Root Mean Square Error (RMSE)
Proposed method7.3872 × 10 5 1.2939 × 10 4
Traditional VMD1.9982 × 10 4 2.2630 × 10 4
Traditional EMD2.6109 × 10 4 2.8105 × 10 4
Wavelet transform1.2901 × 10 4 1.6716 × 10 4
Table 4. Denoising results of four methods at different rotation speed.
Table 4. Denoising results of four methods at different rotation speed.
Rotational SpeedProposed MethodTraditional VMDTraditional EMDWavelet Transform
NIRMSENIRMSENIRMSENIRMSE
5 ( /s)1.2759 × 10 4 1.4448 × 10 4 3.2013 × 10 4 3.2723 × 10 4 4.8796 × 10 4 4.9211 × 10 4 2.3837 × 10 4 2.4785 × 10 4
10 ( /s)1.3437 × 10 4 1.8412 × 10 4 3.2248 × 10 4 3.4618 × 10 4 4.9692 × 10 4 5.1202 × 10 4 2.3852 × 10 4 2.6971 × 10 4

Share and Cite

MDPI and ACS Style

Wang, P.; Gao, Y.; Wu, M.; Zhang, F.; Li, G.; Qin, C. A Denoising Method for Fiber Optic Gyroscope Based on Variational Mode Decomposition and Beetle Swarm Antenna Search Algorithm. Entropy 2020, 22, 765. https://doi.org/10.3390/e22070765

AMA Style

Wang P, Gao Y, Wu M, Zhang F, Li G, Qin C. A Denoising Method for Fiber Optic Gyroscope Based on Variational Mode Decomposition and Beetle Swarm Antenna Search Algorithm. Entropy. 2020; 22(7):765. https://doi.org/10.3390/e22070765

Chicago/Turabian Style

Wang, Pengfei, Yanbin Gao, Menghao Wu, Fan Zhang, Guangchun Li, and Chao Qin. 2020. "A Denoising Method for Fiber Optic Gyroscope Based on Variational Mode Decomposition and Beetle Swarm Antenna Search Algorithm" Entropy 22, no. 7: 765. https://doi.org/10.3390/e22070765

APA Style

Wang, P., Gao, Y., Wu, M., Zhang, F., Li, G., & Qin, C. (2020). A Denoising Method for Fiber Optic Gyroscope Based on Variational Mode Decomposition and Beetle Swarm Antenna Search Algorithm. Entropy, 22(7), 765. https://doi.org/10.3390/e22070765

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop