Bandwidth extension for continuous mode reversal gyroscope

Abstract

A vibratory gyroscope system having a mechanical resonator (proof mass) and drive circuitry for maintaining oscillation in two axes at a small frequency split. Angular rate input is shifted to the frequency split and modulates both the frequency (FM) and the amplitude (AM) of the oscillations. Unlike other gyroscope modulation techniques which derive rate information from only the FM information and are subject to aliasing for rate signals with bandwidth exceeding the modulation frequency, the innovation uses both the FM and AM information. By exploiting their orthogonality, the image frequencies from the modulation cancel, thus removing the bandwidth limitation.

Claims

1. An electromechanical system, comprising: a mechanical resonator having a first mode of vibration and an associated first natural frequency, and a second mode of vibration having an associated second natural frequency, wherein angular rate of motion input couples energy between said first mode of vibration and said second mode of vibration; sensors and actuators for each of said first mode and said second mode for transducing an electrical signal into a mechanical vibration and transducing a mechanical vibration into an electrical signal; sustaining circuitry connected to the actuators to maintain vibrations in the first mode at a first frequency and the second mode at the second frequency; an amplitude and frequency measurement circuit configured for determining the instantaneous amplitude and frequency variation and outputting modulated versions of mode vibration amplitude and frequency as a representation of vibration amplitude and frequency; a phase determination circuit configured for receiving mode vibration signals and generating rate reference signals, and outputting demodulated frequency and amplitude rate signals; and a processing element configured for summing demodulated frequency and amplitude rate signals to generate an output rate in which the image spectra of the demodulation are eliminated from the angular rate determination, by exploiting the opposite polarity of the image spectra during summing.

2. The electromechanical system as recited in claim 1, wherein said electromechanical system comprises a continuous mode reversal gyroscope.

3. The electromechanical system as recited in claim 2, wherein said gyroscope is configured with a proof mass configured to move along at least two orthogonal axes.

4. The electromechanical system as recited in claim 2, wherein said proof mass is suspended by springs.

5. The electromechanical system as recited in claim 1, wherein said mechanical resonator comprises a gyroscope transducer.

6. The electromechanical system as recited in claim 1, wherein said mechanical resonator is symmetric with respect to Coriolis coupled vibration modes.

7. The electromechanical system as recited in claim 1, wherein said first natural frequency and said second natural frequency of said mechanical resonator differs by a range of approximately 1 to 100 Hz.

8. The electromechanical system as recited in claim 7, wherein said mechanical resonator is manufactured or trimmed to achieve said difference between said first natural frequency and said second natural frequency, or a phase frequency control circuit is utilized which is configured to determine these natural frequencies and make adjustments to operation the mechanical resonator.

9. The electromechanical system as recited in claim 1, wherein said system is configured to infer an angular rate of motion from the gains of the sustaining circuitry, the frequency of vibration of the first mode, and the frequency of vibration of the second mode.

10. The electromechanical system as recited in claim 1, wherein said system is configured to infer an angular rate of motion from the sum of the gains of the sustaining circuitry and the frequencies of vibration of the first mode and the frequency of vibration of the second mode.

11. The electromechanical system as recited in claim 1, wherein said system is configured to infer an estimate of the errors of the inferred angular rate from gains of the sustaining circuitry and the frequencies of vibration of the first mode and the second mode, and subtracting this estimate from the inferred output rate.

12. The electromechanical system as recited in claim 11, further comprising a circuit configured for estimating temperature of the mechanical resonator from an estimate of the errors of the inferred angular rate.

13. The electromechanical system as recited in claim 12, further comprising a circuit configured to compensate for errors due to variations in the temperature of the mechanical resonator using the temperature estimate.

14. An electromechanical system, comprising: a mechanical resonator having a first mode of vibration and an associated first natural frequency, and a second mode of vibration having an associated second natural frequency, wherein angular rate of motion input couples energy between said first mode of vibration and said second mode of vibration; sensors and actuators for each of the first mode and second mode of vibration for transducing an electrical signal into a mechanical vibration and transducing a mechanical vibration into an electrical signal; sustaining circuitry with variable gains connected to the actuators to maintain vibrations in the first mode and the second mode with gains adjusted to maintain substantially constant, non-zero velocity amplitude vibrations in the first mode at the first frequency of vibration and the second mode at the second frequency of vibration; an amplitude and frequency measurement circuit configured for determining the instantaneous amplitude and frequency variation and outputting modulated versions of mode vibration amplitude and frequency as a representation of vibration amplitude and frequency; a phase determination circuit configured for receiving mode vibration signals and generating rate reference signals, and outputting demodulated frequency and amplitude rate signals; and output circuitry to infer an angular rate of motion in response to summing demodulated frequency and amplitude rate signals to generate an output rate in which the image spectra of the demodulation are eliminated from the angular rate determination, by exploiting the opposite polarity of the image spectra during summing.

15. The electromechanical system as recited in claim 14, wherein said electromechanical system comprises a continuous mode reversal gyroscope.

16. The electromechanical system as recited in claim 15, wherein said gyroscope is configured with a proof mass configured to move along at least two orthogonal axes.

17. The electromechanical system as recited in claim 15, wherein said proof mass is suspended by springs.

18. The electromechanical system as recited in claim 14, wherein said mechanical resonator comprises a gyroscope transducer.

19. The electromechanical system as recited in claim 14, wherein said mechanical resonator is symmetric with respect to Coriolis coupled vibration modes.

20. The electromechanical system as recited in claim 14, wherein said first natural frequency and said second natural frequency differ by a range of approximately 1 to 100 Hz.

21. The electromechanical system as recited in claim 20, wherein said mechanical resonator is manufactured or trimmed to achieve said difference between said first natural frequency and said second natural frequency, or a phase frequency control circuit is utilized which is configured to determined these natural frequencies and make adjustments to operation the mechanical resonator.

22. The electromechanical system as recited in claim 14, wherein said output circuitry is configured to infer an angular rate of motion from the gains of the sustaining circuitry, and the frequency of vibration of the first mode and the second mode.

23. The electromechanical system as recited in claim 14, wherein said output circuitry is configured to infer an angular rate of motion from a sum of the gains of the sustaining circuitry and the frequencies of vibration of the first mode and second mode.

24. The electromechanical system as recited in claim 23, wherein said output circuitry is configured to infer an estimate of the errors of the inferred angular rate from the gains of the sustaining circuitry and the frequencies of vibration of the first and second mode and subtracting this estimate from the inferred angular rate.

25. The electromechanical system as recited in claim 24, further comprising a circuit configured for estimating the temperature of the electromechanical system from the estimate of the errors of inferred angular rate.

26. The electromechanical system as recited in claim 25, further comprising a circuit configured for compensating for errors due to variations of the temperature of the electromechanical system using the temperature estimate.

Description

BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWING(S)

(1) The technology described herein will be more fully understood by reference to the following drawings which are for illustrative purposes only:

(2) FIG. 1 is a block diagram of a gyroscope having the bandwidth extension circuitry according to the present disclosure.

(3) FIG. 2A and FIG. 2B are plots of intermediary gyroscope signals for a constant input rate as utilized according to an embodiment of the present disclosure.

(4) FIG. 3 is a schematic of a frequency demodulator utilized according to at least one embodiment of the present disclosure.

(5) FIG. 4 is a block diagram of a phase extractor utilized according to at least one embodiment of the present disclosure.

(6) FIG. 5 is a spectrum plot of gyroscope rate input Ω.sub.z, modulated and demodulated rate spectra, and measured rate output, as found according to an embodiment of the present disclosure.

(7) FIG. 6 is a spectrum plot of gyroscope rate input Ω.sub.z with bandwidth greater than the frequency split on the x and y modes, as utilized according to an embodiment of the present disclosure.

(8) FIG. 7 is a spectrum plot of FIG. 6, shown here with additional noise spectra representing errors in practical gyroscope realizations.

DETAILED DESCRIPTION

(9) The FM gyroscope described in this disclosure utilizes a proof mass suspended by springs (of any desired type/implementation), with the proof mass being free to move along at least two orthogonal mechanical modes, as described in applicant's previous work in U.S. Pat. Nos. 9,869,553 and 9,846,055, both of which are incorporated herein by reference in their entireties. These two orthogonal mechanical modes are excited to vibrate simultaneously and are referred to in the description below as x and y modes.

(10) FIG. 1 illustrates an example embodiment 10 of a direct frequency readout gyroscope system according to the presented technology. The gyroscope system 10 is shown with a mechanical resonator element 12 having sensing means 14a and 14b and actuating means 16a and 16b. Sustaining circuitry 20 is shown to mechanically excite the resonator at or near its natural frequencies and overcome resonator damping. Output circuitry 24 converts raw data from the gyroscope to a rate signal 26 to be utilized by the user application circuit.

(11) It should be appreciated that in at least one embodiment of output circuitry 24 is configured 24′ where signals 16a and 16b are digitized and angular rate is computed from these signals using a digital signal processor (DSP) 28. In at least one embodiment the conversion of 16a and 16b is accomplished by phase detection circuitry 22a and 22b (as detailed in FIG. 3) to digitize frequency and phase information (FM) and analog-to-digital converters (ADCs) 25a and 25b to digitize amplitude information (AM). A demodulator 27 is shown also coupled between inputs 16a, 16b and DSP 28 and provides the phase reference for demodulating the frequency and amplitude information.

(12) The present disclosure requires that the natural frequencies differ by a small amount, typically in the range of 1 to 100 Hz. In at least one preferred embodiment, the resonator is symmetric with respect to the Coriolis coupled vibration modes. Accordingly, each Coriolis coupled vibration mode is preferably configured with identical sensing and actuating means, and the natural frequencies of the vibration modes are closely matching. The values of the desired natural frequencies and their difference can be set for example by precision manufacturing, or by trimming, or with a phase frequency control circuit which measures the natural frequencies and their differences and makes adjustments to the transducer 12.

(13) FIG. 2A and FIG. 2B illustrate example signals 30, 50. In FIG. 2A are seen displacements x(t) 32 along with the frequency Ω.sub.FM_m modulation 34 and amplitude modulation Ω.sub.AM_m 36 of 32 induced by angular rate signal Φ.sub.z. In FIG. 2B are shown displacements x(t) 52 and y(t) 54, for example as measured from circuit nodes 16a, 16b from FIG. 1 and the phase difference Φ.sub.xy 56, and frequency Ω.sub.FM_m 58 and amplitude Ω.sub.AM_m 60 modulation of 56 and 58.

(14) Displacements x(t) and y(t) can be understood as sinusoidal signals with varying amplitudes and phase ϕ.sub.x and ϕ.sub.y. The phase difference Φ.sub.xy=ϕ.sub.x−ϕ.sub.y, determines the coupling between the modes: for Φ.sub.xy=0 or 180 degrees, respectively, such that rate in one mode modifies the amplitude of oscillation of the other mode. Likewise, for Φ.sub.xy=90 or 270 degrees, wherein rate in one mode modifies the phase of oscillation of the other mode. For illustration, FIG. 2A shows the displacement x(t) 32 for a constant rate as well as the variations of the frequency Ω.sub.FM_m 34 and amplitude Ω.sub.AM_m 36. These are modulated by the continually changing phase difference Φ.sub.xy, as shown in FIG. 2B. For example, for Φ.sub.xy=0 degrees (i.e., the zero crossings of x(t) and y(t) coincide) the change of amplitude is at a maximum while the frequency is unchanged. Conversely, for Φ.sub.xy=90 degrees (i.e., the zero crossing of y(t) coincides with the peak of x(t)) the change of frequency is at a maximum while the amplitude is unchanged. Rate can therefore be determined by measuring the instantaneous frequency and amplitude change of x(t) and y(t) and demodulating these signals with Φ.sub.xy. For the systems where the displacement amplitude is regulated with a feedback loop, the value of applied feedback forces can be used to extract the amplitude modulation Ω.sub.AM_m instead of using the amplitude of the displacement signals. This implementation is also less sensitive to the Q-factors and natural frequencies of the resonators.

(15) It should also be appreciated that the present disclosure is based on utilizing a mechanical resonator (resonator 12 in FIG. 1) with at least two modes of vibration coupled by the Coriolis Effect. The resonator 12 is shown to have two modes of vibration along the x- and y-directions. A mode of vibration is defined as any number of independent mechanical deformations of the resonator, where the total deformation can be described as the sum of the deformations of the individual modes. A deformation should be understood as a continuous vector field which describes the mechanical displacement of each point of the resonator as a function of its position in three-dimensional space. Deformations are, in general, considered to be time varying. A resonator actuated upon by a forced deformation in a particular vibration mode will store a corresponding amount of strain energy. When released, the resonator will attempt to dissipate this strain energy. For a lossless resonator, the result of this action is a sinusoidal vibration of the resonator in that particular vibration mode. The inverse of the time it takes to complete one cycle of vibration is called the natural frequency of the vibration mode. The amplitude, or envelope, of the displacement vibration is defined as half the peak-to-peak displacement at an antinode of the deformation. An antinode is a position of maximum deformation of the particular vibration mode. The amplitude of the velocity of the vibration is equal to the oscillation frequency multiplied by the amplitude of the displacement vibration.

(16) In order for the resonator to be sensitive to rate, the modes should be coupled by the Coriolis Effect. The Coriolis Effect is an effect observed in the rotating frame, whereby energy from one mode of the resonator is coupled into another mode of the resonator, such as from the x-axis to the y-axis, or between orthogonal modes, including from a first mode of vibration to a second mode of vibration. Modes that are coupled by the Coriolis Effect will occasionally be referred to as axes and given names which correspond to their direction of vibration in the rotating frame, such as the x-axis. The rotating frame is a frame of reference perceived by an observer of the resonator which is also rotating at the same rate as the resonator. Common mechanical resonators include pendulums, lumped mass-spring systems, rings, disks, and hemispheres.

(17) Methods of operation will occasionally be described by the trajectory that a proof mass follows when using the particular mode of operation. These descriptions can only be literally interpreted for pendulum or lumped mass-spring systems, but an extension to rings, disks, and hemispheres can be understood in terms of combinations of independent vibrations of the modes. For example, a pendulum swinging at a 45 degree angle corresponds to equal amplitude, in-phase vibrations on both orthogonal modes. A circular pattern corresponds to equal amplitude oscillations of both orthogonal modes with a constant 90 degree phase difference.

(18) The resonator includes means for transducing mechanical deformations into electrical signals (sensing) and for transducing electrical signals to mechanical deformations (actuating) for each mode of vibration. This is typically accomplished electrostatically through addition of electrodes forming either parallel plate or comb-type capacitors to the resonator. Other transducing means include piezoelectric, magnetic, or optical transducers.

(19) The sustaining circuitry is connected to the resonator through sensing and actuating means. The frequency of oscillation observed and sustained by the sustaining circuitry depends on the input rotation rate, natural frequencies, and unwanted parasitic couplings. The dependence on input rate is introduced because the sensing and actuating means lie in the rotating frame. The frequency observed by the sustaining circuitry will be referred to as the oscillation frequency, which is different than the natural frequency.

(20) In at least one embodiment, the sustaining circuitry for each mode comprises variable gain amplifiers (VGAs). Amplitude controllers adjust the gain of the VGAs in order to maintain a mechanical vibration of a particular amplitude. In this case the amplitude modulation Ω.sub.AM_m is equal to the control inputs the control inputs to the VGAs.

(21) FIG. 3 illustrates an example embodiment 70 of a circuit for determining frequency modulation Ω.sub.FM_m of x and y. It should be appreciated that the circuit is shown by way of example and not limitation, as the present disclosure can be configured to operate with various circuitry which is available for determining frequency information from signals. Two of these circuits are used, one each for x(t) and y(t), as shown in FIG. 1 by 22a and 22b. Referring to FIG. 3 x(t), or y(t), is received at continuous-time comparator 74, which converts it to a square wave FM signal 76 for receipt by a delta-sigma frequency-to-digital converter (ΔΣFDC) 80. FM signal 76 is then processed by the exclusive-OR gate 79, which multiplies it by edge selection signal 78, thus resulting in a frequency-doubled signal 75 which is fed to the phase-frequency detector (PFD) 88. By way of example and not limitation, PFD 88 comprises two D-flip flops 85, 87 receiving signals 75, 83 at their clock inputs, and having their D inputs pulled up to V.sub.dd, with Q outputs coupled through AND gate 83 to drive both reset (R) inputs.

(22) Charge pump 92 includes switches 98 and 100 along with current sources 94 and 96, which outputs a pulse of charge proportional to the phase difference between signals 75 and 83. The output of the phase-frequency detector 88 controls switch 98, which enables current source 94. The fine timer 84 outputs a “done” signal 83 received by edge-to-pulse converter 90, which enables “down” switch 100 for a fixed period of time (e.g., 32 ns). The charge pulse is integrated by the active integrator 102 comprising operational transconductance amplifier 104 and feedback capacitor 106. The voltage across capacitor 106 is quantized by analog-to-digital converter 108, producing the digital output signal 112. This signal is fed to arithmetic block 114, which computes the values to be loaded into coarse timer 82 and fine timer 84. Output signal 112 is first processed by loop filter 115 upon which adder 118 adds a center frequency offset 116. The result is then split, with the least significant bits directed to adder 120, which together with register 122 comprises a digital integrator. The output of adder 120 is fed to fine timer 84 and to register 122. Adder 124 sums the most significant bits from adder 118 with the carry bit from adder 120, with the result coupled to coarse timer 82. Control logic block 86 generates edge selection signal 78, which is toggled after every phase comparison. This signal is used by XOR gate 79 to convert falling edges to rising edges, and is also used by multiplexor (mux) 121 to select one of two values for center frequency 116, which enables the loop to handle signals with an asymmetric duty cycle without requiring a large input frequency range. Output 112 of the ΔΣFDC 80 is a 2nd order sigma-delta modulated version of the input oscillation frequency. It is effectively a digital signal representative of the oscillation frequency which is preferably further processed by a DSP, ASIC, or FPGA to extract input angular rate, based on calculations.

(23) To reduce power, the coarse timer (e.g., an 11-bit counter) 82 runs at a reduced clock rate, with correspondingly low time resolution. The coarse timer is used to enable a fine timer (e.g., a 3-bit counter) 84 near the time of the input edge. As an example, the coarse timer receives a clock with a period of 32 ns (31.25 MHz), while the fine timer receives a clock with a period of 4 ns (250 MHz). Because fine timer 84 effectively only delays the “done” signal of the coarse timer 82, it is necessary to integrate the timer input so that all future edges are delayed by the same amount. This is performed by a modulo integrator comprised of adder 120 and register 122. The carry bit is fed to adder 124, which adds one extra cycle of delay to the coarse timer when the modulo integrator wraps around.

(24) In order to avoid nonlinearity due to mismatch between the up current 94 and down current 96, the down current source 96 is switched on using switch 100 for a constant duration (e.g., 32 ns) during every phase comparison cycle. The on-duration is set by pulse generator 90. The top current source 94 is switched on and off through switch 98. Current source 94 is switched on in response to output 83 from fine counter 84, which causes the output of D flip-flop 85 to go high. The rising edge of signal 75 produced by XOR gate 79 causes the output of D flip-flop 87 to go high; this causes the output of AND gate 83 to go high, resetting flip-flops 87 and 85 and turning off switch 98. Any mismatch between the two current sources 94 and 96 manifests itself as a static phase offset between signals 75 and 83; since the frequency demodulator measures the time derivative of phase, this offset is effectively eliminated and does not appear at the output 112.

(25) While the delta-sigma frequency-to-digital converter (ΔΣFDC) is a preferred embodiment for frequency measurement due to its excellent linearity and scale factor accuracy, numerous other circuit architectures can be used in the present disclosure to extract rate information from the frequency-modulated signal. A conventional analog PLL can be used to demodulate FM, though scale factor stability and linearity is limited by the analog accuracy of the VCO. Another architecture which can be utilized comprises a slope discriminator followed by an envelope detector; this architecture is also subject to analog imperfections. Another possible implementation comprises a high-resolution analog-digital converter to directly digitize the sinusoidal signal from the oscillators, with the FM discrimination and further signal processing performed by a digital signal processor (e.g., comprising a programmable IC and/or a fixed-function application-specific integrated circuit (ASIC)). Another option is to measure the period of the FM signal directly, using either a high-resolution time-to-digital converter or a frequency counter. One of ordinary skill in the art will appreciate that various other methods of FM demodulation, or frequency measurement, may be alternatively utilized or adapted to work with the teachings of the present invention.

(26) FIG. 4 illustrates an example embodiment 170 of a circuit for extracting signals cos Φ.sub.xy and sin Φ.sub.xy needed for demodulating Ω.sub.AM_m and Ω.sub.FM_m. The axis oscillation signals X=x(t) 172 and Y=y(t) 174 are fed into respective band-pass limiters 176, 178 to remove amplitude information. The resulting axis oscillation signals, without amplitude information (cos Φ.sub.xy and sin Φ.sub.xy) are then multiplied at multiplier 180, then multiplied 182 by a constant, seen here as multiplication by the value 2, then passed through a low-pass filter 184, resulting in a signal cos Φ.sub.xy, 185a equal to the cosine of the phase difference between the inputs. Low-pass filter 184 removes the image at the sum of the input frequencies created by the multiplication. The low pass filtered signal is then phase shifted by 90 degrees at a phase shifter 186, to construct the rate reference sin Φ.sub.xy 185b for FM readout. With references cos Φ.sub.xy and sin Φ.sub.xy being known, angular rate can then be determined from Ω.sub.FM_m and Ω.sub.AM_m.

(27) The rate reference 185b is multiplied at multiplier 188 by Ω.sub.FM_m (190 in FIGS. 4 and 112 in FIG. 3) and the result is filtered at low-pass filter 192 to produce a final output 194. The low-pass filter removes wide-band noise and errors. The system used for producing the angular rate from the oscillation frequency and rate reference is referred to as the rate demodulator.

(28) FIG. 5 illustrates an example embodiment 210 showing amplitudes of assumed two-sided spectra of the rotation signal Ω.sub.z (angular rate) 212 and the frequency Ω.sub.FM_m and amplitude Ω.sub.AM_m 214 of the envelope of the modes. The latter are shifted by the relative phase Φ.sub.xy of x and y:
Ω.sub.FM_m=Ω.sub.z sin Φ.sub.xy
Ω.sub.AM_m=Ω.sub.z cos Φ.sub.xy
where Φ.sub.xy≅Δω.Math.t and Δω=2πΔf is the frequency difference between the resonant frequencies of the x and y modes. The signals Ω.sub.FM_m and Ω.sub.AM_m are in quadrature, differing only by their phase, and are therefore not shown individually. The demodulated signals
Ω.sub.FM_d=Ω.sub.FM_m sin Φ.sub.xy
Ω.sub.AM_d=Ω.sub.AM_m cos Φ.sub.xy
comprise the rotation signal Ω.sub.z plus spectra at ±2Δf known in the art as image spectra or spectral images and a well known artifact of AM modulation by those skilled in the art. The unknown rate {circumflex over (Ω)} 220 can be obtained by low-pass filtering either Ω.sub.FM_d 216 or Ω.sub.AM_d 218, or more preferably by summation, thus exploiting the opposite polarity of the image frequencies in Ω.sub.FM_d or Ω.sub.AM_d:
{circumflex over (Ω)}=∩.sub.FM_d+Ω.sub.AM_d=Ω.sub.z

(29) FIG. 6 illustrates an example embodiment 230 which shows the same spectra as FIG. 5, but for the case where the bandwidth BW of Ω.sub.z 232 is greater than Δf. The overlap of the desired signal spectra of Ω.sub.FM_m and Ω.sub.AM_m 234 and their spectral images amounts to signal corruption, also known in the art as aliasing. The demodulated signals Ω.sub.FM_d 236 and Ω.sub.AM_m 238 suffer from the same problem and are therefore not a correct measure of the unknown rate Ω.sub.z.

(30) One solution for addressing this problem is to increase the frequency split between modes x and y. This, however, has adverse effects on gyroscope performance, for example resulting in an increase of the short term noise, also called angular random walk (ARW). Alternatively, if rate is extracted from the sum of Ω.sub.FM_d and Ω.sub.AM_m, as suggested above, the images from Ω.sub.FM_d and Ω.sub.AM_d cancel due to the quadrature nature of Ω.sub.FM_m and Ω.sub.AM_m resulting in a rate estimate {circumflex over (Ω)} 240, where {circumflex over (Ω)}=Ω.sub.z, that is free from aliasing errors, thus removing the bandwidth restriction to rates. Hence rates with frequencies greater than Δf are correctly measured.

(31) Practical implementations suffer from errors that should be eliminated to maximize gyroscope performance. For example, the amplitudes of Ω.sub.FM_d and Ω.sub.AM_d should be matched precisely for the images to be cancelled accurately. Another error results from slowly varying noise components Ω.sub.ε.sub.FM and Ω.sub.ε.sub.AM being added by the sustaining circuits and amplitude and frequency measurement circuitry as well as imperfect extraction of the relative phase Φ.sub.xy used for demodulation. Yet another error results from undesired quadrature motion Ω.sub.q_FM and Ω.sub.d_AM in the x and y modes that are modulated to the same frequency as the signal Ω.sub.z but 90 degrees phase shifted. With these errors added, the equations for the modulated FM and AM rate become:

(32) ${\Omega }_{\mathrm{FM}\_m}=\underset{\mathrm{signal}}{\underbrace{{\Omega }_z\sin {\Phi }_{\mathrm{xy}}}}+\underset{\mathrm{error}}{\underbrace{{\Omega }_{{\epsilon }_{\mathrm{FM}}}+{\Omega }_{q\_\mathrm{FM}}\cos {\Phi }_{\mathrm{xy}}}}{\Omega }_{\mathrm{AM}\_m}=\underset{\mathrm{signal}}{\underbrace{{\Omega }_z\cos {\Phi }_{\mathrm{xy}}}}+\underset{\mathrm{error}}{\underbrace{{\Omega }_{{\epsilon }_{\mathrm{AM}}}+{\Omega }_{q\_\mathrm{AM}}\sin {\Phi }_{\mathrm{xy}}}}$

(33) FIG. 7 illustrates an example embodiment 250 which shows the same spectra as FIG. 6, with bandwidth of Ω.sub.z 252 greater than Δf, an overlap spectra of Ω.sub.FM_m and Ω.sub.AM_m 254, demodulated signals Ω.sub.FM_d 256 and Ω.sub.AM_m 258, and unknown rate {circumflex over (Ω)} 260. However, this figure assumes added noise spectra Ω.sub.ε.sub.FM and Ω.sub.ε.sub.AM. Demodulation shifts these errors to Δf, which for BW>Δf is inside the useful signal band and therefore results in corruption:

(34) $\begin{array}{c}{\Omega }_{\mathrm{FM}\_d}={\Omega }_{\mathrm{FM}\_m}\sin {\Phi }_{\mathrm{xy}}\\ =\frac{{\Omega }_z}{2}-\frac{{\Omega }_z}{2}\cos 2{\Phi }_{\mathrm{xy}}+{\Omega }_{{\epsilon }_{\mathrm{FM}}}\sin {\Phi }_{\mathrm{xy}}+\frac{{\Omega }_{q\_\mathrm{FM}}}{2}\sin 2{\Phi }_{\mathrm{xy}}\end{array}\begin{array}{c}{\Omega }_{\mathrm{AM}\_d}={\Omega }_{\mathrm{AM}\_m}\cos {\Phi }_{\mathrm{xy}}\\ =\frac{{\Omega }_z}{2}+\frac{{\Omega }_z}{2}\cos 2{\Phi }_{\mathrm{xy}}+{\Omega }_{{\epsilon }_{\mathrm{AM}}}\cos {\Phi }_{\mathrm{xy}}+\frac{{\Omega }_{q\_\mathrm{AM}}}{2}\sin 2{\Phi }_{\mathrm{xy}}\end{array}$

(35) Adding these demodulated signals gives an estimate {circumflex over (Ω)}.sub.z of the rate signal Ω.sub.z:

(36) $\begin{array}{cc}\begin{array}{c}{\hat{\Omega }}_z={\Omega }_{\mathrm{FM}\_d}+{\Omega }_{\mathrm{AM}\_d}\\ ={\Omega }_z+\underset{\underset{\mathrm{error}}{︸}}{{\Omega }_{{\epsilon }_{\mathrm{FM}}}\sin {\Phi }_{\mathrm{xy}}+{\Omega }_{{\epsilon }_{\mathrm{AM}}}\cos {\Phi }_{\mathrm{xy}}+\frac{{\Omega }_{q\_\mathrm{FM}}+{\Omega }_{q\_\mathrm{AM}}}{2}\sin 2{\Phi }_{\mathrm{xy}}}\end{array}& \left(1\right)\end{array}$

(37) The error can be calculated as follows:

(38) ${\Omega }_{\mathrm{FM}\_m}\cos {\Phi }_{\mathrm{xy}}=\frac{{\Omega }_z}{2}\sin 2{\Phi }_{\mathrm{xy}}+{\Omega }_{{\epsilon }_{\mathrm{FM}}}\cos {\Phi }_{\mathrm{xy}}+\frac{{\Omega }_{q\_\mathrm{FM}}}{2}+\frac{{\Omega }_{q\_\mathrm{FM}}}{2}\cos 2{\Phi }_{\mathrm{xy}}{\Omega }_{\mathrm{AM}\_m}\sin {\Phi }_{\mathrm{xy}}=\frac{{\Omega }_z}{2}\sin 2{\Phi }_{\mathrm{xy}}+{\Omega }_{{\epsilon }_{\mathrm{AM}}}\sin {\Phi }_{\mathrm{xy}}+\frac{{\Omega }_{q\_\mathrm{AM}}}{2}-\frac{{\Omega }_{q\_\mathrm{AM}}}{2}\cos 2{\Phi }_{\mathrm{xy}}$

(39) Subtracting the two equation yields:

(40) $\begin{array}{cc}{\Omega }_{\mathrm{FM}\_m}\cos {\Phi }_{\mathrm{xy}}-{\Omega }_{\mathrm{AM}\_m}\sin {\Phi }_{\mathrm{xy}}=\underset{\underset{\mathrm{quadrature}\mathrm{leakage}}{︸}}{\frac{{\Omega }_{q\_\mathrm{FM}}-{\Omega }_{q\_\mathrm{AM}}}{2}}+\underset{\mathrm{error}}{\underbrace{{\Omega }_{{\epsilon }_{\mathrm{FM}}}\cos {\Phi }_{\mathrm{xy}}-{\Omega }_{{\epsilon }_{\mathrm{AM}}}\sin {\Phi }_{\mathrm{xy}}+\frac{{\Omega }_{q\_\mathrm{FM}}+{\Omega }_{q\_\mathrm{AM}}}{2}\cos 2{\Phi }_{\mathrm{xy}}}}& \left(2\right)\end{array}$

(41) The error term in Eq. (2) equals the error in Eq. (1) shifted 90 degrees. An error-free estimate {circumflex over (Ω)}.sub.z of the rate signal Ω.sub.z is obtained by high-pass HPF[•] filtering after performing Eq. (2) to remove the quadrature leakage term, applying the Hilbert transform H{ } to shift the remainder by 90 degrees phase, and subtracting the result from Eq. (1):
{circumflex over (Ω)}.sub.z=Ω.sub.FM_d+Ω.sub.AM_d−H{HPF[Ω.sub.FM_m cos Φ.sub.xy−Ω.sub.AM_m sin Φ.sub.xy]}=Ω.sub.z

(42) Although this, in principle, results in an error-free output, practical implementations suffer from residual errors caused by temperature variations. It would therefore be useful to know the precise temperature of the transducer proof-mass to cancel these errors. The relationship between temperature and errors can be established for example based on precise experimental characterization of the device in the factory.

(43) Transducers fabricated from materials such as Silicon or Quartz exhibit well known variations of resonant frequency as a function of device temperature. The error Ω.sub.ε.sub.FM in Eq. (2) above includes the resonant frequency variation and is therefore a precise proxy for device temperature. It can be obtained from Eq. (2) from synchronous demodulation with cos Φ.sub.xy and used in a temperature cancellation circuit that reduces residual gyroscope errors due to temperature variations.

(44) The enhancements described in the presented technology can be readily implemented within various gyroscope hardware configurations, a number of which are implemented to include one or more computer processor devices (e.g., digital signal processor (DSP), central processing unit (CPU), microprocessor, microcontroller, computer enabled application specific integrated circuit (ASIC), etc.) and associated memory storing instructions (e.g., RAM, DRAM, FRAM, NVRAM, FLASH, computer readable media, etc.) whereby programming (instructions) stored in the memory are executed on the processor to perform the steps of the various process methods described herein. The presented technology is non-limiting with regard to the use of processors, memory and computer-readable media, insofar as these are non-transitory, and thus not constituting a transitory electronic signal.

(45) Embodiments of the present technology may be described herein with reference to flowchart illustrations of methods and systems according to embodiments of the technology, and/or procedures, algorithms, steps, operations, formulae, or other computational depictions, which may also be implemented as computer program products. In this regard, each block or step of a flowchart, and combinations of blocks (and/or steps) in a flowchart, as well as any procedure, algorithm, step, operation, formula, or computational depiction can be implemented by various means, such as hardware, firmware, and/or software including one or more computer program instructions embodied in computer-readable program code. As will be appreciated, any such computer program instructions may be executed by one or more computer processors, including without limitation a general purpose computer or special purpose computer, or other programmable processing apparatus to produce a machine, such that the computer program instructions which execute on the computer processor(s) or other programmable processing apparatus create means for implementing the function(s) specified.

(46) Accordingly, blocks of the flowcharts, and procedures, algorithms, steps, operations, formulae, or computational depictions described herein support combinations of means for performing the specified function(s), combinations of steps for performing the specified function(s), and computer program instructions, such as embodied in computer-readable program code logic means, for performing the specified function(s). It will also be understood that each block of the flowchart illustrations, as well as any procedures, algorithms, steps, operations, formulae, or computational depictions and combinations thereof described herein, can be implemented by special purpose hardware-based computer systems which perform the specified function(s) or step(s), or combinations of special purpose hardware and computer-readable program code.

(47) Furthermore, these computer program instructions, such as embodied in computer-readable program code, may also be stored in one or more computer-readable memory or memory devices that can direct a computer processor or other programmable processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory or memory devices produce an article of manufacture including instruction means which implement the function specified in the block(s) of the flowchart(s). The computer program instructions may also be executed by a computer processor or other programmable processing apparatus to cause a series of operational steps to be performed on the computer processor or other programmable processing apparatus to produce a computer-implemented process such that the instructions which execute on the computer processor or other programmable processing apparatus provide steps for implementing the functions specified in the block(s) of the flowchart(s), procedure (s) algorithm(s), step(s), operation(s), formula(e), or computational depiction(s).

(48) It will further be appreciated that the terms “programming” or “program executable” as used herein refer to one or more instructions that can be executed by one or more computer processors to perform one or more functions as described herein. The instructions can be embodied in software, in firmware, or in a combination of software and firmware. The instructions can be stored local to the device in non-transitory media, or can be stored remotely such as on a server or all or a portion of the instructions can be stored locally and remotely. Instructions stored remotely can be downloaded (pushed) to the device by user initiation, or automatically based on one or more factors.

(49) It will further be appreciated that as used herein, that the terms processor, hardware processor, computer processor, central processing unit (CPU), and computer are used synonymously to denote a device capable of executing the instructions and communicating with input/output interfaces and/or peripheral devices, and that the terms processor, hardware processor, computer processor, CPU, and computer are intended to encompass single or multiple devices, single core and multicore devices, and variations thereof.

(50) From the description herein, it will be appreciated that the present disclosure encompasses multiple embodiments which include, but are not limited to, the following:

(51) 1. An electromechanical system, comprising: (a) a mechanical resonator having a first mode of vibration and an associated first natural frequency, and a second mode of vibration having an associated second natural frequency, wherein angular rate of motion input couples energy between said first mode of vibration and said second mode of vibration; (b) sensors and actuators for each of said first mode and said second mode for transducing an electrical signal into a mechanical vibration and transducing a mechanical vibration into an electrical signal; (c) sustaining circuitry connected to the actuators to maintain vibrations in the first mode at a first frequency and the second mode at the second frequency; (d) an amplitude and frequency measurement circuit configured for determining the instantaneous amplitude and frequency variation and outputting modulated versions of mode vibration amplitude and frequency as a representation of vibration amplitude and frequency; (e) a phase determination circuit configured for receiving mode vibration signals and generating rate reference signals, and outputting demodulated frequency and amplitude rate signals; (f) a processing element configured for combining demodulated frequency and amplitude rate signals to generate an output rate in which the image spectra of the demodulation are eliminated from the angular rate determination.

(52) 2. The apparatus, system or method of any preceding or following embodiment, wherein said electromechanical system comprises a continuous mode reversal gyroscope.

(53) 3. The apparatus, system or method of any preceding or following embodiment, wherein said gyroscope is configured with a proof mass configured to move along at least two orthogonal axes.

(54) 4. The apparatus, system or method of any preceding or following embodiment, wherein said proof mass is suspended by springs.

(55) 5. The apparatus, system or method of any preceding or following embodiment, wherein said mechanical resonator comprises a gyroscope transducer.

(56) 6. The apparatus, system or method of any preceding or following embodiment, wherein said mechanical resonator is symmetric with respect to Coriolis coupled vibration modes.

(57) 7. The apparatus, system or method of any preceding or following embodiment, wherein said first natural frequency and said second natural frequency of said mechanical resonator differs by a range of approximately 1 to 100 Hz.

(58) 8. The apparatus, system or method of any preceding or following embodiment, wherein said mechanical resonator is manufactured or trimmed to achieve said difference between said first natural frequency and said second natural frequency, or a phase frequency control circuit is utilized which is configured to determined these natural frequencies and make adjustments to operation the mechanical resonator.

(59) 9. The apparatus, system or method of any preceding or following embodiment, wherein said system is configured to infer an angular rate of motion from the gains of the sustaining circuitry, the frequency of vibration of the first mode, and the frequency of vibration of the second mode.

(60) 10. The apparatus, system or method of any preceding or following embodiment, wherein said system is configured to infer an angular rate of motion from the sum of the gains of the sustaining circuitry and the frequencies of vibration of the first mode and the frequency of vibration of the second mode.

(61) 11. The apparatus, system or method of any preceding or following embodiment, wherein said system is configured to infer an estimate of the errors of the inferred angular rate from gains of the sustaining circuitry and the frequencies of vibration of the first mode and the second mode, and subtracting this estimate from the inferred output rate.

(62) 12. The apparatus, system or method of any preceding or following embodiment, further comprising a circuit configured for estimating temperature of the mechanical resonator from an estimate of the errors of the inferred angular rate.

(63) 13. The apparatus, system or method of any preceding or following embodiment, further comprising a circuit configured to compensate for errors due to variations in the temperature of the mechanical resonator using the temperature estimate.

(64) 14. An electromechanical system, comprising: (a) a mechanical resonator having a first mode of vibration and an associated first natural frequency, and a second mode of vibration having an associated second natural frequency, wherein angular rate of motion input couples energy between said first mode of vibration and said second mode of vibration; (b) sensors and actuators for each of the first mode and second mode of vibration for transducing an electrical signal into a mechanical vibration and transducing a mechanical vibration into an electrical signal; (c) sustaining circuitry with variable gains connected to the actuators to maintain vibrations in the first mode and the second mode with gains adjusted to maintain substantially constant, non-zero velocity amplitude vibrations in the first mode at the first frequency of vibration and the second mode at the second frequency of vibration; (d) an amplitude and frequency measurement circuit configured for determining the instantaneous amplitude and frequency variation and outputting modulated versions of mode vibration amplitude and frequency as a representation of vibration amplitude and frequency; (e) a phase determination circuit configured for receiving mode vibration signals and generating rate reference signals, and outputting demodulated frequency and amplitude rate signals; (f) output circuitry to infer an angular rate of motion in response to combining demodulated frequency and amplitude rate signals to generate an output rate in which the image spectra of the demodulation are eliminated from the angular rate determination.

(65) 15. The apparatus, system or method of any preceding or following embodiment, wherein said electromechanical system comprises a continuous mode reversal gyroscope.

(66) 16. The apparatus, system or method of any preceding or following embodiment, wherein said gyroscope is configured with a proof mass configured to move along at least two orthogonal axes.

(67) 17. The apparatus, system or method of any preceding or following embodiment, wherein said proof mass is suspended by springs.

(68) 18. The apparatus, system or method of any preceding or following embodiment, wherein said mechanical resonator comprises a gyroscope transducer.

(69) 19. The apparatus, system or method of any preceding or following embodiment, wherein said mechanical resonator is symmetric with respect to Coriolis coupled vibration modes.

(70) 20. The apparatus, system or method of any preceding or following embodiment, wherein said first natural frequency and said second natural frequency differ by a range of approximately 1 to 100 Hz.

(71) 21. The apparatus, system or method of any preceding or following embodiment, wherein said mechanical resonator is manufactured or trimmed to achieve said difference between said first natural frequency and said second natural frequency, or a phase frequency control circuit is utilized which is configured to determined these natural frequencies and make adjustments to operation the mechanical resonator.

(72) 22. The apparatus, system or method of any preceding or following embodiment, wherein said output circuitry is configured to infer an angular rate of motion from the gains of the sustaining circuitry, and the frequency of vibration of the first mode and the second mode.

(73) 23. The apparatus, system or method of any preceding or following embodiment, wherein said output circuitry is configured to infer an angular rate of motion from a sum of the gains of the sustaining circuitry and the frequencies of vibration of the first mode and second mode.

(74) 24. The apparatus, system or method of any preceding or following embodiment, wherein said output circuitry is configured to infer an estimate of the errors of the inferred angular rate from the gains of the sustaining circuitry and the frequencies of vibration of the first and second mode and subtracting this estimate from the inferred angular rate.

(75) 25. The apparatus, system or method of any preceding or following embodiment, further comprising a circuit configured for estimating the temperature of the electromechanical system from the estimate of the errors of inferred angular rate.

(76) 26. The apparatus, system or method of any preceding or following embodiment, further comprising a circuit configured for compensating for errors due to variations of the temperature of the electromechanical system using the temperature estimate.

(77) 27. A method of determining angular rate in a gyroscope, comprising: (a) transducing an electrical signal into a mechanical vibration of a mechanical resonator having a first mode of vibration and an associated first natural frequency, and a second mode of vibration having an associated second natural frequency, in which angular rate of motion input couples energy between said first mode of vibration and said second mode of vibration; (b) transducing the mechanical vibration of said mechanical resonator into an electrical signal; (c) maintaining vibrations of said mechanical resonator in the first mode at a first frequency and the second mode at the second frequency; (d) determining instantaneous amplitude and frequency variation and outputting modulated versions of mode vibration amplitude and frequency as a representation of vibration amplitude and frequency; (e) receiving mode vibration signals and generating rate reference signals, and outputting demodulated frequency and amplitude rate signals; and (f) combining demodulated frequency and amplitude rate signals to generate an output rate in which the image spectra of the demodulation are eliminated from the angular rate determination.

(78) 28. An electromechanical system, comprising: (a) a mechanical resonator having a first mode of vibration and an associated first natural frequency, and a second mode of vibration having an associated second natural frequency, wherein angular rate of motion input couples energy between said first mode of vibration and second mode of vibration; (b) sensors and actuators for each of the first mode and second mode for transducing an electrical signal into a mechanical vibration and transducing a mechanical vibration into an electrical signal; (c) sustaining circuitry connected to the actuators to maintain vibrations in the first mode at a first frequency and the second mode at the second frequency; and (d) output circuitry to infer an angular rate of motion from the amplitude of vibration of first mode of vibration and second mode of vibration.

(79) 29. An electromechanical system, comprising: (a) a mechanical resonator having a first mode of vibration and an associated first natural frequency, and a second mode of vibration having an associated second natural frequency, wherein angular rate of motion input couples energy between said first mode of vibration and second mode of vibration; (b) sensors and actuators for each of the first mode and second mode for transducing an electrical signal into a mechanical vibration and transducing a mechanical vibration into an electrical signal; (c) sustaining circuitry with variable gains connected to the actuators to maintain vibrations in the first mode at a first frequency and the second mode at the second frequency with the gains adjusted to maintain substantially constant, non-zero velocity amplitude vibrations in the first mode at a first frequency and the second mode at a second frequency; and (d) output circuitry to infer an angular rate of motion from the gains of the sustaining circuitry.

(80) 30. The apparatus, system or method of any preceding or following embodiment, with output circuitry to infer an angular rate of motion from the gains of the sustaining circuitry, the frequency of vibration of the first mode, and the frequency of vibration of the second mode.

(81) 31. The apparatus, system or method of any preceding or following embodiment, with output circuitry to infer an angular rate of motion from the sum of the gains of the sustaining circuitry and the frequencies of vibration of the first and second mode.

(82) 32. The apparatus, system or method of any preceding or following embodiment, with output circuitry to infer an estimate of the errors of the inferred angular rate from the gains of the sustaining circuitry and the frequencies of vibration of the first and second mode and subtracting this estimate from the inferred angular rate.

(83) 33. The apparatus, system or method of any preceding or following embodiment, with circuitry to estimate the temperature of the mechanical system from the estimate of the errors.

(84) 34. The apparatus, system or method of any preceding or following embodiment, with circuitry to compensate for errors due to variations of the temperature of the mechanical system using the temperature estimate.

(85) 35. In a FM gyroscope of the type having a proof mass that is suspended by springs and is free to move along two orthogonal axes or a ring, or similar structure, wherein vibration of the mass can on occur on both axes simultaneously, and wherein a bias error shifting method is used to modulate bias error to a higher frequency where it can be removed by a filter, the improvement comprising using both AM and FM information from the gyroscope transducer to overcome limitation of useable signal bandwidth to less than the modulation frequency.

(86) 36. The apparatus, system or method of any preceding or following embodiment, wherein modulation with the cos and sin terms, respectively, limit the usable bandwidth of either signal; and wherein combining the two outputs (cos and sin terms) eliminates this bandwidth restriction.

(87) As used herein, the singular terms “a,” “an,” and “the” may include plural referents unless the context clearly dictates otherwise. Reference to an object in the singular is not intended to mean “one and only one” unless explicitly so stated, but rather “one or more.”

(88) As used herein, the term “set” refers to a collection of one or more objects. Thus, for example, a set of objects can include a single object or multiple objects.

(89) As used herein, the terms “substantially” and “about” are used to describe and account for small variations. When used in conjunction with an event or circumstance, the terms can refer to instances in which the event or circumstance occurs precisely as well as instances in which the event or circumstance occurs to a close approximation. When used in conjunction with a numerical value, the terms can refer to a range of variation of less than or equal to ±10% of that numerical value, such as less than or equal to ±5%, less than or equal to ±4%, less than or equal to ±3%, less than or equal to ±2%, less than or equal to ±1%, less than or equal to ±0.5%, less than or equal to ±0.1%, or less than or equal to ±0.05%. For example, “substantially” aligned can refer to a range of angular variation of less than or equal to ±10°, such as less than or equal to ±5°, less than or equal to ±4°, less than or equal to ±3°, less than or equal to ±2°, less than or equal to ±1°, less than or equal to ±0.5°, less than or equal to ±0.1°, or less than or equal to ±0.05°.

(90) Additionally, amounts, ratios, and other numerical values may sometimes be presented herein in a range format. It is to be understood that such range format is used for convenience and brevity and should be understood flexibly to include numerical values explicitly specified as limits of a range, but also to include all individual numerical values or sub-ranges encompassed within that range as if each numerical value and sub-range is explicitly specified. For example, a ratio in the range of about 1 to about 200 should be understood to include the explicitly recited limits of about 1 and about 200, but also to include individual ratios such as about 2, about 3, and about 4, and sub-ranges such as about 10 to about 50, about 20 to about 100, and so forth.

(91) Although the description herein contains many details, these should not be construed as limiting the scope of the disclosure but as merely providing illustrations of some of the presently preferred embodiments. Therefore, it will be appreciated that the scope of the disclosure fully encompasses other embodiments which may become obvious to those skilled in the art.

(92) All structural and functional equivalents to the elements of the disclosed embodiments that are known to those of ordinary skill in the art are expressly incorporated herein by reference and are intended to be encompassed by the present claims. Furthermore, no element, component, or method step in the present disclosure is intended to be dedicated to the public regardless of whether the element, component, or method step is explicitly recited in the claims. No claim element herein is to be construed as a “means plus function” element unless the element is expressly recited using the phrase “means for”. No claim element herein is to be construed as a “step plus function” element unless the element is expressly recited using the phrase “step for”.