Random Error Analysis of MEMS Gyroscope Based on an Improved DAVAR Algorithm

In view that traditional dynamic Allan variance (DAVAR) method is difficult to make a good balance between dynamic tracking capabilities and the confidence of the estimation. And the reason is the use of a rectangular window with the fixed window length to intercept the original signal. So an improved dynamic Allan variance method was proposed. Compared with the traditional Allan variance method, this method can adjust the window length of the rectangular window adaptively. The data in the beginning and terminal interval was extended with the inverted mirror extension method to improve the utilization rate of the interval data. And the sliding kurtosis contribution coefficient and kurtosis were introduced to adjust the length of the rectangular window by sensing the content of shock signal in terminal interval. The method analyzed the window length change factor in different stable conditions and adjusted the rectangular window’s window length according to the kurtosis, sliding kurtosis contribution coefficient. The test results show that the more the kurtosis stability threshold was close to 3, the stronger the dynamic tracking ability of DAVAR would be. But the kurtosis stability threshold was too close to 3, there was a misjudgement in kurtosis analysis of signal stability, resulting in distortion of DAVAR analysis results. When using the improved DAVAR method, the kurtosis stability threshold can be close to 3 to improve the tracking ability and the estimation confidence in stable interval. Therefore, it solved the problem that the dynamic Allan variance tracking ability and confidence level were difficult to take into account, and also solved the problem of misjudgement in the stability analysis of kurtosis.


Introduction
Micro electromechanical system (MEMS) gyroscope has the advantages of small size, light weight, and low cost [1][2][3][4], which is widely used in various fields of modern navigation. However, the accuracy of MEMS gyroscopes is low. So it is important to compensate the errors of MEMS gyroscopes.
MEMS gyroscope error can be divided into deterministic error and random error [5]. Deterministic error can be eliminated by calibrating the MEMS gyroscope [6]. The random error has a strong randomness, so the random error becomes one of the main factors that limit the accuracy of the MEMS gyroscope. At present, the modeling methods of random error about MEMS gyroscopes mainly include autocorrelation analysis, power spectral density method and Allan variance analysis method. Allan variance analysis can characterize and identify various random error sources of MEMS gyroscopes. It is convenient to establish accurate random error model and improve navigation precision. Many scholars such as Gao et al., Chen et al. and Li et al. have made Allan variance, overlapped Allan variance identification for random errors of MEMS gyroscopes and achieved good results [7][8][9]. However, MEMS gyroscope signals are non-stationary, so the Allan variance is difficult to track the dynamic changes of the signal. Wang et al. use a rectangular window to intercept the original Allan variance is defined as shown in Equation (1).
In Equation (1), <> represents the overall average of the data. MEMS gyroscope random error usually contains angle random walk (ARW, A is a coefficient), bias instability (BI, B is a coefficient), rate random walk (RRW, R is a coefficient), rate ramp (RR, K is a coefficient), quantization noise (QN, Q is a coefficient), Markov noise (MN) and sinusoidal noise (SN). Due to the different types of gyroscopes and test environments, random noise of various components may exist in the output data. It is hypothesized that the noise is independent of each other and the Allan variance can be expressed as the sum of squares of the random errors, as follows: Among them, the Allan variance of the five random errors including QN, ARW, BI, RRW, and RR is as following Equations (2)~ (6). Each error source can be expressed as proportional to the −2~2 power of the form, therefore, Allan variance can be expressed as Equation (7).
Micromachines 2018, 9, x 3 of 22 In Equation (1), <> represents the overall average of the data. MEMS gyroscope random error usually contains angle random walk (ARW, A is a coefficient), bias instability (BI, B is a coefficient), rate random walk (RRW, R is a coefficient), rate ramp (RR, K is a coefficient), quantization noise (QN, Q is a coefficient), Markov noise (MN) and sinusoidal noise (SN). Due to the different types of gyroscopes and test environments, random noise of various components may exist in the output data. It is hypothesized that the noise is independent of each other and the Allan variance can be expressed as the sum of squares of the random errors, as follows: Among them, the Allan variance of the five random errors including QN, ARW, BI, RRW, and RR is as following Equations (2)~ (6). Each error source can be expressed as proportional to the −2~2 power of the form, therefore, Allan variance can be expressed as Equation (7).

Allan Variance of Analog MEMS Gyroscope Signal
The output signal of MEMS gyroscope is simulated by setting different amplitude Gaussian white noise at different time. The output signal was set the sampling step to 1 s and the sampling time to 6000 s. In the stationary period, the variance of the analog signal is 1 and the variance of the abrupt signal (3000~4000 s) is 3. The analog signal is obtained as shown in

Allan Variance of Analog MEMS Gyroscope Signal
The output signal of MEMS gyroscope is simulated by setting different amplitude Gaussian white noise at different time. The output signal was set the sampling step to 1 s and the sampling time to 6000 s. In the stationary period, the variance of the analog signal is 1 and the variance of the abrupt signal (3000~4000 s) is 3. The analog signal is obtained as shown in Figure 2. The Allan standard deviation curve is shown in Figure 3, which is obtained by the simulated MEMS gyro signal analysis with Allan variance. A large number of experiments show that different error items usually appear in different time intervals. A certain or a few kinds of random errors play a leading role in different time intervals. Therefore, there is a quantization noise with a slope of −1 in Figure 3. Moreover, angular rate random walk is a kind of random error with long correlation time. Therefore, the random errors such as angular rate random walk and so on are not well represented in Figure 3. Using the least square method to fit the Allan standard deviation curve, the error coefficient of angle random walk can be obtained as 2.36 × 10 −3°/ h 1/2 . The least square fitting results are as shown in Figure 3. The Allan variance has the advantage of dealing with the ideal stationary signal. But the Allan variance has not strong tracking ability and it has the disadvantage of dealing with the non-stationary signal which is shown in Figure 2. Therefore, many scholars intercept the original signal with window function, and then the intercepted signal was respectively processed with Allan variance analysis. The above method is dynamic Allan variance (DAVAR) [19,20].

Dynamic Allan Variance (DAVAR) Analysis
The basic idea of dynamic Allan variance is firstly to design a proper window function to intercept the original signal. Next, the Allan variance of intercepted signal is calculated separately, and then all the results are combined into a three-dimensional graph to represent the dynamic performance of the signal [19,20]. The principle is as follows: (1) t1 is defined as a fixed analysis time t, t = t1.
(2) x1(t) is obtained by capturing the original signal with a rectangular window whose width is T and center is on t1.  The Allan standard deviation curve is shown in Figure 3, which is obtained by the simulated MEMS gyro signal analysis with Allan variance. A large number of experiments show that different error items usually appear in different time intervals. A certain or a few kinds of random errors play a leading role in different time intervals. Therefore, there is a quantization noise with a slope of −1 in Figure 3. Moreover, angular rate random walk is a kind of random error with long correlation time. Therefore, the random errors such as angular rate random walk and so on are not well represented in Figure 3. Using the least square method to fit the Allan standard deviation curve, the error coefficient of angle random walk can be obtained as 2.36 × 10 −3• /h 1/2 . The least square fitting results are as shown in Figure 3. The Allan standard deviation curve is shown in Figure 3, which is obtained by the simulated MEMS gyro signal analysis with Allan variance. A large number of experiments show that different error items usually appear in different time intervals. A certain or a few kinds of random errors play a leading role in different time intervals. Therefore, there is a quantization noise with a slope of −1 in Figure 3. Moreover, angular rate random walk is a kind of random error with long correlation time. Therefore, the random errors such as angular rate random walk and so on are not well represented in Figure 3. Using the least square method to fit the Allan standard deviation curve, the error coefficient of angle random walk can be obtained as 2.36 × 10 −3°/ h 1/2 . The least square fitting results are as shown in Figure 3. The Allan variance has the advantage of dealing with the ideal stationary signal. But the Allan variance has not strong tracking ability and it has the disadvantage of dealing with the non-stationary signal which is shown in Figure 2. Therefore, many scholars intercept the original signal with window function, and then the intercepted signal was respectively processed with Allan variance analysis. The above method is dynamic Allan variance (DAVAR) [19,20].

Dynamic Allan Variance (DAVAR) Analysis
The basic idea of dynamic Allan variance is firstly to design a proper window function to intercept the original signal. Next, the Allan variance of intercepted signal is calculated separately, and then all the results are combined into a three-dimensional graph to represent the dynamic performance of the signal [19,20]. The principle is as follows: (1) t1 is defined as a fixed analysis time t, t = t1.
(2) x1(t) is obtained by capturing the original signal with a rectangular window whose width is T and center is on t1.  The Allan variance has the advantage of dealing with the ideal stationary signal. But the Allan variance has not strong tracking ability and it has the disadvantage of dealing with the non-stationary signal which is shown in Figure 2. Therefore, many scholars intercept the original signal with window function, and then the intercepted signal was respectively processed with Allan variance analysis. The above method is dynamic Allan variance (DAVAR) [19,20].

Dynamic Allan Variance (DAVAR) Analysis
The basic idea of dynamic Allan variance is firstly to design a proper window function to intercept the original signal. Next, the Allan variance of intercepted signal is calculated separately, and then all the results are combined into a three-dimensional graph to represent the dynamic performance of the signal [19,20]. The principle is as follows: (1) t 1 is defined as a fixed analysis time t, t = t 1 . (2) x 1 (t) is obtained by capturing the original signal with a rectangular window whose width is T and center is on t 1 . (3) The Allan standard deviation σ(t 1 , τ) is obtained by calculating the clipped signal x 1 (t). (4) t 2 is defined as another analysis time point t, t = t 2 . The two interception signals overlap through repeat step (2), and then steps (3) and (4) are repeated, respectively. The Allan standard deviation sequence σ(t 1 , τ), σ(t 2 , τ) · · · σ(t n , τ) is finally obtained.
In the window section (t − T/2) ≤ t ≤ (t + T/2), the original signal is captured using the rectangular window function P τ (t), and the result is obtained as follows.
In the above formula, t is the center of the rectangular window and the Allan window h τ (t ) is defined as follows.
Taking Equation (11) into Equation (1), we get: In summary, the expected value of (12) is defined as the dynamic Allan variance.
The square root of (13) is defined as the dynamic Allan standard deviation. Select the window function length of 900. The dynamic Allan variance for the analog signal is shown in Figure 4. ( 1 τ σ t is obtained by calculating the clipped signal x1(t).
(4) t2 is defined as another analysis time point t, t = t2. The two interception signals overlap through repeat step (2), and then steps (3) and (4) are repeated, respectively. The Allan standard deviation sequence is finally obtained.
In the window section , the original signal is captured using the rectangular window function ) (t P τ , and the result is obtained as follows.
In the above formula, t is the center of the rectangular window and the Allan window Taking Equation (11) into Equation (1), we get: (12) In summary, the expected value of (12) is defined as the dynamic Allan variance.
The square root of (13) is defined as the dynamic Allan standard deviation. Select the window function length of 900. The dynamic Allan variance for the analog signal is shown in Figure 4.  Figure 4 shows the dynamic Allan variance of the simulated gyroscope signal. In Figure 4, the dynamic Allan variance can show the dynamic change of the signal fluctuation. Although the dynamic Allan variance can track the dynamic performance of the signal, but the tracking result is not perfect from the analysis of the error coefficient in Figure 5. This is because the dynamic Allan variance intercepts the original signal with a fixed-window rectangular window function. There are some shortcomings in the method of intercepting a signal with a fixed window length rectangular window function. If the window length is short, the amount of data per section is small and the confidence of the analysis results is low. If the length of the window is longer, the dynamic changes  Figure 4 shows the dynamic Allan variance of the simulated gyroscope signal. In Figure 4, the dynamic Allan variance can show the dynamic change of the signal fluctuation. Although the dynamic Allan variance can track the dynamic performance of the signal, but the tracking result is not perfect from the analysis of the error coefficient in Figure 5. This is because the dynamic Allan variance intercepts the original signal with a fixed-window rectangular window function. There are some shortcomings in the method of intercepting a signal with a fixed window length rectangular window function. If the window length is short, the amount of data per section is small and the confidence of the analysis results is low. If the length of the window is longer, the dynamic changes of the signal are sensed very early and at the same time the dynamic changes disappear very late. As a result, the analysis has the disadvantages of early sensitivity to the beginning and sensitive delay to the end of dynamic signal. of the signal are sensed very early and at the same time the dynamic changes disappear very late. As a result, the analysis has the disadvantages of early sensitivity to the beginning and sensitive delay to the end of dynamic signal.

Kurtosis and Sliding Kurtosis Contribution Coefficient
The kurtosis is a normalized fourth-order central moment, dimensionless and can reflect the distribution of the signal. For a signal of length n, the kurtosis coefficient is as in Equation (14).
In Equation (14), μ is the mean value of signal i x , and σ is the standard deviation of signal i x . K = 3 (The signal distribution curve has a normal kurtosis); K > 3 (The peak value of the impact component in the signal is above the normal distribution curve); K < 3 (The impact component amplitude is less than the normal distribution). Wang et al. use rectangular windows to intercept the original signal and then the kurtosis value is calculated to adjust the window length [10,11]. However, the kurtosis of the intercepted signal has the characteristics of a rectangular window intercepting signal, which leads to the shortcomings of early sensitivity to the beginning and sensitive delay to the end of dynamic signal. When using a rectangular window with a window length of 201 and 601, the original signal is intercepted and its kurtosis value is calculated. The kurtosis and the comparison between original signal and the intercepted signal are shown in Figure 6. Through the above analysis, it is difficult to track the signal fluctuations by using the window function intercepting the signal and then calculating the interception signal kurtosis value. The

Kurtosis and Sliding Kurtosis Contribution Coefficient
The kurtosis is a normalized fourth-order central moment, dimensionless and can reflect the distribution of the signal. For a signal of length n, the kurtosis coefficient is as in Equation (14).
In Equation (14), µ is the mean value of signal x i , and σ is the standard deviation of signal x i . K = 3 (The signal distribution curve has a normal kurtosis); K > 3 (The peak value of the impact component in the signal is above the normal distribution curve); K < 3 (The impact component amplitude is less than the normal distribution). Wang et al. use rectangular windows to intercept the original signal and then the kurtosis value is calculated to adjust the window length [10,11]. However, the kurtosis of the intercepted signal has the characteristics of a rectangular window intercepting signal, which leads to the shortcomings of early sensitivity to the beginning and sensitive delay to the end of dynamic signal. When using a rectangular window with a window length of 201 and 601, the original signal is intercepted and its kurtosis value is calculated. The kurtosis and the comparison between original signal and the intercepted signal are shown in Figure 6. of the signal are sensed very early and at the same time the dynamic changes disappear very late. As a result, the analysis has the disadvantages of early sensitivity to the beginning and sensitive delay to the end of dynamic signal.

Kurtosis and Sliding Kurtosis Contribution Coefficient
The kurtosis is a normalized fourth-order central moment, dimensionless and can reflect the distribution of the signal. For a signal of length n, the kurtosis coefficient is as in Equation (14).
In Equation (14), μ is the mean value of signal i x , and σ is the standard deviation of signal i x . K = 3 (The signal distribution curve has a normal kurtosis); K > 3 (The peak value of the impact component in the signal is above the normal distribution curve); K < 3 (The impact component amplitude is less than the normal distribution). Wang et al. use rectangular windows to intercept the original signal and then the kurtosis value is calculated to adjust the window length [10,11]. However, the kurtosis of the intercepted signal has the characteristics of a rectangular window intercepting signal, which leads to the shortcomings of early sensitivity to the beginning and sensitive delay to the end of dynamic signal. When using a rectangular window with a window length of 201 and 601, the original signal is intercepted and its kurtosis value is calculated. The kurtosis and the comparison between original signal and the intercepted signal are shown in Figure 6. Through the above analysis, it is difficult to track the signal fluctuations by using the window function intercepting the signal and then calculating the interception signal kurtosis value. The  Through the above analysis, it is difficult to track the signal fluctuations by using the window function intercepting the signal and then calculating the interception signal kurtosis value. The kurtosis value is more stable with the longer window length, but the worse the tracking results are. In order to accurately track the dynamic characteristics of the signal, the shorter window length of the rectangular window will be selected. When the window length is very short, the kurtosis value is volatile, and it is difficult to track fluctuations of signal. As shown in Figure 6, the kurtosis of window length 11 changes dramatically which make the kurtosis difficult track the dynamic response of the signal. Therefore, this paper proposes the sliding kurtosis contribution coefficient based on the idea of kurtosis contribution coefficient [21]. The principle of sliding kurtosis contribution coefficient is as follows. K is kurtosis of the sampling sequence x i (i = 1, 2, · · · n). x ωp is obtained by intercepting the original signal with a L length rectangular window at the p data point in time t p .
In the equation above, x mp is defined as the sequence in the original signal sequence excluding the x ωp sequence. Though calculation, the kurtosis of the x mp sequence is as follows.
In the equation above, µ mp is the mean value of x mp and σ mp is the standard deviation of x mp . The sliding kurtosis contribution coefficient of sequence x ωp is as follows.
Sliding kurtosis contribution coefficient of different sequences in different analysis points is calculated, which can accurately estimate the abrupt time of the signal. For the above analog gyro signal, when the length of the rectangular window is T = 11, the sliding kurtosis contribution coefficient is shown in Figure 7. kurtosis value is more stable with the longer window length, but the worse the tracking results are. In order to accurately track the dynamic characteristics of the signal, the shorter window length of the rectangular window will be selected. When the window length is very short, the kurtosis value is volatile, and it is difficult to track fluctuations of signal. As shown in Figure 6, the kurtosis of window length 11 changes dramatically which make the kurtosis difficult track the dynamic response of the signal. Therefore, this paper proposes the sliding kurtosis contribution coefficient based on the idea of kurtosis contribution coefficient [21]. The principle of sliding kurtosis contribution coefficient is as follows. K is kurtosis of the sampling sequence i ). xωp is obtained by intercepting the original signal with a L length rectangular window at the p data point in time tp.
In the equation above, xmp is defined as the sequence in the original signal sequence excluding the xωp sequence. Though calculation, the kurtosis of the xmp sequence is as follows.
In the equation above, μmp is the mean value of xmp and σmp is the standard deviation of xmp. The sliding kurtosis contribution coefficient of sequence xωp is as follows.
Sliding kurtosis contribution coefficient of different sequences in different analysis points is calculated, which can accurately estimate the abrupt time of the signal. For the above analog gyro signal, when the length of the rectangular window is T = 11, the sliding kurtosis contribution coefficient is shown in Figure 7. It can be seen from Figure 7 that the sliding kurtosis contribution coefficient which can be applied to the rectangular window with a short length can accurately reflect the dynamic change of the signal.  It can be seen from Figure 7 that the sliding kurtosis contribution coefficient which can be applied to the rectangular window with a short length can accurately reflect the dynamic change of the signal.

Extension
When the performance of signal is analyzed with a rectangular window, the problem of low utilization of the data in the beginning and end of the signal is inevitable. In summary, the beginning and end of the signal data is extended with inverted mirror method of extension, which can enhance the utilization of signal interval data.
When the center of the rectangular window is at 1, L 2 ) ∪ (N − L 2 , N , the standard deviation of the intercepted sequence is calculated as follows.
In the equation above, L is the window length of the rectangular window, m generally take 1/3 the length of the data column, τ 0 is the sampling time.

Improved Window Length Adjusting Adaptively Algorithm
The traditional DAVAR algorithm intercept the original signal using a fixed window length rectangular window and analyze the dynamic Allan variance of the intercepted data, which has disadvantages of early sensitivity to the beginning and sensitive delay to the end of dynamic signal. A rectangular window function is designed by combining slipping kurtosis contribution coefficient and kurtosis. The rectangular windows with longer window lengths are used in the smoothing interval of the signal, which can enhance the confidence and stability of the analysis results. If the signal in the strong changes interval, the signal is dealt with short-window length rectangular window, which can accurately track the dynamic changes of signals and improve the tracking ability of dynamic Allan variance. So the improved window length adjusting adaptively algorithm is as follows.
(1) The L window is defined as a rectangular window with a suitable window length of L. Using the L window to intercept the original signal and analyze the kurtosis of the intercepted signal. The kurtosis threshold for stable signal is [K 1 ,K 2 ], such that at the analysis moment t p , the kurtosis value of this sequence K p is calculated and the following relation is obtained.
(2) The l window is defined as a rectangular window with a suitable window length of l (l < L). The slipping kurtosis contribution coefficient of the signal sequence that is obtained by intercepting the signal from the end of a sequence cut from L window with l window is calculated. Through calculating the kurtosis value K mp of the front L − l length signal in the L window intercepting signal, the sliding kurtosis contribution coefficient S = K p − K mp /K p of the l window is calculated. S reflects the stability of the signal at the end of the L window and achieves an accurate sensitivity to the end of the intercepted signal sequence. The sliding kurtosis contribution coefficient in [S 1 ,S 2 ] is designed to represent the stable signal, which leads to the existence of Equation (17) on the t p moment.
The function f 2 reflects whether there is an impact component at the t p + L−l 2 τ 0 moment. Where τ 0 is the sampling interval.
(3) The above two output function values are combined and then get the following relationship.
When f (t p ) is 0, the characterization of the signal is stationary, and the end of the intercept signal contains no impact components. When f (t p ) is 1, the signal is unstable, and the rising edge of f (t p ) is collected, which is known as the impact component, exists at time t p + L−l 2 τ 0 . (4) When the dynamic Allan analysis of variance is analyzed, the maximum window length L max for a steady signal and the minimum window length L min for an unstable signal are usually designed. When f (t p ) is 1, the signal contains the impact component at the t p + L−l 2 τ 0 moment, and the L window starts to capture the signal at time t p + L−l 2 τ 0 with the minimum window length L min . Assuming that the initial window length of the L window is L , the maximum distance of the center about the L window moving distance is L −L min 2 , and the design of the window length variation factor ∆L 1 is as follows.
When f (t p ) = 1, and at the next point t p+1 , the L window's window length L(t p+1 ) is as follows: When f (t p ) = 0, the intercepted signal is stable and the starting point t p − L 2 τ 0 of the intercepted sequence contains no impact component. In this case, it is required that L window after increasing the window length can not acquire the signal before the t p − L 2 τ 0 time again. At this time, the window length of the L window is increased until it reaches L max . The maximum distance of L window center point moving is L max −L 2 , then the window length variation factor ∆L 0 is designed at this time.
And at the next point t p+1 , the L window's window length L(t p+1 ) is as follows.

Simulation and Test
Taking the rectangular window with a window length of 600 as an example, the adaptive change of window length based on kurtosis is analyzed as follows. The analysis results are as in Figure 8. If the stability threshold of kurtosis is 4, and the sensitive nodes of dynamic signals are 2700 s, 3115 s, 3963 s, and 4300 s. When kurtosis is used to control the change of window length, reduce the length of the window at 2700 s, and the beginning of the rectangle window just collects the 3000 s data; increase the length of window at 3115 s, the end of the rectangular window collects the 2815 s data (3000 s data should be collected), and the signal is in a stable range, so the adjustment of window length is obviously too advanced. When the rectangle window is at 3963 s, the window length begins to decrease to sense the dynamic change of 4000 s data. But the 4263 s data is collected at the beginning of the rectangular window, and the signal is already in a stable state at this time, so there is a serious lag phenomenon in adaptive adjustment of window length.
Taking the rectangular window with a window length of 600 as an example, the adaptive change of window length based on kurtosis is analyzed as follows. The analysis results are as in Figure 8. If the stability threshold of kurtosis is 4, and the sensitive nodes of dynamic signals are 2700 s, 3115 s, 3963 s, and 4300 s. When kurtosis is used to control the change of window length, reduce the length of the window at 2700 s, and the beginning of the rectangle window just collects the 3000 s data; increase the length of window at 3115 s, the end of the rectangular window collects the 2815 s data (3000 s data should be collected), and the signal is in a stable range, so the adjustment of window length is obviously too advanced. When the rectangle window is at 3963 s, the window length begins to decrease to sense the dynamic change of 4000 s data. But the 4263 s data is collected at the beginning of the rectangular window, and the signal is already in a stable state at this time, so there is a serious lag phenomenon in adaptive adjustment of window length. If the stability threshold of kurtosis is 3.45, the sensitive nodes of dynamic signals are 2700 s, 3201 s, 3868 s, 4300 s. When kurtosis is used to control the change of window length, reduce the length of the window at 2700 s, and the beginning of the rectangle window just collects the 3000 s data; increase the length of window at 3201 s, the end of the rectangular window collects the 2901 s data (3000 s data should be collected), and there is still an adjustment ahead of time. When the rectangle window is at 3868 s, the window length begins to decrease to sense the dynamic change of 4000 s data. But the 4168 s data is collected at the beginning of the rectangular window, and there is also a lagging control.
If the stability threshold of kurtosis is 3.3, the rectangle window begins to increase the length of the window at 3228 s, and the end of the rectangular window collects the 2928 s data (3000 s data should be collected); the rectangle window begins to reduce the length of the window at 3840 s, and the 4140 s data is collected at the beginning of the rectangular window.
So the dynamic tracking effect is best when the stability threshold of kurtosis is 3.3. But at this time, the kurtosis value of the stationary interval has exceeded the kurtosis stability threshold, that is, the threshold value of the kurtosis is relatively small. From Figure 8, we can see that the closer the kurtosis threshold is to 3, the more sensitive the adjustment of the rectangular window is, and the better the dynamic tracking effect is. Because in 3000~4000 s, the variance of the noise is increased to 3, that is, there is a sudden increase in the amplitude of noise at 3000 s. According to Equation (14),  If the stability threshold of kurtosis is 3.45, the sensitive nodes of dynamic signals are 2700 s, 3201 s, 3868 s, 4300 s. When kurtosis is used to control the change of window length, reduce the length of the window at 2700 s, and the beginning of the rectangle window just collects the 3000 s data; increase the length of window at 3201 s, the end of the rectangular window collects the 2901 s data (3000 s data should be collected), and there is still an adjustment ahead of time. When the rectangle window is at 3868 s, the window length begins to decrease to sense the dynamic change of 4000 s data. But the 4168 s data is collected at the beginning of the rectangular window, and there is also a lagging control.
If the stability threshold of kurtosis is 3.3, the rectangle window begins to increase the length of the window at 3228 s, and the end of the rectangular window collects the 2928 s data (3000 s data should be collected); the rectangle window begins to reduce the length of the window at 3840 s, and the 4140 s data is collected at the beginning of the rectangular window.
So the dynamic tracking effect is best when the stability threshold of kurtosis is 3.3. But at this time, the kurtosis value of the stationary interval has exceeded the kurtosis stability threshold, that is, the threshold value of the kurtosis is relatively small. From Figure 8, we can see that the closer the kurtosis threshold is to 3, the more sensitive the adjustment of the rectangular window is, and the better the dynamic tracking effect is. Because in 3000~4000 s, the variance of the noise is increased to 3, that is, there is a sudden increase in the amplitude of noise at 3000 s. According to Equation (14), we know that the kurtosis value at 3000 s (the center of the rectangular window is in 2700 s) is positively correlated with the difference of noise variance before and after mutation. When the rectangular window intercepts the abrupt signal, the kurtosis will decrease gradually, no matter how large the variance of the mutation signal is. So the kurtosis value changes slowly between 3000~4000 s. The slow change leads that the adjustment of rectangular windows based on kurtosis is insensitive in this area. The sliding kurtosis contribution coefficient is used to characterize the signal in the kurtosis slowly changing region. So the stability threshold of kurtosis can be achieved closer to 3, and the better the dynamic tracking ability of improved DAVAR will be.
According to the analog MEMS gyroscope designed in Section 2.2, the random error of signal is analyzed by the improved dynamic Allan variance method. Setting parameters: L max = 600, L min = 400; K 1 = 2.5, K 2 = 3.3; S 1 = 0, S 2 = 0.07; l = 100; ∆L 0 = ∆L 1 = 2. The l window's window length and above parameters can be adjusted according to the sensitivity of the impact components in signal. The kurtosis value of the captured signal which is obtained by capturing original signal with rectangular windows with window length of 600, 400, or adaptive window length are obtained, respectively, as shown in Figure 9. It can be found that before the impact signal is intercepted, the kurtosis value of the adaptive window length rectangular window intercepting signal coincides with the kurtosis value of the 600 window length rectangular window intercepting signal. And after the impact signal, the kurtosis value of the adaptive window length rectangular window intercepting signal coincides with the kurtosis value of the 400 window length rectangular window intercepting signal.
how large the variance of the mutation signal is. So the kurtosis value changes slowly between 3000~4000 s. The slow change leads that the adjustment of rectangular windows based on kurtosis is insensitive in this area. The sliding kurtosis contribution coefficient is used to characterize the signal in the kurtosis slowly changing region. So the stability threshold of kurtosis can be achieved closer to 3, and the better the dynamic tracking ability of improved DAVAR will be.
According to the analog MEMS gyroscope designed in Section 2.2, the random error of signal is analyzed by the improved dynamic Allan variance method. Setting parameters: Lmax = 600, Lmin = 400; K1 = 2.5, K2 = 3.3; S1 = 0, S2 = 0.07; l = 100; ∆L0 = ∆L1 = 2. The l window's window length and above parameters can be adjusted according to the sensitivity of the impact components in signal. The kurtosis value of the captured signal which is obtained by capturing original signal with rectangular windows with window length of 600, 400, or adaptive window length are obtained, respectively, as shown in Figure 9. It can be found that before the impact signal is intercepted, the kurtosis value of the adaptive window length rectangular window intercepting signal coincides with the kurtosis value of the 600 window length rectangular window intercepting signal. And after the impact signal, the kurtosis value of the adaptive window length rectangular window intercepting signal coincides with the kurtosis value of the 400 window length rectangular window intercepting signal. The simulation test is carried out according to the above parameters. And the signal kurtosis, f1, sliding kurtosis contribution coefficient, f2, f and adaptive window length is shown in Figure 10. It can be seen from Figure 10 that the adaptive window length becomes the minimum window length at the 2803 s, and the data of the 3003 s is intercepted which is at the end of the intercepted signal and which is the starting point of the shock signal. Compared with the sliding kurtosis contribution coefficient, kurtosis can characterize the signal stability better, but it is insensitive to the end position of the impact component in the signal. As shown in Figure 10, the interval [K1,K2] is relatively small (Some slightly fluctuating signals are shown unstable and f1 jumps to 1). In this case, the shock signal is only perceived when the L window reaches 3840 s. And at this time, the end of the L window intercepting signal is the data at 4140 s. So the design of adaptive window length simply according to kurtosis is difficult to meet the requirements. Compared to the sensitivity of the sliding kurtosis contribution coefficient, the L window is already sensitive to the end of the shock signal when it reaches 3743 s. And at this time, the end of the L window intercepting signal is the data at 4043 s. Therefore, the reduction of window length at this time can improve the sensitivity of the impact component. The simulation test is carried out according to the above parameters. And the signal kurtosis, f 1 , sliding kurtosis contribution coefficient, f 2 , f and adaptive window length is shown in Figure 10. It can be seen from Figure 10 that the adaptive window length becomes the minimum window length at the 2803 s, and the data of the 3003 s is intercepted which is at the end of the intercepted signal and which is the starting point of the shock signal. Compared with the sliding kurtosis contribution coefficient, kurtosis can characterize the signal stability better, but it is insensitive to the end position of the impact component in the signal. As shown in Figure 10, the interval [K 1 ,K 2 ] is relatively small (Some slightly fluctuating signals are shown unstable and f 1 jumps to 1). In this case, the shock signal is only perceived when the L window reaches 3840 s. And at this time, the end of the L window intercepting signal is the data at 4140 s. So the design of adaptive window length simply according to kurtosis is difficult to meet the requirements. Compared to the sensitivity of the sliding kurtosis contribution coefficient, the L window is already sensitive to the end of the shock signal when it reaches 3743 s. And at this time, the end of the L window intercepting signal is the data at 4043 s. Therefore, the reduction of window length at this time can improve the sensitivity of the impact component. According to the improved DAVAR method, the dynamic Allan variance of the analog signal is shown in Figure 11.  Figure 12 shows the random walk analysis results of the signal sequences which are obtained through intercepting the original signal by using the rectangular windows with different window lengths. It can be seen that the rectangular window with the window length of 600 is sensitive to the impact signal very early, and the recovery is later. The rectangular window with the window length of 400 is more sensitive to the impact signal, but the analysis result fluctuates obviously. The rectangular window with adaptive window length can accurately detect the impact component in the signal, and it ensures the stability of the analysis results.  According to the improved DAVAR method, the dynamic Allan variance of the analog signal is shown in Figure 11. According to the improved DAVAR method, the dynamic Allan variance of the analog signal is shown in Figure 11.  Figure 12 shows the random walk analysis results of the signal sequences which are obtained through intercepting the original signal by using the rectangular windows with different window lengths. It can be seen that the rectangular window with the window length of 600 is sensitive to the impact signal very early, and the recovery is later. The rectangular window with the window length of 400 is more sensitive to the impact signal, but the analysis result fluctuates obviously. The rectangular window with adaptive window length can accurately detect the impact component in the signal, and it ensures the stability of the analysis results.   Figure 12 shows the random walk analysis results of the signal sequences which are obtained through intercepting the original signal by using the rectangular windows with different window lengths. It can be seen that the rectangular window with the window length of 600 is sensitive to the impact signal very early, and the recovery is later. The rectangular window with the window length of 400 is more sensitive to the impact signal, but the analysis result fluctuates obviously. The rectangular window with adaptive window length can accurately detect the impact component in the signal, and it ensures the stability of the analysis results. According to the improved DAVAR method, the dynamic Allan variance of the analog signal is shown in Figure 11.  Figure 12 shows the random walk analysis results of the signal sequences which are obtained through intercepting the original signal by using the rectangular windows with different window lengths. It can be seen that the rectangular window with the window length of 600 is sensitive to the impact signal very early, and the recovery is later. The rectangular window with the window length of 400 is more sensitive to the impact signal, but the analysis result fluctuates obviously. The rectangular window with adaptive window length can accurately detect the impact component in the signal, and it ensures the stability of the analysis results.  The angle random walk of gyro is identified by improved DAVAR method. The comparison of confidence level and dynamic tracking ability of the analysis results between improved DAVAR (Adaptive window length DAVAR, A-DAVAR), DAVAR-400 (window length is 400) and DAVAR-600 (window length is 600) is shown in Table 1. In Table 1, the beginning (3000 s) means that the window first senses the 3000 s data, and the window length should decrease at this time. The beginning (3000 s) should be close to DAVAR-600's beginning (3000 s). The end (3000 s) means that the window lost senses the 3000 s data, and the window length should increase at this time. The end (3000 s) should be close to DAVAR-400's end (3000 s). The ARW (0~2000 s) means the mean value from 0 s to 2000 s and the region is in stationary interval, so the ARW (0~2000 s) should be close to the refer value in Section 2.2.
From Table 1, we can see that the improved dynamic Allan variance solves the problem of both confidence and dynamic tracking ability.
The MEMS gyro is sensitive to temperature, so the MEMS gyro is fixed to the stationary preheating 1 h on the turntable to make the temperature reach steady state. The accurate angular velocity can be input to the MEMS gyro using the turntable. But the input angular velocity of the turntable is 0 in the experiment. Since the gyroscope bandwidth is 40 Hz, the sampling frequency is 100 Hz according to sampling theorem.
The gyro data collected for the 600 s is shown in Figure 13. When the turntable input angle speed is 0, the output signal of MEMS gyroscope contains trend and random drift. Usually the mean value method and polynomial fitting method are used to remove the trend terms. So the random drift of MEMS gyro is obtained by removing trend items from data collection with 3-order polynomial fitting method. The random drift is shown in Figure 13. The angle random walk of gyro is identified by improved DAVAR method. The comparison of confidence level and dynamic tracking ability of the analysis results between improved DAVAR (Adaptive window length DAVAR, A-DAVAR), DAVAR-400 (window length is 400) and DAVAR-600 (window length is 600) is shown in Table 1. In Table 1, the beginning (3000 s) means that the window first senses the 3000 s data, and the window length should decrease at this time. The beginning (3000 s) should be close to DAVAR-600's beginning (3000 s). The end (3000 s) means that the window lost senses the 3000 s data, and the window length should increase at this time. The end (3000 s) should be close to DAVAR-400's end (3000 s). The ARW (0~2000 s) means the mean value from 0 s to 2000 s and the region is in stationary interval, so the ARW (0~2000 s) should be close to the refer value in Section 2.2.
From Table 1, we can see that the improved dynamic Allan variance solves the problem of both confidence and dynamic tracking ability.
The MEMS gyro is sensitive to temperature, so the MEMS gyro is fixed to the stationary preheating 1 h on the turntable to make the temperature reach steady state. The accurate angular velocity can be input to the MEMS gyro using the turntable. But the input angular velocity of the turntable is 0 in the experiment. Since the gyroscope bandwidth is 40 Hz, the sampling frequency is 100 Hz according to sampling theorem.
The gyro data collected for the 600 s is shown in Figure 13. When the turntable input angle speed is 0, the output signal of MEMS gyroscope contains trend and random drift. Usually the mean value method and polynomial fitting method are used to remove the trend terms. So the random drift of MEMS gyro is obtained by removing trend items from data collection with 3-order polynomial fitting method. The random drift is shown in Figure 13. The random drift of gyroscope is analyzed with Allan variance, and the Allan standard deviation curve is shown in Figure 14. The Figure 14 shows that the random error includes quantization error, zero bias instability, angular random walk and angular rate random walk. Fitting the Allan standard deviation curve by least square method, the error coefficient of angular random walk is 0.355°/h 1/2 . The random drift of gyroscope is analyzed with Allan variance, and the Allan standard deviation curve is shown in Figure 14. The Figure 14 shows that the random error includes quantization error, zero bias instability, angular random walk and angular rate random walk. Fitting the Allan standard deviation curve by least square method, the error coefficient of angular random walk is 0.355 • /h 1/2 . In the experiment, a slight disturbance is exerted to the gyroscope at 3000~4000 s. The sampling frequency is 100 Hz. The sampling time is 6000 s. According to a rule, one of 100 consecutive sampling data is selected as the research sample. And the gyro signal is shown in Figure 15. The Allan standard deviation curve and the least square fitting curve of the MEMS gyroscope are shown in Figure 16. We can get the error coefficient of random error by least square fitting. However, the analysis results can not reflect the dynamic changes of signals.  The rectangle window length is chosen as 600. The dynamic Allan variance of MEMS gyro is obtained as shown in Figure 17. In the experiment, a slight disturbance is exerted to the gyroscope at 3000~4000 s. The sampling frequency is 100 Hz. The sampling time is 6000 s. According to a rule, one of 100 consecutive sampling data is selected as the research sample. And the gyro signal is shown in Figure 15. The Allan standard deviation curve and the least square fitting curve of the MEMS gyroscope are shown in Figure 16. We can get the error coefficient of random error by least square fitting. However, the analysis results can not reflect the dynamic changes of signals. In the experiment, a slight disturbance is exerted to the gyroscope at 3000~4000 s. The sampling frequency is 100 Hz. The sampling time is 6000 s. According to a rule, one of 100 consecutive sampling data is selected as the research sample. And the gyro signal is shown in Figure 15. The Allan standard deviation curve and the least square fitting curve of the MEMS gyroscope are shown in Figure 16. We can get the error coefficient of random error by least square fitting. However, the analysis results can not reflect the dynamic changes of signals.   In the experiment, a slight disturbance is exerted to the gyroscope at 3000~4000 s. The sampling frequency is 100 Hz. The sampling time is 6000 s. According to a rule, one of 100 consecutive sampling data is selected as the research sample. And the gyro signal is shown in Figure 15. The Allan standard deviation curve and the least square fitting curve of the MEMS gyroscope are shown in Figure 16. We can get the error coefficient of random error by least square fitting. However, the analysis results can not reflect the dynamic changes of signals.  The rectangle window length is chosen as 600. The dynamic Allan variance of MEMS gyro is obtained as shown in Figure 17. The rectangle window length is chosen as 600. The dynamic Allan variance of MEMS gyro is obtained as shown in Figure 17. The quantization noise (QN, Q is a coefficient), angle random walk (ARW, A is a coefficient), bias instability (BI, B is a coefficient), rate random walk (RRW, R is a coefficient), rate ramp (RR, K is a coefficient) are obtained by fitting Allan standard deviation curve with least squares. And the random error is shown in Figure 18. We can see that the random error analysis results show strong dynamic characteristics. According to the improved DAVAR algorithm designed above, the sampled signal is analyzed. The kurtosis, sliding kurtosis contribution coefficient, and adaptive window length are shown in Figure 19. The quantization noise (QN, Q is a coefficient), angle random walk (ARW, A is a coefficient), bias instability (BI, B is a coefficient), rate random walk (RRW, R is a coefficient), rate ramp (RR, K is a coefficient) are obtained by fitting Allan standard deviation curve with least squares. And the random error is shown in Figure 18. We can see that the random error analysis results show strong dynamic characteristics. The quantization noise (QN, Q is a coefficient), angle random walk (ARW, A is a coefficient), bias instability (BI, B is a coefficient), rate random walk (RRW, R is a coefficient), rate ramp (RR, K is a coefficient) are obtained by fitting Allan standard deviation curve with least squares. And the random error is shown in Figure 18. We can see that the random error analysis results show strong dynamic characteristics. According to the improved DAVAR algorithm designed above, the sampled signal is analyzed. The kurtosis, sliding kurtosis contribution coefficient, and adaptive window length are shown in Figure 19. According to the improved DAVAR algorithm designed above, the sampled signal is analyzed. The kurtosis, sliding kurtosis contribution coefficient, and adaptive window length are shown in Figure 19. The quantization noise (QN, Q is a coefficient), angle random walk (ARW, A is a coefficient), bias instability (BI, B is a coefficient), rate random walk (RRW, R is a coefficient), rate ramp (RR, K is a coefficient) are obtained by fitting Allan standard deviation curve with least squares. And the random error is shown in Figure 18. We can see that the random error analysis results show strong dynamic characteristics. According to the improved DAVAR algorithm designed above, the sampled signal is analyzed. The kurtosis, sliding kurtosis contribution coefficient, and adaptive window length are shown in Figure 19.  Figure 19. The analysis result of kurtosis, sliding kurtosis contribution coefficient, and adaptive window length. Figure 19 shows that when the beginning of shock signal is collected by the L window at 3000 s, the length of the rectangular window begins to decrease gradually until the length of the window is reduced to a minimum L min . And when the beginning of the shock signal is collected by the end of the rectangular window, the length of the rectangular window increases gradually to the maximum length L max . In the interval of the end of the shock signal, there is a slight fluctuation in the signal for a short time due to temperature and electromagnetic interference which leads to reduce the length of the rectangular window in the improved DAVAR method. The sliding kurtosis contribution coefficient highlights the more sensitive ability, such as the peak in Figure 19. And the window length decreases in the growth interval after 4000 s, so the tracking ability of dynamic Allan variance has improved by the window length decreasing. The measured signals are analyzed with the improved DAVAR method, and the dynamic Allan variance analysis results of the measured signals are as shown in Figure 20.  Figure 19 shows that when the beginning of shock signal is collected by the L window at 3000 s, the length of the rectangular window begins to decrease gradually until the length of the window is reduced to a minimum Lmin. And when the beginning of the shock signal is collected by the end of the rectangular window, the length of the rectangular window increases gradually to the maximum length Lmax. In the interval of the end of the shock signal, there is a slight fluctuation in the signal for a short time due to temperature and electromagnetic interference which leads to reduce the length of the rectangular window in the improved DAVAR method. The sliding kurtosis contribution coefficient highlights the more sensitive ability, such as the peak in Figure 19. And the window length decreases in the growth interval after 4000 s, so the tracking ability of dynamic Allan variance has improved by the window length decreasing. The measured signals are analyzed with the improved DAVAR method, and the dynamic Allan variance analysis results of the measured signals are as shown in Figure 20. The feasibility of the improved DAVAR algorithm is analyzed by taking the angle random walk as the analysis example. The dynamic changes of the angular random walk error coefficient of the measured MEMS gyroscope signal are shown in Figure 21. We can find that the improved DAVAR algorithm with adaptive window length can be very sensitive to impact signal. And in the stationary interval (0~3000 s, 4000~6000 s), the improved DAVAR analysis results coincide with the DAVAR analysis results obtained from the rectangular window interception signal with a window length of 600. The results of the analysis have good confidence and stability. In the perturbation interval (3000~4000 s), the improved DAVAR analysis results coincide with the DAVAR analysis results obtained from the rectangular window interception signal with a window length of 400. The results of the analysis have good confidence and stability. And the analysis results have good dynamic tracking ability. So the feasibility of the proposed method is verified.   The feasibility of the improved DAVAR algorithm is analyzed by taking the angle random walk as the analysis example. The dynamic changes of the angular random walk error coefficient of the measured MEMS gyroscope signal are shown in Figure 21. We can find that the improved DAVAR algorithm with adaptive window length can be very sensitive to impact signal. And in the stationary interval (0~3000 s, 4000~6000 s), the improved DAVAR analysis results coincide with the DAVAR analysis results obtained from the rectangular window interception signal with a window length of 600. The results of the analysis have good confidence and stability. In the perturbation interval (3000~4000 s), the improved DAVAR analysis results coincide with the DAVAR analysis results obtained from the rectangular window interception signal with a window length of 400. The results of the analysis have good confidence and stability. And the analysis results have good dynamic tracking ability. So the feasibility of the proposed method is verified.
Micromachines 2018, 9, x 16 of 22 Figure 19 shows that when the beginning of shock signal is collected by the L window at 3000 s, the length of the rectangular window begins to decrease gradually until the length of the window is reduced to a minimum Lmin. And when the beginning of the shock signal is collected by the end of the rectangular window, the length of the rectangular window increases gradually to the maximum length Lmax. In the interval of the end of the shock signal, there is a slight fluctuation in the signal for a short time due to temperature and electromagnetic interference which leads to reduce the length of the rectangular window in the improved DAVAR method. The sliding kurtosis contribution coefficient highlights the more sensitive ability, such as the peak in Figure 19. And the window length decreases in the growth interval after 4000 s, so the tracking ability of dynamic Allan variance has improved by the window length decreasing. The measured signals are analyzed with the improved DAVAR method, and the dynamic Allan variance analysis results of the measured signals are as shown in Figure 20. The feasibility of the improved DAVAR algorithm is analyzed by taking the angle random walk as the analysis example. The dynamic changes of the angular random walk error coefficient of the measured MEMS gyroscope signal are shown in Figure 21. We can find that the improved DAVAR algorithm with adaptive window length can be very sensitive to impact signal. And in the stationary interval (0~3000 s, 4000~6000 s), the improved DAVAR analysis results coincide with the DAVAR analysis results obtained from the rectangular window interception signal with a window length of 600. The results of the analysis have good confidence and stability. In the perturbation interval (3000~4000 s), the improved DAVAR analysis results coincide with the DAVAR analysis results obtained from the rectangular window interception signal with a window length of 400. The results of the analysis have good confidence and stability. And the analysis results have good dynamic tracking ability. So the feasibility of the proposed method is verified. The error of gyro angle random walk is identified by improved DAVAR method. The comparison of confidence level and dynamic tracking ability of the analysis results between improved DAVAR (A-DAVAR), DAVAR-400, and DAVAR-600 is as shown in Table 2. The Refer value is obtained by the Allan variance in Figure 16.
In order to further verify the feasibility of the proposed method, change the time of mechanical interference. And the interference is introduced into the gyroscope test at 1000~2000 s and 4000~5000 s. Remove the trend items from the original output and the dynamic random drift is obtained, as shown in Figure 22. The error of gyro angle random walk is identified by improved DAVAR method. The comparison of confidence level and dynamic tracking ability of the analysis results between improved DAVAR (A-DAVAR), DAVAR-400, and DAVAR-600 is as shown in Table 2. The Refer value is obtained by the Allan variance in Figure 16.
In order to further verify the feasibility of the proposed method, change the time of mechanical interference. And the interference is introduced into the gyroscope test at 1000~2000 s and 4000~5000 s. Remove the trend items from the original output and the dynamic random drift is obtained, as shown in Figure 22. If the kurtosis threshold is 3.6, 3.45 and 3.3, the rectangular window with a window length of 600 and 400 is used to intercept signal for dynamic Allan variance analysis. The nodes obtained by the window length adjustment according to kurtosis are shown in Figure 23. As shown in Figure 23, the closer the kurtosis threshold is to 3, the more sensitive the adjustment will be. But if the kurtosis threshold is too close to 3, the judgement of the stationarity of the signal will be distorted according to kurtosis. For example, when the length of the window is 600 and the kurtosis threshold is 3.3, 735~1539 s is judged to be unstable, and there is a big misjudgement. Therefore, we need to combine the sliding kurtosis contribution coefficient and kurtosis together to determine the stable form of the signal.
We can see the dynamic tracking ability of analysis results is best when the stability threshold of kurtosis is 3.3 in Figure 23. So adjust the length of a rectangular window only based on kurtosis (brief account DAVAR-K-3.3) as the same as Wang's method [12,13]. The kurtosis and the window length are shown in Figure 24. In Figure 24, the stability threshold of kurtosis is 3.3 and the window length in the red circle is not perfect. The ideal situation is that the signal should be intercepted by the rectangle window with 600 window length from 0s to 700 s. At 700 s, the beginning of the rectangular window begins to intercept the 1000 s signal, so the length of the window begins to decrease and the beginning of the rectangular window has been at 1000 s. When the rectangle window moves to 800 s, the length of the window is minimized to 400, at this time the beginning of If the kurtosis threshold is 3.6, 3.45 and 3.3, the rectangular window with a window length of 600 and 400 is used to intercept signal for dynamic Allan variance analysis. The nodes obtained by the window length adjustment according to kurtosis are shown in Figure 23. As shown in Figure 23, the closer the kurtosis threshold is to 3, the more sensitive the adjustment will be. But if the kurtosis threshold is too close to 3, the judgement of the stationarity of the signal will be distorted according to kurtosis. For example, when the length of the window is 600 and the kurtosis threshold is 3.3, 735~1539 s is judged to be unstable, and there is a big misjudgement. Therefore, we need to combine the sliding kurtosis contribution coefficient and kurtosis together to determine the stable form of the signal.
We can see the dynamic tracking ability of analysis results is best when the stability threshold of kurtosis is 3.3 in Figure 23. So adjust the length of a rectangular window only based on kurtosis (brief account DAVAR-K-3.3) as the same as Wang's method [12,13]. The kurtosis and the window length are shown in Figure 24. In Figure 24, the stability threshold of kurtosis is 3.3 and the window length in the red circle is not perfect. The ideal situation is that the signal should be intercepted by the rectangle window with 600 window length from 0s to 700 s. At 700 s, the beginning of the rectangular window begins to intercept the 1000 s signal, so the length of the window begins to decrease and the beginning of the rectangular window has been at 1000 s. When the rectangle window moves to 800 s, the length of the window is minimized to 400, at this time the beginning of the rectangular window is still in 1000 s. Then use the rectangular window with a length of 400 to intercept the signal after 1000 s until the movement reaches 1200 s. At 1200 s, the end of the rectangular window begins to intercept the 1000 s signal, the rectangle window length begins to increase, and the end of the rectangular window has been at 1000 s until it moves to 1300 s. And the length of the window grows to the maximum at 1300 s, at this time the end of the rectangular window is still located in 1000 s. The signal should be intercepted by the rectangle window with a 600 window length from 1300 s to 1700 s. This is the most ideal case, with the strongest sensitivity and the highest confidence level at every moment. The process of changing the length of the window is the same as that of the process mentioned above. So the window length in the red circle region influences the analysis results because of the small stability threshold of kurtosis. the rectangular window is still in 1000 s. Then use the rectangular window with a length of 400 to intercept the signal after 1000 s until the movement reaches 1200 s. At 1200 s, the end of the rectangular window begins to intercept the 1000 s signal, the rectangle window length begins to increase, and the end of the rectangular window has been at 1000 s until it moves to 1300 s. And the length of the window grows to the maximum at 1300 s, at this time the end of the rectangular window is still located in 1000 s. The signal should be intercepted by the rectangle window with a 600 window length from 1300 s to 1700 s. This is the most ideal case, with the strongest sensitivity and the highest confidence level at every moment. The process of changing the length of the window is the same as that of the process mentioned above. So the window length in the red circle region influences the analysis results because of the small stability threshold of kurtosis.  The dynamic Allan variance analysis results with DAVAR-K-3.3 are shown in Figure 25. And the angle random walk obtained by the DAVAR-K-3.3 is shown in Figure 26.
Do not change the parameters in the algorithm and the sampled signals are analyzed with the improved DAVAR algorithm. The analysis results of kurtosis, sliding kurtosis contribution coefficient and the window length are shown in the Figure 27. In Figure 27, the window length has been achieved adjust adaptively and the window length adjustment is close to the ideal situation.
The improved dynamic Allan variance analysis of dynamic random drift signals is shown in the Figure 28. the rectangular window is still in 1000 s. Then use the rectangular window with a length of 400 to intercept the signal after 1000 s until the movement reaches 1200 s. At 1200 s, the end of the rectangular window begins to intercept the 1000 s signal, the rectangle window length begins to increase, and the end of the rectangular window has been at 1000 s until it moves to 1300 s. And the length of the window grows to the maximum at 1300 s, at this time the end of the rectangular window is still located in 1000 s. The signal should be intercepted by the rectangle window with a 600 window length from 1300 s to 1700 s. This is the most ideal case, with the strongest sensitivity and the highest confidence level at every moment. The process of changing the length of the window is the same as that of the process mentioned above. So the window length in the red circle region influences the analysis results because of the small stability threshold of kurtosis.  The dynamic Allan variance analysis results with DAVAR-K-3.3 are shown in Figure 25. And the angle random walk obtained by the DAVAR-K-3.3 is shown in Figure 26.
Do not change the parameters in the algorithm and the sampled signals are analyzed with the improved DAVAR algorithm. The analysis results of kurtosis, sliding kurtosis contribution coefficient and the window length are shown in the Figure 27. In Figure 27, the window length has been achieved adjust adaptively and the window length adjustment is close to the ideal situation.
The improved dynamic Allan variance analysis of dynamic random drift signals is shown in the Figure 28. The dynamic Allan variance analysis results with DAVAR-K-3.3 are shown in Figure 25. And the angle random walk obtained by the DAVAR-K-3.3 is shown in Figure 26.
Do not change the parameters in the algorithm and the sampled signals are analyzed with the improved DAVAR algorithm. The analysis results of kurtosis, sliding kurtosis contribution coefficient and the window length are shown in the Figure 27. In Figure 27, the window length has been achieved adjust adaptively and the window length adjustment is close to the ideal situation.
The improved dynamic Allan variance analysis of dynamic random drift signals is shown in the Figure 28.         Taking the angle random walk as an example, the validation of the proposed method is carried out, and the analysis result is shown in Figure 29. As shown in Figure 29, the analysis results of the improved DAVAR algorithm have stronger sensitivity than the DAVAR-600, and they are more stable than the results of DAVAR-400. Taking the angle random walk as an example, the validation of the proposed method is carried out, and the analysis result is shown in Figure 29. As shown in Figure 29, the analysis results of the improved DAVAR algorithm have stronger sensitivity than the DAVAR-600, and they are more stable than the results of DAVAR-400. The comparison of confidence level and dynamic tracking ability of the analysis results between improved DAVAR, DAVAR-400, DAVAR-600, and DAVAR-K-3.3 is shown in Table 3. In Table 3, the beginning (1000 s) means that the window first sense the 1000 s data, and the window length should decrease at this time. The beginning (1000 s) should be close to DAVAR-600's beginning (1000 s). The end (1000 s) means that the window lost sense the 1000 s data, and the window length should increase at this time. The end (1000 s) should be close to DAVAR-400's end (1000 s). The restore 600 (1000 s) means the window length restores to 600 after 1000 s, and the The comparison of confidence level and dynamic tracking ability of the analysis results between improved DAVAR, DAVAR-400, DAVAR-600, and DAVAR-K-3.3 is shown in Table 3. In Table 3, the beginning (1000 s) means that the window first sense the 1000 s data, and the window length should decrease at this time. The beginning (1000 s) should be close to DAVAR-600's beginning (1000 s). The end (1000 s) means that the window lost sense the 1000 s data, and the window length should increase at this time. The end (1000 s) should be close to DAVAR-400's end (1000 s). The restore 600 (1000 s) means the window length restores to 600 after 1000 s, and the window begins to intercept the signal in stable region after 1000 s, so the restore 600 (1000 s) should be close to the DAVAR-600's end (1000 s). Other nodes are similar to the above. ARW (2500~3500 s) means the mean value of ARW from 2500 s to 3500 s and, in the stationary interval, the ARW should be close to the refer value.
From Table 3 and the analysis above, we can see that the improved dynamic Allan variance solves the problem of both confidence and dynamic tracking ability and the analysis results are the closest to the ideal situation in A-DAVAR, DAVAR-400, DAVAR-600, and DAVAR-K-3.3.

Conclusions
This article elaborates the shortcomings of Allan variance and dynamic Allan variance by using these two methods to analyze the analog MEMS gyro signal. Because the dynamic Allan variance analysis using the rectangular window with a fixed window length to capture the original signal has the disadvantage of power leakage and insensitive tracking of the impact signal. In this paper, the DAVAR analysis of signal with a rectangular window with fixed window length is used to analyze the signal. When the window length is long, the dynamic tracking ability of the signal is poor, but the reliability is high. On the contrary, when the window length is small, the dynamic tracking ability is strong but the confidence level is poor. Using kurtosis to adjust the length of the rectangle window, the method can improve the dynamic tracking ability of the analysis results. But by comparing and analyzing the different kurtosis stability threshold, it is found that the more closer the kurtosis stability threshold is close to 3, the stronger the dynamic tracking ability of DAVAR will be. However, when the kurtosis stability threshold is too close to 3, there is a misjudgement in kurtosis analysis of signal stability, resulting in distortion of DAVAR analysis results. So the slipping kurtosis contribution coefficient is proposed to deal with the impact component at the end of sensitive signals. And the paper carry out the beginning and end interval data of the signal with an inverted mirror extension to increase the utilization of beginning and end data. Through the kurtosis and sliding kurtosis contribution coefficient jointly adjusting the window length of the rectangular window, the kurtosis stability threshold can be closer to 3, which improves the dynamic tracking ability of the analysis results and realizes the adaptive change of the length of the window. It solves the problem that the dynamic Allan variance tracking ability and confidence level are difficult to take into account, and also solves the problem of misjudgement in the stability analysis of kurtosis.