CN116576886B - Hemispherical resonator gyro damping non-uniformity identification method - Google Patents

Hemispherical resonator gyro damping non-uniformity identification method Download PDF

Info

Publication number
CN116576886B
CN116576886B CN202310502303.2A CN202310502303A CN116576886B CN 116576886 B CN116576886 B CN 116576886B CN 202310502303 A CN202310502303 A CN 202310502303A CN 116576886 B CN116576886 B CN 116576886B
Authority
CN
China
Prior art keywords
damping
harmonic oscillator
hemispherical
gyro
turntable
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
CN202310502303.2A
Other languages
Chinese (zh)
Other versions
CN116576886A (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 CN202310502303.2A priority Critical patent/CN116576886B/en
Publication of CN116576886A publication Critical patent/CN116576886A/en
Application granted granted Critical
Publication of CN116576886B publication Critical patent/CN116576886B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C25/00Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
    • G01C25/005Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass initial alignment, calibration or starting-up of inertial devices
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C19/00Gyroscopes; Turn-sensitive devices using vibrating masses; Turn-sensitive devices without moving masses; Measuring angular rate using gyroscopic effects
    • G01C19/56Turn-sensitive devices using vibrating masses, e.g. vibratory angular rate sensors based on Coriolis forces
    • G01C19/5776Signal processing not specific to any of the devices covered by groups G01C19/5607 - G01C19/5719
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

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

Abstract

A hemispherical resonator gyro damping non-uniformity identification method belongs to the technical field of inertia. The invention solves the problem of poor gyro output performance caused by harmonic oscillator standing wave drift. The main technical scheme adopted by the invention is as follows: step 1, fixing a hemispherical resonator gyro on a speed turntable, wherein a sensitive axis of the hemispherical resonator gyro is parallel to a rotary axis of the turntable; step 2, driving the hemispherical resonant gyroscope to start vibrating, and maintaining a given amplitude; step 3, driving the turntable to rotate; step 4, cutting off a hemispherical resonant gyro driving loop after the rotating speed of the turntable is stable; step 5, collecting hemispherical resonance gyro vibration signals by adopting an upper computer; step 6, resolving harmonic oscillator energy E; step 7, identifying azimuth angles of damping shafts and damping uneven amplitude values; and 8, compensating harmonic oscillator standing wave drift. The method can be applied to the identification of the damping non-uniformity of the hemispherical resonator gyroscope.

Description

Hemispherical resonator gyro damping non-uniformity identification method
Technical Field
The invention belongs to the technical field of inertia, and particularly relates to a hemispherical resonator gyro damping non-uniformity identification method.
Background
The hemispherical resonant gyroscope works in the principle that when an external angle input exists, the harmonic oscillator vibration standing wave produces precession under the action of coriolis force, and the precession angle is in direct proportion to the input angle. And the harmonic oscillator circumferential damping distribution is uneven due to the limitation of the harmonic oscillator processing technology, and the damping is distributed in a sine way along the harmonic oscillator circumferential direction. The harmonic oscillator standing wave can drift towards the damping shaft, so that the zero-bias stability of the hemispherical resonator gyroscope is affected.
In summary, since the damping distribution is uneven and the harmonic standing wave drift is caused, and the harmonic standing wave drift is caused to have poor gyro output performance, it is necessary to provide a method for reducing gyro drift.
Disclosure of Invention
The invention aims to solve the problem of poor gyro output performance caused by harmonic oscillator standing wave drift, and provides a hemispherical resonator gyro damping non-uniform identification method.
The technical scheme adopted by the invention for solving the technical problems is as follows:
a hemispherical resonator gyro damping non-uniformity identification method specifically comprises the following steps:
step 1, fixing a hemispherical resonator gyro on a speed turntable, wherein a sensitive axis of the hemispherical resonator gyro is parallel to a rotary axis of the turntable;
step 2, driving the hemispherical resonant gyroscope to start vibrating, and maintaining a given amplitude;
step 3, driving the turntable to rotate;
step 4, cutting off a hemispherical resonant gyro driving loop after the rotating speed of the turntable is stable;
step 5, collecting hemispherical resonant gyro vibration signals C by adopting an upper computer x 、S x 、C y And S is y
Step 6, according to the vibration signal C x 、S x 、C y And S is y Resolving harmonic oscillator energy E;
step 7, fitting the obtained harmonic oscillator energy E through a nonlinear least square method according to a harmonic oscillator free attenuation model to obtain an estimated value of the damping axis azimuth angleResistance and method for manufacturing sameEstimated value of Ni-inhomogeneity amplitude +.>
And 8, compensating the standing wave drift of the harmonic oscillator according to the estimated value of the azimuth angle of the damping shaft and the estimated value of the damping non-uniform amplitude.
Further, in the step 3, the turntable is driven to rotate at a speed of 500 °/s.
Further, the vibration signal C x 、S x 、C y And S is y The method comprises the following steps of:
wherein a represents the main wave amplitude of the harmonic oscillator, q represents the orthogonal wave amplitude, delta represents the phase difference between the reference signal and the actual vibration signal, and theta represents the azimuth angle of the standing wave of the harmonic oscillator.
Further, in the step 6, according to the vibration signal C x 、S x 、C y And S is y Resolving harmonic oscillator energy E; the method comprises the following steps:
E=C x 2 +S x 2 +C y 2 +S y 2
further, the harmonic oscillator free-decay model is:
wherein,is the first derivative of E, τ is the time constant, θ τ For damping axis azimuth, delta (τ -1 ) To damp non-uniform amplitude.
Further, an estimate of the damping axis azimuth angleEstimated value of damping non-uniformity amplitudeThe identification process of (2) is as follows:
step S1, setting an initial value theta of a damping axis azimuth angle τ (0) =0, the initial value Δ (τ -1 )(0)=0;
Step S2, calculating a value function of the current k moment:
wherein,for collecting the actual detection value of harmonic oscillator energy at the current k moment after signal processing, the method comprises the steps of +.>The theoretical value of the harmonic oscillator energy at the current k moment is calculated according to the harmonic oscillator free attenuation model;
s3, calculating a jacobian matrix J of the current k moment of the value function r (k):
S4, calculating the increment of the azimuth angle of the damping shaft and the uneven damping amplitude at the current k moment:
wherein, the upper corner mark T represents the transposition of the matrix, the upper corner mark-1 represents the inverse of the matrix, delta theta τ (k) Is the increment of the damping axis azimuth angle at the current k moment, delta [ delta (tau) -1 )](k) Is the increment of the damping non-uniform amplitude at the current k moment;
step S5, updating the azimuth angle of the damping shaft and the uneven damping amplitude at the next moment:
wherein θ τ (k+1) is the damping axis azimuth angle at time k+1, Δ (τ) -1 ) (k+1) is the damping non-uniform amplitude at time k+1;
and S6, judging whether data are input, if so, making k=k+1, and jumping to the step S2, otherwise, completing fitting.
The beneficial effects of the invention are as follows:
the invention provides a hemispherical resonance gyro damping non-uniformity identification method, which can simultaneously identify the azimuth angle of a damping shaft and the damping non-uniformity amplitude by adopting a harmonic oscillator free attenuation model, then compensate according to the azimuth angle of the damping shaft and the damping non-uniformity amplitude in a control loop, and can reduce harmonic oscillator standing wave drift, thereby improving gyro output performance.
Drawings
FIG. 1 is a flow chart of a method for identifying damping non-uniformity of a hemispherical resonator gyroscope.
Detailed Description
The present application is described in further detail below with reference to the accompanying drawings by way of specific embodiments. It will be apparent that the described embodiments are only some, but not all, of the embodiments of the invention. Based on the embodiments of the present invention, other embodiments obtained by those skilled in the art without making any inventive effort are within the scope of the present invention.
Detailed description of the inventionin the first embodiment, this embodiment will be described with reference to fig. 1. The method for identifying damping unevenness of hemispherical resonator gyroscopes specifically comprises the following steps:
step 1, fixing a hemispherical resonator gyro on a speed turntable, wherein a sensitive axis of the hemispherical resonator gyro is parallel to a rotary axis of the turntable;
step 2, driving the hemispherical resonant gyroscope to start vibrating; and maintaining a certain amplitude;
step 3, driving the turntable to rotate;
step 4, cutting off a hemispherical resonant gyro driving loop after the rotating speed of the turntable is stable;
step 5, collecting hemispherical resonant gyro vibration signals C by adopting an upper computer x 、S x 、C y And S is y
Step 6, according to the vibration signal C x 、S x 、C y And S is y Resolving harmonic oscillator energy E;
step 7, fitting the obtained harmonic oscillator energy E through a nonlinear least square method according to a harmonic oscillator free attenuation model to obtain an estimated value of the damping axis azimuth angleEstimated value of damping non-uniform amplitude +.>
The invention can also be realized by adopting other nonlinear identification algorithms such as extended Kalman filtering and the like.
And 8, compensating the standing wave drift of the harmonic oscillator according to the estimated value of the azimuth angle of the damping shaft and the estimated value of the damping non-uniform amplitude.
The second embodiment is as follows: in the present embodiment, the step 3 is different from the specific embodiment in that the turntable is driven to rotate at a speed of 500 °/s.
Other steps and parameters are the same as in the first embodiment.
And a third specific embodiment: this embodiment differs from the first or second embodiments in that the vibration signal C x 、S x 、C y And S is y The method comprises the following steps of:
wherein a represents the main wave amplitude of the harmonic oscillator, q represents the orthogonal wave amplitude, delta represents the phase difference between the reference signal and the actual vibration signal, and theta represents the azimuth angle of the standing wave of the harmonic oscillator.
Other steps and parameters are the same as in the first or second embodiment.
The specific embodiment IV is as follows: this embodiment differs from one to three embodiments in that in the step 6, the vibration signal C is used x 、S x 、C y And S is y Resolving harmonic oscillator energy E; the method comprises the following steps:
E=C x 2 +S x 2 +C y 2 +S y 2
other steps and parameters are the same as in one to three embodiments.
Fifth embodiment: the difference between this embodiment and one to four embodiments is that the harmonic oscillator free-decay model is:
wherein,is the first derivative of E, τ is the time constant, θ τ For damping axis azimuth, delta (τ -1 ) To damp non-uniform amplitude.
Other steps and parameters are the same as in one to four embodiments.
Specific embodiment six: this embodiment differs from one to fifth embodiments in that the estimated value of the damping axis azimuth angleEstimated value of damping non-uniform amplitude +.>The identification process of (2) is as follows:
step S1, setting an initial value theta of a damping axis azimuth angle τ (0) =0, the initial value Δ (τ -1 )(0)=0;
Step S2, calculating a value function of the current k moment:
wherein,for collecting the actual detection value of harmonic oscillator energy at the current k moment after signal processing, the method comprises the steps of +.>The theoretical value of the harmonic oscillator energy at the current k moment is calculated according to the harmonic oscillator free attenuation model;
s3, calculating a jacobian matrix J of the current k moment of the value function r (k):
S4, calculating the increment of the azimuth angle of the damping shaft and the uneven damping amplitude at the current k moment:
wherein, the upper corner mark T represents the transposition of the matrix, the upper corner mark-1 represents the inverse of the matrix, delta theta τ (k) Is the increment of the damping axis azimuth angle at the current k moment, delta [ delta (tau) -1 )](k) Is the increment of the damping non-uniform amplitude at the current k moment;
step S5, updating the azimuth angle of the damping shaft and the uneven damping amplitude at the next moment:
wherein θ τ (k+1) is the damping axis azimuth angle at time k+1, Δ (τ) -1 ) (k+1) is the damping non-uniform amplitude at time k+1;
and S6, judging whether data are input, if so, making k=k+1, and jumping to the step S2, otherwise, completing fitting.
Other steps and parameters are the same as in one of the first to fifth embodiments.
The above examples of the present invention are only for describing the calculation model and calculation flow of the present invention in detail, and are not limiting of the embodiments of the present invention. Other variations and modifications of the above description will be apparent to those of ordinary skill in the art, and it is not intended to be exhaustive of all embodiments, all of which are within the scope of the invention.

Claims (4)

1. A hemispherical resonance gyro damping non-uniformity identification method is characterized by comprising the following steps:
step 1, fixing a hemispherical resonator gyro on a speed turntable, wherein a sensitive axis of the hemispherical resonator gyro is parallel to a rotary axis of the turntable;
step 2, driving the hemispherical resonant gyroscope to start vibrating, and maintaining a given amplitude;
step 3, driving the turntable to rotate;
step 4, cutting off a hemispherical resonant gyro driving loop after the rotating speed of the turntable is stable;
step 5, collecting hemispherical resonant gyro vibration signals C by adopting an upper computer x 、S x 、C y And S is y
Step 6, according to the vibration signal C x 、S x 、C y And S is y Resolving harmonic oscillator energy E;
step 7, fitting the obtained harmonic oscillator energy E through a nonlinear least square method according to a harmonic oscillator free attenuation model to obtain dampingEstimated value of axial azimuthEstimated value of damping non-uniform amplitude +.>
The harmonic oscillator free attenuation model is as follows:
wherein,is the first derivative of E, τ is the time constant, θ τ For damping axis azimuth, delta (τ -1 ) For damping non-uniform amplitude, θ represents harmonic oscillator standing wave azimuth;
estimate of the damping axis azimuthEstimated value of damping non-uniform amplitude +.>The identification process of (2) is as follows:
step S1, setting an initial value theta of a damping axis azimuth angle τ (0) =0, the initial value Δ (τ -1 )(0)=0;
Step S2, calculating a value function of the current k moment:
wherein,for collecting the actual detection value of harmonic oscillator energy at the current k moment after signal processing, the method comprises the steps of +.>The theoretical value of the harmonic oscillator energy at the current k moment is calculated according to the harmonic oscillator free attenuation model;
s3, calculating a jacobian matrix J of the current k moment of the value function r (k):
S4, calculating the increment of the azimuth angle of the damping shaft and the uneven damping amplitude at the current k moment:
wherein, the upper corner mark T represents the transposition of the matrix, the upper corner mark-1 represents the inverse of the matrix, delta theta τ (k) Is the increment of the damping axis azimuth angle at the current k moment, delta [ delta (tau) -1 )](k) Is the increment of the damping non-uniform amplitude at the current k moment;
step S5, updating the azimuth angle of the damping shaft and the uneven damping amplitude at the next moment:
wherein θ τ (k+1) is the damping axis azimuth angle at time k+1, Δ (τ) -1 ) (k+1) is the damping non-uniform amplitude at time k+1;
step S6, judging whether data are input, if so, enabling k=k+1, and jumping to step S2, otherwise, completing fitting;
and 8, compensating the standing wave drift of the harmonic oscillator according to the estimated value of the azimuth angle of the damping shaft and the estimated value of the damping non-uniform amplitude.
2. The method for identifying damping unevenness of a hemispherical resonator gyro according to claim 1, wherein in the step 3, the turntable is driven to rotate at a speed of 500 °/s.
3. The method for identifying damping unevenness of hemispherical resonator gyroscope according to claim 2, wherein the vibration signal C is x 、S x 、C y And S is y The method comprises the following steps of:
wherein a represents the main wave amplitude of the harmonic oscillator, q represents the orthogonal wave amplitude, delta represents the phase difference between the reference signal and the actual vibration signal, and theta represents the azimuth angle of the standing wave of the harmonic oscillator.
4. The method for identifying damping unevenness of hemispherical resonator gyroscope according to claim 3, wherein in step 6, the damping unevenness is identified according to the vibration signal C x 、S x 、C y And S is y Resolving harmonic oscillator energy E; the method comprises the following steps:
E=C x 2 +S x 2 +C y 2 +S y 2
CN202310502303.2A 2023-05-06 2023-05-06 Hemispherical resonator gyro damping non-uniformity identification method Active CN116576886B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310502303.2A CN116576886B (en) 2023-05-06 2023-05-06 Hemispherical resonator gyro damping non-uniformity identification method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310502303.2A CN116576886B (en) 2023-05-06 2023-05-06 Hemispherical resonator gyro damping non-uniformity identification method

Publications (2)

Publication Number Publication Date
CN116576886A CN116576886A (en) 2023-08-11
CN116576886B true CN116576886B (en) 2024-01-26

Family

ID=87544685

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310502303.2A Active CN116576886B (en) 2023-05-06 2023-05-06 Hemispherical resonator gyro damping non-uniformity identification method

Country Status (1)

Country Link
CN (1) CN116576886B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116878477B (en) * 2023-09-06 2023-11-10 湖南二零八先进科技有限公司 Hemispherical resonator gyro damping non-uniformity parameter identification method, equipment and storage medium

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114608612A (en) * 2022-03-11 2022-06-10 中国船舶重工集团公司第七0七研究所 Online compensation system and method for damping non-uniform drift of full-angle mode resonant gyroscope
CN115540903A (en) * 2022-09-14 2022-12-30 哈尔滨工业大学 Method and system for correcting and adjusting non-uniform damping and rigidity of hemispherical resonator gyroscope
CN115950451A (en) * 2023-01-05 2023-04-11 中国船舶集团有限公司第七0七研究所 Damping uneven drift compensation method for full-angle hemispherical resonant gyroscope

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7437253B2 (en) * 2004-07-29 2008-10-14 The Boeing Company Parametrically disciplined operation of a vibratory gyroscope

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114608612A (en) * 2022-03-11 2022-06-10 中国船舶重工集团公司第七0七研究所 Online compensation system and method for damping non-uniform drift of full-angle mode resonant gyroscope
CN115540903A (en) * 2022-09-14 2022-12-30 哈尔滨工业大学 Method and system for correcting and adjusting non-uniform damping and rigidity of hemispherical resonator gyroscope
CN115950451A (en) * 2023-01-05 2023-04-11 中国船舶集团有限公司第七0七研究所 Damping uneven drift compensation method for full-angle hemispherical resonant gyroscope

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Research on the compensation method for damping mismatch of HRG;Fan Bai 等;《EEI 2022:4th International Conference on Electronic Engineering and Informatics》;全文 *
半球谐振陀螺控制及补偿技术;伊国兴 等;《宇航学报》;第41卷(第6期);全文 *
半球谐振陀螺阻尼不均匀误差补偿方法;郭杰 等;《飞控与探测》;第4卷(第1期);全文 *
基于主动驱动旋转的全角模式半球谐振陀螺驻波漂移补偿方法;王奇 等;《中国惯性技术学报》;第29卷(第4期);全文 *

Also Published As

Publication number Publication date
CN116576886A (en) 2023-08-11

Similar Documents

Publication Publication Date Title
CN116576886B (en) Hemispherical resonator gyro damping non-uniformity identification method
CN111536994B (en) Resonant micro gyroscope multi-mode cooperative control method and system and resonant micro gyroscope
CN113670340B (en) Method and system for measuring standing wave azimuth angle of hemispherical resonator gyroscope based on X/Y signal phase difference identification
CN115540903A (en) Method and system for correcting and adjusting non-uniform damping and rigidity of hemispherical resonator gyroscope
CN115876182B (en) Electrode error modeling method of hemispherical resonator gyroscope
CN110686662A (en) Dual-mode differential resonant gyroscope system capable of realizing online self-calibration
CN116086485A (en) Hemispherical resonator gyro error force compensation method and device
US20160202060A1 (en) High bandwidth coriolis vibratory gyroscope (cvg) with in-situ bias self-calibration
CN113899365A (en) Standing wave azimuth angle measurement method based on asymmetric parameter identification of hemispherical resonator gyroscope detection channel
CN113447047B (en) Rate integral gyroscope detection electrode error identification method, device, system and medium
CN114509057A (en) Full-angle mode control method of resonant gyroscope
CN114383590A (en) Phase error identification and compensation method of rate integral gyroscope
CN115077564B (en) Identification and compensation method for unbalanced error of hemispherical resonator gyro detection channel
CN114370886B (en) Full-angle mode vibration gyro measurement error self-calibration method based on virtual rotation
CN116448142A (en) Hemispherical resonator gyro balance mode calibration method
CN114440852A (en) Hemispherical resonator gyroscope X/Y path drive gain asymmetry error identification method
CN116772818A (en) Full angle mode gyroscope damping mismatch compensation method and system based on standing wave pseudo precession
CN115077561B (en) Method and system for adaptively compensating damping anisotropy of hemispherical harmonic oscillator
CN115950451A (en) Damping uneven drift compensation method for full-angle hemispherical resonant gyroscope
Xu et al. Single-channel control for hemispherical resonator gyro based on time division multiplexing and demultiplexing
CN106289209A (en) A kind of gyroscope control method being applicable to wide-range and control system
JP6732453B2 (en) Method for redistribution control of Coriolis vibrating gyroscope (cvg) for performance improvement
CN116576885A (en) Hemispherical resonator gyro starting method and hemispherical resonator gyro starting system based on temperature calibration and compensation
CN112947614B (en) Active vibration control method of variable speed tilting momentum wheel
CN111930154B (en) Nonsingular terminal sliding mode control method and device of MEMS gyroscope

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