CN1908584A - Method for determining initial status of strapdown inertial navigation system - Google Patents
Method for determining initial status of strapdown inertial navigation system Download PDFInfo
- Publication number
- CN1908584A CN1908584A CN 200610112526 CN200610112526A CN1908584A CN 1908584 A CN1908584 A CN 1908584A CN 200610112526 CN200610112526 CN 200610112526 CN 200610112526 A CN200610112526 A CN 200610112526A CN 1908584 A CN1908584 A CN 1908584A
- Authority
- CN
- China
- Prior art keywords
- sins
- angle
- gyroscope
- output
- initial
- Prior art date
Links
- 230000001133 acceleration Effects 0.000 claims abstract description 7
- 230000005484 gravity Effects 0.000 claims abstract description 7
- 239000011159 matrix materials Substances 0.000 claims description 24
- 239000000969 carriers Substances 0.000 claims description 17
- 238000001914 filtration Methods 0.000 description 8
- 235000020127 ayran Nutrition 0.000 description 6
- 238000004364 calculation methods Methods 0.000 description 6
- 238000000034 methods Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 238000010586 diagrams Methods 0.000 description 2
- 241000212893 Chelon labrosus Species 0.000 description 1
- 238000006243 chemical reactions Methods 0.000 description 1
- 230000000875 corresponding Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000002360 preparation methods Methods 0.000 description 1
Abstract
Description
Technical field
The present invention relates to definite method of a kind of strapdown inertial navigation system (SINS) initial attitude, can be used for various in the initial attitude of high-precision strapdown inertial navigation system determine, be particularly suitable for not having indexing mechanism when airborne SINS installs, or the indexing mechanism precision is not high, the corner condition of limited.
Background technology
The ultimate principle of strapdown inertial navigation system (SINS) is the mechanics law according to the relative inertness space of newton's proposition, utilize the line motion and the angular motion parameter in gyroscope, accelerometer measures carrier relative inertness space, under given motion starting condition, carry out integral operation by computing machine, position, speed and attitude information are provided continuously, in real time.SINS relies on its own inertial sensitive element fully, do not rely on any external information and measure navigational parameter, therefore, it has good concealment, be not subjected to the weather condition restriction, advantage such as no signal is lost, and is interference-free, be a kind of complete autonomous type, round-the-clock navigational system, be widely used in fields such as Aeronautics and Astronautics, navigation.According to the ultimate principle of SINS, SINS must obtain initial information before navigator fix resolves, comprise initial position, speed and attitude.Initial position and the velocity information of SINS obtain easily than initial attitude, and the precision after initial attitude is determined is exactly the initial precision of SINS when entering the navigation duty.Carry out fast accurate initial attitude when therefore, SINS starts working and determine it is a crucial step.
Existing strapdown inertial navitation system (SINS) initial attitude is determined to be divided into coarse alignment and two stages of fine alignment.Coarse alignment stage is exactly under quiet pedestal condition, with the direct introducing computing machine of the gyroscope of (fixedly two positions or two positions arbitrarily) on single position or two positions and accelerometer output, calculates the initial attitude of carrier.When using the method, usually ignore the error and the external interference factor of gyroscope and accelerometer, yet these factors can cause error, so the initial attitude computational accuracy is not high.Especially, when the gyroscope on adopting fixing two positions and the output of accelerometer are calculated, require SINS around Z axle Rotate 180 degree or 90 degree, this just need be installed in SINS on the indexing mechanism, utilize indexing mechanism to realize the rotation of 180 degree or 90 degree, extremely inconvenient when engineering is used, and also the precision of indexing mechanism is not high, rotational angle also can't accurately be measured, and has reduced the precision that initial attitude is determined.Because the restriction of SINS installation site in the concrete engineering, also exist rotational angle can't satisfy the situation that 180 degree or 90 degree require, at this moment, existing fixed two positions initial attitude determines that method can't use.In addition, when gyroscope on utilizing any two positions and accelerometer output, need to calculate arcsin function, the error of gyroscope and accelerometer itself and measuring error etc. cause the result of calculation instability easily, the imaginary number phenomenon occurs, therefore often need carry out repeatedly the collection and the calculating of gyroscope and accelerometer output on the two positions arbitrarily, result of calculation is selected, after averaging, as the initial attitude of carrier.So, the time that initial attitude is determined will be multiplied, and the precision that initial attitude calculates also can't be guaranteed.
The fine alignment stage is to carry out on the basis of coarse alignment, utilizes the state-space method of modern control theory, and the data of gyroscope and accelerometer output are handled.When the data of single position are handled, the considerable degree of orientation misalignment is poor, speed of convergence is slower, required time is longer, and gyroscope and accelerometer error are unobservable, therefore, can't estimate preferably, do not reach the purpose that improves the initial attitude precision, can not when initial attitude calculates, realize demarcation gyroscope and accelerometer.When the data of fixing two positions are handled, as above to narrate, existing fixed two positions initial attitude determines that there is the restriction of the anglec of rotation in method, and higher to the indexing mechanism accuracy requirement, often can not use in the concrete engineering.
Summary of the invention
Technology of the present invention is dealt with problems and is: overcome the deficiencies in the prior art, provide a kind of accurately, initial status of strapdown inertial navigation system is determined method easily.
Technical solution of the present invention is: a kind of initial status of strapdown inertial navigation system is determined method,
Concrete steps are as follows:
(1) the SINS preheating is prepared, and concrete setup time is according to different systems and difference;
(2) SINS ready after, keep SINS in initial position (being called primary importance) transfixion, gather gyroscope output and accelerometer and export;
(3) utilize accelerometer output and the relation of acceleration of gravity and the relation of gyroscope output and earth rotation angular speed, calculate the course angle of primary importance 1, pitching angle theta 1With roll angle γ 1
(4) need not indexing mechanism, SINS is rotated to any one position (being called the second place) from primary importance around the space three-dimensional arbitrary axis, keep SINS, gather gyroscope output and accelerometer output in second place transfixion by any means;
(5) adopt Kalman Filter Technology, with 1, θ 1And γ 1As initial parameter, the data of gyroscope and accelerometer output on two positions are handled.Because the change of SINS position, changed the system matrix in the SINS error model, thereby the observability of SINS system improves, filter effect is improved, and the error angle that calculates when better estimating the second place between geographic coordinate system n ' and the true geographic coordinate system n (abbreviates misalignment φ as x, φ yAnd φ z) and gyroscope constant value drift, accelerometer often be worth biasing;
(6) utilize the φ that estimates x, φ yAnd φ zCalculate the transition matrix C between true geographic coordinate system n and the calculating geographic coordinate system n ' N ' nAccording to the angle increment or the angular velocity information of gyroscope output, adopt the hypercomplex number method, calculate the transition matrix C between carrier coordinate system b and the calculating geographic coordinate system n ' b N 'By C N ' nAnd C b N ', calculate the transition matrix C between carrier coordinate system b and the true geographic coordinate system n b n, again by C b nCalculate the course angle of the second place 2, pitching angle theta 2With roll angle γ 2, with its initial attitude as SINS.
Principle of the present invention is: when SINS kept transfixion in some positions, accelerometer output and acceleration of gravity and gyroscope output had following relation with the earth rotation angular speed:
Wherein, f
X1, f
Y1And f
Z1And ω
X1, ω
Y1And ω
Z1The respectively specific force and the angular speed of X-axis, Y-axis and the output of Z axle of SINS on the position for this reason; G is an acceleration of gravity; ω
IeBe the earth rotation angular speed, its in east, north, old name for the Arabian countries in the Middle East make progress be projected as Ω=[0 ω
IeCosL ω
IeSinL]; L is a local latitude; C
n bFor navigation coordinate is tied to the transition matrix of carrier coordinate system,
C n bCan be expressed as this position course angle 1, pitching angle theta 1With roll angle γ 1Expression formula, utilize (1) formula and (2) formula can calculate the of this position 1, β 1And γ 1Computing formula is as follows:
By (1), (2) Shi Kede:
C 11=C 22C 33-C 23C 32
C 21=-C 12C 33+C 13C 32
C 31=C 12C 23-C 13C 22
Course angle 1, pitching angle theta 1With roll angle γ 1The main value computing formula be:
If course angle 1, pitching angle theta 1With roll angle γ 1Span be defined as [0,2 π], [π ,+π], so 1, θ 1And γ 1True value can determine as follows.
θ 1=θ 1 is main(5)
By the definite of (5) formula 1, θ 1And γ 1Be SINS at this locational course angle, the angle of pitch and roll angle.
Owing to do not consider the measuring error of accelerometer biasing, gyroscopic drift and accelerometer and gyroscope output information in (1) formula and (2) formula, the course angle that tries to achieve 1, pitching angle theta 1With roll angle γ 1Can not the accurate description carrier coordinate system and local geographic coordinate system between true angular relationship.Therefore, should further utilize modern estimation theory that the initial attitude misalignment is estimated from stochastic error and random disturbance on this basis.
Utilize Kalman Filter Technology, with 1, θ 1And γ 1As initial parameter, can estimate the error angle that calculates between geographic coordinate system n ' and the true geographic coordinate system n, proofread and correct C b N 'After, can obtain initial attitude more accurately.But, when Kalman filtering is handled the accelerometer of single position and gyroscope output, the incomplete may observe of system, wherein, two accelerometers and a gyrostatic error are unobservable, so estimation effect is poor, do not reach the purpose that improves the initial attitude precision.The present invention is by rotating to any one position from primary importance around the space three-dimensional arbitrary axis with SINS, promptly change the position of SINS, change the system matrix in the SINS error model, thereby improve the observability of SINS system, improve the Kalman filtering effect, estimate the error of misalignment and gyroscope, accelerometer better.The misalignment that utilization estimates is to transition matrix C b N 'Proofread and correct, obtain course angle, the angle of pitch and roll angle more accurately.
The present invention's advantage compared with prior art is:
(1) the present invention broken that conventional fixed two positions initial attitude calculates need be by indexing mechanism with the constraints of SINS around Z axle Rotate 180 degree or 90 degree, but need not to pass through indexing mechanism, SINS is rotated to the optional position around arbitrary axis, is a kind of any two-position initial attitude computing method of space three-dimensional.Avoided because of the indexing mechanism precision is not high, the reduction of the fixedly two positions initial attitude computational accuracy that causes has promptly improved the precision that initial attitude calculates.In addition, SINS being rotated to the optional position around arbitrary axis also greatly facilitates in actual application in engineering.
(2) the present invention utilizes the data of gyroscope and accelerometer output on two positions, improve the observability of SINS system, improve the filter effect of Kalman filter, estimate misalignment and gyroscope constant value drift better, accelerometer often is worth biasing, after the misalignment that utilization estimates is proofreaied and correct attitude matrix, obtain course angle, the angle of pitch and roll angle more accurately.
(3) the present invention can finish to survey when initial attitude calculates and float and the calibration task, has not only improved the precision that system's initial attitude calculates, and calibration precision also is improved, and recompenses at the SINS navigational state, can improve navigation and positioning accuracy effectively.
Description of drawings
Fig. 1 determines the process flow diagram of method for initial attitude of the present invention;
Fig. 2 is the synoptic diagram of course angle , pitching angle theta and roll angle γ, Ox among the figure ny nz nBe navigation coordinate system, i.e. the geographical coordinate system in sky, northeast, Ox by bz bBe carrier coordinate system.Wherein, Fig. 2 a represents that from navigation coordinate be Ox ny nz nBe rotated counterclockwise and carrier coordinate system Ox around the zn axle by bz bOverlap, is course angle; Fig. 2 b represents that from navigation coordinate be Ox ny nz nAround x nAxle is rotated counterclockwise θ and carrier coordinate system Ox by bz bOverlap, θ is the angle of pitch; Fig. 2 c represents that from navigation coordinate be Ox ny nz nAround y nAxle is rotated counterclockwise γ and carrier coordinate system Ox by bz bOverlap, γ is roll angle.
Embodiment
As shown in Figure 1, specific implementation method of the present invention is as follows:
1, the preparation of strapdown inertial navigation system
After the SINS start, enter standby condition.
2, primary importance data acquisition
SINS is ready, keeps SINS in initial position (being called primary importance) transfixion, gathers gyroscope output in 5 minutes and accelerometer output, if but the lower proper extension of the precision of SINS is adopted several times.
3, primary importance Attitude Calculation
Navigation coordinate system is taken as the geographical coordinate system in sky, northeast, the course angle of primary importance place 1, pitching angle theta 1With roll angle γ 1Definition shown in Fig. 2 a, Fig. 2 b and Fig. 2 c.
Accelerometer output and acceleration of gravity and gyroscope output have following relation with the earth rotation angular speed:
In the formula, f X1, f Y1And f Z1And ω X1, ω Y1And ω Z1Be respectively the specific force and the angular speed of X-axis, Y-axis and the output of Z axle of SINS on first position; G is an acceleration of gravity; ω IeBe the earth rotation angular speed, its in east, north, old name for the Arabian countries in the Middle East make progress be projected as Ω=[0 ω IeCosL ω IeSinL], L represents local latitude; C n bFor navigation is tied to the transition matrix that carrier is, can be written as:
By (6), (7) Shi Kede:
C 11=C 22C 33-C 23C 32
C 21=-C 12C 33+C 13C 32
C 31=C 12C 23-C 13C 22
Course angle 1, pitching angle theta 1With roll angle γ 1The main value computing formula be:
θ 1 is main=arcsin (C 23) (9)
If course angle 1, pitching angle theta 1With roll angle γ 1Span be defined as [0,2 π], [π ,+π], so 1, θ 1And γ 1True value can determine as follows.
θ 1=θ 1 is main(10)
By the definite of (10) formula 1, θ 1And γ 1Be course angle, the angle of pitch and the roll angle of SINS on primary importance.
4, second place data acquisition
By any means SINS is rotated to any one position (being called the second place) from initial position around the space three-dimensional arbitrary axis, keep SINS, gather gyroscope output in 5 minutes and accelerometer output in second place transfixion.
5, the data of two positions are carried out Filtering Processing
With 1, θ 1And γ 1As initial parameter, utilize Kalman Filter Technology, the data of gyroscope and accelerometer output on two positions are handled, accurately estimate misalignment φ x, φ yAnd φ z(under the geographical coordinate system in sky, northeast is φ E, φ NAnd φ U) and gyroscope constant value drift, accelerometer often be worth biasing.
(1) the SINS initial attitude is determined the foundation of error model
Navigation coordinate system is taken as the geographical coordinate system in sky, northeast, and position and vertical speed error are omitted, and accelerometer and gyrostatic error are considered as setovering at random and add white-noise process, and this moment, the error model of SINS was:
In the formula, subscript E, N, U represent east, north, sky respectively.Earth rotation angular speed ω IeIn east, north, old name for the Arabian countries in the Middle East make progress be projected as Ω=[0 ω IeCosL ω IeSinL]=[0 Ω NΩ U], L represents local latitude; Subscript x, y, z are carrier coordinate system; T Ij(i=1,2,3; J=1,2,3) be attitude matrix C b nIn element,
(2) the SINS initial attitude is determined the foundation of Kalman filter model
For strapdown inertial navitation system (SINS), consider the random deviation of Gyroscope Random Drift Error and accelerometer, equation (11) is modified to following form:
In the formula, W (t) is N (O, white Gaussian noise Q); State variable X
a=[δ V
Eδ V
Nφ
Eφ
Nφ
U]
T, X
b=[
x
yε
xε
yε
z]
TThe random noise state vector
(12) formula is determined the system equation of Kalman filter model for the strapdown inertial navitation system (SINS) initial attitude.For the application card Thalmann filter carries out the optimal estimation of state vector, also need set up the systematic observation equation.Choosing two horizontal velocity errors is observed quantity, and the systematic observation equation of foundation is
Wherein, η (t) is the systematic observation noise vector, is N (O, white Gaussian noise process R).
(3) foundation of Kalman Filtering for Discrete model
According to said system equation and observation equation, it is as follows to set up the Kalman Filtering for Discrete equation.
Filtering equations:
Gain equation:
Prediction error variance equation:
Filtering error variance equation:
P k=[I-K kH k]P k,k-1 (17)
In the formula, Φ K, k-1Be the state-transition matrix (system matrix) of discretize, Q, R are respectively the covariance matrix of system noise and observation noise.
(4) filtering starting condition
X(0)=0 10×1;
The value of the corresponding medium accuracy SINS of P (0), Q and R is as follows:
P(0)=diag{(0.3m/s) 2,(0.3m/s) 2,(30′) 2,(10″) 2,(10″) 2,
(100μg) 2,(100μg) 2,(0.1°) 2,(0.1°) 2,(0.1°) 2};
Q=diag{(100μg) 2,(100μg) 2,(0.1°) 2,(0.1°) 2,(0.1°) 2,0,0,0,0,0};
R=diag{(0.1m/s) 2,(0.1m/s) 2}。
6, computing gyroscope Random Constant Drift and accelerometer value biasing often at random
The that wave filter estimates x, yBe accelerometer value biasing often at random, ε x, ε y, ε zBe the gyroscope Random Constant Drift.
7, the transition matrix C between calculating carrier coordinate system b and the calculating geographic coordinate system n ' b N '
Can utilize the angle increment or the angular velocity information of gyroscope output, adopt the hypercomplex number method to calculate C b N ', calculation procedure is as follows:
(1) utilizes course angle 1, pitching angle theta 1With roll angle γ 1Attitude during first position of initialization is calculated the initial conversion Matrix C b N 'With hypercomplex number q, computing formula is as follows:
Order
Then have:
(2) upgrade hypercomplex number q 0, q 1, q 2, q 3With transition matrix C b N '
Wherein,
Transition matrix C b N 'More new formula as follows:
8, calculate course angle 2, pitching angle theta 2With roll angle γ 2
The φ that utilizes Kalman Filter Estimation to go out E, φ N, φ UCalculate the transition matrix C between true geographic coordinate system n and the calculating geographic coordinate system n ' N ' nAccording to C N ' n(20) C that calculates of formula b N ', calculate the transition matrix C between carrier coordinate system b and the true geographic coordinate system n b n, again by C b nCalculate the course angle of the second place 2, pitching angle theta 2With roll angle γ 2, with its initial attitude as SINS.Concrete calculation procedure is as follows:
(1) the transition matrix C between true geographic coordinate system n of calculating and the calculating geographic coordinate system n ' n N '
(2) the transition matrix C between calculating carrier coordinate system b and the true geographic coordinate system n b n
(3) calculate course angle 2, pitching angle theta 2With roll angle γ 2
Course angle 2, pitching angle theta 2With roll angle γ 2Definition shown in Fig. 2 a, Fig. 2 b and Fig. 2 c.
The C that (22) formula is calculated b nBe designated as
Again because
Therefore, by (23) formula and (24) formula, can determine course angle 2, pitching angle theta 2With roll angle γ 2Main value, promptly
If course angle 2, pitching angle theta 2With roll angle γ 2Span be defined as respectively [0,2 π], [π ,+π].So, 2, θ 2And γ 2True value can determine by following formula:
θ 2=θ 2 is main(26)
By the definite of (26) formula 2, θ 2And γ 2Be course angle, the angle of pitch and the roll angle of SINS on the second place, it is entered the initial attitude of navigation duty as SINS.
Claims (1)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2006101125264A CN100516775C (en) | 2006-08-23 | 2006-08-23 | Method for determining initial status of strapdown inertial navigation system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2006101125264A CN100516775C (en) | 2006-08-23 | 2006-08-23 | Method for determining initial status of strapdown inertial navigation system |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1908584A true CN1908584A (en) | 2007-02-07 |
CN100516775C CN100516775C (en) | 2009-07-22 |
Family
ID=37699758
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2006101125264A CN100516775C (en) | 2006-08-23 | 2006-08-23 | Method for determining initial status of strapdown inertial navigation system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100516775C (en) |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1983304A2 (en) * | 2007-04-16 | 2008-10-22 | Honeywell International Inc. | Heading stabilization for aided inertial navigation systems |
CN101256080B (en) * | 2008-04-09 | 2010-06-23 | 南京航空航天大学 | Midair aligning method for satellite/inertia combined navigation system |
CN101270993B (en) * | 2007-12-12 | 2011-08-31 | 北京航空航天大学 | Remote high-precision independent combined navigation locating method |
CN101561292B (en) * | 2009-05-12 | 2011-11-16 | 北京航空航天大学 | Method and device for calibrating size effect error of accelerometer |
CN102486377A (en) * | 2009-11-17 | 2012-06-06 | 哈尔滨工程大学 | Method for acquiring initial course attitude of fiber optic gyro strapdown inertial navigation system |
CN102564421A (en) * | 2012-02-09 | 2012-07-11 | 北京机械设备研究所 | Method for prolonging service life of inertia measurement unit |
CN102589546A (en) * | 2012-03-02 | 2012-07-18 | 北京航空航天大学 | Optical-fiber strap-down inertial measurement unit reciprocating-type two-position north finding method for inhibiting slope error influence of devices |
CN102865881A (en) * | 2012-03-06 | 2013-01-09 | 武汉大学 | Quick calibration method for inertial measurement unit |
CN102933937A (en) * | 2010-06-10 | 2013-02-13 | 高通股份有限公司 | Use of inertial sensor data to improve mobile station positioning |
CN102937450A (en) * | 2012-10-31 | 2013-02-20 | 北京控制工程研究所 | Relative attitude determining method based on gyroscope metrical information |
CN103822633A (en) * | 2014-02-11 | 2014-05-28 | 哈尔滨工程大学 | Low-cost attitude estimation method based on second-order measurement update |
EP2239539A3 (en) * | 2009-04-06 | 2014-09-03 | Honeywell International Inc. | Technique to improve navigation performance through carouselling |
CN104919277A (en) * | 2012-12-20 | 2015-09-16 | 大陆-特韦斯贸易合伙股份公司及两合公司 | Method for determining initial data for determining position data of a vehicle |
CN106123921A (en) * | 2016-07-10 | 2016-11-16 | 北京工业大学 | Latitude the unknown Alignment Method of SINS under the conditions of dynamic disturbance |
CN103954298B (en) * | 2014-04-18 | 2017-03-01 | 中国人民解放军国防科学技术大学 | Microthrust test g sensitivity error compensating method in a kind of GNSS and MIMU integrated navigation |
CN109163721A (en) * | 2018-09-18 | 2019-01-08 | 河北美泰电子科技有限公司 | Attitude measurement method and terminal device |
CN110006456A (en) * | 2019-04-24 | 2019-07-12 | 北京星网宇达科技股份有限公司 | A kind of detection vehicle alignment methods, device and equipment |
CN110095115A (en) * | 2019-06-04 | 2019-08-06 | 中国科学院合肥物质科学研究院 | A kind of carrier navigation attitude measurement method updated based on Geomagnetism Information |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101900572B (en) * | 2010-07-09 | 2012-01-04 | 哈尔滨工程大学 | Rapid measuring method for installation error of strapdown inertial system gyroscope based on three-axle rotary table |
-
2006
- 2006-08-23 CN CNB2006101125264A patent/CN100516775C/en active IP Right Grant
Cited By (28)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1983304A3 (en) * | 2007-04-16 | 2014-08-27 | Honeywell International Inc. | Heading stabilization for aided inertial navigation systems |
EP1983304A2 (en) * | 2007-04-16 | 2008-10-22 | Honeywell International Inc. | Heading stabilization for aided inertial navigation systems |
CN101270993B (en) * | 2007-12-12 | 2011-08-31 | 北京航空航天大学 | Remote high-precision independent combined navigation locating method |
CN101256080B (en) * | 2008-04-09 | 2010-06-23 | 南京航空航天大学 | Midair aligning method for satellite/inertia combined navigation system |
US9599474B2 (en) | 2009-04-06 | 2017-03-21 | Honeywell International Inc. | Technique to improve navigation performance through carouselling |
EP2239539A3 (en) * | 2009-04-06 | 2014-09-03 | Honeywell International Inc. | Technique to improve navigation performance through carouselling |
CN101561292B (en) * | 2009-05-12 | 2011-11-16 | 北京航空航天大学 | Method and device for calibrating size effect error of accelerometer |
CN102486377A (en) * | 2009-11-17 | 2012-06-06 | 哈尔滨工程大学 | Method for acquiring initial course attitude of fiber optic gyro strapdown inertial navigation system |
CN102486377B (en) * | 2009-11-17 | 2014-10-22 | 哈尔滨工程大学 | Method for acquiring initial course attitude of fiber optic gyro strapdown inertial navigation system |
CN102933937A (en) * | 2010-06-10 | 2013-02-13 | 高通股份有限公司 | Use of inertial sensor data to improve mobile station positioning |
CN102564421A (en) * | 2012-02-09 | 2012-07-11 | 北京机械设备研究所 | Method for prolonging service life of inertia measurement unit |
CN102589546B (en) * | 2012-03-02 | 2014-09-03 | 北京航空航天大学 | Optical-fiber strap-down inertial measurement unit reciprocating-type two-position north finding method for inhibiting slope error influence of devices |
CN102589546A (en) * | 2012-03-02 | 2012-07-18 | 北京航空航天大学 | Optical-fiber strap-down inertial measurement unit reciprocating-type two-position north finding method for inhibiting slope error influence of devices |
WO2013131471A1 (en) * | 2012-03-06 | 2013-09-12 | 武汉大学 | Quick calibration method for inertial measurement unit |
CN102865881A (en) * | 2012-03-06 | 2013-01-09 | 武汉大学 | Quick calibration method for inertial measurement unit |
CN102865881B (en) * | 2012-03-06 | 2014-12-31 | 武汉大学 | Quick calibration method for inertial measurement unit |
CN102937450B (en) * | 2012-10-31 | 2015-11-25 | 北京控制工程研究所 | A kind of relative attitude defining method based on gyro to measure information |
CN102937450A (en) * | 2012-10-31 | 2013-02-20 | 北京控制工程研究所 | Relative attitude determining method based on gyroscope metrical information |
CN104919277A (en) * | 2012-12-20 | 2015-09-16 | 大陆-特韦斯贸易合伙股份公司及两合公司 | Method for determining initial data for determining position data of a vehicle |
CN103822633B (en) * | 2014-02-11 | 2016-12-07 | 哈尔滨工程大学 | A kind of low cost Attitude estimation method measuring renewal based on second order |
CN103822633A (en) * | 2014-02-11 | 2014-05-28 | 哈尔滨工程大学 | Low-cost attitude estimation method based on second-order measurement update |
CN103954298B (en) * | 2014-04-18 | 2017-03-01 | 中国人民解放军国防科学技术大学 | Microthrust test g sensitivity error compensating method in a kind of GNSS and MIMU integrated navigation |
CN106123921A (en) * | 2016-07-10 | 2016-11-16 | 北京工业大学 | Latitude the unknown Alignment Method of SINS under the conditions of dynamic disturbance |
CN106123921B (en) * | 2016-07-10 | 2019-05-24 | 北京工业大学 | The unknown Alignment Method of the latitude of Strapdown Inertial Navigation System under the conditions of dynamic disturbance |
CN109163721A (en) * | 2018-09-18 | 2019-01-08 | 河北美泰电子科技有限公司 | Attitude measurement method and terminal device |
CN110006456A (en) * | 2019-04-24 | 2019-07-12 | 北京星网宇达科技股份有限公司 | A kind of detection vehicle alignment methods, device and equipment |
CN110095115A (en) * | 2019-06-04 | 2019-08-06 | 中国科学院合肥物质科学研究院 | A kind of carrier navigation attitude measurement method updated based on Geomagnetism Information |
CN110095115B (en) * | 2019-06-04 | 2020-12-25 | 中国科学院合肥物质科学研究院 | Carrier attitude and heading measurement method based on geomagnetic information update |
Also Published As
Publication number | Publication date |
---|---|
CN100516775C (en) | 2009-07-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Li et al. | LIDAR/MEMS IMU integrated navigation (SLAM) method for a small UAV in indoor environments | |
CN106780699B (en) | Visual SLAM method based on SINS/GPS and odometer assistance | |
CN104655131B (en) | Inertial navigation Initial Alignment Method based on ISTSSRCKF | |
US9715234B2 (en) | Multiple rotors aircraft and control method | |
WO2017063387A1 (en) | Method for updating all attitude angles of agricultural machine on the basis of nine-axis mems sensor | |
CN104197927B (en) | Submerged structure detects robot real-time navigation system and method | |
CN105823480B (en) | Underwater moving target location algorithm based on single beacon | |
CN104457754B (en) | SINS/LBL (strapdown inertial navigation systems/long base line) tight combination based AUV (autonomous underwater vehicle) underwater navigation positioning method | |
CN1090314C (en) | Movement detector | |
WO2018184467A1 (en) | Method and device for detecting posture of ball head | |
EP3089267B1 (en) | Method and system for controlling antenna of mobile communication application system based on double quaternions in mems inertial navigation | |
CN106289246B (en) | A kind of flexible link arm measure method based on position and orientation measurement system | |
CN103913181B (en) | A kind of airborne distributed POS Transfer Alignments based on parameter identification | |
US20030164053A1 (en) | Apparatus and method for accurate pipeline surveying | |
US7248967B2 (en) | Autonomous velocity estimation and navigation | |
JP2005115623A (en) | Navigation system using image recognition | |
Ishibashi et al. | The rotation control system to improve the accuracy of an inertial navigation system installed in an autonomous underwater vehicle | |
JP2009019992A (en) | Position detection device and position detection method | |
CN100344938C (en) | Apparatus and method for detecting movable object position in navigation system | |
CN103245360B (en) | Carrier-borne aircraft rotation type strapdown inertial navigation system Alignment Method under swaying base | |
CN103235328B (en) | GNSS (global navigation satellite system) and MEMS (micro-electromechanical systems) integrated navigation method | |
Bryne et al. | Nonlinear observers for integrated INS\/GNSS navigation: implementation aspects | |
CN100476360C (en) | Integrated navigation method based on star sensor calibration | |
Dukan et al. | Dynamic positioning system for a small size ROV with experimental results | |
EP1668381A4 (en) | Inertial gps navigation system using injected alignment data for the inertial system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
C06 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
C10 | Entry into substantive examination | ||
GR01 | Patent grant | ||
C14 | Grant of patent or utility model |