Method for temperature-compensated accelerometer measurement, with at least a device comprising two vibrating resonators

[0001] The present invention relates to the field of accelerometric measurements, and more specifically it relates to improvements made in the field of accelerometric measurements performed with devices having two vibrating resonators working respectively in stress states of opposite signs (in general in tension and compression).

[0002] Thus, the use is already known of devices having two vibrating resonators comprising two respective vibrating beams, parallel and identical, one working in tension and the other in compression, to detect a component of acceleration (see for example documents FR 00 10675 and FR-A-2 784 752 in the name of the Applicant; FR-A-2 685 964, FR-A-2 784 752). As illustrated by way of example in FIG. 1 of the appended drawings, such an accelerometer device 1 comprises, in a very schematic manner, the two miniature accelerometers Acc₁, Acc₂ with nonslaved pendular masses. Each of them comprises a pendulum mass (or seismic mass) M urged toward a rest position by a link with a stand S which consists of a vibrating beam or blade P (the means of vibrational excitation and of detection of the vibration not being illustrated) and of a hinge-forming member C. The two beams P₁, P₂ are identical and parallel and are disposed in a reverse position with respect to each other, so that, under the action of the two respective identical seismic masses M₁, M₂, one P₁ of the beams works in compression and the other P₂ of the beams works in tension. The two individual accelerometers Acc₁, Acc₂ are together enclosed in a leaktight box B in which an intense vacuum is generally established.

[0003] When the accelerometer device 1 is subjected to an acceleration γ directed along the sensitive axis, that is to say for example perpendicular to the span of the vibrating beams P₁, P₂, it results in a variation of the resonant frequencies of the beams related to the acceleration applied and to the stress (compression or tension) undergone by the beams:
f₁=f₀₁+k₁₁γ−k₁₂γ²+k₁₃γ³
f₂=f₀₂−k₂₁γ−k₂₂γ²+k₂₃γ³
where f₀₁ and f₀₂ are the resonant frequencies of the two respective beams P₁ and P₂ in the absence of acceleration, and the coefficients k are constants.

[0004] This type of accelerometer device has several benefits, in particular by reason of its low bulk (flat structure whose thickness may be less than 1 mm for in-plane dimensions of the order of a centimeter) and of its simple and economical fabrication by methods of chemical etching on a crystal wafer, generally piezoelectric, but possibly silicon, with possibility of simultaneous collective fabrication of several accelerometers on one and the same wafer of material.

[0005] However, in the type of accelerometer device considered having double vibrating resonator, the resonant frequencies of the vibrating beams P₁, P₂ depend not only on the acceleration γ to which the device is subjected, but also on the temperature T of the beams. One is therefore led to supplement the above equations with a polynomial in T: $f_1=f_{01}+\prod _1^1\left({\gamma }^n\right)+\prod _1^2\left(T^n\right)$ $f_2=f_{02}+\prod _2^1\left({\gamma }^n\right)+\prod _2^2\left(T^n\right)$

[0006] A simple calculation of the value of the acceleration γ as a function of the measured frequency f₁ and/or f₂ requires knowledge of the temperature of the beams P₁, P₂. However, it is not possible, in practice, to dispose an appropriate sensor, in the box B, in immediate proximity to the vibrating beams, nor to dispose a thermal sensor in actual contact with the vibrating beams P₁, P₂ (which would constitute the ideal solution for obtaining an exact value of the temperature) by reason of the perturbation that the presence of this sensor would cause on the vibration of the beams. Finally, recourse to a method of thermoelectric measurement would be liable to perturb the temperature of the beam itself.

[0007] An object of the invention is therefore to propose a method of calculating the value of the acceleration, temperature-compensated, on the basis of the measurements of the frequencies f₁, f₂ of the two vibrating beams P₁, P₂, without the need for recourse to any thermal detection means disposed in contact with the vibrating beams.

[0008] For these purposes, the invention proposes a method for accelerometer measurement by determination of the value, compensated as a function of temperature, of an acceleration on the basis of measurements of frequencies made on one accelerometer device having two vibrating resonators of similar geometries working respectively in stress states of opposite signs, which process is characterized in that according to the invention it comprises the steps which follow:

• • [0009] the device being subjected to a component of acceleration, the frequencies f₁ and f₂ of the respective vibrations of the two vibrating resonators are determined, given that the respective values of these two frequencies are influenced by the acceleration and by the temperature;
  • [0010] the sum f₁+f₂ of the two measured frequencies is determined, given that this sum is influenced by the temperature;
  • [0011] on the basis of a pre-established correlation of the variation of the quantity f₁+f₂ as a function of temperature, the possible temperature values T_i corresponding to the value of the abovedetermined sum f₁+f₂ are plotted;
  • [0012] the value (G_mes) of a magnitude representative of the temperature in proximity to the device having two vibrating resonators is measured;
  • [0013] with the aid of a previously established transfer function H, an estimated value G_estimated=H(G_mes) of abovesaid magnitude is determined on the basis of said measured value G_mes, this estimated magnitude being representative of the estimated temperature of the vibrating resonators;
  • [0014] on the basis of the estimated value G_estimated of said magnitude and of the possible temperature values T_i, the probability P(T_i) that the value of the temperature of the resonators is T_i, is determined, for each value T_i;
  • [0015] a weighted value T of the temperature of the vibrating resonators is calculated T=ΣP(T_i)×T_i;
  • [0016] finally on the basis of one and/or the measured frequencies f₁ and f₂ and of the weighted value T of the temperature, the temperature-compensated value of the acceleration undergone by the accelerometer device is determined.

[0017] By virtue of the provisions of the invention, it is possible to take account, in the calculation of the value of the acceleration γ, of a temperature value T, or of a magnitude representative of this temperature (for example a magnitude varying linearly as a function of temperature) in the vicinity of the resonators—and hence intrinsically inexact since it is not the value of the temperature of the resonators—and on the basis of a choice, involving a probability calculation, between several possible values of the temperature.

[0018] In practice, the two vibrating resonators cannot be strictly identical from the geometrical point of view: it is then desirable to take account of this geometrical difference by correcting the sum f₁+f₂ with an error term estimated on the basis of the acceleration calculated during the previous calculation cycle, if it exists.

[0019] Again in practice, the implementation of the method of the invention with a polynomial π(Tⁿ) of the high order n turns out to be complex, by reason of the high number n of mathematical solutions, that is to say of the number n of possible values of temperature T_i corresponding to a given value of the frequency f₁+f₂. This also results in a longer calculation of the probabilities P(T_i) and of the value T=ΣP(T_i)×T_i. A value T of sufficient accuracy is obtained, in a concrete manner, with a polynomial of order 2: π(T²) or of order 3: π(T³).

[0020] The mathematical relation between the frequency, the acceleration and the temperature may be written
f=f₀+π(γ³)+a(T−T_o)²
in the case of a polynomial in T of order 2 and
f=f₀+π(γ³)+a(T−T_o)²+b(T−T_o)³
in the case of a polynomial in T of order 3.

[0021] In practice, the terms in γ² and γ³ are small and may be neglected compared with the term in γ. The two equations above may then be written in respectively simplified form:
f=f₀+kγ+a(T−T_o)²
and f=f₀+kγ+a(T−T_o)²+b(T−T_o)³

[0022] The curve corresponding to the function f(T) is a parabola or a cubic, respectively, whose local extremum (situated in the bracket of the implementation temperatures) corresponds to the reversal temperature (or “turnover”) temperature T_o. Thus, with each value f of the frequency may be associated two possible values of temperature T_a and T_b, just one of which is the exact temperature of the vibrating resonator. In the case of a third-order polynomial, the third solution is spurious with respect to the temperature domain considered and may be rejected outright. The determination of the temperature of the resonators calls upon a probabilistic estimation and consists in associating, with the sum f₁+f₂, two values T_a and T_b, in calculating the probability P(T_a) that the temperature of the resonators is T_a and the probability P(T_b) that the temperature of the resonators is T_b, then in calculating a weighted value of the temperature of the resonators
T=P(T_a)×T_a+P(T_b)×T_b.

[0023] In the case of accelerometric and/or gyroscopic measurements in several axes (in practice in two or three axes of a reference system) with the aid of several respective accelerometric and/or gyroscopic and/or gyrometric vibrating devices (inertial sensor block), provision may advantageously be made for a value of the temperature-compensated acceleration to be determined, in each axis, according to one or other of the processes set forth above and, in respect of the determination of one at least of the values of the temperature-compensated acceleration, for the weighted value of temperature determined for the vibrating device associated with this axis together with at least one weighted value of temperature determined for at least one vibrating device associated with at least one other axis to be used as weighted temperature value.

[0024] Such an implementation is most particularly beneficial in the case where all the vibrating resonators of the accelerometric and/or gyroscopic and/or gyrometric devices are machined from one and the same block of crystal.

[0034] The invention is aimed at determining a value, corrected as a function of temperature, of the acceleration γ undergone by an accelerometer device having two vibrating resonators of similar geometries, working respectively in stress states of opposite signs. In the example illustrated, the two resonators respectively comprise two parallel and substantially identical vibrating beams working respectively in particular in tension and in compression, but possibly being combinable with at least one other state of stress such as flexion and/or torsion. In the example illustrated, the acceleration γ or the component of the acceleration which is measured is directed substantially perpendicularly to the beams.

[0035] The frequencies of vibration of the two beams may be written, as a function of temperature T and of acceleration γ: $f_1=f_{01}+\prod _1^1\left({\gamma }^n\right)+\prod _1^2\left(T^n\right)$ $f_2=f_{02}+\prod _2^1\left({\gamma }^n\right)+\prod _2^2\left(T^n\right)$

[0036] According to the invention, the frequencies f₁ and f₂ of the vibration of the two vibrating beams are determined, and then the sum f₁+f₂ is calculated.

[0037] On the basis of a pre-established and stored correlation between the values of frequency and the corresponding values of temperature f→π(Tn), the n possible values Ti of temperature are determined, for the frequency f1+f2.

[0038] An evaluation of the estimated temperature of the resonators is then carried out. For this purpose, the value G_mes of a magnitude representative of the temperature in proximity to the device having two vibrating resonators is measured; then, with the aid of a previously established transfer function H, an estimated value G_estimated=H (G_mes) of abovesaid magnitude G_mes is determined, this estimated magnitude being representative of the estimated temperature of the vibrating resonators. On the basis of the estimated value G_estimated of said representative magnitude of the estimated temperature of the resonators and of the possible temperature values T_i, the probabilities P(T_i) that each possible temperature value T_i is the exact value of the temperature of the resonators are determined.

[0039] A simple way of implementing the above provisions consists in performing a temperature measurement directly, in which case the magnitude G representative of the temperature is the temperature itself. It is in this context that what follows is set forth.

[0040] A temperature sensor that delivers a measured temperature value T_mes is disposed in immediate proximity to the box of the accelerometer device and either inside, or more readily outside the latter. By invoking a transfer function H predetermined through earlier trials, an estimated temperature of the vibrating resonators can be evaluated T_estimated=H(T_mes).

[0041] The probability P(T_i) that the temperature of the resonators is T_i is determined on the basis of the values T_estimated and T_i, for each value T_i. Then a weighted value T of the temperature is determined
T=ΣP(T_i)·T_i
that is taken as temperature value of the two resonators.

[0042] Thereafter, knowing the temperature T, it is possible to calculate the value of the acceleration γ on the basis of one and/or more measured frequencies f₁ and f₂.

[0043] The implementation of a polynomial π(Tn) of high order n entails complex and lengthy calculations, and makes it difficult to establish the correlation f→π(Tn) by reason of the high number of possible values Ti associated with each frequency value.

[0044] Accuracy that is acceptable in practice may be obtained with a polynomial of order 2: π(T²), or preferably of order 3: π(T³).

[0045] The vibration frequency is expressed, as a function of acceleration and temperature, by the relation
f=f₀+π(γ³)+a(T−T_o)²
in the case of a polynomial in T of order 2 or by the relation
f=f₀+π(γ³)+a(T−T_o)²+b(T−T_o)³
in the case of a polynomial in T of order 3.

[0046] The terms in γ² and γ³ being small in practice may be neglected compared with the term in γ so that the two aforesaid relations may be written in simplified form, respectively:
f=f₀+kγ+a(T−T_o)²
and f=f₀+kγ+a(T−T_o)²+b(T−T_o)³

[0047] As represented schematically in FIG. 3, the frequencies f₁ and f₂ of the two vibrating resonators P₁ and P₂ subjected to an acceleration γ are measured at 10 and 11, given that the frequencies f₁ and f₂ are dependent on γ and on the temperature T, the case of a polynomial of order 2 in T:
f₁=f₀₁+k₁γ+a₁(T−T_o1)²
f₂=f₀₂−k₂γ+a₂(T−T_o2)²
(where T_o1 and T_o2 are the respective reversal temperatures (or “turnover” temperatures)), or else, in the case of a polynomial of order 3 in T:
f₁=f₀₁+k₁γ+a₁(T−T_o1)²+b₁(T−T_o1)³
f₂=f₀₂−k₂γ+a₂(T−T_o2)²+b₂(T−T_o2)³

[0048] It is noted that the sum f₁+f₂ is independent of the acceleration γ if k₁=k₂ and depends only on the temperature. Moreover, the difference f₁−f₂ depends essentially on the acceleration γ if the curvature terms a₁ and a₂, as well as the reversal values T_o1 and T_o2 are close to one another, respectively (the terms in T² and T³ dropping out). The method in accordance with the invention and which is set forth hereinbelow with regard to FIG. 3A (case of a polynomial of order 2 in T) and FIG. 3B (case of a polynomial of order 3 in T) is based on the exploitation of the two values f₁+f₂ and f₁−f₂.

[0049] The sum f₁+f₂ is determined at 12 on the basis of the measurements of the frequency f₁ at 10 and of the frequency f₂ at 11.

[0050] The two possible values of temperature T_a, T_b corresponding to the value f₁+f₂ emanating from the measurements f₁ and f₂ are plotted at 13 on the basis of a pre-established correlation, held in memory at 14, between the frequency and the temperature of the resonators (polynomial in T³ and/or T² and T). In the case of a polynomial of order 2 in T (case illustrated in FIG. 3A), the two sought-after values T_a and T_b are the two solutions of the second-degree equation in T. In the case of a polynomial of order 3 in T (case illustrated in FIG. 3B), the third-degree equation in T has three solutions, one of which, however, is spurious with respect to the temperature domain considered and should be rejected; it is thus advisable to hold in memory at 25 (FIG. 3B) limit values, for example T_min; T_max, defining the temperature domain and to compare at 26 the three solutions of the third-degree equation with said limit values so as to eliminate the spurious solution and retain only the two appropriate solutions T_a, T_b.

[0051] The two temperature values T_a, T_b are the two possible values that may correspond to the frequency f₁+f₂, but just one of the two temperature values is exact in practice.

[0052] To perform as accurate an approximation as possible of the exact temperature of the resonators (the two resonators are assumed to be at the same temperature), a thermal sensor CTH is disposed in immediate proximity to the resonators. This thermal sensor may be disposed inside the box B (as shown dashed in FIG. 2), but it may also be disposed (since this is a simpler implementation) outside this box as illustrated in solid line in FIG. 2.

[0053] This thermal sensor CTH provides at 15 (FIGS. 3A and 3B) a measured value of temperature T_mes, which value T_mes is only an approximate value, and not exact, of the temperature of the resonators P₁, P₂.

[0054] Prior to the measurement, a transfer function H is established experimentally between the exact values of temperature of the resonators P₁, P₂ situated inside the box B and the temperature values measured, simultaneously, by the sensor CTH outside the box B. The transfer function H is held in memory at 16.

[0055] To fix matters, the transfer function H may be a simple function of second or third order, such as, with H(p)=T(p)/T_CHT(P)
H(p)=1/(1+αp)²
or H(p)=1/[(1+α₁p)(1+α₂p)(1+α₃p)]
in which p is the Laplace operator and α, α₁, α₂, α₃ are constants depending on the structure of the temperature sensor and on the time for the heat to propagate between the vibrating beams situated inside the box B and the temperature sensor CTH situated outside this box.

[0056] Thus, with the temperature T_mes measured outside the box B, it is possible to associate at 17, through the transfer function H, an estimated value T_est of the temperature of the beams P₁, P₂ inside the box B.

[0057] Then, on the basis of the values T_a, T_b and T_est, probability calculations are performed at 18 giving the probability P(T_a) that the temperature T_a is the temperature of the beams P₁, P₂ and the probability P(T_b) at the temperature T_b is the temperature of the beams P₁, P₂.

[0058] To fix matters, the probabilities P(T_a) and P(T_b) may be calculated in the following manner:
P(T_a)=g(T_a)/(g(T_a)+g(T_b))
P(T_b)=g(T_b)/(g(T_a)+g(T_b))
with
 g(T_a)=1/(|T_a−T_est|)
g(T_b)=1/(|T_b−T_est|)

[0059] The calculations could also be performed with the following functions:
g(T_a)=exp(−α·|T_a−T_est|)
g(T_b)=exp(−α·|T_b−T_est|)
in which α is a constant.

[0060] The formulae indicated here have the advantage of being simple and of leading to fast calculations. However, if necessary, the person skilled in the art may have recourse to formulae that are more complex and more accurate, but also entail lengthier calculations.

[0061] Then, at 19, is calculated a value T of the temperature which is a weighting of the possible values T_a and T_b afforded the previously calculated weighting coefficients P(T_a) and P(T_b):
T=P(T_a)·T_a+P(T_b)·T_b.

[0062] A value T of temperature of the resonators which is statistically the most likely is thus defined.

[0063] Finally, the temperature-compensated value of the acceleration γ undergone by the accelerometer device is determined on the basis of one of the measured frequencies f₁ or f₂, or on the basis of the two measured frequencies f₁ and f₂, and of the weighted value T of the temperature that has just been determined at 19.

[0064] In practice, the quantity f₁−f₂ is calculated at 20. It is known that, in the case of the polynomial of order 2 in T, the mathematical expression therefor is:
f₁−f₂=(f₀₁−f₀₂)+(k₁+k₂)γ+a₁(T−T_o1)²−a₂(T−T_o2)

[0065] Stated otherwise, this expression is of the type:

 f₁−f₂=p+qγ+r T+sT²

[0066] On the basis of the value T of the temperature determined at 19, it is then easy to calculate, at 21, the value γ from T and from f₁−f₂.

[0067] In the case of a polynomial of order 3 in T, the mathematical expression for f₁−f₂ is more complex, but can nevertheless be determined and also allows calculation of the value of γ on the basis of the calculation of the value f₁−f₂.

[0068] A temperature-compensated value of the acceleration γ is ultimately obtained at 22.

[0069] The correlation between the frequency and the temperature (T³ and/or T², T), held in memory at 14, may be stored in the form of a mathematical equation of second or third degree in T whose coefficients are pre-established by prior experimental measurements performed on the accelerometer device 1. The two possible values T_a and T_b are then calculated by solving this equation for each given value of f₁+f₂. An array giving, for all the possible values f₁+f₂, the two corresponding possible values T_a, T_b may also be held in memory at 14.

[0070] In the same way, the transfer function H held in memory at 16 may be stored either in the form of an algebraic relation T_est=H(T_mes), or in the form of an array of pre-established data giving, for each value T_mes, the associated value T_est.

[0071] Experience has shown that the curve that mathematically conveys the relation between the frequency and the temperature (T³ and/or T², T) was prone to variation over time, this variation being conveyed essentially as a drift in the frequency for a given temperature, the shape of the curve f(T³ and/or T², T) remaining substantially unchanged.

[0072] To obtain a method of calculating the acceleration γ which remains accurate over time, it is very desirable to provide for a recalibration as illustrated in FIG. 4 of the appended drawings.

[0073] With this aim, the steps described previously are supplemented with the following steps appearing in FIG. 4 (which is derived more precisely from FIG. 3A). The derivative with respect to time dT_mes/dt is calculated at 22 on the basis of the measured value of the temperature T_mes. If, at 23, this drift is detected to be zero or almost zero, then a recalibration procedure RECAL is instigated at 24, able to act on the transfer function held in memory at 16, on the calculation/selection at 13 of the possible values of temperature T_a, T_b, and on the calculation 21 of the value f₁−f₂=p+q γ+r T leading, in a simple way, to the determination of the acceleration γ.

[0074] Finally, it will be noted that the acceleration γ is calculated iteratively, by successive calculation passes, and that each calculation pass takes into account the value of the acceleration found at the conclusion of the previous calculation pass, in particular for the compensation of the sum f₁+f₂ mentioned above when k₁≠k₂.

[0075] In practice, several vibrating devices may be grouped together to provide accelerometric and/or gyrometric or gyroscopic information in two or three reference axes (for example inertial sensor block), all these vibrating devices calling individually upon a temperature measurement for the determination of the compensated values of the vibration frequencies.

[0076] According to the invention, the temperature-compensated accelerometric and/or gyroscopic and/or gyrometric measurement, performed for each device having two vibrating resonators, takes into account not only its own weighted temperature value, but also at least certain weighted temperatures predetermined for other vibrating devices so as to obtain an interpolated frequency value related to at least a certain number of the measured temperatures.

[0077] By way of example, represented very schematically in FIG. 5 is an inertial sensor 28 grouping together three vibrating devices for accelerometric measurements in three axes x, y, z of a reference frame, each vibrating device possibly being of the type illustrated in FIG. 2 and being referenced 1¹_x, 1¹_yand 1¹_z respectively, as well as three vibrating devices for gyroscopic measurements with respect to the three axes x, y, z of the reference frame, each vibrating device possibly being of the type illustrated in FIG. 2 and being referenced 1²_x, 1²_yand 1²_z respectively.

[0078] With each aforesaid vibrating device is associated an information processing unit that receives the information f₁, f₂ and T_mes provided by each device. As illustrated in FIG. 6, each unit 27 is arranged (by way of example) in accordance with the representation given in FIG. 3A. FIG. 6 shows, by way of example, the unit 27¹_x associated with the vibrating device 1¹_x, it being understood that the other units associated respectively with the other vibrating devices are identical to it.

[0079] All the weighted temperature values determined by all the units 27, T¹_x, T¹_y, T¹_z, T²_x, T²_y, T²_z are sent to a correlation device 29 which, on the basis of a pre- established correlation (for example calculation of the mean), delivers a correlated value of temperature T_c, which is thereafter delivered to each unit 27.

[0080] In the unit 27, the correlated value T_c may then be used directly at 21 to calculate, from f₁−f₂, the value of the component of acceleration (γ¹_x in FIG. 6).

[0081] It is also conceivable, as illustrated in FIG. 7 which shows only a part of the unit 27¹_x for the correlated value of temperature T_c to be combined, according to a pre-established law, at 30, with the weighted value T¹_x the resulting value being used for the calculation of the component γ¹_x from the value f₁−f₂.

[0082] Of course, numerous variants of temperature correlation are conceivable, for example: correlation in respect of the temperatures of the accelerometric vibrating devices alone and correlation in respect of the temperatures of the gyroscopic vibrating devices alone; cross-correlation between the temperatures of the accelerometric and gyroscopic devices; etc.

[0083] An economic advantage may be obtained when all the vibrating devices are machined from one single block of crystal or “wafer”. The vibrating devices are then situated in immediate proximity to one another and there is no longer any need to assign a temperature sensor to each of them: according to the configuration of mutual layout of the devices, a single temperature sensor disposed centrally may suffice, or at the very least two temperature sensors may be envisaged, one associated with the vibrating devices for accelerometric measurements and the other with the vibrating devices for gyroscopic or gyrometric measurements. A gain in dimension in respect of the substrate and a saving as regards the information processing capacity is thus achieved.

[0084] It is in particular in the case of the implementation of several adjacently located vibrating devices that it may be appropriate to process a magnitude G representative of the temperature without calling upon the temperature itself. Such a case may arise for example when the device 1 having two vibrating resonators of frequency f₁ and f₂ is supplemented with an adjacently adjoining vibrating device such as a gyroscopic resonator of the so-called “bowl” type, whose resonant frequency varies linearly as a function of temperature and may be written f=mt+n where m and n are constants. It is then understood that measuring the frequency f of this vibrating device provides a temperature indication which, to within a constant, is equivalent to a direct measurement of the temperature with the aid of a temperature sensor.

[0085] Hence, the processing of aforesaid magnitude G(T), if the latter varies in a simple manner as a function of temperature, makes it possible to proceed with the method of accelerometric measurement according to the invention without it being necessary to implement an actual temperature sensor.

[0086] It is an arrangement of this type that is illustrated in FIG. 8, the diagram of which mimics that of FIG. 3A. At 15 is measured the value G_mes of a magnitude G(T) in particular in a vibrating device adjacent to the device 1 having two vibrating resonators of frequency f₁ and f₂, said magnitude G(T) thus being representative of the temperature in proximity to said vibrating device 1. Then, with the aid of a transfer function H previously established at 16, an estimated value G_estimated=H(G_mes) of the abovementioned magnitude G(T) is determined, the estimated value G_estimated being representative of the estimated temperature of the vibrating resonators of frequency f₁ and f₂. It is this value G_estimated that is used thereafter, in conjunction with the values T_a, T_b, to determine the probabilities P(T_a) and P(T_b) mentioned above, without in fact involving the value of the estimated temperature T_estimated.

[0087] Such a solution has the benefit of avoiding the implementation of one or more temperature sensors, this being advantageous from the economic standpoint. Furthermore, the eliminating of one or more hardware components in an inertial sensor block is always particularly welcome by reason of the gain of space and of bulk thus obtained which goes hand in hand with the search for ever more compact inertial sensor blocks.