JPH0340323B2 - - Google Patents

Info

Publication number
JPH0340323B2
JPH0340323B2 JP56086729A JP8672981A JPH0340323B2 JP H0340323 B2 JPH0340323 B2 JP H0340323B2 JP 56086729 A JP56086729 A JP 56086729A JP 8672981 A JP8672981 A JP 8672981A JP H0340323 B2 JPH0340323 B2 JP H0340323B2
Authority
JP
Japan
Prior art keywords
magnetic
magnetometer
heading
magnetic field
vehicle
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 - Lifetime
Application number
JP56086729A
Other languages
English (en)
Other versions
JPS5752813A (en
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 filed Critical
Publication of JPS5752813A publication Critical patent/JPS5752813A/ja
Publication of JPH0340323B2 publication Critical patent/JPH0340323B2/ja
Granted legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C17/00Compasses; Devices for ascertaining true or magnetic north for navigation or surveying purposes
    • G01C17/38Testing, calibrating, or compensating of compasses

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Measuring Magnetic Variables (AREA)
  • Navigation (AREA)
  • Arrangements For Transmission Of Measured Signals (AREA)
  • Magnetic Ceramics (AREA)

Description

【発明の詳細な説明】
本発明は、乗物(ビークル)の移動中、磁北に
対する乗物の経路もしくは向首方向を任意時点で
指示するように適応された装置に関する。ここで
乗物とは例えば陸上車輌または航空機とすること
ができる。 この種の周知の装置においては、センサとして
該センサが取り付けられている乗物の変位に対し
連続的に地磁界の変化に応答する磁力計を使用す
ることが知られている。実際においては、磁力計
の直交軸線OXおよびOYは乗物の縦軸線ならび
に該縦軸線に対して垂直な軸線にそれぞれ一致す
るように配置されており、これら2つの軸は、乗
物が陸上車輌である場合には水平もしくは実質的
に水平とすることができる乗物の平面P(第1図)
内に位置している。 このような条件下において、磁力計は地磁界H→
の水平成分Hx→およびHy→に対し比例する電圧を発
生することを可能にするものである。理論的に
は、乗物が磁界H中を移動している場合には、該
ベクトルH→の水平面上への射影は端Mが磁力計の
軸線OXおよびOYの原点である中心Oの円上を
移動するベクトルHh→である。この条件下で、乗
物の方向および磁北間に形成される角度である実
際の磁針向首方向φならびに基準軸に対する点M
の磁力計によつて与えられる座標から演繹される
角度ψ′は等しい。 実際上は搬送乗物の所与の磁性形態において、
点Mは乗物のそれ自体を中心とする360゜の回転中
円を描かず楕円を描く。この現象は、地磁界中に
おける乗物の運動過程中角度ψおよびψ′間に程度
の差こそあれ発生しかつ可変であるシフト即ち遷
移に因るものである。このような向首方向誤差
は、乗物の強磁性構造の近接配置ならびに乗物内
を循環する電流によつて誘起される磁界から生ず
る地磁界の振幅および方向における擾乱に起因す
るものである。 このようにして磁力計により与えられる擾乱さ
れた向首方向データを補正するために、一連のよ
く定義されたコンパス主要方位点に対し実際の向
首方向ψと、向首方向測定装置により計算された
向首方向ψ′とを比較して、その比較結果から全て
のコンパス主要点に対しψ′=ψを得るために必要
な補償量を演繹することからなる方法が知られて
いる。 しかしながらこの方法においては時間が長くか
かり、しかも運行レベルにおいて制限が課せられ
る非常に正確な向首方向測定が要求される。 乗物が陸上車輌であつて、水平面上に維持され
ている場合には該乗物は、照準管により非常に正
確な仕方で、360度にわたり分布された12の向
首方向に配置される。このためには、車輌の向首
方向決定装置の磁気補償中該車輌が移動している
陸地のコンパスによる厳密なマーキングが必要と
なる。航空機の場合にも、伝達比ジヤイロを用い
て類似の手順が要求される。 補正を決定するためのこの方法においては専問
家の集まつたチームの場合でも数時間を要するこ
とがあり得る。 本発明による補償装置の目的とするところは、
補償係数の決定を相当に単純化しかつ相当量の時
間を節減することにある。さらに本発明は、計算
後に、従来の方法と比較して一層改善された精度
で実際の磁気向首方向ψを与える係数を得ること
にある。 この目的で本発明は、乗物に搭載された磁針向
首方向を決定するための装置の測定に影響する磁
気擾乱を補償するための装置であつて、基準軸に
沿う擾乱された地磁界の成分に比例するアナログ
電圧を発生し、その場合にいわゆる補償係数を磁
力計と組合わせて設けられている計算機モジユー
ルのメモリに導入して磁力計により発生される情
報を補正し実際の向首方向を求めることを可能に
する装置に関するものである。なお補償係数は磁
力計から供給される擾乱地磁界電圧から直接的
に、即ち向首方向基準を参照することなしに決定
されるものである。 陸上車輌の場合には、車輌はそれ程大きな精度
を伴なわずに、360度を実質的に等分した任意8
つの向首方向に配位される。各位置において、磁
力計(第2図)により供給される擾乱電圧Uxpお
よびUypは、擾乱された地磁界ベクトルの端
M′の水平面内の座標XpおよびYpを表わす。こ
れら8つの点は車輌のそれ自体を中心とする1完
全回転中に点M′によつて描かれる楕円上に在る
ものと考えてもよい。計算機はこの楕円のパラメ
ータを求めて磁力計軸OXおよびOY上に心出し
された理論円にこの楕円を変換せしめる補償係数
を算出する。なお該理論円は擾乱が存在しなかつ
た場合には点Mの位置に対応するものである。 向首方向決定装置と固定的に連結していても、
あるいは連結していなくてもよいこの計算機は、
楕円のパラメータを定めて、擾乱もしくは妨害を
受けた磁針向首方向ψ′から実際の磁針向首方向ψ
に移すために、磁力計から供給される電圧Uxp,
Uypに印加することが必要とされる補正係数を計
算する。陸上車輌が水平位でない場合には、地磁
界の垂直成分Hvは、横揺れ角ならびに縦揺れ
角θの値に従い地磁界の唯一の水平成分Hhの成
分HxおよびHyに影響を及ぼす。 後に述べた補正は、センサ、補償モジユールお
よび表示装置を有する向首方向決定装置の段階で
行なわれる。計算機によつて求められた補償係数
は補償過程中向首方向決定装置の補償モジユール
のメモリに装入される。 このモジユールの役割りは、磁力計によつて供
給される電圧Uxp=kHxpおよびUyp=kHypな
らびに傾斜計から供給される電圧から、補償され
た向首方向即ち実際の向首方向ψの正弦および余
弦に対応するアナログ電圧形態の電気信号を発生
することにある。これら信号は回転目盛カードお
よび固定インデツクスによつて構成される表示モ
ジユールに送られて、磁化に対する乗物の実際の
向首方向ψの表示を可能にする。 磁力計の2つの軸OXおよびOYは第1図の面
P内に位置する。軸OZは、面Pが水平であると
した場合に垂線を表わす。このような条件下で、
地理上の位置もしくは場所に依存して変化する角
度だけ水平線に対し傾斜している磁界H→は軸
OX,OY,OZに沿つて3つのベクトル成分に分
解することができる。即ち垂直成分Hv→と、面P
におけるH→の投影であるHh→の成分Hhx―→hxおよび
Hhy―→である。 乗物の近傍に地磁界Hの外乱もしくは擾乱が存
在しないとした場合には、角度XOMは実際の磁
針向首方向ψを表わす。磁界H→が擾乱されると、
磁力計は点M′を、XOMが偽磁針向首方向に対応
する角度ψ′となるように定める面P内の成分Hhp―→
hpを有する磁界Hp→を測定する。そこで磁力計は
2つの擾乱もしくは外乱されたアナログ電圧Uxp
=kHxpおよびUyp=kHypを発生する。ここで
kは磁力計の換算係数に対応する係数である。 航空機の場合には補償は次のようにして行なわ
れる。 即ち2つの相続く測定間において地磁界のモジ
ユールが一定であると仮定し、擾乱もしくは外乱
を受けた地磁界の空間におけるモジユールの展開
(楕円体)を検査し、反復的な仕方で、地磁界の
モジユールの展開を楕円体から磁力計の軸に心出
しされた理論的球体に戻して擾乱された測定値と
実際もしくは現実の測定値との間の対応を定める
ことを可能にする補償係数を計算することにより
行なわれる。 陸上車輌の強磁性が地磁界に何ら影響を有しな
いとした場合には、この地磁界の水平成分Hhの
端Mは磁力計の軸OXおよびOYの原点である中
心Oを有する円C1を描くことになろう。 しかしながら現実には製造の過程中に磁化され
たり、あるいは電気的な起源に由る磁界の影響下
で磁化された車輌もしくは乗物の強磁性構造に起
因する永久磁界によるいわゆる「硬鉄」擾乱もし
くは外乱が存在する。これら永久磁界は車輌もし
くは乗物に結合されて該乗物と共に回転する。 この結果磁力計はその軸の原点に対してXo,
Yoだけ偏差している中心O′を有する円C2上を移
動する端M″を有するベクトルの成分を測定する
ことになる。 磁力計に結合された空間および基準マークにお
いて、これら永久磁界はベクトルH→o=Xoi→+
Yoj→+Zov→によつて表わすことができよう。上
式中i,jおよびvは陸上車輌に取り付けられた
三軸磁力計のマークされた直交軸を表わす。 磁力計の測定結果に影響を与える第2の擾乱も
しくは外乱がある。これはいわゆる「軟鉄」外乱
と称されるものであつて、地磁界の作用により誘
起される磁界に起因するものである。これら磁界
は車輌と共に回転せず、その極性および振幅は車
輌に対する地磁界の方向および振幅の関数として
変化する。 これら誘起される磁界Hi→は次のような式に書
くことができよう。 Hi→=(kxxHx+kxyHy+kxzHz) +i→+(kyxHx +kyyHy+kyzHz)j→+(kzxHx+kzyHy +kzxHz)v→ またマトリツクス形態として、 Hi=K.R.H 上式中Rは地理基準(北、東、垂直)の磁力計
の基準における通過マトリツクスである。 全ての磁気擾乱もしくは外乱が含まれている場
合には、磁力計は擾乱された磁界Hpを次のよう
に測定する。 Hp→=Ho→+K.R.H. (1) 陸上車輌の向首方向を探知する用途に対して
は、二軸磁力計の面内で方程式(1)を解けば充分で
ある。即ち、 H→=R-1K-1(Hp→−Ho→) 角度θおよびが小さい大きさである場合に
は、補償された磁界の水平成分は次のように表わ
すことができる。 Hh′x=1xx(Hxp−Hxo) +1xy(Hyp−Hyo)−θHv Hh′y=1yx(Hxy−Hxo) +1yy(Hyp−Hyo)+Hv 上式中〔1xx 1xy 1yx 1yy〕は「軟鉄」補償マトリ ツクス(行列)であり、HxoおよびHyoは「硬
鉄」擾乱もしくは外乱に起因する磁界の成分であ
り、そして Hvは地磁界の垂直成分である。 補償された(即ち現実のもしくは実際の)磁針
向首方向は次式から求められる。 ψ=arc tanHh′y/Hh′x パラメータHxo,Hyo,1xx,1xy,1yxおよ
び1yyは次の手順によつて決定される。 二軸磁力計によりその軸OXおよびOYによつ
て測定される擾乱された磁界Hpの式は次のよう
に書くことができる。 Hpx=Hxo+kxxHx+kxyHy+kxzHz Hpy=Hyo+kyxHx+kyyHy+kyzHz (2) 車輌が補償時点で水平位にある場合にはHz=
Hv(垂直成分)であり、所与の軌道に対して一定
である。 Hx=Hh cosψおよびHy=Hh sinψ(Hhは地磁
界H→の水平成分である)。 HxoおよびHyoは「硬鉄」に起因する2つの
一定の磁界である。 点M″が水平面上に円C2を描く場合には、点
M′は次の中心を有する楕円C3を描く。 Xo=(Hxo+kxzHz) Yo=(Hyo+kyzHz) この楕円の軸の長さはkxx,kyyに依存しそし
て、kxyおよびKyxに依存する基準軸に対し傾き
Φを有する。 磁力計は周知の仕方で2つの直交探針もしくは
プローブによつて構成されており、これらプロー
ブの軸は搬送担体の軸と一致するかあるいはそれ
に平行にするのが有利である。 このような条件下で、磁力計は次式で表わされ
る2つの擾乱された電圧UxpおよびUypの形態に
あるアナログ情報を発生する。 式中ψは向首方向決定装置で補正後に読み取ら
れた磁針向首方向角度である。 Ux非擾乱電圧=Uh cosψ,但しUh=kHh Uy非擾乱電圧=Uh sinψ,但しUh=kHh Uz=Uo,垂直成分kHvに対応する一定電圧。 補償方法は、水平面上における乗物の360゜にわ
たつて分布されたいろいろな配位に関するいくつ
かの測定結果から次の係数を計算することからな
る。 Xo=Uxo+kxzUz Yo=Uyo+kyzUz 1xx,1xy,1yx,1yy, は次のような関係にある。 Ux=Uxp1xx+Uyp1xy−(Xo+XT+Xv)−θ
U Uy=Uyp1yy+Uxp1yx−(Yo+YT+Yv)+U (4) は、実際のもしくは現実の向首方向ψを得るのに
必要な補正された電圧を表わし、ここでXTおよ
びYTは磁性体の変位に対応する補正係数であり、
そしてXvおよびYvは乗物を一時的に流れる電流
に対応する補正係数である。 これら係数を求めるためには、次のような処理
が必要である。 1) 楕円C3の決定 このためには、車輌の各配位に対して磁力計
により供給される電圧Uxp,Uypは楕円のいろ
いろな点M′の軸OX,OYに沿う座標に対応す
るものであると見なす。即ち、M′1(Ux1p,
Uy1p)、M′2(Ux2p,Uy2p)……M′n(Uxnp,
Uynp) 実際にはこの楕円を定義するのに5ないし8
の点で充分である。 2) 楕円C3の円C1への変換 この楕円C3を軸OXおよびOYを有する同心
円C1に変換するために、これら電圧が受ける
変換は計算機または図表処理により行なわれ
る。ベクトルO′O―→(O′は楕円の中心である)の
変換ならびに軸OXに沿う比R/aおよび軸
OYに沿うR/bの二重アフイン変換(Rは円
C1の半径でありそしてaおよびbは楕円C3
短軸である)により楕円C3上の点M′をこの円
C1上の対応の点Mに移す。半径Rは角度を変
更することなく任意に選択することができる。 これら2つの変換は「硬鉄(ハード・アイア
ン)」および「軟鉄(リフト・アイアン)」作用
の補償に対応する。 3) 楕円C3の点と円C1上の対応の点との間の
対応の決定 楕円C3の任意点M′には理論円C1上の対応の
点Mが、次のように影響を及ぼす。即ち角度
OX,OMが、擾乱もしくは外乱によつて偽化
された向首方向ψ′に対応する測定角OX,
OM′から演繹される乗物の非擾乱磁針向首方
向ψを表わすように影響を受ける。 楕円の1点と円の1点との間の対応が最後
に、中心Oの囲りの円C1の点の回転の選択に
よつて作られる。 本方法では次のような選択が用いられる。変換
O′Oおよび二重アフイン変換後に、回転は恒等元
即ち零角の回転に還元される。このことは軟鉄外
乱係数が対称マトリツクスを構成する、即ち1xy
=1yxであることを意味する。従つてA.B.C.D.の
変換された点は第2図に示されている点A′,B′,
C′,D′であり、O′A,O′B,O′CおよびO′Dはそ
れぞれOA′,OB′,OC′およびOD′に平行である。 測定は例えば1回転を等分された水平位にある
乗物の5つの配位において行なわれる。 このようにして磁力計の測定値は、 (Ux1,Uy1),(Ux2,Uy2)……(Ux5
Uy5)として描かれる。 磁力計の基準Xoyにおける楕円の方程式は次の
ように表わされる。 αUx2+βUy2+γUx.Uy+δUx+εUy=1 楕円を定義するパラメータ(α,β,γ,δ,
ε)=Vは方程式W・V=Zを検証する。ここで Z=1 1 1 1 1 W=UX1 2,Uy1 2,Ux1・Uy1,Ux1,Uy1 Ux2 2,Uy2 2,Ux2・Uy2,Ux2,Uy2 〓 〓 〓 〓 〓 Ux2 5,Uy5 2,Ux5・Uy5,Ux5,Uy5 計算機は行列もしくはマトリツクスWを反転し
て次の積を求める。 V=W-1・z 実際には測定は8つのオリエンテーシヨンもし
くは方位(配位)に沿つて行なうのが有利であ
り、パラメータVは測定ノイズを部分的に除去す
るために慣用のいわゆる最小自乗法を用いて計算
することができる。 次いで第2図に示すように、楕円のパラメータ
a,b,Φ,Xo,Yoが計算される。 a=楕円C3の長軸 b=楕円C3の短軸 Φ=磁力計の軸線に対する楕円の軸の傾き角 Xo,Yo=楕円C3の中心の座標 Xo=γε−2βδ/4αβ−γ2 Yo=γδ−2αε/4αβ−γ2 tg=β−α+k√(α−β)2+r2/γ
【式】 但しα+β>0ならばk=1 α+β<0ならばk=1 楕円C3のパラメータが既知となつたならば該
楕円を円C1に変換するための補正量が求められ
る。この目的で変換 Ux−Xo Uy−Yo ならびに二重アフイン変換が行なわれる。 このアフイン変換を行なうために楕円の次のよ
うな基準式が用いられる。 Uex=cosΦ(Ux−Xo)+sinΦ(Uy−Yo) Uey=sinΦ(Ux−Xo)+cosΦ(Uy−Yo) アフイン変換は次のように表わすことができ
る。 Xe=R/a〔cosΦ(Ux−Xo)+sinΦ(Uy−Yo)〕 Ye=R/b〔−sinΦ(Ux−Xo)+cosΦ(Uy−Yo)〕 また、磁力計の基準については、 X=cosΦXe−sinΦYe Y=sinΦXe−cosΦYe また、 X=(X−Xo)(R/acos2Φ+R/bsin2Φ) +(y−yo)(R/asinΦcosΦ−R/bsinΦcos
Φ) Y=(X−Xo)(R/asinΦcosΦ−R/bsinΦcos
Φ +(Y−Yo)(R/asin2Φ+R/bcos2Φ) 従つて補償パラメータを得るための最終処理は
次のように表わされる。 「硬鉄」の補償−Xo −Yo 上式中Rは円C1の半径であつて、任意の値で
ありそしてRyx=Rxyである。何故ならばマトリ
ツクスもしくは行列が対称であるからである。 次に航空機の事例について検討する。 航空機の強磁性が地磁界に何ら影響を与えない
とした場合には、磁界を表わすベクトルの端M
は、航空機があらゆる可能な姿勢を採つたとして
も、磁力計の軸OX,OY,OZの原点である中心
Oを有する球(第3図)を描くことになろう。実
際には、磁力計はいわゆる永久硬鉄磁界に起因し
て磁力計の軸の原点に対しB1,B2,B3から変換
された中心O′を有する球Sz内で移動する端M″を
有するベクトルの成分を測定する。 磁力計に結合された空間および基準内で、これ
ら永久磁界は次のベクトルによつて表わすことが
できる。 B→=B1i→+B2j→+B3k→ i→,j→,k→はそれぞれ航空機に取り付けられた
三軸磁力計のマークされた直交軸を表わす。 誘起される軟鉄(ソフト・アイアン)磁界Hi
は次のように表わすことができよう。 H→i=(k11Ht1+K12Ht2+k13Ht3)i+(k21Ht1 +k22Ht2+k23Ht3)j→+(k31Ht1+(k31Ht1 +k32Ht2+k33Ht3)k→ また行列式ではHi→:KHt―→a H→taは航空機の中心線における地磁界である。 全ての擾乱もしくは外乱を含めて、磁力計は擾
乱された磁界Hmを次のように測定する。 Hna=Hta+B+K・Hta=B+(I+K)Hta 向首方向を得るには次式 H→ta=(1+K)-1(H→na−B→) を計算し、しかる後に垂直基準による補償された
地磁界の水平成分を計算すれば充分である。 垂直基準および航空機内の地磁界の測定値から
磁針向首方向の計算を可能にする方程式は周知の
ように次のとおりである。 軸(X,Y,Z)の〔A〕、航空機基準および
軸(Nm,Em,V)の〔Gm〕、極部地磁基準。 任意時点tにおいて〔A〕の〔Gm〕への遷移
は連続回転(−,−θ,−Cm)によつて定義さ
れる第4図参照。 なお :横揺れ角 θ:縦揺れ角 Cm:磁針向首方向 地磁界ベクトルは次式によつて定義される。 〔A〕ではHt=(Ht1+Ht2+Ht3) 〔Gm〕ではHt=(Kn,O,Hv) 基準〔Gm〕内のパラメータは次のようにして
モジユールおよび傾斜データから演繹される。 Ho=Htcosα Ht:磁界ベクトルのモジユール Hv=Htsinα α:水平面に対する磁界の傾き 従つて次の関係が成り立つ。 (11) Hn=Ht cosα=(cosθcos Cm)Ht1 +(cos Cm sinθsin −sin Cm cos)Ht2 +(cos Cm sinθcos +sin Cm sin)Ht3 (12) O=sin Cm cosθHt1 +(sin Cm sinθsin +cos Cm cos)Ht2 +(sin Cm sinθcos −cos Cm sin)Ht3 (13) Hv=Ht sinα=−sinθHt1 +sincosθHt2 +coscosθHt3 第2番目の式から磁針向首方向Cmは次のよう
にして計算できる。 Cm=Arc tan(Ht3sin−Ht2cos/Ht1cos
θ+Ht2sinθsin+Ht3sinθcos)+kπ アーク・タンジエント関数に由る180゜の不確定
性は式(11)および(12)から得られるcos Cmの記号を
チエツクすることにより取り払われる。 即ち cos Cm=Ht1cosθ+Ht2sinθsin+Ht3sinθcos/Ht
cosα 上式中αは定義により−(π/2)および+
(π/2)間に含まれるものであり、従つてHt
cosα>0である。 cos Cmは従つて式tan Cmの分母である分子
の記号であり、この理由からCmの不確定性を取
り除くシーケンスは付加的な代数計算を必要とし
ない。 硬鉄および軟鉄補償量を計算するためには、磁
力計によつて発生される測定量だけが用いられ
る。 或るベクトルに対してはそのモジユールからの
情報だけが測定用三面体から独立している。従つ
て地磁界ベクトルのモジユールは測定三面体には
依存せず仮説により無視し得る程に充分弱い勾配
を有する。時間変動は零とみなしてよい。 上に述べたことは次のように定式化することが
できる。 任意のtに対しdHt/dt=0ならば、 t=√t1 2t2 2t3 2 上式中Ht1,Ht2,Ht3は航空機の中心線におけ
る磁界の成分であり、Htは地磁界のモジユール
である。 硬鉄および軟鉄擾乱をモデル化すると(測定ベ
クトルが楕円体上を移動する)、測定成分は次式
で表わされる。 Hni=Hti+Bi3j=1 Kij Htj 測定磁界のモジユールの方程式は、(擾乱もし
くは外乱が測定成分に対し一次の無限小であると
言う仮定の基に)第一次の項だけを保留して次の
ように書くことができる。n 2t 2+23i=1 BiHti+23j=1 Hti3j=1 kijHtj+ε 上式中εは二次および三次の無限小の和であ
る。 この式を時間tに関して微分すると次式が得ら
れる。 dHn 2/dt=dHt 2/dt+23i=1 BidHti/dt +3i=1 3j=1 (kij+kji)d(Hti+Htj)/dt Htが定数であるので次式が得られる。 dHn 2/dt=23i=1 BidHti/dt3i=1 3j=1 (kij+kji)d(Hti+Htj)/dt この方程式は線形であり、9つのパラメータを
有する。この方程式は軟鉄(ソフト・アイアン)
の対称部分だけの推定を可能にする(kij+kji)。 さらにパラメータの係数は地磁界の推定成分に
よつてのみ近似することが可能な地磁界の成分の
導関数である。9つのパラメータは慣用の帰納的
最小自乗法によつて推定される。 各標本化時点において、磁力計からの情報がサ
ンプリングされ、これら新しい測定値はパラメー
タの推定値を改善するために考慮される。 従つて時点nにおいては d(n 2)=2 no2 no-1 かつH→to=〔I+Ko-1-1(H→no−B→o-1) これを解くと、 dH2 n=23 〓 〓i=1 Bid→Htei+ΣΣ(kij+kji)α(Htei・Htej)/d
t 上式中H→teは推定地磁界である。 測定の結果が示すところによれば、〔K〕の項
は、磁力計がよく調和されておつて、搬送体のか
なり大きな擾乱もしくは外乱から相当に離れてい
る場合には非常に弱く、従つて次式が成り立つ。 〔I+K〕-1=I−K さらにマトリツクス(行列)が対称であると言
う仮説は実際上、磁力計の軸が搬送体の軸と一致
している場合に確証される。 第3図は航空機が空間であらゆる方向に移動し
得る時に端M′が楕円体Eを描く磁界ベクトルを
示す。補償方法は、第1番目の段階において、こ
の楕円体をO′上に心出しされた球体S′に変換し、
しかる後に該球体S′を変換O′O―→により磁力計の軸
OX,OYおよびOZの原点O上に心出しすること
からなる。 このような条件下においてはは擾乱がない
とした場合の理論ベクトルH→tである。平面P
(OX,OY)に関しては次のように記述すること
ができる。 角度OX,OH′は擾乱された磁針向首方向ψ′を
表わし、そして 角度OX,OHは磁針向首方向ψを表わす。 第4図は航空機中心線の軸X,YおよびZから
局部地磁基準Nm,Emおよび垂線Vの軸への遷
移を可能にするいろいろな回転を示す。即ち、 Zを中心とする磁針向首方向Cmの回転、 Y″を中心とする縦揺れ角θの回転、 X′を中心とする横揺れ角の回転、 なお角度αは磁界の傾きを表わす。 本発明は、添付図面を参照して以下の説明から
一層容易に理解されるであろう。 図面を参照するに、陸上車輛(第5図)の実際
の向首方向装置は非磁性材料から作られた本体1
0内に次のような要素を有するセンサ・アツセン
ブリを備えている。 電源50。 2つの外乱もしくは擾乱を受けたアナログ電圧
UxpおよびUypを発生する好ましくは2つの軸に
縮小された静的磁力計20。この磁力計としては、
例えば本出願人のフランス国特許願第73−00084
号明細書に開示されている型のものがある。 2つの電解レベル・センサ301,302から
なり、所定の緯度に対して一定である地磁界の垂
直成分を現わす電圧Uoから乗物の傾きを補償す
るための2つの電圧θUoおよびUoを供給する傾
斜計30(第6図)。 この傾斜計30には電子回路が組合わされてお
り、この電子回路は、例えば400Hzの発振器30
5とそれに続いて設けられて各レベル・センサ
(水準センサ)に逆位相の2つの信号を供給する
と共に復調器307および307′に線路306
を介して同期信号を供給する整形段304から構
成されている。それぞれ復調器307および30
7′の上流側に接続されている抵抗器R1およびR2
に現われる信号の振幅は、各センサの傾き角の絶
対値に比例し、そしてその位相は傾き角の記号を
表わす。復調器の出力側に接続された同期復調能
動フイルタ308は、センサの傾き角に直接比例
する振幅および同じ記号を有する出力信号を得る
ことを可能にする。 磁力計および傾斜計からのデータは、補償モジ
ユールと称される電子装置40に供給される。該
電子装置40は補償の時点で計算機80によりそ
のメモリに記憶されている軟鉄(ソフト・アイア
ン)係数1xx,1xy,1yyから補償された電圧Ux
およびUyを計算する。 これら電圧は接続ボツクス70を介してコンパ
ス・カード型の表示装置60に送られる。 さらに乗物に設けられているセンサでスイツチ
33および34が作動されてそれによりモジユー
ル40には、磁性体の明確に定義された変位
(XT,YT)および/または乗物の搭載回路を一時
的に流れる電流(Xv,Yv)に対応し前もつて計
算されている補完用硬鉄(ハード・アイアン)補
正係数が供給される。 補償モジユール40は、向首方向装置が実際の
向首方向を決定するために地上で用いられている
かあるいは補償過程中測定量を得るために用いら
れているかに依存して、2つの仕方で動作する。 乗物が移動しつつある時には、補償モジユール
(第8図)は次の2つの補償された電圧を計算す
る。 Ux=Uxp1xx+Uyp1xy−(Xo+XT+Xv)−θUo Uy=Uyp1yy+Uxp1yx−(Yo+YT+Yv)−Uo (4) 上式中Uhは実際の磁針向首方向読取り量であ
り、 Ux=Uh cosψ Uy=Uh sinψ、但しUh=K Hh 補償過程中、モジユール40は磁力計および傾
斜計の測定量を計算器80に伝送し、該計算機8
0はそれから補償係数を演繹する。補償過程の終
時に、計算機80はこれら係数をモジユール40
のメモリ2に装入する。そこでモジユール40は
電圧UxpおよびUypに補正を行ない、実際の磁針
向首方向ψを求めるに当つて乗物の傾斜を考慮す
ることができる。 電子装置40(第8図)は2つのサブアツセン
ブリを有している。第1のサブアツセンブリは例
えば論理C−Mosでデイジタル・モードで機能す
るものであつて、クロツク4、シーケンス発生器
5およびアドレス計数器を備えたメモリ2を有す
る。 このサブアセンブリは補償の時点で計算機80
から到来する制御信号および情報を受けて、計算
機80を切離した時にはメモリ2によるデータ・
バンクとしての働きをなす。 第2のサブアツセンブリはアナログ・モード
(可変電圧モード)で動作するものであつて、入
力マルチプレクサ1、基準電圧発生手段7、乗算
器として用いられるデイジタル−アナログ変換器
18、出力マルチプレクサ6および出力増幅加算
器3,3′および3″を備えている。 このサブアツセンブリは磁力計20からデータ
UxpおよびUypを受けると共に磁性体の相対変位
または電動機の起動に対応する搬送乗物の磁気構
造形態の変更中にのみ考慮されるXoおよびYoな
らびに他の補完用電圧XT,YTおよびXv,Yvの
ような硬鉄補正電圧を受ける。 向首方向は次のようにして計算される。 クロツク4からのパルスで、シーケンス発生器
5は、例えば電圧Uxの第1の補正項1xxをメモ
リ2に読込むためにアドレス計数器9を増分す
る。同時にマルチプレクサ制御装置15を介して
シーケンス発生器5は磁力計20からの入力項
Uxpならびにマルチプレクサ6の出力チヤンネル
1xxを選択する。乗算「Uxp1xx」がデイジタル
−アナログ変換器18で、レジスタ8を介して実
施されると、その結果が、アナログ電圧の形態
で、マルチプレクサ6の出力端に到るチヤンネル
(1xxUxp)に接続されたコンデンサによつて記
憶される(第7図)。 このようにして各式UxおよびUyの項が1つず
つ計算され、その部分結果はそれぞれ、各計算項
目の各線路、即ち1xx,1xy,Xo,XTおよびXv
ならびに1yy,1yx,Yo,YTおよびYvのマルチ
プレクサ6の出力端に設けられているコンデンサ
によつてアナログ形態で記憶される。 全べての項が計算されたならば、方程式(4)にし
たがつてUxおよびUyを求めるための加算が2つ
の増幅加算器3および3′によつて行なわれる。
なお、これら加算器は、その入力端および出力端
にそれぞれコンデンサCx,Cyならびに抵抗器Rx
およびRyが接続されておつてろ過機能をも行な
うように構成されている(第7図)。 地磁界の垂直成分Hvを表わす電圧Uoはメモリ
2に記憶されている対応の情報もしくはデイジタ
ル量から同じ仕方で計算される。この情報はまた
回路点7で供給される基準電圧Vrefに乗ぜられ
て傾斜計30の入力端に印加され(第8図)、そ
れにより電圧θUoおよびUoが得られ、これら電
圧はチヤンネルUxおよびUyに対応する増幅加算
器3および3′に供給される。硬鉄補正Yo,Xo
および補完補正XT,YTならびにXv,Xvも、補
償の時点でメモリ2に記憶されている硬鉄(ハー
ド・アイアン)パラメータのデイジタル量を乗ぜ
られた基電圧Vrefから同じ仕方で計算される。 補償の過程中、即ち補償計数を決定する時点
で、向首方向装置は補償計算機と称される計算機
80は接続されて、この計算機により動作を制御
される。計算機80は接続ボツクス70を介して
向首方向装置から到来するデータUxp,Uyp,
θUo,UoおよびUoを考慮する。計算機は補償
パラメータを計算してこれらパラメータをデイジ
タル信号の形態でモジユール40のメモリ2に供
給する。そこでチヤンネルa(第8図)によりク
ロツク4を切離してチヤンネルbによりメモリ2
のアドレス計数器9を零にリセツトし、発生器5
を介して入力直列チヤンネルcにより補償データ
をデイジタル形態でレジスタ8に装入し、そこで
メモリ2において全べての補償パラメータの記録
が行なわれてしまうまで、チヤンネルdに現われ
る有効化信号を用いて選択されたアドレスで、メ
モリ2への上記レジスタ8の内容の転送を制御す
る。このメモリ2は、搬送乗物の搭載回路の給電
が遮断した場合でもデタの損失を阻止するため
に、例えばリチウム電池によつて発生される安全
電圧で非揮発性にされている。 向首方向装置の補償モジユール42より発生さ
れる前述の2つの電圧Ux=R cosψおよびUy=
R sinψは、2つの役割を有する表示装置60
(第9図)に伝送される。この表示装置60は、
固定の基準マークの前で回転するコンパス・カー
ドにより乗物が辿る磁針向首方向を永久的に表わ
すことを可能にすると共に、他方またはインデツ
クスにより既に表示されている辿るべき向首方向
と乗物が実際上取る向首方向との間の角度偏差を
表わす電圧Qの形態でアナログ情報を発生する働
きをなす。 この表示手段は、連続的に360゜回転して正弦波
発生機能を行なう2つのポテンシヨメータ91,
92を有している。該ポテンシヨメータの摺動子
91aおよび92aは直流サーボ制御モータ93
の軸100に連結されている。これら摺動子はこ
の軸100上に互に90゜の間隔をあけて固定され
ており、それにより一方のポテンシヨメータは
sinψ′機能を行ない、他方のポテンシヨメータは
cosinψ′機能を行なう(ここでψ′は機械的零基準
に対しポテンシヨメータの摺動子によつて作られ
る角度である)。 モータ93の軸にやはり固定されているダイヤ
ルおよびコンパス・カードは該軸と共に回転す
る。 このサーボ制御の目的はψ′−ψ=0を実行する
ことである。 このためには次の方程式を解かねばならない。 R sin(ψ′−ψ)=R sinψ′cosψ −R sinψcosψ′=0 2つのインバータ増幅器110が反転された電
圧(−R sinψ)および(−R cosψ)を供給
する。 sinψ′機能を行なう第1番目のポテンシヨメー
タ91はその出力端91bに、R sinψ′cosψを
表わす正弦波電圧を発生し、他方cosψ′機能を有
する第2のポテンシヨメータの出力端92bはR
sinψcosψ′の形態にある正弦波信号を発生する。
2のインピーダンス・アダプタ94および95に
よりこれら2つの電圧のレベルの調整後に、これ
ら電圧は高利得の差動増幅器96に印加される。
該増幅器96はその差を求めて得られた電圧をサ
ーボ制御モータ93に供給する。このサーボ制御
モータ93の回転で固定のインデツクス(指標)
97の前でコンパス・カードのダイヤル98が駆
動されて向首方向を表わす。 増幅器96から発生する電圧Tの極性は、アツ
センブリの機械的形態が、ψ′−ψ<0またはψ′−
ψ>0のいずれに対応するかに依存して変わり、
それによりψ′−ψ=0となるまでポテンシヨメー
タと協働してモータを対応の1つの方向に回転せ
しめる。ψ′−ψ=0になる時点で、基準マーク9
7の前に設けられているコンパス・カードによつ
て表わされる角度は修正された向首方向の角度ψ
に対応する。 最後に、表示装置はノブによつて示される向首
方向に対する向首方向偏差に比例する電圧を発生
することができる。インデツクスが設けられてい
るこのノブは、コンパス・カードのダイヤル98
に摩擦結合されておつて、辿るべき向首方向の表
示に依存して摺動子99を位置決めすることを可
能にする(第7図)。 この摺動子は円形のポテンシヨメータ軌道98
a上を摺動し、その場合の変化は360゜にわたり線
形である。この軌道には電圧Qoが印加される。
摺動子99から得られる電圧はしたがつて、向首
方向偏差に対応する角度に比例する。 磁力計および傾斜計の直交軸XおよびYは組立
時に、乗物もしくは車輛の軸に調和される。この
目的でそれぞれのX軸は乗物もしくは車輛の縦軸
線と一致するかあるいは該縦軸線に平行にしなけ
ればならない。 既に述べたように、乗物の補償係数の計算は計
算機80を用いて行なわれる。この計算機80
は、補償過程に必要な時間中だけ、接続ボツクス
70を介して向首方向装置に接続し得るようにす
るのが有利である。この計算機には周辺装置が取
付けられており、その機能は乗物の電源電圧によ
つて確保される。軽量で取扱い易い周辺装置ケー
ス(外箱)を得ることができる。 このケースもしくは段はリンク70を介して向
首方向装置のセンサ電圧を測定して補償係数の計
算を保証する。算出された補償係数は電子装置4
0のメモリ2に装入される。電圧Uoの形態で補
償場所の地磁界の傾斜を導入することが可能であ
る。この情報はケースの蓋に設けられている地図
から見付けられる。 補償手順は次のようにして行なわれる。乗物も
しくは車輛を扁平な水平面上に設定する。この表
面は例えば金属製格納庫を設けることにより磁気
外乱から免がれていなければならない。 乗物の可動磁性体の位置を明確に定義したなら
ば、補償ケースを向首装置に接続する。傾斜計か
ら到来する情報で、乗物の水平位を約±0.5度内
に制御することができる。しかる後に、制御ノブ
により、オペレータはセンサで行なわれる測定値
の自動的標本化を行なう。 この測定値標本化もしくは採取は360゜にわたり
近似的に分布された車輛もしくは乗物の7つの他
の方位について繰返えされる。 ケースは自動的に、補償係数の算出を可能にす
る計算を行なう。計算の終時に、ケースは向首装
置のメモリに対して、乗物に固有の磁気外乱を補
正するのに必要な補償係数を供給する。表示ラン
プで補償が終つたこと、したがつてケースを向首
装置から切離すことができることが表示される。 補償ケース(第10図)は次のような構成を有
し機能を行なう。この補償ケースは読出し専用メ
モリ23に格納されているそのプログラムを実行
するためにクロツク421で同期されるマイクロ
プロセツサ21を有している。 このプログラムは実際上一度作成されると、変
更し得ないものであつて、一連の命令を有してい
る。これら命令は、磁力計から到来する測定量
Uxp,Uypならびに向首方向装置で算出された
Uo,Uo,θUoの採取もしくは標本化を可能に
する。 さらに補償係数の計算を行なうことを可能にす
る。 さらに補償過程の終時に向首方向装置のモジユ
ール40に補償係数を装入することを可能にする
命令である。 計算の展開にまた、向首方向装置の測定標本値
を、変換器29でアナログ−デイジタル変換した
後に一時的に記憶すると共に計算過程における中
間結果を一時的に記憶するための主メモリ24が
必要とされる。ケースはさらにオペレータによつ
て操作される制御モジユール25を有しており、
この制御モジユール25は本質的に次の要素から
なる。即ち、動作モードを選択するための回転接
触器もしくはチヤンネル・セレクタと標本化形態
に変更があつた後にも自動標本計算を可能にする
「連続」制御用プツシユ式接触器(コンタクト)
と、電圧Uoの形態で地磁気の傾斜を導入するた
めのポテンシヨメータである。 マイクロプロセツサ21はアドレス母線212
により、そのプログラム・メモリに格納されてい
る命令の選択を行ない、データ母線213を介し
て読出し専用メモリ23にこれら命令を読込む。 この母線はまた向首方向装置から到来するデー
タ、即ちマルチプレクサ27、フイルタ28、ア
ナログ−デイジタル変換器29およびカプラ30
を介してデイジタル形態で到来するUxp,Uyp,
Uo,θUoおよびUoのデータを上記マイクロプ
ロセツサが受けることを可能にする。 最後に、母線213はまた、プログラムの展開
中に必要とされる中間記憶のために、主メモリ2
4に対するアクセス(読込みおよび読出し)をも
可能にする。 こゝで、「母線」とは複数の回路に対して共通
の導体束を意味する。 ケースはまた補償過程中に異常が生じた場合
に、例えば警報のような光メツセージを例えば3
つのランプを用いて表示することを可能にするイ
ンターフエース126を備えた表示モジユール6
0を有している。横揺れ角および縦揺れ角の値
は、標本採集中、例えば発光ダイオードによつて
構成された7セグメント型の表示モジユールで表
示される。 サービス・モジユール216は、サービス母線
214を介して、いろいろなモジユールおよびメ
モリに対する情報を有効化する働きをなす。カプ
ラ211によりケースもしくは外箱はまたデイジ
タル線路上に、向首方向装置に対し補正パラメー
タを供給するのに必要とされるサービス信号を発
生する。 マイクロプロセツサ21は、マルチプレクサ2
7を介して、補償の計算に介入する第1のパラメ
ータ(例えばUxp)に対応するチヤンネルを選択
する。 このパラメータはアナログ形態で与えられて、
28でろ過され、変換器29によりデイジタル形
態に符号化され、そしてカプラ30を介してデー
タ母線213に供給される。マイクロプロセツサ
21はこの情報を読み取り、記録しそしてアドレ
ス母線217および選択線路215を介してマル
チプレクサ27の他のチヤンネルを選択すること
により次のパラメータの採取に移行する。 5つのパラメータUxp,Uyp,Uo,θUoおよ
びUoが符号化され記録されると、ケースはオ
ペレータに対し、指示ランプを介して、新しい磁
気パターンならびに乗物の新しい方位もしくは配
位を示すように要求し、これら5つのパラメータ
をこの新しいパターンで符号化しかつ記録する。 標本採取が乗物のいろいろな可能な磁気状態に
おいて8つの配位もしくは方位で行なわれたなら
ば、マイクロプロセツサ21は補償係数の計算を
行なう。この計算が終ると結果はデイジタル形態
で、カプラ211を介して向首方向装置の補償モ
ジユール40のメモリ2に供給される。 上に述べた過程中に計算機によつて計算された
係数から、磁力計により供給される電圧を修正す
るためのこの装置は、磁力計の環境に影響を与え
る強磁性から生ずる磁気外乱もしくは擾乱を迅速
にかつ効果的に補償することを可能にするもので
ある。 この装置は特に、乗物に搭載されて、任意時点
で該乗物の実際の磁針向首方向を算出する「向首
方向装置」を構成するのに適している。 最後に本発明の範囲から逸脱することなく、上
に述べた補償方法を、(向首方向誤差式にsin3ψ,
cos3ψ等々を考慮することにより)楕円以外の数
学的に一層展開されたモデルを用いて、擾乱され
た磁力計測定値を補償する用途にも拡張し得るこ
とは言う迄もない。この場合本発明による補償方
法においては、向首方向基準を用いることなく、
360°にわたつて描かれた充分な数の磁力計測定値
から採用される幾何学的モデルの特性であるパラ
メータを同定し、しかる後に既述の例の場合と同
様に該採用されたモデルを表わす幾何学的形態を
基準軸の原点に心出しされた円に変換することが
できるように補償係数を決定することからなる。 第11図に示した航空機の向首方向装置に関し
て述べると、この装置は、3つの磁力計センサが
内部に固定して設けられている非磁性ボツクスに
よつて構成される三軸静的磁力計401を有して
いる。これらセンサの検知軸は、組立時に、航空
機のピツチ(縦揺れ)姿勢、横揺れおよび片揺れ
軸と調和されている三直角三面角を画定する。各
センサは周知の原理に従つて動作し、高い透磁率
を有する棒をAC信号によるヒステレシス・サイ
クルに曝す。 外部磁界は上記ヒステレシス・サイクルに非対
称を惹起し、この非対称は検出されて、コイル内
を循環する電流によつて発生される反対磁界によ
り相殺される。の電流から外部磁界に比例する電
圧が求められ、それにより該外部磁界が測定され
ることになる。 さらに垂直ジヤイロスコープ402もしくは慣
性航行装置が設けられ、この装置は航空機の横揺
れおよび縦揺れ姿勢もしくは角度θを表わすア
ナログ電圧かまたはデイジタル信号を発生する。 さらに計算装置410が設けられ、この計算装
置は実時間で補償係数の計算および実際の磁針向
首方向の計算を行なう。 この計算装置は、例えば、固定小数点で16ビツ
トでの計算かまたは浮動小数点で24ビツトでの計
算のいずれかで用いられるN−MOS型式の8ビ
ツト・マイクロプロセツサMC6800とすることが
できる。 この計算装置410はまたデータ母線414お
よびアドレス母線415に接続されているメモリ
を備えている。これらメモリには、例えば8Kオ
クテツトのプログラム・メモリ411(N−
MOS型式のREPROM)、512オクテツトの第1
の主メモリ412(例えばC−MOS型のRAM)
および、航空機の電源が1つの飛行から他の飛行
に切り換わる時に遮断される際に自動補償結果を
安全保護するためのデータを格納しておくのに用
いられる1Kビツトの第2の主メモリ413を含
む。後者のメモリはC−MOS型のRAMメモリに
よつて構成される。 入/出力インターフエースは到来するデータを
取り込み、デイジタル的に符号化することを可能
にする。これらインターフエースには、マルチプ
レクサ403および404、ジヤイロスコープ4
02のためのピーク値標本化回路405、磁力計
401から発生される3つの電圧を逐次受け取る
標本化保持回路406が含まれる。 標本化回路は制御パルスを発生するクロノメー
タ407によつて制御される。 このようにして取り込まれた入力信号は、マル
チプレクサ408に供給されかつアナログ−デイ
ジタル変換器409に供給され、しかる後に計算
装置410で利用することができる。この計算装
置410はデイジタル形態で、航行に必要ないろ
いろなデータ即ちピツチ横揺れおよび自動補償さ
れた磁針向首方向を算出する。 実施に当つては、ジヤイロスコープおよび磁力
計の軸が航空機の対応の軸と調和するように、磁
気擾乱もしくは外乱補償装置を配列するのが有利
であろう。 補償係数は、計算機が搭載されている場合には
実時間で計算することができるし、また計算機が
地上に接地されている場合には飛行中に取り込ま
れた情報を記憶しておいて補償係数を計算するこ
とができよう。
【図面の簡単な説明】
第1図は第5図の乗物の基準軸を示し、第2図
は第5図の装置の磁力計から供給される電圧を図
解し、第3図は第11図の装置の磁力計から供給
される電圧を図解し、第4図は航空機の中心線な
らびに第11図に示した装置に設けられる局部地
磁基準を示し、第5図は陸上車輌の向首方向装置
をブロツクダイヤグラムで略示し、第6図は第5
図の装置の傾斜計およびそれに関連の回路を示す
ブロツクダイヤグラム、第7図は第5図の装置の
電子系統に含まれる出力マルチプレクサの出力回
路を示し、第8図は第5図の装置の電子系統の詳
細図であり、第9図は第5図の装置の表示回路を
示し、第10図は第5図の装置の補償ケース(外
箱)のブロツクダイヤグラムを示し、そして第1
1図は航空機の向首方向装置を略示するブロツク
ダイヤグラムである。 1……入力マルチプレクサ、2……メモリ、3
……増幅加算器、4……クロツク、5……シーケ
ンス発生器、6……出力マルチプレクサ、7……
電圧発生手段、8……レジスタ、10……本体、
20……磁力計、21……マイクロプロセツサ、
23……読出し専用メモリ、24……主メモリ、
25……制御モジユール、27,408……マル
チプレクサ、28,308……フイルタ、29,
409……アナログ−デイジタル変換器、30,
211……傾斜計、33,34……スイツチ、4
0……電子装置、50……電源、60……表示装
置、70……接続ボツクス、80……計算機、2
12……アドレス母線、213,414……デー
タ母線、216……サービス・モジユール、30
1,302……レベル・センサ、304……整形
段、305……発振器、306……線路、307
……復調器、401……磁力計、402……ジヤ
イロスコープ、405……標本化回路、406…
…標本化保持回路、407……クロノメータ、4
10……計算装置、411……プログラム・メモ
リ、412,413……主メモリ。

Claims (1)

  1. 【特許請求の範囲】 1 陸上車輌に搭載されている磁針向首方向を決
    定する装置の測定に影響を及ぼす磁気外乱を補償
    するための装置において、擾乱された地磁界の基
    準軸線に沿う成分に比例するアナログ電圧を供給
    する磁力計20と、いわゆる補償係数を該磁力計
    と組み合わせて装入されている計算機モジユール
    40のメモリとを備えていて、該磁力計により供
    給される情報を補正して実際の向首方向を求める
    もので更に装置は乗物が水平面上での360゜にわた
    る車輌の自転に際して、擾乱された磁界を表わす
    第1のベクトルHhp―→の端M′によつて描かれる閉
    じた平面の幾何学的形状C3を同定する手段と、
    該閉じた幾何学的形状を、前記磁力計の基準軸の
    原点に心出しされて、擾乱されていない磁界を表
    わす第2のベクトルHh→の端Mの軌道を表わす半
    径Rの円C1に変換する手段と、そして前記閉じ
    た幾何学的形状C3の全ての点と前記円C1上の該
    点のホモロジーとの間の対応を定める係数を計算
    する手段とを備える磁針向首方向測定装置に影響
    する磁気外乱を補償する装置。 2 前記閉じた幾何学的形状が楕円C3である特
    許請求の範囲第1項記載の磁針向首方向測定装置
    に影響する磁気外乱を補償する装置。 3 前記車輌が補償過程中水平面上にあることを
    チエツクするため該磁力計と組み合わせて設けら
    れている傾斜計30と、そして擾乱されていない
    地磁界H→の成分(Hx→およびHy→)を表わす2つの
    電圧UxおよびUyを計算するのに必要な補償係数
    を装入する向首方向決定装置のメモリとを備える
    特許請求の範囲第1項記載の磁針向首方向測定装
    置に影響する磁気外乱を補償する装置。 4 前記磁力計から供給される擾乱された電圧
    Uxp,Uypを含む情報をアナログ形態で供給され
    る手段と、前記地磁界の垂直成分を表わす一定電
    圧Uoおよび前記傾斜計から供給される車輌の勾
    配を表わす電圧Uo,θUoと、前記情報を多重化
    する手段27と、前記情報をろ過する手段28
    と、前記情報をデイジタル符号化する手段29
    と、補償係数を計算する計算機80のデータ母線
    213に前記情報を与える手段と、該補償係数を
    伝送する向首方向決定装置の計算機モジユール4
    0のメモリ2とを備える特許請求の範囲第3項記
    載の磁針向首方向測定装置に影響する磁気外乱を
    補償する装置。 5 磁力計の軸に沿う擾乱された地磁界の成分に
    比例するアナログ電圧を発生する磁力計20と、
    電源50と、前記磁力計によつて行なわれる測定
    ならびにメモリに予め記憶されている補償係数か
    ら実際の向首方向を計算する該計算機モジユール
    40とを備えるもので更に装置は磁力計の軸と同
    じ軸に沿う乗物の傾斜に感応して対応する電圧を
    発生するセンサを含む傾斜計30を備えていて、
    計算機モジユール40が磁方計および傾斜計の軸
    に沿う擾乱されていない地磁界の成分に応ずる電
    圧を計算し乗物の横揺れ及び縦揺れ角の影響を考
    慮する特許請求の範囲第3項或は第4項のいずれ
    か記載の磁針向首方向測定装置に影響する磁気外
    乱を補償する装置。 6 前記磁力計20および傾斜計30の軸が搬送
    乗物の軸と調和している特許請求の範囲第5項記
    載の磁針向首方向測定装置に影響する磁気外乱を
    補償する装置。 7 擾乱されていない地磁界の成分に対応する電
    圧が、固定の基準マークの前で回転するコンパ
    ス・カード60を有し、乗物の磁針向首方向を与
    える表示手段によつて利用される特許請求の範囲
    第5項記載の磁針向首方向測定装置に影響する磁
    気外乱を補償する装置。 8 航空機に搭載されている磁針向首方向を決定
    する装置の測定に影響を及ぼす磁気外乱を補償す
    るための装置において擾乱された地磁界の基準軸
    線に沿う成分に比例するアナログ電圧を供給する
    磁力計と、いわゆる補償係数を該磁力計と組み合
    わせて装入されている計算機モジユールのメモリ
    とを備えていて、該磁力計により供給される情報
    を補正して実際の向首方向を求めるもので、更に
    装置は地磁界のモジユールが2つの相続く測定間
    において一定であると仮定し、擾乱された地磁界
    のモジユールの空間における展開を検査する手段
    と、外乱を受けた測定と実際の測定との間の対応
    を定めることにより、前記磁力計の軸に心出しさ
    れた理論球体に楕円体から地磁界のモジユールの
    展開を戻すことを可能にする補償係数を反復的な
    仕方で計算する手段とを備える磁針向首方向測定
    装置に影響する磁気外乱を補償する装置。 9 前記磁力計により測定され、該磁力計の3つ
    の軸OX,OY,OZに沿う磁界成分を表わす3つ
    のアナログ電圧を、デイジタル符号化する手段
    と、しかる後に実時間で補償係数を計算する計算
    機のデータ母線とを備える特許請求の範囲第8項
    記載の磁針向首方向測定装置に影響する磁気外乱
    を補償する装置。 10 前記計算機はまた、垂直ジヤイロスコープ
    または慣性測定装置から供給される補完横揺れお
    よび縦揺れ姿勢データから実際の磁針向首方向を
    計算するように配される特許請求の範囲第9項記
    載の磁針向首方向測定装置に影響する磁気外乱を
    補償する装置。 11 前記表示手段が360゜回転する2つのポテン
    シヨメータ91,92を備え、該ポテンシヨメー
    タの摺動子は90゜偏位しておつて、前記コンパ
    ス・カード60にも接続されているサーボ制御モ
    ータ93の軸に連結されており、さらに前記ポテ
    ンシヨメータ91,92によつて供給される電圧
    を受け、前記サーボ制御モータ93に接続されて
    いる出力端を有する差動増幅器96を備えている
    特許請求の範囲第10項記載の磁針向首方向測定
    装置に影響する磁気外乱を補償する装置。 12 磁力計の軸OXおよびOYに沿う擾乱され
    た地磁界Hp→の成分(HxpおよびHyp)に比例す
    るアナログ電圧UxpおよびUypを発生する磁力計
    と、電源50と、前記磁力計によつて行なわれる
    測定ならびにメモリ2に予め記憶されている補償
    係数から実際の向首方向を計算する該計算機モジ
    ユール40とを備えるもので更に装置は3つの軸
    に沿つて地磁界の成分に比例する3つの電気量を
    うるための手段401と、搬送体の横揺れおよび
    縦揺れ姿勢を表わす電気量の形態で垂直基準を供
    給する手段402と、該他の手段から供給される
    電気量から補償係数および実際の磁針向首方向を
    計算するための手段410とを備えた特許請求の
    範囲第8項或は第9項のいづれか記載の磁針向首
    方向測定装置に影響する磁気外乱を補償する装
    置。 13 地磁界の前記成分に比例する電気量を得る
    ための手段が磁力計401であり、そして前記垂
    直基準を与える手段が垂直ジヤイロスコープ40
    2または慣性測定装置を有している特許請求の範
    囲第12項記載の磁針向首方向測定装置に影響す
    る磁気外乱を補償する装置。 14 磁力計が、その軸OX,OY,OZに沿う擾
    乱された地磁界の成分に比例する3つのアナロゲ
    電圧Hm1,Hm2,Hm3を発生する3軸の静的磁
    力計401であり、前記垂直ジヤイロスコープ4
    02は乗物の横揺れ姿勢および縦揺れ姿勢θを
    表わす信号を発生し、さらに2つのマルチプレク
    サ403,404と、第1のピーク値標本化手段
    405と、第2の標本化保時回路406と、制御
    パルスを発生するクロノメータ407と、アナロ
    グ−デイジタル変換器409にアナログ情報を供
    給するための第3のマルチプレクサ408と、擾
    乱されていない磁界の成分を求めて該成分から実
    際の磁針向首方向を導出するのに必要な補償係数
    を求めるための、メモリ411,412,413
    ならびにデータおよびアドレス母線414,41
    5を備えた計算機装置410とを有している特許
    請求の範囲第13項記載の磁針向首方向測定装置
    に影響する磁気外乱を補償する装置。 15 前記磁力計およびジヤイロスコープの軸が
    航空機の軸と調和している特許請求の範囲第13
    項記載の磁針向首方向測定装置に影響する磁気外
    乱を補償する装置。
JP56086729A 1980-06-05 1981-06-05 Method of compensating for magnetic disturbance affecting magnetic needle head direction gauge and apparatus for determining magnetic needle head direction Granted JPS5752813A (en)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
FR8012725A FR2484079A1 (fr) 1980-06-05 1980-06-05 Procede de compensation des perturbations magnetiques dans la determination d'un cap magnetique, et dispositif pour la mise en oeuvre de ce procede

Publications (2)

Publication Number Publication Date
JPS5752813A JPS5752813A (en) 1982-03-29
JPH0340323B2 true JPH0340323B2 (ja) 1991-06-18

Family

ID=9242845

Family Applications (1)

Application Number Title Priority Date Filing Date
JP56086729A Granted JPS5752813A (en) 1980-06-05 1981-06-05 Method of compensating for magnetic disturbance affecting magnetic needle head direction gauge and apparatus for determining magnetic needle head direction

Country Status (8)

Country Link
US (1) US4414753A (ja)
EP (1) EP0041892B1 (ja)
JP (1) JPS5752813A (ja)
AT (1) ATE6443T1 (ja)
AU (1) AU553800B2 (ja)
CA (1) CA1194119A (ja)
DE (1) DE3162423D1 (ja)
FR (1) FR2484079A1 (ja)

Families Citing this family (76)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4672565A (en) * 1981-03-10 1987-06-09 Nippon Soken, Inc. Direction detecting system for vehicles
FR2507770A1 (fr) * 1981-06-11 1982-12-17 Crouzet Sa Procede de compensations des perturbations magnetiques dans la determination d'un cap magnetique et dispositif pour la mise en oeuvre de ce procede
JPS59155714A (ja) * 1982-10-12 1984-09-04 ロケ マナ リサーチ リミテッド 輸送機関用電子コンパス
US4546551A (en) * 1983-03-24 1985-10-15 Prince Corporation Electrical control system
JPS601510A (ja) * 1983-06-17 1985-01-07 Honda Motor Co Ltd 車載地磁気センサの補正方法
US4593359A (en) * 1983-06-24 1986-06-03 Ilan Sadeh Vehicle navigation system
US4797841A (en) * 1983-11-28 1989-01-10 Magnavox Government And Industrial Electronics Company Method and apparatus for automatic calibration of magnetic compass
US4611293A (en) * 1983-11-28 1986-09-09 Magnavox Government And Industrial Electronics Company Method and apparatus for automatic calibration of magnetic compass
DE3345818A1 (de) * 1983-12-17 1985-06-27 Robert Bosch Gmbh, 7000 Stuttgart Verfahren zur ermittlung von stoerfeldern in fahrzeugen mit einem elektronischen kompass
DE3405933A1 (de) * 1984-02-18 1985-08-22 Teldix Gmbh, 6900 Heidelberg Verfahren zur ermittlung des kurses eines bewegbaren objekts
DE3422491A1 (de) * 1984-06-16 1985-12-19 Robert Bosch Gmbh, 7000 Stuttgart Verfahren zur ermittlung der fahrtrichtung eines fahrzeuges mit elektronischem kompass
US4677381A (en) * 1984-10-19 1987-06-30 Prince Corporation Flux-gate sensor electrical drive method and circuit
GB8504949D0 (en) * 1985-02-26 1985-03-27 Shell Int Research Determining azimuth of borehole
DE3509548A1 (de) * 1985-03-16 1986-09-18 Robert Bosch Gmbh, 7000 Stuttgart Verfahren zur ermittlung der fahrtrichtung eines fahrzeuges mit elektronischem kompass
EP0214817B1 (en) * 1985-09-03 1991-12-04 British Aerospace Public Limited Company Calibration of magnetic compasses
DE3534480A1 (de) * 1985-09-27 1987-04-02 Bosch Gmbh Robert Verfahren zur bestimmung des einbauortes eines magnetfeldsensors in kraftfahrzeugen
US4698912A (en) * 1985-12-11 1987-10-13 The Laitram Corporation Magnetic compass calibration
DE3568374D1 (en) * 1985-12-20 1989-03-30 Litef Gmbh Method to determine heading by using and automatically calibrating a 3-axis magnetometer rigidly fitted in an aircraft
JPH0629732B2 (ja) * 1986-10-08 1994-04-20 三菱電機株式会社 移動体用方位検出装置
JPS6394108A (ja) * 1986-10-08 1988-04-25 Mitsubishi Electric Corp 移動体用方位検出装置
DE3644682A1 (de) * 1986-12-30 1988-07-14 Bosch Gmbh Robert Navigationsverfahren fuer fahrzeuge mit elektronischem kompass
DE3644681A1 (de) * 1986-12-30 1988-07-14 Bosch Gmbh Robert Navigationsverfahren fuer fahrzeuge mit elektronischem kompass
JP2642649B2 (ja) * 1987-12-28 1997-08-20 本田技研工業株式会社 車体の進行方位検出装置
DE3833798A1 (de) * 1988-10-05 1990-04-12 Bosch Gmbh Robert Verfahren zur ermittlung des stoereinflusses von verbrauchern auf navigationssysteme
JPH07101174B2 (ja) * 1989-01-06 1995-11-01 日産自動車株式会社 車両用方位計
US5046031A (en) * 1989-02-06 1991-09-03 Magnavox Government And Industrial Electronics Company Method and apparatus for automatic flux-gate compass calibration
FR2643709B1 (fr) * 1989-02-27 1993-11-19 Fabrication Instruments Mesure Procede d'aide a la navigation et navigateur terrestre perfectionne
IL92239A (en) * 1989-11-07 1993-04-04 Israel State Magnetic compass
DE59204446D1 (de) * 1991-09-17 1996-01-04 Siemens Ag Verfahren zur kompensation eines magnetischen störfeldes in einem fahrzeug.
US5828984A (en) * 1991-12-27 1998-10-27 Chrysler Corporation Data processing method for an electronic compass system
US5297065A (en) * 1991-12-27 1994-03-22 Chrysler Corporation Magnetic transient detection and calibration technique for an auto-calibrating compass
CA2134191C (en) * 1993-11-17 2002-12-24 Andrew Goodwin Brooks Method of correcting for axial and transverse error components in magnetometer readings during wellbore survey operations
US5452518A (en) * 1993-11-19 1995-09-26 Baker Hughes Incorporated Method of correcting for axial error components in magnetometer readings during wellbore survey operations
US5517430A (en) * 1994-06-20 1996-05-14 Directional Robotics Research, Inc. Real-time digital orientation device
US5477470A (en) * 1994-06-20 1995-12-19 Lewis; W. Stan Real-time digital orientation device
US5511319A (en) * 1994-08-29 1996-04-30 Prince Corporation Vehicle compass correction circuit
US5737226A (en) * 1995-06-05 1998-04-07 Prince Corporation Vehicle compass system with automatic calibration
US5878370A (en) * 1995-12-01 1999-03-02 Prince Corporation Vehicle compass system with variable resolution
US20030158699A1 (en) * 1998-12-09 2003-08-21 Christopher P. Townsend Orientation sensor
US6356851B1 (en) 1999-02-05 2002-03-12 Delphi Technologies, Inc. Accelerated calibration for electronic compass module
US6301794B1 (en) * 1999-05-27 2001-10-16 Johnson Controls, Inc. Vehicle compass system with continuous automatic calibration
US6577976B1 (en) 1999-09-17 2003-06-10 Hrl Laboratories, Llc Method for dynamic autocalibration of a multi-sensor tracking system and apparatus incorporating it therein
US6536123B2 (en) * 2000-10-16 2003-03-25 Sensation, Inc. Three-axis magnetic sensor, an omnidirectional magnetic sensor and an azimuth measuring method using the same
US6539639B2 (en) 2000-12-06 2003-04-01 Honeywell International Inc. Monitoring accuracy of an electronic compass
US6543146B2 (en) 2000-12-06 2003-04-08 Honeywell International, Inc. Electronic compass and compensation of large magnetic errors for operation over all orientations
KR100533106B1 (ko) * 2002-08-06 2005-12-05 삼성전자주식회사 지자계 센서의 자세 오차 보상장치 및 방법
KR100491597B1 (ko) * 2002-11-28 2005-05-27 삼성전자주식회사 지자기 센서의 보정완료를 자동으로 추출하는 장치 및방법
US20050081334A1 (en) * 2003-10-15 2005-04-21 Tai Wen H. Hinge structure
WO2006035505A1 (ja) * 2004-09-29 2006-04-06 C & N Inc 磁気センサの制御方法、制御装置、および携帯端末装置
US7437257B2 (en) * 2004-10-07 2008-10-14 Yamaha Corporation Geomagnetic sensor and geomagnetic sensor correction method, temperature sensor and temperature sensor correction method, geomagnetism detection device
CN101726291B (zh) * 2004-10-07 2012-02-15 雅马哈株式会社 地磁检测装置
KR100799536B1 (ko) * 2005-11-08 2008-01-31 한국전자통신연구원 2축 지자계 센서의 경사각 오차를 보상하기 위한 가상축 지자계 데이터 추정 장치 및 그 방법과, 그를 이용한 방위각 산출 시스템
KR20080026395A (ko) * 2006-09-20 2008-03-25 삼성전자주식회사 자기 환경을 고려한 컴퍼스 센서의 교정 방법 및 장치와이를 이용한 방위각 측정 방법 및 장치
FR2931248B1 (fr) * 2008-05-14 2012-07-27 Movea Sa Systeme de mesure de champ magnetique comprenant un capteur triaxial de mesure de champ magnetique mobile conjointement avec un element porteur perturbant les mesures et procede associe
US8370097B2 (en) * 2009-06-05 2013-02-05 Apple Inc. Calibration techniques for an electronic compass in portable device
US8108171B2 (en) * 2009-09-14 2012-01-31 Honeywell International, Inc. Systems and methods for calibration of gyroscopes and a magnetic compass
US8825426B2 (en) 2010-04-09 2014-09-02 CSR Technology Holdings Inc. Method and apparatus for calibrating a magnetic sensor
JP5656775B2 (ja) * 2011-08-31 2015-01-21 日立オートモティブシステムズ株式会社 電力変換装置
US9677889B2 (en) 2012-06-21 2017-06-13 Innovative Solutions & Support, Inc. Method and system for compensating for soft iron magnetic disturbances in multiple heading reference systems
US9207079B2 (en) 2012-06-21 2015-12-08 Innovative Solutions & Support, Inc. Method and system for compensating for soft iron magnetic disturbances in a heading reference system
CN104237822B (zh) * 2013-06-20 2018-10-19 意法半导体(中国)投资有限公司 用于电子磁力计传感器的补偿磁干扰
JP6983565B2 (ja) 2016-07-29 2021-12-17 イノベイティブ・ソリューションズ・アンド・サポート・インコーポレイテッド 乗物の機首方位基準システムにおける軟鉄磁気擾乱を補償するための方法とシステム
US10209078B2 (en) * 2016-08-31 2019-02-19 Yost Labs Inc. Local perturbation rejection using time shifting
CN108802839B (zh) * 2018-06-08 2024-06-21 北京桔灯地球物理勘探股份有限公司 基于固定翼无人机的铯光泵磁测方法
FR3082611B1 (fr) * 2018-06-13 2020-10-16 Sysnav Procede de calibration de magnetometres equipant un objet
US10914795B2 (en) * 2019-02-20 2021-02-09 Crocus Technology Inc. Apparatus and method for magnetic sensor output compensation based upon ambient temperature
CN112630807A (zh) * 2020-11-24 2021-04-09 洛阳师范学院 一种地磁效应及卫星导航失灵时的自修正导航方法
CN112630808A (zh) * 2020-11-26 2021-04-09 洛阳师范学院 一种用于地磁效应及卫星导航失灵时的修正导航系统
CN112577477B (zh) * 2020-12-26 2025-02-28 南方科技大学 一种海底柱状取样技术的磁补偿测试装置及其磁补偿方法
CN112858961B (zh) * 2021-03-01 2023-02-07 哈尔滨工业大学 一种针对航空平台多源磁干扰的补偿方法
CN113447993B (zh) * 2021-07-26 2023-09-15 中国人民解放军61540部队 磁力矢量测量的补偿飞行方法、系统及磁补偿方法、系统
CN114184988B (zh) * 2021-11-11 2022-10-11 北京大学 包含补偿平台电流磁干扰的航磁补偿方法及其装置
CN116203640A (zh) * 2023-03-03 2023-06-02 吉林大学 一种基于地磁总场垂向差值的磁性目标体定位方法
CN117537792B (zh) * 2024-01-03 2024-04-30 西南应用磁学研究所(中国电子科技集团公司第九研究所) 电子罗盘自适应方位角矫正方法
CN118426063B (zh) * 2024-04-10 2025-07-18 北京自动化控制设备研究所 抑制地质磁干扰的航空磁补偿方法及系统
CN120949131B (zh) * 2025-10-14 2025-12-12 南方科技大学 一种飞行器航磁补偿方法、系统、设备和介质

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR1085870A (fr) * 1953-05-30 1955-02-08 Perfectionnements apportés aux appareils asservis au champ magnétique terrestre et munis d'un système de correction de quadrantale
GB1202533A (en) * 1966-11-19 1970-08-19 Elliott Brothers London Ltd Improvements relating to magnetically sensitive arrangements
US3621584A (en) * 1969-12-10 1971-11-23 Wayne E Burt Magnetometer compass
US3683668A (en) * 1971-01-26 1972-08-15 Sperry Rand Corp Compass calibrator
US4006631A (en) * 1974-12-11 1977-02-08 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Magnetic heading reference
US3943763A (en) * 1974-12-11 1976-03-16 The United States Of America As Represented By The United States National Aeronautics And Space Administration Office Of General Counsel-Code Gp Magnetic heading reference
NO772126L (no) * 1976-06-17 1977-12-20 Laitram Corp Fremgangsm}te og apparat for frembringelse av deviasjonskorreksjons-signaler for et magnetisk kompass
US4031630A (en) * 1976-06-17 1977-06-28 The Laitram Corporation Calibration apparatus for automatic magnetic compass correction
US4091543A (en) * 1976-06-17 1978-05-30 The Laitram Corporation Automatic magnetic compass correction
IT1073691B (it) * 1976-08-26 1985-04-17 Importazioni Esportazioni Rapp Perfezionamento nelle bussole statiche

Also Published As

Publication number Publication date
AU7140181A (en) 1981-12-10
ATE6443T1 (de) 1984-03-15
DE3162423D1 (en) 1984-04-05
US4414753A (en) 1983-11-15
EP0041892A1 (fr) 1981-12-16
EP0041892B1 (fr) 1984-02-29
CA1194119A (en) 1985-09-24
FR2484079B1 (ja) 1982-06-11
AU553800B2 (en) 1986-07-31
FR2484079A1 (fr) 1981-12-11
JPS5752813A (en) 1982-03-29

Similar Documents

Publication Publication Date Title
JPH0340323B2 (ja)
JPS59155714A (ja) 輸送機関用電子コンパス
Caruso Applications of magnetic sensors for low cost compass systems
CA2431722C (en) Electronic compass and compensation of large magnetic errors
CN114674301B (zh) 一种强磁干扰环境下电子罗盘主动补偿方法及系统
EP0225493B1 (en) Magnetic compass calibration
CN108507553A (zh) 电子罗盘的校正方法
CN115524762B (zh) 基于三维亥姆赫兹线圈的地磁矢量测量系统补偿方法
US3683668A (en) Compass calibrator
CN104316037B (zh) 一种电子罗盘的校正方法及装置
Caruso Applications of magnetoresistive sensors in navigation systems
EP0436794A1 (en) Magnetic compass
US5321631A (en) Method and self-contained system for harmonizing equipments on board a vehicle, using means of measuring the earth's gravitational and magnetic fields
KR101210394B1 (ko) 지자기 검지장치
Liu et al. Compensation of geomagnetic vector measurement system with differential magnetic field method
GB2128749A (en) Electronic compass with tilt compensation
JP4026937B2 (ja) 電子磁気コンパス
JPS6345043B2 (ja)
Fu et al. Research on control algorithm of strong magnetic interference compensation for MEMS electronic compass
US2887873A (en) Method of measuring and compensating for deviation errors in earth's field sensitivedirection indicators
CN108710145A (zh) 一种无人机定位系统及方法
CN116400424B (zh) 基于最小二乘两步法的地磁矢量测量多源误差校正方法
JP2007163388A (ja) 方位センサおよび記録媒体
JPH08278137A (ja) 方位出力装置
NO872526L (no) Innretning til styring av et magnetisk egenbeskyttelses(meb-)-anlegg.