Mechanical and Electrical Noise in Sense Channel of MEMS Vibratory Gyroscopes

This paper presents a theoretical analysis of mechanical and electrical noise in the sense channel of micro-electromechanical systems (MEMS) vibratory gyroscopes. Closed-form expressions for the power spectral density (PSD) of the noise equivalent rate (NER) of gyroscopes in the open-loop and the force-rebalance operations are derived by using an averaged PSD model and an equivalent transfer function. The obtained expressions are verified through numerical simulations, demonstrating close agreements between the analytic and the numerical models. Based on the derived expressions for the PSD of the NER, the impacts of the modal frequency split, quality factor, and the gain of the feedback forcer, as well as the gain of the signal conditioning circuit, on the gyroscope noise characteristics are theoretically analyzed. In addition, the angle random walk (ARW) and the standard deviation of the NER are also discussed through the PSD models. Finally, the effects of the loop closing, the mode matching, and the gain of the feedback forcer on the PSD of the NER were verified via a MEMS vibratory gyroscope with a tunable modal frequency split.


Introduction
Gyroscopes based on micro-electromechanical systems (MEMS) technologies have great advantages of small size, low power consumption, low cost, and batch fabrication over their traditional counterparts. MEMS vibratory gyroscopes have been widely adopted in various consumer electronics, as well as in industrial fields. However, the noise performance is one of the essential limitations of MEMS vibratory gyroscopes operating on high-precision occasions. Noise analysis of the gyroscope sense channel is of great importance to system designs. The impacts of different operations and system parameters on the noise performance can provide constructive guidance on determining a sufficient quality factor of the resonator, a reasonable modal frequency split, and adopting proper feedback forcers.
There have been several classic studies focusing on the noise analysis of MEMS vibratory gyroscopes. A comprehensive review of noise in MEMS devices is given in [1]. Noise in MEMS vibratory gyroscopes typically contains mechanical-thermal noise and electrical noise. The mechanical noise has a fundamental intrinsic mechanism, namely Brownian motion, which is due to dynamic unbalanced forces caused by random impacts of molecules on a structure. Leland investigated the mechanical-thermal noise in MEMS gyroscopes via the method of stochastic averaging [2]. The expressions for the impacts of the mechanical noise, including the angle random walk (ARW), the noise equivalent rate (NER), and the spectral density of the rate noise, are derived in [2]. Common electrical noise sources are summarized in [3], including electrical-thermal noise (Johnson noise), Besides noise, bias drift is another important performance of gyroscopes. The bias drift is a result of the unstable behavior of the zero rate output of the gyroscopes due to environmental changes, especially to the change of temperature. The zero rate output of MEMS vibratory gyroscopes is studied in detail in [16]. Different operations and parameters of gyroscopes will change the zero rate output and, therefore, will certainly change the bias drift. However, considering the bias drift and the noise have different formation mechanisms, we will not discuss the effects of different operations and parameters on the bias drift in this paper.
This paper is organized as follows. Section 2 gives a general background of the sense channel model in MEMS vibratory gyroscopes, including the open-loop operation and the force-rebalance operation. The expressions for the frequency response of the gyroscope and the scale-factor are briefly introduced in both operations. The open-loop noise model is analyzed in Section 3. The noise propagation in the open-loop is demonstrated, and the closed-form expression for the PSD of the NER are derived. Section 4 presents the noise analysis of the force-rebalance model. The validity of the obtained closed-form expressions is verified through numerical simulations for both the mode-split and mode-matched gyroscopes. Through the obtained PSD models, the ARW and the standard deviation of the NER are discussed in Section 5. Some conclusions in Section 5 were verified via a MEMS vibratory gyroscope with a tunable frequency split, and the results are shown in Section 6. Finally, Section 7 concludes this paper with a summary.

Model of Sense Channel in MEMS Vibratory Gyroscopes
Typically a MEMS vibratory gyroscope can be modeled as a two-dimensional oscillator where m is the mass of the oscillator; x and y are displacements along the two orthogonal dimensions, that is x-and y-axes; η is an angular gain factor, which is determined by the structural geometry of the gyroscope; Ω is the input angular rate; c x and c y represent the damping coefficients; k x and k y are stiffness; and F x and F y are external forces exerted on the x-and y-axes. ω x = √ k x /m and ω y = k y /m are natural frequencies of the x-axis and the y-axis, respectively. If ω x = ω y , gyroscopes are in the mode-matched situation; otherwise, gyroscopes are in the mode-split situation. The damping and the stiffness couplings are ignored in Equations (1) and (2) because these non-ideal components only introduce output bias and, accordingly, have no significant impacts on the gyroscope noise performance.
We refer to the x-axis of vibratory gyroscopes as the drive axis and the y-axis as the sense axis. The drive axis is excited into stable resonance by an external electrostatic force, F x , whose frequency and amplitude are respectively regulated by a phase-locked loop (PLL) and an automatic gain control (AGC) loop. In the open-loop operation, there is no external force applied along the sense axis, which means F y = 0, and the sense axis is only driven by the Coriolis force. The vibration amplitude of the sense axis is proportional to the angular rate input and, therefore, can be used as a measure of the angular rate. In the force-rebalance operation, namely the closed-loop operation, an electrostatic force is applied to exactly balance the Coriolis force and to null the vibration of the sense axis. In this operation, the voltages applied to the feedback electrodes, which generate the electrostatic force, can be used as an estimation of the angular rate.
We assume the drive axis has a steady vibration under the control of the PLL and the AGC loop where A x is the vibration amplitude of the drive axis in meters, and ω d is the drive frequency and always tracks ω x . Equation (3) indicates that the Coriolis force exerted on the sense axis is We refer to the sense axis, together with the corresponding readout and control circuits, as the sense channel. Figure 1 demonstrates the general block diagrams of the sense channel in both the open-loop operation and the force-rebalance operation. Symbols in Figure 1 are explained as follows. G y (s) is the transfer function of the sense axis from applied forces to pickoff voltages where K xv is the pickoff circuit gain from displacement to voltage in V/m; and Q y = mω y /c y is the quality factor. K a is the gain of the signal conditioning circuit in V/V. L(s) is the transfer function of a low pass filter (LPF) with a cutoff frequency of f l0 . α o and α c are scale-factors of the open-loop operation and the force-rebalance operation, respectively. Module PI in Figure 1b is the feedback controller, which is commonly implemented by a proportional-integral controller. The gain of the feedback forcer, K f , in N/V, is given by where C f b is the capacitance of the forcer electrodes, K DA is the gain of the buffer of the digital-to-analog converter, and V D is the bias direct current (DC) voltage applied to the forcer electrodes. ϕ d is a demodulation phase, whose optimal value is determined by In practice, ϕ d will be set as 0 for mode-split gyroscopes and −π/2 for mode-matched gyroscopes. Considering the modulation-demodulation process shown in Figure 2, the equivalent transfer function can be expressed as [17] where the approximation is taken when Q y >> 1/2 and ω y,d >> 2π f l0 .
and the open-loop scale-factor as where ∆ω ω y − ω d . Similarly, the block diagram of the force-rebalance operation in Figure 1b can be equivalent to the one shown in Figure 3. For convenience, we denote where K p and K i are parameters of the PI controller. Considering the overall transfer function of the force-rebalance sense channel is and the force-rebalance scale-factor is Equations (9) and (13) respectively give the gyroscope frequency responses in the open-loop operation and the force-rebalance operation to the applied angular rate.

Open-Loop Noise Model
To evaluate the noise characteristics in the open-loop sense channel, we ignore the angular rate input and introduce the mechanical and the electrical noise, as depicted in Figure 4. N m (t) is mechanical-thermal noise in Newton. N ep (t) is the electrical noise brought in by the pickoff and the signal conditioning circuits, including the flicker noise and the Johnson noise. Ω N (t) is the measured NER in the gyroscope output. (1) The vibration of the drive axis under the electrostatic force excitation is sufficiently greater than the mechanical-thermal disturbance. Therefore, the NER coupling from the amplitude noise of the drive axis through Coriolis effect is insignificant and can be ignored. (2) Considering the drift of ω x is extremely slow, the control parameters of the PLL in the drive axis can be set to achieve a very low phase noise for the frequency tracking. Consequently, the phase noise in the demodulation process is neglected in this paper. (3) Generally, there are multiple electrical noise sources in the system. However, in linear systems, multi stage noise sources can be equivalently converted into the one at the first stage [18]. In analog systems, the noise sources after the demodulation, typically introduced by the LPF, are trivial and thus can be neglected. In digital systems, no additional electrical noise will be presented after the analog-to-digital conversion. Taking these facts into account, we only model one electrical noise source at the pickoff node. (4) In addition, digitally-based noise, namely spurious, which is usually introduced by clocks and power supplies, is also ignored in this paper.
The PSD of the mechanical noise and the electrical noise can be respectively described by where f is the frequency in Hertz, k B is the Boltzmann constant, T is the temperature in Kelvin, R p represents an equivalent resistance for the Johnson noise, and k 1 is a constant for the flicker noise. Double-side spectral models are adopted in Equations (15) and (16) for the convenience of theoretical analysis. It should be noted that the equivalent resistance R p is not the resistance of a real resistor but a variable representing the Johnson noise equivalently. The noise performance of the pickoff circuit can be evaluated by the ratio of S N ep ( f )/K xv , which is always preferred as low as possible.
The PSD of N (t) in Figure 4 can be readily obtained as where To scrutinize the PSD of the modulated noise N c (t) = N (t) · cos (ω d t + ϕ d ), we examine its autocorrelation function (ACF) via where E[·] denotes the expectation operator, and R N (τ) represents the ACF of N (t). Equation (19) states that N c (t) is not wide-sense stationary due to its time dependence. However, because R N c (τ, t) is a fast-changing variable with respect to time, we can find an averaged ACF of N c (t) as over a period of π/ω d . Accordingly, applying Fourier transform to Equation (20) and denoting where The approximations in Equations (22) and (23) hold under the conditions that ∆ω << ω y and 2π f << ω y . Figure 5 visualizes the noise PSD described by Equations (15)- (17) and (21). For the mechanical noise in N (t), the PSD at f y = ω y /2π is larger than those in low frequencies by a factor of Q 2 y , predicted by Equation (18). The flicker noise is pronounced in low frequencies, and the Johnson noise prevails in high frequencies. The crossover frequency, at which the PSD of the flicker noise equals that of the Johnson noise, is generally from several Hertz to tens of Hertz. As depicted in Figure 5c, the PSD feature of N (t) can be utilized to identify the mechanical noise and the electrical noise for situations where the Johnson noise PSD is below the peak of the mechanical noise PSD. Once the peak magnitude of the mechanical noise PSD at the resonant frequency is obtained, the overall PSD of the mechanical noise can be determined via Equation (18). An example of this procedure will be demonstrated in Section 6. Considering the low-pass feature of L(s), the modulated flicker noise can be sufficiently suppressed in the gyroscope output, yielding the PSD of the NER as where L( f ) = L(s)| s=j2π f . For mode-split gyroscopes, provided that ∆ω >> ω y /2Q y , the open-loop scale-factor (Equation (10)) can be simplified as where the modulation phase ϕ d is taken as 0. By combining Equations (24) and (25), we obtain the expression of the PSD of the NER for open-loop mode-split gyroscopes as (26) Figure 6 shows the influences of the modal frequency split ∆ f and the quality factor Q y on the NER. The negative frequency parts of the plots are not shown here for the even property of S Ω N ( f ). The default parameter values used in this paper are listed in Table 1 unless some of them are re-specified in the context or figures. In Figure 6a, the PSD of the mechanical NER is not notably affected by ∆ f in low frequencies, but its peak shifts with ∆ f ; the PSD of the electrical NER increases as ∆ f grows. The PSD of the mechanical NER is inversely proportional to Q y at low frequencies but proportional to Q y at the frequency of ∆ f , Figure 6b. The PSD of the electrical NER will not be affected by Q y , revealed by Equation (26), and hence their relationship is not presented in Figure 6b for the sake of clarity.

Parameters Values Parameters Values
For mode-matched gyroscopes, ∆ω = 0, hence ϕ d should be set as −π/2, and the open-loop scale-factor (Equation (10)) can be simplified as Equations (24) and (27) give the PSD of the NER for open-loop mode-matched gyroscopes The relationships between the NER and Q y are shown in Figure 7, which clearly demonstrates that the boost of Q y can reduce both the mechanical NER and the electrical NER. High Q y improves the electrical NER in open-loop mode-matched gyroscopes because the sensitivity of these gyroscopes is directly proportional to Q y , Equation

Force-Rebalance Noise Model
In the force-rebalance sense channel, besides the noise sources introduced in the open-loop, we model another electrical noise source N e f (t) at the feedback node, as demonstrated in Figure 8. The reason is that N e f (t) is separated from N ep (t) by the modulation-demodulation process in the loop and, as a result, cannot be equivalent to N ep (t). The electrical noise source at the output stage is still ignored for the same reason discussed in Section 3. Different from the open-loop operation, noise sources of the force-rebalance operation are in the closed-loop, making the noise propagate circularly. The basic idea is to isolate the noise sources from the closed-loop to make the analysis straightforward. Consider the system dynamics of the model demonstrated in Figure 8 f where * denotes the convolve operation, and g y (t) and c(t) are inverse Laplace transforms of G y (s) and C(s), respectively. Through simple algebraic manipulations of above equations, we get which can be translated into the model presented in Figure 9. Via this approach, the noise sources are all separated from the closed-loop. Figure 9. The equivalent noise model of the force-rebalance sense channel.
By applying the linearization method presented in Figure 2 to the model presented in Figure 9, we can easily obtain In Figure 9, the noise propagation in the force-rebalance operation is close to that in the open-loop operation S where and Combining Equations (15), (16) and (36)-(39) yields where the modulated flicker noise is ignored in consideration of the low-pass feature of T N (s).
It should be noted that thanks to the modulation and the followed filtering process, the flicker noise in the loop, either open-loop or closed-loop, does not contribute to the NER in the output of gyroscopes. In practice, the NER that resulted from the flicker noise is believed to be introduced by the electrical noise source at the output stage.
Equations (12) and (35) share a common denominator, suggesting that the PSD of the NER and the frequency response of the gyroscope have magnitude peaks at the same frequencies.
In mode-split cases, Equation (8) can be further simplified to H(s) ≈ K xv 4mω y ∆ω s 2 + ω y s/Q y + ∆ω 2 . (42) In Figure 10, the numerically simulated PSD curves are acquired via the noise model in Figure 8, while the analytic PSD curves are obtained through the closed-form Equation (41), along with Equations (14), (22), (23), (35) and (42). The consistency between the numerical and analytic curves manifests that the closed-form expressions can effectively describe the force-rebalance noise model in Figure 8. The mechanical NER and the feedback electrical NER are both determined by T N (s) and G y (s), leading to similar PSD shapes in Figure 10a,c. By contrast, the pickoff electrical NER is determined by T N (s) only.  Figure 11 illustrates the impacts of ∆ f , Q y , K a , and K f on the NER. Similar to the relationships exhibited in Figure 6, in low frequencies, the PSD of the mechanical NER is not influenced by the modal frequency split, and the PSD of the pickoff electrical NER increases as the split increases, Figure 11a. The PSD of the mechanical NER is inversely proportional to Q y while the PSD of the pickoff electrical NER is not affected by Q y , Figure 11b. The PSD shapes of the feedback electrical NER and those of the mechanical NER are very close except that (1) the PSD of the feedback electrical NER does not depend on Q y in low frequencies whereas that of the mechanical NER does, Figure 11b, and (2) the PSD of the feedback electrical NER will be changed by K f whereas that of the mechanical NER is relatively irrelevant to K f in low frequencies, Figure 11d. Here, we pay more attention to the PSD values in low frequencies because those values are directly related to the ARW of gyroscopes, which we will discuss in Section 5.
In mode-matched cases, Equation (8) will degenerate to Equations (14), (22), (23), (35), (41) and (43) give the closed-form expressions of the PSD of the NER in force-rebalance mode-matched gyroscopes. Figure 12 presents the corresponding analytic and numerically simulated PSD curves. In this operation, the controller parameters are set as K p = 1 and K i = 120. The analytic curves are in close agreement with the numerical ones, indicating the validity of the equivalent noise model.  Comparisons of the PSD of the NER between different Q y , K a , and K f are demonstrated in Figure 13. For mode-matched gyroscopes working in the force-rebalance operation, the enhancement of Q y can significantly improve the mechanical and the pickoff electrical noise. Similar to mode-split gyroscopes in the force-rebalance operation, in low frequencies, K f only affects the PSD of the feedback electrical NER. From Figures 11c and 13b, the gain of the signal conditioning circuit will not influence the PSD of the NER in low frequencies for both the mode-split and the mode-matched cases.

Angle Random Walk and Standard Deviation of Noise Equivalent Rate
The random noise of gyroscopes is typically characterized by ARW, which can be derived from the Allan Variance of the gyroscope output data, at an integration time of 1 s. The ARW is related to the PSD of the NER via [2,8] According to the PSD of the NER discussed in Sections 3 and 4, the values of S Ω N (0) for gyroscopes in different operations are summarized in Table 2.  (1) The increase in A x can improve the ARW attributable to all noise components, which suggests an effective approach to improve the noise performance. (2) The boost of Q y can reduce the ARW resulted from the mechanical noise in all of the operation cases but only decreases the ARW, to which the pickoff electrical noise contributes, in mode-matched cases. (3) Mode matching will not influence the ARW from the mechanical noise but will affect the ARW from the pickoff electrical noise by a factor of ∆ω/(ω y /2Q y ).
(4) The enhancement of the pickoff circuit, namely a decrease in the ratio of S N ep ( f )/K xv , can reduce the ARW introduced by the pickoff electrical noise but does not affect the mechanically induced ARW. (5) The gain of the feedback forcer only affects the ARW from the feedback electrical noise. In other words, in a well-designed system, if the feedback electrical noise is insignificant, the design value of K f is non-essential for the noise consideration.
(6) For a given gyroscope, loop closing of the sense channel, namely changing the sense channel from open-loop to force-rebalance, will not change the ARW from the mechanical noise and the pickoff electrical noise, neither in the mode-split case nor the mode-matched case.
In some applications, the standard deviation of the gyroscope output, σ Ω , in a specified frequency band, is utilized to specify random drift [7]. The standard deviation of the NER, gives the lower limit of σ Ω , assuming that the output bias does not drift. B in Equation (45) is the frequency band of interest and is generally comparable to the gyroscope bandwidth. Signals with frequencies far beyond the gyroscope bandwidth do not contain effective rate information and, consequently, can be smoothed or averaged. Based on Equation (45), the calculation of the standard deviation of the NER is very straightforward in open-loop sense channels. In mode-split cases, and in mode-matched cases, where the LPF is assumed as Considering the complexity of the expressions in force-rebalance channels, the corresponding σ Ω N can be discretely estimated via where f res is the frequency resolution.

Experiments and Results
The noise characteristics of MEMS vibratory gyroscopes were evaluated by a silicon tuning folk micro-gyroscope with a tunable modal frequency split, previously reported in [15]. Due to the negative stiffness effect, the resonant frequency of the sense axis can be lowered by DC voltages over the frequency tuning electrodes shown in Figure 14, which presents the schematic of the gyroscope design. For the gyroscope under the test, the natural frequency of the drive axis was 3563 Hz with Q x = 2284, and the natural frequency of the sense axis was 3629 Hz with Q y = 187, indicating a modal frequency split of 66 Hz. The modal frequency split could be tuned to 0 with a DC voltage of 8.2 V. The control system of the gyroscope was implemented by a field-programmable gate array (FPGA).  Firstly, the method to identify the mechanical noise and the electrical noise mentioned in Section 3 was verified through the primary resonator (drive axis) of the gyroscope. The output of the pickoff circuit was recorded by a dynamic signal analyzer at a sampling frequency of 102.4 kHz for 100 s when the resonator was kept stationary. The resultant frequency resolution was 0.01 Hz. Figure 15 demonstrates the PSD of the acquired signal, which is similar to the right half PSD shown in Figure 5c. The crossover frequency of the electrical noise was measured as 12 Hz, and the PSD of the electrical-thermal noise was 7.793 × 10 −13 V 2 /Hz. A magnitude peak at the natural frequency of the resonator could be clearly observed, suggesting the mechanical noise PSD at the natural frequency was 6.968 × 10 −12 V 2 /Hz. According to Equation (18), the PSD of the mechanical noise in low frequencies could be estimated as 1.336 × 10 −18 V 2 /Hz, which was much smaller than that of the electrical-thermal noise. By this approach, we can separate the mechanical noise from the very noisy electrical noise, as long as the magnitude peak of the mechanical noise PSD can be notably observed.
Then, the gyroscope was configured to operate in the open-loop mode-split case. Figure 16a shows the corresponding spectral density of the NER, along with the measured frequency response of the gyroscope. The output rate of the gyroscope was 1 kHz and the recording time was 100 s. The frequency response of the gyroscope was measured by an angular vibration table from 5 Hz to 80 Hz. The spectral density of the NER and the gyroscope frequency response had magnitude peaks at the same frequency, as explained in Section 4. In Figure 16a, S Ω N (0) ≈ 1.307 × 10 −3 ( • /s)/ √ Hz, which was an average value between 0.01 Hz to 1 Hz. The approximated ARW was consistent with the Allan Variance in Figure 16b. To calculate the Allan Variance, the gyroscope output rate was averaged to 1 Hz, and the data was recorded for 5400 s.  Next, a DC voltage of 5 V was applied to the frequency tuning electrodes, resulting in a decrease in ∆ f from 66 Hz to 43 Hz. Figure 17 demonstrates the NER spectral density of the open-loop gyroscope and the measured gyroscope frequency response. By comparing Figures 16a and 17, it can be found that the open-loop scale-factor increased by 1.546 times as the frequency split decreased by a factor of 1.535, in agreement with Equation (25). However, the ARW only improved by 1.19 times, suggesting that for the gyroscope under test, the ARW was significantly affected by both the mechanical noise and the electrical noise. Finally, Figure 18 presents the NER spectral density of the gyroscopes in the force-rebalance operation. Revealed from Figures 16a and 18a, the loop closing of the sense channel did not improve the ARW, as predicted in Table 2. Figure 18a,b manifest that the mode matching did not change the scale-factor in the force-rebalance operation but improved the ARW by 1.52 times. The unremarkable improvement of the ARW was due to the considerable mechanical noise. Figure 18c was obtained by reducing V D in Equation (6) from 5 V to 2.5 V, resulting in a half K f . The corresponding scale-factor was increased by a factor of 2, which was in agreement with Equation (14). The ARW was not significantly affected by the change of K f . From the discussion in Section 5, K f would not influence the ARW attributable to the mechanical noise or the pickoff electrical noise. If the impact of the feedback electrical noise is negligible, the design of the feedback forcer is irrelevant to the noise performance. In summary, in both the open-loop operation and the force-rebalance operation, the decrease of the modal frequency split improved the ARW due to the decrease of the electrical noise, as presented in Table 2. However, the mechanical noise was not affected by the modal frequency split, resulting in a moderate improvement of the ARW through mode matching. To further improve the ARW of the gyroscopes under the mode-matched situations, the boost of quality factor is essential. The closing of the sense loop did not significantly change the ARW according to the experimental data, showing an agreement with the models demonstrated in previous sections. The gain of the feedback forcer did not distinctly affect the overall ARW of the gyroscopes, as discussed in Section 5. This result can give very important guidance on the structural design of the feedback forcers, implying that the gain of the forcers is irrelevant to the noise consideration.

Conclusions
In this paper, we have presented a theoretical analysis of mechanical and electrical noise in the sense channel of MEMS vibratory gyroscopes. The closed-form expressions for the PSD of the NER of gyroscopes in the open-loop and the force-rebalance operations have been obtained. In the open-loop noise model, the features of the noise after the demodulation are investigated by the averaged PSD model. In the force-rebalance noise model, the noise sources are separated from the closed-loop, and the modulation-demodulation process in the closed-loop is linearized by the equivalent transfer function. Via this approach, the closed-loop noise model can be equivalent to an open-loop one and, hence, the analysis method used in the open-loop model can be applied. The analytic models are in close agreement with the numerical models, verifying the validity of the closed-form expressions. The influences of the structure parameters, such as modal frequency split, quality factor, and the gain of the feedback forcer, as well as the gain of the signal conditioning circuit, on the gyroscope noise features have been theoretically analyzed by using the PSD model. In addition, the ARW and the standard deviation of the NER also have been discussed. The impacts of the modal frequency split, the loop closing, the mode matching, and the gain of the feedback forcer on the gyroscope noise characteristics were verified through a MEMS vibratory gyroscope with a tunable modal frequency split.
In both the open-loop and the force-rebalance operations, the decrease of the modal frequency split improved the ARW due to the decrease of the electrical noise. However, the mechanical noise was not affected by the modal frequency split. To further improve the ARW of the gyroscopes under the mode-matched situations, it is important to boost the quality factor. The loop closing did not significantly change the ARW according to the experimental data, which is in agreement with the models proposed in this paper. The gain of the feedback forcer will not distinctly affect the ARW as long as the electrical noise at the feedback node is negligible.
Author Contributions: X. Ding and H. Li analyzed the theoretical model; X. Ding and Y. Gao designed the experiments; J. Jia performed the experiments; X. Ding and J. Jia analyzed the data; X. Ding wrote the paper.

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations
The following abbreviations are used in this manuscript:  σ Ω N standard deviation of NER, rad/s ϕ d demodulation phase, rad ϕ od optimal demodulation phase, rad ω d vibration frequency of the x-axis, rad/s ω x natural frequency of the x-axis, rad/s ω y natural frequency of the y-axis, rad/s Ω(s) Laplace transform of Ω(t) Ω(s) Laplace transform of Ω(t) Ω(t) input angular rate, rad/s Ω(t) measured angular rate, rad/s Ω ARW angle random walk, (rad/s)/ √ Hz