CN108120452A - The filtering method of MEMS gyroscope dynamic data - Google Patents
The filtering method of MEMS gyroscope dynamic data Download PDFInfo
- Publication number
- CN108120452A CN108120452A CN201711396612.7A CN201711396612A CN108120452A CN 108120452 A CN108120452 A CN 108120452A CN 201711396612 A CN201711396612 A CN 201711396612A CN 108120452 A CN108120452 A CN 108120452A
- Authority
- CN
- China
- Prior art keywords
- mtd
- mtr
- msub
- mtable
- gyroscope
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
- G01C25/005—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass initial alignment, calibration or starting-up of inertial devices
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
- G06F2218/04—Denoising
Landscapes
- Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Gyroscopes (AREA)
Abstract
A kind of filtering method of MEMS gyroscope dynamic data, belongs to field of signal processing.The purpose of the present invention is being handled mainly for the dynamical output data of inexpensive MEMS gyroscope system, so as to improve the filtering method of the MEMS gyroscope dynamic data of gyroscope dynamical output data precision.The method of the present invention for improving MEMS gyroscope dynamical output data precision comprises the following steps:It determines MEMS gyroscope constant value drift, establishes MEMS gyroscope output data model, determines Kalman filter process-noise varianceWith measurement noise variance
Description
Technical field
The invention belongs to field of signal processing.
Background technology
MEMS gyroscope is for measuring the sensor of object angular speed, due to small, light-weight, at low cost, reliable
Property it is high the advantages that, be widely used in the various fields such as military, business, civilian.But due to system structure deviation, perturbed force
How the influence of square, ambient noise etc. there are random error in the output data of MEMS gyroscope, therefore improves MEMS gyroscope
Precision be the hot issue studied both at home and abroad at present.For this problem, there are two types of main solutions, and one kind is from hard
Improve gyroscope arrangement in terms of part, gyroscope precision is improved by reducing structural failure;Another kind is set about from software aspects, is led to
It crosses filtering algorithm elimination random error and reaches raising gyroscope precision purpose.Since the construction cycle is short, cost is few, current second
The research work of method is relatively more, and conventional method is first to establish gyroscope output data model based on data with existing, is then based on
This model application filtering algorithm is filtered gyroscope output data.A current research work part is based on MEMS tops
Spiral shell instrument static data is carried out, and filtering algorithm is only applicable to the Static output data of gyroscope;Other filtering algorithm can
Realize dynamic data filtering, but calculation amount is excessive, it is relatively difficult to achieve for inexpensive real-time system.
The content of the invention
The purpose of the present invention is being handled mainly for the dynamical output data of inexpensive MEMS gyroscope system, so as to
Improve the filtering method of the MEMS gyroscope dynamic data of gyroscope dynamical output data precision.
The method of the present invention for improving MEMS gyroscope dynamical output data precision comprises the following steps:
(1) MEMS gyroscope constant value drift is determined:When MEMS gyroscope steady operation, one hour quiet is gathered respectively in four times
State output data calculates constant value drift of the average of four MEMS gyroscope static datas as gyroscope, and what is gathered later is every
A data are required for cutting this constant value drift;
(2) MEMS gyroscope output data model is established:Gyroscope dynamical output data sequence is Xk, dynamic data difference sequence
For Dk, wherein Dk=Xk+1-Xk;
(3) Kalman filter process-noise variance Q and measurement noise variance R are determined:Process noise and measurement noise are and the time
Unrelated white Gaussian noise, their variance Q and R are constant, and the value needs of Q and R are determined by testing;
(4) Kalman filter is carried out to dynamic data:The Kalman filter model of dynamic data is the MEMS established in step 2
Gyroscope output data model.
MEMS gyroscope output data model of the present invention is:
Wherein state vector is MEMS gyroscope Output speed signal XkWith angular velocity difference sub-signal Dk, corresponding state is defeated
Go out respectively xkAnd dk, Wk, VkRespectively systematic procedure noise and measurement noise act on differential signal DkOn.
The present invention to dynamic data carry out Kalman filter the step of be:
(1) state one-step prediction:This when, is predicted by the optimal State Estimation value of last moment
The state at quarter;System mode vectorState-transition matrix
(2) one-step prediction varivance matrix is calculated:It uses
One moment estimation error variance and process-noise variance calculate one-step prediction error variance;Systematic procedure noise inputs matrix
(3) filtering gain matrix are calculated:Pass through one-step prediction variance and survey
It measures noise variance and calculates Kalman filter gain;Systematic survey matrix
(4) state estimation:It is corrected using the reality output data of gyroscope
One-step prediction is as a result, so as to obtain optimal estimation;ZkFor the angular velocity signal and angular speed difference of k moment gyroscope reality outputs
The vector of signal composition;
(5) estimation error variance matrix is calculated:Pk=[I-KkHk]PK, k-1;This estimation error variance is calculated, is next time
Estimation is prepared;Unit diagonal matrix
(6) state output:
The invention discloses a kind of filtering methods of MEMS gyroscope dynamic data, are determined that constant value is floated using static data
It moves, the data using angular velocity signal and angular velocity difference sub-signal as quantity of state is established based on MEMS gyroscope data output characteristic
Output model carries out Kalman filter to dynamic data based on this model, improves the precision of gyroscope dynamical output data.This
The Kalman filter method that invention uses is a kind of iterative estimate method, only needs the estimate of last moment and the measurement at this moment
Value can provide the angular speed optimal estimation value at this moment, and calculation amount is small, is easy to low-cost system Project Realization.
Description of the drawings
Fig. 1 is the flow chart of gyroscope dynamic data filtering algorithm of the present invention;
Fig. 2 is Kalman filter process, I, II, III in figure, IV, V, and correspond to respectively in technical scheme steps 4 six of VI walk
Suddenly;
Fig. 3 a are first group of static data;
Fig. 3 b are second group of static data;
Fig. 3 c are the 3rd group of static data;
Fig. 3 d are the 4th group of static data;
Fig. 4 a are first group of data mean value;
Fig. 4 b are second group of data mean value;
Fig. 4 c are the 3rd group of data mean value;
Fig. 4 d are the 4th group of data mean value;
Fig. 5 a are first group of data variance;
Fig. 5 b are second group of data variance;
Fig. 5 c are the 3rd group of data variance;
Fig. 5 d are the 4th group of data variance;
Fig. 6 a are first group of data auto-correlation function;
Fig. 6 b are second group of data auto-correlation function;
Fig. 6 c are the 3rd group of data auto-correlation function;
Fig. 6 d are the 4th group of data auto-correlation function;
Fig. 7 is comparison diagram before and after dynamic data filtering.
Specific embodiment
The method of the present invention for improving MEMS gyroscope dynamical output data precision comprises the following steps:
(1) MEMS gyroscope constant value drift is determined:When MEMS gyroscope steady operation, one hour quiet is gathered respectively in four times
State output data calculates constant value drift of the average of four MEMS gyroscope static datas as gyroscope, and what is gathered later is every
A data are required for cutting this constant value drift;
(2) MEMS gyroscope output data model is established:Gyroscope dynamical output data sequence is Xk, dynamic data difference sequence
For Dk, wherein Dk=Xk+1-Xk;
(3) Kalman filter process-noise variance Q and measurement noise variance R are determined:Process noise and measurement noise are and the time
Unrelated white Gaussian noise, their variance Q and R are constant, and the value needs of Q and R are determined by testing;
(4) Kalman filter is carried out to dynamic data:The Kalman filter model of dynamic data is the MEMS established in step 2
Gyroscope output data model.
MEMS gyroscope output data model of the present invention is:
Wherein state vector is MEMS gyroscope Output speed signal XkWith angular velocity difference sub-signal Dk, corresponding state is defeated
Go out respectively xkAnd dk, Wk, VkRespectively systematic procedure noise and measurement noise act on differential signal DkOn.
The present invention to dynamic data carry out Kalman filter the step of be:
(1) state one-step prediction:This when, is predicted by the optimal State Estimation value of last moment
The state at quarter;System mode vectorState-transition matrix
(2) one-step prediction varivance matrix is calculated:It uses
One moment estimation error variance and process-noise variance calculate one-step prediction error variance;Systematic procedure noise inputs matrix
(3) filtering gain matrix are calculated:Pass through one-step prediction variance and survey
It measures noise variance and calculates Kalman filter gain;Systematic survey matrix
(4) state estimation:It is corrected using the reality output data of gyroscope
One-step prediction is as a result, so as to obtain optimal estimation;ZkFor the angular velocity signal and angular speed difference of k moment gyroscope reality outputs
The vector of signal composition;
(5) estimation error variance matrix is calculated:Pk=[I-KkHk]PK, k-1;This estimation error variance is calculated, is next time
Estimation is prepared;Unit diagonal matrix
(6) state output:
Said program is described further with reference to specific embodiments and the drawings, the present invention includes but are not limited to down
State embodiment:
1st, MEMS gyroscope constant value drift is determined:When MEMS gyroscope steady operation, one hour quiet is gathered respectively in four times
State output data, Fig. 3 are the four groups of MEMS gyroscope static datas collected.Calculate the equal of four MEMS gyroscope static datas
It is worth the constant value drift as gyroscope.The average that the present embodiment collects static data is -0.8375 °/s.What is gathered later is every
A data are required for cutting this constant value drift.Judge whether the method for steady operation is to examine static data to MEMS gyroscope
Stationarity, stationarity refer to Random time sequence { xtThe statistical property of (t=1,2,3...) do not change over time, leads here
Examine mean μ=E (xt), the variances sigma of static data2=Var (xt) and covariance function γk=Cov (xt, xt+k) at any time
Situation of change, when meet the time increase when, average tends to definite value and does not change over time, variance tend to definite value not at any time and
Change, when covariance is unrelated with the time, is only related with time interval, judge Random time sequence for stationary sequence.If Fig. 4 is four
Group data mean value changes with time, and Fig. 5 changes with time for four groups of data variances, Fig. 6 for four groups of data covariances at any time
Between variation, as seen from the figure the four of the present embodiment group static data meet stationarity.When the static data of this step acquisition has
During stationarity, the constant value drift calculated could be closer to real constant value drift.
2nd, MEMS gyroscope output data model is established:Gyroscope dynamical output data sequence is Xk, dynamic data difference
Sequence is Dk, wherein Dk=Xk+1-Xk.Since gyroscope actual angular speed is a signal become slowly, from k-1 moment to k
The angular speed increment at moment and the angular speed increment approximately equal from the k moment to the k+1 moment, i.e. Dk=Dk-1.So MEMS gyro
Instrument output data model is:
Wherein state vector is MEMS gyroscope Output speed signal XkWith angular velocity difference sub-signal Dk, corresponding state is defeated
Go out respectively xkAnd dk, Wk, VkRespectively system noise and measurement noise act on differential signal DkOn.What this step was established
Output data model is the model that Kalman filter is used in the 4th step.
3rd, Kalman filter process-noise variance Q and measurement noise variance R are determined:Process noise and measurement noise be and when
Between unrelated white Gaussian noise, their variance Q and R is constant, and the value needs of Q and R are determined by testing.It sets respectively
The value of multigroup Q and R takes the value of the preferable one group of Q and R of static data filter effect as process-noise variance and measurement noise side
Difference.Here judge the value of Q and R with static data filter effect and be without the reason for dynamic data, since top can not be obtained
The actual value of spiral shell instrument dynamic data, therefore the quality of dynamic filter effect can not be analyzed, and static data actual value is 0, filtering
Effect can represent that variance is smaller, and filter effect is better with the variance of filtered data and static data actual value.This implementation
When taking in exampleWhen, Static Filtering effect is preferable, if Fig. 7 is four groups of static numbers
According to filter effect figure.The variance of static data such as table 1 before and after filtering.
Table 1
Data | First group | Second group | 3rd group | 4th group |
Poor in front of filtering (°/s)2 | 1.9745 | 1.9604 | 1.9450 | 1.9612 |
Variance after filtering (°/s)2 | 0.1643 | 0.1708 | 0.1689 | 0.1719 |
4th, Kalman filter is carried out to dynamic data:
The Kalman filter model of dynamic data is the MEMS gyroscope output data model established in step 2:
It is to the step of dynamic data progress Kalman filter:
(1) state one-step prediction:It is predicted by the optimal State Estimation value of last moment
The state at this moment.System mode vectorState-transition matrix
(2) one-step prediction varivance matrix is calculated:With
Last moment estimation error variance and process-noise variance calculate one-step prediction error variance.System noise input matrixSystematic procedure noise variance
(3) filtering gain matrix are calculated:Pass through one-step prediction variance
Kalman filter gain is calculated with measurement noise variance.Systematic survey matrixMeasurement noise variance
(4) state estimation:Come using the reality output data of gyroscope
One-step prediction is corrected as a result, so as to obtain optimal estimation.ZkFor the angular velocity signal and angular speed of k moment gyroscope reality outputs
The vector of differential signal composition.
(5) estimation error variance matrix is calculated:Pk=[I-KkHk]PK, k-1.This estimation error variance is calculated, for next time
Estimation prepare.Unit diagonal matrix
(6) state output:
It is above Kalman filter step, since Kalman filter is an iterative process, gives ability after filtering initial value
Start filtering, it is therefore desirable to determine state initial valueWith estimation error variance initial value P0.Selection state initial valueWhen should not
It is excessive to deviate true initial state, estimation error variance initial value P0It is chosen in the range of (0,1).The present embodiment chooses state initial valueChoose estimation error variance initial valueFig. 7 is that one group of dynamic data of the present embodiment filters
Front and rear comparison diagram, as seen from the figure filtering method proposed by the present invention effectively inhibit in gyroscope dynamical output data with
Chance error is poor, improves the precision of gyroscope dynamical output data.
Fig. 3 (a), (b), (c), four groups of static datas that (d) is the MEMS gyroscope in the present embodiment, every group of static data
Length be one hour.The average of four groups of data is calculated to obtain the constant value drift of the present embodiment MEMS gyroscope.When acquisition
When data volume is bigger, obtained constant value drift is closer to actual constant value drift.The present embodiment thinks that gathering four groups of data obtains
The close enough truthful data of constant value drift constant value drift.
Fig. 4, Fig. 5, Fig. 6 are respectively the mean function, variance function and auto-correlation function of four groups of static datas, by can in figure
To obtain, when increasing the time, the average of four groups of numbers tends to definite value, and variance tends to definite value, and auto-correlation function is unrelated with time point,
It is only related with time interval, it is taken as that the four groups of static datas collected in the present embodiment meet stationarity.
Fig. 7 is comparison diagram before and after the filtering of MEMS gyroscope dynamic data in the present embodiment, and fine line is initial data in figure,
Heavy line is filtered data, and filtered dynamic data is more smooth as can be seen from Figure, preferably eliminates and makes an uproar at random
Sound, hence it is demonstrated that the validity of the method for the present invention.
Claims (3)
1. a kind of filtering method of MEMS gyroscope dynamic data, it is characterised in that:Its step is:
(1) MEMS gyroscope constant value drift is determined:When MEMS gyroscope steady operation, one hour quiet is gathered respectively in four times
State output data calculates constant value drift of the average of four MEMS gyroscope static datas as gyroscope, and what is gathered later is every
A data are required for cutting this constant value drift;
(2) MEMS gyroscope output data model is established:Gyroscope dynamical output data sequence is Xk, dynamic data difference sequence
For Dk, wherein Dk=Xk+1-Xk;
(3) Kalman filter process-noise variance Q and measurement noise variance R are determined:Process noise and measurement noise are and the time
Unrelated white Gaussian noise, their variance Q and R are constant, and the value needs of Q and R are determined by testing;
(4) Kalman filter is carried out to dynamic data:The Kalman filter model of dynamic data is the MEMS established in step 2
Gyroscope output data model.
2. the filtering method of MEMS gyroscope dynamic data according to claim 1, it is characterised in that:MEMS gyroscope is defeated
Going out data model is:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>X</mi>
<mi>k</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>D</mi>
<mi>k</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>X</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>D</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>+</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<msub>
<mi>W</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>x</mi>
<mi>k</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>d</mi>
<mi>k</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>X</mi>
<mi>k</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>D</mi>
<mi>k</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>+</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<msub>
<mi>V</mi>
<mi>k</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
Wherein state vector is MEMS gyroscope Output speed signal XkWith angular velocity difference sub-signal Dk, corresponding state is defeated
Go out respectively xkAnd dk, Wk, VkRespectively systematic procedure noise and measurement noise act on differential signal DkOn.
3. the filtering method of MEMS gyroscope dynamic data according to claim 1, it is characterised in that:To dynamic data into
The step of row Kalman filter is:
(1) state one-step prediction:This when, is predicted by the optimal State Estimation value of last moment
The state at quarter;System mode vectorState-transition matrix
(2) one-step prediction varivance matrix is calculated:Use one
Moment estimation error variance and process-noise variance calculate one-step prediction error variance;Systematic procedure noise inputs matrix
(3) filtering gain matrix are calculated:Pass through one-step prediction variance and measurement
Noise variance calculates Kalman filter gain;Systematic survey matrix
(4) state estimation:One is corrected using the reality output data of gyroscope
Prediction result is walked, so as to obtain optimal estimation;ZkBelieve for the angular velocity signal and angular speed difference of k moment gyroscope reality outputs
Number composition vector;
(5) estimation error variance matrix is calculated:Pk=[I-KkHk]PK, k-1;This estimation error variance is calculated, for estimating next time
Meter is prepared;Unit diagonal matrix
(6) state output:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711396612.7A CN108120452B (en) | 2017-12-21 | 2017-12-21 | Filtering method for dynamic data of MEMS gyroscope |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711396612.7A CN108120452B (en) | 2017-12-21 | 2017-12-21 | Filtering method for dynamic data of MEMS gyroscope |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108120452A true CN108120452A (en) | 2018-06-05 |
CN108120452B CN108120452B (en) | 2021-12-24 |
Family
ID=62231139
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711396612.7A Active CN108120452B (en) | 2017-12-21 | 2017-12-21 | Filtering method for dynamic data of MEMS gyroscope |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108120452B (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110485222A (en) * | 2019-07-29 | 2019-11-22 | 中国铁路总公司 | A kind of dynamic data inversion method and device |
CN111339494A (en) * | 2020-01-16 | 2020-06-26 | 四川建筑职业技术学院 | Gyroscope data processing method based on Kalman filtering |
CN112632454A (en) * | 2020-12-17 | 2021-04-09 | 长光卫星技术有限公司 | MEMS gyro filtering method based on adaptive Kalman filtering algorithm |
US11680798B2 (en) | 2020-08-24 | 2023-06-20 | Invensense, Inc. | Digital demodulator and complex compensator for MEMS gyroscope |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102662083A (en) * | 2012-03-28 | 2012-09-12 | 北京航空航天大学 | Accelerometer calibration method based on GPS velocity information |
CN103217176A (en) * | 2013-05-06 | 2013-07-24 | 南京智真电子科技有限公司 | Method for eliminating null shift of micro-electromechanical system (MEMS) gyroscope |
CN103808331A (en) * | 2014-03-05 | 2014-05-21 | 北京理工大学 | MEMS (micro-electromechanical system) three-axis gyroscope error calibration method |
CN104215262A (en) * | 2014-08-29 | 2014-12-17 | 南京航空航天大学 | On-line dynamic inertia sensor error identification method of inertia navigation system |
CN104880201A (en) * | 2015-03-26 | 2015-09-02 | 武汉大学 | Automatic calibration method of MEMS gyroscopes |
EP3007429A1 (en) * | 2014-10-10 | 2016-04-13 | Parrot | Mobile apparatus, in particular a rotary-wing drone, provided with a video camera supplying image sequences with dynamic correction of the wobble effect |
CN107014367A (en) * | 2017-04-10 | 2017-08-04 | 中国船舶重工集团公司第七0七研究所 | A kind of method for optimizing dynamic property of optical fiber gyroscope |
US20170254832A1 (en) * | 2011-11-04 | 2017-09-07 | Google Inc. | Calibrating Inertial Sensors Using an Image |
-
2017
- 2017-12-21 CN CN201711396612.7A patent/CN108120452B/en active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20170254832A1 (en) * | 2011-11-04 | 2017-09-07 | Google Inc. | Calibrating Inertial Sensors Using an Image |
CN102662083A (en) * | 2012-03-28 | 2012-09-12 | 北京航空航天大学 | Accelerometer calibration method based on GPS velocity information |
CN103217176A (en) * | 2013-05-06 | 2013-07-24 | 南京智真电子科技有限公司 | Method for eliminating null shift of micro-electromechanical system (MEMS) gyroscope |
CN103808331A (en) * | 2014-03-05 | 2014-05-21 | 北京理工大学 | MEMS (micro-electromechanical system) three-axis gyroscope error calibration method |
CN104215262A (en) * | 2014-08-29 | 2014-12-17 | 南京航空航天大学 | On-line dynamic inertia sensor error identification method of inertia navigation system |
EP3007429A1 (en) * | 2014-10-10 | 2016-04-13 | Parrot | Mobile apparatus, in particular a rotary-wing drone, provided with a video camera supplying image sequences with dynamic correction of the wobble effect |
CN104880201A (en) * | 2015-03-26 | 2015-09-02 | 武汉大学 | Automatic calibration method of MEMS gyroscopes |
CN107014367A (en) * | 2017-04-10 | 2017-08-04 | 中国船舶重工集团公司第七0七研究所 | A kind of method for optimizing dynamic property of optical fiber gyroscope |
Non-Patent Citations (4)
Title |
---|
XIA YAN,ETC: "Research on Random Drift Modeling and A Kalman Filter Based on The Differential Signal of MEMS Gyroscope", 《2013 25TH CHINESE CONTROL AND DECISION CONFERENCE (CCDC)》 * |
XUE LIANG,ETC: "A novel Kalman filter for combining outputs of MEMS gyroscope array", 《MEASUREMENT》 * |
王励扬,等: "低成本MEMS陀螺实时滤波方法", 《电子技术应用》 * |
蒙涛,等: "MEMS陀螺误差建模与滤波方法", 《***工程与电子技术》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110485222A (en) * | 2019-07-29 | 2019-11-22 | 中国铁路总公司 | A kind of dynamic data inversion method and device |
CN110485222B (en) * | 2019-07-29 | 2020-11-27 | 中国铁路总公司 | Dynamic data inversion method and device |
CN111339494A (en) * | 2020-01-16 | 2020-06-26 | 四川建筑职业技术学院 | Gyroscope data processing method based on Kalman filtering |
US11680798B2 (en) | 2020-08-24 | 2023-06-20 | Invensense, Inc. | Digital demodulator and complex compensator for MEMS gyroscope |
CN112632454A (en) * | 2020-12-17 | 2021-04-09 | 长光卫星技术有限公司 | MEMS gyro filtering method based on adaptive Kalman filtering algorithm |
Also Published As
Publication number | Publication date |
---|---|
CN108120452B (en) | 2021-12-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108120452A (en) | The filtering method of MEMS gyroscope dynamic data | |
CN105222809B (en) | The method of the boat magnetic disturbance penalty coefficient estimation of a kind of ground magnetic gradient robust | |
CN101464152B (en) | Adaptive filtering method for SINS/GPS combined navigation system | |
CN107832575A (en) | Band feedback maneuvering target Asynchronous Track Fusion based on pseudo-measurement | |
CN109883392B (en) | Strapdown inertial navigation heave measurement method based on phase compensation | |
CN102540252B (en) | High-precision median stacking method on basis of cross-correlation | |
WO2009105964A1 (en) | Method of pre-stack two-dimension-like transformation of three-dimensional seismic record | |
CN104062642B (en) | A kind of method that laser radar waveform data is carried out Gauss echo decomposition | |
CN106840211A (en) | A kind of SINS Initial Alignment of Large Azimuth Misalignment On methods based on KF and STUPF combined filters | |
CN109724602A (en) | A kind of attitude algorithm system and its calculation method based on hardware FPU | |
CN107688554A (en) | Frequency domain identification method based on adaptive Fourier decomposition | |
CN103968910B (en) | Method for improving motor vehicle fuel quantity estimation precision | |
CN109827571A (en) | A kind of dual acceleration meter calibration method under the conditions of no turntable | |
CN102722102B (en) | H-infinity feedback and iterative learning coordination control method based on error analysis | |
CN104613966B (en) | A kind of cadastration off-line data processing method | |
CN104597460A (en) | Beidou satellite navigation receiver based carrier wave tracking loop crystal oscillator acceleration speed sensitivity coefficient calibration method | |
CN109489987A (en) | Fanjet measurement biases fault-tolerant gas circuit performance distributed and filters estimation method | |
CN101929862A (en) | Method for determining initial attitude of inertial navigation system based on Kalman filtering | |
CN108037317B (en) | Dynamic decoupling method and system of accelerometer | |
CN107168359A (en) | A kind of control method of large-scale flexible satellite high precision high stability degree | |
CN106197376B (en) | Car body obliqueness measurement method based on uniaxial MEMS inertial sensor | |
CN109067381A (en) | A kind of Real-Time Filtering system and method for MEMS gyroscope random noise | |
CN106053884B (en) | Dual-axis resonance type accelerometer digital control circuit device based on Coordinate Rotating Algorithm | |
CN110677140A (en) | Random system filter containing unknown input and non-Gaussian measurement noise | |
CN105698799B (en) | A kind of optimal FIR filter of pretreatment for improving Strapdown Inertial Navigation System attitude accuracy |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |