CN114780903B - Hemispherical resonator gyroscope X/Y channel drive gain asymmetry and electrode non-orthogonal error identification method - Google Patents

Hemispherical resonator gyroscope X/Y channel drive gain asymmetry and electrode non-orthogonal error identification method Download PDF

Info

Publication number
CN114780903B
CN114780903B CN202210480919.XA CN202210480919A CN114780903B CN 114780903 B CN114780903 B CN 114780903B CN 202210480919 A CN202210480919 A CN 202210480919A CN 114780903 B CN114780903 B CN 114780903B
Authority
CN
China
Prior art keywords
electrode
orthogonal
gain
asymmetry
driving
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202210480919.XA
Other languages
Chinese (zh)
Other versions
CN114780903A (en
Inventor
解伟男
王奇
奚伯齐
孙一为
伊国兴
王常虹
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN202210480919.XA priority Critical patent/CN114780903B/en
Publication of CN114780903A publication Critical patent/CN114780903A/en
Application granted granted Critical
Publication of CN114780903B publication Critical patent/CN114780903B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C19/00—Gyroscopes; Turn-sensitive devices using vibrating masses; Turn-sensitive devices without moving masses; Measuring angular rate using gyroscopic effects
    • G01C19/56—Turn-sensitive devices using vibrating masses, e.g. vibratory angular rate sensors based on Coriolis forces
    • G01C19/567—Turn-sensitive devices using vibrating masses, e.g. vibratory angular rate sensors based on Coriolis forces using the phase shift of a vibration node or antinode
    • G01C19/5691—Turn-sensitive devices using vibrating masses, e.g. vibratory angular rate sensors based on Coriolis forces using the phase shift of a vibration node or antinode of essentially three-dimensional vibrators, e.g. wine glass-type vibrators

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Algebra (AREA)
  • Mathematical Optimization (AREA)
  • Remote Sensing (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Gyroscopes (AREA)

Abstract

A method for identifying the asymmetry of X/Y channel drive gain and non-orthogonal error of electrode of hemispherical resonator gyroscope belongs to the technical field of identification of error parameters of hemispherical resonator gyroscopes. The invention solves the problems of standing wave azimuth drift and control loop error caused by inconsistent gains of two X/Y driving signals and non-orthogonality of driving electrodes. The method is based on the existence of a drive gain asymmetric coefficient k and an electrode non-orthogonal deflection angle
Figure DDA0003627818440000011
The time precession speed model identifies the error parameters by using an identification method, so that the driving gain asymmetric error coefficient and the electrode non-orthogonal deflection angle are calculated, the problems of standing wave azimuth drift and control loop errors are solved, and the performance of the hemispherical resonator gyroscope is improved. The method can be applied to identifying the asymmetry of the X/Y channel drive gain of the half-sphere resonance gyroscope and the non-orthogonal error of the electrode.

Description

Hemispherical resonator gyroscope X/Y channel drive gain asymmetry and electrode non-orthogonal error identification method
Technical Field
The invention belongs to the technical field of hemispherical resonator gyroscope error parameter identification, and particularly relates to a hemispherical resonator gyroscope X/Y path drive gain asymmetry and electrode non-orthogonal error identification method.
Background
The hemispherical resonator gyroscope is a new generation high-precision gyroscope developed on the basis of the traditional mechanical rotor gyroscope and optical gyroscope, is a high-precision gyroscope with inertial navigation level performance in a Coriolis vibration gyroscope, is one of the mainstream high-precision inertial devices at present, and is widely applied to the fields of aviation, aerospace, navigation and the like. The hemispherical resonator gyroscope detects the position of the standing wave in real time through the X/Y signals, and then the external input angle and the angular speed can be measured. When the hemispherical resonator gyroscope works normally, an amplitude control loop is needed to stabilize the vibration amplitude of the harmonic oscillator, and an orthogonal control loop is needed to reduce the amplitude of orthogonal waves. The gyro parameters are influenced by factors such as parameter matching conditions, mechanical errors and environment of the X/Y driving circuits, so that gains of the X/Y driving signals are inconsistent and driving electrodes are not orthogonal, problems such as standing wave azimuth drift and control loop errors are caused, and the performance of the gyro is reduced. Therefore, it is very significant to provide a method for identifying the asymmetry of the driving gain of the X/Y path of the hemispherical resonator gyroscope and the non-orthogonal error of the electrode.
Disclosure of Invention
The invention aims to solve the problems of standing wave azimuth drift and control loop error caused by inconsistent gain of two X/Y driving signals and non-orthogonality of driving electrodes, and provides a method for identifying the asymmetry of the X/Y channel driving gain of a hemispherical resonator gyroscope and the non-orthogonality of the electrodes.
The technical scheme adopted by the invention for solving the technical problems is as follows:
a method for identifying the asymmetry of X/Y channel drive gain and non-orthogonal error of an electrode of a hemispherical resonator gyroscope specifically comprises the following steps:
step 1, solving out precession speed omega according to standing wave azimuth angle theta and vibration amplitude A r (ii) a The specific process comprises the following steps:
Figure BDA0003627818420000011
wherein, F x 、F y Electrostatic forces applied to the x and y axis electrodes, respectively, when k is x ≠k y And the y-axis electrode has an angle of deflection in the direction orthogonal to the x-axis
Figure BDA0003627818420000012
When F is present x 、F y The expression of (a) is:
Figure BDA0003627818420000013
wherein k is x 、k y Respectively x-axis and y-axis drive circuit gain coefficients, G is harmonic oscillator kinematic model, omega 0 A given precessional rotation speed;
step 2, define H = k 0 Gω 0 ，k 0 The gain coefficients of the ideal x-axis and y-axis driving circuits;
and giving an initial value k (0) of the driving gain asymmetry coefficient k and the non-orthogonal deflection angle of the electrode
Figure BDA0003627818420000021
Initial value of
Figure BDA0003627818420000022
And an initial value of H, H (0);
step 3, according to P r ＝Aω r To obtain
Figure BDA0003627818420000023
Then according to P r Expression of the calc function r (i):
r(i)＝P d (i)-P r (i)
wherein, P r (i) Is P r Calculated value at time i, P d (i) R (i) is a value function of the moment i, and is an actual detection value obtained after processing a signal collected at the moment i;
step 4, calculating the Jacobian matrix J of the value function r (i) r (i) Then according to the Jacobian matrix J r (i) Calculating increment of driving gain asymmetry coefficient k at moment i and non-orthogonal deflection angle of electrode
Figure BDA0003627818420000024
And increments of H;
step 5, driving gain asymmetry coefficient k and electrode non-orthogonal deflection angle of each increment pair i +1 moment calculated in the step 4
Figure BDA0003627818420000025
And H, updating;
and 6, repeatedly executing the processes from the step 3 to the step 5 until no control system variable is input, and taking the driving gain asymmetric coefficient and the electrode non-orthogonal deflection angle obtained by the last iteration as final outputs, namely obtaining the identification results of the driving gain asymmetric coefficient and the electrode non-orthogonal deflection angle.
Further, the Jacobian matrix J r (i) The expression of (a) is:
Figure BDA0003627818420000026
further, the function is according to a Jacobian matrix J r (i) Calculating increment of driving gain asymmetry coefficient k at moment i and non-orthogonal deflection angle of electrode
Figure BDA0003627818420000027
The increment of (c) and the increment of H are specifically:
Figure BDA0003627818420000028
wherein, the delta k (i) is the increment of the driving gain asymmetry coefficient k at the moment i,
Figure BDA0003627818420000029
for non-orthogonal deflection angles of electrodes
Figure BDA00036278184200000210
The increment at time i, Δ H (i), is the increment of H at time i, the superscript T represents the transpose of the matrix, and the superscript-1 represents the inverse of the matrix.
Further, the specific process of step 5 is as follows:
Figure BDA0003627818420000031
wherein k (i + 1),
Figure BDA0003627818420000032
H (i + 1) is the driving gain asymmetry coefficient k and the non-orthogonal deflection angle of the electrode
Figure BDA0003627818420000033
And the value of H at time i + 1.
Further, an initial value k (0) =1 of the driving gain asymmetry coefficient k.
Further, the electrodes are non-orthogonally off-angled
Figure BDA0003627818420000034
Initial value of
Figure BDA0003627818420000035
Further, an initial value H (0) =0 of H.
The invention has the beneficial effects that:
the method is based on the existence of a drive gain asymmetric coefficient k and an electrode non-orthogonal deflection angle
Figure BDA0003627818420000036
And the time precession speed model identifies the error parameters by using an identification method, so that the asymmetric error coefficient of the driving gain and the non-orthogonal deflection angle of the electrode are calculated, the problems of standing wave azimuth drift and control loop error are solved, and the performance of the hemispherical resonant gyroscope is improved.
Drawings
FIG. 1 is a graph of a driving gain asymmetry coefficient k identification;
FIG. 2 shows non-orthogonal deflection angles of electrodes
Figure BDA0003627818420000037
The graph is identified.
Detailed Description
The first embodiment is as follows: the method for identifying the X/Y channel drive gain asymmetry and the electrode non-orthogonal error of the hemispherical resonator gyroscope in the embodiment comprises the following specific processes:
step 1, solving out precession speed omega according to standing wave azimuth angle theta and vibration amplitude A r (ii) a The specific process comprises the following steps:
Figure BDA0003627818420000038
wherein, F x 、F y Electrostatic forces applied to the x and y axis electrodes, respectively, when k is x ≠k y And the y-axis electrode has a deflection angle with respect to the x-axis orthogonal direction
Figure BDA0003627818420000039
When F is present x 、F y The expression of (a) is:
Figure BDA00036278184200000310
wherein k is x 、k y Respectively, the gain coefficients of the driving circuit in x and y axes, G is a harmonic oscillator kinematic model, omega 0 A given precessional rotation speed;
step 2, define H = k 0 Gω 0 ，k 0 The gain coefficients of the ideal x-axis and y-axis driving circuits;
and giving an initial value k (0) of the driving gain asymmetry coefficient k and the non-orthogonal deflection angle of the electrode
Figure BDA0003627818420000041
Initial value of
Figure BDA0003627818420000042
And an initial value of H, H (0);
usually, the gain deviation value of the electrode is less than 10 percent, the non-orthogonal deflection angle of the electrode is less than 10 degrees, and H is invariable in the precession process;
wherein G is harmonic oscillator kinematic model, k 0 Is the gain coefficient, omega, of an ideal x-axis and y-axis drive circuit 0 A given precessional rotation speed;
step (ii) of3. According to P r ＝Aω r To obtain
Figure BDA0003627818420000043
Then according to P r Expression of the calc function r (i):
r(i)＝P d (i)-P r (i)
wherein, P r (i) Is P r Calculated value at time i, P d (i) R (i) is a value function of the i moment, and is an actual detection value obtained after the signal acquired at the i moment is processed;
step 4, calculating the Jacobian matrix J of the value function r (i) r (i) Then according to the Jacobian matrix J r (i) Calculating increment of driving gain asymmetry coefficient k at moment i and non-orthogonal deflection angle of electrode
Figure BDA0003627818420000045
And increments of H;
step 5, according to the driving gain asymmetry coefficient k of each increment calculated in the step 4 to the i +1 moment and the non-orthogonal deflection angle of the electrode
Figure BDA0003627818420000046
And H, updating;
and 6, repeatedly executing the processes from the step 3 to the step 5 until no control system variable is input, and taking the driving gain asymmetric coefficient and the electrode non-orthogonal deflection angle obtained by the last iteration as final outputs, namely obtaining the identification results of the driving gain asymmetric coefficient and the electrode non-orthogonal deflection angle.
In the embodiment, firstly, the hemispherical harmonic oscillator is subjected to parameter excitation, so that the vibration amplitude of the harmonic oscillator is stable; applying a constant electrostatic force along the direction vertical to the amplitude axis of the ideal standing wave to make the standing wave precess; then collecting parameters C related to harmonic oscillator standing wave vibration x 、S x 、C y 、S y And performing secondary combination to obtain a hemispherical resonator gyro control system variable E, Q, S, R signal.
Solving a standing wave azimuth angle theta and a vibration amplitude A according to the control system variables E, Q, S and R;
Figure BDA0003627818420000044
Figure BDA0003627818420000051
the second embodiment, which is different from the first embodiment, is: the Jacobian matrix J r (i) The expression of (a) is:
Figure BDA0003627818420000052
other steps and parameters are the same as those in the first embodiment.
The third embodiment, which is different from the first or second embodiment, is: the Jacobian matrix J r (i) Calculating increment of driving gain asymmetry coefficient k at moment i and non-orthogonal deflection angle of electrode
Figure BDA0003627818420000053
The increment of (c) and the increment of H are specifically:
Figure BDA0003627818420000054
wherein, the delta k (i) is the increment of the driving gain asymmetry coefficient k at the moment i,
Figure BDA0003627818420000055
for non-orthogonal deflection angles of electrodes
Figure BDA0003627818420000056
The increment at time i, Δ H (i), is the increment of H at time i, the superscript T represents the transpose of the matrix, and the superscript-1 represents the inverse of the matrix.
Other steps and parameters are the same as those in the first or second embodiment.
The fourth embodiment and the differences between this embodiment and the first to the third embodiments are: the specific process of the step 5 is as follows:
Figure BDA0003627818420000057
wherein k (i + 1),
Figure BDA0003627818420000058
H (i + 1) is the driving gain asymmetry coefficient k and the non-orthogonal deflection angle of the electrode
Figure BDA0003627818420000059
And the value of H at time i + 1.
Other steps and parameters are the same as those in one of the first to third embodiments.
The fifth embodiment is different from the first to the fourth embodiments in that: an initial value k (0) =1 of the driving gain asymmetry coefficient k.
Other steps and parameters are the same as in one of the first to fourth embodiments.
Sixth embodiment, the difference between this embodiment and one of the first to fifth embodiments, is: the electrode non-orthogonal deflection angle
Figure BDA00036278184200000510
Initial value of (2)
Figure BDA00036278184200000511
Other steps and parameters are the same as those in one of the first to fifth embodiments.
The seventh embodiment and the differences between the first to sixth embodiments are as follows: an initial value H (0) =0 of the H.
Other steps and parameters are the same as those in one of the first to sixth embodiments.
Examples
The process of the invention implemented by taking a simulation experiment as an example is as follows:
step 1, performing parameter excitation on a hemispherical harmonic oscillator to stabilize the vibration amplitude of the harmonic oscillator;
step 2, setting the sampling frequency to be f s =100Hz, sampling time t =100s, drive electrode gain ratio k =1.03, and electrode non-orthogonal deflection angle
Figure BDA0003627818420000061
Applying a constant electrostatic force along the direction vertical to the amplitude axis of the ideal standing wave to make the standing wave precess;
step 3, collecting parameters C related to harmonic oscillator standing wave vibration x 、S x 、C y 、S y And performing secondary combination to obtain a hemispherical resonator gyro control system variable E, Q, S, R, and calculating a standing wave azimuth angle theta and a precession speed omega r And vibration amplitude A:
Figure BDA0003627818420000062
Figure BDA0003627818420000063
Figure BDA0003627818420000064
step 4, since the nonlinear least square identification method is a recursive algorithm, an initial estimation value must be given first
Figure BDA0003627818420000065
Figure BDA0003627818420000066
Normally, the gain deviation value of the electrode is less than 10%, the non-orthogonal deviation angle of the electrode is less than 10 degrees, and H is an invariant in the precession process. Therefore, the least square method is given to identify the initial value
Figure BDA0003627818420000067
H＝k 0 Gω 0
Figure BDA0003627818420000068
Wherein G is harmonic oscillator kinematic model, k 0 Is the gain coefficient, omega, of an ideal x-axis and y-axis drive circuit 0 For a given precessional speed, k x 、k y Respectively x-axis and y-axis drive circuit gain coefficients, k is a drive gain asymmetry coefficient,
Figure BDA0003627818420000069
is the non-orthogonal deflection angle of the electrode;
step 5, calculating a current time value function:
when the included angle between the amplitude axis of the standing wave and the x-axis electrode is theta, the parameter equation for driving the static electricity force of the standing wave precession is as follows:
Figure BDA00036278184200000610
wherein, F x 、F y Electrostatic forces applied to the x and y axis electrodes, respectively.
When k is 0 ＝k x ≠k y And the y-axis electrode has an angle of deflection in the direction orthogonal to the x-axis
Figure BDA0003627818420000071
The parameter equation of the static electricity force for driving the standing wave precession is as follows:
Figure BDA0003627818420000072
the precession speed perpendicular to the amplitude axis direction of the standing wave is as follows:
Figure BDA0003627818420000073
simplifying to obtain:
Figure BDA0003627818420000074
the value function at this time: r (i) = P d (i)-P r (i)
Wherein, P = A ω r ，P d For acquiring actual detection values after signal processing, P r Calculating a value according to a theoretical model;
step 6, calculating a Jacobian matrix of the function at the current moment:
Figure BDA0003627818420000075
step 7, calculating the increment of the identification parameter at the current moment:
Figure BDA0003627818420000076
step 8, updating the identification parameters at the next moment:
Figure BDA0003627818420000077
step 9, judging whether a hemispherical resonator gyro control system variable signal E, Q, S, R is input, if so, jumping to step 5, and if not, jumping to step 10;
step 10, outputting the driving gain coefficient
Figure BDA0003627818420000078
In conclusion, the asymmetry of the X/Y channel driving gain of the half-sphere resonance gyroscope and the identification of the non-orthogonal error of the electrode are realized.
The simulation results are shown in FIG. 1 and FIG. 2, from which the final identification can be obtained
Figure BDA0003627818420000079
The method of the invention has high identification precision.
The above-described calculation examples of the present invention are merely to describe the calculation model and the calculation flow of the present invention in detail, and are not intended to limit the embodiments of the present invention. It will be apparent to those skilled in the art that other variations and modifications of the present invention can be made based on the above description, and it is not intended to be exhaustive or to limit the invention to the precise form disclosed, and all such modifications and variations are possible and contemplated as falling within the scope of the invention.

Claims (5)

1. A method for identifying the asymmetry of X/Y channel drive gain and non-orthogonal error of an electrode of a hemispherical resonator gyroscope is characterized by comprising the following steps:
step 1, solving out precession speed omega according to standing wave azimuth angle theta and vibration amplitude A r (ii) a The specific process comprises the following steps:
Figure FDA0003905212500000011
wherein, F x 、F y Electrostatic forces applied to the x and y axis electrodes, respectively, when k is x ≠k y And the non-orthogonal off-angle of the electrode is
Figure FDA0003905212500000012
When F is present x 、F y The expression of (a) is:
Figure FDA0003905212500000013
wherein k is x 、k y Respectively x-axis and y-axis drive circuit gain coefficients, G is harmonic oscillator kinematic model, omega 0 A given precessional rotation speed;
step (ii) of2. Definition H = k 0 Gω 0 ，k 0 The gain coefficients of the ideal x-axis and y-axis driving circuits;
and giving an initial value k (0) of the driving gain asymmetry coefficient k and the non-orthogonal deflection angle of the electrode
Figure FDA0003905212500000014
Initial value of
Figure FDA0003905212500000015
And an initial value of H, H (0);
step 3, according to P r ＝Aω r To obtain
Figure FDA0003905212500000016
Then according to P r Expression of the calc function r (i):
r(i)＝P d (i)-P r (i)
wherein, P r (i) Is P r Calculated value at time i, P d (i) R (i) is a value function of the moment i, and is an actual detection value obtained after processing a signal collected at the moment i;
step 4, calculating the Jacobian matrix J of the value function r (i) r (i) Then according to the Jacobian matrix J r (i) Calculating increment of driving gain asymmetry coefficient k at moment i and non-orthogonal deflection angle of electrode
Figure FDA0003905212500000017
And increments of H;
the Jacobian matrix J r (i) The expression of (a) is:
Figure FDA0003905212500000018
the Jacobian matrix J r (i) Calculating increment of driving gain asymmetry coefficient k at moment i and non-orthogonal deflection angle of electrode
Figure FDA0003905212500000021
The increment of (c) and the increment of H are specifically:
Figure FDA0003905212500000022
wherein Δ k (i) is the increment of the driving gain asymmetry coefficient k at the moment i,
Figure FDA0003905212500000023
for non-orthogonal deflection angles of electrodes
Figure FDA0003905212500000024
Delta H (i) is the increment of H at the moment i, an upper superscript T represents the transposition of the matrix, and an upper superscript-1 represents the inverse of the matrix;
step 5, according to the driving gain asymmetry coefficient k of each increment calculated in the step 4 to the i +1 moment and the non-orthogonal deflection angle of the electrode
Figure FDA0003905212500000025
And H, updating;
and 6, repeatedly executing the processes from the step 3 to the step 5 until no control system variable is input, and taking the driving gain asymmetric coefficient and the electrode non-orthogonal deflection angle obtained by the last iteration as final outputs, namely obtaining the identification results of the driving gain asymmetric coefficient and the electrode non-orthogonal deflection angle.
2. The method for identifying the asymmetry of the driving gain of the X/Y path of the hemispherical resonator gyroscope and the non-orthogonal error of the electrode as claimed in claim 1, wherein the specific process of the step 5 is as follows:
Figure FDA0003905212500000026
wherein k (i + 1),
Figure FDA0003905212500000027
H (i + 1) is the driving gain asymmetry coefficient k and the non-orthogonal deflection angle of the electrode
Figure FDA0003905212500000028
And the value of H at time i + 1.
3. The method for identifying the X/Y channel driving gain asymmetry and the electrode non-orthogonality error of the hemispherical resonator gyroscope of claim 2, wherein an initial value k (0) =1 of the driving gain asymmetry coefficient k.
4. The method as claimed in claim 3, wherein the electrode non-orthogonal declination is determined by the asymmetry of the X/Y channel driving gain of the hemispherical resonator gyroscope
Figure FDA0003905212500000029
Initial value of
Figure FDA00039052125000000210
5. The method as claimed in claim 4, wherein the initial value of H is H (0) =0.
CN202210480919.XA 2022-05-05 2022-05-05 Hemispherical resonator gyroscope X/Y channel drive gain asymmetry and electrode non-orthogonal error identification method Active CN114780903B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210480919.XA CN114780903B (en) 2022-05-05 2022-05-05 Hemispherical resonator gyroscope X/Y channel drive gain asymmetry and electrode non-orthogonal error identification method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210480919.XA CN114780903B (en) 2022-05-05 2022-05-05 Hemispherical resonator gyroscope X/Y channel drive gain asymmetry and electrode non-orthogonal error identification method

Publications (2)

Publication Number Publication Date
CN114780903A CN114780903A (en) 2022-07-22
CN114780903B true CN114780903B (en) 2023-02-07

Family

ID=82435897

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210480919.XA Active CN114780903B (en) 2022-05-05 2022-05-05 Hemispherical resonator gyroscope X/Y channel drive gain asymmetry and electrode non-orthogonal error identification method

Country Status (1)

Country Link
CN (1) CN114780903B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116625410B (en) * 2023-07-20 2023-09-19 中国船舶集团有限公司第七〇七研究所 Hemispherical resonator gyro excitation electrode gain calibration method

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109104904B (en) * 2014-12-10 2016-12-28 上海新跃仪表厂 A kind of hemispherical reso nance gyroscope high-precision driving circuit based on dynamic balance mode
CN108489512A (en) * 2018-02-28 2018-09-04 北京控制工程研究所 A kind of the compensation scaling method and device of hemispherical reso nance gyroscope constant multiplier
CN113670340A (en) * 2021-09-13 2021-11-19 哈尔滨工业大学 Method and system for measuring standing wave azimuth angle of hemispherical resonator gyroscope based on X/Y signal phase difference identification
CN113899365A (en) * 2021-11-22 2022-01-07 哈尔滨工业大学 Standing wave azimuth angle measurement method based on asymmetric parameter identification of hemispherical resonator gyroscope detection channel

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2993738A1 (en) * 2012-07-23 2014-01-24 St Microelectronics Grenoble 2 METHOD AND APPARATUS FOR ESTIMATING GAIN AND PHASE CORRECTIVE PARAMETERS WHEN RECEIVING A MODULE SIGNAL ACCORDING TO OFDM TYPE MODULATION
US20210348503A1 (en) * 2020-05-11 2021-11-11 Institute Of Geology And Geophysics, Chinese Academy Of Sciences Self-calibration method and system of solid-state resonator gyroscope
CN112506040B (en) * 2020-11-16 2024-01-30 哈尔滨工程大学 Single-channel control system and method for full-angle hemispherical resonant gyroscope
CN113587954B (en) * 2021-08-06 2023-05-23 大连海事大学 Compensation control method and system for damping non-uniformity of full-angle hemispherical resonant gyroscope

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109104904B (en) * 2014-12-10 2016-12-28 上海新跃仪表厂 A kind of hemispherical reso nance gyroscope high-precision driving circuit based on dynamic balance mode
CN108489512A (en) * 2018-02-28 2018-09-04 北京控制工程研究所 A kind of the compensation scaling method and device of hemispherical reso nance gyroscope constant multiplier
CN113670340A (en) * 2021-09-13 2021-11-19 哈尔滨工业大学 Method and system for measuring standing wave azimuth angle of hemispherical resonator gyroscope based on X/Y signal phase difference identification
CN113899365A (en) * 2021-11-22 2022-01-07 哈尔滨工业大学 Standing wave azimuth angle measurement method based on asymmetric parameter identification of hemispherical resonator gyroscope detection channel

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
半球谐振陀螺控制及补偿技术;伊国兴;《宇航学报》;20200630;第41卷(第6期);780-789 *

Also Published As

Publication number Publication date
CN114780903A (en) 2022-07-22

Similar Documents

Publication Publication Date Title
CN114440852B (en) Hemispherical resonator gyroscope X/Y path drive gain asymmetry error identification method
CN113670340B (en) Method and system for measuring standing wave azimuth angle of hemispherical resonator gyroscope based on X/Y signal phase difference identification
Wu et al. Generalized linear quaternion complementary filter for attitude estimation from multisensor observations: An optimization approach
CN113899365B (en) Standing wave azimuth angle measurement method based on detection of asymmetric parameter identification of channel of hemispherical resonator gyroscope
US11015934B2 (en) Vibrating-mass gyroscope systems and method
CN105783899B (en) High bandwidth Coriolis Vibration Gyroscope (CVG) with in-situ bias self-calibration
CN111578923A (en) Closed-loop control method and system for resonant gyroscope
CN113686356B (en) Resonant gyro zero bias on-line self-compensation system and method based on RBF network
CN114780903B (en) Hemispherical resonator gyroscope X/Y channel drive gain asymmetry and electrode non-orthogonal error identification method
CN112504298A (en) GNSS-assisted DVL error calibration method
CN108489485B (en) Error-free strapdown inertial navigation value updating method
Hu et al. Control and damping imperfection compensation for a rate integrating MEMS gyroscope
CN115407657A (en) Intelligent control method of hemispherical resonant gyroscope under input saturation
CN113137977B (en) SINS/polarized light combined navigation initial alignment filtering method
Zhang et al. SINS initial alignment based on fifth-degree cubature Kalman filter
Chikovani et al. Digital rate MEMS vibratory gyroscope modeling, tuning and simulation results
Liu et al. Design of a digital closed control loop for the sense mode of a mode-matching MEMS vibratory gyroscope
CN116698084A (en) Real-time zero-bias drift suppression and zero-bias stability improvement method for symmetrical gyroscope
Ma et al. Application and comparison of two methods for alignment of FOG SINS
Shi et al. Event-triggered output feedback control for MEMS gyroscope with prescribed performance
CN114964309A (en) Temperature compensation method for MEMS inertial navigation system
CN113777350A (en) Acceleration sensor data processing method based on stabilized numerical integration
Bo et al. Research on a new decoupled dual-mass micro-gyroscope
CN111256731A (en) Ground static alignment method capable of fast convergence without being influenced by historical data
CN109342767B (en) Method, device and system for testing accelerometer scale factor

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant