CN106908060A - A kind of high accuracy indoor orientation method based on MEMS inertial sensor - Google Patents
A kind of high accuracy indoor orientation method based on MEMS inertial sensor Download PDFInfo
- Publication number
- CN106908060A CN106908060A CN201710079910.7A CN201710079910A CN106908060A CN 106908060 A CN106908060 A CN 106908060A CN 201710079910 A CN201710079910 A CN 201710079910A CN 106908060 A CN106908060 A CN 106908060A
- Authority
- CN
- China
- Prior art keywords
- zero
- inertial sensor
- error
- mems inertial
- speed
- 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
- 238000000034 method Methods 0.000 title claims abstract description 43
- 230000033001 locomotion Effects 0.000 claims abstract description 31
- 238000001514 detection method Methods 0.000 claims abstract description 19
- 238000012937 correction Methods 0.000 claims abstract description 15
- 230000001133 acceleration Effects 0.000 claims description 23
- 238000001914 filtration Methods 0.000 claims description 13
- 238000002474 experimental method Methods 0.000 claims description 10
- 238000005259 measurement Methods 0.000 claims description 9
- 230000008859 change Effects 0.000 claims description 7
- 230000005540 biological transmission Effects 0.000 claims description 6
- 230000003044 adaptive effect Effects 0.000 claims description 4
- 238000013461 design Methods 0.000 claims description 4
- 230000009194 climbing Effects 0.000 claims description 2
- 238000003672 processing method Methods 0.000 claims description 2
- 230000008569 process Effects 0.000 description 11
- 230000008901 benefit Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000005021 gait Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/10—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration
- G01C21/12—Navigation; 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/16—Navigation; 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/20—Instruments for performing navigational calculations
- G01C21/206—Instruments for performing navigational calculations specially adapted for indoor navigation
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Automation & Control Theory (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Navigation (AREA)
Abstract
The invention discloses a kind of high accuracy indoor orientation method based on MEMS inertial sensor, comprise the following steps:(1) MEMS inertial sensor is connected firmly on pedestrian's pin, MEMS inertial sensor is sensed the motion state of foot, foot navigation information is obtained in real time, and realize transmitting by bluetooth;(2) the hand-held Android mobile phone of pedestrian, Android client real-time reception and preserve MEMS inertial sensor offer data;(3) denoising is carried out to data;(4) zero velocity interval is obtained using zero-speed detection algorithm first, error correction is secondly carried out using zero-velocity curve algorithm bonding state algorithm for estimating;(5) Android client is shown by the data after error correction and compensation in real time by user interface.The present invention does not need extra auxiliary foundation to set;Good positioning precision can be kept to various complicated motion states in indoor positioning;Corrected in real time using mobile terminal and compensate error, run trace is shown by User Interface in real time.
Description
Technical Field
The invention relates to an indoor positioning method, in particular to a high-precision indoor positioning method based on an MEMS inertial sensor without installing additional infrastructure.
Background
Lbs (location based services) is becoming more and more important due to the development of the internet, the popularity of mobile devices and personal devices, and users obtain location information and use it for services such as navigation, tracking, monitoring, information push, and the like. The GPS can conveniently provide outdoor personal positioning information, but because the GPS needs to receive at least 4 satellites to realize positioning, the positioning effect is greatly affected by the received satellite signals, and the positioning effect is seriously affected due to the fact that the GPS cannot receive the satellite signals under an indoor environment.
The inertial navigation system has the advantages of strong autonomy, high output frequency, high short-time precision and the like, and particularly, the MEMSIMU is rapidly developed in recent years, so that the inertial navigation system becomes small in size and low in cost. The advantage of autonomous navigation can be well played by using a mode of Pedestrian Dead Reckoning (PDR) for indoor positioning of the MEMS IMU. From the current indoor positioning research condition, the MEMS indoor positioning is a big research hotspot already abroad, the condition of error accumulation of the MEMS inertial sensor can be well inhibited by judging the motion state of a carrier firstly and then detecting a zero-speed interval through information output by the MEMS inertial sensor, and further performing zero-speed correction (zero velocity Update), meanwhile, the initial alignment can not be performed according to the high-precision inertial navigation system initial alignment method due to the low precision of the low-cost MEMS IMU gyro, and the magnetometer is stable compared with the gyro, so that the magnetometer is introduced to assist course alignment and provide course compensation.
In the process of indoor positioning of the MEMS inertial sensor, zero-speed correction is adopted to correct errors accumulated in a walking period, so that the process has important significance in the process of improving the positioning accuracy of the MEMS inertial sensor, the zero-speed correction accuracy is directly determined by detecting a zero-speed interval in a walking process, a threshold value is set for detection methods of data searched at home and abroad about the zero-speed interval, corresponding detection quantity is compared with the threshold value, the selection of the threshold value is frequently one-sided and can only aim at a certain specific motion state, indoor footsteps can have various complex motion states such as walking, running, going upstairs, going downstairs and the like, and therefore the detection accuracy needs to be improved urgently.
Disclosure of Invention
The purpose of the invention is as follows: in order to overcome the defects of the prior art, the invention aims to provide a high-precision indoor positioning method based on a MEMS inertial sensor, which is suitable for autonomous positioning in an indoor environment and does not need additional installation of infrastructure.
The technical scheme is as follows: in order to realize the purpose of the invention, the following technical scheme is adopted: based on MEMS inertial sensor
The high-precision indoor positioning method comprises the following steps:
(1) the MEMS inertial sensor is fixedly connected to the foot of the pedestrian, so that the MEMS inertial sensor senses the motion state of the foot and acquires navigation information of the foot in real time, and the navigation information is transmitted through the Bluetooth transmission module;
(2) the pedestrian holds the android mobile phone, receives data provided by the MEMS inertial sensor for navigating the foot motion state in real time at an android client, and stores the data;
(3) denoising the data;
(4) aiming at the problem of low positioning accuracy of the MEMS inertial sensor, firstly, a zero-speed detection algorithm is adopted to obtain a zero-speed interval, and secondly, a zero-speed correction algorithm is adopted to be combined with a state estimation algorithm to perform error correction so as to improve the positioning accuracy; the zero-speed interval detection method comprises the following steps of acquiring data transmitted by the MEMS inertial sensor through an android client, and detecting the zero-speed interval by adopting a sliding window method:
(a) according to the noise characteristics of an accelerometer and a gyroscope of the MEMS inertial sensor, repeated experiments are carried out in different motion states to obtain zero-speed interval thresholds in the corresponding motion states, wherein the different motion states comprise normal walking, running, climbing stairs or descending stairs;
(b) according to the acceleration information of the accelerometer acquired by the android client, firstly, the motion state of the footstep at the moment is determined by analyzing the characteristics of the acceleration in the sliding window and comparing the change of the acceleration at the adjacent moments and the conditions of the maximum value and the minimum value of the acceleration in the window, and the threshold value of the zero-speed detection corresponding to the motion state is automatically selected;
(c) respectively judging whether the specific force module value, the specific force variance and the angular velocity module value are within a threshold range in the sliding window interval, if the specific force module value, the specific force variance and the angular velocity module value simultaneously meet a threshold condition, judging that the interval is a zero-velocity interval, and after obtaining the zero-velocity interval, determining that the interval is an interval in which the foot is at zero velocity;
at the moment, zero-speed correction is adopted in combination with an adaptive Kalman filtering algorithm, and errors accumulated in a walking interval are eliminated.
(5) And the android client displays the data subjected to error correction and compensation in real time through the user interface display module, and presents the motion state to the user in real time.
The data provided by the MEMS inertial sensor comprises data measured by a three-axis accelerometer, a three-axis gyroscope and a three-axis magnetometer. The android client has the function of calibrating the data, and can record the data and process the data in real time.
The android client has the function of displaying the walking track of the user in real time.
In the step (3), aiming at the noise problem generated when the feet of the MEMS inertial sensor are severely shaken, the denoising processing method is wavelet denoising, and can remove noise and interference signals contained in high-frequency signals.
In the step (4), aiming at the problem that the error of the low-cost inertial sensor is easy to disperse along with the time, the state estimation algorithm is an adaptive Kalman filtering algorithm, and the step of correcting the error by adopting a zero-speed correction algorithm and combining the state estimation algorithm is as follows:
the accurate zero-speed interval can be obtained through the zero-speed detection algorithm, the speed is zero in the zero-speed interval process, and therefore the speed output by navigation calculation is the speed error. The relationship between the speed error and the gyro drift, the acceleration zero offset and the attitude angle error can be described by an error differential equation. A Kalman filtering equation is established on the basis of an error equation, the speed error is used as an observed quantity, and an attitude angle error, a position error and other error quantities are estimated and fed back to a navigation resolving module, so that the purpose of correcting the error of the system is achieved.
Error equation of attitude angle
Equation of speed error
Equation of position error
A Kalman filtering equation is established on the basis of an error equation, the speed error is used as an observed quantity, an attitude angle error, a position error and other differential quantities are estimated and fed back to a navigation resolving module,
selecting the state quantity of the Kalman filter as follows:
whereinIndicating attitude angle error, vnIndicating a speed error; p is a radical ofnIndicating a position error; the gyro is represented by a constant drift of the gyro,the system equation representing the accelerometer constant zero offset, Kalman is expressed as
X(t)=FX(t)+W(t)
Measuring velocity errors of equations, i.e.
Zk=[vn]=HkXk
The measurement can be obtained only in the zero-speed interval, the measurement is related to the specific motion and interference conditions, and after the zero-speed interval is detected by the zero-speed detection algorithm, the Kalman filtering updates the time and the measurement, and feeds the estimated error back to the system for error compensation.
The MEMS inertial sensor and the Bluetooth module adopt a patch type hardware design method, and have the advantage of small volume.
Has the advantages that: according to the human body walking dynamics model, the MEMS inertial sensor with low cost is adopted, the system has the advantages of small volume and light weight, and compared with a wireless positioning mode, the system structure can realize autonomous positioning without additional auxiliary foundation arrangement; the invention adopts a self-adaptive threshold matching method, can keep good positioning precision for various complicated motion states of footsteps in the indoor positioning process, and has stronger popularity; aiming at the current popular mobile terminal, the mobile terminal is designed to receive the measurement data of the inertial sensor transmitted by the Bluetooth in real time, has the functions of storing and processing the data, corrects and compensates errors of the data in real time, and simultaneously designs a user interaction interface to display the walking track of a user in real time.
Drawings
FIG. 1 is a diagram of a MEMS inertial sensor indoor positioning system;
FIG. 2 is a graph of a MEMS inertial sensor before denoising and after denoising smoothing;
FIG. 3 is a flow chart of a zero-velocity interval detection method;
FIG. 4 is a diagram of the detection effect of the zero-speed interval;
FIG. 5 is an android client interface;
FIG. 6 is a two-dimensional planar walking trajectory diagram of the present invention;
fig. 7 is a three-dimensional up and down stair walking track diagram of the invention.
Detailed Description
The technical solution of the present invention will be further described in detail with reference to the following examples and accompanying drawings.
As shown in fig. 1, the high-precision indoor positioning apparatus employed in the present invention includes: low-cost MEMS inertial sensors such as MPU9260, MPU6050 and the like, a Bluetooth transmission module and an android client; the android client comprises a data calibration module, a data storage module, a data processing module and a user interface display module. The MEMS inertial sensor is fixedly connected to the foot, the change of the step state is sensed, the measured parameters are transmitted through the Bluetooth transmission module, and the data transmitted by the Bluetooth transmission module are received through the android client. The android client firstly calibrates the three-axis accelerometer, the three-axis gyroscope and the three-axis magnetometer respectively through the data calibration module, stores and processes data in real time, and simultaneously displays the data on a user interface in real time. The MEMS inertial sensor and the Bluetooth module both adopt a surface mount hardware design method, and have the advantages of compact structure and small volume.
The method comprises the following concrete steps:
(1) the low-cost inertial sensor is fixedly connected to the foot, so that the inertial sensor can sense the motion state of the foot, the MEMS inertial sensor can normally output information of the three-axis accelerometer, the three-axis gyroscope and the three-axis magnetometer, navigation information of the foot can be acquired in real time, and the navigation information is transmitted through the Bluetooth transmission module.
(2) The android mobile phone is held by a pedestrian, data transmitted by a foot inertial sensor is received at an android client in real time, a data calibration module of the android client is firstly opened to respectively calibrate a three-axis accelerometer, a three-axis gyroscope and measurement data of the three-axis magnetometer are calibrated, the accelerometer is removed for the main purpose of calibration, the gyroscope and the zero offset of the magnetometer, the calibration method of the accelerometer and the gyroscope is to level the MEMS inertial sensor, when the zero offset displayed by the calibration module does not change any more, the zero offset at the moment is stored, namely the calibration of the accelerometer and the gyroscope is completed, the MEMS inertial sensor is repeatedly rotated around different directions for the calibration of a magnetic field, and the calibration of the magnetic field is completed after the zero offset of the magnetic field is displayed to change no more by multiple movements. And after the calibration of the data is finished, whether the accelerometer and the gyroscope are zero or not is judged when the MEMS inertial sensor is in a static state, and the data is stored.
(3) The wavelet denoising processing is performed on the obtained data, as shown in fig. 2, which is a graph comparing the data measured by the MEMS inertial sensor and the curve of the data after wavelet denoising.
(4) The data of the MEMS inertial sensor is detected in a zero-speed interval by adopting a sliding window method, and the detection method of the zero-speed interval is shown in figure 3.
The zero-speed detection algorithm plays an important role in the whole MEMS positioning process and directly determines the positioning precision. The identification method comprises the following steps:
calculating the triaxial acceleration module value at each sampling moment i:
selecting a sliding window with the size of 2N +1, wherein N represents the window width, and the acceleration variance in the window is as follows:
wherein,
scanning acceleration data at regular intervals, identifying the state that the pace in the period of time is normal walking, slow speed, fast walking or going upstairs and downstairs and the like through the difference between the maximum value and the minimum value of the scanned acceleration, and determining the acceleration threshold values of the normal walking, the slow speed, the fast walking or the going upstairs and downstairs by multiple tests to respectively take tha1、tha2、tha3、tha4(ii) a Threshold angular velocity is taken thw1、thw2、thw3、thw4(ii) a Covariance th of accelerationσ1、thσ2、thσ3、thσ4。
Above mainly because the motion state identification method of acceleration modulus value also can discern based on angular velocity or acceleration covariance method, after discerning the motion state, carries out zero velocity interval and detects, in order to improve the precision that detects, according to the motion state of discerning above, relies on single element to detect must have the error, and this embodiment adopts the method that the three combines to detect:
the first condition is as follows: modulus of proportionalityWhen the feet are at the zero-speed moment, the output result of the acceleration measured by the MEMS inertial sensor fixed on the feet theoretically has only the local gravity acceleration value, the acceleration accords with chi-square distribution, the judgment can be carried out only by setting a proper confidence interval, and the assumption is that the judgment can be carried out on the condition that the acceleration accords with chi-square distributionThe minimum value and the maximum value of the confidence interval are respectively;
where 1 indicates a stationary state is detected and 0 indicates a moving state.
And a second condition: variance of specific force. The sudden change stage of the signal in the gait cycle can be effectively detected by observing the variance of the specific force through the sliding window. The variance variation of the proportion in the zero-velocity interval is weak and hardly fluctuates, and the step variation in the non-zero-velocity interval has strong fluctuation, so that the zero-velocity interval can be judged by the specific force variance being smaller than a given threshold value. The calculation is as follows:
wherein,is a specific force average value; m is the size of the sliding window and is related to the output frequency;
where 1 indicates a stationary state is detected and 0 indicates a moving state.
And (3) carrying out a third condition: the angular velocity modulus. Because the sole part keeps a complete contact state with the ground in the zero-speed interval, the angular speed and the change value thereof tend to zero in the zero-speed interval, and the angular speed has stronger fluctuation at the time outside the zero-speed interval, so that the corresponding angular speed threshold th can be determined through the identified walking statewTo determine a zero-speed interval in which the angular velocity module value is
Where 1 indicates a stationary state is detected and 0 indicates a moving state.
And when the condition I, the condition II and the condition III are simultaneously satisfied, recording that the state is 1 at the moment, the zero speed interval is required to be detected, and if the state is 0, the motion state is indicated.
Fig. 4 shows the zero velocity interval detected by the above-mentioned zero velocity interval detection method.
When the zero-speed interval is detected by the method, the error accumulated by the MEMS inertial sensor in a walking period is corrected in the zero-speed interval so as to compensate the whole positioning system.
The accurate zero-speed interval can be obtained through the zero-speed detection algorithm, the speed is zero in the zero-speed interval process, and therefore the speed output by navigation calculation is the speed error. The relationship between the speed error and the gyro drift, the acceleration zero offset and the attitude angle error can be described by an error differential equation. A Kalman filtering equation can be established on the basis of an error equation, the speed error is used as an observed quantity, an attitude angle error is estimated, and a position error and other error quantities are fed back to a navigation resolving module, so that the purpose of correcting the error of the system is achieved.
(a) Equation of attitude error
(b) Error in velocity
(c) Equation of position error
Selecting the state quantity of the Kalman filter as follows:
whereinIndicating the attitude error angle, vnIndicating a speed error; p is a radical ofnIndicating a position error; the gyro is represented by a constant drift of the gyro,representing that the accelerometer has a constant value of zero offset, theoretically, the influence factors of the states should be considered as many as possible, the more the factors are considered, the higher the accuracy of the navigation system estimation is, but the order of the system model is increased[11]The system of Kalman equations may be expressed as
X(t)=FX(t)+W(t)
Measuring velocity errors of equations, i.e.
Zk=[vn]=HkXk
Measurement of v quantitykThe Kalman filtering method is characterized in that the Kalman filtering method only can be obtained in a zero-speed interval, the measurement is related to specific motion and interference conditions, after the zero-speed interval is detected through the zero-speed detection algorithm, time updating and measurement updating are carried out through Kalman filtering, estimated errors are fed back to a system, and error compensation is carried out
(5) And displaying the corrected data in real time on a user interface display module, so that a user can clearly see the motion state of the user.
The invention relates to an independent autonomous high-precision indoor positioning method, and the beneficial effects of the invention are verified through partial experiments.
In order to verify the effectiveness of the method, a normal two-dimensional plane walking experiment and a three-dimensional stair-climbing walking experiment are respectively carried out by using a low-cost MEMS MPU9250, data transmitted by an MEMS sensor through Bluetooth are received in real time by a handheld android mobile phone, and an android client interface of the mobile phone is shown in fig. 5.
Experiment one: two-dimensional plane normal walking experiment
In the experimental environment, a corridor of a student building is selected to be reserved in the south-east university, wherein the corridor is in a rectangular shape with the length of 5m and the width of 1 m.
Fig. 6 is a walking trajectory diagram processed by the algorithm of the present invention for data received by an android client, where the walking trajectory is a closed rectangular curve under an ideal situation without error, and a start point and an end point are coincident, but due to the error, the start point and the end point are often misaligned. The positioning accuracy is judged according to the misalignment of the two, and the positioning accuracy of the two-dimensional walking experiment is within 0.5 m.
Experiment two: three-dimensional stair-climbing walking experiment
The experimental environment is selected in the central east-south university building, the starting point is the 2 nd th floor of the central building, the walking route is the position from the right side stair on the 2 nd floor of the central building to the 4 th floor of the central building through the 3 rd floor, and then the walking route goes from the left side of the corridor on the 4 th floor to the left side stair on the central building to the 2 nd floor through the 3 rd floor and returns to the starting point of the 2 nd floor, and the walking route also runs in a closed route. The walking process goes through going upstairs, normal walking, going downstairs and normal walking.
FIG. 7 is a walking trajectory graph of data obtained by an android client processed by the algorithm of the present invention, and the positioning error is calculated to be within 2 m.
Through experimental analysis, the invention can keep higher positioning precision for normal two-dimensional plane walking and still has higher positioning precision for complex upstairs and downstairs movement states. The corrected data are displayed to the user through a user interface display module of the android client, and autonomous navigation of the user can be achieved in various unknown environments.
Claims (6)
1. A high-precision indoor positioning method based on an MEMS inertial sensor is characterized by comprising the following steps:
(1) the MEMS inertial sensor is fixedly connected to the foot of the pedestrian, so that the MEMS inertial sensor senses the motion state of the foot and acquires navigation information of the foot in real time, and the navigation information is transmitted through the Bluetooth transmission module;
(2) the pedestrian holds the android mobile phone, receives data provided by the MEMS inertial sensor for navigating the foot motion state in real time at an android client, and stores the data;
(3) denoising the data;
(4) firstly, obtaining a zero-speed interval by adopting a zero-speed detection algorithm, and secondly, carrying out error correction by adopting a zero-speed correction algorithm in combination with a state estimation algorithm; the zero-speed interval detection method comprises the following steps of acquiring data transmitted by the MEMS inertial sensor through an android client, and detecting the zero-speed interval by adopting a sliding window method:
(a) according to the noise characteristics of an accelerometer and a gyroscope of the MEMS inertial sensor, repeated experiments are carried out in different motion states to obtain zero-speed interval thresholds in the corresponding motion states, wherein the different motion states comprise normal walking, running, climbing stairs or descending stairs;
(b) according to the acceleration information of the accelerometer acquired by the android client, firstly, the motion state of the footstep at the moment is determined by analyzing the characteristics of the acceleration in the sliding window and comparing the change of the acceleration at the adjacent moments and the conditions of the maximum value and the minimum value of the acceleration in the window, and the threshold value of the zero-speed detection corresponding to the motion state is automatically selected;
(c) respectively judging whether the specific force module value, the specific force variance and the angular velocity module value are within a threshold range in the sliding window interval, if the specific force module value, the specific force variance and the angular velocity module value simultaneously meet a threshold condition, judging that the interval is a zero-velocity interval, and after obtaining the zero-velocity interval, determining that the interval is an interval in which the foot is at zero velocity;
(5) and the android client displays the data subjected to error correction and compensation in real time through the user interface display module.
2. The high-precision indoor positioning method based on the MEMS inertial sensor of claim 1, wherein: the data provided by the MEMS inertial sensor comprises data measured by a three-axis accelerometer, a three-axis gyroscope and a three-axis magnetometer.
3. The high-precision indoor positioning method based on the MEMS inertial sensor of claim 1, wherein: the android client has the function of displaying the walking track of the user in real time.
4. The high-precision indoor positioning method based on the MEMS inertial sensor of claim 1, wherein: in the step (3), the denoising processing method is wavelet denoising.
5. The high-precision indoor positioning method based on the MEMS inertial sensor of claim 1, wherein: in the step (4), the state estimation algorithm is an adaptive kalman filtering algorithm, and the step of performing error correction by adopting a zero-speed correction algorithm and combining the adaptive kalman filtering algorithm is as follows:
and (3) describing the relationship between the speed error and the gyro drift, the acceleration zero offset and the attitude angle error by using an error differential equation:
error equation of attitude angle
Equation of speed error
Equation of position error
Establishing a Kalman filtering equation on the basis of an error equation, taking a speed error as an observed quantity, estimating an attitude angle error, a position error, a gyro drift and an acceleration zero offset, feeding back to a navigation resolving module,
selecting the state quantity of the Kalman filter as follows:
wherein the attitude angle error, v, is expressednIndicating a speed error; p is a radical ofnIndicating a position error; the gyro is represented by a constant drift of the gyro,the system equation representing the accelerometer constant zero offset, Kalman is expressed as
X(t)=FX(t)+W(t)
Measuring velocity errors of equations, i.e.
Zk=[vn]=HkXk
And after the zero-speed interval is detected, the Kalman filtering updates time and measurement, and feeds the estimated error back to the system for error compensation.
6. The high-precision indoor positioning method based on the MEMS inertial sensor of claim 1, wherein: the MEMS inertial sensor and the Bluetooth module adopt a patch type hardware design method.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710079910.7A CN106908060A (en) | 2017-02-15 | 2017-02-15 | A kind of high accuracy indoor orientation method based on MEMS inertial sensor |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710079910.7A CN106908060A (en) | 2017-02-15 | 2017-02-15 | A kind of high accuracy indoor orientation method based on MEMS inertial sensor |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106908060A true CN106908060A (en) | 2017-06-30 |
Family
ID=59207504
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710079910.7A Pending CN106908060A (en) | 2017-02-15 | 2017-02-15 | A kind of high accuracy indoor orientation method based on MEMS inertial sensor |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106908060A (en) |
Cited By (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107490378A (en) * | 2017-07-17 | 2017-12-19 | 北京工业大学 | A kind of indoor positioning based on MPU6050 and smart mobile phone and the method for navigation |
CN107843256A (en) * | 2017-10-11 | 2018-03-27 | 南京航空航天大学 | Adaptive zero-velocity curve pedestrian navigation method based on MEMS sensor |
CN107990895A (en) * | 2017-11-08 | 2018-05-04 | 北京工商大学 | A kind of building floor gap pedestrian track tracking and system based on wearable IMU |
CN108362282A (en) * | 2018-01-29 | 2018-08-03 | 哈尔滨工程大学 | A kind of inertia pedestrian's localization method based on the adjustment of adaptive zero-speed section |
CN108507572A (en) * | 2018-05-28 | 2018-09-07 | 清华大学 | A kind of attitude orientation error correcting method based on MEMS Inertial Measurement Units |
CN108593965A (en) * | 2018-05-02 | 2018-09-28 | 福州大学 | A kind of accelerometer moored condition scaling method based on specific force mould and stable inertia |
CN108600033A (en) * | 2018-05-15 | 2018-09-28 | 哈尔滨理工大学 | A kind of wireless sensor network universal nodes and method for diagnosing faults |
CN108680184A (en) * | 2018-04-19 | 2018-10-19 | 东南大学 | A kind of zero-speed detection method based on Generalized Likelihood Ratio statistic curve geometric transformation |
CN109186603A (en) * | 2018-08-16 | 2019-01-11 | 浙江树人学院 | 3-D positioning method in a kind of fireman room based on multisensor |
CN109361767A (en) * | 2018-12-06 | 2019-02-19 | 苏州思必驰信息科技有限公司 | Optimize server-side, client process method and the server of client data display error, be able to carry out the client that data are shown |
CN109489694A (en) * | 2019-01-02 | 2019-03-19 | 中国船舶重工集团公司第七0七研究所 | A kind of individual soldier's navigation system zero-speed detection method of voltage sensitive sensor auxiliary |
CN109579838A (en) * | 2019-01-14 | 2019-04-05 | 湖南海迅自动化技术有限公司 | The localization method and positioning system of AGV trolley |
CN109579832A (en) * | 2018-11-26 | 2019-04-05 | 重庆邮电大学 | A kind of personnel's height autonomous positioning algorithm |
WO2019084918A1 (en) * | 2017-11-03 | 2019-05-09 | Beijing Didi Infinity Technology And Development Co., Ltd. | System and method for determining a trajectory |
CN110121149A (en) * | 2019-04-22 | 2019-08-13 | 西安邮电大学 | A kind of indoor orientation method based on bluetooth and PDR data fusion |
CN110133692A (en) * | 2019-04-18 | 2019-08-16 | 武汉苍穹电子仪器有限公司 | The high-precision GNSS dynamic tilt measuring system and method for inertial navigation technique auxiliary |
CN110207704A (en) * | 2019-05-21 | 2019-09-06 | 南京航空航天大学 | A kind of pedestrian navigation method based on the identification of architectural stair scene intelligent |
CN110274592A (en) * | 2019-07-18 | 2019-09-24 | 北京航空航天大学 | A kind of zero-speed section of waist foot Inertial Measurement Unit information fusion determines method |
CN110455284A (en) * | 2019-07-03 | 2019-11-15 | 中国人民解放军战略支援部队信息工程大学 | A kind of pedestrian movement patterns' recognition methods and device based on MEMS-IMU |
CN110657807A (en) * | 2019-09-30 | 2020-01-07 | 西安电子科技大学 | Indoor positioning displacement measurement method for detecting discontinuity based on wavelet transformation |
CN111929689A (en) * | 2020-07-22 | 2020-11-13 | 杭州电子科技大学 | Object imaging method based on sensor of mobile phone |
CN112435423A (en) * | 2019-08-22 | 2021-03-02 | 杭州海康威视数字技术股份有限公司 | Monitoring method and device |
CN113092819A (en) * | 2021-04-14 | 2021-07-09 | 东方红卫星移动通信有限公司 | Dynamic zero-speed calibration method and system for foot accelerometer |
CN113111480A (en) * | 2021-02-22 | 2021-07-13 | 同济大学 | Method and device for diagnosing and detecting running state of drainage pipe network |
CN113295158A (en) * | 2021-05-14 | 2021-08-24 | 江苏大学 | Indoor positioning method fusing inertial data, map information and pedestrian motion state |
CN113776525A (en) * | 2021-09-01 | 2021-12-10 | 东南大学 | Inertia/single sound source passive combination navigation method based on slope distance difference matching |
CN114440867A (en) * | 2021-12-17 | 2022-05-06 | 际络科技(上海)有限公司 | Zero-speed detection method and device for heavy truck |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040107072A1 (en) * | 2002-12-03 | 2004-06-03 | Arne Dietrich | Ins-based user orientation and navigation |
CN103616030A (en) * | 2013-11-15 | 2014-03-05 | 哈尔滨工程大学 | Autonomous navigation system positioning method based on strapdown inertial navigation resolving and zero-speed correction |
CN104296750A (en) * | 2014-06-27 | 2015-01-21 | 大连理工大学 | Zero speed detecting method, zero speed detecting device, and pedestrian navigation method as well as pedestrian navigation system |
CN104713554A (en) * | 2015-02-01 | 2015-06-17 | 北京工业大学 | Indoor positioning method based on MEMS insert device and android smart mobile phone fusion |
CN104931049A (en) * | 2015-06-05 | 2015-09-23 | 北京信息科技大学 | Movement classification-based pedestrian self-positioning method |
CN105043385A (en) * | 2015-06-05 | 2015-11-11 | 北京信息科技大学 | Self-adaption Kalman filtering method for autonomous navigation positioning of pedestrians |
-
2017
- 2017-02-15 CN CN201710079910.7A patent/CN106908060A/en active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040107072A1 (en) * | 2002-12-03 | 2004-06-03 | Arne Dietrich | Ins-based user orientation and navigation |
CN103616030A (en) * | 2013-11-15 | 2014-03-05 | 哈尔滨工程大学 | Autonomous navigation system positioning method based on strapdown inertial navigation resolving and zero-speed correction |
CN104296750A (en) * | 2014-06-27 | 2015-01-21 | 大连理工大学 | Zero speed detecting method, zero speed detecting device, and pedestrian navigation method as well as pedestrian navigation system |
CN104713554A (en) * | 2015-02-01 | 2015-06-17 | 北京工业大学 | Indoor positioning method based on MEMS insert device and android smart mobile phone fusion |
CN104931049A (en) * | 2015-06-05 | 2015-09-23 | 北京信息科技大学 | Movement classification-based pedestrian self-positioning method |
CN105043385A (en) * | 2015-06-05 | 2015-11-11 | 北京信息科技大学 | Self-adaption Kalman filtering method for autonomous navigation positioning of pedestrians |
Non-Patent Citations (2)
Title |
---|
张嗣瀛,王福利主编: "《2006中国控制与决策学术年会论文集》", 31 May 2006, 东北大学出版社 * |
张晓东: "基于MEMS惯性器件的个人导航系统研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (45)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107490378A (en) * | 2017-07-17 | 2017-12-19 | 北京工业大学 | A kind of indoor positioning based on MPU6050 and smart mobile phone and the method for navigation |
CN107490378B (en) * | 2017-07-17 | 2020-05-22 | 北京工业大学 | Indoor positioning and navigation method based on MPU6050 and smart phone |
CN107843256A (en) * | 2017-10-11 | 2018-03-27 | 南京航空航天大学 | Adaptive zero-velocity curve pedestrian navigation method based on MEMS sensor |
US11692829B2 (en) | 2017-11-03 | 2023-07-04 | Beijing Didi Infinity Technology And Development Co., Ltd. | System and method for determining a trajectory of a subject using motion data |
WO2019084918A1 (en) * | 2017-11-03 | 2019-05-09 | Beijing Didi Infinity Technology And Development Co., Ltd. | System and method for determining a trajectory |
TWI706295B (en) * | 2017-11-03 | 2020-10-01 | 大陸商北京嘀嘀無限科技發展有限公司 | System and method for determining a trajectory |
CN107990895A (en) * | 2017-11-08 | 2018-05-04 | 北京工商大学 | A kind of building floor gap pedestrian track tracking and system based on wearable IMU |
CN107990895B (en) * | 2017-11-08 | 2020-12-18 | 北京工商大学 | Building inter-floor pedestrian trajectory tracking method and system based on wearable IMU |
CN108362282A (en) * | 2018-01-29 | 2018-08-03 | 哈尔滨工程大学 | A kind of inertia pedestrian's localization method based on the adjustment of adaptive zero-speed section |
CN108680184A (en) * | 2018-04-19 | 2018-10-19 | 东南大学 | A kind of zero-speed detection method based on Generalized Likelihood Ratio statistic curve geometric transformation |
CN108680184B (en) * | 2018-04-19 | 2021-09-07 | 东南大学 | Zero-speed detection method based on generalized likelihood ratio statistical curve geometric transformation |
CN108593965A (en) * | 2018-05-02 | 2018-09-28 | 福州大学 | A kind of accelerometer moored condition scaling method based on specific force mould and stable inertia |
CN108600033A (en) * | 2018-05-15 | 2018-09-28 | 哈尔滨理工大学 | A kind of wireless sensor network universal nodes and method for diagnosing faults |
CN108507572B (en) * | 2018-05-28 | 2021-06-08 | 清华大学 | Attitude positioning error correction method based on MEMS inertial measurement unit |
CN108507572A (en) * | 2018-05-28 | 2018-09-07 | 清华大学 | A kind of attitude orientation error correcting method based on MEMS Inertial Measurement Units |
CN109186603A (en) * | 2018-08-16 | 2019-01-11 | 浙江树人学院 | 3-D positioning method in a kind of fireman room based on multisensor |
CN109186603B (en) * | 2018-08-16 | 2021-07-30 | 浙江树人学院 | Multi-sensor-based firefighter indoor three-dimensional positioning method |
CN109579832B (en) * | 2018-11-26 | 2022-12-27 | 重庆邮电大学 | Personnel height autonomous positioning algorithm |
CN109579832A (en) * | 2018-11-26 | 2019-04-05 | 重庆邮电大学 | A kind of personnel's height autonomous positioning algorithm |
CN109361767B (en) * | 2018-12-06 | 2021-11-02 | 思必驰科技股份有限公司 | Processing method for optimizing client data display error, server and client |
CN109361767A (en) * | 2018-12-06 | 2019-02-19 | 苏州思必驰信息科技有限公司 | Optimize server-side, client process method and the server of client data display error, be able to carry out the client that data are shown |
CN109489694A (en) * | 2019-01-02 | 2019-03-19 | 中国船舶重工集团公司第七0七研究所 | A kind of individual soldier's navigation system zero-speed detection method of voltage sensitive sensor auxiliary |
CN109579838A (en) * | 2019-01-14 | 2019-04-05 | 湖南海迅自动化技术有限公司 | The localization method and positioning system of AGV trolley |
CN110133692A (en) * | 2019-04-18 | 2019-08-16 | 武汉苍穹电子仪器有限公司 | The high-precision GNSS dynamic tilt measuring system and method for inertial navigation technique auxiliary |
CN110133692B (en) * | 2019-04-18 | 2023-01-31 | 武汉苍穹融新科技有限公司 | Inertial navigation technology-assisted high-precision GNSS dynamic inclination measurement system and method |
CN110121149A (en) * | 2019-04-22 | 2019-08-13 | 西安邮电大学 | A kind of indoor orientation method based on bluetooth and PDR data fusion |
CN110207704A (en) * | 2019-05-21 | 2019-09-06 | 南京航空航天大学 | A kind of pedestrian navigation method based on the identification of architectural stair scene intelligent |
CN110207704B (en) * | 2019-05-21 | 2021-07-13 | 南京航空航天大学 | Pedestrian navigation method based on intelligent identification of building stair scene |
CN110455284A (en) * | 2019-07-03 | 2019-11-15 | 中国人民解放军战略支援部队信息工程大学 | A kind of pedestrian movement patterns' recognition methods and device based on MEMS-IMU |
CN110274592B (en) * | 2019-07-18 | 2021-07-27 | 北京航空航天大学 | Zero-speed interval determination method for information fusion of waist and foot inertia measurement units |
CN110274592A (en) * | 2019-07-18 | 2019-09-24 | 北京航空航天大学 | A kind of zero-speed section of waist foot Inertial Measurement Unit information fusion determines method |
CN112435423A (en) * | 2019-08-22 | 2021-03-02 | 杭州海康威视数字技术股份有限公司 | Monitoring method and device |
CN110657807B (en) * | 2019-09-30 | 2022-09-23 | 西安电子科技大学 | Indoor positioning displacement measurement method for detecting discontinuity based on wavelet transformation |
CN110657807A (en) * | 2019-09-30 | 2020-01-07 | 西安电子科技大学 | Indoor positioning displacement measurement method for detecting discontinuity based on wavelet transformation |
CN111929689A (en) * | 2020-07-22 | 2020-11-13 | 杭州电子科技大学 | Object imaging method based on sensor of mobile phone |
CN111929689B (en) * | 2020-07-22 | 2023-04-07 | 杭州电子科技大学 | Object imaging method based on sensor of mobile phone |
CN113111480B (en) * | 2021-02-22 | 2022-07-29 | 同济大学 | Method and device for diagnosing and detecting running state of drainage pipe network |
CN113111480A (en) * | 2021-02-22 | 2021-07-13 | 同济大学 | Method and device for diagnosing and detecting running state of drainage pipe network |
CN113092819B (en) * | 2021-04-14 | 2022-11-18 | 东方红卫星移动通信有限公司 | Dynamic zero-speed calibration method and system for foot accelerometer |
CN113092819A (en) * | 2021-04-14 | 2021-07-09 | 东方红卫星移动通信有限公司 | Dynamic zero-speed calibration method and system for foot accelerometer |
CN113295158A (en) * | 2021-05-14 | 2021-08-24 | 江苏大学 | Indoor positioning method fusing inertial data, map information and pedestrian motion state |
CN113295158B (en) * | 2021-05-14 | 2024-05-14 | 江苏大学 | Indoor positioning method integrating inertial data, map information and pedestrian motion state |
CN113776525A (en) * | 2021-09-01 | 2021-12-10 | 东南大学 | Inertia/single sound source passive combination navigation method based on slope distance difference matching |
CN113776525B (en) * | 2021-09-01 | 2023-12-05 | 东南大学 | Inertia/single sound source passive combined navigation method based on oblique distance difference matching |
CN114440867A (en) * | 2021-12-17 | 2022-05-06 | 际络科技(上海)有限公司 | Zero-speed detection method and device for heavy truck |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106908060A (en) | A kind of high accuracy indoor orientation method based on MEMS inertial sensor | |
CN109827577B (en) | High-precision inertial navigation positioning algorithm based on motion state detection | |
Bai et al. | A high-precision and low-cost IMU-based indoor pedestrian positioning technique | |
CN104296750B (en) | Zero speed detecting method, zero speed detecting device, and pedestrian navigation method as well as pedestrian navigation system | |
CN110553646B (en) | Pedestrian navigation method based on inertia, magnetic heading and zero-speed correction | |
EP2951529B1 (en) | Inertial device, method, and program | |
CN111024126B (en) | Self-adaptive zero-speed correction method in pedestrian navigation positioning | |
AU781848B2 (en) | Pedestrian navigation method and apparatus operative in a dead reckoning mode | |
US9752879B2 (en) | System and method for estimating heading misalignment | |
JP7023234B2 (en) | How to estimate pedestrian movement | |
US10267646B2 (en) | Method and system for varying step length estimation using nonlinear system identification | |
US20110106487A1 (en) | Moving body positioning device | |
JP5586994B2 (en) | POSITIONING DEVICE, POSITIONING METHOD OF POSITIONING DEVICE, AND POSITIONING PROGRAM | |
US20130338961A1 (en) | Method and system for estimating a path of a mobile element or body | |
KR101394984B1 (en) | In-door positioning apparatus and method based on inertial sensor | |
CN101151508A (en) | Traveling direction measuring apparatus and traveling direction measuring method | |
JP2010210634A (en) | Method for calculating trajectory of geographical track | |
CN110398245A (en) | The indoor pedestrian navigation Attitude estimation method of formula Inertial Measurement Unit is worn based on foot | |
CN111854752B (en) | Dead reckoning by determining a misalignment angle between a direction of movement and a direction of sensor travel | |
CN106403952A (en) | Method for measuring combined attitudes of Satcom on the move with low cost | |
CN107014388A (en) | A kind of pedestrian track projectional technique and device detected based on magnetic disturbance | |
KR20150106004A (en) | Method and apparatus for handling vertical orientations of devices for constraint free portable navigation | |
JP2002139340A (en) | Walking navigation device and navigation system using the same | |
US10466054B2 (en) | Method and system for estimating relative angle between headings | |
US10914793B2 (en) | Method and system for magnetometer calibration |
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: 20170630 |
|
RJ01 | Rejection of invention patent application after publication |