CN113984090B - Wheel type robot IMU error online calibration and compensation method and device - Google Patents

Wheel type robot IMU error online calibration and compensation method and device Download PDF

Info

Publication number
CN113984090B
CN113984090B CN202111242117.7A CN202111242117A CN113984090B CN 113984090 B CN113984090 B CN 113984090B CN 202111242117 A CN202111242117 A CN 202111242117A CN 113984090 B CN113984090 B CN 113984090B
Authority
CN
China
Prior art keywords
imu
correction data
accelerometer
gyroscope
zero offset
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
CN202111242117.7A
Other languages
Chinese (zh)
Other versions
CN113984090A (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.)
University of Science and Technology Beijing USTB
Original Assignee
University of Science and Technology Beijing USTB
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 University of Science and Technology Beijing USTB filed Critical University of Science and Technology Beijing USTB
Priority to CN202111242117.7A priority Critical patent/CN113984090B/en
Publication of CN113984090A publication Critical patent/CN113984090A/en
Application granted granted Critical
Publication of CN113984090B publication Critical patent/CN113984090B/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
    • 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)
  • Manufacturing & Machinery (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Gyroscopes (AREA)

Abstract

The invention discloses a wheel robot IMU error online calibration and compensation method and device, wherein the method comprises the following steps: acquiring first correction data; the first correction data are a section of continuous data collected by the IMU when the wheeled robot provided with the IMU is in a standing state on a horizontal plane, and the continuous data comprise acceleration and angular velocity in the X, Y, Z direction; acquiring second correction data; the second correction data is a section of continuous data collected by the IMU in the in-situ rotation process of the wheeled robot, and comprises the angular speed in the X, Y, Z direction; establishing an accelerometer and gyroscope online calibration error model, and calculating calibration parameters; and compensating the output of the IMU through the calculated parameters. According to the invention, the condition that the IMU has an installation error is considered, and the IMU of the wheeled robot can be calibrated and compensated on line through simple operation, so that the measurement accuracy of the IMU is improved.

Description

Wheel type robot IMU error online calibration and compensation method and device
Technical Field
The invention relates to the technical field of automatic control of wheeled robots, in particular to a method and a device for calibrating and compensating an IMU error of a wheeled robot on line.
Background
Today, where service robots are becoming more and more popular, wheeled robots take up a significant role, while during their travel inertial measurement units (Inertial Measurement Unit, IMU) are a key sensor for the wheeled robots to achieve positioning, speed measurement, and attitude estimation. The IMU is used as a measuring element, and is sensitive to interference generated by unavoidable external factors such as impact, vibration, stress, temperature and the like in the processes of factory installation, transportation and user use, so that time-varying and non-negligible IMU errors can be caused. If the accumulation of the errors of the IMU is not processed, the gesture and navigation precision of the robot can be seriously reduced, so that the performance of the robot is affected. Therefore, on-line calibration and compensation are required for the IMU error, so that the measurement accuracy of the IMU is ensured, and the performance of the robot is ensured.
Before the IMU is installed on the wheeled robot, the IMU is generally fixed on a turntable for factory calibration, and calibration parameters comprise a scale factor and zero offset. Through testing, after the IMU leaves the factory and installs on the robot, the scale factors of the accelerometer and gyroscope are very small in change caused by factors such as installation, transportation, temperature change and aging, and the like, so that online calibration is not needed; however, the zero offset factors of gyroscopes and accelerometers can be greatly changed due to the factors, and online calibration and compensation are needed.
Generally, in the prior art, in order to perform on-line calibration on zero offset of a gyroscope and an accelerometer installed in a wheeled robot, the zero offset can be calculated by only allowing the robot to stand on a horizontal ground, and since the theoretical output of the gyroscope is 0 and the theoretical output of the accelerometer is the gravitational acceleration g.
However, considering that the mounting of the IMU on the wheeled robot is not horizontal, there is an installation error angle between the IMU and the robot, so even if the robot is placed on a horizontal ground, it is difficult for the theoretical horizontal axis in the IMU to remain parallel to the ground. And because of the existence of the gravity acceleration g, components of the gravity acceleration g exist on 3 axes of the accelerometer, which brings difficulty to the calculation of the zero offset of the accelerometer. It is difficult to estimate the accelerometer zero offset simply by resting on a level ground. Therefore, the existing method for calibrating the zero offset of the IMU installed in the wheeled robot cannot accurately realize the calibration and compensation of the IMU error of the wheeled robot.
Disclosure of Invention
The invention provides a method and a device for on-line calibration and compensation of an IMU error of a wheeled robot, which are used for solving the technical problem that the existing method for calibrating the zero offset of the IMU installed in the wheeled robot cannot accurately realize the calibration and compensation of the IMU error of the wheeled robot.
In order to solve the technical problems, the invention provides the following technical scheme:
in one aspect, the invention provides a wheel robot IMU error online calibration and compensation method, which comprises the following steps:
acquiring first correction data; the first correction data are a section of continuous data collected by the Inertial Measurement Unit (IMU) in a first preset time period when the wheeled robot provided with the IMU is in a standing state on a horizontal plane, and the first correction data comprise acceleration and angular velocity in the X, Y, Z direction;
acquiring second correction data; the second correction data is a section of continuous data collected by the IMU in the process that the wheeled robot rotates in situ for a preset angle on a horizontal plane at a preset speed within a second preset time period, and the second correction data comprises an angular speed in the X, Y, Z direction;
establishing an accelerometer and gyroscope online calibration error model;
based on the first correction data and the second correction data, calculating zero offset of the gyroscope and the accelerometer through an online calibration error model of the accelerometer and the gyroscope, and completing online calibration of the IMU error;
and compensating the output data of the IMU through the calculated zero offset of the gyroscope and the accelerometer.
Further, the expression of the accelerometer and gyroscope online calibration error model is as follows:
Figure BDA0003319574870000021
wherein,,
Figure BDA0003319574870000022
represents the angular velocity measurement value, ω represents the angular velocity theoretical value, w g Representing gyro noise of the type Gaussian white noise +.>
Figure BDA0003319574870000023
Represents the measured acceleration value, a represents the theoretical acceleration value, w a Representing accelerometer noise of the type Gaussian white noise, b g Indicating zero bias of gyroscope, b a Indicating the accelerometer zero offset.
Further, the calculating, based on the first correction data and the second correction data, zero offset of the gyroscope and the accelerometer through the accelerometer and the gyroscope on-line calibration error model includes:
calculating an average value of accelerations in the first correction data
Figure BDA0003319574870000024
And the average value of the angular velocity in the first correction data +.>
Figure BDA0003319574870000025
Since the IMU is in a rest state and gaussian white noise mathematical expectation is 0, ω=0, gyroscope noise average +.>
Figure BDA0003319574870000026
Accelerometer noise mean>
Figure BDA0003319574870000027
This gives:
Figure BDA0003319574870000028
wherein,,
Figure BDA0003319574870000031
representing the average value of zero offset of the gyroscope within a first preset duration;
calculating an average value of the angular velocity in the second correction data
Figure BDA0003319574870000032
Since the robot rotates on the horizontal plane, the rotation axis direction is the same as the gravity direction, and there are:
Figure BDA0003319574870000033
and because:
Figure BDA0003319574870000034
so that:
Figure BDA0003319574870000035
further, the compensating the output data of the IMU by the calculated zero offset of the gyroscope and the accelerometer includes:
based on the calculated zero offset of the gyroscope and the accelerometer and the online calibration error model of the accelerometer and the gyroscope, the output data of the IMU is compensated by using the following formula:
Figure BDA0003319574870000036
wherein omega 2 A, outputting data for the compensated gyroscope 2 And outputting data for the compensated accelerometer.
On the other hand, the invention also provides a wheel type robot IMU error on-line calibration and compensation device, which comprises:
the first correction data acquisition module is used for acquiring first correction data; the first correction data are a section of continuous data collected by the Inertial Measurement Unit (IMU) in a first preset time period when the wheeled robot provided with the IMU is in a standing state on a horizontal plane, and the first correction data comprise acceleration and angular velocity in the X, Y, Z direction;
a second correction data acquisition module configured to acquire second correction data; the second correction data is a section of continuous data collected by the IMU in the process that the wheeled robot rotates in situ for a preset angle on a horizontal plane at a preset speed within a second preset time period, and the second correction data comprises an angular speed in the X, Y, Z direction;
the calibration error model construction module is used for establishing an accelerometer and gyroscope online calibration error model;
the on-line calibration module is used for calculating zero offset of the gyroscope and the accelerometer through the on-line calibration error model of the accelerometer and the gyroscope established by the calibration error model establishment module based on the first correction data acquired by the first correction data acquisition module and the second correction data acquired by the second correction data acquisition module, so as to complete on-line calibration of the IMU error;
and the compensation module is used for compensating the output data of the IMU through the zero offset of the gyroscope and the accelerometer calculated by the online calibration module.
Further, the expression of the online calibration error model established by the calibration error model establishing module is as follows:
Figure BDA0003319574870000041
wherein,,
Figure BDA0003319574870000042
represents the angular velocity measurement value, ω represents the angular velocity theoretical value, w g Representing gyro noise of the type Gaussian white noise +.>
Figure BDA0003319574870000043
Represents the measured acceleration value, a represents the theoretical acceleration value, w a Representing accelerometer noise of the type Gaussian white noise, b g Indicating zero bias of gyroscope, b a Indicating the accelerometer zero offset.
Further, the online calibration module is specifically configured to:
calculating an average value of accelerations in the first correction data
Figure BDA0003319574870000044
And the average value of the angular velocity in the first correction data +.>
Figure BDA0003319574870000045
Since the IMU is in a rest state and gaussian white noise mathematical expectation is 0, ω=0, gyroscope noise average +.>
Figure BDA0003319574870000046
Accelerometer noise mean>
Figure BDA0003319574870000047
This gives:
Figure BDA0003319574870000048
wherein,,
Figure BDA0003319574870000049
representing the average value of zero offset of the gyroscope within a first preset duration;
calculating an average value of the angular velocity in the second correction data
Figure BDA00033195748700000410
Since the robot rotates on the horizontal plane, the rotation axis direction is the same as the gravity direction, and there are:
Figure BDA00033195748700000411
and because:
Figure BDA00033195748700000412
so that:
Figure BDA00033195748700000413
further, the compensation module is specifically configured to:
based on the calculated zero offset of the gyroscope and the accelerometer and the online calibration error model of the accelerometer and the gyroscope, the output data of the IMU is compensated by using the following formula:
Figure BDA00033195748700000414
wherein omega 2 A, outputting data for the compensated gyroscope 2 And outputting data for the compensated accelerometer.
In yet another aspect, the present invention also provides an electronic device including a processor and a memory; wherein the memory stores at least one instruction that is loaded and executed by the processor to implement the above-described method.
In yet another aspect, the present invention also provides a computer readable storage medium having at least one instruction stored therein, the instruction being loaded and executed by a processor to implement the above method.
The technical scheme provided by the invention has the beneficial effects that at least:
1. according to the invention, the on-line calibration of the IMU of the wheeled robot can be completed through real-time acquisition calculation, the calibration result is good, and the positioning and attitude-fixing performance can be effectively improved.
2. The mathematical formula in the invention is simple and easy to understand, and the program is convenient to realize; the implementation and operation are simple, the implementation is easy, and the operation time is short. Only one horizontal ground is needed for operation, and other specialized tools are not needed.
3. The wheel type robot IMU error online calibration and compensation method provided by the invention is suitable for the condition of larger IMU installation error and has good practicability.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings required for the description of the embodiments will be briefly described below, and it is apparent that the drawings in the following description are only some embodiments of the present invention, and other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
Fig. 1 is a schematic diagram of an execution flow of a method for on-line calibration and compensation of an IMU error of a wheeled robot according to an embodiment of the present invention.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the present invention more apparent, the embodiments of the present invention will be described in further detail with reference to the accompanying drawings.
First embodiment
The embodiment provides a wheel robot IMU error online calibration and compensation method, which can be realized by electronic equipment. The execution flow of the method is shown in fig. 1, and comprises the following steps:
s1, acquiring first correction data; the first correction data are a section of continuous data collected by the Inertial Measurement Unit (IMU) in a first preset time period when the wheeled robot provided with the IMU is in a standing state on a horizontal plane, and the first correction data comprise acceleration and angular velocity in the X, Y, Z direction;
in this embodiment, the first correction data is obtained by: placing the wheeled robot provided with the IMU on a horizontal plane, enabling the wheeled robot to be in a standing state, and then acquiring data output by the IMU; the data output by the IMU is IMU data corrected by the factory calibration coefficient.
S2, acquiring second correction data; the second correction data is a section of continuous data collected by the IMU in the process that the wheeled robot rotates in situ for a preset angle on a horizontal plane at a preset speed within a second preset time period, and the second correction data comprises an angular speed in the X, Y, Z direction;
in this embodiment, the second correction data is obtained by: placing a wheeled robot provided with an IMU on a horizontal plane, driving the robot to rotate in situ at a certain angle on the horizontal plane clockwise (or anticlockwise) at a certain speed, and acquiring a section of continuous data collected by the IMU in a in-situ rotation state; the data output by the IMU is IMU data corrected by the factory calibration coefficient.
S3, establishing an accelerometer and gyroscope online calibration error model;
in the embodiment, when modeling the error models of the accelerometer and the gyroscope, the calibration factors of the gyroscope and the accelerometer are considered to be calibrated when leaving the factory, and the changes caused by factors such as installation, transportation, temperature change, aging and the like after leaving the factory are very small, so that online calibration is not needed; however, the zero offset factors of gyroscopes and accelerometers can cause large changes due to the factors, and online calibration is needed.
Based on the above, the present embodiment builds an accelerometer and gyroscope error model as follows:
Figure BDA0003319574870000061
wherein,,
Figure BDA0003319574870000062
represents the angular velocity measurement value, ω represents the angular velocity theoretical value, w g Representing gyro noise of the type Gaussian white noise +.>
Figure BDA0003319574870000063
Represents the measured acceleration value, a represents the theoretical acceleration value, w a Representing accelerometer noise of the type Gaussian white noise, b g Indicating zero bias of gyroscope, b a Indicating the accelerometer zero offset.
S4, based on the first correction data and the second correction data, calculating zero offset of the gyroscope and the accelerometer through an online calibration error model of the accelerometer and the gyroscope, and completing online calibration of IMU errors;
specifically, in the present embodiment, the error parameter solving process of the accelerometer and the gyroscope is as follows:
and respectively averaging the acceleration and the angular velocity in the first correction data to obtain: average value of acceleration in first correction data
Figure BDA0003319574870000064
And an average value +.of the angular velocity in the first correction data>
Figure BDA0003319574870000065
Since the IMU is in a rest state and gaussian white noise mathematical expectation is 0, ω=0, gyroscope noise average +.>
Figure BDA0003319574870000066
Accelerometer noise average
Figure BDA0003319574870000067
This gives:
Figure BDA0003319574870000068
Figure BDA0003319574870000069
is the local gravitational acceleration.
Wherein,,
Figure BDA00033195748700000610
and representing the average value of the zero offset of the gyroscope in the first preset time period.
Averaging the angular velocity of the second correction data to obtain an average value of the angular velocity in the second correction data
Figure BDA00033195748700000611
Since the robot rotates on the horizontal plane, the rotation axis direction is the same as the gravity direction, and there are:
Figure BDA00033195748700000612
and because gaussian white noise is mathematically expected to be 0,
Figure BDA0003319574870000071
the error model is brought into the above equation:
Figure BDA0003319574870000072
and because:
Figure BDA0003319574870000073
so that:
Figure BDA0003319574870000074
so far, the zero offset of the gyroscope and the accelerometer is solved, and the online calibration of the IMU error of the robot is completed.
And S5, compensating the output data of the IMU through the calculated zero offset of the gyroscope and the accelerometer.
Specifically, when the robot operates in practical application, the IMU error compensation method is as follows:
Figure BDA0003319574870000075
wherein omega 2 A, outputting data for the compensated gyroscope 2 Outputting data for the compensated accelerometer; b g B is zero offset error of gyroscope a For zero offset error of the accelerometer, w g Data noise of gyroscope is Gaussian white noise, w a The data noise of the accelerometer is Gaussian white noise;
Figure BDA0003319574870000076
output data for the factory corrected gyroscope and accelerometer, respectively, can be represented by the following formulas:
Figure BDA0003319574870000077
wherein K is g Is 3-order square matrix, B g The 3-dimensional column vectors are respectively a factory calibrated gyroscope scale coefficient and zero offset; k (K) a Is 3-order square matrix, B a And 3-dimensional column vectors are respectively the calibration coefficient and zero offset of the factory calibrated accelerometer.
Figure BDA0003319574870000078
Raw data output by the gyroscope and the accelerometer respectively.
In summary, in the method for calibrating and compensating the wheel robot IMU error online in the embodiment, zero offset is solved by using the average value of the output of the gyroscope when the robot stands still; the zero offset of the accelerometer is solved by utilizing the characteristic that the direction of the rotation axis is the same as the direction of gravity when the robot rotates on the horizontal plane. The method can complete on-line calibration and compensation of the IMU of the wheeled robot through simple operation and real-time acquisition calculation, has good result and can effectively improve the positioning and attitude-fixing performance; the mathematical formula is simple and easy to understand, the program is convenient to realize, the implementation and operation are simple and easy to realize, the operation time is short, the operation can be performed only by a horizontal ground, and other professional tools are not needed. Moreover, the calibration and compensation method is suitable for the condition of large IMU installation error, and has good practicability.
Second embodiment
The embodiment provides an on-line calibration and compensation device for IMU errors of a wheeled robot, which comprises:
the first correction data acquisition module is used for acquiring first correction data; the first correction data are a section of continuous data collected by the Inertial Measurement Unit (IMU) in a first preset time period when the wheeled robot provided with the IMU is in a standing state on a horizontal plane, and the first correction data comprise acceleration and angular velocity in the X, Y, Z direction;
a second correction data acquisition module configured to acquire second correction data; the second correction data is a section of continuous data collected by the IMU in the process that the wheeled robot rotates in situ for a preset angle on a horizontal plane at a preset speed within a second preset time period, and the second correction data comprises an angular speed in the X, Y, Z direction;
the calibration error model construction module is used for establishing an accelerometer and gyroscope online calibration error model;
the on-line calibration module is used for calculating zero offset of the gyroscope and the accelerometer through the on-line calibration error model of the accelerometer and the gyroscope established by the calibration error model establishment module based on the first correction data acquired by the first correction data acquisition module and the second correction data acquired by the second correction data acquisition module, so as to complete on-line calibration of the IMU error;
and the compensation module is used for compensating the output data of the IMU through the zero offset of the gyroscope and the accelerometer calculated by the online calibration module.
The wheel type robot IMU error online calibration and compensation device of the embodiment corresponds to the wheel type robot IMU error online calibration and compensation method of the first embodiment; the functions realized by the functional modules in the wheel robot IMU error online calibration and compensation device in the embodiment are in one-to-one correspondence with the flow steps in the wheel robot IMU error online calibration and compensation method; therefore, the description is omitted here.
Third embodiment
The embodiment provides an electronic device, which comprises a processor and a memory; wherein the memory stores at least one instruction that is loaded and executed by the processor to implement the method of the first embodiment.
The electronic device may vary considerably in configuration or performance and may include one or more processors (central processing units, CPU) and one or more memories having at least one instruction stored therein that is loaded by the processors and performs the methods described above.
Fourth embodiment
The present embodiment provides a computer-readable storage medium having stored therein at least one instruction that is loaded and executed by a processor to implement the method of the first embodiment described above. The computer readable storage medium may be, among other things, ROM, random access memory, CD-ROM, magnetic tape, floppy disk, optical data storage device, etc. The instructions stored therein may be loaded by a processor in the terminal and perform the methods described above.
Furthermore, it should be noted that the present invention can be provided as a method, an apparatus, or a computer program product. Accordingly, embodiments of the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, embodiments of the invention may take the form of a computer program product on one or more computer-usable storage media having computer-usable program code embodied therein.
Embodiments of the present invention are described with reference to flowchart illustrations and/or block diagrams of methods, terminal devices (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, embedded processor, or other programmable data processing terminal device to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing terminal device, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks. These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
It is further noted that relational terms such as first and second, and the like are used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions. The terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or terminal that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or terminal. Without further limitation, an element defined by the phrase "comprising one … …" does not exclude the presence of other like elements in a process, method, article or terminal device comprising the element.
It is finally pointed out that the above description of the preferred embodiments of the invention, it being understood that although preferred embodiments of the invention have been described, it will be obvious to those skilled in the art that, once the basic inventive concepts of the invention are known, several modifications and adaptations can be made without departing from the principles of the invention, and these modifications and adaptations are intended to be within the scope of the invention. It is therefore intended that the following claims be interpreted as including the preferred embodiment and all such alterations and modifications as fall within the scope of the embodiments of the invention.

Claims (4)

1. The method for calibrating and compensating the wheel type robot IMU error on line is characterized by comprising the following steps:
acquiring first correction data; the first correction data are a section of continuous data collected by the Inertial Measurement Unit (IMU) in a first preset time period when the wheeled robot provided with the IMU is in a standing state on a horizontal plane, and the first correction data comprise acceleration and angular velocity in the X, Y, Z direction;
acquiring second correction data; the second correction data is a section of continuous data collected by the IMU in the process that the wheeled robot rotates in situ for a preset angle on a horizontal plane at a preset speed within a second preset time period, and the second correction data comprises an angular speed in the X, Y, Z direction;
establishing an accelerometer and gyroscope online calibration error model;
based on the first correction data and the second correction data, calculating zero offset of the gyroscope and the accelerometer through an online calibration error model of the accelerometer and the gyroscope, and completing online calibration of the IMU error;
compensating the output data of the IMU through the calculated zero offset of the gyroscope and the accelerometer;
the expression of the accelerometer and gyroscope on-line calibration error model is as follows:
Figure FDA0004262079460000011
wherein,,
Figure FDA0004262079460000012
represents the angular velocity measurement value, ω represents the angular velocity theoretical value, w g Representing gyro noise of the type Gaussian white noise +.>
Figure FDA0004262079460000013
Represents the measured acceleration value, a represents the theoretical acceleration value, w a Representing accelerometer noise of the type Gaussian white noise, b g Indicating zero bias of gyroscope, b a Indicating the accelerometer zero offset;
based on the first correction data and the second correction data, calculating zero offset of the gyroscope and the accelerometer through the accelerometer and the gyroscope on-line calibration error model, including:
calculating an average value of accelerations in the first correction data
Figure FDA0004262079460000014
And the average value of the angular velocity in the first correction data +.>
Figure FDA0004262079460000015
Since the IMU is in a rest state and gaussian white noise mathematical expectation is 0, ω=0, gyroscope noise average +.>
Figure FDA0004262079460000016
Accelerometer noise mean>
Figure FDA0004262079460000017
This gives:
Figure FDA0004262079460000018
wherein,,
Figure FDA0004262079460000019
representing the average value of zero offset of the gyroscope within a first preset duration;
calculating an average value of the angular velocity in the second correction data
Figure FDA00042620794600000110
Since the robot rotates on the horizontal plane, the rotation axis direction is the same as the gravity direction, and there are:
Figure FDA00042620794600000111
and because:
Figure FDA0004262079460000021
so that:
Figure FDA0004262079460000022
2. the method for on-line calibration and compensation of an IMU error of a wheeled robot according to claim 1, wherein said compensating the output data of said IMU by the calculated zero offset of the gyroscope and the accelerometer comprises:
based on the calculated zero offset of the gyroscope and the accelerometer and the online calibration error model of the accelerometer and the gyroscope, the output data of the IMU is compensated by using the following formula:
Figure FDA0004262079460000023
wherein omega 2 A, outputting data for the compensated gyroscope 2 To compensate forThe latter accelerometer outputs data.
3. The utility model provides a wheeled robot IMU error is demarcating and compensation arrangement on line which characterized in that includes:
the first correction data acquisition module is used for acquiring first correction data; the first correction data are a section of continuous data collected by the Inertial Measurement Unit (IMU) in a first preset time period when the wheeled robot provided with the IMU is in a standing state on a horizontal plane, and the first correction data comprise acceleration and angular velocity in the X, Y, Z direction;
a second correction data acquisition module configured to acquire second correction data; the second correction data is a section of continuous data collected by the IMU in the process that the wheeled robot rotates in situ for a preset angle on a horizontal plane at a preset speed within a second preset time period, and the second correction data comprises an angular speed in the X, Y, Z direction;
the calibration error model construction module is used for establishing an accelerometer and gyroscope online calibration error model;
the on-line calibration module is used for calculating zero offset of the gyroscope and the accelerometer through the on-line calibration error model of the accelerometer and the gyroscope established by the calibration error model establishment module based on the first correction data acquired by the first correction data acquisition module and the second correction data acquired by the second correction data acquisition module, so as to complete on-line calibration of the IMU error;
the compensation module is used for compensating the output data of the IMU through the zero offset of the gyroscope and the accelerometer calculated by the online calibration module;
the expression of the online calibration error model established by the calibration error model construction module is as follows:
Figure FDA0004262079460000024
wherein,,
Figure FDA0004262079460000031
represents the angular velocity measurement value, ω represents the angular velocity theoretical value, w g Representing gyro noise of the type Gaussian white noise +.>
Figure FDA0004262079460000032
Represents the measured acceleration value, a represents the theoretical acceleration value, w a Representing accelerometer noise of the type Gaussian white noise, b g Indicating zero bias of gyroscope, b a Indicating the accelerometer zero offset;
the on-line calibration module is specifically used for:
calculating an average value of accelerations in the first correction data
Figure FDA0004262079460000033
And the average value of the angular velocity in the first correction data +.>
Figure FDA0004262079460000034
Since the IMU is in a rest state and gaussian white noise mathematical expectation is 0, ω=0, gyroscope noise average +.>
Figure FDA0004262079460000035
Accelerometer noise mean>
Figure FDA0004262079460000036
This gives:
Figure FDA0004262079460000037
wherein,,
Figure FDA0004262079460000038
representing the average value of zero offset of the gyroscope within a first preset duration;
calculating the angle in the second correction dataAverage value of speed
Figure FDA0004262079460000039
Since the robot rotates on the horizontal plane, the rotation axis direction is the same as the gravity direction, and there are:
Figure FDA00042620794600000310
and because:
Figure FDA00042620794600000311
so that:
Figure FDA00042620794600000312
4. the on-line calibration and compensation device for the error of the IMU of the wheeled robot according to claim 3, wherein the compensation module is specifically configured to:
based on the calculated zero offset of the gyroscope and the accelerometer and the online calibration error model of the accelerometer and the gyroscope, the output data of the IMU is compensated by using the following formula:
Figure FDA00042620794600000313
wherein omega 2 A, outputting data for the compensated gyroscope 2 And outputting data for the compensated accelerometer.
CN202111242117.7A 2021-10-25 2021-10-25 Wheel type robot IMU error online calibration and compensation method and device Active CN113984090B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111242117.7A CN113984090B (en) 2021-10-25 2021-10-25 Wheel type robot IMU error online calibration and compensation method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111242117.7A CN113984090B (en) 2021-10-25 2021-10-25 Wheel type robot IMU error online calibration and compensation method and device

Publications (2)

Publication Number Publication Date
CN113984090A CN113984090A (en) 2022-01-28
CN113984090B true CN113984090B (en) 2023-07-04

Family

ID=79741093

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111242117.7A Active CN113984090B (en) 2021-10-25 2021-10-25 Wheel type robot IMU error online calibration and compensation method and device

Country Status (1)

Country Link
CN (1) CN113984090B (en)

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101788305A (en) * 2010-03-26 2010-07-28 中北大学 Method for rapid field calibration of micro inertial measurement unit
CN101949710A (en) * 2010-07-28 2011-01-19 北京泰豪联星技术有限公司 Rapid online dynamic calibration method for zero offset of GNSS (Global Navigation Satellite System) auxiliary MEMS (Micro Electro Mechanical Systems) inertial sensor
CN103278177A (en) * 2013-04-27 2013-09-04 中国人民解放军国防科学技术大学 Calibration method of inertial measurement unit based on camera network measurement
CN103776451A (en) * 2014-03-04 2014-05-07 哈尔滨工业大学 High-precision three-dimensional posture inertia measurement system and method based on MEMS (Micro Electro Mechanical Systems)
CN107576334A (en) * 2016-07-04 2018-01-12 北京合众思壮科技股份有限公司 The scaling method and device of Inertial Measurement Unit
CN107727114A (en) * 2017-09-20 2018-02-23 千寻位置网络有限公司 Acceleration scaling method and system, service terminal, memory based on gyroscope
WO2018129335A1 (en) * 2017-01-06 2018-07-12 Idhl Holdings, Inc. Methods, apparatus and systems for gyroscope scale calibration
CN108534800A (en) * 2018-03-09 2018-09-14 中国科学院长春光学精密机械与物理研究所 A kind of MEMS-IMU warm population parameter calibration compensation method entirely
CN109596143A (en) * 2018-11-07 2019-04-09 深圳市欧盛自动化有限公司 Test method, system, equipment and the computer readable storage medium of gyroscope
WO2019196093A1 (en) * 2018-04-13 2019-10-17 深圳市固胜智能科技有限公司 Pan tilt head-based imu calibration method and device and storage medium
CN111272199A (en) * 2020-03-23 2020-06-12 北京爱笔科技有限公司 Method and device for calibrating installation error angle of IMU (inertial measurement Unit)
CN112792821A (en) * 2021-04-01 2021-05-14 北京科技大学 Method and system for deblurring of vision of moving robot under inertia assisted facilitation exposure

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101788305A (en) * 2010-03-26 2010-07-28 中北大学 Method for rapid field calibration of micro inertial measurement unit
CN101949710A (en) * 2010-07-28 2011-01-19 北京泰豪联星技术有限公司 Rapid online dynamic calibration method for zero offset of GNSS (Global Navigation Satellite System) auxiliary MEMS (Micro Electro Mechanical Systems) inertial sensor
CN103278177A (en) * 2013-04-27 2013-09-04 中国人民解放军国防科学技术大学 Calibration method of inertial measurement unit based on camera network measurement
CN103776451A (en) * 2014-03-04 2014-05-07 哈尔滨工业大学 High-precision three-dimensional posture inertia measurement system and method based on MEMS (Micro Electro Mechanical Systems)
CN107576334A (en) * 2016-07-04 2018-01-12 北京合众思壮科技股份有限公司 The scaling method and device of Inertial Measurement Unit
WO2018129335A1 (en) * 2017-01-06 2018-07-12 Idhl Holdings, Inc. Methods, apparatus and systems for gyroscope scale calibration
CN107727114A (en) * 2017-09-20 2018-02-23 千寻位置网络有限公司 Acceleration scaling method and system, service terminal, memory based on gyroscope
CN108534800A (en) * 2018-03-09 2018-09-14 中国科学院长春光学精密机械与物理研究所 A kind of MEMS-IMU warm population parameter calibration compensation method entirely
WO2019196093A1 (en) * 2018-04-13 2019-10-17 深圳市固胜智能科技有限公司 Pan tilt head-based imu calibration method and device and storage medium
CN109596143A (en) * 2018-11-07 2019-04-09 深圳市欧盛自动化有限公司 Test method, system, equipment and the computer readable storage medium of gyroscope
CN111272199A (en) * 2020-03-23 2020-06-12 北京爱笔科技有限公司 Method and device for calibrating installation error angle of IMU (inertial measurement Unit)
CN112792821A (en) * 2021-04-01 2021-05-14 北京科技大学 Method and system for deblurring of vision of moving robot under inertia assisted facilitation exposure

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
A calibration and error analysis method for IMU using centrifuge;Zhao yan;《transactions of beijing institute of technology》;第35卷(第4期);414-420 *
A method of low-cost IMU calibration and alignment;JIXIN LV;《2016 IEEE/SICE international symposium on system integration(SII)》;373-378 *
IMU在卧式三轴转台上的标定方法研究;崔凯鑫;《中国优秀硕士学位论文全文数据库信息科技辑》(第2期);I136-2097 *
单轴旋转式光纤惯导系统技术研究;张思;《中国优秀硕士学位论文全文数据库信息科技辑》(第3期);I136-1504 *

Also Published As

Publication number Publication date
CN113984090A (en) 2022-01-28

Similar Documents

Publication Publication Date Title
KR101988786B1 (en) Initial alignment of inertial navigation devices
CN107655493B (en) SINS six-position system-level calibration method for fiber-optic gyroscope
CN109459061B (en) Micro-inertia measurement unit calibration method, device and computer readable storage medium
CN108592952A (en) The method for demarcating more MIMU errors simultaneously with positive and negative times of rate based on lever arm compensation
CN105806367B (en) Gyro free inertia system error calibrating method
CN111712688A (en) Calibration method, calibration apparatus, stabilizer, and computer-readable storage medium
CN111551174A (en) High-dynamic vehicle attitude calculation method and system based on multi-sensor inertial navigation system
CN106153069B (en) Attitude rectification device and method in autonomous navigation system
KR101106048B1 (en) Method for calibrating sensor errors automatically during operation, and inertial navigation using the same
CN109084806A (en) Scalar domain MEMS inertia system scaling method
CN110887481A (en) Carrier dynamic attitude estimation method based on MEMS inertial sensor
EP3527948B1 (en) Air data aided inertial measurement unit
CN110887480A (en) Flight attitude estimation method and system based on MEMS sensor
CN110361031A (en) A kind of IMU population parameter error quick calibrating method theoretical based on backtracking
CN114413934A (en) Vehicle positioning system correction method and device
CN113984090B (en) Wheel type robot IMU error online calibration and compensation method and device
EP2853901B1 (en) Telematics system and associated method
CN117268431A (en) Automatic calibration method, device and system for angle sensor
JPH04231813A (en) Method for measuring angle and angular characteristic curve
CN114964214B (en) Extended Kalman filtering attitude calculation method of attitude heading reference system
CN113959464B (en) Gyroscope-assisted accelerometer field calibration method and system
CN115876225A (en) MEMS IMU calibration method and system based on two-degree-of-freedom turntable
CN108716925A (en) A kind of scaling method and device of nine axle sensors
US11371848B2 (en) Method for characterising an inertial measurement unit
RU2727344C1 (en) Method for increasing calibration accuracy of unit of angular velocity micromechanical sensors

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