CN101241009B - Magneto- electronic compass error compensation method - Google Patents

Magneto- electronic compass error compensation method Download PDF

Info

Publication number
CN101241009B
CN101241009B CN2007103045358A CN200710304535A CN101241009B CN 101241009 B CN101241009 B CN 101241009B CN 2007103045358 A CN2007103045358 A CN 2007103045358A CN 200710304535 A CN200710304535 A CN 200710304535A CN 101241009 B CN101241009 B CN 101241009B
Authority
CN
China
Prior art keywords
angle
magnetic
alpha
compensation
compass
Prior art date
Application number
CN2007103045358A
Other languages
Chinese (zh)
Other versions
CN101241009A (en
Inventor
李希胜
王磊
舒雄鹰
王立锦
Original Assignee
北京科技大学
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 北京科技大学 filed Critical 北京科技大学
Priority to CN2007103045358A priority Critical patent/CN101241009B/en
Publication of CN101241009A publication Critical patent/CN101241009A/en
Application granted granted Critical
Publication of CN101241009B publication Critical patent/CN101241009B/en

Links

Abstract

The invention provides an electronic magnetic compass error compensation method which bases on deformation round distribution and cyclical assumptions. The measurement error caused by around ferromagnetic materials and other factors must be compensated to enhance electronic magnetic compass measurement accuracy. At present, domestic and foreign popular methods of compensation are based on assumptions which relate to the effect of the environmental interference magnetic field causing the magnetic field strength level projection component distribution from round into oval. But because of more and complicated influence factors, in fact the magnetic field strength level projection component distribution can not be guarantee as perfect oval. In view of this situation, the invention provides error compensation algorithmic method which bases on deformation round distribution and cyclical assumptions, the objective is to compensate the measurement error caused by around ferromagnetic materials and other factors, and achieve magnetic azimuth high-precision measurement.

Description

A kind of magneto-electronic compass error compensation method

Technical field

Patent of the present invention relates to the sub-compass of magnetoelectricity that utilizes the terrestrial magnetic field to realize orientating function, has provided the measurement Error Compensation method of the sub-compass of a kind of magnetoelectricity.

Background technology

After the sub-compass of magnetoelectricity was fixed on the carrier, the magnetic azinuth measurement result can be subjected to the influence of ferromagnetic material on every side.Disturbing magnetic field is made up of hard iron magnetic field and soft iron magnetic field two parts.When having disturbing magnetic field on every side, the distribution of magnetic field intensity horizontal projection component no longer is positive circle, and distortion (influence in soft iron magnetic field) has taken place its shape, and skew (influence in hard iron magnetic field) has also taken place its center relative coordinate initial point.The existence of above factor causes the sub-compass of magnetoelectricity measuring error to occur in the magnetic azinuth measuring process.Distribution that the popular backoff algorithm influence that is based on environmental interference magnetic field both at home and abroad causes magnetic field intensity horizontal projection component becomes oval-shaped hypothesis by positive circle.Promptly fit magnetic field intensity horizontal projection component and distribute, ask for relevant coefficient, in use carry out the compensation of magnetic-field measurement value according to relevant coefficient according to ellipse.Above backoff algorithm is supposed the ellipse that is distributed as rule of magnetic field intensity horizontal projection component, and in fact is difficult to guarantee the strictness establishment of this condition owing to the complicacy of concrete condition, thereby has limited the raising of measuring accuracy.

Summary of the invention

The present invention proposes based on being out of shape the Error Compensation Algorithm that circle distributes and periodically supposes, with a kind of method that can be used for the magneto-electronic compass error compensation of caused by factors measuring error such as ferromagnetic material around the compensation.

The present invention is a kind of method that can be used for magneto-electronic compass error compensation, when the sub-compass of magnetoelectricity uniform rotation multi-turn in the environment around, and the azimuth angle alpha of azimuth angle alpha, the not compensated that records so m, azimuth measurement error Δ α is with the generating period variation, can be expressed as the form of the periodic function of the sub-compass corner of magnetoelectricity β, [0,2 π] in the scope, the form of the function generate convergent series behind periodic extension can be similar to this function with the algebraic sum of finite term, promptly

Δα = a 0 + Σ n = 1 j ( a n cos nα + b n sin nα )

In the following formula, j is the finite value positive integer, and a0, an, bn are coefficient, n=1 wherein, and 2,3 ...,

Calibration process comprises following each step:

(1) the no magnetic turntable of sub-compass of magnetoelectricity and caliberating device is equipped with in rotation, makes it point to positive north, determines that this moment, the position angle was 0 degree; (2) according to the actual measurement magnetic field data and as data computation position angle, geomagnetic declination approximate value α m

(3) the corner data α computer azimuth angle error delta α=α that exports according to caliberating device m-α;

(4) rotate no magnetic turntable and make the position angle increase a certain angle, repeat second step and the 3rd step, the sub-compass azimuth angle error of the magnetoelectricity when obtaining the respective party parallactic angle.Reach 360 degree (identical) until the position angle with 0 degree;

(5) calculate compensation formula according to azimuth angle error

Δα = a 0 + Σ n = 1 j ( a n cos nα + b n sin nα )

In penalty coefficient a0, an, bn, n=1 wherein, 2,3 ...;

During use, the compensation formula coefficient calculations that will obtain according to demarcation more on the occasion of, its step is as follows:

(1) according to the actual measurement magnetic field data and as data computation position angle, geomagnetic declination approximate value α m

(2) according to compensation formula Calculate more on the occasion of-Δ α, at first use position angle approximate value α during calculating mReplace position angle actual value α, then carry out repeatedly iteration as the need degree of precision;

(3) calculate compensation back bearing measured value α '=α m+ (Δ α).

Advantage of the present invention:

Proposed with caused by factors measuring error such as ferromagnetic materials around the compensation, to realize the high-acruracy survey of magnetic azinuth based on being out of shape the Error Compensation Algorithm that circle distributes and periodically supposes.

Embodiment

Provide an instantiation that adopts patent of the present invention to realize the magneto-electronic compass error compensation below.

In this example, be the magnetic azimuth angle transducer with high performance three magnetic resistance sensors, measure the terrestrial magnetic field component on carrier x axle, y axle, z axle respectively With the double-axel acceleration sensor is the angle of pitch (θ) and roll angle (γ) survey sensor.The information of comprehensive above a plurality of sensors, the high-acruracy survey of realization magnetic azinuth.

The realization of the magneto-electronic compass error compensation method that patent of the present invention proposes is divided into two stages, i.e. the error compensation stage of the coefficient calculations stage of calibration process and use.In calibration process, calculate the compensation formula coefficient according to measurement error value; In use, utilize compensation formula and relevant coefficient calculations more on the occasion of, realize error compensation.

When the sub-compass of magnetoelectricity uniform rotation multi-turn in the environment around, the azimuth angle alpha of the not compensated that records so m, azimuth measurement error Δ α is with the generating period variation, can be expressed as the form of the periodic function of azimuth angle alpha, [0,2 π] in the scope, the form of the function generate convergent series behind periodic extension can be similar to this function with the algebraic sum of finite term, promptly

Δα = a 0 + Σ n = 1 j ( a n cos nα + b n sin nα )

In the following formula, j is the finite value positive integer, and a0, an, bn are coefficient, n=1 wherein, and 2,3 ...;

Timing signal is installed in the sub-compass of magnetoelectricity on the no magnetic turntable, and the turntable corner is measured by high-precision photoelectric encoder.Calibration process is made up of following steps.

The first step is rotated no magnetic turntable, makes its forward direction refer to north, and the output of the photoelectric encoder that resets determines that the position angle is 0 degree.

Second step is according to measured data computer azimuth angle approximate value.At first, according to following formula calculating magnetic field strength level projection H x, H y

In the formula, For Magnetic Sensor in θ ≠ 0,3 geomagnetic field components that record under γ ≠ 0 state.

Then, according to H x, H yCalculate the magnetic azinuth approximate value

α ~ M = arctan ( H y H x )

At last, according to the magnetic declination value D computer azimuth angle approximate value of locality

α m = α ~ M + D

The 3rd step is according to the corner data α computer azimuth angle error of photoelectric encoder output.

Δα=α m

The 4th step, rotate no magnetic turntable and make the position angle increase by 15 degree, repeat second step and the 3rd step, the sub-compass azimuth angle error of the magnetoelectricity when obtaining the respective party parallactic angle.Reach 360 degree (identical) until the position angle with 0 degree.

In the 5th step, calculate compensation formula according to azimuth angle error

Δα = a 0 + Σ n = 1 j ( a n cos nα + b n sin nα )

In penalty coefficient a0, an, bn, n=1 wherein, 2,3 ....

In the above calibration process, we only need to rotate no magnetic turntable and get final product to respective angles.

During use, the compensation formula coefficient calculations that will obtain according to demarcation more on the occasion of, its step is as follows.

The first step is according to measured data computer azimuth angle approximate value.At first, according to following public affairs

Formula calculating magnetic field strength level projection H x, H y

In the formula, For Magnetic Sensor in θ ≠ 0,3 geomagnetic field components that record under γ ≠ 0 state.

Then, according to H x, H yCalculate the magnetic azinuth approximate value

α ~ M = arctan ( H y H x )

At last, according to the magnetic declination value D computer azimuth angle approximate value of locality

α m = α ~ M + D

Second step is according to compensation formula Calculate more on the occasion of-Δ α, at first use position angle approximate value α during calculating mReplace position angle actual value α, then repeatedly change as the need degree of precision.

In the 3rd step, calculate compensation back bearing measured value α '=α m+ (Δ α).

Claims (2)

1. method that is used for magneto-electronic compass error compensation is characterized in that when the sub-compass of magnetoelectricity uniform rotation multi-turn in the environment around, the azimuth angle alpha of the not compensated that records so m, azimuth measurement error Δ α is with the generating period variation, all is expressed as the form of the periodic function of azimuth angle alpha, in [0,2 π] scope, the form of the function generate convergent series behind periodic extension is similar to this function with the algebraic sum of finite term, promptly
Δα = a 0 + Σ n = 1 j ( a n cos nα + b n sin nα )
In the following formula, j is the finite value positive integer, a 0, a n, b nBe coefficient, n=1 wherein, 2,3 ...; Calibration process comprises following each step:
(1) the no magnetic turntable of sub-compass of magnetoelectricity and caliberating device is equipped with in rotation, makes it point to positive north, determines that this moment, the position angle was 0 degree;
(2) according to the actual measurement magnetic field data and as data computation position angle, geomagnetic declination approximate value α m
(3) the corner data α computer azimuth angle error delta α=α that exports according to caliberating device m-α;
(4) rotate no magnetic turntable and make the position angle increase a certain angle, repeat second step and the 3rd step, the sub-compass azimuth angle error of the magnetoelectricity when obtaining the respective party parallactic angle reaches 360 degree until the position angle;
(5) calculate compensation formula according to azimuth angle error
Δα = a 0 + Σ n = 1 j ( a n cos nα + b n sin nα )
In penalty coefficient a 0, a n, b n, n=1 wherein, 2,3 ...;
During use, the compensation formula coefficient calculations that will obtain according to demarcation more on the occasion of, its step is as follows:
(1) according to the actual measurement magnetic field data and as data computation position angle, geomagnetic declination approximate value α m
(2) according to compensation formula Calculate more on the occasion of-Δ α, at first use position angle approximate value α during calculating mReplace position angle actual value α, then carry out repeatedly iteration as the need degree of precision;
(3) calculate compensation back bearing measured value α '=α m+ (Δ α).
2. a kind of method that is used for the magneto-electronic compass error compensation according to claim 1 is characterized in that the sub-compass of described magnetoelectricity is based on the measuring sensor of fluxgate technology or anisotropic magnetoresistance technology or giant magnetoresistance technology or giant magnetic impedance technology.
CN2007103045358A 2007-12-28 2007-12-28 Magneto- electronic compass error compensation method CN101241009B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2007103045358A CN101241009B (en) 2007-12-28 2007-12-28 Magneto- electronic compass error compensation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2007103045358A CN101241009B (en) 2007-12-28 2007-12-28 Magneto- electronic compass error compensation method

Publications (2)

Publication Number Publication Date
CN101241009A CN101241009A (en) 2008-08-13
CN101241009B true CN101241009B (en) 2010-06-09

Family

ID=39932712

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2007103045358A CN101241009B (en) 2007-12-28 2007-12-28 Magneto- electronic compass error compensation method

Country Status (1)

Country Link
CN (1) CN101241009B (en)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101587132B (en) * 2009-06-26 2011-07-20 中国科学院合肥物质科学研究院 Field weakening direction sensor calibration method
US8374816B2 (en) * 2010-06-08 2013-02-12 Honeywell International Inc. Automatic data collection algorithm for 3D magnetic field calibration with reduced memory requirements
CN102313543B (en) * 2011-07-11 2013-07-17 上海大学 Magnetic azimuth measuring system based on giant magneto-resistance sensor, measurement method and perpendicular compensation method
CN102506837A (en) * 2011-10-10 2012-06-20 海能达通信股份有限公司 Terminal positioning method and terminal with electronic compass
CN104111066B (en) * 2014-07-29 2016-06-15 中船重工中南装备有限责任公司 A kind of electronic compass testing-calibrating instrument
CN104792301A (en) * 2015-04-22 2015-07-22 华中科技大学 Method and device for correcting azimuth of borehole television probe under ferromagnetic interference
CN104965102B (en) * 2015-05-25 2018-04-27 厦门大学 A kind of method for effectively improving acoustic Doppler fluid velocity profile instrument flow directional detection precision
CN105675905B (en) * 2016-01-28 2018-10-09 北京理工大学 A kind of rotary bullet tachometric survey error compensating method based on Geomagnetism Information
CN105865492B (en) * 2016-05-31 2019-02-05 清华大学 The online error compensating method of two axis magnetometers and system
CN106052633B (en) * 2016-07-01 2018-10-30 西安华讯天基通信技术有限公司 The azimuthal device and method in sector in a kind of determining radio station
CN107340001B (en) * 2017-05-23 2020-02-28 中国人民解放军军械工程学院 Geomagnetic measurement error compensation test device
CN108398124B (en) * 2018-02-05 2020-10-16 无锡北微传感科技有限公司 Electronic compass calibration method

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1181808A (en) * 1996-03-13 1998-05-13 莱卡公开股份有限公司 Process for determining the direction of the earth's magnetic field
CN1216104A (en) * 1997-02-10 1999-05-05 莱卡地球系统公开股份有限公司 Method for determining correction parameters
CN1323979A (en) * 2000-05-12 2001-11-28 精工电子有限公司 Calibration mechanism, electron compass having said mechanism and electron clock
CN1627040A (en) * 2003-12-13 2005-06-15 三星电机株式会社 Method for searching for dip angle using tilt-compensated electric compass

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1181808A (en) * 1996-03-13 1998-05-13 莱卡公开股份有限公司 Process for determining the direction of the earth's magnetic field
CN1216104A (en) * 1997-02-10 1999-05-05 莱卡地球系统公开股份有限公司 Method for determining correction parameters
CN1323979A (en) * 2000-05-12 2001-11-28 精工电子有限公司 Calibration mechanism, electron compass having said mechanism and electron clock
CN1627040A (en) * 2003-12-13 2005-06-15 三星电机株式会社 Method for searching for dip angle using tilt-compensated electric compass

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
JP特开2000-9468A 2000.01.14
JP特开2005-106569A 2005.04.21

Also Published As

Publication number Publication date
CN101241009A (en) 2008-08-13

Similar Documents

Publication Publication Date Title
CN101672665B (en) Encoder output signal correction apparatus and method
CN103337066B (en) 3D obtains the calibration steps of system
CN1323282C (en) Apparatus and method of compensating for the attitude change of an earth magnetic sensor
CN101551240B (en) Large-scale gear measuring method based on laser tracking technology
US20030079358A1 (en) Method of an apparatus for correction of coordinate measurement errors due to vibrations in coordinate measuring machines (cmms)
Van Camp et al. Uncertainty of absolute gravity measurements
EP3249347B1 (en) Magnetoresistive angle sensor and corresponding strong magnetic field error correction and calibration methods
CN101694373B (en) Antenna deformation measuring method
Jan et al. Large variability of the Kuroshio at 23.75 N east of Taiwan
CN103630705B (en) A kind of measuring method utilizing solid state two dimensional wind speed and direction measuring instrument
CN102818564B (en) Calibration method of three-dimensional electronic compass
CN102313543B (en) Magnetic azimuth measuring system based on giant magneto-resistance sensor, measurement method and perpendicular compensation method
CN104121905A (en) Course angle obtaining method based on inertial sensor
US7834616B2 (en) Magnetic speed, direction, and/or movement extent sensor
CN100549731C (en) Demodulator probe secondary localization method
CN101265804B (en) Well drilling high precision gradient meter sensor perpendicular installation error compensation process
Cea et al. Bathymetric error estimation for the calibration and validation of estuarine hydrodynamic models
WO2010088119A3 (en) Reducing error contributions to gyroscopic measurements from a wellbore survey system
US20100176799A1 (en) Sensor including two code rings and a magnetic field sensor between the code rings
CN103940332B (en) Magnetic grating displacement transducer based on Hall magnetic sensitive element array
CN102607521A (en) Accelerometer based tilt sensor and method for using same
CN102171628A (en) Pointer with motion sensing resolved by data merging
Duffy et al. Application of spatial cross correlation to detection of migration of submarine sand dunes
CN101975872B (en) Method for calibrating zero offset of quartz flexible accelerometer component
CN103925933B (en) A kind of multi-turn absolute magnetic encoder

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: 20100609

Termination date: 20161228