CN103674000B - Electronic compass real time calibration algorithm - Google Patents
Electronic compass real time calibration algorithm Download PDFInfo
- Publication number
- CN103674000B CN103674000B CN201310531661.2A CN201310531661A CN103674000B CN 103674000 B CN103674000 B CN 103674000B CN 201310531661 A CN201310531661 A CN 201310531661A CN 103674000 B CN103674000 B CN 103674000B
- Authority
- CN
- China
- Prior art keywords
- vectorial
- electronic compass
- vector
- real time
- value
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C17/00—Compasses; Devices for ascertaining true or magnetic north for navigation or surveying purposes
- G01C17/38—Testing, calibrating, or compensating of compasses
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Traffic Control Systems (AREA)
- Navigation (AREA)
Abstract
The invention belongs to the correlative technology field of aviation aircraft, be specially a kind of electronic compass real time calibration algorithm.The present invention can make aviation aircraft constantly carry out real time calibration to electronic compass when not affecting normal use, thus is user-friendly to, improves the environmental suitability of electronic compass and meet the little requirement of algorithm calculated amount.
Description
Technical field
The invention belongs to the correlative technology field of aviation aircraft, be specially a kind of electronic compass real time calibration algorithm.
Background technology
At present, aviation aircraft obtains self direction, and the parameter generally all relying on electronic compass to provide carrys out identification.
Electronic compass, also referred to as three axle magnetic resistance meters, can be obtained the field strength value of x, y, z three axles in theory, thus extrapolate geomagnetic field intensity and the direction of environment residing for electronic compass by the magnetometer read value of its three axles.But owing to not calibrating zero point by the reading of magnetic resistance meter, using the reasons such as the magnetic interference of the aviation aircraft internal current generation of electronic compass, if its electronic compass three axle reading is directly used in calculate direction, terrestrial magnetic field, be there is huge deviation in the direction, terrestrial magnetic field making to calculate.
It causes the reason of deviation can be divided into three major types principle: the interference of zero deviation, Hard Magnetic and soft magnetism interference.Wherein, zero deviation is droop, is determined by electronic compass self individual difference; Hard Magnetic interference is affected by environment, and under same environmental baseline, this deviation is relatively fixing; Soft magnetism interference causes primarily of the ferromagnetic type objects near electronic compass, both can substantially avoid this to disturb as long as avoid being arranged on by the part having soft magnetism effect near electronic compass.
Summary of the invention
The object of the invention is to overcome above-mentioned deficiency, a kind of electronic compass calibration algorithm is provided, aviation aircraft can be made constantly to carry out real time calibration to electronic compass when not affecting normal use, thus be user-friendly to, improve the environmental suitability of electronic compass and meet the little requirement of algorithm calculated amount.
For realizing above-mentioned technical purpose, scheme provided by the invention is: a kind of electronic compass real time calibration algorithm, comprises the steps.
Step one, setting offset vector F, reads configuration information, obtains to run last gained last time and offset to value and give F, and some f is the corresponding point of vectorial F, and f equals initial point and adds vectorial F.
Step 2, reads the value of real time readouts as vectorial A of the first moment electronic compass, and some a is the corresponding point of vectorial A, and a equals initial point and adds vectorial A.
Step 3, reads the value of real time readouts as vectorial B of the second moment electronic compass, and some b is the corresponding point of vectorial B, and b equals initial point and adds vectorial B.
Step 4, if | A-B| is less than the mould of the surface level projection vector of earth-magnetic field vector, then return step 3; If | A-B| is greater than the mould of the surface level projection vector of earth-magnetic field vector, then enter next step.
Step 5, makes plane D, and D is perpendicular to the straight line C through a, b 2, and D is through the mid point c of a, b 2.
Step 6, some f is to D vertical projection, and obtain some e, initial point is E to the vector of e.
Step 7, makes vectorial F equal vectorial E, makes vectorial A equal vectorial B.
Step 8, skips to step 3, circulation above-mentioned steps.
And described earth-magnetic field vector is that the vector that electronic compass reading is corresponding deducts vectorial F gained.
The present invention can make aviation aircraft constantly carry out real time calibration to electronic compass when not affecting normal use, thus is user-friendly to, improves the environmental suitability of electronic compass and meet the little requirement of algorithm calculated amount.
Embodiment
Below in conjunction with embodiment, the invention will be further described.
The present embodiment provides a kind of electronic compass real time calibration algorithm, if the vector value of zero deviation is M, Hard Magnetic interference vector value is N, and earth-magnetic field vector is T.Vector M can not change substantially, vector N only has the Shi Caihui that change such as the inner airborne equipment of aviation aircraft to change, direction the changing towards with attitudes vibration with aviation aircraft of vector T, the modulus of vector T then changes with latitude and longitude coordinates change residing for aircraft, electronic compass reading R=M+N+T can be established according to this, have | R-M-N|=|T|, that is: all electronic compass readings are all at a periphery, this circle home position is M+N, and radius is | T|.There is this to obtain, to obtain correct magnetic direction, then needing the numerical value (i.e. the position in the circle center of circle) first obtaining M+N.
Appoint and get two electronic compass readings and have R1, R2, there is plane, this plane orthogonal is in R1R2 connecting line, and this plane is through R1R2 mid point, and this plane is through the center of circle of circle.Now establish arbitrfary point, space A1, make A1 perpendicular to this plane projection, obtain the some A2 in plane, the distance that can obtain A2 and the circle center of circle is less than or equal to the distance of A1 to the circle center of circle, that is: Equation f (a, b, c) has A2=f(A1, R1, R2).This is had to obtain, M+N=0 can be established temporarily when value the unknown of M+N, Equation f (a is together substituted into as a and the electronic compass reading newly obtained, b, c), the result of acquisition together brings f(a into as a and the electronic compass reading newly obtained again, b, c), in, continuous iteration and so forth, the result obtained will constantly approach the actual value of M+N.Even if change appears in the value of M+N, after several iterative computation, the result that obtains will constantly approach M+N actual value again.
The concrete scheme that the present embodiment provides comprises the steps.
First need to do following setting.
1. set electronic compass zero deviation and suffered Hard Magnetic to disturb the actual value of sum as vector
, wherein
,
,
for constant.
2. set electronic compass zero deviation and suffered Hard Magnetic to disturb the calculated value of sum as vector
, wherein
,
,
for variable.
3. set aircraft latitude and longitude real time readouts as
,
.
4. known tables, record all parts of the world geomagnetic field intensity and angle between direction, terrestrial magnetic field and surface level, so can establish function
, wherein
for earth-magnetic field vector to surface level projection gained vector field homoemorphism,
for latitude,
for longitude.
5. set electronic compass real time readouts as
, two reading buffer memorys are respectively
with
.
Step one, reads configuration information, obtains and runs last gained and offset to value last time, for
assignment.
Step 2, obtains the real time readouts of the first moment electronic compass
, order
.
Step 3, obtains the real time readouts of the second moment electronic compass
, order
.
Step 4, if
, then step 3 is returned.
Step 5, makes plane, perpendicular to vector
, and through point
,
Can plane equation be obtained thus:
。
Step 6, does straight line, the plane of this linear vertical in step 5, and through point
, can straight-line equation be obtained
, the intersection point of this straight line and step 5 midplane can be asked according to this, by straight-line equation and plane equation simultaneous, have
Namely have
, by gained
the intersection point that value substitution straight-line equation can obtain both is
。
Step 7, can obtain newly according to content in step 6
,
,
value be
, order
,
,
,
,
,
.
Step 8, skips to step 3, circulation above-mentioned steps.
Constantly circulated by above flow process, can make
infinitely approach
, then just can obtain earth-magnetic field vector value accurately
, thus obtain geomagnetic field intensity and azimuth information.
The above is only the preferred embodiment of the present invention; it should be pointed out that for those skilled in the art, under the premise without departing from the principles of the invention; can also make some improvement or distortion, these improve or distortion also should be considered as protection scope of the present invention.
Claims (2)
1. an electronic compass real time calibration algorithm, comprises the steps:
Step one, setting offset vector F, reads configuration information, obtains to run last gained last time and offset to value and give F, and some f is the corresponding point of vectorial F, and f equals initial point and adds vectorial F;
Step 2, reads the value of real time readouts as vectorial A of the first moment electronic compass, and some a is the corresponding point of vectorial A, and a equals initial point and adds vectorial A;
Step 3, reads the value of real time readouts as vectorial B of the second moment electronic compass, and some b is the corresponding point of vectorial B, and b equals initial point and adds vectorial B;
Step 4, if | A-B| is less than the mould of the surface level projection vector of earth-magnetic field vector, then return step 3; If | A-B| is greater than the mould of the surface level projection vector of earth-magnetic field vector, then enter next step;
Step 5, makes plane D, and D is perpendicular to the straight line C through a, b 2, and D is through the mid point c of a, b 2;
Step 6, some f is to D vertical projection, and obtain some e, initial point to the vector of e is E;
Step 7, makes vectorial F equal vectorial E, makes vectorial A equal vectorial B;
Step 8, skips to step 3, circulation above-mentioned steps.
2. a kind of electronic compass real time calibration algorithm according to claim 1, is characterized in that: described earth-magnetic field vector is that the vector that electronic compass reading is corresponding deducts vectorial F gained.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310531661.2A CN103674000B (en) | 2013-11-01 | 2013-11-01 | Electronic compass real time calibration algorithm |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310531661.2A CN103674000B (en) | 2013-11-01 | 2013-11-01 | Electronic compass real time calibration algorithm |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103674000A CN103674000A (en) | 2014-03-26 |
CN103674000B true CN103674000B (en) | 2015-10-21 |
Family
ID=50312272
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310531661.2A Expired - Fee Related CN103674000B (en) | 2013-11-01 | 2013-11-01 | Electronic compass real time calibration algorithm |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103674000B (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107063218B (en) * | 2017-04-21 | 2020-09-01 | 北京京东尚科信息技术有限公司 | Data processing method and device |
CN107144271A (en) * | 2017-05-16 | 2017-09-08 | 北京京东尚科信息技术有限公司 | Data processing method and device |
CN109059960B (en) * | 2018-07-18 | 2021-08-31 | 福州大学 | Calibration method of three-dimensional electronic compass |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0436794A1 (en) * | 1989-11-07 | 1991-07-17 | State of Israel Ministry of Defense Armaments Development Authority, Rafael | Magnetic compass |
CN1624425A (en) * | 2003-12-03 | 2005-06-08 | 三星电机株式会社 | Automatic calibration method for use in electronic compass |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6513252B1 (en) * | 1999-04-08 | 2003-02-04 | Donnelly Corporation | Vehicle compass compensation |
JP4528342B2 (en) * | 2006-03-30 | 2010-08-18 | 京セラ株式会社 | Portable electronic device and geomagnetic sensor calibration method |
-
2013
- 2013-11-01 CN CN201310531661.2A patent/CN103674000B/en not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0436794A1 (en) * | 1989-11-07 | 1991-07-17 | State of Israel Ministry of Defense Armaments Development Authority, Rafael | Magnetic compass |
CN1624425A (en) * | 2003-12-03 | 2005-06-08 | 三星电机株式会社 | Automatic calibration method for use in electronic compass |
Non-Patent Citations (1)
Title |
---|
《高精度磁电子罗盘的研制》;李希胜等;《传感技术学报》;20061231;第19卷(第6期);第2441-2444页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103674000A (en) | 2014-03-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104406610B (en) | A kind of magnetometer real time correction device and method | |
CN102252689A (en) | Electronic compass calibration method based on magnetic sensor | |
CN106323334B (en) | A kind of magnetometer calibration method based on particle group optimizing | |
CN105737858B (en) | A kind of Airborne Inertial Navigation System attitude parameter calibration method and device | |
CN110146839A (en) | A kind of mobile platform magnetic gradient tensor system compensation method | |
CN106959105B (en) | Method for calibrating compass and device | |
CN109725360A (en) | One-point positioning method based on magnetic gradient tensor invariant | |
CN104697523A (en) | Inertia/terrestrial magnetism matching and positioning method based on iterative computation | |
CN104459728B (en) | A kind of magnetic declination calibration steps positioned based on GNSS | |
CN104913777A (en) | Error compensation algorithm of geomagnetic sensor application system | |
CN110736484B (en) | Background magnetic field calibration method based on fusion of gyroscope and magnetic sensor | |
CN104142485A (en) | Magnetic sensor calibration method and device | |
CN103674000B (en) | Electronic compass real time calibration algorithm | |
CN109725361B (en) | Magnetic target positioning method based on invariant of magnetic gradient tensor | |
CN114674301A (en) | Electronic compass active compensation method and system under strong magnetic interference environment | |
CN104748714B (en) | A kind of spaceborne rope net deployable antenna surface accuracy measuring method | |
Chen et al. | A novel calibration method for tri-axial magnetometers based on an expanded error model and a two-step total least square algorithm | |
CN103913179A (en) | Ship-based single-axial rotation modulation strapdown inertial navigation attitude angle correction method | |
CN109031439A (en) | A kind of geomagnetic diurnal variations numerical value based on difference of latitude and distance determines method and system | |
Jiang et al. | Scalar calibration of aeromagnetic data using BPANN and LS algorithms based on fixed-wing UAV platform | |
Liu et al. | Simplified ellipsoid fitting-based magnetometer calibration for pedestrian dead reckoning | |
CN105865444A (en) | Affine transformation-based inertial/geomagnetic matching iterative positioning method | |
CN108061922B (en) | A kind of distribution magnetic anomaly detection system dynamic magnetic compensation method | |
Ji et al. | Integrated compensation method of three-axis magnetometer in geomagnetic navigation | |
CN108200547A (en) | Rigid body localization method based on measurement distance |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20151021 Termination date: 20211101 |
|
CF01 | Termination of patent right due to non-payment of annual fee |