CN101788305A - Method for rapid field calibration of micro inertial measurement unit - Google Patents
Method for rapid field calibration of micro inertial measurement unit Download PDFInfo
- Publication number
- CN101788305A CN101788305A CN201010134734A CN201010134734A CN101788305A CN 101788305 A CN101788305 A CN 101788305A CN 201010134734 A CN201010134734 A CN 201010134734A CN 201010134734 A CN201010134734 A CN 201010134734A CN 101788305 A CN101788305 A CN 101788305A
- Authority
- CN
- China
- Prior art keywords
- measurement unit
- inertial measurement
- micro inertial
- accelerometer
- micro
- 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.)
- Pending
Links
Landscapes
- Navigation (AREA)
Abstract
The invention relates to a method for calibration of a micro inertial measurement unit, in particular to a method for rapid field calibration of a micro inertial measurement unit, which solves the problems that the operation of the existing method for calibration of the micro inertial measurement unit is complex and time-consuming and the field calibration lacks direction and position criteria. The method for rapid field calibration of the micro inertial measurement unit is realized through the following steps: (1) establishing a calibration model; (2) accurately measuring the local gravitational acceleration of the fixed position point to be calibrated; (3) randomly rotating the micro inertial measurement unit in the fixed position point to be calibrated; (4) fitting an ellipsoid equation; (5) collating the fitted ellipsoid equation into a standard ellipsoid equation; (6) calculating the axial scale factor and null bias of an accelerometer for the micro inertial measurement unit; and (7) compensating the output information of the micro inertial measurement unit. The invention is applicable to the field calibration of the micro inertial measurement unit.
Description
Technical field
The present invention relates to the scaling method of Micro Inertial Measurement Unit, specifically is a kind of method for quick field calibration of micro inertial measurement unit.
Background technology
Micro Inertial Measurement Unit is the core component of micro-inertial measuring system.Before using Micro Inertial Measurement Unit, must carry out test calibration to it, with determine wherein inertia device (accelerometer) zero partially, demarcate parameters such as factor and between centers alignment error angle, so that the output information by Micro Inertial Measurement Unit accurately calculates carrier coordinate system each axial linear acceleration and angular velocity of rotation, and then provide input information really and accurately for the attitude algorithm and the navigation calculating of system.Correlation test shows, as time passes and environmental change, the between centers alignment error angle of the inertia device in the Micro Inertial Measurement Unit can not change, but inertia device zero partially, demarcate factor and can change.Wherein, zero of inertia device changes particularly evident partially, all can be different when working on power at every turn, and As time goes on, zero changes constantly aggravation partially.As everyone knows, zero of the inertia device error effect that changes partially navigation and location increases by secondary and the cube of time respectively.This shows, if the parameter of demarcating in Micro Inertial Measurement Unit testing laboratory is applied in the actual test in later stage, can produce than mistake, reduce the measuring accuracy of Micro Inertial Measurement Unit, and Micro Inertial Measurement Unit, particularly be applied to the Micro Inertial Measurement Unit of conventional weapons such as rocket projectile, shell equipment, need have higher measuring accuracy in a short time, is badly in need of a kind of can be before Micro Inertial Measurement Unit be used the scene method of demarcation quick and precisely.The scaling method of existing Micro Inertial Measurement Unit all needs microtest equipment that directional reference and position reference are provided, leave these accurate experimental facilitiess and carry out accurately demarcation of scene, can only carry out simple on-site and demarcate with regard to directly causing Micro Inertial Measurement Unit before reality is used, to be difficult to have ready conditions.And simple calibrating often lack direction benchmark and position reference, demarcating partially with zero of accelerometer is example, only be seek two fully the positions of symmetry just be difficult to accomplish.In a word, the scaling method of existing Micro Inertial Measurement Unit need rely on laboratory equipment, and complicated operation, time-consuming and stated accuracy is low is difficult to satisfy the application demand of micro-inertial measuring system.Based on this, be necessary to invent a kind of simple to operate, save time and method for quick field calibration of micro inertial measurement unit that stated accuracy is high.
Summary of the invention
The present invention is for the scaling method complicated operation that solves existing Micro Inertial Measurement Unit, time-consuming and on-site proving is lacked direction and position reference and the low problem of simple calibrating precision provides a kind of method for quick field calibration of micro inertial measurement unit.
The present invention adopts following technical scheme to realize: a kind of method for quick field calibration of micro inertial measurement unit, this method are to adopt following steps to realize:
(1) sets up peg model; The concrete form of peg model is as follows:
Wherein, k
Ij(i=x, y, z; J=x, y, z; I=j) be the axial constant multiplier of accelerometer,
(i=x, y, z; J=x, y, z; I ≠ j) is an accelerometer between centers coupling constant multiplier; u
I0(i=x, y, z) be x axle, y axle, z direction of principal axis accelerometer zero partially;
(i=x, y, z; J=x, y, z; The value of i ≠ j) is demarcated by Micro Inertial Measurement Unit testing laboratory;
(2) accurately measure local gravitational acceleration in fixed position point to be calibrated;
(3) rotate Micro Inertial Measurement Unit at random in fixed position point to be calibrated, make its attitude angle (angle between the coordinate axis of Micro Inertial Measurement Unit place coordinate system and the coordinate axis of quiet ground coordinate system) span cover Micro Inertial Measurement Unit place three dimensions scope (promptly making the coordinate axis end points motion of Micro Inertial Measurement Unit place coordinate system form the ellipsoid track) as far as possible, and then obtain a series of accelerometer output voltage measured values;
(4) simulate the ellipsoid equation by a series of accelerometer output voltage measured values;
(5) the ellipsoid equation arrangement that simulates is the standard spheroid equation, obtains spheroid-like matrix and ellipsoid centre coordinate;
(6) by spheroid-like matrix and ellipsoid centre coordinate, the axial constant multiplier of accelerometer that calculates Micro Inertial Measurement Unit is inclined to one side with zero;
(7), the output information of Micro Inertial Measurement Unit is compensated with the axial constant multiplier of the accelerometer that calculates and zero substitution peg model partially.
Usually, the general peg model of Micro Inertial Measurement Unit is:
In the formula (1): u
i(i=x, y z) are x axle, the y axle of Micro Inertial Measurement Unit, the output voltage measured value of z direction of principal axis accelerometer; k
Ij(i=x, y, z; J=x, y, z; I ≠ j) is an accelerometer between centers coupling constant multiplier; k
Ij(i=x, y, z; J=x, y, z; I=j) be the axial constant multiplier of accelerometer; a
i(i=x, y z) are the output accekeration of x axle, y axle, z direction of principal axis accelerometer; u
I0(i=x, y, z) be x axle, y axle, z direction of principal axis accelerometer zero partially.
Consider accelerometer between centers coupling constant multiplier k
Ij(i=x, y, z; J=x, y, z; I ≠ j) itself is very little, and variable quantity is very little in a short time, and consequent influence is very little, so uses recently accurately calibration result during on-site proving
The approximate k that replaces
Ij(i=x, y, z; J=x, y, z; I ≠ j) (promptly
Value demarcate by Micro Inertial Measurement Unit testing laboratory), the only bigger axial constant multiplier k of accelerometer during on-site proving to influence
Ij(i=x, y, z; J=x, y, z; I ≠ j) and zero is u partially
I0(i=x, y z) demarcate, and therefore, the peg model in the described step (1) specifically is expressed as follows for the peg model after approximate, the peg model after being similar to:
Since the vector of the output accekeration of x axle, y axle, z direction of principal axis accelerometer and be always local gravitational acceleration, that is:
a
Ta=[k
-1(u-u
0)]
T[k
-1(u-u
0)]=G
2 (3)
In the formula (3): G is a local gravitational acceleration;
Can obtain following second normal form equation by formula (3):
In the formula (3):
By formula (4) as can be known, the output voltage measured value of x axle, y axle, z direction of principal axis accelerometer satisfies a secondary Ellipsoidal Surface equation, and its geometric meaning is that the coordinate points of the output voltage measured value of x axle, y axle, z direction of principal axis accelerometer is positioned in measurement coordinate system on the Ellipsoidal Surface of being determined by formula (4).
In the described step (4), the process of match ellipsoid equation is based on the least squares theory of ellipsoid constraint, and fit procedure relies on computing machine to finish.
In the described step (5), the standard spheroid equation that the ellipsoid equation arrangement through simulating draws is as follows:
In the formula (5): S is the spheroid-like matrix; X
0Center point coordinate for ellipsoid.
Formula (5) can draw with formula (4) comparison:
Centre coordinate X according to spheroid-like matrix S and ellipsoid
0, can calculate the axial constant multiplier k of x axle, y axle, z direction of principal axis accelerometer
Ij(i=x, y, z; J=x, y, z; I ≠ j) and zero is u partially
I0(i=x, y, z).Concrete computation process is as follows:
The group of solving an equation calculates the axial constant multiplier k of accelerometer of Micro Inertial Measurement Unit
Ij(i=x, y, z; J=x, y, z; I ≠ j) and zero is u partially
I0(i=x, y, z) as follows:
With the axial constant multiplier k of the accelerometer that calculates
Ij(i=x, y, z; J=x, y, z; I ≠ j) and zero is u partially
I0(z) peg model in the described step of substitution (1) is in conjunction with the accelerometer between centers coupling constant multiplier of accurately being demarcated by Micro Inertial Measurement Unit testing laboratory in the recent period for i=x, y
(i=x, y, z; J=x, y z), promptly finishes the quick field calibration to Micro Inertial Measurement Unit, thereby makes Micro Inertial Measurement Unit obtain measurement result more accurately.Compare with the scaling method of existing Micro Inertial Measurement Unit, that method for quick field calibration of micro inertial measurement unit of the present invention has is easy to operate, save time, need not characteristics such as benchmark location, stated accuracy height, can satisfy the demand of micro-inertial navigation system practical application.
The present invention is according to the Micro Inertial Measurement Unit peg model of setting up, least square fitting theory in conjunction with the ellipsoid constraint, efficiently solve the scaling method complicated operation, time-consuming and on-site proving is lacked direction and position reference and the low problem of simple calibrating precision is widely used in the on-site proving of Micro Inertial Measurement Unit in each field of existing Micro Inertial Measurement Unit.
Embodiment
A kind of method for quick field calibration of micro inertial measurement unit, this method are to adopt following steps to realize:
(1) sets up peg model; The concrete form of peg model is as follows:
Wherein, k
Ij(i=x, y, z; J=x, y, z; I=j) be the axial constant multiplier of accelerometer,
(i=x, y, z; J=x, y, z; I ≠ j) is an accelerometer between centers coupling constant multiplier; u
I0(i=x, y, z) be x axle, y axle, z direction of principal axis accelerometer zero partially;
(i=x, y, z; J=x, y, z; The value of i ≠ j) is demarcated by Micro Inertial Measurement Unit testing laboratory;
(2) accurately measure local gravitational acceleration in fixed position point to be calibrated;
(3) rotate Micro Inertial Measurement Unit at random in fixed position point to be calibrated, make its attitude angle span cover Micro Inertial Measurement Unit place three dimensions scope as far as possible, and then obtain a series of accelerometer output voltage measured values;
(4) simulate the ellipsoid equation by a series of accelerometer output voltage measured values;
(5) the ellipsoid equation arrangement that simulates is the standard spheroid equation, obtains spheroid-like matrix and ellipsoid centre coordinate;
(6) by spheroid-like matrix and ellipsoid centre coordinate, the axial constant multiplier of accelerometer that calculates Micro Inertial Measurement Unit is inclined to one side with zero;
(7), the output information of Micro Inertial Measurement Unit is compensated with the axial constant multiplier of the accelerometer that calculates and zero substitution peg model partially.
Claims (1)
1. method for quick field calibration of micro inertial measurement unit is characterized in that: this method is to adopt following steps to realize:
(1) sets up peg model; The concrete form of peg model is as follows:
Wherein, k
Ij(i=x, y, z; J=x, y, z; I=j) be the axial constant multiplier of accelerometer,
(i=x, y, z; J=x, y, z; I ≠ j) is an accelerometer between centers coupling constant multiplier; u
I0(i=x, y, z) be x axle, y axle, z direction of principal axis accelerometer zero partially;
(i=x, y, z; J=x, y, z; The value of i ≠ j) is demarcated by Micro Inertial Measurement Unit testing laboratory;
(2) accurately measure local gravitational acceleration in fixed position point to be calibrated;
(3) rotate Micro Inertial Measurement Unit at random in fixed position point to be calibrated, make its attitude angle span cover Micro Inertial Measurement Unit place three dimensions scope as far as possible, and then obtain a series of accelerometer output voltage measured values;
(4) simulate the ellipsoid equation by a series of accelerometer output voltage measured values;
(5) the ellipsoid equation arrangement that simulates is the standard spheroid equation, obtains spheroid-like matrix and ellipsoid centre coordinate;
(6) by spheroid-like matrix and ellipsoid centre coordinate, the axial constant multiplier of accelerometer that calculates Micro Inertial Measurement Unit is inclined to one side with zero;
(7), the output information of Micro Inertial Measurement Unit is compensated with the axial constant multiplier of the accelerometer that calculates and zero substitution peg model partially.
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201010134734A CN101788305A (en) | 2010-03-26 | 2010-03-26 | Method for rapid field calibration of micro inertial measurement unit |
CN2010102004189A CN101832782B (en) | 2010-03-26 | 2010-06-01 | Method for quick field calibration of micro inertial measurement unit |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201010134734A CN101788305A (en) | 2010-03-26 | 2010-03-26 | Method for rapid field calibration of micro inertial measurement unit |
Publications (1)
Publication Number | Publication Date |
---|---|
CN101788305A true CN101788305A (en) | 2010-07-28 |
Family
ID=42531611
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201010134734A Pending CN101788305A (en) | 2010-03-26 | 2010-03-26 | Method for rapid field calibration of micro inertial measurement unit |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101788305A (en) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101832782A (en) * | 2010-03-26 | 2010-09-15 | 中北大学 | Method for quick field calibration of micro inertial measurement unit |
CN103782180A (en) * | 2011-09-08 | 2014-05-07 | 莱卡地球系统公开股份有限公司 | Method and system of recalibrating an inertial sensor |
CN105242070A (en) * | 2014-07-09 | 2016-01-13 | 北京自动化控制设备研究所 | Accelerometer unit calibration method without vector standard |
CN105424040A (en) * | 2016-01-15 | 2016-03-23 | 极翼机器人(上海)有限公司 | Novel MEMS (micro-electromechanical systems) inertial sensor array redundancy configuration method |
CN105606846A (en) * | 2015-09-18 | 2016-05-25 | 北京理工大学 | Accelerometer calibration method based on attitude information |
CN107044860A (en) * | 2017-04-17 | 2017-08-15 | 吉林化工学院 | A kind of measurement scaling method of inertial sensor |
CN110988400A (en) * | 2019-12-11 | 2020-04-10 | 北京自动化控制设备研究所 | MEMS accelerometer combination calibration method and calibration device |
CN113984090A (en) * | 2021-10-25 | 2022-01-28 | 北京科技大学 | Online calibration and compensation method and device for IMU (inertial measurement Unit) error of wheeled robot |
-
2010
- 2010-03-26 CN CN201010134734A patent/CN101788305A/en active Pending
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101832782B (en) * | 2010-03-26 | 2011-11-09 | 中北大学 | Method for quick field calibration of micro inertial measurement unit |
CN101832782A (en) * | 2010-03-26 | 2010-09-15 | 中北大学 | Method for quick field calibration of micro inertial measurement unit |
CN103782180A (en) * | 2011-09-08 | 2014-05-07 | 莱卡地球系统公开股份有限公司 | Method and system of recalibrating an inertial sensor |
CN105242070A (en) * | 2014-07-09 | 2016-01-13 | 北京自动化控制设备研究所 | Accelerometer unit calibration method without vector standard |
CN105242070B (en) * | 2014-07-09 | 2018-09-18 | 北京自动化控制设备研究所 | A kind of accelerometer combination scaling method of scalar potential standard |
CN105606846B (en) * | 2015-09-18 | 2018-06-19 | 北京理工大学 | A kind of accelerometer calibration method based on attitude information |
CN105606846A (en) * | 2015-09-18 | 2016-05-25 | 北京理工大学 | Accelerometer calibration method based on attitude information |
CN105424040B (en) * | 2016-01-15 | 2019-09-13 | 极翼机器人(上海)有限公司 | A kind of novel MEMS inertial sensor array redundant configurations method |
CN105424040A (en) * | 2016-01-15 | 2016-03-23 | 极翼机器人(上海)有限公司 | Novel MEMS (micro-electromechanical systems) inertial sensor array redundancy configuration method |
CN107044860A (en) * | 2017-04-17 | 2017-08-15 | 吉林化工学院 | A kind of measurement scaling method of inertial sensor |
CN110988400A (en) * | 2019-12-11 | 2020-04-10 | 北京自动化控制设备研究所 | MEMS accelerometer combination calibration method and calibration device |
CN113984090A (en) * | 2021-10-25 | 2022-01-28 | 北京科技大学 | Online calibration and compensation method and device for IMU (inertial measurement Unit) error of wheeled robot |
CN113984090B (en) * | 2021-10-25 | 2023-07-04 | 北京科技大学 | Wheel type robot IMU error online calibration and compensation method and device |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101788305A (en) | Method for rapid field calibration of micro inertial measurement unit | |
JP6613236B2 (en) | High-precision optical fiber gyro inertial surveying method | |
CN101290326B (en) | Parameter identification calibration method for rock quartz flexibility accelerometer measuring component | |
CN101839719A (en) | Inertial measurement unit based on gyroscope and geomagnetic sensor | |
CN101706284B (en) | Method for increasing position precision of optical fiber gyro strap-down inertial navigation system used by ship | |
CN104344836B (en) | Posture observation-based redundant inertial navigation system fiber-optic gyroscope system level calibration method | |
CN103852085B (en) | A kind of fiber strapdown inertial navigation system system for field scaling method based on least square fitting | |
CN103323625B (en) | Error calibration compensation method of accelerometers in MEMS-IMU under dynamic environment | |
CN105180968A (en) | IMU/magnetometer installation misalignment angle online filter calibration method | |
CN101246023A (en) | Closed-loop calibration method of micro-mechanical gyroscope inertial measuring component | |
CN101187568A (en) | Multi-position strapping north-seeking system direction effect calibration method | |
CN102937450B (en) | A kind of relative attitude defining method based on gyro to measure information | |
CN103389114B (en) | 3 axis MEMS gyro non-orthogonal errors scaling method based on g sensitivity | |
CN104316079B (en) | Drop point precision estimation method for inertia measurement system based on rocket sled test | |
CN103743413A (en) | Installation error online estimation and north-seeking error compensation method for modulating north seeker under inclined state | |
CN103983276A (en) | Three-frame four-axis inertial platform error calibration method based on navigation datum system | |
CN101738203A (en) | Optimal position calibration method of static drifting zero and primary acceleration related term error model of flexible gyroscope | |
CN110361031A (en) | A kind of IMU population parameter error quick calibrating method theoretical based on backtracking | |
CN106017452A (en) | Dual gyro anti-disturbance north-seeking method | |
CN103808349A (en) | Error correction method and device for vector sensors | |
CN109931952A (en) | The direct analytic expression coarse alignment method of inertial navigation under the conditions of unknown latitude | |
CN103900566A (en) | Method for eliminating influence on accuracy of rotating modulation strapdown inertial navigation system caused by earth rotation angular velocity | |
CN102636184B (en) | Specific force-sensitive term calibration method for flexible gyroscope based on centrifuge in environment without angular movement | |
CN103344252B (en) | A kind of Airborne Hyperspectral imaging system analysis of Positioning Error method | |
CN101832782B (en) | Method for quick field calibration of micro inertial measurement unit |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication |