CN100588906C - Carrier posture measuring method suitable for optical fiber gyroscope - Google Patents
Carrier posture measuring method suitable for optical fiber gyroscope Download PDFInfo
- Publication number
- CN100588906C CN100588906C CN200710144846A CN200710144846A CN100588906C CN 100588906 C CN100588906 C CN 100588906C CN 200710144846 A CN200710144846 A CN 200710144846A CN 200710144846 A CN200710144846 A CN 200710144846A CN 100588906 C CN100588906 C CN 100588906C
- Authority
- CN
- China
- Prior art keywords
- omega
- carrier
- optical fibre
- angle
- posture renewal
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Gyroscopes (AREA)
Abstract
The invention discloses a measuring method for a carrier gesture of an optical fiber gyro, which comprises determining an initial position parameter and an initial speed value of a carrier through anexternal device, initially aligning strapdown inertial navigation system of the optical fiber gyro, determining an initial gesture of the carrier which is relative to a navigation coordinate system, obtaining an initial value of a gesture hypercomplex number, determining a gesture renewal cycle H = tm-tm-1, calculating an increment delta phi of rotating vector through collecting a carrier which isoutput by the optical fiber gyro which is corresponded to an angle speed of the inertial coordinate system, obtaining a gesture renewal hypercomplex number q(H) of the gesture renewal cycle H througha relation of the rotating vector and the hypercomplex number, renewing the gesture hypercomplex number through a renewal equation of the gesture hypercomplex number, calculating a strapdown matrix Tthat a carrier coordinate system b is relative to a navigation coordinate system n, and obtaining a gesture angle of the carrier which is relative to the navigation coordinate system and the like. The invention solves the problem that a cone effect influences a measuring accuracy of the carrier gesture under a high dynamic environment of a carrier or a high-frequency vibration environment.
Description
(1) technical field
The present invention is to provide a kind of measuring method, specifically a kind of to use optical fibre gyro to measure carrier be the method for attitude parameter with respect to navigation coordinate.
(2) background technology
In strapdown inertial navigation system, inertial sensor (gyro and accelerometer) directly connects firmly on carrier.Their sensitive carriers change and velocity variations information with respect to the attitude of inertial coordinates system, the processing output navigational parameter of these information via navigational computers, thus finish the navigational guidance task.Because the mode that adopts inertial sensor to connect firmly, the vibration of the direct sensitive carrier of inertial sensor meeting and the interference of environment bring certain influence to the measuring accuracy of strapdown inertial navigation system.In strapdown inertial navigation system, employing be mathematical platform, mainly rely on the output information of gyro to make up mathematical platform, to describe the attitude of the relative navigation coordinate of carrier system.And gyro is output as angle increment or angular velocity, can not be directly used in the relativeness of describing coordinate system.In the attitude measurement method of the strapdown inertial navigation system of routine, always think that the conversion between carrier coordinate system and navigation coordinate system realizes by a series of rotation.So,, need try to achieve hypercomplex number or the rotating vector of describing the coordinate system relation by corresponding angular velocity or angle increment in order to measure the attitude of carrier.But under high dynamic environment, be different in the position of each gyro sampling corresponding vectors coordinate system of posture renewal in the cycle.In addition, be that the order of the rotation of hypothesis between from the carrier coordinate system to the reference frame can not considered in the classic method, this is based on infinitely small rotation is that the principle of vector obtains.But in the practical project, particularly under high dynamic mobility situation of carrier or abominable vibration environment, the rotation of carrier usually is limited rotation, and limited rotation is not a vector, and it rotates order and can not exchange.Thereby traditional attitude measurement method of finding the solution hypercomplex number after the angle increment summation with all gyro outputs of posture renewal in the cycle again will bring bigger error, i.e. coning error.Coning error is similar to gyroscopic drift, brings negative effect for the attitude measurement of carrier.The complement taper error, the attitude of measuring carrier in high quality just must the high performance attitude measurement method of design.
In traditional attitude measurement method, utilize the angle increment of gyro output to go to estimate the circular cone compensation term, thereby estimate the increment of rotating vector.These basic ideas are based in the output of traditional gyro, and obtaining of angle increment is very easily.The increment that promptly will expect rotating vector must be known the angle increment of gyro output.And in present widely used fiber optic gyro strapdown inertial navigation system, the output signal of optical fibre gyro is an angular velocity information, and for this reason, common way is to adopt piece-wise linearization, thinks that the angular velocity of carrier is constant in each sampling period.Set up when this hypothesis and angular velocity varies carrier enough little in the sampling period is slower, considered but just must add in addition for the dynamic occasion of height.Available research achievements shows, under the situation of optical fibre gyro output angle speed, for traditional attitude measurement method based on rotating vector, performance can't be embodied in the superiority of reply cone effect not as the conventional attitude measurement method based on hypercomplex number.
(3) summary of the invention
The object of the present invention is to provide a kind of can the solution in carrier high dynamic environment or high-frequency vibration environment, the carrier posture measuring method that is suitable for optical fibre gyro of the problem that cone effect exerts an influence for the attitude of carrier measuring accuracy.
The object of the present invention is achieved like this:
Step 1, determine the initial position parameters and the initial velocity value of carrier by external unit;
The carrier of step 4, the output of collection optical fibre gyro is with respect to the increment Delta φ of the angular speed calculation rotating vector of inertial coordinates system;
Step 5, by the relation of rotating vector and hypercomplex number, obtain posture renewal hypercomplex number q (H) in the posture renewal cycle H,
Wherein posture renewal hypercomplex number q (H) is obtained by step 5, Q (t
m), Q (t
M-1) represent that respectively carrier is at t
m, t
M-1Attitude quaternion constantly;
Step 7, the t that utilizes step 6 to obtain
mMoment attitude quaternion Q (t
m)=[q
0q
1q
2q
3]
TCalculating carrier coordinate system b system is the strapdown matrix T of n system with respect to navigation coordinate,
θ
0=arcsinT
32
Wherein the field of definition of grid heading angle ψ is (0 a °, 360 °); The field of definition of pitch angle θ is (90 °, 90 °); The field of definition of roll angle γ is (90 °, 90 °),
Obtain their true value by the field of definition of attitude angle,
The true value at grid heading angle is
The true value of pitch angle is
θ=θ
0
The true value of roll angle is
The present invention can also comprise:
1, the initial position at the carrier described in the step 1 is provided by GPS device or outside high-precision integrated navigation equipment, and the initial velocity value of described carrier is provided by DVL Doppler log or outside high-precision integrated navigation equipment.
2, determine to be divided into two kinds of situations at the initial attitude of the relative navigation coordinate of the carrier described in the step 2 system: when carrier is in quiet pedestal, adopt Alignment Method based on classic control theory; When carrier is in moving pedestal, adopt combination alignment methods based on kalman filtering theory.
3, the calculating concrete steps at the rotating vector increment Delta φ described in the step 4 are as follows:
Step 41: can survey item by the inertia that the output of optical fibre gyro output is calculated among the Δ φ, i.e. angle increment α; The H=Nh that concerns by posture renewal cycle and optical fibre gyro sampling period; Known N+1 optical fibre gyro sampling angular velocity omega in each posture renewal cycle H
0, ω
1, L ω
NThen the angle increment α in the posture renewal cycle H is
Wherein
Step 42: can survey item by the non-inertia that the output of optical fibre gyro output is calculated among the Δ φ, i.e. circular cone compensation term β; The H=Nh that concerns by posture renewal cycle and optical fibre gyro sampling period; Known N+1 optical fibre gyro sampling angular velocity omega in each posture renewal cycle H
0, ω
1, L ω
N
In the posture renewal cycle H, circular cone compensation term β is calculated by the linear combination of optical fibre gyro sampling angular velocity multiplication cross item
K wherein
N-iFor optimizing coefficient, in the typical cone environment, determine;
Step 43:,, obtain rotating vector increment Delta φ in the posture renewal cycle H with the posture renewal cycle H female cone compensation term β addition that obtains in the step 41 with angle increment α in the posture renewal cycle H that obtains in the step 41;
Δφ=α+β
N=3 is set, and gyro is at the t angular velocity omega of sampling constantly
0, and each posture renewal cycle H is at t+H/3, t+2H/3, the t+H ω that samples constantly
1, ω
2, ω
3Obtaining the interior rotating vector increment Delta φ of posture renewal cycle H is
Wherein
N=2 is set, and gyro is at the t angular velocity omega of sampling constantly
0, and in each computation period H, at t+H/2, the t+H ω that samples constantly twice
1, ω
2Obtaining the interior rotating vector increment Delta φ of posture renewal cycle H is
Wherein
The present invention compare with classic method advantage be mainly reflected in:
(1) be output as the situation of angular velocity at optical fibre gyro, be different from traditional measurement method with angle increment as input, inventive method is directly with the angular velocity input, fits the circular cone compensation term with the multiplication cross item of angular velocity.When having avoided the use optical fibre gyro, classic method is for the measured deviation of sculling compensation item.
(2) with respect to the another kind of approach that improves strapdown inertial navigation system Attitude Calculation precision: select high performance optical fibre gyro for use.The measuring method of invention only need be done the navigation software design, does not need to increase the manufacturing cost of strapdown inertial navigation system.
(4) description of drawings
Fig. 1 is the attitude measurement method process flow diagram of suitable fiber optic gyro strapdown inertial navigation system of the present invention.
Fig. 2 is in the typical cone environment, the attitude measurement method of invention; The empirical curve of traditional attitude measurement method.The typical cone environment is defined as: have same frequency on two orthogonal axes of carrier, the angular oscillation of phase differential 90 degree.The sample frequency of optical fibre gyro is 100Hz.
(5) embodiment
For example the present invention is done description in more detail below in conjunction with accompanying drawing:
In conjunction with Fig. 1, the attitude measurement method that is suitable for fiber optic gyro strapdown inertial navigation system of the present invention comprises the steps:
Step 1, determine the initial position parameters and the initial velocity value of carrier by external unit.
The carrier of step 4, the output of collection optical fibre gyro is with respect to the increment Delta φ of the angular speed calculation rotating vector of inertial coordinates system.
Step 5, by the relation of rotating vector and hypercomplex number, obtain posture renewal hypercomplex number q (H) in the posture renewal cycle H.
Wherein posture renewal hypercomplex number q (H) is tried to achieve by step 5.Q (t
m), Q (t
M-1) represent that respectively carrier is at t
m, t
M-1Attitude quaternion constantly.
Step 7, the t that utilizes step 6 to obtain
mMoment attitude quaternion Q (t
m)=[q
0q
1q
2q
3]
TCalculating carrier coordinate system b system is the strapdown matrix T of n system with respect to navigation coordinate.
θ
0=arcsinT
32(4)
Wherein the field of definition of grid heading angle ψ is (0 a °, 360 °); The field of definition of pitch angle θ is (90 °, 90 °); The field of definition of roll angle γ is (90 °, 90 °).
Can determine their true value by the field of definition of attitude angle.The true value at grid heading angle is
The true value of pitch angle is
θ=θ
0
The true value of roll angle is
Initial position at the carrier described in the step 1 is provided by GPS device or outside high-precision integrated navigation equipment, and the initial velocity value of described carrier is provided by DVL Doppler log or outside high-precision integrated navigation equipment.
Initial attitude in the relative navigation coordinate of the carrier described in the step 2 system determines to be divided into two kinds of situations: when carrier is in quiet pedestal, can adopt the Alignment Method based on classic control theory; When carrier is in moving pedestal, can adopt combination alignment methods based on kalman filtering theory.
Calculating concrete steps at the rotating vector increment Delta φ described in the step 4 are as follows:
Step 41: can survey item by the inertia that the output of optical fibre gyro output is calculated among the Δ φ, i.e. angle increment α.The H=Nh that concerns by posture renewal cycle and optical fibre gyro sampling period.Known N+1 optical fibre gyro sampling angular velocity omega in each posture renewal cycle H
0, ω
1, L ω
NThen the angle increment α in the posture renewal cycle H is
Wherein
Step 42: can survey item by the non-inertia that the output of optical fibre gyro output is calculated among the Δ φ, i.e. circular cone compensation term β.The H=Nh that concerns by posture renewal cycle and optical fibre gyro sampling period.Known N+1 optical fibre gyro sampling angular velocity omega in each posture renewal cycle H
0, ω
1, L ω
N
In the posture renewal cycle H, circular cone compensation term β can be calculated by the sample linear combination of angular velocity multiplication cross item of optical fibre gyro
K wherein
N-iFor optimizing coefficient, can in the typical cone environment, determine.
Step 43:,, obtain rotating vector increment Delta φ in the posture renewal cycle H with the posture renewal cycle H female cone compensation term β addition that obtains in the step 41 with angle increment α in the posture renewal cycle H that obtains in the step 41.
Δφ=α+β(8)
N=3 is set, and gyro is at the t angular velocity omega of sampling constantly
0, and each posture renewal cycle H is at t+H/3, t+2H/3, the t+H ω that samples constantly
1, ω
2, ω
3Obtaining the interior rotating vector increment Delta φ of posture renewal cycle H is
Wherein
N=2 is set, and gyro is at the t angular velocity omega of sampling constantly
0, and in each computation period H, at t+H/2, the t+H ω that samples constantly twice
1, ω
2Obtaining the interior rotating vector increment Delta φ of posture renewal cycle H is
Wherein
In the posture renewal cycle H, the employed optical fibre gyro hits of attitude measurement method is many more, and the effect of attitude measurement method reply circular cone compensation is just good more, and the attitude accuracy of fiber optic gyro strapdown inertial navigation system output is just high more.
Claims (3)
1, a kind of carrier posture measuring method that is suitable for optical fibre gyro is characterized in that:
Step 1, determine the initial position parameters and the initial velocity value of carrier by external unit;
Step 2, fiber optic gyro strapdown inertial navigation system carry out initial alignment, determine the initial attitude of the relative navigation coordinate of carrier system, obtain the initial value of attitude quaternion;
Step 3, determine posture renewal cycle H=t
m-t
M-1, described posture renewal cycle H equals N optical fibre gyro sampling period h=t doubly
l-t
L-1, described N is the integer greater than zero;
The carrier of step 4, the output of collection optical fibre gyro is with respect to the increment Delta φ of the angular speed calculation rotating vector of inertial coordinates system, and concrete steps are as follows:
Step 41: can survey item by the inertia that the output of optical fibre gyro output is calculated among the Δ φ, i.e. angle increment α; The H=Nh that concerns by posture renewal cycle and optical fibre gyro sampling period; Known N+1 optical fibre gyro sampling angular velocity omega in each posture renewal cycle H
0, ω
1... ω
NThen the angle increment α in the posture renewal cycle H is
Wherein
Step 42: can survey item by the non-inertia that the output of optical fibre gyro output is calculated among the Δ φ, i.e. circular cone compensation term β; The H=Nh that concerns by posture renewal cycle and optical fibre gyro sampling period; Known N+1 optical fibre gyro sampling angular velocity omega in each posture renewal cycle H
0, ω
1... ω
N
In the posture renewal cycle H, circular cone compensation term β is calculated by the linear combination of optical fibre gyro sampling angular velocity multiplication cross item
K wherein
N-iFor optimizing coefficient, in the typical cone environment, determine;
Step 43:,, obtain rotating vector increment Delta φ in the posture renewal cycle H with the posture renewal cycle H female cone compensation term β addition that obtains in the step 41 with angle increment α in the posture renewal cycle H that obtains in the step 41;
Δφ=α+β
N=3 is set, and gyro is at the t angular velocity omega of sampling constantly
0, and each posture renewal cycle H is at t+H/3, t+2H/3, the t+H ω that samples constantly
1, ω
2, ω
3Obtaining the interior rotating vector increment Delta φ of posture renewal cycle H is
Wherein
N=2 is set, and gyro is at the t angular velocity omega of sampling constantly
0, and in each computation period H, at t+H/2, the t+H ω that samples constantly twice
1, ω
2Obtaining the interior rotating vector increment Delta φ of posture renewal cycle H is
Wherein
Step 5, by the relation of rotating vector and hypercomplex number, obtain posture renewal hypercomplex number q (H) in the posture renewal cycle H,
Step 6, upgrade attitude quaternion by the attitude quaternion renewal equation
Wherein posture renewal hypercomplex number q (H) is obtained by step 5, Q (t
m), Q (t
M-1) represent that respectively carrier is at t
m, t
M-1Attitude quaternion constantly;
Step 7, the t that utilizes step 6 to obtain
mMoment attitude quaternion Q (t
m)=[q
0q
1q
2q
3]
TCalculating carrier coordinate system b system is the strapdown matrix T of n system with respect to navigation coordinate,
Step 8, the strapdown matrix T of utilizing step 7 to obtain are tried to achieve the attitude angle of the relative navigation coordinate of carrier system, navigation coordinate system is chosen for the azimuthal coordinates system of moving about, then the grid heading angle ψ of carrier, pitch angle θ, roll angle γ can usually be represented by the unit in the strapdown matrix T
θ
0=arcsin?T
32
Wherein the field of definition of grid heading angle ψ is (0 a °, 360 °); The field of definition of pitch angle θ is (90 °, 90 °); The field of definition of roll angle γ is (90 °, 90 °),
Obtain their true value by the field of definition of attitude angle, the true value at grid heading angle is
The true value of pitch angle is
θ=θ
0
The true value of roll angle is
2, the carrier posture measuring method that is suitable for optical fibre gyro according to claim 1, it is characterized in that: the initial position at the carrier described in the step 1 is provided by GPS device or outside high-precision integrated navigation equipment, and the initial velocity value of described carrier is provided by DVL Doppler log or outside high-precision integrated navigation equipment.
3, the carrier posture measuring method that is suitable for optical fibre gyro according to claim 1 and 2, it is characterized in that: the initial attitude in the relative navigation coordinate of the carrier described in the step 2 system determines to be divided into two kinds of situations: when carrier is in quiet pedestal, adopt the Alignment Method based on classic control theory; When carrier is in moving pedestal, adopt combination alignment methods based on kalman filtering theory.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200710144846A CN100588906C (en) | 2007-12-18 | 2007-12-18 | Carrier posture measuring method suitable for optical fiber gyroscope |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200710144846A CN100588906C (en) | 2007-12-18 | 2007-12-18 | Carrier posture measuring method suitable for optical fiber gyroscope |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101187561A CN101187561A (en) | 2008-05-28 |
CN100588906C true CN100588906C (en) | 2010-02-10 |
Family
ID=39480052
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200710144846A Expired - Fee Related CN100588906C (en) | 2007-12-18 | 2007-12-18 | Carrier posture measuring method suitable for optical fiber gyroscope |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100588906C (en) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102997932B (en) * | 2011-09-15 | 2015-08-05 | 北京自动化控制设备研究所 | A kind of method eliminating high accuracy inertial navigation system demarcation intermediate station effect of jitter |
CN102768043B (en) * | 2012-06-14 | 2015-03-04 | 辽宁工程技术大学 | Integrated attitude determination method without external observed quantity for modulated strapdown system |
CN103148868B (en) * | 2013-02-22 | 2015-12-09 | 哈尔滨工程大学 | Based on the combination alignment methods that Doppler log Department of Geography range rate error is estimated under at the uniform velocity sailing through to |
JP2015094631A (en) * | 2013-11-11 | 2015-05-18 | セイコーエプソン株式会社 | Position calculation device, and position calculation method |
CN103792843B (en) * | 2014-01-24 | 2016-05-04 | 北京航天控制仪器研究所 | A kind of inertial platform quick rotation control method |
CN103900608B (en) * | 2014-03-21 | 2016-08-17 | 哈尔滨工程大学 | A kind of low precision inertial alignment method based on quaternary number CKF |
CN106997053B (en) * | 2017-02-24 | 2019-11-12 | 北京航天自动控制研究所 | A kind of navigation Consistency test method of three strapdown inertial measurement unit of stand base altogether |
CN109297487A (en) * | 2017-07-25 | 2019-02-01 | 北京信息科技大学 | A kind of attitude decoupling method under the conditions of turning rate input |
CN110345942B (en) * | 2019-07-17 | 2023-05-30 | 哈尔滨工程大学 | Interpolation three subsampled cone error compensation algorithm based on angular rate input |
CN111351483B (en) * | 2020-03-31 | 2021-12-07 | 北京控制工程研究所 | Recursive multi-subsample large dynamic inertial navigation method |
CN112857366B (en) * | 2021-01-14 | 2022-03-18 | 东南大学 | Optical fiber strapdown inertial navigation system attitude calculation method based on compression structure |
CN113720206B (en) * | 2021-09-02 | 2023-04-11 | 重庆零壹空间科技集团有限公司 | Rocket ground aiming method, system, computer equipment and storage medium |
-
2007
- 2007-12-18 CN CN200710144846A patent/CN100588906C/en not_active Expired - Fee Related
Non-Patent Citations (9)
Title |
---|
Navigation Computation in Terrestrial Strapdown InertialNavigation Systems. ITZHACK Y.BAR-ITZHACK.IEEE TRANSACTIONS ON AEROSPACE AND ELECTRONIC SYSTEMS,Vol.AES-13 No.6. 1977 |
Navigation Computation in Terrestrial Strapdown InertialNavigation Systems. ITZHACK Y.BAR-ITZHACK.IEEE TRANSACTIONS ON AEROSPACE AND ELECTRONIC SYSTEMS,Vol.AES-13 No.6. 1977 * |
光纤陀螺及其船用航姿***的姿态解算. 曹晓平等.光学与光电技术,第1卷第3期. 2003 |
光纤陀螺及其船用航姿***的姿态解算. 曹晓平等.光学与光电技术,第1卷第3期. 2003 * |
基于多普勒和光纤陀螺水下机器人导航***研究. 俞建成等.航海工程,第2期. 2006 |
基于多普勒和光纤陀螺水下机器人导航***研究. 俞建成等.航海工程,第2期. 2006 * |
惯性导航. 秦永元,237,科学出版社. 2007 惯性导航初始对准. 万德钧,48-82,东南大学出版社. 1998 |
捷联惯导算法中螺旋理论的应用. 孙枫等.***工程与电子技术,第29卷第9期. 2007 |
捷联惯导算法中螺旋理论的应用. 孙枫等.***工程与电子技术,第29卷第9期. 2007 * |
Also Published As
Publication number | Publication date |
---|---|
CN101187561A (en) | 2008-05-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100588906C (en) | Carrier posture measuring method suitable for optical fiber gyroscope | |
CN100547352C (en) | The ground speed testing methods that is suitable for fiber optic gyro strapdown inertial navigation system | |
CN101173858B (en) | Three-dimensional posture fixing and local locating method for lunar surface inspection prober | |
US10670424B2 (en) | Method for initial alignment of an inertial navigation apparatus | |
CN105300379B (en) | A kind of Kalman filtering Attitude estimation method and system based on acceleration | |
CN107560613B (en) | Robot indoor track tracking system and method based on nine-axis inertial sensor | |
CN101701825A (en) | High-precision laser gyroscope single-shaft rotating inertial navigation system | |
CN101629826A (en) | Coarse alignment method for fiber optic gyro strapdown inertial navigation system based on single axis rotation | |
CN107063254B (en) | Gesture resolving method for gyros and geomagnetic combination | |
CN103697878B (en) | A kind of single gyro list accelerometer rotation modulation north finding method | |
CN103900608A (en) | Low-precision inertial navigation initial alignment method based on quaternion CKF | |
CN103499348A (en) | High-precision attitude data calculation method for AHRS (Attitude and Heading Reference System) | |
CN109931952A (en) | The direct analytic expression coarse alignment method of inertial navigation under the conditions of unknown latitude | |
CN103900566A (en) | Method for eliminating influence on accuracy of rotating modulation strapdown inertial navigation system caused by earth rotation angular velocity | |
CN103955005B (en) | A kind of rocket sledge track gravity method for real-time measurement | |
Vaknin et al. | Coarse leveling of gyro-free INS | |
CN113551669B (en) | Combined navigation positioning method and device based on short base line | |
Wongwirat et al. | A position tracking experiment of mobile robot with inertial measurement unit (imu) | |
CN104297525A (en) | Accelerometer calibration method for inertia measurement system on basis of rocket sled test | |
CN104101345A (en) | Multisensor attitude fusion method based on complementary reconstruction technology | |
Zhe et al. | Adaptive complementary filtering algorithm for imu based on mems | |
CN111649738B (en) | Method for calculating initial attitude of accelerometer under microgravity field | |
TW498169B (en) | Interruption-free hand-held positioning method and system thereof | |
WO2018126911A1 (en) | Positioning method and apparatus, service processing system and computer readable storage medium | |
CN106705954B (en) | Floating ball platform carrier posture measuring method |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20100210 Termination date: 20171218 |