CN1740746A - Micro-dynamic carrier attitude measuring apparatus and measuring method thereof - Google Patents
Micro-dynamic carrier attitude measuring apparatus and measuring method thereof Download PDFInfo
- Publication number
- CN1740746A CN1740746A CN 200510011763 CN200510011763A CN1740746A CN 1740746 A CN1740746 A CN 1740746A CN 200510011763 CN200510011763 CN 200510011763 CN 200510011763 A CN200510011763 A CN 200510011763A CN 1740746 A CN1740746 A CN 1740746A
- Authority
- CN
- China
- Prior art keywords
- state variable
- micro
- carrier
- attitude
- analog
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims description 28
- 238000006243 chemical reaction Methods 0.000 claims abstract description 15
- 238000004891 communication Methods 0.000 claims abstract description 8
- 230000001133 acceleration Effects 0.000 claims description 26
- 238000005259 measurement Methods 0.000 claims description 18
- 239000011159 matrix material Substances 0.000 claims description 12
- 238000005070 sampling Methods 0.000 claims description 10
- 230000033001 locomotion Effects 0.000 claims description 9
- 230000003068 static effect Effects 0.000 claims description 4
- 238000012937 correction Methods 0.000 claims description 3
- 238000009434 installation Methods 0.000 claims 1
- 210000000056 organ Anatomy 0.000 claims 1
- 238000000691 measurement method Methods 0.000 abstract 1
- 230000005484 gravity Effects 0.000 description 4
- 238000001228 spectrum Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 2
- 239000000969 carrier Substances 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- -1 aircraft Substances 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000009191 jumping Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
Images
Landscapes
- Gyroscopes (AREA)
- Navigation (AREA)
Abstract
The measuring equipment includes three-axis rate gyroscope, three-axis magnetic field meter, uniaxial accelerometer, temperature sensor, A/D conversion circuit, microprocessor and memory, serial communication interface and its correspondent software. Said invention also provides the concrete steps of its measurement method.
Description
Technical field
The present invention relates to a kind of microminiature carrier attitude measuring apparatus, be specially adapted to dynamic carrier attitude and measure, belong to signal Processing, measurement, navigation field.
The invention still further relates to a kind of measuring method according to this device.
Background technology
The measurement of attitude information has crucial meaning for motion carriers such as aircraft, vehicle, boats and ships.In the classic method, people use the attitude angle of gyroscope survey motion carrier.Good gyroscope has quite high precision, but simultaneously it also exist cost an arm and a leg, complex structure, volume are big, quality heavily waits shortcoming, this makes its application be subjected to considerable restraint.The rate gyro price is lower, volume is little, light weight, but its precision is lower, and drift is serious, therefore can not directly be used for measuring attitude angle.
3 components that utilize acceleration transducer and geomagnetic sensor to measure acceleration of gravity and terrestrial magnetic field also can calculate the attitude angle (as United States Patent (USP) 20020188416) of carrier.The volume of this class sensor can be accomplished very little, but because therefore also responsive acceleration of motion of acceleration transducer while is subjected to acceleration of motion easily and disturbs, can not be used under the nonequilibrium condition.
Another kind method is on the basis of acceleration and geomagnetic sensor, adds rate gyro, and structure is estimated (as the United States Patent (USP) 6647352 of having authorized) based on the Kalman filter of hypercomplex number in real time to attitude angle.The method is utilized 3 rate gyro output update mode amounts, and the measured value that adopts 3 axle accelerations and 3 magnetometers obtains the optimal estimation of state variable as observed reading under the least square meaning, and then obtains the estimated value of attitude angle.The advantage of the method is the random disturbance of filtering sensor signal preferably, and the drift of compensation gyro has improved precision effectively.But it requires the measured value of acceleration of gravity in long-time is no inclined to one side, that is to say, when carrier is in when being in acceleration (perhaps slow down, turn) state for a long time, the precision of the method will seriously reduce, and will be even not available.In order to address this problem, people have proposed to adjust automatically according to different conditions the method for wave filter noise.This method increases the degree of dependence to gyro signal when being in acceleration, deceleration or turn condition for a long time.Its prerequisite is that gyro signal drifts about in the period at this section and do not have significant change, and this performance to gyro has proposed higher requirement, but uses high accuracy gyroscope, the price and the weight of meeting increase system, and along with the growth of time, error is accumulation constantly, therefore still can not address this problem.
In addition, utilizing the method measurement attitude of carrier angle of differential GPS or GPS array also is a kind of common method, this method is subjected to the influence that satellite-signal is lost easily, and effect is bad under more situation of blocking, and the position of GPS receiving antenna and distance also can influence the precision of attitude measurement.
Summary of the invention
The present invention has introduced a kind of micro-dynamic carrier attitude measuring apparatus.This device comprises 3 rate gyros (101), 3 magnetometers (102), single-axis accelerometer (103), temperature sensor (104), individual axis velocity sensor (105), analog to digital conversion circuit (106), microprocessor and storer (107), serial communication interface (108), and the corresponding software measuring method.Wherein individual axis velocity sensor (105) is optional device, does not have its system still can operate as normal, helps to improve measuring accuracy but add individual axis velocity sensor (105).
3 rate gyros (101) are connected respectively to each road input end of analog to digital conversion circuit (106) to the output voltage signal of (105) 5 kinds of sensors of individual axis velocity sensor.The control signal of analog to digital conversion circuit (106) is connected with the corresponding port of microprocessor with storer (107) respectively with data-signal.The output port of microprocessor and storer (107) is connected with serial communication unit (108).Microprocessor and storer (107) are sampled to each road sensor signal according to default sampling rate control analog to digital conversion circuit (106), and read corresponding data, handle, and calculate attitude angle, store or output to serial communication unit (108).
Utilize above-mentioned measurement mechanism to measure attitude of carrier information, step is as follows:
The first step is connected measurement mechanism and tested carrier.Definition geographic coordinate system and carrier coordinate system are as shown in Figure 2.Fig. 2 (a) is a geographic coordinate system, three orthogonal axis of N, E, D respectively energized north, east, 3 directions.Fig. 2 (b) is a carrier coordinate system, x, y, three orthogonal axes of z respectively with N, E, when three of D align, each attitude angle is defined as zero, direction is as the criterion with the right-hand rule.The symbolic representation of attitude angle and scope can be defined as: the angle of pitch (θ) scope-90 °~90 °, roll angle (φ) scope-180 °~180 °, 0 °~360 ° of course angle (ψ) scopes.
According to the order of course → pitching → lift-over, the direction cosine matrix R to the carrier coordinate can be expressed as by inertial coordinate:
Wherein s and c are respectively writing a Chinese character in simplified form of function sin and cos.
The present invention constructs a vector v with first row and the 3rd row of direction cosine matrix R,
v=[r
11?r
21?r
31?r
13?r
23?r
33]
R wherein
IjBe illustrated in the capable j row of the i element among the R.
Second step, system initialization.Set the sampling period, the initial value of set condition variable, the initial value of the covariance matrix of state estimation error is measured noise, process noise.With the static placement of measurement mechanism, set the zero point of 3 rate gyros.
The initial value of the covariance matrix of state estimation error can be given arbitrarily, do not influence the convergence of system.The size of measuring noise and process noise then will be according to concrete sensor and applied environment decision.The initial value of state variable need be determined by initial alignment, is zero position such as carrier being positioned over each attitude angle, and the initial value with state variable can be set to v then
0=[1 0000 1].
The 3rd step, the pick-up transducers signal.Gather each road sensing data by microprocessor control analog to digital conversion circuit according to the sampling period, read analog-to-digital result in microprocessor.
In the 4th step, sensing data is compensated and demarcates.According to the sensor temperature family curve sensor sample result is carried out temperature compensation.Because the error that causes is installed also to be needed to demarcate.
The temperature characteristics of 3 rate gyros of sensor (101), 3 magnetometers (102), single-axis accelerometer (103), individual axis velocity sensor (105) will be through measuring in advance, and be kept in the storer.Temperature value during according to sampling, microprocessor reads the temperature characteristics respective value from storer, the sensor data are compensated.Be installed in the center of rotation of carrier if can not guarantee accelerometer (103), then according to current rotational angular velocity with go out to the distance calculation of center of rotation because the centripetal acceleration that rotation causes compensates.
X wherein
AxBe the distance of x axis accelerometer to center of rotation, a
x' be the measured value of accelerometer, a
xFor compensating owing to the accekeration behind rotation and the eccentric acceleration that causes.
In the 5th step, utilize gyro data update mode variable.Utilize the estimated value of the state variable of 3 rate gyro measured values and previous moment can calculate the step forecast estimated value of current time state variable.
More new formula is as follows:
Wherein F is the function of the current angular speed of gyro,
Be the estimated value of kh moment state variable v,
For the step forecast of state variable is estimated.
The 6th step is according to observed reading and constraint condition correction state variable.The constraint condition of 3 magnetometers, single-axis accelerometer and two state variables as measured value, is revised a step predicted value of state variable.
According to measured value and constraint condition, go on foot the quantity of state that forecast obtains to one
Revise:
P(kh+h|kh)=F·P(kh)·F
T+Φ
ev
K(kh+h)=P(kh+h|kh)·C
T[C·P(kh+h|kh)·C
T+Φ
ey]
-1
P(kh+h)=[I-K(kh+h)·C]·P(kh+h|kh)。(3)
In the 7th step, calculate and export current attitude angle.Calculate current attitude angle information according to state variable, and storage or output as required.
V calculates attitude angle by state variable, and concrete formula is as follows:
The angle of pitch:
θ=-arcsin(r
13) (4)
Roll angle:
Course angle:
In the 8th step, jump to the 3rd and go on foot or withdraw from.Jump to the continuation of the 3rd step and measure, perhaps stop to withdraw from.
The 5th step and the 6th goes on foot based on following discrete state spatial model:
v(kh+h)=F·v(kh)+e
v(kh) (7)
y(kh+h)=C·v(kh+h)+e
y(kh+h) (8)
Wherein
y(kh+h)=[m
x(kh+h)?m
y(kh+h)?m
z(kh+h)?a
x(kh+h)?1?1]
T
e
vAnd e
yBe that power spectrum density is Ф
EvAnd Ф
EyThe zero-mean white noise.
With state equation (7), measure equation (8) structure Kalman Filtering for Discrete device, wherein measure Matrix C and constitute by a last moment state variable.Because the measured value y in the equation (8) only contains a
xSo, be not subjected to the influence of y, z direction of principal axis acceleration of motion item.As previously mentioned, if utilize
Replace a
xThe influence of x direction of principal axis acceleration of motion will be eliminated, under high dynamic environment, still high-acruracy survey can be guaranteed like this.
The discrete state space equation is obtained according to sampling period h discretize by following continuity equation.
R gets differential to the direction cosine matrix, and it is as follows to obtain equation:
Wherein
Because v is made up of first row and the 3rd row of R, so can draw by equation (9):
Wherein 0 is 3 * 3 null matrix, e
vBe approximately the zero-mean white noise, its power spectrum density is approximately Ф
Ev
The coordinate representation of acceleration of gravity vector in inertial coordinates system is a
0=[0 0 g]
TMeasured value is a=[a in carrier coordinate system
xa
ya
z]
TThe acceleration of gravity vector is transformed into carrier coordinate system by inertial coordinates system:
a=Ra
0+e
a (11)
E wherein
aFor measuring noise.When carrier is in static state or low dynamic environment, e
aCan be approximated to be power spectrum density is Ф
aThe zero-mean white noise, but when carrier was in high dynamic environment, acceleration measurement and gravitational acceleration component differed greatly, when particularly having long period low frequency movement acceleration, as spiraling e for a long time
aCan not think the zero-mean white noise, but this acceleration is less relatively for the influence of x axle acceleration component.According to equation (6), acceleration at the axial component of carrier coordinate system x is:
a
x=r
13g+e
ax (12)
E wherein
AxCan be approximately power spectrum density is Ф
AxThe zero-mean white noise.If the x direction of principal axis adds speed pickup, use
Replace a in the following formula
xWill help further to improve precision.
The coordinate representation of ground magnetic vector in inertial coordinates system is m
0=[M cos β 0 M sin β]
T, wherein M is the mould (in the statement of back M being omitted as the unit length amount) of ground magnetic vector, β is local earth's magnetic dip angle, revises after ignoring the geomagnetic declination or calculating course angle.Measured value is m=[m in carrier coordinate system
xm
ym
z]
TThe ground magnetic vector is transformed into carrier coordinate system by inertial coordinates system:
m=Rm
0+e
m (13)
E wherein
mFor measuring noise, can be approximately independently zero-mean white noise.
Six elements noticing vector v are not separate, have following two constraint conditions:
With " pseudo-measure equation " of constraint condition (14), constitute the measurement equation of vector v with equation (12), (13) as vector v.Equation contains nonlinear terms in (14), therefore need do linearization process.One in quadratic term in (14) is write in the output equation, is about to equation (12), (13), (14) and is write as following form:
E wherein
y=[e
Mxe
Mye
Mze
Axe
C1e
C2]
T, e
C1, e
C2Be the error of calculation of constraint condition, can be approximately white noise.
Equation (10) and equation (15) will be obtained the discrete system equation that constitutes by equation (7) and equation (8) by sampling period h discretize.
Description of drawings
Fig. 1 is the hardware composition frame chart of this attitude measuring.
Fig. 2 describes inertial coordinates system and carrier coordinate system, wherein (a) inertial coordinates system (NED), (b) carrier coordinate system (xyz).
Fig. 3 is the measuring method process flow diagram.
Fig. 4 is the aircraft coordinate system.
Embodiment
Be example to measure attitude of flight vehicle below, introduce the concrete real-time process of this attitude measuring and method.
Fig. 1 is the hardware composition frame chart of this attitude measuring.3 rate gyros (101), 3 magnetometers (102), single-axis accelerometer (103), temperature sensor (104), individual axis velocity sensor (105), analog to digital conversion circuit (106), microprocessor and storer (107), serial communication interface (108), wherein speed pickup (105) is optional device.3 rate gyros (101), 3 magnetometers (102), single-axis accelerometer devices such as (103) are selected the sensor based on the MEMS technology for use, and volume is little, in light weight.
The output voltage signal of 3 rate gyros (101), 3 magnetometers (102), single-axis accelerometer (103), temperature sensor (104), (105) 5 kinds of sensors of individual axis velocity sensor is connected respectively to each road input end of analog to digital conversion circuit (106).The control signal of analog to digital conversion circuit (106) is connected with the corresponding port of microprocessor with storer (107) respectively with data-signal.The output port of microprocessor and storer (107) is connected with serial communication unit (108).
Fig. 2 describes inertial coordinates system and carrier coordinate system, wherein (a) inertial coordinates system (NED), (b) carrier coordinate system (xyz).
Fig. 3 is the measuring method process flow diagram, and method and the step of utilizing said apparatus to measure attitude of flight vehicle have been described.
(301) this measurement mechanism and aircraft are connected, Fig. 4 illustrated coordinate axis x, y, z respectively with each corresponding axially aligning of aircraft.X axle and aircraft heading are to it, and y axle and aircraft wing direction are to it, and the z axle is vertical with the aircraft wing direction.
(302) system initialization.The setting sampling period is 25Hz; Carrier is positioned over each attitude angle is zero position, the initial value of state variable is set to then; The initial value of the covariance matrix of state estimation error is set at unit matrix arbitrarily; The setting measurement noise, the value of process noise.With the static placement of measurement mechanism, the measured value of setting 3 rate gyros is zero point.
(303) pick-up transducers signal.Gather each road sensing data by microprocessor control analog to digital conversion circuit according to the sampling period, read analog-to-digital result in microprocessor.
(304) temperature compensation and demarcation according to temperature sensor measurements, are read from storer and are demarcated good modified value in advance, and sensing data is carried out temperature compensation.Below used sensor signal be through the data after the temperature compensation.
Measure the distance X of x axis accelerometer to center of rotation
Ax, according to current rotational angular velocity and X
Ax, calculate owing to rotate the centripetal acceleration that causes, compensate.
A wherein
x' be the measured value of accelerometer, a
xFor compensating owing to the accekeration behind rotation and the eccentric acceleration that causes.
(305) utilize gyro data update mode variable.Utilize the estimated value of the state variable of 3 rate gyro measured values and previous moment can calculate the step forecast estimated value of current time state variable.
More new formula is as follows:
Wherein F is the function of the current angular speed of gyro,
Be the estimated value of kh moment state variable v,
For the step forecast of state variable is estimated.
(306) according to observed reading and constraint condition correction state variable.The constraint condition of 3 magnetometers, single-axis accelerometer and two state variables as measured value, is revised a step predicted value of state variable.
According to measured value and constraint condition, go on foot the quantity of state that forecast obtains to one
Revise:
P(kh+h|kh)=F·P(kh)·F
T+Φ
ev
K(kh+h)=P(kh+h|kh)·C
T[C·P(kh+h|kh)·C
T+Φ
ey]
-1
P(kh+h)=[I-K(kh+h)·C]·P(kh+h|kh);
(307) calculate and export current attitude angle.Calculate current attitude angle information according to state variable, and storage or output as required.
V calculates attitude angle by state variable, and concrete formula is as follows:
The angle of pitch:
θ=-arcsin(r
13)
Roll angle:
Course angle:
(308) jumping to the 3rd goes on foot or withdraws from.Jump to the continuation of the 3rd step and measure, perhaps stop to withdraw from.
This device volume is little, in light weight, compare with the conventional inertia guider, less demanding to sensor performance, can select the very little MEMS sensor of volume for use, and pass through special algorithm, therefore eliminated the interference of acceleration of motion effectively, attitude measurement accuracy can be accomplished higher, is specially adapted in the high dynamic environment the dynamic attitude to aircraft, vehicle, boats and ships or other carriers and measures.
Claims (6)
1, a kind of micro-dynamic carrier attitude measuring apparatus, it is characterized in that this device is made up of 3 rate gyros (101), 3 magnetometers (102), single-axis accelerometer (103), temperature sensor (104), analog to digital conversion circuit (106), microprocessor and storer (107), serial communication interface (108);
3 rate gyros (101) are connected respectively to each road input end of analog to digital conversion circuit (106) to the output voltage signal of (104) 4 kinds of sensors of temperature sensor;
The control signal of analog to digital conversion circuit (106) is connected with the corresponding port of microprocessor with storer (107) respectively with data-signal;
The output port of microprocessor and storer (107) is connected with serial communication unit (108).
2, a kind of micro-dynamic carrier attitude measuring apparatus according to claim 1 is characterized in that, this device also has an individual axis velocity sensor (105); The output voltage signal of individual axis velocity sensor (105) is connected to the input end of analog to digital conversion circuit (106).
3, a kind of micro-dynamic carrier attitude measuring apparatus according to claim 1 and 2 is characterized in that, 3 axial organ spares align with each respective coordinates axle in the body axis system respectively.
4, a kind of micro-dynamic carrier attitude measuring apparatus according to claim 1 and 2 is characterized in that, the single-axis accelerometer installation direction is that acceleration of motion is disturbed less one.
5, a kind of micro-dynamic carrier attitude measuring apparatus according to claim 2 is characterized in that, speed pickup and accelerometer are installed in the same way.
6, a kind of micro-dynamic carrier attitude measuring apparatus according to claim 1 carries out a kind of measuring method of attitude of carrier information measurement, it is characterized in that, this measuring method comprises the steps:
The first step: measurement mechanism and tested carrier are connected; Definition geographic coordinate system and carrier coordinate system;
According to the order of course → pitching → lift-over, construct by the direction cosine matrix R of inertial coordinate to the carrier coordinate:
Wherein, s and c are respectively writing a Chinese character in simplified form of function sin and cos;
With first row and vector v of the 3rd row structure of direction cosine matrix R,
v=[r
11?r
21?r
31?r
13?r
23?r
33]
R wherein
IjBe illustrated in the capable j row of the i element among the R;
Second step, system initialization
Set the sampling period, the initial value of set condition variable, the initial value of the covariance matrix of state estimation error is measured noise, process noise;
With the static placement of measurement mechanism, set the zero point of 3 rate gyros;
The 3rd step, the pick-up transducers signal
Microprocessor control analog to digital conversion circuit is gathered each road sensing data according to the sampling period, reads analog-to-digital result in microprocessor;
In the 4th step, sensing data is compensated and demarcates
According to the sensor temperature family curve sensor sample result is carried out temperature compensation;
Because the error that causes is installed also to be needed to demarcate;
In the 5th step, utilize gyro data update mode variable
Utilize the estimated value of the state variable of 3 rate gyro measured values and previous moment to calculate the step forecast estimated value of current time state variable; More new formula is:
Wherein F is the function of the current angular speed of gyro,
Be the estimated value of kh moment state variable v,
For the step forecast of state variable is estimated;
The 6th step is according to observed reading and constraint condition correction state variable
The constraint condition of 3 magnetometers, single-axis accelerometer and two state variables as measured value, is revised a step predicted value of state variable;
According to measured value and constraint condition, go on foot the quantity of state that forecast obtains to one
Revise:
P(kh+h|kh)=F·P(kh)·F
T+Φ
ev
k(kh+h)=P(kh+h|kh)·C
T[C·P(kh+h|kh)·C
T+Ф
ey]
-1
P(kh+h)=[I-K(kh+h)·C]·P(kh+h|kh)
In the 7th step, calculate and export current attitude angle
Calculate current attitude angle information according to state variable, and storage or output as required;
V calculates attitude angle by state variable, and concrete formula is as follows:
The angle of pitch:
θ=-arcsin(r
13)
Roll angle:
Course angle:
In the 8th step, jump to the 3rd and go on foot or withdraw from;
Jump to the continuation of the 3rd step and measure, perhaps stop, withdrawing from, finish measurement.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 200510011763 CN1740746B (en) | 2005-05-23 | 2005-05-23 | Micro-dynamic carrier attitude measuring apparatus and measuring method thereof |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 200510011763 CN1740746B (en) | 2005-05-23 | 2005-05-23 | Micro-dynamic carrier attitude measuring apparatus and measuring method thereof |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1740746A true CN1740746A (en) | 2006-03-01 |
CN1740746B CN1740746B (en) | 2010-08-04 |
Family
ID=36093206
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 200510011763 Expired - Fee Related CN1740746B (en) | 2005-05-23 | 2005-05-23 | Micro-dynamic carrier attitude measuring apparatus and measuring method thereof |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN1740746B (en) |
Cited By (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100423038C (en) * | 2006-10-24 | 2008-10-01 | 中国科学院合肥物质科学研究院 | Sailboard attitude testing device and analysis method |
CN102289306A (en) * | 2011-08-30 | 2011-12-21 | 江苏惠通集团有限责任公司 | Attitude sensing equipment and positioning method thereof as well as method and device for controlling mouse pointer |
CN102313546A (en) * | 2011-04-14 | 2012-01-11 | 南京航空航天大学 | Motion platform gesture sensing method based on polarized electromagnetic wave information chain |
CN102529850A (en) * | 2012-01-16 | 2012-07-04 | 华南理工大学 | Safe state monitoring method of motor vehicle based on wheel load type intelligent sensing |
CN102706360A (en) * | 2012-06-11 | 2012-10-03 | 北京航空航天大学 | Method utilizing optical flow sensors and rate gyroscope to estimate state of air vehicle |
CN102853834A (en) * | 2012-01-09 | 2013-01-02 | 北京信息科技大学 | High-precision scheme of IMU for rotating carrier and denoising method |
CN103090870A (en) * | 2013-01-21 | 2013-05-08 | 西北工业大学 | Spacecraft attitude measurement method based on MEMS (micro-electromechanical systems) sensor |
CN101617199B (en) * | 2007-01-09 | 2013-06-19 | C2Sat通讯股份公司 | Method for compensating drift in a position measuring device |
CN103712598A (en) * | 2013-12-31 | 2014-04-09 | 渤海大学 | Attitude determination system and method of small unmanned aerial vehicle |
CN104567787A (en) * | 2013-10-12 | 2015-04-29 | 北京航天计量测试技术研究所 | Method for calibrating measurement accuracy of dynamic angle measuring system |
CN105352502A (en) * | 2015-11-26 | 2016-02-24 | 湖北三江航天红峰控制有限公司 | Attitude obtaining method of micro-inertia sailing attitude reference system |
CN105573117A (en) * | 2014-11-04 | 2016-05-11 | 霍尼韦尔国际公司 | Configurable inferential sensor for vehicle control systems |
CN105866459A (en) * | 2016-03-25 | 2016-08-17 | 中国人民解放军国防科学技术大学 | Method for estimating restraining angular velocity of gyro-free inertial measurement system |
CN106246168A (en) * | 2016-08-29 | 2016-12-21 | 中国科学院地质与地球物理研究所 | A kind of nearly bit drilling tool attitude measurement while drilling device and measuring method |
US9970288B2 (en) | 2016-08-31 | 2018-05-15 | Institute of geology and geophysics, Chinese Academy of Science | Receiving apparatus for downhole near-bit wireless transmission |
US10025003B1 (en) | 2017-01-19 | 2018-07-17 | Institute Of Geology And Geophysics, Chinese Academy Of Sciences | Calibration method under near-bit wireless short-transmission ground envrionment based on electric field theory |
US10030504B2 (en) | 2016-11-21 | 2018-07-24 | Institute Of Geology And Geophysics, Chinese Academy Of Sciences | Receiving apparatus suitable for azimuthally acoustic logging while drilling |
US10082021B2 (en) | 2016-11-21 | 2018-09-25 | Institute Of Geology And Geophysics, Chinese Academy Of Sciences | Azimuthally acoustic while drilling signal receiving transducer encapsulating apparatus |
US10120795B2 (en) | 2016-09-21 | 2018-11-06 | Institute Of Geology And Geophysics, Chinese Academy Of Sciences | Wear-leveling nandflash memory reading/writing method |
US10317204B2 (en) | 2016-08-31 | 2019-06-11 | Institute Of Geology And Geophysics, Chinese Academy Of Sciences | Near-bit dynamic well deviation angle measurement method and apparatus |
CN110192001A (en) * | 2016-11-09 | 2019-08-30 | 萨兰达有限公司 | Sensor for rotatable element |
US10428646B2 (en) | 2016-08-31 | 2019-10-01 | Institute Of Geology And Geophysics, Chinese Academy Of Sciences | Apparatus for downhole near-bit wireless transmission |
US10578754B2 (en) | 2017-01-19 | 2020-03-03 | Institute Of Geology And Geophysics, Chinese Academy Of Sciences | Sinusoidal excitation method and apparatus for multi-pole acoustic logging while drilling |
US10662764B2 (en) | 2016-08-31 | 2020-05-26 | Institute Of Geology And Geophysics, Chinese Academy Of Sciences | Near-bit constant-power wireless short-distance transmission method and apparatus |
CN111248922A (en) * | 2020-02-11 | 2020-06-09 | 中国科学院半导体研究所 | Human body respiration condition acquisition paste based on accelerometer and gyroscope and preparation method thereof |
CN111856074A (en) * | 2020-06-04 | 2020-10-30 | 中国人民解放军国防科技大学 | Combined type atmospheric data measurement experiment cabin section and flight data measurement method thereof |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102175244A (en) * | 2011-03-16 | 2011-09-07 | 公安部沈阳消防研究所 | In-building person positioning system and positioning method thereof based on inertial sensor |
CN110954103B (en) * | 2019-12-18 | 2022-02-08 | 无锡北微传感科技有限公司 | Method and system for estimating dynamic attitude of vehicle body based on MEMS sensor |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4318063A (en) * | 1979-05-03 | 1982-03-02 | The United States Of America As Represented By The Secretary Of The Air Force | Crystal oscillator compensated for g-sensitivity |
US4696112A (en) * | 1986-09-05 | 1987-09-29 | Condor Pacific Industries, Inc. | Bore hole navigator |
US6008731A (en) * | 1997-07-30 | 1999-12-28 | Union Switch & Signal, Inc. | Detector for sensing motion and direction of a railway vehicle |
US6647352B1 (en) * | 1998-06-05 | 2003-11-11 | Crossbow Technology | Dynamic attitude measurement method and apparatus |
CN1152237C (en) * | 2001-03-30 | 2004-06-02 | 清华大学 | Miniature navigation system based on micro electromechanical techn. |
-
2005
- 2005-05-23 CN CN 200510011763 patent/CN1740746B/en not_active Expired - Fee Related
Cited By (37)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100423038C (en) * | 2006-10-24 | 2008-10-01 | 中国科学院合肥物质科学研究院 | Sailboard attitude testing device and analysis method |
CN101617199B (en) * | 2007-01-09 | 2013-06-19 | C2Sat通讯股份公司 | Method for compensating drift in a position measuring device |
CN102313546B (en) * | 2011-04-14 | 2012-12-26 | 南京航空航天大学 | Motion platform gesture sensing method based on polarized electromagnetic wave information chain |
CN102313546A (en) * | 2011-04-14 | 2012-01-11 | 南京航空航天大学 | Motion platform gesture sensing method based on polarized electromagnetic wave information chain |
CN102289306A (en) * | 2011-08-30 | 2011-12-21 | 江苏惠通集团有限责任公司 | Attitude sensing equipment and positioning method thereof as well as method and device for controlling mouse pointer |
CN102853834A (en) * | 2012-01-09 | 2013-01-02 | 北京信息科技大学 | High-precision scheme of IMU for rotating carrier and denoising method |
CN102529850A (en) * | 2012-01-16 | 2012-07-04 | 华南理工大学 | Safe state monitoring method of motor vehicle based on wheel load type intelligent sensing |
CN102529850B (en) * | 2012-01-16 | 2014-06-25 | 华南理工大学 | Safe state monitoring method of motor vehicle based on wheel load type intelligent sensing |
CN102706360A (en) * | 2012-06-11 | 2012-10-03 | 北京航空航天大学 | Method utilizing optical flow sensors and rate gyroscope to estimate state of air vehicle |
CN102706360B (en) * | 2012-06-11 | 2014-11-26 | 北京航空航天大学 | Method utilizing optical flow sensors and rate gyroscope to estimate state of air vehicle |
CN103090870A (en) * | 2013-01-21 | 2013-05-08 | 西北工业大学 | Spacecraft attitude measurement method based on MEMS (micro-electromechanical systems) sensor |
CN103090870B (en) * | 2013-01-21 | 2015-07-01 | 西北工业大学 | Spacecraft attitude measurement method based on MEMS (micro-electromechanical systems) sensor |
CN104567787B (en) * | 2013-10-12 | 2017-05-17 | 北京航天计量测试技术研究所 | Method for calibrating measurement accuracy of dynamic angle measuring system |
CN104567787A (en) * | 2013-10-12 | 2015-04-29 | 北京航天计量测试技术研究所 | Method for calibrating measurement accuracy of dynamic angle measuring system |
CN103712598A (en) * | 2013-12-31 | 2014-04-09 | 渤海大学 | Attitude determination system and method of small unmanned aerial vehicle |
CN105573117B (en) * | 2014-11-04 | 2021-09-24 | 盖瑞特交通一公司 | Configurable inferential sensor for vehicle control systems |
CN105573117A (en) * | 2014-11-04 | 2016-05-11 | 霍尼韦尔国际公司 | Configurable inferential sensor for vehicle control systems |
CN105352502A (en) * | 2015-11-26 | 2016-02-24 | 湖北三江航天红峰控制有限公司 | Attitude obtaining method of micro-inertia sailing attitude reference system |
CN105866459A (en) * | 2016-03-25 | 2016-08-17 | 中国人民解放军国防科学技术大学 | Method for estimating restraining angular velocity of gyro-free inertial measurement system |
CN105866459B (en) * | 2016-03-25 | 2018-10-26 | 中国人民解放军国防科学技术大学 | Non-gyro inertial measurement system restriction Attitude rate estimator method |
CN106246168A (en) * | 2016-08-29 | 2016-12-21 | 中国科学院地质与地球物理研究所 | A kind of nearly bit drilling tool attitude measurement while drilling device and measuring method |
CN106246168B (en) * | 2016-08-29 | 2017-09-19 | 中国科学院地质与地球物理研究所 | A kind of nearly bit drilling tool posture measurement while drilling device and measuring method |
US10202841B2 (en) | 2016-08-29 | 2019-02-12 | Institute Of Geology And Geophysics, Chinese Academy Of Sciences | Near-bit tool attitude measurement while drilling apparatus and method |
US10317204B2 (en) | 2016-08-31 | 2019-06-11 | Institute Of Geology And Geophysics, Chinese Academy Of Sciences | Near-bit dynamic well deviation angle measurement method and apparatus |
US9970288B2 (en) | 2016-08-31 | 2018-05-15 | Institute of geology and geophysics, Chinese Academy of Science | Receiving apparatus for downhole near-bit wireless transmission |
US10662764B2 (en) | 2016-08-31 | 2020-05-26 | Institute Of Geology And Geophysics, Chinese Academy Of Sciences | Near-bit constant-power wireless short-distance transmission method and apparatus |
US10428646B2 (en) | 2016-08-31 | 2019-10-01 | Institute Of Geology And Geophysics, Chinese Academy Of Sciences | Apparatus for downhole near-bit wireless transmission |
US10120795B2 (en) | 2016-09-21 | 2018-11-06 | Institute Of Geology And Geophysics, Chinese Academy Of Sciences | Wear-leveling nandflash memory reading/writing method |
CN110192001A (en) * | 2016-11-09 | 2019-08-30 | 萨兰达有限公司 | Sensor for rotatable element |
US11402205B2 (en) | 2016-11-09 | 2022-08-02 | Salunda Limited | Sensor for a rotatable element |
US10082021B2 (en) | 2016-11-21 | 2018-09-25 | Institute Of Geology And Geophysics, Chinese Academy Of Sciences | Azimuthally acoustic while drilling signal receiving transducer encapsulating apparatus |
US10030504B2 (en) | 2016-11-21 | 2018-07-24 | Institute Of Geology And Geophysics, Chinese Academy Of Sciences | Receiving apparatus suitable for azimuthally acoustic logging while drilling |
US10578754B2 (en) | 2017-01-19 | 2020-03-03 | Institute Of Geology And Geophysics, Chinese Academy Of Sciences | Sinusoidal excitation method and apparatus for multi-pole acoustic logging while drilling |
US10025003B1 (en) | 2017-01-19 | 2018-07-17 | Institute Of Geology And Geophysics, Chinese Academy Of Sciences | Calibration method under near-bit wireless short-transmission ground envrionment based on electric field theory |
CN111248922A (en) * | 2020-02-11 | 2020-06-09 | 中国科学院半导体研究所 | Human body respiration condition acquisition paste based on accelerometer and gyroscope and preparation method thereof |
CN111248922B (en) * | 2020-02-11 | 2022-05-17 | 中国科学院半导体研究所 | Human body respiration condition acquisition paste based on accelerometer and gyroscope and preparation method thereof |
CN111856074A (en) * | 2020-06-04 | 2020-10-30 | 中国人民解放军国防科技大学 | Combined type atmospheric data measurement experiment cabin section and flight data measurement method thereof |
Also Published As
Publication number | Publication date |
---|---|
CN1740746B (en) | 2010-08-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1740746A (en) | Micro-dynamic carrier attitude measuring apparatus and measuring method thereof | |
CN111878056B (en) | Gyro measurement while drilling system and method | |
CN110031882B (en) | External measurement information compensation method based on SINS/DVL integrated navigation system | |
KR101320035B1 (en) | Location and path-map generation data acquisition and analysis systems | |
CN110792430B (en) | While-drilling inclination measurement method and device based on multi-sensor data fusion | |
CN103245360A (en) | Autocollimation method of carrier aircraft rotating type strapdown inertial navigation system under shaking base | |
CN110007354B (en) | Device and method for measuring flight parameters of semi-aviation transient electromagnetic receiving coil of unmanned aerial vehicle | |
CN104848859A (en) | Three-axis inertial stabilization platform and self-positioning and orientation control method thereof | |
CN109579836B (en) | Indoor pedestrian orientation calibration method based on MEMS inertial navigation | |
CN101044416A (en) | Unmanned airborne vehicle for geophysical surveying | |
US20210348502A1 (en) | Gyro measurement while drilling system and method therefor | |
CN106017452B (en) | Double tops disturbance rejection north finding method | |
WO2006063432A1 (en) | Gyroscopically-oriented survey tool | |
CN102087110B (en) | Miniature underwater moving vehicle autonomous attitude detecting device and method | |
CN104374388A (en) | Flight attitude determining method based on polarized light sensor | |
CN111189474A (en) | Autonomous calibration method of MARG sensor based on MEMS | |
Xing et al. | Quaternion-based Complementary Filter for Aiding in the Self-Alignment of the MEMS IMU | |
CN101581221B (en) | Measurement-while-drilling system | |
CN110736484B (en) | Background magnetic field calibration method based on fusion of gyroscope and magnetic sensor | |
Sun et al. | Coarse alignment based on IMU rotational motion for surface ship | |
Zheng et al. | Compensation for stochastic error of gyros in a dual-axis rotational inertial navigation system | |
CN112729335A (en) | Inertial/starlight combined navigation system calibration method suitable for shaking base | |
CN110988393A (en) | Unmanned aerial vehicle wind speed and direction measurement and correction algorithm based on ultrasonic anemoscope | |
CN116660579A (en) | Wind speed data correction method, system and device | |
CN115727842A (en) | Rapid alignment method and system for unmanned aerial vehicle, computer equipment and storage medium |
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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20100804 Termination date: 20110523 |