CN102788598A - Error suppressing method of fiber strap-down inertial navigation system based on three-axis rotation - Google Patents

Error suppressing method of fiber strap-down inertial navigation system based on three-axis rotation Download PDF

Info

Publication number
CN102788598A
CN102788598A CN201210305216XA CN201210305216A CN102788598A CN 102788598 A CN102788598 A CN 102788598A CN 201210305216X A CN201210305216X A CN 201210305216XA CN 201210305216 A CN201210305216 A CN 201210305216A CN 102788598 A CN102788598 A CN 102788598A
Authority
CN
China
Prior art keywords
integral
epsiv
rightarrow
imu
order
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
Application number
CN201210305216XA
Other languages
Chinese (zh)
Other versions
CN102788598B (en
Inventor
孙伟
徐爱功
徐宗秋
车莉娜
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Liaoning Technical University
Original Assignee
Liaoning Technical University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Liaoning Technical University filed Critical Liaoning Technical University
Priority to CN201210305216.XA priority Critical patent/CN102788598B/en
Publication of CN102788598A publication Critical patent/CN102788598A/en
Application granted granted Critical
Publication of CN102788598B publication Critical patent/CN102788598B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Navigation (AREA)
  • Gyroscopes (AREA)

Abstract

The invention provides an error suppressing method of a fiber strap-down inertial navigation system based on three-axis rotation. The error suppressing method comprises the following steps: determining an initial position parameter of a carrier by using a GPS (global positioning system); collecting data output by a fiber-optic gyroscope and a quartz accelerometer; determining attitude information of the carrier according to relationship of output of the accelerometer and gravity acceleration as well as relationship of output of the gyroscope and the rotational angular velocity of the earth, and accomplishing initial alignment of the system; employing an indexing scheme of a rotational period of twelve rotation-stop orders by IMU (inertial measurement unit); converting data generated by the fiber-optic gyroscope and the quartz accelerometer after rotation of the IMU to be in navigation coordinate system to obtain a modulation format of constant deviation of an inertial unit; updating a strap-down matrix by utilizing an output value of the fiber-optic gyroscope; and calculating position information of the carrier after the IMU is rotated and modulated. With the adoption of the error suppressing method disclosed by the invention, constant deviation of the inertial unit in three-axis direction can be modulated, thereby navigation positioning precision is improved.

Description

Fiber strapdown inertial navigation system error inhibition method based on three rotations
(1) technical field
What the present invention relates to is a kind of measuring method, in particular a kind of fiber strapdown inertial navigation system error inhibition method based on three rotations.
(2) background technology
Inertial navigation is to utilize inertia sensitive element (gyroscope and accelerometer) to measure the line motion and the angular motion in carrier relative inertness space, and under known starting condition, goes out the navigational parameters such as speed, position and attitude of carrier with COMPUTER CALCULATION.It relies on the Sensitive Apparatus of self to accomplish navigation task fully, need not to rely on any external information, also not to any energy of external radiation, is a kind of autonomous navigational system fully, therefore has good concealment, anti-interference, the advantage that do not receive any meteorological condition restriction.In addition, inertial navigation system also has the characteristics of data updating rate height, short-term accuracy height and good stability.Just because of above advantage, it has obtained widespread use at Aeronautics and Astronautics, navigation and a lot of civil area.In strapdown inertial navigation system; All inertance elements are directly installed on the carrier; What inertance element was exported is exactly angular velocity and the acceleration of carrier with respect to inertial space; By computing machine the acceleration information that records under the carrier coordinate system is transformed into navigation coordinate system and carries out navigation calculation again, be equivalent to utilize the gyroscope output data in computing machine, to make up of the reference of a mathematical platform as navigation calculating.
Optical fibre gyro is as a kind of novel angular rate sensor; Compare with traditional gyroscope (liquid floated gyroscope, dynamically tuned gyro, DTG, electrostatic gyro); Has significant advantage: 1) owing to have no rotary part; Thereby firm in structure, anti-vibration, shock resistance, anti-big overload, reliability is high.Simultaneity factor is low in energy consumption, does not need preheating, and start-up time is short, and need not keep in repair, and the life-span is long; 2) because optical fiber is nonmetallic materials, so radiation resistance, strong interference immunity, stable performance can work in the comparatively abominable electromagnetic environment; 3) owing to the area of sensitivity with fiber optic loop is directly proportional, can increase the area of fiber optic loop through the way that increases the fiber optic loop number of turns, the sensitivity that improves gyro, so volume is little, simple in structure, processing technology is simple and cost is low; 4) dynamic range is big, the latch-up phenomenon in the time of can not low rate occurring, and can directly export digital signal, be convenient to utilize computing machine to carry out system in combination.
In the strapdown inertial navigation system, people have promoted the fast development of inertia device for the lasting research of inertia devices such as gyroscope that constitutes Inertial Measurement Unit and accelerometer.But the device precision is high more, and further the cost of boost device precision is just big more.After the inertia device precision reached certain requirement, the performance that adopts the next further improvement of the method system of compensation inertia device deviation was to realize a more realization approach of high precision navigation.The compensation method of inertance element has two kinds: the one, and utilize external information to compensate correction; Another kind method is the self compensation of inertia device deviation; The rotation modulation technique is a kind of method of self compensation; Through around an axle or a plurality of rotator inertia measuring units (IMU), navigation error is modulated, reach the purpose that navigation accuracy is dispersed, improved to the control navigation error.
The single shaft rotation only can compensate the normal value deviation of inertia device on two sensitive axes directions; Though the twin shaft rotation can compensate the normal value deviation of inertia device on three sensitive axes directions, can't avoid the negative effect of carrier angular motion to the rotation modulation technique.Therefore, how reasonable in design three rotation compensation modes have important meaning for the navigation accuracy of further raising fiber strapdown inertial navigation system.
(3) summary of the invention
Technology of the present invention is dealt with problems and is: overcome the deficiency of prior art, a kind of fiber strapdown inertial navigation system error inhibition method based on three sensitive axes rotations of Inertial Measurement Unit is provided.
Technical solution of the present invention is: a kind of fiber strapdown inertial navigation system error inhibition method based on three rotations; It is characterized in that adopting three transposition schemes of Inertial Measurement Unit to come to isolate fully the carrier angular motion; Make the relative geographic coordinate system of Inertial Measurement Unit static; Avoid the carrier angular motion for the negative influence of adopting Inertial Measurement Unit rotation modulation technique, can confirm that inertia device often is worth the inhibition form of deviation, to realize more precise navigation.Its concrete steps are following:
(1) confirms the initial position parameters of carrier through GPS, they are bound to navigational computer;
(2) SINS carries out preheating and prepares, and gathers data that fibre optic gyroscope and quartz accelerometer export and data are handled;
(3) IMU adopts 12 commentaries on classics to stop the transposition scheme that order is a swing circle (like accompanying drawing 2);
Order 1, IMU rotates counterclockwise 180 ° of in-position B, stand-by time T from the A point sOrder 2, IMU rotates counterclockwise 180 ° of in-position C, stand-by time T from the B point sOrder 3, IMU rotates counterclockwise 180 ° of in-position A, stand-by time T from the C point sOrder 4, IMU rotates counterclockwise 180 ° of in-position C, stand-by time T from the A point sOrder 5, IMU rotates counterclockwise 180 ° of in-position B, stand-by time T from the C point sOrder 6, IMU rotates counterclockwise 180 ° of in-position A, stand-by time T from the B point sOrder 7, IMU clockwise rotates 180 ° of in-position B, stand-by time T from the A point sOrder 8, IMU clockwise rotates 180 ° of in-position C, stand-by time T from the B point sOrder 9, IMU clockwise rotates 180 ° of in-position A, stand-by time T from the C point sOrder 10, IMU clockwise rotates 180 ° of in-position C, stand-by time T from the A point sOrder 11, IMU clockwise rotates 180 ° of in-position B, stand-by time T from the C point sOrder 12, IMU clockwise rotates 180 ° of in-position A, stand-by time T from the B point sIMU rotates sequential loop according to this to carry out.
(4) data-switching that Inertial Measurement Unit rotation back gyroscope is generated obtains the modulation format that inertia device often is worth deviation under carrier coordinate system;
Suppose that the gyroscope constant value drift on the IMU horizontal direction is respectively ε xAnd ε yUnder the carrier quiescent conditions; Because three positions of A, B, C that IMU pauses are with respect to navigation coordinate system symmetry; Therefore on three fixed positions in one three transposition cycles, three gyroscope constant value drifts are fastened the attitude error that projection causes at navigation coordinate and must be satisfied:
3 ( ∫ 0 T s ϵ x n dt ) A + 3 ( ∫ 0 T s ϵ x n dt ) B + 3 ( ∫ 0 T s ϵ x n dt ) C = 0
3 ( ∫ 0 T s ϵ y n dt ) A + 3 ( ∫ 0 T s ϵ y n dt ) B + 3 ( ∫ 0 T s ϵ y n dt ) C = 0
3 ( ∫ 0 T s ϵ z n dt ) A + 3 ( ∫ 0 T s ϵ z n dt ) B + 3 ( ∫ 0 T s ϵ z n dt ) C = 0
Exist the symmetry problem of rotation according to the rotation in three scheme of rotation of IMU, ignore the influence of carrier movement and with local geographic coordinate system as a reference, 12 order transposition schemes can be expressed as:
Process 1: order 1,6,7,12, in the rotation period of formation, the gyroscopic drift of x, y axle is ox at navigation coordinate ny nDemonstrate the Changing Pattern in positive and negative each week in the plane, the normal value deviation that therefore in the integral process of complete cycle, produces is zero, that is:
( ∫ 0 T z ϵ x n dt ) A → + B + ( ∫ 0 T z ϵ x n dt ) B → + A + ( ∫ 0 T z ϵ x n dt ) A → - B + ( ∫ 0 T z ϵ x n dt ) B → - A = 0
( ∫ 0 T z ϵ y n dt ) A → + B + ( ∫ 0 T z ϵ y n dt ) B → + A + ( ∫ 0 T z ϵ y n dt ) A → - B + ( ∫ 0 T z ϵ y n dt ) B → - A = 0
Wherein, the time of each rotation process is counted T z, around the responsive coordinate axis of Inertial Measurement Unit rotate counterclockwise into+, clockwise rotate into-.
Process 2: order 2,5,8,11, in the rotation period of formation, the gyroscopic drift of y, z axle is oy at navigation coordinate nz nDemonstrate the Changing Pattern in positive and negative each week in the plane, the normal value deviation that therefore in the integral process of complete cycle, produces is zero, that is:
( ∫ 0 T z ϵ y n dt ) B → + C + ( ∫ 0 T z ϵ y n dt ) C → + B + ( ∫ 0 T z ϵ y n dt ) B → - C + ( ∫ 0 T z ϵ y n dt ) C → - B = 0
( ∫ 0 T z ϵ z n dt ) B → + C + ( ∫ 0 T z ϵ z n dt ) C → + B + ( ∫ 0 T z ϵ z n dt ) B → - C + ( ∫ 0 T z ϵ z n dt ) C → - B = 0
Process 3: order 3,4,9,10, in the rotation period of formation, the gyroscopic drift of x, z axle is ox at navigation coordinate nz nDemonstrate the Changing Pattern in positive and negative each week in the plane, the normal value deviation that therefore in the integral process of complete cycle, produces is zero, that is:
( ∫ 0 T z ϵ x n dt ) C → + A + ( ∫ 0 T z ϵ x n dt ) A → + C + ( ∫ 0 T z ϵ x n dt ) C → - A + ( ∫ 0 T z ϵ x n dt ) A → - C = 0
( ∫ 0 T z ϵ z n dt ) C → + A + ( ∫ 0 T z ϵ z n dt ) A → + C + ( ∫ 0 T z ϵ z n dt ) C → - A + ( ∫ 0 T z ϵ z n dt ) A → - C = 0
It is exactly the value that periodically changes the strapdown matrix that 12 order change the process of stopping, make three gyrostatic sensitive axes in a rotation period along center of rotation be symmetrically distributed (like accompanying drawing 3).Having proved intuitively that one 12 order changes stops in the process, and gyroscope often is worth the relative navigation coordinate of deviation system by modulation fully, and the navigation accuracy of system is not exerted an influence.In like manner stop in the cycle a complete commentaries on classics because the symmetrical distribution of three fixed positions and rotation process, can obtain that Inertial Measurement Unit stops and the transposition process in the accelerometer zero drift in the similar effect effect of navigation coordinate system.
(5) output valve of gyroscope under the IMU coordinate system brought in the strapdown inertial navigation system, adopts the equivalent rotating vector method that strapdown matrix
Figure BSA00000768634100041
is upgraded:
ω ns s = ω is s - ( C s n ) T ( ω ie n + ω en n )
Wherein:
Figure BSA00000768634100043
is the component of rotational-angular velocity of the earth under navigation system;
Figure BSA00000768634100044
is the component of motion angular velocity under navigation system of spherical coordinate system relatively for navigation coordinate;
Figure BSA00000768634100045
is the component of motion angular velocity on the IMU coordinate system of the relative navigation coordinate of IMU system.
If the equivalent rotating vector differential equation of the relative navigation coordinate of IMU coordinate system system is:
Φ · = ω ns s + 1 2 Φ × ω ns s + 1 12 Φ × ( Φ × ω ns s )
Solve the rotating vector of equivalence and replace hypercomplex number to separate according to angular velocity
Figure BSA00000768634100047
q = cos Φ 2 + Φ | Φ | sin Φ 2
Because q=q 0+ q 1I+q 2J+q 3K, i, j, k are direction vector.Therefore the renewal process of attitude matrix
Figure BSA00000768634100049
is:
C s n = q 0 2 + q 1 2 - q 2 2 - q 3 2 2 ( q 1 q 2 - q 0 q 3 ) 2 ( q 1 q 3 + q 0 q 2 ) 2 ( q 1 q 2 + q 0 q 3 ) q 0 2 - q 1 2 + q 2 2 - q 3 2 2 ( q 2 q 3 - q 0 q 1 ) 2 ( q 1 q 3 - q 0 q 2 ) 2 ( q 2 q 3 + q 0 q 1 ) q 0 2 - q 1 2 - q 2 2 + q 3 2
(6) the use of quartz accelerometer output value
Figure BSA000007686341000411
and Step (5) the calculated attitude matrix
Figure BSA000007686341000412
calculated after IMU rotation modulated carrier position.
1) calculates the acceleration f down of navigation system n:
f n = C s n f is s
2) position of calculating carrier:
According to t 1Carrier east orientation horizontal velocity V constantly x(t 1) and north orientation horizontal velocity V y(t 1), ask for t 2Carrier positions is constantly:
Figure BSA000007686341000414
3) calculate the carrier positions error:
Figure BSA000007686341000415
Wherein:
Figure BSA00000768634100051
λ 0Longitude and the latitude of representing initial time carrier present position respectively; Δ λ representes the latitude of carrier, the variable quantity of longitude respectively; R N, R MThe radius-of-curvature of representing earth meridian circle, prime vertical respectively; t 1, t 2Two the adjacent time points in the process that resolve for inertial navigation system.
The present invention's advantage compared with prior art is: the present invention has broken in traditional SINS IMU and carrier and has been connected and causes system's navigation accuracy to receive the constraint of inertia device effects; Propose a kind of IMU often to be worth the deviation modulation scheme around the inertia device that three fixing position rotatings of the sensitive axes of three directions of carrier stop; This method can often be worth deviation with all inertia devices modulates, and improves navigation and positioning accuracy effectively.
Effect to the present invention is useful is explained as follows:
Under the VC++ simulated conditions, this method is carried out emulation experiment:
Carrier remains static, the error model parameters of IMU three positions 12 order rotation-stop schemes:
The dead time of three positions: T s=5 minutes;
The time that consumes when rotating 180 °: T z=12 seconds;
Rotate in 180 ° the process, the acceleration and deceleration time in each transposition respectively is 4 seconds;
Carrier initial position: 45.7796 ° of north latitude, 126.6705 ° of east longitudes;
The initial attitude error angle: three initial attitude error angles are zero;
Equatorial radius: R e=6378393.0 meters;
Ellipsoid degree: e=3.367e-3;
The earth surface acceleration of gravity that can get by universal gravitation: g 0=9.78049;
Rotational-angular velocity of the earth (radian per second): 7.2921158e-5;
The gyroscope constant value drift: 0.01 degree/hour;
Accelerometer bias: 10 -4g 0
Constant: π=3.1415926;
It is as shown in Figure 4 to utilize the said method of invention to obtain the carrier positions graph of errors.The result shows that the 12 order commentaries on classics of IMU three positions stops under the condition, adopts the inventive method can obtain high orientation precision.
(4) description of drawings
Fig. 1 is a kind of fiber strapdown inertial navigation system error inhibition method process flow diagram based on three rotations of the present invention;
Fig. 2 is the fiber strapdown inertial navigation system IMU rotation-stop scheme detailed step figure based on three rotations of the present invention;
Fig. 3 is that the fiber strapdown inertial navigation system IMU based on three rotations of the present invention changes stopping time constant value drift orientation distribution;
Contrast experiment's curve of carrier positioning error when Fig. 4 is carrier positions error and the IMU stationary state of the fiber strapdown inertial navigation system based on three rotations of the present invention.
(5) embodiment
Describe in detail below in conjunction with the accompanying drawing specific embodiments of the invention:
(1) confirms the initial position parameters of carrier through GPS, they are bound to navigational computer;
(2) SINS carries out preheating and prepares, and gathers data that fibre optic gyroscope and quartz accelerometer export and data are handled;
(3) IMU adopts 12 commentaries on classics to stop the transposition scheme that order is a swing circle (like accompanying drawing 2);
Order 1, IMU rotates counterclockwise 180 ° of in-position B, stand-by time T from the A point sOrder 2, IMU rotates counterclockwise 180 ° of in-position C, stand-by time T from the B point sOrder 3, IMU rotates counterclockwise 180 ° of in-position A, stand-by time T from the C point sOrder 4, IMU rotates counterclockwise 180 ° of in-position C, stand-by time T from the A point sOrder 5, IMU rotates counterclockwise 180 ° of in-position B, stand-by time T from the C point sOrder 6, IMU rotates counterclockwise 180 ° of in-position A, stand-by time T from the B point sOrder 7, IMU clockwise rotates 180 ° of in-position B, stand-by time T from the A point sOrder 8, IMU clockwise rotates 180 ° of in-position C, stand-by time T from the B point sOrder 9, IMU clockwise rotates 180 ° of in-position A, stand-by time T from the C point sOrder 10, IMU clockwise rotates 180 ° of in-position C, stand-by time T from the A point sOrder 11, IMU clockwise rotates 180 ° of in-position B, stand-by time T from the C point sOrder 12, IMU clockwise rotates 180 ° of in-position A, stand-by time T from the B point sIMU rotates sequential loop according to this to carry out.
(4) data-switching that Inertial Measurement Unit rotation back gyroscope is generated obtains the modulation format that inertia device often is worth deviation under carrier coordinate system;
Suppose that the gyroscope constant value drift on the IMU horizontal direction is respectively ε xAnd ε yUnder the carrier quiescent conditions; Because three positions of A, B, C that IMU pauses are with respect to navigation coordinate system symmetry; Therefore on three fixed positions in one three transposition cycles, three gyroscope constant value drifts are fastened the attitude error that projection
Figure BSA00000768634100061
causes at navigation coordinate and must be satisfied:
3 ( ∫ 0 T s ϵ x n dt ) A + 3 ( ∫ 0 T s ϵ x n dt ) B + 3 ( ∫ 0 T s ϵ x n dt ) C = 0
3 ( ∫ 0 T s ϵ y n dt ) A + 3 ( ∫ 0 T s ϵ y n dt ) B + 3 ( ∫ 0 T s ϵ y n dt ) C = 0 - - - ( 1 )
3 ( ∫ 0 T s ϵ z n dt ) A + 3 ( ∫ 0 T s ϵ z n dt ) B + 3 ( ∫ 0 T s ϵ z n dt ) C = 0
Exist the symmetry problem of rotation according to the rotation in three scheme of rotation of IMU, ignore the influence of carrier movement and with local geographic coordinate system as a reference, 12 order transposition schemes can be expressed as:
Process 1: order 1,6,7,12, in the rotation period of formation, the gyroscopic drift of x, y axle is ox at navigation coordinate ny nDemonstrate the Changing Pattern in positive and negative each week in the plane, the normal value deviation that therefore in the integral process of complete cycle, produces is zero, that is:
( ∫ 0 T z ϵ x n dt ) A → + B + ( ∫ 0 T z ϵ x n dt ) B → + A + ( ∫ 0 T z ϵ x n dt ) A → - B + ( ∫ 0 T z ϵ x n dt ) B → - A = 0
( 2 )
( ∫ 0 T z ϵ y n dt ) A → + B + ( ∫ 0 T z ϵ y n dt ) B → + A + ( ∫ 0 T z ϵ y n dt ) A → - B + ( ∫ 0 T z ϵ y n dt ) B → - A = 0
Wherein, the time of each rotation process is counted T z, around the responsive coordinate axis of Inertial Measurement Unit rotate counterclockwise into+, clockwise rotate into-.
Process 2: order 2,5,8,11, in the rotation period of formation, the gyroscopic drift of y, z axle is oy at navigation coordinate nz nDemonstrate the Changing Pattern in positive and negative each week in the plane, the normal value deviation that therefore in the integral process of complete cycle, produces is zero, that is:
( ∫ 0 T z ϵ y n dt ) B → + C + ( ∫ 0 T z ϵ y n dt ) C → + B + ( ∫ 0 T z ϵ y n dt ) B → - C + ( ∫ 0 T z ϵ y n dt ) C → - B = 0
( 3 )
( ∫ 0 T z ϵ z n dt ) B → + C + ( ∫ 0 T z ϵ z n dt ) C → + B + ( ∫ 0 T z ϵ z n dt ) B → - C + ( ∫ 0 T z ϵ z n dt ) C → - B = 0
Process 3: order 3,4,9,10, in the rotation period of formation, the gyroscopic drift of x, z axle is ox at navigation coordinate nz nDemonstrate the Changing Pattern in positive and negative each week in the plane, the normal value deviation that therefore in the integral process of complete cycle, produces is zero, that is:
( ∫ 0 T z ϵ x n dt ) C → + A + ( ∫ 0 T z ϵ x n dt ) A → + C + ( ∫ 0 T z ϵ x n dt ) C → - A + ( ∫ 0 T z ϵ x n dt ) A → - C = 0
( 4 )
( ∫ 0 T z ϵ z n dt ) C → + A + ( ∫ 0 T z ϵ z n dt ) A → + C + ( ∫ 0 T z ϵ z n dt ) C → - A + ( ∫ 0 T z ϵ z n dt ) A → - C = 0
It is exactly the value that periodically changes the strapdown matrix that 12 order change the process of stopping, make three gyrostatic sensitive axes in a rotation period along center of rotation be symmetrically distributed (like accompanying drawing 3).Having proved intuitively that one 12 order changes stops in the process, and gyroscope often is worth the relative navigation coordinate of deviation system by modulation fully, and the navigation accuracy of system is not exerted an influence.In like manner stop in the cycle a complete commentaries on classics because the symmetrical distribution of three fixed positions and rotation process, can obtain that Inertial Measurement Unit stops and the transposition process in the accelerometer zero drift in the similar effect effect of navigation coordinate system.
(5) output valve of gyroscope under the IMU coordinate system brought in the strapdown inertial navigation system, adopts the equivalent rotating vector method that strapdown matrix
Figure BSA000007686341000711
is upgraded:
ω ns s = ω is s - ( C s n ) T ( ω ie n + ω en n ) - - - ( 5 )
Wherein:
Figure BSA000007686341000713
is the component of rotational-angular velocity of the earth under navigation system;
Figure BSA000007686341000714
is the component of motion angular velocity under navigation system of spherical coordinate system relatively for navigation coordinate;
Figure BSA000007686341000715
is the component of motion angular velocity on the IMU coordinate system of the relative navigation coordinate of IMU system.
If the equivalent rotating vector differential equation of the relative navigation coordinate of IMU coordinate system system is:
Φ · = ω ns s + 1 2 Φ × ω ns s + 1 12 Φ × ( Φ × ω ns s ) - - - ( 6 )
Solve the rotating vector of equivalence and replace hypercomplex number to separate according to angular velocity
Figure BSA00000768634100082
q = cos Φ 2 + Φ | Φ | sin Φ 2 - - - ( 7 )
Because q=q 0+ q 1I+q 2J+q 3K, i, j, k are direction vector.Therefore the renewal process of attitude matrix is:
C s n = q 0 2 + q 1 2 - q 2 2 - q 3 2 2 ( q 1 q 2 - q 0 q 3 ) 2 ( q 1 q 3 + q 0 q 2 ) 2 ( q 1 q 2 + q 0 q 3 ) q 0 2 - q 1 2 + q 2 2 - q 3 2 2 ( q 2 q 3 - q 0 q 1 ) 2 ( q 1 q 3 - q 0 q 2 ) 2 ( q 2 q 3 + q 0 q 1 ) q 0 2 - q 1 2 - q 2 2 + q 3 2 - - - ( 8 )
(6) the use of quartz accelerometer output value
Figure BSA00000768634100086
and Step (5) the calculated attitude matrix
Figure BSA00000768634100087
calculated after IMU rotation modulated carrier position.
1) calculates the acceleration f down of navigation system n:
f n = C s n f is s - - - ( 9 )
2) position of calculating carrier:
According to t 1Carrier east orientation horizontal velocity V constantly x(t 1) and north orientation horizontal velocity V y(t 1), ask for t 2Carrier positions is constantly:
Figure BSA00000768634100089
3) calculate the carrier positions error:
Figure BSA000007686341000810
Wherein: λ 0Longitude and the latitude of representing initial time carrier present position respectively;
Figure BSA000007686341000812
Δ λ representes the latitude of carrier, the variable quantity of longitude respectively; R N, R MThe radius-of-curvature of representing earth meridian circle, prime vertical respectively; t 1, t 2Two the adjacent time points in the process that resolve for inertial navigation system.

Claims (5)

1. fiber strapdown inertial navigation system error inhibition method based on three rotations is characterized in that may further comprise the steps:
(1) confirms the initial position parameters of carrier through GPS, they are bound to navigational computer;
(2) SINS carries out preheating and prepares, and gathers data that fibre optic gyroscope and quartz accelerometer export and data are handled;
(3) IMU adopts 12 commentaries on classics to stop the transposition scheme that order is a swing circle (like accompanying drawing 2);
Order 1, IMU rotates counterclockwise 180 ° of in-position B, stand-by time T from the A point sOrder 2, IMU rotates counterclockwise 180 ° of in-position C, stand-by time T from the B point sOrder 3, IMU rotates counterclockwise 180 ° of in-position A, stand-by time T from the C point sOrder 4, IMU rotates counterclockwise 180 ° of in-position C, stand-by time T from the A point sOrder 5, IMU rotates counterclockwise 180 ° of in-position B, stand-by time T from the C point sOrder 6, IMU rotates counterclockwise 180 ° of in-position A, stand-by time T from the B point sOrder 7, IMU clockwise rotates 180 ° of in-position B, stand-by time T from the A point sOrder 8, IMU clockwise rotates 180 ° of in-position C, stand-by time T from the B point sOrder 9, IMU clockwise rotates 180 ° of in-position A, stand-by time T from the C point sOrder 10, IMU clockwise rotates 180 ° of in-position C, stand-by time T from the A point sOrder 11, IMU clockwise rotates 180 ° of in-position B, stand-by time T from the C point sOrder 12, IMU clockwise rotates 180 ° of in-position A, stand-by time T from the B point sIMU rotates sequential loop according to this to carry out.
(4) data-switching that Inertial Measurement Unit rotation back gyroscope is generated obtains the modulation format that inertia device often is worth deviation under carrier coordinate system;
Suppose that the gyroscope constant value drift on the IMU horizontal direction is respectively ε xAnd ε yUnder the carrier quiescent conditions; Because three positions of A, B, C that IMU pauses are with respect to navigation coordinate system symmetry; Therefore on three fixed positions in one three transposition cycles, three gyroscope constant value drifts are fastened the attitude error that projection causes at navigation coordinate and must be satisfied:
3 ( ∫ 0 T s ϵ x n dt ) A + 3 ( ∫ 0 T s ϵ x n dt ) B + 3 ( ∫ 0 T s ϵ x n dt ) C = 0
3 ( ∫ 0 T s ϵ y n dt ) A + 3 ( ∫ 0 T s ϵ y n dt ) B + 3 ( ∫ 0 T s ϵ y n dt ) C = 0
3 ( ∫ 0 T s ϵ z n dt ) A + 3 ( ∫ 0 T s ϵ z n dt ) B + 3 ( ∫ 0 T s ϵ z n dt ) C = 0
Exist the symmetry problem of rotation according to the rotation in three scheme of rotation of IMU, ignore the influence of carrier movement and with local geographic coordinate system as a reference, 12 order transposition schemes can be expressed as:
Process 1: order 1,6,7,12, in the rotation period of formation, the gyroscopic drift of x, y axle is ox at navigation coordinate ny nDemonstrate the Changing Pattern in positive and negative each week in the plane, the normal value deviation that therefore in the integral process of complete cycle, produces is zero, that is:
( ∫ 0 T z ϵ x n dt ) A → + B + ( ∫ 0 T z ϵ x n dt ) B → + A + ( ∫ 0 T z ϵ x n dt ) A → - B + ( ∫ 0 T z ϵ x n dt ) B → - A = 0
( ∫ 0 T z ϵ y n dt ) A → + B + ( ∫ 0 T z ϵ y n dt ) B → + A + ( ∫ 0 T z ϵ y n dt ) A → - B + ( ∫ 0 T z ϵ y n dt ) B → - A = 0
Wherein, the time of each rotation process is counted T z, around the responsive coordinate axis of Inertial Measurement Unit rotate counterclockwise into+, clockwise rotate into-.
Process 2: order 2,5,8,11, in the rotation period of formation, the gyroscopic drift of y, z axle is oy at navigation coordinate nz nDemonstrate the Changing Pattern in positive and negative each week in the plane, the normal value deviation that therefore in the integral process of complete cycle, produces is zero, that is:
( ∫ 0 T z ϵ y n dt ) B → + C + ( ∫ 0 T z ϵ y n dt ) C → + B + ( ∫ 0 T z ϵ y n dt ) B → - C + ( ∫ 0 T z ϵ y n dt ) C → - B = 0
( ∫ 0 T z ϵ z n dt ) B → + C + ( ∫ 0 T z ϵ z n dt ) C → + B + ( ∫ 0 T z ϵ z n dt ) B → - C + ( ∫ 0 T z ϵ z n dt ) C → - B = 0
Process 3: order 3,4,9,10, in the rotation period of formation, the gyroscopic drift of x, z axle is ox at navigation coordinate nz nDemonstrate the Changing Pattern in positive and negative each week in the plane, the normal value deviation that therefore in the integral process of complete cycle, produces is zero, that is:
( ∫ 0 T z ϵ x n dt ) C → + A + ( ∫ 0 T z ϵ x n dt ) A → + C + ( ∫ 0 T z ϵ x n dt ) C → - A + ( ∫ 0 T z ϵ x n dt ) A → - C = 0
( ∫ 0 T z ϵ z n dt ) C → + A + ( ∫ 0 T z ϵ z n dt ) A → + C + ( ∫ 0 T z ϵ z n dt ) C → - A + ( ∫ 0 T z ϵ z n dt ) A → - C = 0
It is exactly the value that periodically changes the strapdown matrix that 12 order change the process of stopping, make three gyrostatic sensitive axes in a rotation period along center of rotation be symmetrically distributed (like accompanying drawing 3).Having proved intuitively that one 12 order changes stops in the process, and gyroscope often is worth the relative navigation coordinate of deviation system by modulation fully, and the navigation accuracy of system is not exerted an influence.In like manner stop in the cycle a complete commentaries on classics because the symmetrical distribution of three fixed positions and rotation process, can obtain that Inertial Measurement Unit stops and the transposition process in the accelerometer zero drift in the similar effect effect of navigation coordinate system.
(5) output valve
Figure FSA00000768634000027
of gyroscope under the IMU coordinate system brought in the strapdown inertial navigation system, adopts the equivalent rotating vector method that strapdown matrix
Figure FSA00000768634000028
is upgraded:
ω ns s = ω is s - ( C s n ) T ( ω ie n + ω en n )
Wherein:
Figure FSA000007686340000210
is the component of rotational-angular velocity of the earth under navigation system;
Figure FSA000007686340000211
is the component of motion angular velocity under navigation system of spherical coordinate system relatively for navigation coordinate; is the component of motion angular velocity on the IMU coordinate system of the relative navigation coordinate of IMU system.
If the equivalent rotating vector differential equation of the relative navigation coordinate of IMU coordinate system system is:
Φ · = ω ns s + 1 2 Φ × ω ns s + 1 12 Φ × ( Φ × ω ns s )
Solve the rotating vector of equivalence and replace hypercomplex number to separate according to angular velocity
Figure FSA00000768634000031
q = cos Φ 2 + Φ | Φ | sin Φ 2
Because q=q 0+ q 1I+q 2J+q 3K, i, j, k are direction vector.Therefore the renewal process of attitude matrix
Figure FSA00000768634000033
is:
C s n = q 0 2 + q 1 2 - q 2 2 - q 3 2 2 ( q 1 q 2 - q 0 q 3 ) 2 ( q 1 q 3 + q 0 q 2 ) 2 ( q 1 q 2 + q 0 q 3 ) q 0 2 - q 1 2 + q 2 2 - q 3 2 2 ( q 2 q 3 - q 0 q 1 ) 2 ( q 1 q 3 - q 0 q 2 ) 2 ( q 2 q 3 + q 0 q 1 ) q 0 2 - q 1 2 - q 2 2 + q 3 2
(6) the use of quartz accelerometer output value
Figure FSA00000768634000035
and steps (5) calculation of attitude matrix
Figure FSA00000768634000036
calculated after IMU rotation modulated carrier position.
1) calculates the acceleration f down of navigation system n:
f n = C s n f is s
2) position of calculating carrier:
According to t 1Carrier east orientation horizontal velocity V constantly x(t 1) and north orientation horizontal velocity V y(t 1), ask for t 2Carrier positions is constantly:
Figure FSA00000768634000038
3) calculate the carrier positions error:
Figure FSA00000768634000039
Wherein:
Figure FSA000007686340000310
λ 0Longitude and the latitude of representing initial time carrier present position respectively;
Figure FSA000007686340000311
Δ λ representes the latitude of carrier, the variable quantity of longitude respectively; R N, R MThe radius-of-curvature of representing earth meridian circle, prime vertical respectively; t 1, t 2Two the adjacent time points in the process that resolve for inertial navigation system.
2. the fiber strapdown inertial navigation system error inhibition method based on three rotations according to claim 1 is characterized in that adopting 12 commentaries on classics to stop the transposition scheme that order is a swing circle IMU, specifically comprises the steps:
Order 1, IMU rotates counterclockwise 180 ° of in-position B, stand-by time T from the A point sOrder 2, IMU rotates counterclockwise 180 ° of in-position C, stand-by time T from the B point sOrder 3, IMU rotates counterclockwise 180 ° of in-position A, stand-by time T from the C point sOrder 4, IMU rotates counterclockwise 180 ° of in-position C, stand-by time T from the A point sOrder 5, IMU rotates counterclockwise 180 ° of in-position B, stand-by time T from the C point sOrder 6, IMU rotates counterclockwise 180 ° of in-position A, stand-by time T from the B point sOrder 7, IMU clockwise rotates 180 ° of in-position B, stand-by time T from the A point sOrder 8, IMU clockwise rotates 180 ° of in-position C, stand-by time T from the B point sOrder 9, IMU clockwise rotates 180 ° of in-position A, stand-by time T from the C point sOrder 10, IMU clockwise rotates 180 ° of in-position C, stand-by time T from the A point sOrder 11, IMU clockwise rotates 180 ° of in-position B, stand-by time T from the C point sOrder 12, IMU clockwise rotates 180 ° of in-position A, stand-by time T from the B point sIMU rotates sequential loop according to this to carry out.
3. the fiber strapdown inertial navigation system error inhibition method based on three rotations according to claim 1; The data-switching that it is characterized in that Inertial Measurement Unit rotation back gyroscope is generated is under carrier coordinate system; Obtain the modulation format that inertia device often is worth deviation, specifically comprise the steps:
Suppose that the gyroscope constant value drift on the IMU horizontal direction is respectively ε xAnd ε yUnder the carrier quiescent conditions; Because three positions of A, B, C that IMU pauses are with respect to navigation coordinate system symmetry; Therefore on three fixed positions in one three transposition cycles, three gyroscope constant value drifts are fastened the attitude error that projection causes at navigation coordinate and must be satisfied:
3 ( ∫ 0 T s ϵ x n dt ) A + 3 ( ∫ 0 T s ϵ x n dt ) B + 3 ( ∫ 0 T s ϵ x n dt ) C = 0
3 ( ∫ 0 T s ϵ y n dt ) A + 3 ( ∫ 0 T s ϵ y n dt ) B + 3 ( ∫ 0 T s ϵ y n dt ) C = 0
3 ( ∫ 0 T s ϵ z n dt ) A + 3 ( ∫ 0 T s ϵ z n dt ) B + 3 ( ∫ 0 T s ϵ z n dt ) C = 0
Exist the symmetry problem of rotation according to the rotation in three scheme of rotation of IMU, ignore the influence of carrier movement and with local geographic coordinate system as a reference, 12 order transposition schemes can be expressed as:
Process 1: order 1,6,7,12, in the rotation period of formation, the gyroscopic drift of x, y axle is ox at navigation coordinate ny nDemonstrate the Changing Pattern in positive and negative each week in the plane, the normal value deviation that therefore in the integral process of complete cycle, produces is zero, that is:
( ∫ 0 T z ϵ x n dt ) A → + B + ( ∫ 0 T z ϵ x n dt ) B → + A + ( ∫ 0 T z ϵ x n dt ) A → - B + ( ∫ 0 T z ϵ x n dt ) B → - A = 0
( ∫ 0 T z ϵ y n dt ) A → + B + ( ∫ 0 T z ϵ y n dt ) B → + A + ( ∫ 0 T z ϵ y n dt ) A → - B + ( ∫ 0 T z ϵ y n dt ) B → - A = 0
Wherein, the time of each rotation process is counted T z, around the responsive coordinate axis of Inertial Measurement Unit rotate counterclockwise into+, clockwise rotate into-.
Process 2: order 2,5,8,11, in the rotation period of formation, the gyroscopic drift of y, z axle is oy at navigation coordinate nz nDemonstrate the Changing Pattern in positive and negative each week in the plane, the normal value deviation that therefore in the integral process of complete cycle, produces is zero, that is:
( ∫ 0 T z ϵ y n dt ) B → + C + ( ∫ 0 T z ϵ y n dt ) C → + B + ( ∫ 0 T z ϵ y n dt ) B → - C + ( ∫ 0 T z ϵ y n dt ) C → - B = 0
( ∫ 0 T z ϵ z n dt ) B → + C + ( ∫ 0 T z ϵ z n dt ) C → + B + ( ∫ 0 T z ϵ z n dt ) B → - C + ( ∫ 0 T z ϵ z n dt ) C → - B = 0
Process 3: order 3,4,9,10, in the rotation period of formation, the gyroscopic drift of x, z axle is ox at navigation coordinate nz nDemonstrate the Changing Pattern in positive and negative each week in the plane, the normal value deviation that therefore in the integral process of complete cycle, produces is zero, that is:
( ∫ 0 T z ϵ x n dt ) C → + A + ( ∫ 0 T z ϵ x n dt ) A → + C + ( ∫ 0 T z ϵ x n dt ) C → - A + ( ∫ 0 T z ϵ x n dt ) A → - C = 0
( ∫ 0 T z ϵ z n dt ) C → + A + ( ∫ 0 T z ϵ z n dt ) A → + C + ( ∫ 0 T z ϵ z n dt ) C → - A + ( ∫ 0 T z ϵ z n dt ) A → - C = 0
It is exactly the value that periodically changes the strapdown matrix that 12 order change the process of stopping, make three gyrostatic sensitive axes in a rotation period along center of rotation be symmetrically distributed (like accompanying drawing 3).Having proved intuitively that one 12 order changes stops in the process, and gyroscope often is worth the relative navigation coordinate of deviation system by modulation fully, and the navigation accuracy of system is not exerted an influence.In like manner stop in the cycle a complete commentaries on classics because the symmetrical distribution of three fixed positions and rotation process, can obtain that Inertial Measurement Unit stops and the transposition process in the accelerometer zero drift in the similar effect effect of navigation coordinate system.
4. the fiber strapdown inertial navigation system error inhibition method based on three rotations according to claim 1; It is characterized in that the output valve
Figure FSA00000768634000055
of gyroscope under the IMU coordinate system brought in the strapdown inertial navigation system, adopt the equivalent rotating vector method that strapdown matrix
Figure FSA00000768634000056
is upgraded:
ω ns s = ω is s - ( C s n ) T ( ω ie n + ω en n )
Wherein:
Figure FSA00000768634000058
is the component of rotational-angular velocity of the earth under navigation system;
Figure FSA00000768634000059
is the component of motion angular velocity under navigation system of spherical coordinate system relatively for navigation coordinate;
Figure FSA000007686340000510
is the component of motion angular velocity on the IMU coordinate system of the relative navigation coordinate of IMU system.
If the equivalent rotating vector differential equation of the relative navigation coordinate of IMU coordinate system system is:
Φ · = ω ns s + 1 2 Φ × ω ns s + 1 12 Φ × ( Φ × ω ns s )
Solve the rotating vector of equivalence and replace hypercomplex number to separate according to angular velocity
Figure FSA000007686340000512
q = cos Φ 2 + Φ | Φ | sin Φ 2
Because q=q 0+ q 1I+q 2J+q 3K, i, j, k are direction vector.Therefore the renewal process of attitude matrix
Figure FSA000007686340000514
is:
C s n = q 0 2 + q 1 2 - q 2 2 - q 3 2 2 ( q 1 q 2 - q 0 q 3 ) 2 ( q 1 q 3 + q 0 q 2 ) 2 ( q 1 q 2 + q 0 q 3 ) q 0 2 - q 1 2 + q 2 2 - q 3 2 2 ( q 2 q 3 - q 0 q 1 ) 2 ( q 1 q 3 - q 0 q 2 ) 2 ( q 2 q 3 + q 0 q 1 ) q 0 2 - q 1 2 - q 2 2 + q 3 2
5. the fiber strapdown inertial navigation system error inhibition method based on three rotations according to claim 1 is characterized in that utilizing the output valve
Figure FSA00000768634000061
of quartz accelerometer and the attitude matrix
Figure FSA00000768634000062
of step (5) calculating to calculate the position through IMU rotation modulation back carrier.
1) calculates the acceleration f down of navigation system n:
f n = C s n f is s
2) position of calculating carrier:
According to t 1Carrier east orientation horizontal velocity V constantly x(t 1) and north orientation horizontal velocity V y(t 1), ask for t 2Carrier positions is constantly:
3) calculate the carrier positions error:
Figure FSA00000768634000065
Wherein:
Figure FSA00000768634000066
λ 0Longitude and the latitude of representing initial time carrier present position respectively;
Figure FSA00000768634000067
Δ λ representes the latitude of carrier, the variable quantity of longitude respectively; R N, R MThe radius-of-curvature of representing earth meridian circle, prime vertical respectively; t 1, t 2Two the adjacent time points in the process that resolve for inertial navigation system.
CN201210305216.XA 2012-08-16 2012-08-16 Error suppressing method of fiber strap-down inertial navigation system based on three-axis rotation Expired - Fee Related CN102788598B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210305216.XA CN102788598B (en) 2012-08-16 2012-08-16 Error suppressing method of fiber strap-down inertial navigation system based on three-axis rotation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210305216.XA CN102788598B (en) 2012-08-16 2012-08-16 Error suppressing method of fiber strap-down inertial navigation system based on three-axis rotation

Publications (2)

Publication Number Publication Date
CN102788598A true CN102788598A (en) 2012-11-21
CN102788598B CN102788598B (en) 2014-12-03

Family

ID=47154074

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210305216.XA Expired - Fee Related CN102788598B (en) 2012-08-16 2012-08-16 Error suppressing method of fiber strap-down inertial navigation system based on three-axis rotation

Country Status (1)

Country Link
CN (1) CN102788598B (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102997919A (en) * 2012-11-22 2013-03-27 北京理工大学 Method for improving error inhibition effect of rotary type strapdown inertial navigation by insulation of carrier movement
CN103090865A (en) * 2013-01-06 2013-05-08 哈尔滨工程大学 Method for restraining attitude errors of modulation type strapdown inertial navigation system
CN103792843A (en) * 2014-01-24 2014-05-14 北京航天控制仪器研究所 Inertial platform fast rotating control method
CN106843293A (en) * 2017-01-06 2017-06-13 西北工业大学 A kind of gyro-stabilized platform drift rejection method for high speed rolling carrier
CN107631737A (en) * 2016-11-03 2018-01-26 重庆市皓颖科技发展有限公司 Indexable level meter and calibration method
CN109882157A (en) * 2019-04-09 2019-06-14 中油奥博(成都)科技有限公司 The optical fiber inertial navigation system and its data processing method of underground multi-component measurements instrument
CN113551688A (en) * 2021-05-27 2021-10-26 北京航天发射技术研究所 Quick non-support disassembly-free calibration method and device for vehicle-mounted positioning and directional navigation equipment
CN116858280A (en) * 2023-05-28 2023-10-10 中国人民解放军国防科技大学 Full-parameter error comprehensive modulation method for laser gyro biaxial rotation inertial navigation system

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0783673A (en) * 1993-09-10 1995-03-28 Nippon Sharyo Seizo Kaisha Ltd Azimuth estimation device of traveling body
CN101514899A (en) * 2009-04-08 2009-08-26 哈尔滨工程大学 Optical fibre gyro strapdown inertial navigation system error inhibiting method based on single-shaft rotation
CN101629830A (en) * 2009-08-20 2010-01-20 北京航空航天大学 Calibration method and device of three-axis integrative high precision fiber optic gyro
CN101718560A (en) * 2009-11-20 2010-06-02 哈尔滨工程大学 Strapdown system error inhibition method based on uniaxial four-position rotation and stop scheme

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0783673A (en) * 1993-09-10 1995-03-28 Nippon Sharyo Seizo Kaisha Ltd Azimuth estimation device of traveling body
CN101514899A (en) * 2009-04-08 2009-08-26 哈尔滨工程大学 Optical fibre gyro strapdown inertial navigation system error inhibiting method based on single-shaft rotation
CN101629830A (en) * 2009-08-20 2010-01-20 北京航空航天大学 Calibration method and device of three-axis integrative high precision fiber optic gyro
CN101718560A (en) * 2009-11-20 2010-06-02 哈尔滨工程大学 Strapdown system error inhibition method based on uniaxial four-position rotation and stop scheme

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
孙枫,孙伟: ""基于单轴旋转的光纤捷联惯导系统误差特性与实验分析"", 《宇航学报》 *
孙枫,孙伟: ""旋转自动补偿捷联惯导系统技术研究"", 《系统工程与电子技术》 *
车莉娜,孙伟,徐爱功: ""IMU转动角速度对旋转SINS的精度影响分析"", 《仪器仪表学报》 *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102997919A (en) * 2012-11-22 2013-03-27 北京理工大学 Method for improving error inhibition effect of rotary type strapdown inertial navigation by insulation of carrier movement
CN102997919B (en) * 2012-11-22 2015-07-15 北京理工大学 Method for improving error inhibition effect of rotary type strapdown inertial navigation by insulation of carrier movement
CN103090865A (en) * 2013-01-06 2013-05-08 哈尔滨工程大学 Method for restraining attitude errors of modulation type strapdown inertial navigation system
CN103090865B (en) * 2013-01-06 2015-08-12 哈尔滨工程大学 A kind of modulation type strapdown inertial navigation system attitude error suppressing method
CN103792843A (en) * 2014-01-24 2014-05-14 北京航天控制仪器研究所 Inertial platform fast rotating control method
CN107631737A (en) * 2016-11-03 2018-01-26 重庆市皓颖科技发展有限公司 Indexable level meter and calibration method
CN106843293A (en) * 2017-01-06 2017-06-13 西北工业大学 A kind of gyro-stabilized platform drift rejection method for high speed rolling carrier
CN109882157A (en) * 2019-04-09 2019-06-14 中油奥博(成都)科技有限公司 The optical fiber inertial navigation system and its data processing method of underground multi-component measurements instrument
CN113551688A (en) * 2021-05-27 2021-10-26 北京航天发射技术研究所 Quick non-support disassembly-free calibration method and device for vehicle-mounted positioning and directional navigation equipment
CN116858280A (en) * 2023-05-28 2023-10-10 中国人民解放军国防科技大学 Full-parameter error comprehensive modulation method for laser gyro biaxial rotation inertial navigation system
CN116858280B (en) * 2023-05-28 2024-03-01 中国人民解放军国防科技大学 Full-parameter error comprehensive modulation method for laser gyro biaxial rotation inertial navigation system

Also Published As

Publication number Publication date
CN102788598B (en) 2014-12-03

Similar Documents

Publication Publication Date Title
CN102788598B (en) Error suppressing method of fiber strap-down inertial navigation system based on three-axis rotation
CN101514900B (en) Method for initial alignment of a single-axis rotation strap-down inertial navigation system (SINS)
CN101718560B (en) Strapdown system error inhibition method based on uniaxial four-position rotation and stop scheme
CN101514899B (en) Optical fibre gyro strapdown inertial navigation system error inhibiting method based on single-shaft rotation
Sun et al. MEMS-based rotary strapdown inertial navigation system
CN103090867A (en) Error restraining method for fiber-optic gyroscope strapdown inertial navigation system rotating relative to geocentric inertial system
CN101629826A (en) Coarse alignment method for fiber optic gyro strapdown inertial navigation system based on single axis rotation
CN101701825A (en) High-precision laser gyroscope single-shaft rotating inertial navigation system
CN101493008A (en) Strapping inertial navigation gyroscope clinometer based on MEMS device
CN101713666B (en) Single-shaft rotation-stop scheme-based mooring and drift estimating method
CN102829781B (en) Implementation method of rotation type strapdown optical-fiber compass
CN101419080B (en) Mini quick-connecting inertia measurement system zero speed correcting method
CN101963512A (en) Initial alignment method for marine rotary fiber-optic gyroscope strapdown inertial navigation system
CN103344226B (en) A kind of north-seeking system based on MEMS rotation technique and method
CN102798399A (en) SINS error inhibiting method based on biaxial rotation scheme
CN105973271A (en) Self-calibration method of hybrid type inertial navigation system
CN101162147A (en) Marine fiber optic gyroscope attitude heading reference system mooring extractive alignment method under the large heading errors
CN103148854A (en) Attitude measurement method of micro-electro mechanical system (MEMS) inertial navigation system based on single-shaft forward revolution and reverse revolution
CN103256943A (en) Compensation method for scale factor error in single-axial rotating strapdown inertial navigation system
CN102207386A (en) North-finding method based on orientation effect error compensation
CN103743413A (en) Installation error online estimation and north-seeking error compensation method for modulating north seeker under inclined state
CN102788597B (en) Error suppressing method of rotary strap-down inertial navigation system based on space stabilization
CN109752000A (en) A kind of MEMS dual-axis rotation modulation type strapdown compass Initial Alignment Method
CN103090866A (en) Method for restraining speed errors of single-shaft rotation optical fiber gyro strapdown inertial navigation system
CN103925930B (en) A kind of compensation method of gravimeter biax gyrostabilized platform course error effect

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20141203

Termination date: 20150816

EXPY Termination of patent right or utility model