Research Article  Open Access
DelayFree Tracking Differentiator Design Based on Variational Mode Decomposition: Application on MEMS Gyroscope Denoising
Abstract
This paper presents a delayfree tracking differentiator based on variational mode decomposition (VMD) for extracting the useful signal from a noisy measurement of gyroscope. Sigmoid functionbased tracking differentiator (STD) is a novel tracking differentiator with the advantages of noiseattenuation ability and dynamical performance. However, there is a contradiction in STD; i.e., selecting a larger acceleration factor may cause faster convergence but bad random noise reduction whereas selecting a smaller acceleration factor may lead to signal delay but effective random noise reduction. Here, multiscale transformation is introduced to overcome the contradiction of STD. VMD is selected to decompose the noisy signal into multiscale components, and the correlation coefficients between each component and original signal are calculated, then the component with biggest correlation coefficient is reserved and other components are filtered by the proposed adaptive STD algorithm based on the correlation coefficient of each component, and finally the denoising result is obtained after reconstruction. The prominent advantages of the proposed algorithm are as follows: (i) compared to traditional tracking differentiators, better noise suppression ability can be achieved with suppression of time delay; (ii) compared to other widely used denoising methods, a simpler structure but better denoising ability can be obtained.
1. Introduction
In the early of 1990s, MEMS gyroscope was introduced by the Draper Laboratory. Due to the extensive research and advancements in fabrication technologies and readout electronics, MEMS gyroscope’s performance has been improved over the last twenty years [1]. For further MEMS gyroscope research, a main motivation is the promise of achieving even better performance that competes with conventional bulky and expensive counterparts, but at small size and low cost [2, 3].
In the application of MEMS gyroscope, the noise becomes the main bottleneck which degrades of the signal accuracy. Therefore, it is very important to study the denoising technique to improve the performance of MEMS gyroscope. Many literatures have been dedicating great effort to remove MEMS gyroscope noises. Multiscale transformation method is a widely used technique for gyroscope denoising. In the mentioned waveletbased techniques [4, 5], the proposed method combines a multidimensional gray model (GM) and a wavelet compression technique based on a novel IFOG structure with quadrature demodulation for noise suppression and sensitivity enhancement. In addition, multiscale transformation methods can decompose signal from timedomain into frequencydomain, and the high frequency noise can be eliminated for gyroscope denoising. In mode decompositionbased techniques [6–8], a hybrid algorithm of an optimized local mean decompositionkernel principal component analysis (OLMDKPCA) method is proposed for mechanism structure improvement. The optimal estimation theory is another widely used technique for gyroscope denoising, in which Kalman filter (KF) is the most representative algorithm [9, 10]. Besides, to solve the divergent problem of KF, adaptive Kalman filters (AKF) have been proposed [11–13], such as adaptive sampling strong tracking algorithm (ASSTA), residual based adaptive estimation (RAE), and multiple models based adaptive estimation (MMAE). Furthermore, AKF and variant one are applied to minimize the random noise in interferometric fiberoptic gyroscope (IFOG). ASSTA and scaled unscented Kalman filter (SUKF) algorithm are integrated for denoising the IFOG signal. Experimental results in literatures show that the adjusting KF or AKF mentioned before are suitable for gyroscope denoising in both static and dynamic conditions. However, the specific parameters in systems for the realtime filtering application still need to be improved. Moreover, many other denoising methods have been reported for gyroscope, such as forward linear predication (FLP) [14], timefrequency peak filtering (TFPF) [15], and hybrid denoising algorithms which can describe that an adaptive multiscale method based on the combination generalized morphological filter (CGMF) is presented for denoising [16]. These algorithms are available for gyroscope denoising but still with the problems of complex model or time delay.
In this paper, to develop a novel denoising algorithm for MEMS gyroscope, a Sigmoid functionbased tracking differentiator (STD) based on multiscale decomposition is proposed with advantages of strong noise suppression ability and delayfree. In addition, through both simulations and experiments, the superiority and effectiveness of the proposed denoising algorithm in significantly reducing the noise are verified.
2. Problem Formulation and Innovation Solution
2.1. Sigmoid FunctionBased TD
In order to provide the accurate estimation of derivative of virtual control, STD was first developed with the advantages of being simple structure, global fast convergence, and chatteringfree in differential estimation. Here only a brief introduction of STD will be given. In [18], the detailed description can be found.
Lemma 1 (see [19]). The following system is considered:where stands for the sigmoid function. If there exist , then the system is considered as globally asymptotically stable.
Lemma 2 (see [18]). The following novel tracking differentiator is considered:where R is referred to as acceleration factor and . If there exist , then for any arbitrary bounded and integral function and a constant , the solution of (2) satisfies
It can be concluded that , which is the tracking estimation, averagely converges to the input signal and , which is the differential estimation, converges to the generalized derivative of .
2.2. Problem Formulation and Analysis
Although STD has already shown advantages of dynamical performance and noiseattenuation ability, note that STD is still constrained in supplying a relatively smooth denoising results without signal delay. This statement motivates us to pursue for an improved algorithm model which can radically relax the contradictory and give a better solution in suppressing noises. Here, a simulation noisy signal is utilized to depict the influence of acceleration factor R on the denoising performance of STD, where is random white noise. Figure 1 is the denoising results by using STD with different acceleration factor , in which we can see that bigger means smoother denoising result but severer time delay, while smaller R means worse denoising result but weaker signal delay. Therefore, it is really hard to select an exact R for STD just by experience to guarantee both smoother denoising result and smaller time delay especially for strong random noise.
(a)
(b)
(c)
(d)
(e)
(f)
In Figure 1, we know that noisy signal is composed of three components, which are , , and , respectively. When STD is employed for filtering the noisy signal , it can be considered that STD is worked on the three components. The results of STD applied on the three components as are depicted in Figure 2, Figures 2(a), 2(c), and 2(e) are components , , and . Figures 2(b), 2(d), and 2(f) are the filtered components by STD where the acceleration factor R is set as 0.5. From Figure 1(f) we can see that a favourable noise component suppression result is obtained. But in Figures 1(b) and 1(d), it can be seen that the time delay has occurred. Actually, the two components in Figures 2(a) and 2(c) do not need to be filtered, or in other words the STD with is inappropriate for filtering the two components, which leads to time delay. In traditional application, this phenomenon happens inevitably due to the fact that the components are filtered by STD integrity as one signal. If we can decompose noisy signal into different components with different filtering ways by STD, the time delay phenomenon is expected to be restrained.
(a)
(b)
(c)
(d)
(e)
(f)
2.3. Innovation Solution
After analysis, how to decompose the noisy signal into different component and filter the components with different but appropriate acceleration factors become the key problem that needs to be resolved. As a novel nonrecursive signal processing technique, VMD can adaptively decompose a real valued signal into discrete set of bandlimited subsignals intrinsic mode functions (BLIMFs) owing specific sparsity properties [20]. Thus, in our study, VMD is selected for decomposing noisy signal before the application of STD.
There is another important issue that needs to be resolved, which is how to determine the acceleration factor (R) of STD for each BLIMF after VMD decomposition. The traditional method is determining R just by experience; however, the denoising effect cannot be guaranteed. The correlation coefficient (CC) can be used to reflect the correlation degree between BLIMF and original signal. Hence, CC is used here as a determination criterion for each R of STD. The steps of VMDSTD algorithm are given below [21].
Step 1 (decomposition). The VMD is utilized to decompose the signal into BLIMFs, which are , respectively.
VMD is a novel method of signal decomposition; multicomponent signals are decomposed into BLIMFs, so as to minimize the sum of the bandwidth estimation of each mode. VMD can be represented as a constrained variational problem which is given bywhere is the kth mode, Z(t) is the original signal to be decomposed, is the corresponding center frequency, represents gradient with respect to t, and is the impulse function.
Step 2. Calculate CC for each BLIMF, where m is the length of .
Step 3. Judge the biggest and select the corresponding . Set an origin value of R; then STD is used to track and the tracking result is obtained. Calculate the difference between and ; if the difference is smaller than the defined threshold θ, then R is selected as ; if the difference is bigger than θ, adjust the value of R until is reached. The threshold is defined as (6) by experience, where q is the length of .
After is determined, the other can be calculated by (7), where .
Step 4. Design n STDs for all the BLIMFs by using the calculated R in (7). Then each BLIMF is denoised by the corresponding STD.
Step 5 (reconstruction). The last step is to get the final denoising result by reconstruction. Note that with the biggest is the most useful component of original signal according to the definition of correlation coefficient; therefore will be reserved and this step can make sure that the timedelay would be avoid. The are added together and the final denoising result is obtained.
Since the proposed denoising algorithm is the combination of STD and VMD, it is named as VMDSTD algorithm. The flowchart of VMDSTD algorithm is shown as Figure 3.
3. Simulation
3.1. Simulated Signal Process
In order to present the superiority of the proposed VMDSTD denoising algorithm, the simulated signal is employed for verification. According to the steps of VMDSTD algorithm, firstly is decomposed into BLIMFs as shown in Figure 4. From Figure 4 it can be seen that low frequency components are extracted mainly in BLIMF1 and high frequency noise components are mainly in other BLIMFs. Therefore, STDs for each BLIMF should be designed separately.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
The decomposition is the first step of the proposed VMDSTD algorithm, and the second step is to calculate the CC of each BLIMF. The calculation results are shown in Table 1, from where it can be seen that BLIMF1 has the biggest CC. The third step is applying STD on BLIMF1 to determine the first R according to criterion in Step 3, and the other acceleration factors R can be obtained by (5).

The next step is the application of STD on BLIMFs with determined R. It is noted that BLIMF1 will be reserved without any processing, which can guarantee there is no signal delay of low frequency component. Figure 5 is the denoising results of each BLIMF by using STD with determined R. It can be seen that the details of BLIMF1 are reserved totally; BLIMF2 is filtered by STD with R = 0. 36 which makes part of the details reserved, and signal delay is occurring compared to Figure 4(b); BLIMF3 to BLIMF9 are filtered by STD with bigger R which make a set of smooth denoising results obtained, and it goes without saying that the signal delays are occurring.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
The last step is reconstruction. All the filtered components in Figure 5 are added together and then the reconstruction result is obtained, which is shown as Figure 6(c). Comparison between Figures 6(b) and 6(c) show that the noise is filtered effectively; comparison between Figures 6(a) and 6(c) show that there is almost no signal delay occurring after denoising by the proposed VMDSTD algorithm due to the fact that the low frequency component BLIMF1 is reserved directly; from Figure 6(d) we can see that the residual is still existing, which means the noiseattenuation ability or signal distortion problem of STD still needs to be considered which would be our future work.
(a)
(b)
(c)
(d)
But beyond that, the prior work conducts a simulation process with the single VMDSTD algorithm. To verify our proposed algorithm feasible indeed, some other advanced denoising methods are applied for comparison, like adaptive robust Kalman filter (ARKF), detrended fluctuation analysisVMD (DFAVMD), and empirical mode decomposition–forward linear prediction (EMDFLP). These methods can be also well applied for noise suppression with gyroscope output signal. In Table 2, we added the corresponding noise of the different signaltonoise power ratios (SNRs) to verify the noise reduction capability. Comparing DFAVMD, EMDFLP, VMDSTD, and ARKF methods, the DFAVMD and EMDFLP have similarity in function for noise reduction. The ARKF method has good suppression effect, but the SNR from denoising gyroscope signal of VMDSTD is relatively higher. Consequently, our proposed VMDSTD method is superior to other denoising algorithms.

3.2. Computational Complexity of VMDSTD Denoising Algorithm
In order to assess the complexity of the proposed VMDSTD algorithm, the computation time and space complexity are analysed. Firstly, we assume that the time cost of each operator is the same; therefore only the performance and running hardware are concerned. All required operations, such as addition (ADD), subtraction (SUB), multiplication (MUL), definition (DEF), and division (DIV) are considered. S is introduced to stand for the length of input signal, i.e., the data scale, and N is introduced to stand for the maximum of loops and iterations. The VMD algorithm has already been analysed by [20], as shown in Table 3. The simulated signal in this paper is decomposed into k modes, and the variate N in Table 2 stands for the maximum number of iterations.

In our VMD algorithm, the initialization parameters are set as follows: alpha = 2000, tau = 0, tol = 1e7, and N = 1000. In Table 3, the time and space complexity of VMD are of logarithmic order and linear order O(N), respectively.
In the STD algorithm, the latter relates specifically to the parameter initialization and system function output. From the detailed computation, the STD’s time and space complexity are both of linear order O (N) and linear order O (N).
As shown in Table 4, the time and space complexities for the STD based on the VMD algorithm are listed. The time complexity is of logarithmic order O (, whereas the space complexity is of linear order O (N). In general, we evaluate the magnitude of the complexity . Hence, the results of the analysis demonstrate VMDSTD method is a valid algorithm solvable in polynomial time.

To compare the execution time of different denoising algorithms, a simulation with signals of lengths ranging from to with SNR = 5 dB was performed on a personal computer (Intel® Core™ @2.80 GHz and 8 GB RAM memory) running Windows 10.
Comparing these current algorithms with VMDSTD, the same experimental conditions were applied in the tests. The actual execution times are listed in Table 5. By complexity analysis, EMDFLP time complexity is of O ( and DFAVMD is of O(. The ARKF algorithm has the largest magnitude of O (. The STDVMD time complexity is equal to DFAVMD and EMDFLP whereas the execution times indicate that VMDSTD increases more in time than DFAVMD and less than EMDFLP. It can be explained that the influences of algorithm factors are different. The actual VMDSTD algorithm is with a larger factor than FLP but smaller than DFA. Thus, the proposed algorithm has a simple and rational construction as a whole.

4. Experimental and Verification
4.1. Denoising Test
In this section, the output of gyroscope is employed for verifying the effectiveness of VMDSTD denoising algorithm. MEMS Ssprings vibrating ring gyroscope (MSVRG) [17, 22] is employed in this paper, and the structure is shown in Figure 7; its mode analysis is shown in Figure 8.
(a)
(b)
(c)
(a)
(b)
As Figure 7 shows, a glass substrate with patterned electrode leads, a ring resonator, and twentyfour silicon capacitor electrodes which are evenly distributed around the ring are used to form the MSVRG. On the upper surface of the glass, the patterned metal leads connected to the silicon electrodes are provided to realize the electrical conduction with the measuring and controlling circuit of the gyroscope. In Figure 8, the resonant frequencies in the operating modes of the ring resonator are 8.8078 kHz and 8.8050 kHz, respectively, and the frequency splitting value is 2.8 Hz. Figure 9 is the experimental setup of gyroscope. The equipment is mainly included: temperature control cabinet, a single axis MSVRG, data acquisition system (Multimeter), computer, and Power Supply.
In this work, one set of data is collected from gyroscope with temperature changing from +10°C to 10°C to +10°C, where the temperature change rate is less than 1°C/min. From Figure 10 we can see that, during the temperature variation, there is an obvious drift trend which is submerged in large noises. To extract the precise drift, it is necessary to eliminate the noises effectively. The denoising procedure by using VMDSTD algorithm for gyroscope output signal is the same as section “Simulation.” At the same time, we also apply other advanced denoising methods for signal denoising. ARKF, DFAVMD, EMDFLP, and traditional STD are employed for comparison as shown in Figure 10.
From Figure 10 it can be obviously seen that the noises of original signal are removed. However, the denoising effects are different among different denoising algorithms. When STD algorithm is employed for denoising, a smooth denoising result can be obtained but the signal delay is happening. When DFAVMD, EMDFLP, ARKF, and VMDSTD algorithms are employed for denoising, the signal delay is avoided, and more effective denoising results are obtained. From the principle of DFAAMD we can see that the relevant modes composed of VMD can be selected by DFA and then used for the construction of denoised signal, which means that the BLIMFs with smaller DFAs would be abandon, and the rest of BLIMFs are remaining for reconstruction; therefore the noise in the remaining BLIMFs will be reserved. While in the process of VMDSTD, almost all the BLIMFs would be processed by STD except the BLIMF with biggest correlation coefficient, and the smaller correlated BLIMFs would be suppressed strictly; therefore the denoising performance of VMDSTD algorithm is better. The principle of EMDFLP algorithm is decomposing the signal into IMFs by EMD and then filtering each IMF by FLP, which is similar as our proposed VMDSTD algorithm. However, there is no selection criterion for the parameters of FLP for each IMF which means that it is hard to select the best parameters just by experience; therefore the denoising result of EMDFLP algorithm is not optimal. The denoising result of ARKF is not as smooth as VMDSTD algorithm; that is because ARKF is denoising the signal in timedomain directly, while the VMDSTD algorithm is decomposing the timedomain signal into different BLIMFs firstly and then each BLIMF is denoised independently which makes the denoising process more specific. Besides, compared to ARKF, it does not need to predesign any model in the application of VMDSTD, which makes the structure of the proposed algorithm simpler.
To evaluate the denoising ability of the proposed algorithm, Allan variance analysis is introduced to quantitative comparison. By using Allan variance analysis, the noise coefficients of MEMS gyroscope can be identified and evaluated. Normally, the identified noise coefficients are Q (quantification noise), N (angle random walk), B (bias instability), K (rate random walk), and R (angular rate ramp), respectively, where, N and B are the most important parameters to characterize the noise performance of gyroscope. Simply speaking, N stands for the white noise and B means the 1/f noise or other noises induced by environment. From Figure 11 and Table 6 we can see that the N is 6.84°/h^{1/2} to 0.10°/h^{1/2}, and the B is reduced from 6.84°/h to 0.10°/h after denoising by the proposed VMDSTD algorithm, which performs better than other denoising algorithms.

4.2. Other Discussions
In order to demonstrate the effectiveness of the proposed VMDSTD denoising algorithm, dynamic (step rotation) and bandwidth tests are carried out. Figure 12 is the dynamic test result. The dynamic test is carried out on rotation plat under constant temperature (25°C). It can be concluded that the noises are eliminated, and the useful signals are preserved effectively.
Figure 13 is the bandwidth test result. The bandwidth of our gyroscope is 100 Hz. The input angular rate signal is a sine signal (frequency is 100 Hz) with 2°/s amplitude which is employed for verification. From Figures 13 and 14 we can see that the amplitude of gyroscope output signal keeps 2°/s, and the curve is smooth which indicates that the noise component is decreased obviously. And the results show that the bandwidth is not reduced by the proposed VMDSTD denoising algorithm.
5. Conclusions
To minimize the random noise of MEMS gyroscope, a VMDSTD denoising algorithm is proposed. The main contributions of this paper include the following: firstly, a novel denoising method, named as VMDSTD algorithm, is studied for gyroscope signal processing; secondly the signal delay of STD is solved by combining VMD and STD together; thirdly a VMDbased denoising method is given which makes VMDbased denoising method expanded. The proposed VMDSTD algorithm is verified by the collected MEMS gyroscope data, and the experimental results show that the best denoising result is obtained by using VMDSTD compared to other advanced denoising algorithms, and the temperature drift of MEMS gyroscope can be extracted without signal delay. Additionally, it can be concluded by stationary simulated data and nonstationary experimental output that the proposed denoising algorithm is effective for both stationary and nonstationary signals.
Data Availability
The [Original.mat] data used to support the findings of this study are included within the supplementary information file.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Authors’ Contributions
Xi Zhang and Huiliang Cao contributed equally to this paper.
Acknowledgments
This work was supported in part by the National Natural Science Foundation of China (61603353, 51705477), the PreResearch Field Foundation (6140518010201), the Scientific and Technology Innovation Programs of Higher Education Institutions in Shanxi (201802084), the Program for the Top Young Academic Leaders of Higher Learning Institutions of Shanxi, the Young Academic Leaders Foundation in North University of China, Science Foundation of North University of China (XJJ201822), the Fund for Shanxi “1331 Project” Key Subjects Construction, and the Shanxi Province Outstand Researcher (2016M180018).
Supplementary Materials
The original data are collected from gyroscope with temperature changing from +10°C to 10°C to +10°C, where the temperature change rate is less than 1°C/min. (Supplementary Materials)
References
 H. Cao, Y. Zhang, Z. Han et al., “Polezerotemperature compensation circuit design and experiment for dualmass mems gyroscope bandwidth expansion,” IEEE/ASME Transactions on Mechatronics, 2019. View at: Publisher Site  Google Scholar
 H. L. Cao, Y. Zhang, C. Shen et al., “Temperature energy influence compensation for MEMS vibration gyroscope based on RBF NNGAKF method,” Shock and Vibration, vol. 2018, Article ID 2830686, 10 pages, 2018. View at: Google Scholar
 S. Chong, S. Rui, L. Jie et al., “Temperature drift modeling of MEMS gyroscope based on geneticElman neural network,” Mechanical Systems and Signal Processing, vol. 7273, pp. 897–905, 2016. View at: Publisher Site  Google Scholar
 X. Xu, X. S. Xu, Y. Q. Yao et al., “Optimizationbased alignment for inertial navigation systems: theory and algorithm,” Review of Scientific Instruments, vol. 88, Article ID 035001, 2017. View at: Google Scholar
 X. Guo, J. Tang, J. Li, C. Wang, C. Shen, and J. Liu, “Determine turntable coordinate system considering its nonorthogonality,” Review of Scientific Instruments, vol. 90, no. 3, Article ID 033704, 2019. View at: Publisher Site  Google Scholar
 Q. Zeng, S. Gu, J. Liu et al., “A gyro signal characteristics analysis method based on empirical mode decomposition,” Journal of Sensors, vol. 2016, Article ID 3613878, 7 pages, 2016. View at: Publisher Site  Google Scholar
 C. Shen, J. Yang, J. Tang, J. Liu, and H. Cao, “Note: Parallel processing algorithm of temperature and noise error for microelectromechanical system gyroscope based on variational mode decomposition and augmented nonlinear differentiator,” Review of Scientific Instruments, vol. 89, no. 7, Article ID 076107, 2018. View at: Publisher Site  Google Scholar
 Z. Wang, J. Wang, and W. Du, “Research on fault diagnosis of gearbox with improved variational mode decomposition,” Sensors, vol. 18, no. 10, Article ID 3510, 2018. View at: Publisher Site  Google Scholar
 R. Peesapati, S. L. Sabat, K. P. Karthik, J. Nayak, and N. Giribabu, “Efficient hybrid Kalman filter for denoising fiber optic gyroscope signal,” Optik  International Journal for Light and Electron Optics, vol. 124, no. 20, pp. 4549–4556, 2013. View at: Publisher Site  Google Scholar
 R. Peesapati, S. L. Sabat, K. P. Karthik, M. Narasimhappa, N. Giribabu, and J. Nayak, “FPGAbased embedded platform for fiber optic gyroscope signal denoising,” International Journal of Circuit Theory and Applications, vol. 42, no. 7, pp. 744–757, 2014. View at: Publisher Site  Google Scholar
 M. Narasimhappa, S. L. Sabat, R. Peesapati et al., “Fiberoptic gyroscope signal denoising using an adaptive robust kalman filter,” IEEE Sensors Journal, vol. 16, no. 10, pp. 3711–3718, 2016. View at: Publisher Site  Google Scholar
 M. Narasimhappa, S. L. Sabat, and J. Nayak, “Adaptive sampling strong tracking scaled unscented Kalman filter for denoising the fiber optic gyroscope drift signal. IET Science,” Measurement & Technology, vol. 9, no. 3, pp. 241–249, 2015. View at: Publisher Site  Google Scholar
 H. Guo, Q. Zhu, J. Tang et al., “A temperature and humidity synchronization detection method based on microwave coupledresonator,” Sensors and Actuators B: Chemical, vol. 261, pp. 434–440, 2018. View at: Publisher Site  Google Scholar
 J. Cheng, J. Fang, W. Wu, and J. Li, “Temperature drift modeling and compensation of RLG based on PSO tuning SVM,” Measurement, vol. 55, pp. 246–254, 2014. View at: Publisher Site  Google Scholar
 Y. Wu, C. Shen, H. Cao et al., “Improved morphological filter based on variational mode decomposition for MEMS gyroscope denoising,” Micromachines, vol. 9, no. 5, p. 246, 2018. View at: Publisher Site  Google Scholar
 C. Shen, H. Cao, J. Li et al., “Hybrid denoising approach for fiber optic gyroscopes combing improved empirical mode decomposition and forward linear prediction algorithms,” Review of Scientific Instruments, vol. 87, no. 3, pp. 1–8, 2016. View at: Publisher Site  Google Scholar
 Z. Kou, J. Liu, H. Cao, Y. Shi, J. Ren, and Y. Zhang, “A novel MEMS Ssprings vibrating ring gyroscope with atmosphere package,” AIP Advances, vol. 7, no. 12, Article ID 125301, 2017. View at: Publisher Site  Google Scholar
 X. Shao and H. Wang, “Backstepping robust trajectory linearization control for hypersonic reentry vehicle via novel tracking differentiator,” Journal of The Franklin Institute, vol. 353, no. 9, pp. 1957–1984, 2016. View at: Publisher Site  Google Scholar  MathSciNet
 H. Guo, Y. Chen, D. Wu et al., “Plasmonenhanced sensitivity of spinbased sensors based on a diamond ensemble of nitrogen vacancy color centers,” Optics Expresss, vol. 42, no. 3, pp. 403–406, 2017. View at: Publisher Site  Google Scholar
 K. Dragomiretskiy and D. Zosso, “Variational mode decomposition,” IEEE Transactions on Signal Processing, vol. 62, no. 3, pp. 531–544, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 X. Zhang, H. Cao, X. Shao, J. Liu, and C. Shen, “FOG denoising algorithm based on augmented nonlinear differentiator and singular spectrum analysis,” Applied Sciences, vol. 8, no. 10, Article ID 1710, 2018. View at: Publisher Site  Google Scholar
 Z. Kou, J. Liu, H. Cao et al., “Investigation, modeling, and experiment of an MEMS Ssprings vibrating ring gyroscope,” Journal of Micro/Nanolithography, MEMS, and MOEMS, vol. 17, no. 01, Article ID 015001, 2018. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2019 Xi Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.