CN109579878B - Inertial element error model rapid identification method based on frequency scanning excitation signal - Google Patents

Inertial element error model rapid identification method based on frequency scanning excitation signal Download PDF

Info

Publication number
CN109579878B
CN109579878B CN201910079721.9A CN201910079721A CN109579878B CN 109579878 B CN109579878 B CN 109579878B CN 201910079721 A CN201910079721 A CN 201910079721A CN 109579878 B CN109579878 B CN 109579878B
Authority
CN
China
Prior art keywords
excitation signal
frequency
error model
inertial
test
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
CN201910079721.9A
Other languages
Chinese (zh)
Other versions
CN109579878A (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 CN201910079721.9A priority Critical patent/CN109579878B/en
Publication of CN109579878A publication Critical patent/CN109579878A/en
Application granted granted Critical
Publication of CN109579878B publication Critical patent/CN109579878B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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

Landscapes

  • Engineering & Computer Science (AREA)
  • Manufacturing & Machinery (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
  • Gyroscopes (AREA)

Abstract

The invention provides a method for quickly identifying an inertial element error model based on a frequency scanning excitation signal, and belongs to the field of inertial test. The application starts from the application of an inertial element in aerospace tasks or weaponry and other aspects, and provides an excitation signal capable of fitting the working state of the inertial element by analyzing the working environment, wherein the excitation signal excites high-order terms which are often ignored in a parameter model compared with the existing signal, so that the aim of more accurately and truly obtaining the high-order term coefficient in an error model is fulfilled, and a cushion is laid for the subsequent compensation work; the time cost of the existing testing method is greatly reduced for the design of the whole testing method and the design of the error model parameter identification method, and the identification precision and the identification comprehensiveness of each order of parameters of the error model are greatly improved compared with the existing quick calibration method.

Description

Inertial element error model rapid identification method based on frequency scanning excitation signal
Technical Field
The invention relates to a method for quickly identifying an inertial element error model based on a frequency scanning excitation signal, and belongs to the field of inertial test.
Background
The inertia element test means that output data of the inertia element is obtained on given input special test equipment, and experimental data is filtered, fitted and the like to obtain performance parameters of the inertia element or parameter items of different orders are separated, so that an error model of the inertia element is perfected, and theoretical support is provided for subsequent compensation. Currently, the error model of the inertial element includes a static error model, a dynamic error model, and a random error model. The static error model reflects the relation between the error and the specific force, the dynamic error model reflects the relation between the error and the angular motion, the static error model and the dynamic error model can compensate the error and the angular motion to a large extent through mathematical processing, and the part which cannot be compensated belongs to the random error model along with experimental environment factors, excitation coupling and the like. Therefore, the random error is an important index representing the performance of the inertial element.
At present, the test items for the gyroscope mainly include a drift test, a static test and a dynamic test. The static test mainly comprises a force feedback method, polar axis rolling and the like, and the dynamic test mainly completes error model identification by processing response data of an inertial element to a multi-axis turntable excited single-frequency dynamic excitation test signal. The test of the accelerometer mainly comprises a gravity field rolling experiment, a centrifuge test, a vibration table test and the like. In many test methods, static test development is mature, and currently, the random error model is determined by using a model in a static state as a result, but the adaptability verification in dynamic test is difficult to complete. The dynamic test excitation signal is generally a uniform angular velocity or a coupling angular velocity that makes the uniform angular velocity in a specific direction, and only represents the response of the inertial element at a certain frequency in a frequency domain. In order to obtain error model parameters, a series of static and dynamic combined experiments need to be carried out in batches, and under the condition that the inertial elements need to be produced, tested and calibrated in a large scale, the consumption of manpower, financial resources, material resources and time cost is large; meanwhile, static tests or single-frequency dynamic tests performed on the inertia element cannot completely simulate the dynamic excitation borne by the inertia element in actual operation.
The fast self-calibration method for the error coefficient of the high-precision inertial gyroscope is provided by the Yan-hai swallow in research on the self-calibration method for the error coefficient of the high-precision inertial gyroscope, the calibration method is carried out by means of a platform provided with the gyroscope, external auxiliary equipment is not needed, the attitude of the platform is continuously changed relative to the earth under the action of excitation torque, the inertial navigation platform generates an attitude angle error of the platform under the driving of gyroscope drift, and therefore error information of the gyroscope can be separated. The specific test method is that the static gravity component excitation is applied by changing the position of the gyroscope (namely, the position of the gyroscope is turned over), so that the feedback output (the output of the platform attitude angle sensor) of the gyroscope at different positions is measured, and a simultaneous equation set is established, thereby solving each drift coefficient of the inertial instrument.
Wangjia proposes a rapid calibration method for an inertial measurement unit in an IMU440 inertial measurement unit rapid calibration method and experiment based on a biaxial rate turntable, and establishes a deterministic error model of a gyroscope and an accelerometer. The inertial measurement unit is integrated with three accelerometers and three gyroscopes, wherein the rate excitation applied to the gyroscopes is positive and negative bidirectional and the excitation applied to the accelerometers is static position excitation.
The general error model of the gyroscope is as follows: omegagRepresenting the steady state output of the gyroscope; kIRepresenting a gyroscope scale factor; omegaeIThe method is characterized in that the rotation angular velocity component of the earth on an input shaft is represented, and the input of the component in the direction of a sensitive shaft of a gyroscope can be 0 under the laboratory condition according to the change of the direction of a shaft system; dERepresenting zero offset of important index parameters of the gyroscope; dS、DI、DORespectively representing the sensitivity coefficients of the gyroscope rotating shaft, the input shaft and the output shaft; dIO、DOS、DSIRespectively representing the cross coupling coefficients of all the shafts under the coupling action of specific force components under the gravity field; dSS、DIIRespectively representing second-order sensitivity coefficients of the gyroscope in the directions of the rotation axis and the input axis; gI、gO、gSRespectively representing the components of specific force in the input shaft, the output shaft and the rotation shaft under the gravity field; epsilongRepresenting random errors in the test procedure.
ωg-K1ωeI=DF+DIgI+DOgO+DSgS+DIOgIgS+DOSgOgS
+DSIgSgI+DIIgi 2+DSSgs 2g
The existing error model parameter obtaining method needs to carry out a series of tests continuously, wherein the tests comprise a static gravity excitation method, a multi-position rolling method and a speed calibration method, each method can obtain partial parameters in a parameter model through corresponding calculation according to different side emphasis points and excitation mechanisms, the time for carrying out a gyroscope error model calibration test for one time is as long as dozens of hours, and a zero searching process is needed in each test. As described above, for the method widely used in static test, static test development is mature for testing inertial components, and currently, the random error model is determined by using the model in a static state as a result, but the adaptive verification of the static error model in dynamic test is difficult to perform. At present, the dynamic test excitation signal is generally a uniform angular velocity or a coupling angular velocity of the uniform angular velocity in a specific direction, and the excitation signal in the frequency domain mostly represents a peak value at a single frequency, so that the frequency domain response only represents the response of the inertia element at a certain frequency. In order to obtain error model parameters, a series of static and dynamic combined experiments need to be carried out in batches, and under the condition that the inertial elements need to be produced, tested and calibrated in a large scale, the consumption of manpower, financial resources, material resources and time cost is large; meanwhile, because the working environment of the inertia test element is a high-G environment or an environment with high overload and severe dynamic process, the static test or single-frequency dynamic test of the inertia element cannot truly restore the dynamic excitation borne by the inertia element in the actual work, and therefore, the judgment on the reliability of the obtained parameter model has deviation.
Disclosure of Invention
The invention aims to solve the problems in the prior art, and further provides a method for quickly identifying an inertial element error model based on a frequency scanning excitation signal.
The purpose of the invention is realized by the following technical scheme:
a method for quickly identifying an inertial element error model based on a frequency scanning excitation signal comprises the following specific steps: the method comprises the following steps: the excitation signal is designed to be expressed as
Figure BDA0001959992890000031
a=2π(fc-f0)/(n+1)Tn
b=2πf0
In the formula (f)cRepresenting a frequency of interest; f. of0Representing an initial frequency; t represents the excitation signal duration; n is0Represents the attenuation coefficient; g represents amplitude gain, t represents a time variable, and n is the order of the time variable;
step two: and (3) solving a first derivative and a second derivative of the excitation signal to obtain the velocity and the acceleration of the excitation signal:
real-time velocity information of the excitation signal:
Figure BDA0001959992890000033
acceleration information of real-time signals:
Figure BDA0001959992890000034
step three: by the formula
Figure BDA0001959992890000041
The method comprises the steps that an input signal obtains time of any frequency and obtains frequency of any time;
step four: and identifying parameters by adopting a harmonic analysis method and an entropy analysis method.
The invention relates to a method for quickly identifying an inertial element error model based on a frequency scanning excitation signal, wherein n is0、n、f0、fcAll can carry out the parameter value setting according to actual need equipment bearing range and components and parts operational environment promptly.
The invention relates to a method for quickly identifying an inertial element error model based on a frequency scanning excitation signal, which comprises the following steps of: and respectively obtaining the output of the inertia test element under each attitude condition by adopting a test mode of combining multiple attitudes, carrying out Fourier transform analysis on the output, obtaining a constant term coefficient, a first harmonic term coefficient and a second harmonic term coefficient by comparison, and identifying the error parameter model.
The invention relates to a method for quickly identifying an inertial element error model based on a frequency scanning excitation signal, which takes the application of an inertial element in the aspects of space missions or weaponry and the like as a starting point, provides an excitation signal capable of fitting the working state of the inertial element practically through the analysis of the working environment, and compared with the existing signal, the excitation signal excites the high-order terms which are often ignored in a parameter model, so that the purpose is to more accurately and truly obtain the high-order term coefficients in the error model and lay a cushion for the subsequent compensation work; the time cost of the existing testing method is greatly reduced for the design of the whole testing method and the design of the error model parameter identification method, and the identification precision and the identification comprehensiveness of each order of parameters of the error model are greatly improved compared with the existing quick calibration method.
Drawings
FIG. 1 is a waveform diagram of an excitation signal of the fast inertial element error model identification method based on a frequency-swept excitation signal according to the present invention.
FIG. 2 is a second order position signal plot of the excitation signal of the present invention.
FIG. 3 is a second order rate signal plot of the excitation signal of the present invention.
FIG. 4 is a second order acceleration signal curve of the excitation signal of the present invention.
Detailed Description
The invention will be described in further detail below with reference to the accompanying drawings: the present embodiment is implemented on the premise of the technical solution of the present invention, and a detailed implementation is given, but the scope of the present invention is not limited to the following embodiments.
The first embodiment is as follows: as shown in fig. 1 to 4, the method for quickly identifying an error model of an inertial component based on a frequency-scan excitation signal according to the present embodiment includes the following specific steps: the method comprises the following steps: the excitation signal is designed to be expressed as
Figure BDA0001959992890000051
a=2π(fc-f0)/(n+1)Tn
b=2πf0
In the formula (f)cRepresenting a frequency of interest; f. of0Representing an initial frequency; t represents the excitation signal duration; n is0Represents the attenuation coefficient; g represents amplitude gain, t represents a time variable, and n is the order of the time variable;
step two: and (3) solving a first derivative and a second derivative of the excitation signal to obtain the velocity and the acceleration of the excitation signal:
real-time velocity information of the excitation signal:
Figure BDA0001959992890000053
acceleration information of real-time signals:
Figure BDA0001959992890000054
step three: by passing
Figure BDA0001959992890000055
The method comprises the steps that an input signal obtains time of any frequency and obtains frequency of any time;
step four: and identifying parameters by adopting a harmonic analysis method and an entropy analysis method.
The invention relates to a method for quickly identifying an inertial element error model based on a frequency scanning excitation signal, n0、n、f0、fcAll can carry out the parameter value setting according to actual need equipment bearing range and components and parts operational environment promptly.
The input excitation adopted by the existing inertia test element test method is a gravity component, constant speed, multidirectional constant speed coupling or single frequency sinusoidal signal. The signals have different expression forms in the time domain, but are uniformly expressed as single spectrum content in the frequency domain, so that the full-coverage dynamic excitation of the inertia test element in the working environment frequency band of the actual task execution equipment cannot be met. To improve the existing testing methods, the present application is in exciting the signalThe design of the number firstly ensures that the signal has enough abundant frequency spectrum content; secondly, based on the analysis simulation of the working environment of the inertia test element, the amplitude and the frequency of the excitation signal are adjusted to enable the amplitude and the frequency to be within an allowable range so as to excite the inertia test element to the maximum extent. Considering the above conditions comprehensively, the excitation signal expression is designed as follows:
Figure BDA0001959992890000061
a=2π(fc-f0)/(n+1)Tn
b=2πf0
in the formula (f)cRepresenting a frequency of interest; f. of0Representing an initial frequency; t represents the excitation signal duration; n is0Represents the attenuation coefficient; g represents amplitude gain, t represents time variable, and n is the order of the time variable. n is0、n、f0、fcParameter value setting can be carried out according to actual needs (equipment bearing range and component working environment). Fig. 1 is a frequency domain characteristic of an excitation signal, and the amplitude of the excitation signal is continuous in a frequency domain of interest in a frequency domain characteristic diagram, which can satisfy excitation in a wide frequency domain.
Considering the performance upper limit of the inertia instrument test equipment in design, such as the limit of factors of bearable maximum overload, load and the like, the excitation waveform is introduced with n0And G to limit the magnitude margins of acceleration and velocity, respectively. The adjusting parameters enable the inertial instrument and the inertial instrument testing equipment not to be damaged when the input speed and the acceleration value are large by adjusting the waveform, and can effectively limit the acceleration and the speed position instruction within the testing capacity range of hardware. What the excitation signal represents is position information for which the first and second derivatives can be taken for its speed and acceleration, respectively:
Figure BDA0001959992890000063
Figure BDA0001959992890000064
y 'represents the real-time velocity information of the excitation signal, Y' represents the acceleration information of the real-time signal, and n is 2 in consideration of the performance of the inertial meter test equipment, and the waveform curves of the position, velocity and acceleration of the excitation signal are shown in fig. 2-4. The waveform transformation rate can analyze that the dynamic process starts more slowly with the increase of the value of the representation time order, but the dynamic process is accelerated quickly, so the time order can be determined according to the requirement. The frequency spectrum of the time domain continuous signal in the frequency domain is represented as a covering f0To fcA rectangular region of (a); the frequency band threshold value can be freely set according to the actual working characteristics of the inertial instrument. For a more complete test, the following expression was introduced:
Figure BDA0001959992890000071
Figure BDA0001959992890000072
the above formula enables acquisition of time of any frequency and acquisition of frequency of any time for an input signal, and provides convenience for actual test operation and data analysis. The practical significance of the individual parameters in the signal is explained below for a complete experimental design. In the testing process of the inertial instrument, the working frequency range of the inertial instrument can be obtained by knowing the working mode and the working performance of the inertial instrument, if the performance near a certain frequency is mined, the working frequency range can be solved according to the formula, the initial frequency can be randomly set in the effective working range, wherein the initial frequency is fc<f0And fc>f0Two cases. When f isc>f0In the test with the duration of T, the time for reaching any set frequency is TiThe working frequency of the inertia instrument at the T moment after the test is finished is the cut-off frequency fc. When f isc<f0In the case of tests of duration T, due to the presence of the attenuation factorThe second derivative will reach the peak value, i.e. the frequency is 0, and the time to reach the frequency 0 is tzIn this case, it is taken from 0 to tzAnd analyzing the test data at any moment.
Example two: as shown in fig. 1, in the present embodiment, a method for quickly identifying an error model of an inertial element based on a frequency-scan excitation signal includes: and respectively obtaining the output of the inertia test element under each attitude condition by adopting a test mode of combining multiple attitudes, carrying out Fourier transform analysis on the output, obtaining a constant term coefficient, a first harmonic term coefficient and a second harmonic term coefficient by comparison, and identifying the error parameter model.
The identification of the parameters of the error model of the inertial instrument refers to separating interested harmonic terms from output data of the inertial instrument by processing such as harmonic analysis, filtering and the like. The current main research method is to determine the system state variable by combining the characteristics of an inertial navigation test equipment error model and an inertial instrument error parameter model, and establish a specific expression of a system state equation according to the input compensated for an error source. The parameter identification mainly considers a harmonic analysis method and an entropy analysis method.
Harmonic analysis method:
the dynamic error model of a classical rotor gyroscope is shown as follows:
Figure BDA0001959992890000081
the accelerometer error model is shown as follows:
U=KF+KIaI+KIIaI 2+Koq|aI|aI+KIIIaI 3+KPaP+KPPaP 2+KPPPap 3+KOaO+KOOaO 2
+KOOOaO 3+KIPaIaP+KIOaIaO+KPOaPaO+KIPPaIaP 2+KIOOaIaO 2+KOPPaOaP 2
+KOIIaOaI 2+KPIIaPaI 2+KPOOaPaO 2+KIPOaIaPaO
the amplitude and the influence degree of the constant term and the primary term in the error model are larger than those of the high-order term, so that the influence of the second-order term and the high-order term in the traditional ship-based inertia test element can be ignored, but high overload is common in mission environments such as aerospace, strategic missiles and the like, so that the excitation degree of the high-order term is higher, and the identification of the second-order term is also of great significance in the identification of the system error model. In the experimental process, the output of the inertia test element under each attitude condition is respectively obtained by adopting a test mode of combining multiple attitudes, Fourier transform analysis is carried out on the output, a constant term coefficient, a first harmonic term coefficient and a second harmonic term coefficient can be obtained by comparison, and the identification of an error parameter model is realized.
Entropy analysis method:
entropy analysis methods and methods of estimation or model fitting based on isentropic optimization criteria such as maximum entropy and minimum mutual entropy have been widely used in various mathematical statistical analyses, information and signal processing, system analysis and identification, error analysis and data processing. The estimation principle of the maximum entropy method is a method which does not make any subjective assumption except the actual measurement data obtained, namely, the estimation is carried out according to the principle of the least certainty; the minimum mutual entropy method is to apply the existing data and also require the data to be in accordance with the constraint of the prior condition, and is closest to the prior information. In the application, the entropy analysis method is considered preferentially for data processing and error model coefficient separation, and more information related to model parameters is obtained from the output data of the inertial instrument from the information aspect.
Example three: as shown in fig. 1-4, the present embodiment relates to a method for fast identifying an error model of an inertial element based on a frequency-scan excitation signal, and the method has the following advantages:
compared with the prior art, the excitation principle designed by the application is greatly innovative in excitation mode, is not limited to excitation under static conditions (for example, gravity components under different postures are used as input or uniform motion is used as input), and is greatly innovative in input signal form and content when the user looks at dynamic excitation under high overload.
Compared with the prior art, the excitation signal designed by the application is more comprehensively improved on the excitation mechanism of the system parameters, and usually, the static error model coefficient of the gyroscope is calibrated through the turntable test under gravity. However, the input acceleration in the gravitational field is small and the second order error coefficients of the gyroscope cannot be excited. In practice, gyroscopes typically operate under high acceleration conditions, in which case high-precision navigation cannot ignore the effect of second-order errors on the accuracy of the gyroscope. Therefore, it is very important to calibrate the second order error coefficients of the gyroscope through a high acceleration device. In order to more truly obtain the dynamic characteristics of the system and more vividly restore the working environment of the inertia test element, a frequency scanning signal with high overload and wider frequency domain is selected on the excitation of the system, and high-order parameters in the system are more comprehensively excited, so that the reliability and the authenticity of the identification of high-order term coefficients in an error model of the inertia test element are greatly improved.
Compared with the prior art, the testing method designed by the application has the advantages that the testing operation redundancy, the initial state error introduction and the testing time cost are greatly improved. In the prior art, error model parameters of an inertia test element are often determined through a series of tests combining a multi-position method, a polar axis rolling method and a geographic coordinate axis rolling method under a static condition, the test time can be dozens of hours, and the test method designed by the application can be realized by repeatedly exciting for dozens of seconds. Secondly, to all can have error factors such as measurement noise, alignment deviation and installation error to experimental every time, have certain influence to error model parameter identification precision, and this application can reach the purpose through repeated short-time excitation signal, has realized the purpose of high-efficient test, quick test.
Compared with the prior art, the output data processing mode applied by the method has the advantages that harmonic analysis is carried out on the system model parameter acquisition through the result of Fourier analysis to determine each order parameter of the system at one time or least square fitting is carried out on each attitude data to obtain each order parameter, and the calculation complexity is greatly optimized compared with the traditional method.
The above description is only a preferred embodiment of the present invention, and these embodiments are based on different implementations of the present invention, and the scope of the present invention is not limited thereto, and any changes or substitutions that can be easily conceived by those skilled in the art within the technical scope of the present invention are included in the scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.

Claims (3)

1. The method for quickly identifying the inertial element error model based on the frequency scanning excitation signal is characterized by comprising the following specific steps of:
the method comprises the following steps: the excitation signal is designed to be expressed as
Figure FDA0003082489980000011
Figure FDA0003082489980000012
b=2πf0
In the formula (f)cRepresenting a frequency of interest; f. of0Representing an initial frequency; t represents the excitation signal duration; n is0Represents the attenuation coefficient; g represents amplitude gain, t represents a time variable, and n is the order of the time variable;
step two: and (3) solving a first derivative and a second derivative of the excitation signal to obtain the velocity and the acceleration of the excitation signal:
real-time velocity information of the excitation signal:
Figure FDA0003082489980000013
acceleration information of real-time signals:
Figure FDA0003082489980000014
step three: by the formula
Figure FDA0003082489980000015
The method comprises the steps that an input signal obtains time of any frequency and obtains frequency of any time;
step four: and identifying parameters by adopting a harmonic analysis method and an entropy analysis method.
2. The method for fast identification of an inertial element error model based on a swept-frequency excitation signal according to claim 1, wherein n is0、n、f0、fcAll carry out the parameter value setting according to actual need equipment bearing range and components and parts operational environment promptly.
3. The method for rapidly identifying an error model of an inertial element based on a frequency-swept excitation signal according to claim 1, wherein the harmonic analysis method specifically comprises: and respectively obtaining the output of the inertia test element under each attitude condition by adopting a test mode of combining multiple attitudes, carrying out Fourier transform analysis on the output, obtaining a constant term coefficient, a first harmonic term coefficient and a second harmonic term coefficient by comparison, and identifying the error parameter model.
CN201910079721.9A 2019-01-28 2019-01-28 Inertial element error model rapid identification method based on frequency scanning excitation signal Active CN109579878B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910079721.9A CN109579878B (en) 2019-01-28 2019-01-28 Inertial element error model rapid identification method based on frequency scanning excitation signal

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910079721.9A CN109579878B (en) 2019-01-28 2019-01-28 Inertial element error model rapid identification method based on frequency scanning excitation signal

Publications (2)

Publication Number Publication Date
CN109579878A CN109579878A (en) 2019-04-05
CN109579878B true CN109579878B (en) 2021-07-02

Family

ID=65918210

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910079721.9A Active CN109579878B (en) 2019-01-28 2019-01-28 Inertial element error model rapid identification method based on frequency scanning excitation signal

Country Status (1)

Country Link
CN (1) CN109579878B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114324977B (en) * 2021-12-03 2022-09-02 伸瑞科技(北京)有限公司 Method for calibrating accelerometer of inertial navigation on centrifugal machine based on full-error analysis

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103558843A (en) * 2013-11-05 2014-02-05 中国航空工业集团公司西安飞机设计研究所 Automatic amplitude modulation frequency scanning method for airplane servo elasticity frequency response test
CN104215262A (en) * 2014-08-29 2014-12-17 南京航空航天大学 On-line dynamic inertia sensor error identification method of inertia navigation system
CN106844880A (en) * 2016-12-29 2017-06-13 中国空间技术研究院 A kind of satellite momentum wheel disturbs shake experiment and data interpretation method
CN107063295A (en) * 2016-10-26 2017-08-18 中国矿业大学(北京) A kind of method for analyzing stability of resonant mode gyro
CN107764257A (en) * 2017-09-14 2018-03-06 中国电子科技集团公司第五十四研究所 A kind of inertia device method for numerical simulation

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10071824B2 (en) * 2012-02-08 2018-09-11 Alan D Reth Method and apparatus for spacecraft gyroscope scale factor calibration

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103558843A (en) * 2013-11-05 2014-02-05 中国航空工业集团公司西安飞机设计研究所 Automatic amplitude modulation frequency scanning method for airplane servo elasticity frequency response test
CN104215262A (en) * 2014-08-29 2014-12-17 南京航空航天大学 On-line dynamic inertia sensor error identification method of inertia navigation system
CN107063295A (en) * 2016-10-26 2017-08-18 中国矿业大学(北京) A kind of method for analyzing stability of resonant mode gyro
CN106844880A (en) * 2016-12-29 2017-06-13 中国空间技术研究院 A kind of satellite momentum wheel disturbs shake experiment and data interpretation method
CN107764257A (en) * 2017-09-14 2018-03-06 中国电子科技集团公司第五十四研究所 A kind of inertia device method for numerical simulation

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于参数激励的MEMS环式陀螺驱动方法与实现;钱磊等;《计算机测量与控制》;20181231;第26卷(第11期);全文 *

Also Published As

Publication number Publication date
CN109579878A (en) 2019-04-05

Similar Documents

Publication Publication Date Title
CN101290326B (en) Parameter identification calibration method for rock quartz flexibility accelerometer measuring component
Zhang et al. Coupled microvibration analysis of a reaction wheel assembly including gyroscopic effects in its accelerance
Chen et al. On-chip capacitive sensing and tilting motion estimation of a micro-stage for in situ MEMS gyroscope calibration
Wei et al. High-precision synchronous test method of vibration performance parameters for fused quartz hemispherical resonator
EP1976760B1 (en) System for and method of monitoring free play of aircraft control surfaces
CN103323625A (en) Error calibration compensation method of accelerometers in MEMS-IMU under dynamic environment
KR101658473B1 (en) Compensation method of acceleration sensitivity of mems gyroscope
Patil et al. A multi-view digital image correlation for extracting mode shapes of a tire
CN109579878B (en) Inertial element error model rapid identification method based on frequency scanning excitation signal
CN113865583B (en) Accelerometer combination dynamic installation deviation matrix determining and compensating method
Yingbo et al. Calibration method of quartz accelerometer on dynamic centrifuge
Ma et al. Analysis and compensation of mass imperfection effects on 3-D sensitive structure of bell-shaped vibratory gyro
Sun et al. Sequential calibration method of nonlinear errors of PIGA on counter-rotating platform centrifuge
CN104655876A (en) Method for calibrating linear accelerometer under situations of constant acceleration and vibration composite input
Zhou et al. A novel ring vibrating gyroscope based on side piezo-electrodes
CN110988400A (en) MEMS accelerometer combination calibration method and calibration device
Ren et al. A new method for calibrating nonlinear coefficients of PIGA on linear vibrator
CN113945230B (en) Identification method for high-order error coefficient of inertial device
Xing et al. Offline calibration for MEMS gyroscope G-sensitivity error coefficients based on the newton iteration and least square methods
Yang et al. The static unbalance analysis and its measurement system for gimbals axes of an inertial stabilization platform
Fakhari et al. A theoretical and experimental disturbance analysis in a product of inertia measurement system
Ermakov et al. Development of a vibrational error model of a hemispherical resonator gyroscope
Šiaudinytė et al. Modal analysis and experimental research into improved centering–leveling devices
Tang et al. An improved scale factor calibration model of MEMS gyroscopes
Pi et al. Implementation of effective scale factor matching and axial centripetal error compensation for a rotating accelerometer gravity gradiometer

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