CN115060259A - Multi-IMU fusion positioning system and method - Google Patents

Multi-IMU fusion positioning system and method Download PDF

Info

Publication number
CN115060259A
CN115060259A CN202210996020.3A CN202210996020A CN115060259A CN 115060259 A CN115060259 A CN 115060259A CN 202210996020 A CN202210996020 A CN 202210996020A CN 115060259 A CN115060259 A CN 115060259A
Authority
CN
China
Prior art keywords
imu
fusion
virtual
fused
imus
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
Application number
CN202210996020.3A
Other languages
Chinese (zh)
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.)
Jiangsu Deyijia Security Technology Co ltd
Original Assignee
Jiangsu Deyijia Security Technology Co ltd
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 Jiangsu Deyijia Security Technology Co ltd filed Critical Jiangsu Deyijia Security Technology Co ltd
Priority to CN202210996020.3A priority Critical patent/CN115060259A/en
Publication of CN115060259A publication Critical patent/CN115060259A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/10Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration
    • G01C21/12Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning
    • G01C21/16Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/20Instruments for performing navigational calculations

Abstract

The invention provides a positioning system and a method for fusing multiple IMUs, which comprises the following steps: the system comprises a multi-IMU fusion module and a position and speed updating module; a multi-IMU fusion module: the virtual fusion IMU is obtained according to the measurement data of the IMU; a location speed update module: and the IMU is used for positioning according to the virtual fusion IMU. Compared with the prior art, the method improves the reliability and robustness of positioning through the integrated fusion of multiple IMUs.

Description

Multi-IMU fusion positioning system and method
Technical Field
The invention relates to the technical field of positioning, in particular to a positioning system and a method with multiple IMUs integrated.
Background
The fire fighters have complex and dangerous operation environments in the fire fighting process, and the importance of accurate positioning and tracking of the fire fighters is increasingly prominent for guaranteeing the life safety and property safety of the fire fighters and the national people. At present, there are various methods for positioning fire-fighting operators, and common positioning methods include an Inertial Measurement Unit (IMU), a Global Navigation Satellite System (GNSS), and various indoor positioning technologies.
An Inertial Measurement Unit (IMU) is a low cost motion sensor that can measure angular velocity and gravity-compensated linear acceleration of a moving platform and is widely used in modern positioning systems. The indoor positioning system is independent of external information, and is calculated based on sensors such as acceleration and angular velocity of the indoor positioning system, so that information such as movement speed, direction and position is determined, and indoor positioning of people is basically based on IMU (inertial measurement unit), but the indoor positioning system has the following defects: errors accumulate over time and eventually become large over time, requiring correction by other external sources of navigation, which inevitably increases the complexity and cost of the system.
To date, existing inertial positioning methods utilize only one IMU. While a single IMU positioning may provide acceptable accuracy and robustness for different use cases, the following challenges remain for personnel positioning: the information which can be obtained by the judgment logic of zero-speed updating is limited, so that an accurate zero-speed updating algorithm cannot be implemented, and the error of position calculation is greatly reduced; the system redundancy of a single IMU is low, and the reliability is low for a complex operation environment; the error accumulation of a single IMU is fast, and the effective working time is short.
Patent document CN109827593A discloses an error self-calibration method, system and storage medium based on multiple IMUs, which includes the steps of: the method comprises the steps of conducting front and back surface pasting on n IMUs according to the principle that the variables with the same attribute are parallel, simplifying all errors of the n IMUs into parameter error models of three types of comprehensive errors according to output variables, removing the maximum and minimum output quantities of the IMUs based on statistical information, conducting zero offset error compensation on the IMUs with gross errors removed, directly adding the zero offset error quantities in output quantities, conducting output averaging of the compensated zero offset error values in the same property and direction, and outputting corresponding accelerometers, angular velocities and magnetic field velocities. However, this method does not effectively solve the problem of low reliability of the IMU in a complex working environment.
Disclosure of Invention
Aiming at the defects in the prior art, the invention aims to provide a positioning system and a positioning method for multi-IMU fusion.
The invention provides a positioning system with multiple IMUs fused, which comprises the following modules: the system comprises a multi-IMU fusion module and a position and speed updating module.
A multi-IMU fusion module: the virtual fusion IMU is obtained according to the measurement data of the IMU;
a location speed update module: and the positioning processing is carried out according to the virtual fusion IMU.
Preferably, the multi-IMU fusion module includes: multiple IMUs and multiple IMUs integrated processing units.
Preferably, the positioning process comprises: location update and velocity update.
Preferably, the plurality of IMUs are provided on the device in a preset manner.
Preferably, the apparatus is provided at a plurality of locations on the worker.
The positioning method for fusing multiple IMUs provided by the invention is suitable for a positioning system for fusing multiple IMUs, and comprises the following steps:
step 1: acquiring measurement data of each IMU;
step 2: obtaining a relative relation matrix according to the relative position relation of the IMUs;
and step 3: obtaining a virtual fusion IMU according to the measurement data and the relative relation matrix;
and 4, step 4: and performing positioning processing according to the virtual fusion IMU.
Preferably, the measurement data comprises angular velocity data and acceleration data;
step 1, comprising:
step 101: obtaining angular velocity data according to the gyroscope of each IMU;
step 102: acceleration data is obtained from the accelerometer of each IMU.
Preferably, step 4, comprises:
step 401: obtaining a virtual fusion gyroscope model according to the virtual fusion IMU;
step 402: obtaining a virtual fusion accelerometer model according to the virtual fusion IMU;
step 403: and carrying out positioning processing according to the virtual fusion gyroscope model and the virtual fusion accelerometer model.
Preferably, the virtual fused IMU is optimized by a least squares method.
Preferably, the positioning process comprises: location update and velocity update.
Compared with the prior art, the invention has the following beneficial effects:
1. the invention improves the reliability and robustness of the system through the integrated fusion of multiple IMUs.
2. According to the method, the updating precision of the IMU position speed is improved through the optimization solution of the probability error.
3. The invention constructs the virtual integrated IMU, so that the calculation amount of position and speed updating is not changed, and the development and design of a system are facilitated.
4. The positioning system with the multiple IMUs integrated can be mounted or worn on various parts and equipment of people, different combination schemes can be customized, and design and use of different scenes are facilitated.
5. The positioning system with the multiple IMUs integrated can avoid the situation that a single IMU can judge the walking and running gaits only by means of a single motion mode, so that the positioning accuracy of the position and the speed is further improved.
Drawings
Other features, objects and advantages of the invention will become more apparent upon reading of the detailed description of non-limiting embodiments with reference to the following drawings:
FIG. 1 is a schematic structural view of the present invention;
FIG. 2 is a flow diagram of a multi-IMU integrated processing unit of the present invention;
FIG. 3 is a schematic flow chart of the present invention;
FIG. 4 is a diagram illustrating the location and velocity update of the present invention.
Detailed Description
The present invention will be described in detail with reference to specific examples. The following examples will assist those skilled in the art in further understanding the invention, but are not intended to limit the invention in any way. It should be noted that it would be obvious to those skilled in the art that various changes and modifications can be made without departing from the spirit of the invention. All falling within the scope of the present invention.
Fig. 1 is a schematic structural diagram of the present invention, and as shown in fig. 1, the present invention provides a positioning system with multiple IMU fusion, which includes the following modules: the system comprises a multi-IMU fusion module and a position and speed updating module. The multi-IMU fusion module includes a plurality of IMUs, as shown in fig. 1, which may be represented as: IMU1, IMU2, …, IMUn; the multi-IMU fusion module also comprises a multi-IMU integrated processing unit, and the multi-IMU integrated processing unit receives and processes the measurement information of all IMUs, performs fusion optimization on the measurement information of the IMUs and constructs a virtual fusion IMU. The position and speed updating module comprises position and speed updating, and performs position updating and speed updating on the virtual fusion IMU established by the multi-IMU fusion module.
Wherein, many IMU fuse the module: the virtual fusion IMU is obtained according to the measurement data of the IMU; a location speed update module: and the positioning processing is carried out according to the virtual fusion IMU.
Specifically, in the present invention, the measurement probabilities of all physical IMUs are mapped onto one virtual fusion IMU and randomly estimated by using a least squares estimator. Accordingly, a classical filter-based or optimization-based sensor fusion algorithm may be used for the localization process.
Preferably, the multi-IMU fusion module includes: multiple IMUs and multiple IMUs integrated processing units.
Specifically, a plurality of IMUs are set on the device in a preset manner.
The preset mode is not limited in the present invention, and may be set according to specific situations, for example, the preset mode includes, but is not limited to, an array, a straight line, or a stack.
Preferably, the apparatus is provided at a plurality of locations on the worker.
Preferably, the positioning process comprises: location update and velocity update.
Specifically, the multi-IMU fusion module may include two or more IMUs. Moreover, the IMUs may be of the same model, or may be configured with different capabilities or different models.
Wherein the physical locations of the plurality of IMUs may be arranged on the same device, in particular, by means including, but not limited to, an array, a line, or a stack.
Further, the multi-IMU fusion module may be decomposed into a plurality of devices as needed, and each device may have a certain number of IMUs arranged thereon in a manner including, but not limited to, an array, a line, or a stack. And each device performs transmission of measurement data and characteristic parameters of the IMU and integration of multiple IMUs through a certain communication mechanism.
It will be appreciated that these IMU-lined devices may be carried independently or distributed about the body part of the firefighter, such as the feet, limbs, head or torso.
Specifically, the multi-IMU integrated processing unit is used for receiving and processing the measurement data of all IMUs, performing fusion optimization on multi-IMU information, and constructing a virtual fusion IMU for positioning and tracking of fire fighters.
In the invention, the multi-IMU integrated processing unit can be integrated with one, a plurality of or all IMUs on the same device, or can independently reside on a special device.
Each device transmits IMU measurement data and characteristic parameters through a certain communication mechanism and is used for supporting the integration and fusion of multiple IMUs.
Fig. 2 is a schematic flow diagram of a multi-IMU integrated processing unit according to the present invention, and as shown in fig. 2, firstly, the measurement data of each IMU is received, then, the validity of the measurement data of each IMU is verified, further, the measurement data of invalid IMU is removed, and finally, multi-IMU fusion optimization is performed according to the measurement data of IMU.
Fig. 3 is a schematic flow chart of the present invention, and as shown in fig. 3, the present invention provides a positioning method for multi-IMU fusion, which is suitable for the above positioning system for multi-IMU fusion, and includes the following steps:
step 1: measurement data for each IMU is acquired.
Wherein the measurement data comprises angular velocity data and acceleration data.
Specifically, step 1 may include: step 101: obtaining angular velocity data according to the gyroscope of each IMU; step 102: acceleration data is obtained from the accelerometer of each IMU.
Illustratively, there are n IMUs, and the measurement equation for the gyroscope of each IMU is shown in equation (1):
Figure 579189DEST_PATH_IMAGE001
(1)
wherein the content of the first and second substances,
Figure 965171DEST_PATH_IMAGE002
is shown as
Figure DEST_PATH_IMAGE003
Angular speed output values of a gyroscope of each IMU in three directions of an x axis, a y axis and a z axis under an IMU body coordinate system;
Figure 666280DEST_PATH_IMAGE004
denotes the first
Figure 530331DEST_PATH_IMAGE005
Angular velocity ideal values of a gyroscope of each IMU in three directions of an x axis, a y axis and a z axis under an IMU body coordinate system;
Figure 298698DEST_PATH_IMAGE006
representing random walk errors conforming to a normal distribution, i.e. Gaussian noise, in which the covariance
Figure 590002DEST_PATH_IMAGE007
Figure 60298DEST_PATH_IMAGE008
A zero matrix representing the corresponding dimension, typically provided by the IMU manufacturer;
Figure 728039DEST_PATH_IMAGE009
representing static measurement error;
Figure 115027DEST_PATH_IMAGE010
the noise representing the static measurement error, i.e. the dynamic measurement error, can be seen as a scalar normal markov process.
Wherein the measurement equation for the accelerometer of each IMU is shown in equation (2):
Figure 577232DEST_PATH_IMAGE011
;(2)
wherein the content of the first and second substances,
Figure 800403DEST_PATH_IMAGE012
is shown as
Figure 494340DEST_PATH_IMAGE013
Acceleration output values of an accelerometer of each IMU in three directions of an x axis, a y axis and a z axis under an IMU body coordinate system;
Figure 486566DEST_PATH_IMAGE014
is shown as
Figure 385252DEST_PATH_IMAGE013
Acceleration ideal values of an accelerometer of each IMU in three directions of an x axis, a y axis and a z axis under an IMU body coordinate system;
Figure 95719DEST_PATH_IMAGE015
representing random walk errors following a normal distribution, i.e. Gaussian noise, in which the covariance
Figure 88952DEST_PATH_IMAGE016
Figure 201265DEST_PATH_IMAGE017
The zero matrix representing the corresponding dimension is typically provided by the IMU manufacturer;
Figure 739693DEST_PATH_IMAGE018
representing static measurement error;
Figure 937456DEST_PATH_IMAGE019
the noise representing the static measurement error, i.e., the dynamic measurement error, can also be considered a scalar normal Markov process.
Step 2: and obtaining a relative relation matrix according to the relative position relation of the IMUs.
Specifically, determining
Figure 235845DEST_PATH_IMAGE020
Relative relationship matrix of position between IMUs
Figure 202664DEST_PATH_IMAGE021
Specifically, as shown in formula (3):
Figure 443152DEST_PATH_IMAGE022
;(3)
wherein the content of the first and second substances,
Figure 862632DEST_PATH_IMAGE023
is as follows
Figure 728826DEST_PATH_IMAGE024
IMU relative to
Figure 550152DEST_PATH_IMAGE025
The body coordinate transformation matrix of each IMU meets the formula (4):
Figure 695962DEST_PATH_IMAGE026
;(4)
wherein the content of the first and second substances,
Figure 602738DEST_PATH_IMAGE027
it should be noted that if the physical locations of multiple IMUs are arranged on the same device by way of including, but not limited to, an array, a line, or a stack, the relative relationship matrix is
Figure 974420DEST_PATH_IMAGE028
Is a constant matrix. If these physical IMUs are distributed across multiple devices or mounted at different locations on a firefighter, the relative relationship matrix
Figure 915832DEST_PATH_IMAGE028
Is a dynamic matrix.
And step 3: and obtaining the virtual fusion IMU according to the measurement data and the relative relation matrix.
In particular, based on this
Figure 498123DEST_PATH_IMAGE029
A physical IMU for constructing a virtual IMU with subscript defined as
Figure 157774DEST_PATH_IMAGE030
And defining a virtual body coordinate system, the virtual fusion IMU having a coordinate transformation matrix with respect to each physical IMU
Figure 834612DEST_PATH_IMAGE031
Figure 364950DEST_PATH_IMAGE032
. Here, the body coordinate system of the virtual fused IMU may be selected to be independent of all the physical IMUs, or may be consistent with the body coordinate system of one of the physical IMUs.
In general, the virtual fused IMU may be defined relative to a selected second IMU
Figure 852564DEST_PATH_IMAGE033
The relative position of the individual physical IMUs,
Figure 750243DEST_PATH_IMAGE034
. Thereby obtaining the virtual fused IMU relative to the second
Figure 247084DEST_PATH_IMAGE033
Coordinate transformation matrix of physical IMU
Figure 897508DEST_PATH_IMAGE035
If the virtual fusion IMU is consistent with the body coordinate system of the physical IMU, the coordinate transformation matrix
Figure 556023DEST_PATH_IMAGE035
Can be expressed by equation (5):
Figure 173955DEST_PATH_IMAGE036
;(5)
further, equation (6) is obtained:
Figure 208907DEST_PATH_IMAGE037
;(6)
and 4, step 4: and performing positioning processing according to the virtual fusion IMU.
Preferably, step 4, comprises: step 401: obtaining a virtual fusion gyroscope model according to the virtual fusion IMU;
step 402: obtaining a virtual fusion accelerometer model according to the virtual fusion IMU; step 403: and positioning processing is carried out according to the virtual fusion gyroscope model and the virtual fusion accelerometer model.
Preferably, the virtual fused IMU is optimized by a least squares method.
Wherein the positioning process comprises: location update and velocity update.
Specifically, the measurement equation of the gyroscope in the virtual fusion IMU is shown in equation (7):
Figure 979417DEST_PATH_IMAGE038
;(7)
wherein, the first and the second end of the pipe are connected with each other,
Figure 808832DEST_PATH_IMAGE039
representing angular velocity output values of a gyroscope with a virtual fusion IMU in three directions of an x axis, a y axis and a z axis under a body coordinate system;
Figure 141157DEST_PATH_IMAGE040
a vector of gyroscope output values representing all physical IMUs;
Figure 714221DEST_PATH_IMAGE041
to represent
Figure 73658DEST_PATH_IMAGE042
The generalized inverse matrix of (a), i.e.,
Figure 588822DEST_PATH_IMAGE043
Figure 932078DEST_PATH_IMAGE044
can be expressed by equation (8):
Figure 308833DEST_PATH_IMAGE045
;(8)
similarly, the measurement equation for the accelerometers of the virtual fusion IMU is shown in equation (9):
Figure 788356DEST_PATH_IMAGE046
;(9)
wherein the content of the first and second substances,
Figure 710307DEST_PATH_IMAGE047
representing the acceleration output values of the accelerometer of the virtual fusion IMU in three directions of an x axis, a y axis and a z axis under a body coordinate system of the accelerometer;
Figure 275280DEST_PATH_IMAGE048
a vector of output values representing the accelerometers of all physical IMUs.
As in the above, the above-mentioned,
Figure 721305DEST_PATH_IMAGE049
is composed of
Figure 55335DEST_PATH_IMAGE050
The generalized inverse matrix of (a), i.e.,
Figure 912301DEST_PATH_IMAGE051
Figure 964571DEST_PATH_IMAGE050
represented by formula (8).
Further, the virtual fusion gyroscope model is solved.
In particular, according to the virtual converged IMU, the institute can be integratedOutput value vector for a gyroscope with a physical IMU
Figure 683128DEST_PATH_IMAGE052
Rewritten as formula (10):
Figure 137243DEST_PATH_IMAGE053
;(10)
wherein the content of the first and second substances,
Figure 398067DEST_PATH_IMAGE054
then, the gaussian noise in the angular velocity measurement equation for each physical IMU is rewritten to be expressed on the respective axes as shown in equation (11):
Figure 937632DEST_PATH_IMAGE055
;(11)
wherein the content of the first and second substances,
Figure 459881DEST_PATH_IMAGE056
further, define
Figure 768502DEST_PATH_IMAGE057
Figure 967271DEST_PATH_IMAGE057
Expressed by equation (12):
Figure 994133DEST_PATH_IMAGE058
;(12)
then, the best estimate of angular velocity of the virtually fused IMU
Figure 320072DEST_PATH_IMAGE059
Satisfies formula (13):
Figure 217621DEST_PATH_IMAGE060
;(13)
further, equation (13) is solved by the least square method. For simplicity, the Gaussian noise in the angular velocity measurement equation for each physical IMU is assumed to be consistent across the various axes and is defined as
Figure 88756DEST_PATH_IMAGE061
That is to say that,
Figure 602914DEST_PATH_IMAGE062
similarly, if all the physical IMUs are of the same type, it may be further assumed that the gaussian noise of each IMU is consistent and defined as
Figure 466965DEST_PATH_IMAGE063
That is to say that,
Figure 484600DEST_PATH_IMAGE064
by solving equation (13) by the least square method and writing it in conformity with the measurement equation form of the physical IMU, it can be expressed as equation (14):
Figure 290750DEST_PATH_IMAGE065
;(14)
wherein the content of the first and second substances,
Figure 26625DEST_PATH_IMAGE066
and the gyroscope represents angular speed output values of the virtual fusion IMU in the directions of x, y and z under the IMU body coordinate system.
Figure 959946DEST_PATH_IMAGE067
The method comprises the following steps of virtually fusing ideal angular velocity values of a gyroscope of an IMU in three directions of an x axis, a y axis and a z axis under an IMU body coordinate system;
Figure 832087DEST_PATH_IMAGE068
representing random walk errors following a normal distribution, i.e. Gaussian noiseWherein the covariance
Figure 47955DEST_PATH_IMAGE069
Is obtained by solving through a least square method,
Figure 5546DEST_PATH_IMAGE070
a zero matrix of corresponding dimensions;
Figure 211400DEST_PATH_IMAGE071
representing static measurement errors and solving by a least square method;
Figure 452894DEST_PATH_IMAGE072
the noise representing the static measurement error, i.e., the dynamic measurement error, can be considered as a scalar normal markov process, and is solved by the least square method.
Then, the virtual fusion accelerometer model is solved according to the method.
In particular, according to the virtual fused IMU, the output values of the accelerometers of all physical IMUs can be vectorized
Figure 617159DEST_PATH_IMAGE073
Rewritten as formula (15):
Figure 62047DEST_PATH_IMAGE074
;(15)
wherein, the first and the second end of the pipe are connected with each other,
Figure 71591DEST_PATH_IMAGE075
then, the gaussian noise in the acceleration measurement equation for each physical IMU is rewritten for expression on the respective axis, as shown in equation (16):
Figure 183904DEST_PATH_IMAGE076
;(16)
further, define
Figure 473065DEST_PATH_IMAGE077
Figure 670828DEST_PATH_IMAGE078
Can be expressed as formula (17):
Figure 952905DEST_PATH_IMAGE079
;(17)
optimal estimation of virtual fused IMU acceleration
Figure 434571DEST_PATH_IMAGE080
Satisfies formula (18):
Figure 409480DEST_PATH_IMAGE081
;(18)
equation (18) will be solved by the least square method below. For simplicity, the Gaussian noise in the acceleration measurement equation for each physical IMU is assumed to be consistent across the axes and is defined as
Figure 94539DEST_PATH_IMAGE082
That is to say that,
Figure 445886DEST_PATH_IMAGE083
similarly, if all the physical IMUs are of the same type, it may be further assumed that the gaussian noise of each IMU is consistent and defined as
Figure 15014DEST_PATH_IMAGE084
That is to say that,
Figure 160825DEST_PATH_IMAGE085
solving by the least squares method, and writing it to be consistent with the form of the measurement equation of the physical IMU, equation (19) can be obtained:
Figure 67601DEST_PATH_IMAGE086
;(19)
wherein the content of the first and second substances,
Figure 957060DEST_PATH_IMAGE087
representing acceleration output values of an accelerometer virtually fused with the IMU in three directions of an x axis, a y axis and a z axis under an IMU body coordinate system;
Figure 882159DEST_PATH_IMAGE088
representing the ideal acceleration values of the accelerometer virtually fused with the IMU in the three directions of the x axis, the y axis and the z axis under the IMU body coordinate system;
Figure 198871DEST_PATH_IMAGE089
representing random walk errors following a normal distribution, i.e. Gaussian noise, in which the covariance
Figure 592943DEST_PATH_IMAGE090
Is obtained by solving through a least square method,
Figure 36825DEST_PATH_IMAGE091
a zero matrix of corresponding dimensions;
Figure 832743DEST_PATH_IMAGE092
representing static measurement errors and solving by a least square method;
Figure 585935DEST_PATH_IMAGE093
the noise representing the static measurement error, i.e. the dynamic measurement error, can also be regarded as a scalar normal Markov process, and is obtained by solving with the least square method.
It can be known that, through the above steps, a virtual fused IMU obtained based on fusion of a plurality of physical IMUs can be obtained, the form of the virtual fused IMU conforms to the corresponding characteristic parameters of a single physical IMU, and the position and the speed of the person who is equipped with the plurality of IMUs can be calculated through a conventional IMU propagation equation.
Preferably, the position and speed updating module updates the position and speed of the virtual fused IMU established by the multi-IMU fusion module.
Fig. 4 is a schematic diagram of position and velocity update according to the present invention, as shown in fig. 4, acceleration coordinate conversion and posture update are performed by an accelerometer and a gyroscope, velocity update and position cross flow are performed after acceleration coordinate conversion, gravity update and direction change rate update are performed after position update, direction change rate update is performed after velocity update, direction update and posture update are performed after direction change rate update, further, posture update data may be used for acceleration coordinate conversion, direction update data and gravity update data may be used for velocity update, and velocity update data and position update data are used for positioning processing.
The invention provides a positioning system and a positioning method for fusing multiple IMUs (inertial measurement units) to improve the positioning accuracy of firefighters.
It should be noted that, in the present invention, by introducing the dynamic relative relationship matrix of multiple IMUs, multiple IMUs may be fixedly loaded in the same positioning device, or may be distributed and loaded in each related positioning sub-device. The real-time dynamic relative relationship matrix among the IMUs is determined, so that on one hand, the precision of the constructed virtual integrated IMU can be improved; on the other hand, the gait, the movement posture and the like of the fire fighter in the operation process can be judged; therefore, the accuracy of the IMU zero-speed updating algorithm is improved, the precision of calculating the positions of the firefighters is improved, and auxiliary monitoring is carried out on the task execution state and the states of the firefighters. The multi-IMU combined navigation method is mostly seen in a common fusion algorithm, but a method for constructing a virtual integrated IMU based on a plurality of IMUs is not involved for positioning fire fighters, the influence of the common fusion algorithm on error accumulation is avoided, and the accuracy of updating the speed and the position of the virtual integrated IMU for a long time is ensured through the optimal solution of probability errors.
The technical principle of the invention is as follows:
the present invention fuses measurement data from multiple IMUs to enable significant performance improvements for a positioning system without incurring additional computational costs. The invention maps the measurement probabilities of all physical IMUs onto one virtual fusion IMU and performs stochastic estimation by using a least squares estimator. Accordingly, classical filter-based or optimization-based sensor fusion algorithms can be used for localization.
The invention mainly solves the technical problems that:
1. the information which can be obtained by the judgment logic of zero-speed updating is limited, so that an accurate zero-speed updating algorithm cannot be implemented, and the error of position calculation is greatly reduced.
2. The system redundancy of a single IMU is low, and the reliability is low for a complex operation environment.
3. The error accumulation of a single IMU is fast, and the effective working time is short.
4. The overall performance of the system may be further improved by using multiple IMUs, resulting in higher reliability and slower accumulation of errors over time.
Compared with the prior art, the invention has the following beneficial effects:
1. the invention improves the reliability and robustness of the system through the integrated fusion of multiple IMUs.
2. According to the method, the updating precision of the IMU position speed is improved through the optimization solution of the probability error.
3. The invention constructs the virtual integrated IMU, so that the calculation amount of position and speed updating is unchanged, and the development and design of the system are facilitated.
4. The positioning system with the multiple IMUs integrated can be mounted or worn on various parts and equipment of people, different combination schemes can be customized, and design and use of different scenes are facilitated.
5. The positioning system with the multiple IMUs integrated can avoid the situation that a single IMU can judge the walking and running gaits only by means of a single motion mode, so that the positioning accuracy of the position and the speed is further improved.
Those skilled in the art will appreciate that, in addition to implementing the systems, apparatus, and various modules thereof provided by the present invention in purely computer readable program code, the same procedures can be implemented entirely by logically programming method steps such that the systems, apparatus, and various modules thereof are provided in the form of logic gates, switches, application specific integrated circuits, programmable logic controllers, embedded microcontrollers and the like. Therefore, the system, the device and the modules thereof provided by the present invention can be considered as a hardware component, and the modules included in the system, the device and the modules thereof for implementing various programs can also be considered as structures in the hardware component; modules for performing various functions may also be considered to be both software programs for performing the methods and structures within hardware components.
The foregoing description of specific embodiments of the present invention has been presented. It is to be understood that the present invention is not limited to the specific embodiments described above, and that various changes or modifications may be made by one skilled in the art within the scope of the appended claims without departing from the spirit of the invention. The embodiments and features of the embodiments of the present application may be combined with each other arbitrarily without conflict.

Claims (10)

1. A multi-IMU fused localization system, the system comprising: the system comprises a multi-IMU fusion module and a position and speed updating module;
a multi-IMU fusion module: the virtual fusion IMU is obtained according to the measurement data of the IMU;
a location speed update module: and the IMU is used for positioning according to the virtual fusion IMU.
2. The multi-IMU fused localization system of claim 1, wherein the multi-IMU fusion module comprises: a plurality of the IMUs and a multi-IMU integrated processing unit.
3. The multi-IMU fused localization system of claim 1, wherein the localization process comprises: location update and velocity update.
4. The multi-IMU fused positioning system of claim 2, wherein a plurality of said IMUs are disposed on a device in a predetermined manner.
5. The multi-IMU fused positioning system of claim 4, wherein the apparatus is provided on a plurality of locations on a worker.
6. A multi-IMU fused localization method adapted to the multi-IMU fused localization system of claim 1, the method comprising:
step 1: acquiring measurement data of each IMU;
step 2: obtaining a relative relation matrix according to the relative position relation of the IMUs;
and 3, step 3: obtaining a virtual fusion IMU according to the measurement data and the relative relation matrix;
and 4, step 4: and performing positioning processing according to the virtual fusion IMU.
7. The multi-IMU fused localization method of claim 6, wherein the measurement data comprises angular velocity data and acceleration data;
the step 1 comprises the following steps:
step 101: obtaining the angular velocity data according to the gyroscope of each IMU;
step 102: and obtaining the acceleration data according to the accelerometer of each IMU.
8. The method for multi-IMU fused localization according to claim 6 or 7, wherein the step 4 comprises:
step 401: obtaining a virtual fusion gyroscope model according to the virtual fusion IMU;
step 402: obtaining a virtual fusion accelerometer model according to the virtual fusion IMU;
step 403: and positioning according to the virtual fusion gyroscope model and the virtual fusion accelerometer model.
9. The multi-IMU fused localization method of claim 6, wherein the virtual fused IMU is optimized by a least squares method.
10. The multi-IMU fused localization method of claim 6, wherein the localization process comprises: location update and velocity update.
CN202210996020.3A 2022-08-19 2022-08-19 Multi-IMU fusion positioning system and method Pending CN115060259A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210996020.3A CN115060259A (en) 2022-08-19 2022-08-19 Multi-IMU fusion positioning system and method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210996020.3A CN115060259A (en) 2022-08-19 2022-08-19 Multi-IMU fusion positioning system and method

Publications (1)

Publication Number Publication Date
CN115060259A true CN115060259A (en) 2022-09-16

Family

ID=83207560

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210996020.3A Pending CN115060259A (en) 2022-08-19 2022-08-19 Multi-IMU fusion positioning system and method

Country Status (1)

Country Link
CN (1) CN115060259A (en)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040036650A1 (en) * 2002-08-26 2004-02-26 Morgan Kenneth S. Remote velocity sensor slaved to an integrated GPS/INS
CN104677355A (en) * 2015-03-06 2015-06-03 九江飞恩微电子有限公司 Multi-sensor fusion based virtual gyroscope and method
CN106482734A (en) * 2016-09-28 2017-03-08 湖南优象科技有限公司 A kind of filtering method for IMU Fusion
US20190145776A1 (en) * 2017-11-15 2019-05-16 Baidu Online Network Technology (Beijing) Co., Ltd Integrated positioning method and system
CN109827593A (en) * 2018-09-11 2019-05-31 广东星舆科技有限公司 A kind of error self-calibrating method, system and storage medium based on more IMU
CN112611380A (en) * 2020-12-03 2021-04-06 燕山大学 Attitude detection method based on multi-IMU fusion and attitude detection device thereof
CN112815937A (en) * 2020-12-31 2021-05-18 中国人民解放军战略支援部队航天工程大学 Data fusion optimal weight estimation method for redundant inertial measurement unit
CN114526729A (en) * 2022-01-14 2022-05-24 重庆邮电大学 Course optimization method of MEMS inertial positioning system based on redundancy technology

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040036650A1 (en) * 2002-08-26 2004-02-26 Morgan Kenneth S. Remote velocity sensor slaved to an integrated GPS/INS
CN104677355A (en) * 2015-03-06 2015-06-03 九江飞恩微电子有限公司 Multi-sensor fusion based virtual gyroscope and method
CN106482734A (en) * 2016-09-28 2017-03-08 湖南优象科技有限公司 A kind of filtering method for IMU Fusion
US20190145776A1 (en) * 2017-11-15 2019-05-16 Baidu Online Network Technology (Beijing) Co., Ltd Integrated positioning method and system
CN109827593A (en) * 2018-09-11 2019-05-31 广东星舆科技有限公司 A kind of error self-calibrating method, system and storage medium based on more IMU
CN112611380A (en) * 2020-12-03 2021-04-06 燕山大学 Attitude detection method based on multi-IMU fusion and attitude detection device thereof
CN112815937A (en) * 2020-12-31 2021-05-18 中国人民解放军战略支援部队航天工程大学 Data fusion optimal weight estimation method for redundant inertial measurement unit
CN114526729A (en) * 2022-01-14 2022-05-24 重庆邮电大学 Course optimization method of MEMS inertial positioning system based on redundancy technology

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
MING ZHANG等: "A Lightweight and Accurate Localization Algorithm Using Multiple Inertial Measurement Units", 《IEEE ROBOTICS AND AUTOMATION LETTERS》 *
王玲玲等: "基于低成本IMU的捷联航姿系统软件设计与实现", 《中国惯性技术学报》 *

Similar Documents

Publication Publication Date Title
Simanek et al. Evaluation of the EKF-based estimation architectures for data fusion in mobile robots
US10247576B2 (en) Method and system for verifying measured data
JP5258362B2 (en) Fault detection, isolation, and reconfiguration of inertial measurement devices using parity logic
US8756995B2 (en) Device and method for combining samples from an inertial measurement sensor cluster
CA2510714C (en) System for navigation redundancy
Youn et al. Combined quaternion-based error state Kalman filtering and smooth variable structure filtering for robust attitude estimation
US8260477B2 (en) Method and apparatus for tracking center of gravity of air vehicle
CN105606096B (en) A kind of posture of carrier movement status information auxiliary and course calculate method and system
JP2010032398A (en) Location detecting apparatus and method of navigation system
CN107014376A (en) A kind of posture inclination angle method of estimation suitable for the accurate operation of agricultural machinery
CN102506875B (en) The air navigation aid of unmanned plane and device
CN110567457A (en) Inertial navigation self-detection system based on redundancy
Malleswaran et al. Integration of INS and GPS using radial basis function neural networks for vehicular navigation
KR101642286B1 (en) Heading Orientation Estimation Method Using Pedestrian Characteristics in Indoor Environment
TW202016506A (en) Satellite attitude data fusion system and method thereof applicable to the space satellite environment for estimating attitude data of the satellite
Wu et al. Wheel-INS2: multiple MEMS IMU-based dead reckoning system with different configurations for wheeled robots
CN111141286A (en) Unmanned aerial vehicle flight control multi-sensor attitude confidence resolving method
CN116608853B (en) Carrier dynamic posture estimation method, device and storage medium
Zhao et al. Learning to compensate for the drift and error of gyroscope in vehicle localization
CN115060259A (en) Multi-IMU fusion positioning system and method
Buhmann et al. A GPS aided full linear accelerometer based gyroscope-free navigation system
CN101949702A (en) Quick self-testing method for GNSS PVT quality by using MEMS accelerometer
Jao et al. Prio-imu: Prioritizable imu array for enhancing foot-mounted inertial navigation accuracy
Shu-zhi et al. Optimization design and calibration of installation error coefficients for gyroscope-free strapdown inertial measurement unit
CN113449248B (en) Data fusion method and device for integrated SINS/GNSS system

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20220916

RJ01 Rejection of invention patent application after publication