(57)【特許請求の範囲】
【請求項１】
少なくともＮ個の三軸磁力計（Ｍｉ，ｊ）であって、前記磁力計（Ｍｉ，ｊ）の既知の相対位置を維持するために自由度なしで互いに機械的に接続される少なくともＮ個の三軸磁力計（Ｍｉ，ｊ）の配列に対する少なくとも１個の可動磁性体（ＯＭＭｋ）の位置を特定する方法であって、Ｎは少なくとも２に等しい整数であり、前記方法は、
−  現時点（ｔ）における前記磁力計（Ｍｉ，ｊ）の測定値（Ｂｉ，ｊｍｅｓ）をロードするステップと、
−  平均地球磁場および前記可動磁性体（ＯＭＭｋ）に起因する平均磁場を表す前記磁力計（Ｍｉ，ｊ）の前記測定値（Ｂｉ，ｊｍｅｓ）の加重平均（＜Ｂｍｅｓ＞）を計算するステップと、
−  修正測定値（Ｂｉ，ｊｍｅｓ−＜Ｂｍｅｓ＞）を取得するために、前記加重平均（＜Ｂｍｅｓ＞）を前記測定値（Ｂｉ，ｊｍｅｓ）の各々から減算するステップと、
−  前記１個以上の可動磁性体（ＯＭＭｋ）の位置を特定するためのフィルタリング演算（ＦＬ）への入力として、前記修正測定値（Ｂｉ，ｊｍｅｓ−＜Ｂｍｅｓ＞）および前記現時点（ｔ）における前記１個以上の可動磁性体（ＯＭＭｋ）の位置（Ｐｏｓｔ，Ｍｔ）をロードするステップと、
−  前記位置フィルタリング演算（ＦＬ）を実行するステップを含み、前記実行するステップが、
−  前記磁力計（Ｍｉ）により伝送されたデータの推定値（Ｂｉ，ｊｅｓｔ）の加重平均（＜Ｂｅｓｔ＞）を計算することを含み、推定値（Ｂｉ，ｊｅｓｔ）が、推定アルゴリズムを用いて計算され、推定アルゴリズムが、磁石の位置、向き、および磁気モーメントを空間のあらゆる点で生じた磁場に関連付ける測定モデルの使用に基づき、前記実行するステップがさらに、
−  前記加重平均（＜Ｂｅｓｔ＞）を前記推定データ（Ｂｉ，ｊｅｓｔ）の各々から減算することと、
−  後続時点（ｔ＋１）における前記１個以上の可動磁性体（ＯＭＭｋ）の位置（Ｐｏｓｔ＋１，Ｍｔ＋１）を出力として伝送することとを含む、方法。
【請求項２】
前記加重平均は、前記Ｎ個の磁力計（Ｍｉ，ｊ）の各々に対して同一である重み係数を使用する、請求項１に記載の方法。
【請求項３】
前記加重平均は、前記Ｎ個の磁力計（Ｍｉ，ｊ）のうち、前記１個以上の可動磁性体（ＯＭＭｋ）に最も近いＰ個の磁力計に対してより低い重み係数を使用し、ＰはＮよりも小さい整数である、請求項１に記載の方法。
【請求項４】
前記加重平均は、前記Ｐ個の磁力計に対してゼロの重み係数を使用する、請求項３に記載の方法。
【請求項５】
前記加重平均は、前記Ｎ個の磁力計（Ｍｉ，ｊ）のうち、前記１個以上の可動磁性体（ＯＭＭｋ）から最も遠い前記磁力計以外に対してゼロの重み係数を使用する、請求項１に記載の方法。
【請求項６】
少なくとも１個の可動磁性体（ＯＭＭｋ）の位置を特定するシステムであって、
−  少なくともＮ個の三軸磁力計（Ｍｉ，ｊ）の配列を含み、前記磁力計（Ｍｉ，ｊ）は、前記磁力計（Ｍｉ，ｊ）の既知の相対位置を維持するために自由度なしで互いに機械的に接続され、Ｎは少なくとも２に等しい整数であり、前記システムはさらに、
−  前記配列の前記磁力計（Ｍｉ，ｊ）および位置フィルタ（ＦＬ）により伝送された測定値（Ｂｉｍｅｓ）から前記１個以上の可動磁性体（ＯＭＭｋ）の位置（Ｐｏｓｔ，Ｍｔ）を判定するために好適な電子処理ユニット（１１）を含み、
前記電子処理ユニット（１１）は、
−  現時点（ｔ）における前記磁力計（Ｍｉ，ｊ）の測定値（Ｂｉ，ｊｍｅｓ）をロードする手段と、
−  平均地球磁場および前記可動磁性体（ＯＭＭｋ）に起因する平均磁場を表す前記磁力計（Ｍｉ，ｊ）の前記測定値（Ｂｉ，ｊｍｅｓ）の加重平均（＜Ｂｍｅｓ＞）を計算する手段と、
−  修正測定値（Ｂｉ，ｊｍｅｓ−＜Ｂｍｅｓ＞）を取得するために、前記加重平均（＜Ｂｍｅｓ＞）を前記測定値（Ｂｉ，ｊｍｅｓ）の各々から減算する手段とを含み、
前記位置フィルタ（ＦＬ）は、
−  前記可動磁性体（ＯＭＭ）の位置を特定するためのフィルタリング演算（ＦＬ）への入力として、前記修正測定値（Ｂｉ，ｊｍｅｓ−＜Ｂｍｅｓ＞）および前記現時点（ｔ）における前記１個以上の可動磁性体（ＯＭＭｋ）の位置（Ｐｏｓｔ，Ｍｔ）をロードする手段と、
−  前記磁力計（Ｍｉ，ｊ）により伝送されたデータの推定値（Ｂｉ，ｊｅｓｔ）の加重平均（＜Ｂｅｓｔ＞）を計算する手段とを含み、推定値（Ｂｉ，ｊｅｓｔ）が、推定アルゴリズムを用いて計算され、推定アルゴリズムが、磁石の位置、向き、および磁気モーメントを空間のあらゆる点で生じた磁場に関連付ける測定モデルの使用に基づき、前記位置フィルタ（ＦＬ）はさらに、
−  前記加重平均（＜Ｂｅｓｔ＞）を前記推定データ（Ｂｉ，ｊｅｓｔ）の各々から減算する手段と、
−  後続時点（ｔ＋１）における前記１個以上の可動磁性体（ＯＭＭｋ）の位置（Ｐｏｓｔ＋１，Ｍｔ＋１）を出力として伝送する手段とを含む、システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、少なくとも１個の可動磁性体の位置を特定する方法および関連するシステムに関する。
【背景技術】
【０００２】
磁性体の位置を特定する方法および装置に関する文献、例えば文献仏国特許第２９８８８６２号明細書等が知られている。
【０００３】
文献仏国特許第２９８８８６２号明細書は、Ｎ個の三軸磁力計であって、これらの磁力計の各々間の既知の距離を維持するために自由度なしで互いに機械的に接続されるＮ個の三軸磁力計を含む、三軸磁力計の配列であって、Ｎは５よりも大きい整数である、配列により取得された測定値から磁性体の位置を特定する方法を開示している。磁力計の配列は、任意選択的に重ね合わされた一枚以上の用紙等の筆記媒体用の荷重面を含む装置内に配置される。
【０００４】
磁性体の位置は、その配置（磁力計の配列に対して固定された座標系内での（ｘ、ｙ、ｚ）座標）と、磁性体の向きを表す磁気モーメント（同一の座標系内での（Ｍｘ，Ｍｙ，Ｍｚ））の値からなると考えられる。
【０００５】
この種の方法で使用される測定値へのバイアスの導入は、位置を不正確にするか、または不安定にさえする。
【０００６】
このような位置特定装置の動的な使用に関連して、すなわち筆記媒体を搭載した位置特定装置をユーザーが動かすことができる場合、その向きを修正することができ、これは、磁力計の配列により測定される座標系内での地球磁場の回転と、従って全ての磁力計の測定におけるバイアスとに直接つながる。
【０００７】
従って、均一な磁場内での位置特定装置の移動により、このような位置特定装置を使用することがほぼ不可能になる。
【０００８】
このような装置が手書き文字をデジタル化するシステムで使用される場合、そのようなシステムによる位置の判定に要する精度は１ｍｍ未満であるため、磁力計の座標系内での地球磁場におけるわずかな変化であっても、求める精度よりも大きい測定誤差に直接つながる。
【０００９】
従来技術において、起動時に、全てのセンサにおいてまたは一切の磁性体が存在しない状態で測定バイアスを除去することが知られている。
【００１０】
このような方法では、磁性体の位置特定実行中の移動が不可能であり、従って、移動的使用に関して述べた潜在的な問題を解決しない。特に、筆記デジタル化アプリケーションに関連して、ユーザーは、自らの姿勢および動きを修正しながら可能な筆記角度を得るために位置特定装置を動かす場合が多い。
【００１１】
また、磁力計を可動磁性体から十分な距離に置くことにより、この磁力計が可動磁性体により生じた磁場の影響を受けないようにする方法が知られている。この遠隔磁力計により測定される密封は、地球磁場により誘導される磁性成分を抑制するのに役立つ。
【００１２】
上述のような方法は、小型の装置が使用されることを許容せず、従って携帯機器に適していない。更に、上述の方法は、この遠隔磁力計の測定誤差、特に地球以外の磁場による遠隔磁力計の誤動作または撹乱に極めて敏感である。
【先行技術文献】
【特許文献】
【００１３】
【特許文献１】仏国特許第２９８８８６２号明細書
【発明の概要】
【発明が解決しようとする課題】
【００１４】
本発明の目的の１つは、これらの問題を緩和することである。
【課題を解決するための手段】
【００１５】
本発明の一態様によれば、少なくともＮ個の三軸磁力計であって、これらの磁力計の既知の相対位置を維持するために自由度なしで互いに機械的に接続される少なくともＮ個の三軸磁力計の配列に対する少なくとも１個の可動磁性体の位置を特定する方法が提供され、Ｎは少なくとも２に等しい整数であり、本方法は、
−  現時点における磁力計の測定値をロードすることと、
−  平均地球磁場および可動磁性体に起因する平均磁場を表す磁力計の前記測定値の加重平均を計算することと、
−  修正測定値を取得するために、前記測定値の前記加重平均を前記測定値の各々から減算することと、
−  １個以上の可動磁性体の位置を特定するためのフィルタリング演算への入力として、前記修正測定値および現時点における１個以上の可動磁性体の位置をロードすることと、
−  位置フィルタリング演算を実行することであって、
−  前記磁力計により伝送されたデータの推定値の前記加重平均を計算することと、
−  前記推定値の前記加重平均を前記推定データの各々から減算することと、
−  後続時点における１個以上の可動磁性体の位置を出力として伝送することと
を含むステップを含む、実行することと
を含むステップを含む。
【００１６】
本方法により、上述の磁気位置特定装置を、装置の向きを修正しながらおよび／または移動しながら使用することが可能になる。
【００１７】
また、本方法により、装置に対する均一な磁場（地球磁場がその一例である）の影響を抑制して、全ての磁力計に均一な磁場をもたらす撹乱の発生源を除去することが可能になる。大部分の発生源は磁気双極子に例えることができる。装置と撹乱の発生源との距離が十分に大きい場合、検出される磁場は均一な磁場となり、従って本方法により修正される。
【００１８】
磁気撹乱の発生源としてＤＣモーターが知られており、例えば１ガウスよりも強い強磁場を生成する。従って、本方法により、これらの装置を列車または路面電車に使用することができる。
【００１９】
「位置を特定する」とは、例えばデカルト座標（３本の直交軸に沿った成分）で表された磁性体の位置および磁気モーメントを判定することにより、可動磁性体の位置および向きを判定することを意味する。一変形形態として、可動磁性体の位置および向きは、磁気モーメントを球面座標系、すなわちノルムおよび２個の角度で表すことにより判定することができる。この代替案により、磁気モーメントを既知の値に初期化することにより可動磁性体の磁気モーメントを設定でき、従って位置および２個の角度のみを推定すればよく、利点として位置フィルタにより推定すべき未知数の個数を減らすことができる。
【００２０】
一実施形態において、前記加重平均は、Ｎ個の磁力計の各々に対して同一である重み係数を使用する。この方法は、必要とされる計算能力が異なる重みを使用する場合よりも低くて済む利点がある。具体的には、３ｎ回の乗算（Ｎはセンサの個数である）が冗長になる。
【００２１】
一実施形態によれば、前記加重平均は、Ｎ個の磁力計のうち、各々の可動磁性体に最も近いＰ個の磁力計に対してより低い重み係数を使用し、ＰはＮよりも小さい整数である。
【００２２】
従って、可動磁性体により最も影響を受け易いかまたは飽和することもあるＰ個の磁力計の影響は無視され、限定的である。
【００２３】
複数の可動部が存在する場合、磁力計と各可動磁性体との距離が配列の各磁力計について判定され、次いで、判定された距離が全ての磁力計および可動磁性体について昇順に分類され、分類された距離が最も短いＰ個の磁力計に対して低い相関係数が使用される。
【００２４】
例えば、前記加重平均は、前記Ｐ個の磁力計に対してゼロの重み係数を使用する。
【００２５】
従って、可動磁性体により最も影響を受け易いかまたは飽和することもあるＰ個の磁力計は無視される。
【００２６】
一実施形態において、前記加重平均は、Ｎ個の磁力計のうち、１個以上の可動磁性体から最も遠い磁力計以外に対してゼロの重み係数を使用する。
【００２７】
従って、測定された地球磁場の平均は、１個以上の可動磁性体から最も遠い磁力計により測定された磁場に等しく設定され、従って１個以上の可動磁性体により撹乱される度合いが最も少ない磁場に設定される。
【００２８】
本発明の別の態様によれば、少なくとも１個の可動磁性体の位置を特定するシステムも提供され、前記システムは、
−  少なくとも三軸磁力計であって、これらの磁力計の既知の相対位置を維持するために自由度なしで互いに機械的に接続される少なくとも三軸磁力計の配列であって、Ｎは少なくとも２に等しい整数である、配列と、
−  配列の磁力計および位置フィルタにより伝送された測定値から１個以上の可動磁性体の位置を判定するために好適な電子処理ユニットと
を含み、
前記電子処理ユニットは、
−  現時点における磁力計の測定値をロードする手段と、
−  平均地球磁場および可動磁性体に起因する平均磁場を表す磁力計の前記測定値の加重平均を計算する手段と、
−  修正測定値を取得するために、前記測定値の前記加重平均を前記測定値の各々から減算する手段と
を含み、
位置フィルタは、
−  可動磁性体の位置を特定するためのフィルタリング演算への入力として、前記修正測定値および現時点における１個以上の可動磁性体の位置をロードする手段と、
−  前記磁力計により伝送されたデータの推定値の加重平均を計算する手段と、
−  前記推定値の前記加重平均を前記推定データの各々から減算する手段と、
−  後続時点における１個以上の可動磁性体の位置を出力として伝送する手段と
を含む。
【００２９】
本発明は、添付図面に例示する完全に非限定的な例を使用して記述するいくつかの実施形態を調べることでより良く理解されるであろう。
【図面の簡単な説明】
【００３０】
【図１】公知の従来技術による位置特定装置を概略的に示す。
【図２】公知の従来技術による位置特定装置を概略的に示す。
【図３】本発明の一態様による方法を概略的に示す。
【発明を実施するための形態】
【００３１】
全ての図面を通じて同一の参照符号で示す要素は同一である。
【００３２】
以下の記述において当業者に公知の特徴および機能は詳述しない。
【００３３】
図１に、従来技術による、例えば文献仏国特許出願公開第２９８８８６２Ａ１号明細書に記述されている可動磁性体ＯＭＭの位置を特定する装置１を示す。
【００３４】
この場合、位置特定装置１を使用して、電子装置２、例えば画面２上の画像の表示を制御することができる制御装置３に接続された画面２を制御する。
【００３５】
図示する例において、装置３は画面２上のカーソル４の位置および向きを制御する。例えば、カーソル４は、直方体等の三次元形状を有する。
【００３６】
可動磁性体ＯＭＭｋは、外部磁場および非磁性筆記具６が存在しない場合であっても非ゼロの磁気モーメントを有する永久磁石５を含む。非磁性とは、測定可能な磁気特性を有していない材料で作られた筆記具を意味する。例えば、筆記具６は木またはプラスチック製の鉛筆、ペン、またはインク消しである。筆記具６の位置は永久磁石５の位置から特定することができる。筆記具６により、特に永久磁石５のサイズが小さい場合に永久磁石５を扱い易くする。典型的に、永久磁石５および筆記具６の組み合わせは、装置１と自由度なしで固定された直交ＸＹＺ座標系に存在する人間の手により直接自由に動かすことができる。この場合、ＸおよびＹ方向は水平（筆記媒体の平面）であり、Ｚ方向は垂直である。この目的のため、この組み合わせの重量は１キロ未満、好適には２００ｇ未満である。この組み合わせの寸法は、ユーザーの片手でつまみ上げて動かすのに十分小さい。１個以上の可動磁性体ＯＭＭｋは、直交ＸＹＺ座標系内で自由に動かすことができる。
【００３７】
例えば、磁石の保磁力は１００Ａ．ｍ−１または５００Ａ．ｍ−１よりも高い。例えば、それは強またはフェリ磁性体で作られている。永久磁石５は長方形の形状をなしている。図１において、永久磁石５の磁気モーメントの方向を、この物体の長手方向に平行な矢印で示す。一変形形態として、永久磁石５は、筆記具６を囲むリングの形状をなしていてよい。永久磁石５で最長寸法を以下でＬと表記する。
【００３８】
永久磁石５の強度は、典型的に０．０１Ａ．ｍ２または０．１Ａ．ｍ２よりも強い。本実施形態において、永久磁石５は、自由度なしで筆記具６に固定されている。
【００３９】
位置特定装置１により、永久磁石５のＸＹＺ座標系内での位置を特定することができる。ここで位置を特定するとは、ＸＹＺ座標系内での永久磁石５のｘ、ｙ、ｚ位置を判定し、かつＸＹＺ座標系のＸ、ＹおよびＺ軸に関する永久磁石５の向きを判定することを意味する。例えば、永久磁石５の向きは、座標系のＸ、ＹおよびＺ軸に対する、ＸＹＺ座標系内での永久磁石５の磁気モーメントのデカルト座標Ｍｘ，Ｍｙ、およびＭｚにより、またはＸＹＺ座標系内での永久磁石５の磁気モーメントの球面座標Ｍ、θ、φにより表すことができる。
【００４０】
位置特定装置１は、Ｎ個の三軸磁力計Ｍｉｊの配列を含む。図１において、波形の垂直線は、いくつかの位置特定装置１が示されていないことを示す。
【００４１】
典型的に、Ｎは５よりも大きく、好適には１６もしくは３２よりも大きく、または６４に等しくてもよい。
【００４２】
本実施形態において、磁力計Ｍｉｊは、行および列に配置されてマトリクス配列または単に配列を形成する。この場合、このマトリクス配列は８行および８列を含む。添え字ｉおよびｊは、それぞれこのマトリクス配列の行および列を識別し、その交差部に磁力計Ｍｉｊが見出されれる。
【００４３】
図１において、行ｉの磁力計Ｍｉ１、Ｍｉ２、Ｍｉ３、Ｍｉ４およびＭｉ８のみを示す。互いに相対的な磁力計Ｍｉｊの位置を図２に更に詳細に示す。
【００４４】
各磁力計Ｍｉｊは、他の磁力計と自由度なしで接合されている。この目的のため、磁力計Ｍｉｊは、剛性パネル８の背面７に自由度なしで固定されている。この剛性パネル８は、永久磁石５に向けられた前面９を有する。パネル８は、固い非磁性材料で作られている。例えば、パネル８はガラス製であってよい。
【００４５】
各磁力計Ｍｉｊは、永久磁石５により生じた磁場の方向および強度を測定する。これを行うため、各磁力計Ｍｉｊは、この磁力計の３本の測定軸に沿ってこの磁力計Ｍｉｊの永久磁石５により生じた磁場の直交射影のノルムを測定する。この場合、これらの３本の測定軸は互い直交している。例えば、磁力計Ｍｉｊの各々の測定軸は座標系のＸ、Ｙ、およびＺ軸に平行である。
【００４６】
各磁力計Ｍｉｊは、データバス１０を介して処理ユニット１１に接続されている。
【００４７】
処理ユニット１１は、磁力計Ｍｉｊの測定値から、ＸＹＺ座標系内での永久磁石５の位置および向きを判定することができる。この目的のため、装置１１は、データ記憶媒体に格納された命令を実行可能なプログラム可能電子プロセッサ１２を含む。装置１１は、従って、プロセッサ１２による図３の方法の実行に必要な命令を含むメモリ１３も含む。
【００４８】
特に、装置１１は、複数の磁力計Ｍｉｊの測定値に永久磁石５、すなわち同じくＸＹＺ座標系内での筆記具６の位置および向きを表すパラメータを関連付ける物理数学的モデルを実装している。
【００４９】
このモデルは、非線形推定値フィルタ、例えば拡張カルマンフィルタ形式で実装される。
【００５０】
このモデルは、典型的に静磁場物理方程式から構築される。このモデルを構築するために永久磁石５を磁気双極子で近似している。この近似により、永久磁石５と磁力計Ｍｉｊとの距離が２Ｌよりも大きく、好適には８Ｌよりも大きい場合、生じる誤差が極めて小さい（Ｌは永久磁石５の最大寸法である）。典型的には、Ｌは、２０ｃｍよりも小さく、好適には１０または５ｃｍ未満である。
【００５１】
装置１１はまた、永久磁石５の位置および測定された向きをインターフェース１４に渡すことができる。
【００５２】
制御装置３は、このインターフェース１４を介して装置１１に接続されている。
【００５３】
図２は、位置特定装置１のいくつかの磁力計Ｍｉｊを示す。これらの磁力計Ｍｉｊは、Ｘ方向に平行な行ｉに整列配置されている。これらの磁力計はまた、マトリクス配列または単に配列を形成するためにＹ方向に平行な列ｊに整列配置されている。行ｉおよび列ｊは、添え字の昇順に配置されている。
【００５４】
磁力計Ｍｉｊの中心は、行ｉおよび列ｊの交差箇所に位置している。磁力計の中心は、磁場がこの磁力計により測定される点に対応している。この場合、添え字ｉおよびｊは、区間［１；８］に属している。
【００５５】
行ｉに沿って直接隣接する２個の磁力計Ｍｉｊ、Ｍ１，ｊ＋１の中心は、既知の距離ｄｉ，ｊ，ｊ＋１だけ離れている。同様に、所与の列ｊに沿って直接隣接する２個の磁力計Ｍｉｊ、Ｍｉ＋１，ｊの中心は、既知の距離ｄｊ，ｉ，ｉ＋１だけ離れている。
【００５６】
ここに記述する特定の場合において、いずれの行ｉについても距離ｄｉ，ｊ，ｊ＋１は同一である。従って、この距離をｄｊと表記する。同様に、いずれの列ｊについても２個の磁力計間の距離ｄｊ，ｉ，ｉ＋１は同一である。従って、この距離をｄｉと表記する。記述する例において、距離ｄｉおよびｄｊは共にｄに等しい。
【００５７】
典型的に、距離ｄが１〜４ｃｍであるのは、１
−  永久磁石の強度が０．５Ａ．ｍ２であり、
−  磁力計の感度が４×１０−７Ｔであり、および
−  磁力計Ｍｉｊの個数が６４であるときである。
【００５８】
図３は、本発明の一態様による方法の一実施形態を示す。本方法のステップは、電子処理ユニット１１により実装される。
【００５９】
各測定ステップにおいて、平均地球磁場と可動磁性体ＯＭＭｋに起因する平均磁場との和を表す各磁力計Ｍｉ，ｊの前記測定値Ｂｉ，ｊｍｅｓの加重平均＜Ｂｍｅｓ＞を計算する。
【００６０】
加重平均＜Ｂｍｅｓ＞は、非限定的に、重み係数が全て等しい場合に算術平均であってよい。
【００６１】
従って、地球磁場は、１個以上の可動磁性体ＯＭＭｋにより生成された平均磁場と同時に除去される。
【００６２】
地球磁場またはシステムに対する均一な磁気撹乱を抑制することで求める修正が得られるのに対し、強磁性物質により生じた磁場の平均を抑制することで誤差が生じる。
【００６３】
次に、生じた誤差を考慮に入れて位置特定アルゴリズムの精度を下げないようにする。
【００６４】
従って、平均地球磁場および可動磁性体ＯＭＭｋに起因する平均磁場を表す磁力計Ｍｉｊの測定値
【数１】
0006487050.tif 000002
の加重平均
【数２】
0006487050.tif 000003
を計算する。
【００６５】
次に、測定値Ｂｉ，ｊｍｅｓの加重平均
【数３】
0006487050.tif 000004
を各々の測定値
【数４】
0006487050.tif 000005
から減算して修正測定値
【数５】
0006487050.tif 000006
を得る。
【００６６】
例えば、加重平均が算術平均である場合、
【数６】
0006487050.tif 000007
が成立し、ここで、ａ×ｂ＝Ｎである。
【００６７】
次いで位置フィルタリング演算（ＦＬ）は、入力として、修正測定値
【数７】
0006487050.tif 000008
および現時点ｔでの１個以上の可動磁性体（ＯＭＭｋ）の位置
【数８】
0006487050.tif 000009
を使用する。
【００６８】
位置フィルタリング演算は、前記磁力計（Ｍｉｊ）により伝送したデータの推定値
【数９】
0006487050.tif 000010
の推定加重平均
【数１０】
0006487050.tif 000011
の計算を使用する。この場合、使用する加重平均は、測定された加重平均（すなわち、測定値の重み付けおよび推定値の重み付けにおいて磁力計に各々関連付けられた同一の重み係数であり、１個の重み係数が各磁力計に関連付けられている）に使用したものと同一でなければならない。
【００６９】
推定アルゴリズムは、磁石の位置、向き、および磁気モーメントを空間のあらゆる点で生じた磁場に関連付ける測定モデルの使用に基づいている。
【００７０】
推定加重平均の計算は、本実施形態では、同一の加重平均、この場合には算術平均
【数１１】
0006487050.tif 000012
を使用し、ここで、ａ×ｂ＝Ｎである。
【００７１】
次に、前記加重平均
【数１２】
0006487050.tif 000013
を前記推定データＢｉ，ｊｅｓｔの各々から減算する。
【００７２】
位置フィルタＦＬは、次いで、出力として後続時点ｔ＋１における１個以上の可動磁性体ＯＭＭｋの位置
【数１３】
0006487050.tif 000014
を伝送する。
【００７３】
一変形形態として、前記加重平均は、Ｎ個の磁力計Ｍｉ，ｊのうち、各々の可動磁性体ＯＭＭｋに最も近いＰ個の磁力計に対してより低い重み係数を使用し、ＰはＮよりも小さい整数である。当然のことながら、測定値の加重平均は依然として推定値の加重平均と同一である。
【００７４】
従って、可動磁性体により最も影響を受け易いかまたは飽和することもあるＰ個の磁力計の影響は無視され、限定的である。
【００７５】
加重平均は、可動部磁性体により最も影響を受け易いかまたは飽和することもあるＰ個の磁力計の影響を完全に除去して無視できるように、前記Ｐ個の磁力計に対してゼロの重み係数を使用してもよい。
【００７６】
一変形形態として、加重平均は、地球磁場同等の平均を１個以上の可動磁性体により最も影響を受けにくい磁力計の測定値に設定するために、Ｎ個の磁力計のうち、１個以上の可動磁性体から最も遠い磁力計以外に対してゼロの重み係数を使用してもよい。
【００７７】
本実施形態において、計算された平均値は磁力計の測定値と厳密に一致する。
【００７８】
上述の方法のステップは、入力データに作用して出力データを生成することにより本発明の機能を実行するために、コンピュータプログラムを実行する１個以上のプログラム可能プロセッサにより実行され得る。
【００７９】
コンピュータプログラムは、コンパイラまたはインタープリタ言語を含む任意の形式のプログラミング言語で記述さていてよく、コンピュータプログラムは、スタンドアローンプログラムまたはサブプログラムとして、コンピュータ環境での使用に適した要素または他の装置を含む任意の形式で提供することができる。コンピュータプログラムは、１台のコンピュータ、または１箇所もしくは複数箇所にわたり分散されて通信ネットワークにより互いに接続された複数台のコンピュータ上で実行するために提供することができる。
【００８０】
本発明の好適な実施形態について記述してきた。本発明の趣旨および範囲から逸脱することなく各種の変更形態がなされ得る。従って、他の実施形態は添付の請求項の範囲内に含まれる。
【図１】
0006487050.tif 000015
【図２】
0006487050.tif 000016
【図３】
0006487050.tif 000017