CN112197767A - Filter design method for improving filtering error on line - Google Patents
Filter design method for improving filtering error on line Download PDFInfo
- Publication number
- CN112197767A CN112197767A CN202011078829.5A CN202011078829A CN112197767A CN 112197767 A CN112197767 A CN 112197767A CN 202011078829 A CN202011078829 A CN 202011078829A CN 112197767 A CN112197767 A CN 112197767A
- Authority
- CN
- China
- Prior art keywords
- matrix
- noise
- error
- filter
- simplified
- 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
- 238000000034 method Methods 0.000 title claims abstract description 52
- 238000013461 design Methods 0.000 title claims abstract description 31
- 238000001914 filtration Methods 0.000 title claims description 39
- 239000011159 matrix material Substances 0.000 claims abstract description 113
- 238000012937 correction Methods 0.000 claims abstract description 17
- 238000005259 measurement Methods 0.000 claims description 33
- 230000007704 transition Effects 0.000 claims description 21
- 239000013598 vector Substances 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 8
- 238000012546 transfer Methods 0.000 claims description 3
- 238000003491 array Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
Images
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
- G01C21/165—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 combined with non-inertial navigation instruments
-
- 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
-
- 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
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/42—Determining position
- G01S19/45—Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement
- G01S19/47—Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement the supplementary measurement being an inertial measurement, e.g. tightly coupled inertial
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Manufacturing & Machinery (AREA)
- Computer Networks & Wireless Communication (AREA)
- Navigation (AREA)
Abstract
A filter design method for improving filter error on line, according to the application requirement of inertial navigation system and GPS receiver, adopting indirect closed-loop correction under position and speed loose combination mode to process the filter estimated value obtained by indirect method, then simplifying the design of filter estimated value, substituting the coefficient matrix after simplification into indirect closed-loop, and continuously carrying out on-line estimation or correction on system model parameter, noise statistical characteristic and state gain array, thereby realizing on-line improvement of filter design parameter to reduce actual filter error.
Description
Technical Field
The invention relates to the technical field of filter design, in particular to a filter design method for improving a filtering error on line.
Background
In a GPS/INS integrated navigation system, navigation parameters are estimated by Kalman filtering, and an important prerequisite for the application of a classical Kalman filter is that the statistical characteristics of noise must be accurately known, but because of the instability of components (an accelerometer and a gyroscope) of an actual system, the estimation environment is not constant all the time, such as multipath errors and SA errors of a GPS, the accurate description of system noise and observation noise is difficult.
Disclosure of Invention
The technical problem to be solved by the present invention is to provide a filter design method for improving the filtering error on line, so as to solve the above problems in the background art.
The technical problem solved by the invention is realized by adopting the following technical scheme:
a filter design method for improving filter error on line, according to the application requirement of inertial navigation system and GPS receiver, adopting indirect closed-loop correction under position and speed loose combination mode to process the filtering estimated value obtained by indirect method, then simplifying the design of the filtering estimated value, substituting the coefficient matrix after simplification into indirect closed-loop, and continuously carrying out on-line estimation or correction on system model parameter, noise statistical characteristic and state gain array, thereby realizing on-line improvement of filter design parameter, so as to reduce the actual filter error and achieve the purpose of improving the accuracy and reliability of the navigation system; the method comprises the following specific steps:
1) processing an indirect filtering estimated value by indirect closed-loop correction in a position and speed loose combination mode according to the application requirements of an inertial navigation system and a GPS receiver;
selecting a state vector of a system state equation
Selecting 3 position error components, 3 speed error components and 3 attitude angle error components from the state quantities of the position and speed loose combination mode, solving colored noise by adopting a state expansion method, and expanding the errors of a gyroscope and an accelerometer into an error equation as the state quantities:
calculating to obtain a three-dimensional platform error angle, a three-dimensional speed error and a three-dimensional position error, then taking the three-dimensional platform error angle, the three-dimensional speed error and the three-dimensional position error as state vectors of a system state equation, obtaining respective estimated values by Kalman filtering, and correcting the output of an inertial navigation system through the estimated values, thereby improving the navigation precision;
② selecting system white noise
Taking the white noise drift (omega) of the gyroscopegx、ωgy、ωgz) Accelerometer random white noise drift (omega)ax、ωay、ωaz) And gyroscope first order Markov driven white noise (ω)rx、ωry、ωrz) As system white noise w (t):
W(t)=[ωgx ωgy ωgz ωrx ωry ωrz ωax ωay ωaz]T;
setting system state transfer matrix and system noise coefficient matrix
Assigning values to a system state transition matrix F (t) and a system noise coefficient matrix G (t);
selecting the measurement vector of the system measurement equation
Wherein the subscript I denotes information of the inertial navigation system, the subscript G denotes information of the GPS,
because the combined mode of position and speed is adopted, two groups of observed values are shared: one group is position observation values, namely longitude and latitude and height information given by an inertial navigation system and corresponding information difference values given by a GPS receiver, and the other group is speed difference values given by the inertial navigation system and the GPS receiver respectively;
taking a speed error in the northeast direction and a position error of the longitude and latitude height of the GPS receiver as measurement white noise;
setting measurement noise coefficient matrix
And assigning a value to the measurement noise coefficient matrix N (t);
2) simplifying the design of the system state transition matrix F (t), the system noise coefficient matrix G (t) and the measurement noise coefficient matrix N (t) according to the engineering state and the requirement
The kalman filter equation is as follows:
k (k) -the filter gain array at time k;
z (k) -the measured value at time k;
p (k/k-1) -a covariance matrix of the prediction error estimates at time k-1 versus time k;
p (k) -a covariance matrix of error estimates at time k;
q (k) -system noise variance matrix;
r (k) -system observation noise variance matrix;
a) design simplification system state transition matrix phik,k-1
Because the engineering object is generally a continuous system, but the Kalman filtering is commonly used for describing a system by a discretization model, the state transition matrix phi of the system is firstly determinedk,k-1Carrying out discretization treatment:
Φk,k-1=eFT
wherein, T is a filtering period, and F is a system state transition matrix F (T);
considering that the filter period T is generally short, phik,k-1It can also be simplified as follows:
Φk,k-1=I+FT+O(T2)≈I+FT(5)
The system noise variance required in the filtering calculation is in the form ofTherefore, Q is not generally calculated separatelykInstead, a continuous system of G (t) Q (t) GT(t) direct calculation of
Wherein, G (t) is a system noise coefficient matrix, Q (t) is a variance matrix corresponding to system white noise W (t), Q (t) is a diagonal matrix, and the diagonal elements of the diagonal matrix are the variances of the white noise;
wherein T is the filter period, phik,k-1To simplify the system state transition matrix;
c) design simplification measurement noise error equation matrix Rk
RkFor measuring the variance matrix corresponding to the noise coefficient matrix N (t), and RkIs a diagonal matrix whose diagonal elements are the variance of the process noise;
Initial state of filteringTaking different values as required, taking zero vector for convenient calculation, filteringInitial error variance matrix P0The method is characterized in that the diagonal matrix is a diagonal matrix, and diagonal elements of the diagonal matrix are determined according to the range of navigation parameter initial errors, attitude misalignment angle errors, gyro drift errors and accelerometer deviations in combination with engineering;
4) design of Kalman filter using simplified coefficient matrix
The simplified system state transition matrix phi obtained by the processing in the step 2)k,k-1Simplified process noise variance matrixSimplified measurement noise error equation matrix RkAnd the initial state determined in step 3)And the initial error variance matrix P0The method substitutes indirect closed loop to carry out correction, and continuously carries out online estimation or correction on unknown or unknown system model parameters, noise statistical characteristics and state gain arrays while carrying out filtering by utilizing observation data, thereby realizing online improvement of filter design parameters, reducing actual filtering errors and achieving the purpose of improving the precision and reliability of a navigation system.
Has the advantages that: the method comprises the steps of processing a filtering estimated value obtained by an indirect method by adopting an indirect method closed loop in a position and speed loose combination mode, then simply designing the filtering estimated value, determining an initial error range by combining engineering experience in the process of developing a missile system, substituting a simplified coefficient matrix and the initial error into the indirect method closed loop, and continuously carrying out online estimation or correction on system model parameters, noise statistical characteristics and a state gain array, thereby realizing online improvement of filter design parameters and reducing the actual filtering error; the filter designed by the invention has a simple structure, is convenient for the concrete realization of actual engineering, is beneficial to the implementation of calculation by a computer, has a good filtering effect, and can effectively improve the filtering precision.
Drawings
FIG. 1 is a schematic diagram of indirect closed loop in the preferred embodiment of the present invention.
FIG. 2 is a schematic diagram of an inertial navigation system according to a preferred embodiment of the present invention.
Detailed Description
In order to make the technical means, the creation characteristics, the achievement purposes and the effects of the invention easy to understand, the invention is further explained below by combining the specific drawings.
A filter design method for improving filter error on line, firstly, according to the application requirements of an inertial navigation system and a GPS receiver, processing the filter estimation value obtained by an indirect method by indirect method closed loop correction under a position and speed loose combination mode, then simplifying the design of the filter estimation value, substituting the coefficient matrix after simplification into the indirect method closed loop, and continuously carrying out on-line estimation or correction on system model parameters, noise statistical characteristics and a state gain array, thereby realizing the on-line improvement of filter design parameters to reduce the actual filter error; the method comprises the following specific steps:
1) according to the application requirements of an inertial navigation system and a GPS receiver, indirect closed-loop correction in a position and speed loose combination mode is adopted to process indirect filtering estimated values, and the principle is shown in figure 1;
selecting a state vector of a system state equation
Selecting 3 position error components, 3 speed error components and 3 attitude angle error components from the state quantities of the position and speed loose combination mode, solving colored noise by adopting a state expansion method, and expanding the errors of a gyroscope and an accelerometer into an error equation as the state quantities:
calculating to obtain a three-dimensional platform error angle, a three-dimensional speed error and a three-dimensional position error, then taking the three-dimensional platform error angle, the three-dimensional speed error and the three-dimensional position error as state vectors of a system state equation, obtaining respective estimated values by Kalman filtering, and correcting the output of an inertial navigation system through the estimated values, thereby improving the navigation precision;
② selecting system white noise
Taking the white noise drift (omega) of the gyroscopegx、ωgy、ωgz) Accelerometer random white noise drift (omega)ax、ωay、ωaz) And gyroscope first order Markov driven white noise (ω)rx、ωry、ωrz) As system white noise w (t):
W(t)=[ωgx ωgy ωgz ωrx ωry ωrz ωax ωay ωaz]T;
setting system state transfer matrix and system noise coefficient matrix
Assigning values to a system state transition matrix F (t) and a system noise coefficient matrix G (t);
selecting the measurement vector of the system measurement equation
Wherein the subscript I denotes information of the inertial navigation system, the subscript G denotes information of the GPS,
because this embodiment adopts the position, speed combination mode, so there are two sets of observed values in total: one group is position observation values, namely longitude and latitude and height information given by an inertial navigation system and corresponding information difference values given by a GPS receiver, and the other group is speed difference values given by the inertial navigation system and the GPS receiver respectively;
taking a speed error in the northeast direction and a position error of the longitude and latitude height of the GPS receiver as measurement white noise;
setting measurement noise coefficient matrix
And assigning a value to the measurement noise coefficient matrix N (t);
2) simplifying the design of the system state transition matrix F (t), the system noise coefficient matrix G (t) and the measurement noise coefficient matrix N (t) according to the engineering state and the requirement
The kalman filter equation is designed as follows:
k (k) -the filter gain array at time k;
z (k) -the measured value at time k;
p (k/k-1) -a covariance matrix of the prediction error estimates at time k-1 versus time k;
p (k) -a covariance matrix of error estimates at time k;
q (k) -system noise variance matrix;
r (k) -system observation noise variance matrix;
a) design simplification system state transition matrix phik,k-1
Because the engineering object is generally a continuous system, but the Kalman filtering is commonly used for describing a system by a discretization model, the state transition matrix phi of the system is firstly determinedk,k-1Carrying out discretization treatment:
Φk,k-1=eFT
wherein, T is a filtering period, and F is a system state transition matrix F (T);
considering that the filter period T is generally short, phik,k-1It can also be simplified as follows:
Φk,k-1=I+FT+O(T2)≈I+FT (5)
The system noise variance required in the filtering calculation is in the form ofTherefore, Q is not generally calculated separatelykInstead, a continuous system of G (t) Q (t) GT(t) direct calculation of
Wherein, G (t) is a system noise coefficient matrix, Q (t) is a variance matrix corresponding to system white noise W (t), Q (t) is a diagonal matrix, and the diagonal elements of the diagonal matrix are the variances of the white noise;
wherein T is the filter period, phik,k-1To simplify the system state transition matrix;
because the measurement information of the GPS in the SINS/GPS integrated navigation system is provided at the time interval of about 1 second, the measurement equation is discrete, and the discretization processing is not needed;
c) design simplification measurement noise error equation matrixRk
RkFor measuring the variance matrix corresponding to the noise coefficient matrix N (t), and RkIs a diagonal matrix whose diagonal elements are the variance of the process noise;
Initial state of filteringTaking different values as required, taking the values as zero vectors for convenient calculation, and taking the initial error variance matrix P of filtering0The method is characterized in that the diagonal matrix is a diagonal matrix, and diagonal elements of the diagonal matrix are determined according to the range of navigation parameter initial errors, attitude misalignment angle errors, gyro drift errors and accelerometer deviations in combination with engineering;
4) design of Kalman filter using simplified coefficient matrix
As shown in fig. 2, the simplified system state transition matrix Φ obtained by the processing in step 2) is usedk,k-1Simplified process noise variance matrixSimplified measurement noise error equation matrix RkAnd the initial state determined in step 3)And the initial error variance matrix P0The method substitutes indirect closed loop to carry out correction, and continuously carries out online estimation or correction on unknown or unknown system model parameters, noise statistical characteristics and state gain arrays while carrying out filtering by utilizing observation data, thereby realizing online improvement of filter design parameters, reducing actual filtering errors and achieving the purpose of improving the precision and reliability of a navigation system.
Claims (8)
1. A filter design method for improving filter error on line is characterized in that according to application requirements of an inertial navigation system and a GPS receiver, indirect closed-loop correction under a position and speed loose combination mode is adopted to process filter estimated values obtained by an indirect method, then simplified design is carried out on the filter estimated values, simplified coefficient matrixes are substituted into indirect closed-loop, and on-line estimation or correction is carried out on system model parameters, noise statistical characteristics and state gain matrixes continuously, so that filter design parameters are improved on line, and actual filter errors are reduced.
2. The method for designing the filter for improving the filtering error online according to claim 1, comprising the following steps:
1) processing an indirect filtering estimated value by indirect closed-loop correction in a position and speed loose combination mode according to the application requirements of an inertial navigation system and a GPS receiver;
selecting a state vector of a system state equation
Selecting 3 position error components, 3 speed error components and 3 attitude angle error components from the state quantities of the position and speed loose combination mode, solving colored noise by adopting a state expansion method, and expanding the errors of a gyroscope and an accelerometer into an error equation as the state quantities:
taking the obtained three-dimensional platform error angle, three-dimensional speed error and three-dimensional position error as state vectors of a system state equation, and obtaining respective estimated values by Kalman filtering;
② selecting system white noise
Taking the white noise drift (omega) of the gyroscopegx、ωgy、ωgz) Accelerometer random white noise drift (omega)ax、ωay、ωaz) And gyroscope first order Markov driven white noise (ω)rx、ωry、ωrz) As system white noise w (t):
W(t)=[ωgx ωgy ωgz ωrx ωry ωrz ωax ωay ωaz]T;
setting system state transfer matrix and system noise coefficient matrix
Assigning values to a system state transition matrix F (t) and a system noise coefficient matrix G (t);
selecting the measurement vector of the system measurement equation
Because of the position and speed combination mode, two groups of observed values are shared: one group is position observation values, namely longitude and latitude and height information given by an inertial navigation system and corresponding information difference values given by a GPS receiver, and the other group is speed difference values given by the inertial navigation system and the GPS receiver respectively;
taking a speed error in the northeast direction and a position error of the longitude and latitude height of the GPS receiver as measurement white noise;
setting measurement noise coefficient matrix
And assigning a value to the measurement noise coefficient matrix N (t);
2) simplifying the system state transition matrix F (t), the system noise coefficient matrix G (t) and the measurement noise coefficient matrix N (t) according to the engineering state and the requirement to obtain a simplified system state transition matrix phik,k-1Simplified process noise variance matrixSimplified measurement noise error equation matrix Rk;
4) Design of Kalman filter using simplified coefficient matrix
The simplified system state transition matrix phi obtained by the processing in the step 2)k,k-1Simplified process noise variance matrixSimplified measurement noise error equation matrix RkAnd the initial state determined in step 3)And the initial error variance matrix P0And substituting the closed loop of an indirect method, and continuously carrying out online estimation or correction on unknown or unknown system model parameters, noise statistical characteristics and a state gain array so as to realize online improvement of filter design parameters.
4. The method for designing the filter for improving the filtering error online according to claim 2, wherein in step 2), the kalman filtering equation is as follows:
k (k) -the filter gain array at time k;
z (k) -the measured value at time k;
p (k/k-1) -a covariance matrix of the prediction error estimates at time k-1 versus time k;
p (k) -a covariance matrix of error estimates at time k;
q (k) -system noise variance matrix;
r (k) -system observation noise variance matrix;
derived from equation (4): the simplified system state transition matrix phik,k-1:
Φk,k-1=I+FT+O(T2)≈I+FT (5)
Wherein, T is a filtering period, and F is a system state transition matrix F (T);
Wherein T is the filter period, phik,k-1To simplify the system state transition matrix.
5. The method as claimed in claim 4, wherein the simplified process noise variance matrix is a simplified process noise variance matrixIs calculated asThe following:
the system noise variance required in the filtering calculation is in the form ofTherefore, Q is not generally calculated separatelykInstead, a continuous system of G (t) Q (t) GT(t) direct calculation of
Wherein, g (t) is a system noise coefficient matrix, q (t) is a variance matrix corresponding to the system white noise w (t), q (t) is a diagonal matrix, and the diagonal element thereof is the variance of the white noise, and the formula (7) is obtained by simplifying the formula (6).
6. The method as claimed in claim 2, wherein in step 2), the simplified measurement noise error equation matrix R is used for designing the filter for improving the filtering error onlinekThe corresponding variance matrix of the noise coefficient matrix N (t) is measured.
7. The method as claimed in claim 6, wherein the simplified measurement noise error equation matrix R is a simplified measurement noise error equation matrixkIs a diagonal matrix.
8. The method of claim 7, wherein the simplified measurement noise error equation matrix R is used for designing a filter for improving filtering errors on linekThe diagonal element is the variance of the process noise.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011078829.5A CN112197767B (en) | 2020-10-10 | 2020-10-10 | Filter design method for improving filtering error on line |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011078829.5A CN112197767B (en) | 2020-10-10 | 2020-10-10 | Filter design method for improving filtering error on line |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112197767A true CN112197767A (en) | 2021-01-08 |
CN112197767B CN112197767B (en) | 2022-12-23 |
Family
ID=74013308
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011078829.5A Active CN112197767B (en) | 2020-10-10 | 2020-10-10 | Filter design method for improving filtering error on line |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112197767B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116165885A (en) * | 2022-11-29 | 2023-05-26 | 华东交通大学 | Model-free adaptive robust control method and system for high-speed train |
Citations (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050060092A1 (en) * | 2003-08-05 | 2005-03-17 | The Boeing Company | Laser range finder closed-loop pointing technology of relative navigation, attitude determination, pointing and tracking for spacecraft rendezvous |
FR2906893A1 (en) * | 2006-10-06 | 2008-04-11 | Thales Sa | METHOD AND DEVICE FOR MONITORING THE INTEGRITY OF INFORMATION DELIVERED BY AN INS / GNSS HYBRID SYSTEM |
FR2927705A1 (en) * | 2008-02-19 | 2009-08-21 | Thales Sa | HYBRIDIZATION NAVIGATION SYSTEM BY PHASE MEASUREMENTS |
US20120326926A1 (en) * | 2011-06-24 | 2012-12-27 | Mayflower Communications Company, Inc. | High sensitivity gps/gnss receiver |
WO2013080183A1 (en) * | 2011-11-30 | 2013-06-06 | Applanix Corporation | A quasi tightly coupled gnss-ins integration process |
CN104132662A (en) * | 2014-07-25 | 2014-11-05 | 辽宁工程技术大学 | Closed-loop Kalman filter inertial positioning method based on zero velocity update |
CN104181572A (en) * | 2014-05-22 | 2014-12-03 | 南京理工大学 | Missile-borne inertia/ satellite tight combination navigation method |
US20140372063A1 (en) * | 2012-03-06 | 2014-12-18 | Wuhan University | Quick calibration method for inertial measurement unit |
CN106501832A (en) * | 2016-12-16 | 2017-03-15 | 南京理工大学 | A kind of fault-tolerant vector tracking GNSS/SINS deep integrated navigation methods |
CN106646575A (en) * | 2016-11-15 | 2017-05-10 | 南京航空航天大学 | FPGA double-core-based integrated navigation system and construction method thereof |
CN107621264A (en) * | 2017-09-30 | 2018-01-23 | 华南理工大学 | The method for adaptive kalman filtering of vehicle-mounted micro- inertia/satellite combined guidance system |
CN108426574A (en) * | 2018-02-02 | 2018-08-21 | 哈尔滨工程大学 | A kind of MEMS pedestrian navigation methods of the course angle correction algorithm based on ZIHR |
CN109000642A (en) * | 2018-05-25 | 2018-12-14 | 哈尔滨工程大学 | A kind of improved strong tracking volume Kalman filtering Combinated navigation method |
CN109556632A (en) * | 2018-11-26 | 2019-04-02 | 北方工业大学 | INS/GNSS/polarization/geomagnetic integrated navigation alignment method based on Kalman filtering |
US20190129044A1 (en) * | 2016-07-19 | 2019-05-02 | Southeast University | Cubature Kalman Filtering Method Suitable for High-dimensional GNSS/INS Deep Coupling |
CN110165707A (en) * | 2019-02-26 | 2019-08-23 | 国网吉林省电力有限公司 | Light-preserved system optimal control method based on Kalman filtering and Model Predictive Control |
CN110567455A (en) * | 2019-09-25 | 2019-12-13 | 哈尔滨工程大学 | tightly-combined navigation method for quadrature updating of volume Kalman filtering |
CN110702144A (en) * | 2019-10-25 | 2020-01-17 | 江西洪都航空工业集团有限责任公司 | Method for verifying strapdown inertia and GPS satellite integrated navigation system by mounting aircraft |
CN110864688A (en) * | 2019-11-28 | 2020-03-06 | 湖南率为控制科技有限公司 | Attitude heading method for vehicle-mounted azimuth open-loop horizontal attitude angle closed loop |
CN110926465A (en) * | 2019-12-11 | 2020-03-27 | 哈尔滨工程大学 | MEMS/GPS loose combination navigation method |
CN110988926A (en) * | 2019-12-20 | 2020-04-10 | 福建海峡北斗导航科技研究院有限公司 | Method for realizing position accurate fixed point deception migration in loose GNSS/INS combined navigation mode |
CN111006659A (en) * | 2019-12-06 | 2020-04-14 | 江西洪都航空工业集团有限责任公司 | Navigation system with multi-navigation-source information fusion function |
CN111156994A (en) * | 2019-12-31 | 2020-05-15 | 西安航天华迅科技有限公司 | INS/DR & GNSS loose integrated navigation method based on MEMS inertial component |
CN111623779A (en) * | 2020-05-20 | 2020-09-04 | 哈尔滨工程大学 | Time-varying system adaptive cascade filtering method suitable for unknown noise characteristics |
-
2020
- 2020-10-10 CN CN202011078829.5A patent/CN112197767B/en active Active
Patent Citations (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050060092A1 (en) * | 2003-08-05 | 2005-03-17 | The Boeing Company | Laser range finder closed-loop pointing technology of relative navigation, attitude determination, pointing and tracking for spacecraft rendezvous |
FR2906893A1 (en) * | 2006-10-06 | 2008-04-11 | Thales Sa | METHOD AND DEVICE FOR MONITORING THE INTEGRITY OF INFORMATION DELIVERED BY AN INS / GNSS HYBRID SYSTEM |
FR2927705A1 (en) * | 2008-02-19 | 2009-08-21 | Thales Sa | HYBRIDIZATION NAVIGATION SYSTEM BY PHASE MEASUREMENTS |
US20120326926A1 (en) * | 2011-06-24 | 2012-12-27 | Mayflower Communications Company, Inc. | High sensitivity gps/gnss receiver |
WO2013080183A1 (en) * | 2011-11-30 | 2013-06-06 | Applanix Corporation | A quasi tightly coupled gnss-ins integration process |
US20140372063A1 (en) * | 2012-03-06 | 2014-12-18 | Wuhan University | Quick calibration method for inertial measurement unit |
CN104181572A (en) * | 2014-05-22 | 2014-12-03 | 南京理工大学 | Missile-borne inertia/ satellite tight combination navigation method |
CN104132662A (en) * | 2014-07-25 | 2014-11-05 | 辽宁工程技术大学 | Closed-loop Kalman filter inertial positioning method based on zero velocity update |
US20190129044A1 (en) * | 2016-07-19 | 2019-05-02 | Southeast University | Cubature Kalman Filtering Method Suitable for High-dimensional GNSS/INS Deep Coupling |
CN106646575A (en) * | 2016-11-15 | 2017-05-10 | 南京航空航天大学 | FPGA double-core-based integrated navigation system and construction method thereof |
CN106501832A (en) * | 2016-12-16 | 2017-03-15 | 南京理工大学 | A kind of fault-tolerant vector tracking GNSS/SINS deep integrated navigation methods |
CN107621264A (en) * | 2017-09-30 | 2018-01-23 | 华南理工大学 | The method for adaptive kalman filtering of vehicle-mounted micro- inertia/satellite combined guidance system |
CN108426574A (en) * | 2018-02-02 | 2018-08-21 | 哈尔滨工程大学 | A kind of MEMS pedestrian navigation methods of the course angle correction algorithm based on ZIHR |
CN109000642A (en) * | 2018-05-25 | 2018-12-14 | 哈尔滨工程大学 | A kind of improved strong tracking volume Kalman filtering Combinated navigation method |
CN109556632A (en) * | 2018-11-26 | 2019-04-02 | 北方工业大学 | INS/GNSS/polarization/geomagnetic integrated navigation alignment method based on Kalman filtering |
CN110165707A (en) * | 2019-02-26 | 2019-08-23 | 国网吉林省电力有限公司 | Light-preserved system optimal control method based on Kalman filtering and Model Predictive Control |
CN110567455A (en) * | 2019-09-25 | 2019-12-13 | 哈尔滨工程大学 | tightly-combined navigation method for quadrature updating of volume Kalman filtering |
CN110702144A (en) * | 2019-10-25 | 2020-01-17 | 江西洪都航空工业集团有限责任公司 | Method for verifying strapdown inertia and GPS satellite integrated navigation system by mounting aircraft |
CN110864688A (en) * | 2019-11-28 | 2020-03-06 | 湖南率为控制科技有限公司 | Attitude heading method for vehicle-mounted azimuth open-loop horizontal attitude angle closed loop |
CN111006659A (en) * | 2019-12-06 | 2020-04-14 | 江西洪都航空工业集团有限责任公司 | Navigation system with multi-navigation-source information fusion function |
CN110926465A (en) * | 2019-12-11 | 2020-03-27 | 哈尔滨工程大学 | MEMS/GPS loose combination navigation method |
CN110988926A (en) * | 2019-12-20 | 2020-04-10 | 福建海峡北斗导航科技研究院有限公司 | Method for realizing position accurate fixed point deception migration in loose GNSS/INS combined navigation mode |
CN111156994A (en) * | 2019-12-31 | 2020-05-15 | 西安航天华迅科技有限公司 | INS/DR & GNSS loose integrated navigation method based on MEMS inertial component |
CN111623779A (en) * | 2020-05-20 | 2020-09-04 | 哈尔滨工程大学 | Time-varying system adaptive cascade filtering method suitable for unknown noise characteristics |
Non-Patent Citations (6)
Title |
---|
ANTON RENNER等: "Loosely Coupled Integrated Navigation System for Excavators", 《2020 EUROPEAN CONTROL CONFERENCE (ECC)》 * |
LI XIANG等: "External Acceleration Elimination for Complementary Attitude Filter", 《2017 IEEE INTERNATIONAL CONFERENCE ON INFORMATION AND AUTOMATION (ICIA)》 * |
MATTIA DE AGOSTINO & CHIARA PORPORATO等: "Italian low cost GNSS/INS system suitable for mobile mapping", 《IEEE AEROSPACE AND ELECTRONIC SYSTEMS MAGAZINE》 * |
周艳丽等: "基于GPS/MEMS微惯性导航***的混合校正方法研究", 《导航与控制》 * |
李翔等: "航姿参考***姿态算法性能优化仿真研究", 《计算机仿真》 * |
王均晖等: "无人机组合导航直接法与间接法滤波方式比较", 《北京航空航天大学学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116165885A (en) * | 2022-11-29 | 2023-05-26 | 华东交通大学 | Model-free adaptive robust control method and system for high-speed train |
CN116165885B (en) * | 2022-11-29 | 2023-11-14 | 华东交通大学 | Model-free adaptive robust control method and system for high-speed train |
Also Published As
Publication number | Publication date |
---|---|
CN112197767B (en) | 2022-12-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109211276B (en) | SINS initial alignment method based on GPR and improved SRCKF | |
CN110487301B (en) | Initial alignment method of radar-assisted airborne strapdown inertial navigation system | |
CN101949703B (en) | Strapdown inertial/satellite combined navigation filtering method | |
CN112097763B (en) | Underwater vehicle combined navigation method based on MEMS IMU/magnetometer/DVL combination | |
CN100541135C (en) | Fiber-optic gyroscope strapdown inertial navigation system initial attitude based on Doppler is determined method | |
CN103344259B (en) | A kind of INS/GPS integrated navigation system feedback correction method estimated based on lever arm | |
CN102538792B (en) | Filtering method for position attitude system | |
CN106500693B (en) | A kind of AHRS algorithm based on adaptive extended kalman filtering | |
CN104764467B (en) | Re-entry space vehicle inertial sensor errors online adaptive scaling method | |
CN105091907B (en) | DVL orientation alignment error method of estimation in SINS/DVL combinations | |
CN109945895B (en) | Inertial navigation initial alignment method based on fading smooth variable structure filtering | |
Xue et al. | In-motion alignment algorithm for vehicle carried SINS based on odometer aiding | |
CN103674064B (en) | Initial calibration method of strapdown inertial navigation system | |
CN110849360B (en) | Distributed relative navigation method for multi-machine collaborative formation flight | |
CN113340298B (en) | Inertial navigation and dual-antenna GNSS external parameter calibration method | |
CN113959462B (en) | Quaternion-based inertial navigation system self-alignment method | |
CN105606093B (en) | Inertial navigation method and device based on gravity real-Time Compensation | |
CN110567455A (en) | tightly-combined navigation method for quadrature updating of volume Kalman filtering | |
CN114877915B (en) | Device and method for calibrating g sensitivity error of laser gyro inertia measurement assembly | |
CN114777812B (en) | Inter-advancing alignment and attitude estimation method for underwater integrated navigation system | |
CN109489661A (en) | Gyro constant value drift estimation method when a kind of satellite is initially entered the orbit | |
CN111912427B (en) | Method and system for aligning motion base of strapdown inertial navigation assisted by Doppler radar | |
CN109084755B (en) | Accelerometer zero offset estimation method based on gravity apparent velocity and parameter identification | |
CN105300407B (en) | A kind of marine dynamic starting method for single axis modulation laser gyro inertial navigation system | |
CN112197767B (en) | Filter design method for improving filtering error on line |
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 |