JPH11505620A - 地磁場の方向の算定のための方法 - Google Patents

Abstract

(57)【要約】 装置に固着した磁気的な材料及び電流によって乱されている可能性のある地磁場の方向を、三つの磁場センサーと二つの傾斜感知装置とをもっている電子工学的な磁気コンパスを用いて測定するための方法が提供される。当該電子工学的な磁気コンパスは、Ｎ個の異なる空間的な位置に配置される。その際、これらのＮ個の位置のそれぞれにおいて、傾斜センサー信号及び磁場センサー信号が測定され、これらの信号から傾斜値及び磁場値が算定される。これらの傾斜値及び磁場値に基づいて、地磁場ベクトルの量がベクトル方程式 を用いて算定される。ここで、 である。また、Ｎが当該ベクトル方程式の決定すべきパラメータの数に少なくとも等しくなければならない。

Description

【発明の詳細な説明】 地磁場の方向の算定のための方法 本発明は、装置に固着した磁気的な材料及び電流の磁場により乱されている可 能性のある地磁場の方向を、三つの磁場センサーと二つの傾斜感知装置とをもっ ている電子工学的な磁気コンパスを用いて算定するための方法に関する。 米国特許第４６８６７７２号明細書により、例えば戦車(Ｐanzer)の進行方向 を算定するといわれる、電子工学的な磁気コンパスが知られて 鉛直の軸線のまわりで回転可能である。磁気コンパスは、車両ボディーに配置さ れており、振り子のない三軸の磁力計(pendelfreies dreiachsiges Ｍagnetomet er)を備える。当該磁力計は、当該磁力計の位置における三つの磁場成分を示す 電気的な磁場信号を出力する。二つの傾斜センサーが設けられており、当該傾斜 センサーが、当該傾斜センサーが取り付けられている車両ボディーのピッチ角( Ｎickwinkel)及びロール角(Ｒollwinkel)を示す電気的な信号を出力する。さら に、角度測定装置が二つの鉄製ボディーの間の角度についての回転角度信号を与 える。記憶装置には、複数の予め較正された修正ファクター(Ｋorrekturfaktore n)が記憶されている。これらの修正ファクターが、両方の鉄製ボディーの間の複 数の回転角についての磁力計の測定への車両によって誘導される磁場の影響を修 正するとされている。コンピューターを使って、リ アルタイムで戦車の進行方向が、電気的な磁場信号と、対応する回転角について の記憶された修正ファクターによって修正された傾斜センサーの信号とから算出 される。記載された磁気コンパスが使用できるようになる前には、較正測定を行 うことが必要不可欠である。このために、車両が少なくとも二つの異なる互いに 平行でない平面上に配置される。これらの位置で、車両ボディーの異なる方向づ けでしかし装甲砲塔の等しい相対的な方向づけの場合に、方位角(Ａzimutwinkel )、ピッチ角、及びロール角の測定が行われる。これらの角度を測定できるよう に、例えば経緯儀(Ｔheodolit)が必要不可欠である。較正の際には、測定箇所に おける地磁場が知られていることが前提とされる。このために、そのときどきの 測定箇所についての地磁場の磁気的な傾角及び偏角(magnetische Ｉnklination und Ｄeklination)の値がカード資料から読み取られる。 米国特許第４５３９７６０号明細書には、三つの磁気センサー(Ｍagnetsensor en)を有する、車両のための電子工学的な磁気コンパスが記載されている。これ らの磁気センサーは、地磁場及び車両と結びつい を示す。当該磁気センサーは、これらの成分に対応する電気的な信号を発生させ る。さらに、傾斜センサーが水平な平面に関しての車両の傾斜に反応を示す。車 両が磁気コンパスの較正のために一つの円上を通って回転させられると、データ 処理装置と記憶装置とを用いて、較正修正値(Ｋalibrierkorrekturwerte)として 測定センサーから導出される信号が記憶される。磁気的なノイズ場の影響を排除 するために、データ処理装置が較正過程の終了後に較正修正値を考慮にいれて車 両の位置で の地磁場についての修正された値を算出する。それから、これらの修正された値 を出発点として、水平線(Ｈorizont)の算定に用いる傾斜センサーの値を使用し て、車両進行方向の方位が算出される。これらの算出の際には、修正行列(Ｋorr ekturmatrix)が対称であることを出発点とする。このことが当てはまるケースは 、そもそも、まれにしかない。 本発明の課題は、比較的少ない経費で(weniger aufwendig)実行され得る、電 子工学的な磁気コンパスを用いて地磁場の方向を算定するための方法を提供する ことである。 本発明により、この課題は、請求項１の特徴的構成によって解決される。 本発明の対象物の有利な別の構成は、従属項に記載されている。 本発明に係る方法は付加的な測定装置による特殊な較正測定を必要としないこ とが有利である。当該電子工学的な磁気コンパスの場合には、使用するためにそ れを準備するための、特別なデータの入力も行われない。当該磁気コンパスをさ まざまな任意の空間位置に配置することだけが必要である。それぞれの空間位置 において、三つの磁場成分が算定されることが有利である。望まれるならば、こ れらの磁場成分が互いに直交しているとよいが、しかしそうでなければならない わけではない。それぞれの空間位置で付加的に測定された傾斜量、すなわちピッ チ角及びロール角を用いて、そのときどきの磁場成分から地磁場の実際の磁場ベ クトルの方向が算定され得る。 本発明に係る方法は、電子工学的な磁気コンパスがその取り付けのまえに較正 されていなければならないということを要求しない。本発明に係る方法によって 、磁気的なノイズ場だけでなく製造公差、セン サーの異なった感度の高さ等も考慮に入れられる。従って、すでに製造者によっ て較正された磁気コンパスを使用することが必要不可欠ではない。 本発明に係る方法の場合には、それぞれの固定した箇所での重力ベクトルと地 磁場ベクトルとの間の傾角がシステムの一時的な位置に依存せずに一定のままで あるという事実が使用される。 本発明の理解を容易にするために図を示す。 図１は、電子工学的な磁気コンパスとソフトマグネティックな及びハードマグ ネティックなノイズ場発生源との配置を図式的に示している。 図２は、地磁場の測定に関連したベクトルを示す図である。 ポアソン(Ｐoisson)はすでに１９世紀に以下のような問題と取り組んでいた。 すなわち、磁力計(Ｍagnetometer)がそれ自体に磁気的な構成要素(magnetische Ｋomponenten)を有するシステムに取り付けられているときに、どのようにすれ ば実際に存在する磁場を算定できるか、という問題である。この状況を記述する ポアソンの公式は、そのような場合には測定された磁場は実際に存在する場の一 次関数(lineare Ｆunktion)であることをあらわしている。従って、アフィン・ マッピング(affine Ａbbildung)の問題である。これに関しては、すでに言及し た米国特許第４６８６７７２号明細書、第２欄、第２６行〜第３０行を参照のこ と。 一般的な場合には、測定された磁場が測定箇所でのソフトマグネティック(wei chmagnetisch)にひずまされた地磁場から及びハードマグネティックな成分(hart magnetische Ｋomponente)から構成されることが意味され得る。ソフトマグネテ ィックにひずまされた地磁場には、地磁場によ って誘導された磁力（磁気、Ｍagnetismus）がかかわっている。ハードマグネテ ィックな成分は、例えば当該システムにおける電流あるいは永久磁石によって発 生させられる磁力計の位置での一定の磁場を含む。ハードマグネティックな成分 は外場の変化によって影響を及ぼされ得ない。 数学的には、上述のポアソンの公式は、いくらか変形されて以下のように示さ れる： ここで、 前記の記述において並びに以下の説明において、ベクトル及び行列は太字で記 載されている。ベクトルは、通常デカルト座標系に適用される。行列は一般的に ３×３行列である。 図式的に示されている。当該装置ではハードマグネティックなノイズ 料とが協働し、ＤＭＣの測定値に影響を及ぼす。図２には、地磁場の測定に関連 したベクトルと水平面への関連した射影が示されている。車両に固定して配置さ れた磁気コンパスの場合には、目標方向 (Ｚielrichtung)が進行方向に対応する。 電子工学的な磁気コンパスは、三つの磁場センサーと二つの傾斜センサーとを 有する。その際、前者が三つの磁場成分を算定する。それらの三つの磁場成分は 必ずしも互いに直交している必要はない。 永久磁石及び電流は、磁場センサーの位置に固定した場を発生させる。この場 は、前記磁場成分によって構成される座標系の零点変位（ゼロシフト、Ｎullpun ktverschiebung）として現われる。 ソフトマグネティックな材料は、存在する場から、場方向における弱められた あるいは強められた場及びそれに対して垂直な方向における付加的な場成分を発 生させる。これは、場方向ｘ、ｙ、及びｚの 例えば製造者による較正(Ｋalibrierung)なしの「未処理の」(roh)磁力計のよ うな、感度の異なる三つの非直交の磁場センサーと零点変位とを備える磁力計を 観察するときもまた上記の同一の方程式が生じる。 そのような「未処理の」磁力計での磁場の測定は、数学的に次のように表現さ れ得る： わち単位ベクトル、ｏ_iはオフセットである。 ｆ_iｅ_i＝(Ｍ_i1,Ｍ_i2,Ｍ_i3)及びｏ_i＝ｂ_Oiと書くと、再び上記の方程式(1)が得 られる。 実際の地磁場の大きさを算定するために、上記のベクトル方程式(1)をｂ_Eに関 して解くことが必要不可欠である。逆行列(Ｉnversion)及び減法(Ｓubtraktion) によって以下の関係が得られる： 未知の量Ｍまたはｍ＝Ｍ^-1及びｂ_Oを方程式(1)あるいは(2)から算定するため に、すでに言及した米国特許第４６８６７７２号明細書により 能である。そこでは、測定された磁場ベクトルｂ_gemに対して常にさらに地磁場 ベクトルｂ_Eも明確に知られていることを出発点とする。ベクトル方程式(1)及び (2)は、それぞれ、未知数Ｍ_ijまたはｍ_ij及びｂ_Oiについての線形方程式システ ム(lineares Ｇleichungssystem)である。少なくとも四つの幾何学的に異なる位 置での測定によって、地磁場が知られている場合に、Ｍ_ij及びｂ_Oiが、線形方程 式システムを解くための基本的な方法を用いてすぐにわかる。 前述の米国特許第４５３９７６０号明細書に記載されている解決策の場合には 、地磁場の大きさ(Ｂetrag)は磁力計の位置に依存しないという事実を出発点と する。量ｍ_ij及びｂ_Oが正しく算定されていれば、磁力計のそれぞれの位置にお いて同一の長さを有する地磁場ベクトルｂ_Eが生じる。従って： ここで、Ｕ^T＝Ｕは、Ｕが対称行列であることを意味する。 この式(3)によって、求められる行列ｍの積Ｕ＝ｍ^Tｍだけが算出され得ること がすぐにわかる。この行列が対称であることが受け入れられるとき及び対角成分 (Ｄiagonale)が正であるときにだけこの行列の成分が算出され得る。しかし、最 初に述べた仮定(Ａnnahme)は、まれなケー スでしか該当しない。というのは、それが、ソフトマグネティックな対称を意味 したからである。このソフトマグネティックな対称は、例えば飛行機あるいは動 力車のような技術的な装置の場合には、与えられることが極めてありそうにない 。 本発明により、量ｂ_Eを決めるために、すなわち前記のベクトル方程式を解く ために、地理学上の同一箇所での測定システムのそれぞれの位置で水平面と地磁 場との間の角度、すなわち傾角(Ｉnklinationswinkel)が一定のままであるとい う事実を利用した式を出発点とする。そのとき、このことは明らかに重力ベクト ルｇの方向と地磁場ベクトルｂ_Eとの間の角度にも当てはまる。従って、以下の ように書くことができる： ここで、 この方程式は、重力場の方向における、従って水平面に対して垂直な磁場ベク トルの成分がすべてのシステム位置において同一のままであることを意味する。 この方程式において、量ｍはリニア（一次）であり、積としては存在しない。 重力ベクトルｇは、傾斜センサーを使って算定され得る。それ故に、この量ｍは 、前述の米国特許第４６８６７７２号明細書に示されているような場の測定を実 行する必要なしに、あるいは前述の米国特許第４５３９７６０号明細書に示され ているような特別な対称条件を前提とする必要なしに直接に算出され得る。 決定すべきパラメータの数は、以下のようになる： ｍ ＝3×3＝9 ｂ_O＝3 ｂ_g＝1 したがって、全部で9+3+1=13のパラメータが生じる。それらの決定のために少な くとも１３個の方程式が必要とされる。 例えば前述の米国特許第４５３９７６０号明細書（第４欄、第３行以下）から 知られているように、さらに、任意のスケールファクター(ＭaBstabsfaktor)が 選ばれ得る。可能なスケール固定(ＭaBstabsfestlegung)は、例えば以下の一欄 、すなわち ｂ_g ＝const. ｍ₁₁ ＝const. ｍ₁₁＋ｍ₂₂＋ｍ₃₃＝const. ｍ₁₁ ²＋ｍ₂₂ ²＋ｍ₃₃ ²＝const. ｍ₁₁ ²＋ｍ₁₂ ²＋...+ｍ₃₃ ²＝const. det ｍ＝const. の一つ、あるいは別の適当なものであるだろう。その際、定数としては１を選ん でもよい。スケール固定によって、パラメータの数は１だけ減らされる。その結 果、なお１２個のパラメータが残っているままである。これらのパラメータに対 して対応する数の方程式が必要とされる。そのために、磁気コンパスの最初の較 正の際に、１２個の異なる幾何学的な位置ｊ＝１，．．．，１２が占められる。 それらの位置において三つの磁場成分と両方の傾斜角(Ｎeigungswinkel)とが測 定される。 以下のような線形方程式システムが得られる： ここで、ｕ_O＝ｍｂ_Oであり、例えばｂ_g＝1である。その場合、明示的に以下のよ うになる： １２より多くの方程式が与えられれば、ここで例えば最も良い適合が、一般的に 知られている最小二乗偏差の方法(das Ｖerfahren der kleinsten quadratische n Ａbweichung)によって決定され得る。 本発明のフレームにおいて、同一の固定した箇所での傾角またはｂ_gの不変性 の事実を磁場ベクトルの長さ（ｂ_E＝｜ｂ_E｜）の存在する不変性の事実とともに 利用することの可能性も考慮に入れられている。換言すると、このことは、上述 の方程式(5)だけでなく方程式(3)も当該方法のために使用され得ることを意味す る。そのとき、付け加えてさらにパラメータｂ_Eが決定される必要がある。その 結果、パラメータの数Ｎ、従って必要不可欠な測定の数ＮがＮ＝１３へ１だけ増 加する。 有利には、方程式(3)及び(5)の使用によって、与えられているデータのより良 い利用が得られる。なぜならば、それぞれの測定が二つの方程式において使用さ れるからである。それゆえに、必要不可欠な測定の数は半分になる。この方法の 実行のためには、統計学上の補整計算(statistische Ａusgleichsrechnung)が要 求される。この補整計算は、例えば最小二乗偏差の方法に基づいてされ得る。 以下に記述される方法では、上述の方程式(3)及び(5)が統計学上の適切な方法 で使用される。 次に、磁場センサーシグナルμ_ij；ｉ＝１，２，３と地磁場ベクトルｂ_Ejとの 間にある関係を詳細に説明する。これは、以下の方程式によっ で表わされ得る： 地磁場ベクトルの成分ｂ_Ejは、磁気コンパスのさまざまの位置においては知ら れていない。しかし、それらを傾斜センサーを使って部分的に示すことは可能で ある。その際、二つの傾斜センサーが、それぞれ、重力場の１に規格化されたベ クトルの一つの成分を測定する。このベクトルは鉛直方向(Ｌotrichtung)に延び る。 その際、α_jは方位角、すなわち、水平面における空間に固定された座標 系に対してのセンサー固有の座標系の回転角である。 パラメータＭ₁₁...Ｍ₃₃,μ₀₁,μ₀₂,μ₀₃,α₁...α_N及びｉの値はそれ自体公知の 最適化計算(Ｏptimierungsrechnung)の使用のもとで決定され得る。その際、例 えば、統計学上の和 についてミニマムが要求される。本発明のフレームにおいては、例えばカルマン フィルター(Ｋalmanfilter)、ファジー法(Ｆuzzy-Ｖerfahren)、あるいはニュー ロン・ネットワーク(neuronale Ｎetze)によってのようにさまざまな別の解法が 可能である。 量ΔＭ＝０であると、本質的に数学的な簡単化が生じることがわかる。これは 、ソフトマグネティックなノイズ場が考慮に入れられる必要がない場合である。 そのとき、Ｍは単位行列に合致する。 以上では、傾斜センサーを引き合いにだして説明した。これらの代わりに、二 つの直交して取り付けられたエンコーダも使用され得る。これらによって、角度 が基準点(Ｂezugspunkt)に関して測定されるだろう。 しかしそのとき、実際的な実行の場合には、固定取り付けされた取り付け具が設 けられねばならないだろう。当該取り付け具に対して、磁気コンパス及び妨害す る系が相対回転するだろう。 望ましい傾斜角を、システムに取り付けられた二つの回転率センサー(Ｄrehra tensensoren)、すなわちジャイロ(Ｇyros)を用いて決定してもよいだろう。回転 率情報(Ｄrehrateninformationen)の統合によって回転角が導出され得る。

───────────────────────────────────────────────────── 【要約の続き】 である。また、Ｎが当該ベクトル方程式の決定すべきパ ラメータの数に少なくとも等しくなければならない。

Claims (1)

1. 【特許請求の範囲】 １．装置に固定された磁気的な材料及び電流の磁場によって乱されている可能性 のある地磁場の方向を、三つの磁場センサー及び二つの傾斜感知装置をもってい る一つの電子工学的な磁気コンパスを用いて算定するための方法にして、 当該電子工学的な磁気コンパスがＮ個の異なる空間的な位置に配置され、 これらのＮ箇所のそれぞれにおいて、前記傾斜感知装置の信号及び前記磁場 センサーの信号が測定され、且つこれらの信号から傾斜値及び磁場値が算定され 、 これらの傾斜値及び磁場値に基づいて地磁場ベクトルの大きさが、ベクトル 方程式 その際、 を用いて算出され、 その際、Ｎが当該ベクトル方程式の決定すべきパラメータの数に少なくとも 等しくなければならない、前記の方法。 ２．以下のスケール固定、すなわち ｂ_g ＝const.；ｍ₁₁＝const.；ｍ₁₁＋ｍ₂₂＋ｍ₃₃＝const.； ｍ₁₁ ²＋ｍ₂₂ ²＋ｍ₃₃ ²＝const.； ｍ₁₁ ²＋ｍ₁₂ ²＋...＋ｍ₃₃ ²＝const.；det ｍ＝Const.； のうちの一つが選択されることを特徴とする、請求項１に記載の方法。 ３．ソフトマグネティックなノイズ場を考慮に入れる必要がない場合に、ΔＭ＝ ０とおかれることを特徴とする、請求項１に記載の方法。 ４．異なる位置での測定の数Ｎが与えられている方程式の数よりも大きいときに 、前記磁場ベクトルが最適化計算を用いて算定されることを特徴とする、請求項 １に記載の方法。 ５．地磁場ベクトルの大きさが、前記算定された傾斜値及び磁場値に基づいて、 以下の方程式 を付加的に使用して、統計学的な最適化法を用いて算定されることを特徴とす る、請求項１に記載の方法。 ６．傾斜感知装置として傾斜センサーが使用されることを特徴とする、請求項１ 〜５のいずれか一項に記載の方法。 ７．傾斜感知装置として回転率センサー（ジャイロ）が使用され、回転率情報の 統合によって傾斜角が導出されることを特徴とする、請求項１〜５のいずれか一 項に記載の方法。 ８．傾斜感知装置としてエンコーダが使用され、当該エンコーダによって基準点 に関して角度が測定されることを特徴とする、請求項１〜５のいずれか一項に記 載の方法。