(57)【要約】      （修正有）【課題】擬似的距離測定値の選定による位置推定の修正方法を提供する。【解決手段】改善された位置Ｘｈｕｂｅｒが、ナビゲーション受信機とＮ基の衛星との間の測定距離に相当する、擬似的距離ρｉのＮ個の測定値（１０２）、及び前記受信機によってなされた前記受信機の位置の推定値Ｘｐｒｉｍ（１０３）を用いて、ロバスト推定アルゴリズムの適用（１０４）により決定される。測定値ρｉからの擬似的距離の剰余Δρｉｈｕｂｅｒよりのノルム剰余値Δｒｉｈｕｂｅｒを決定するステップと、Ｎ−ｋ個のノルム剰余値Δｒｉｈｕｂｅｒを含むシグマ（ギリシア文字σの大文字）個の部分集合（１０５、１０６）を決定するステップ（１１１）と、最小の標準偏差σＳＥＯを伴う部分集合ＳＥＯを選定するステップと、異常でない測定値を選定するステップ（１０７）と、選定された測定値からの位置の修正された推定値Ｘｃｌｎを決定するステップとを含む。【選択図】図１
2012073246.tif 000002
【特許請求の範囲】
【請求項１】
位置の推定を修正するための方法であって、
改善された位置Ｘｈｕｂｅｒが、ナビゲーション受信機とＮ基の衛星との間の測定距離に相当する、擬似的距離ρｉのＮ個の測定値（１０２）、及び前記受信機によってなされた前記受信機の位置の推定値Ｘｐｒｉｍ（１０３）を用いて、ロバスト推定アルゴリズムの適用（１０４）により決定され、少なくとも以下のステップ：
―前記測定値ρｉからの擬似的距離の剰余Δρｉｈｕｂｅｒよりのノルム剰余値Δｒｉｈｕｂｅｒの決定ステップと、
―ｋが１よりも厳密に大きい整数である、Ｎ−ｋ個のノルム剰余値Δｒｉｈｕｂｅｒを含むシグマ（ギリシア文字σの大文字）個の部分集合（１０５、１０６）の決定ステップ（１１１）と、
―最小の標準偏差σＳＥＯを伴う部分集合ＳＥＯの選定ステップと、
―前記ノルム剰余Δｒｉｈｕｂｅｒと前記部分集合ＳＥＯのノルム剰余の平均μＳＧＯとの間の差が、所定のしきい値Ｔ１よりも小さい場合に選定される、異常でない測定値の選定ステップ（１０７）と、
―前記選定された測定値からの、前記位置の修正された推定値Ｘｃｌｎの決定ステップとを含むことを特徴とする方法。
【請求項２】
ロバスト推定アルゴリズムがＨｕｂｅｒ（フーバー）アルゴリズムであることを特徴とする、請求項１に記載の方法。
【請求項３】
前記Ｈｕｂｅｒ（フーバー）アルゴリズムが、前記位置Ｘｐｒｉｍを用いて初期化されることを特徴とする、請求項２に記載の方法。
【請求項４】
前記剰余Δρｉｈｕｂｅｒが次の式：
Δρｉｈｕｂｅｒ＝δ（Ｘｉｓａｔ−Ｘｈｕｂｅｒ）
を用いることにより決定され、ここで、
Ｘｉｓａｔがｉ番目の衛星位置を表わし、前記位置が例えば信号伝達チャンネルを用いて前記受信機へと伝達され、
δ（Ｘｉｓａｔ−Ｘｈｕｂｅｒ）が前記差（Ｘｉｓａｔ−Ｘｈｕｂｅｒ）の推定値と前記測定された擬似的距離ρｉとの間の偏差を表わすことを特徴とする、請求項１〜３のいずれか一項に記載の方法。
【請求項５】
位置Ｘｃｌｎが最小二乗法又はロバスト推定法を用いることにより決定されることを特徴とする、請求項１〜４のいずれか一項に記載の方法。
【請求項６】
前記ノルム剰余Δｒｉｈｕｂｅｒが次の式：
Δｒｉｈｕｂｅｒ＝Δρｉｈｕｂｅｒ／σｉ 
を用いることにより決定され、ここで、
σｉが測定誤差の先験的分散を表わすことを特徴とする、請求項１〜５のいずれか一項に記載の方法。
【請求項７】
部分集合の数シグマ（ギリシア文字σの大文字）が、Ｎ個の内で（Ｎ−ｋ）個の測定値の組合せ数に相当することを特徴とする、請求項１〜６のいずれか一項に記載の方法。
【請求項８】
擬似的距離のＮ個の剰余Δｒｉｃｌｎを決定するためのステップを含み、次の式：
Δｒｉｃｌｎ＝δ（Ｘｉｓａｔ−Ｘｃｌｎ）／σｉ
が使用され、ここで、
δ（Ｘｉｓａｔ−Ｘｃｌｎ）が前記差（Ｘｉｓａｔ−Ｘｃｌｎ）の推定値と前記測定された擬似的距離ρｉとの間の偏差を表わすことを特徴とする、請求項１〜７のいずれか一項に記載の方法。
【請求項９】
既に選定された測定値から、測定値ρｉを選定するための第２のステップ（１０８）を含み、前記測定値は次の式：
｜Δｒｉｃｌｎ−Δｒｉｈｕｂｅｒ｜＜Ｔ２ 
が満足される場合に選定され、ここでＴ２が所定のしきい値であることを特徴とする、請求項７に記載の方法。
【請求項１０】
ロバスト位置推定Ｘｒｏｂが、第２の選定ステップ（１０８）の間に選定された測定値に対する前記最小二乗法の適用により推定されることを特徴とする、請求項８に記載の方法。
【請求項１１】
検出半径の値（１１０）が、前記第２の選定ステップ（１０８）の間に選定された前記測定値に基づいて決定されることを特徴とする、請求項９あるいは１０のいずれか一項に記載の方法。
【請求項１２】
請求項１〜１１のいずれか一項に記載の方法を実施することを特徴とする、ナビゲーション受信機。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、擬似的距離測定値の選定による位置推定の修正方法、及びその方法を実施する受信機に関する。それはとりわけ衛星ナビゲーションシステムの分野に適用される。
【背景技術】
【０００２】
衛星測位システムは通常、「全地球的航法衛星システム」（“Ｇｌｏｂａｌ  Ｎａｖｉｇａｔｉｏｎ  Ｓａｔｅｌｌｉｔｅ  Ｓｙｓｔｅｍ”）を表わす略語ＧＮＳＳによって示される。このタイプのシステムは、組み込まれたナビゲーション受信機を用いて端末の位置を推定できる。この受信機はそれに到達する複数の人工衛星から送信される信号に対する測定を行なう。これらの測定値は例えば擬似的距離に相当し、すなわち信号の送信時点における衛星と、信号の受信時点における受信機との間の距離の測定値に相当する。
【０００３】
ナビゲーション受信機により推定される位置は、必ずしも正確ではない。位置測定における精度の欠如を考慮に入れるため、推定される位置の周りに保護半径を定義することが一般的である。
【０００４】
様々なタイプの誤差が測定値の精度を悪化させ得る。これらの誤差は公称誤差及び非公称誤差の２つの範疇に分類され得る。
【０００５】
公称誤差は、システムの正常な動作において生じる妨害に起因する測定誤差である。一例として、衛星に組み込まれた送信機によって使用されるクロックは、完全には予測出来ない挙動を有する。推定の精度はまた、受信機の近傍における環境によっても影響され得る。公称誤差の別の形態は、受信機の近傍における衛星信号の地上又は建物での反射に起因する。また言及する価値があるのは、測定される信号の伝播へ遅れを導入し得る気象又は電離層現象である。
【０００６】
非公称誤差はシステムの機能不良の結果である。そのような誤差によって影響される測定値は、以下の記述において異常測定値とも呼ばれる。
【０００７】
位置推定の精度に対するこれらの測定誤差の影響を制限するため、誤差を取り除くためにそれらの識別を可能にし、利用可能な測定値に従って位置誤差に対する制限を計算出来るようにする手段が通常実施され、「受信機自律型システム完全性監視」“Ｒｅｃｅｉｖｅｒ  Ａｕｔｏｎｏｍｏｕｓ  Ｉｎｔｅｇｒｉｔｙ  Ｍｏｎｉｔｏｒｉｎｇ”を表わす頭字語ＲＡＩＭにより示される。
【０００８】
ＲＡＩＭの機能性を含む現行の位置測定機器は、多くの問題に悩まされている。
【０００９】
第１の問題は、これらの機器が完全に一体化されており、それはナビゲーション信号を取得する機器と、完全性の機能を確保する機器の位置を計算する機器とを別々に選定出来ないことを意味する。
【００１０】
第２の問題はこれらの機器が通常、最小二乗タイプのアルゴリズムに基づいていることである。これは、これらの機器が非公称誤差に由来する誤った測定値の存在により、及び／又はこれらの測定値に影響する誤差の振幅にかかわらず、誤差の不完全なモデル化により不安定にされることを意味する。従って提案された解は、そのような誤差に直面して信頼性がない。
【００１１】
Ｈｕｂｅｒ（フーバー）アルゴリズムの使用のような、ロバスト推定法は、異常測定値が存在する中で推定精度の改善を可能にするが、これらの誤差の影響は、それでもなお無視出来ない。
【発明の概要】
【発明が解決しようとする課題】
【００１２】
本発明の１つの目的は、特に上述の欠点を克服することである。
【課題を解決するための手段】
【００１３】
このため、本発明の主題は位置の推定を修正するための方法であり、改善された位置Ｘｈｕｂｅｒは、ナビゲーション受信機とＮ基の衛星との間の測定距離に相当する、擬似的距離ρｉのＮ個の測定値、及び前記受信機によってなされた前記受信機の位置の推定値Ｘｐｒｉｍを用いて、ロバスト推定アルゴリズムの適用により決定される。本方法は少なくとも以下のステップ：
―測定値ρｉからの擬似的距離の剰余Δρｉｈｕｂｅｒよりのノルム剰余値Δｒｉｈｕｂｅｒの決定ステップと、
―ｋが１よりも厳密に大きい整数である、Ｎ−ｋ個のノルム剰余値Δｒｉｈｕｂｅｒを含むシグマ（ギリシア文字σの大文字）個の部分集合の決定ステップと、
―最小の標準偏差σＳＥＯを伴う部分集合ＳＥＯの選定ステップと、
―ノルム剰余Δｒｉｈｕｂｅｒと部分集合ＳＥＯのノルム剰余の平均μＳＧＯとの間の差が、所定のしきい値Ｔ１よりも小さい場合に選定される、異常でない測定値の選定ステップと、
―選定された測定値からの、位置の修正された推定値Ｘｃｌｎの決定ステップと
を含む。
【００１４】
本発明の１つの態様によれば、ロバスト推定アルゴリズムはＨｕｂｅｒ（フーバー）アルゴリズムである。
【００１５】
Ｈｕｂｅｒ（フーバー）アルゴリズムは、例えば位置Ｘｐｒｉｍを用いることにより初期化される。
【００１６】
１つの実施形態において、剰余Δρｉｈｕｂｅｒは次の式を用いて決定される。
Δρｉｈｕｂｅｒ＝δ（Ｘｉｓａｔ−Ｘｈｕｂｅｒ）
ここで、
Ｘｉｓａｔはｉ番目の衛星位置を表わし、前記位置は例えば信号伝達チャンネルを用いて受信機へと伝達され、
δ（Ｘｉｓａｔ−Ｘｈｕｂｅｒ）は差（Ｘｉｓａｔ−Ｘｈｕｂｅｒ）の推定値と測定された擬似的距離ρｉとの間の偏差を表わす。
【００１７】
本発明の有利な実施形態によれば、位置Ｘｃｌｎは最小二乗法又はロバスト推定法を用いて決定される。
【００１８】
ノルム剰余値Δｒｉｈｕｂｅｒは、次の式を用いて決定され得る。
Δｒｉｈｕｂｅｒ＝Δρｉｈｕｂｅｒ／σｉ
ここで、σｉは測定誤差の先験的分散を表わす。
【００１９】
本発明の１つの態様によれば、部分集合の数シグマ（ギリシア文字σの大文字）はＮ個の内で（Ｎ−ｋ）個の測定値の組合せ数に相当する。
【００２０】
さらに、本発明は擬似的距離のＮ個の剰余Δｒｉｃｌｎを決定するためのステップを提供し、次の式が用いられる。
Δｒｉｃｌｎ＝δ（Ｘｉｓａｔ−Ｘｃｌｎ）／σｉ
ここで、
δ（Ｘｉｓａｔ−Ｘｃｌｎ）は差（Ｘｉｓａｔ−Ｘｃｌｎ）の推定値と測定された擬似的距離ρｉとの間の偏差を表わす。
【００２１】
本方法は例えば、既に選定された測定値から、測定値ρｉを選定するための第二のステップを含み、前記測定値は次の式が満足される場合に選定される。
｜Δｒｉｃｌｎ−Δｒｉｈｕｂｅｒ｜＜Ｔ２
ここで、Ｔ２は所定のしきい値である。
【００２２】
１つの実施モードにおいて、ロバスト位置推定値Ｘｒｏｂは第二の選定ステップの間に選定された測定値に対し、最小二乗法を適用することにより推定される。
【００２３】
検出半径の値は、例えば第二の選定ステップの間に選定された測定値に基づき決定される。
【００２４】
本発明はまた前述の方法を実施するナビゲーション受信機にも関する。
【００２５】
有利なことに、本方法はＲＡＩＭ処理を行なうハードウェアとは無関係に、ナビゲーション信号を取得するためのハードウェアの選定を最適化できる。さらに、ロバストＲＡＩＭ法は、標準のＲＡＩＭと比較して検出に関する完全性の性能を改善しながら、位置推定の信頼性を高めることができる。
【００２６】
本発明の別の特徴及び利点は、限定されない例証として与えられる以下の記述から、及び添付図を踏まえて明らかになるであろう。
【図面の簡単な説明】
【００２７】
【図１】擬似的距離を選定することにより、ナビゲーション受信機の位置を推定する方法の単純化された図を示す。
【発明を実施するための形態】
【００２８】
本方法は、一次受信機１００、すなわち前記方法を実施しない従来のナビゲーション受信機によって送信されるデータに基づく。これらのデータは一方で一次受信機により推定される解１０１に相当し、他方でＮ基の衛星から生じる信号を処理することにより得られる、擬似的距離の測定値１０２に相当する。一次受信機により推定される解は、Ｘｐｒｉｍで表わされる受信機の推定位置、及び推定されるクロックのずれ（時間差）に相当する。ρｉで表わされる擬似的距離のＮ個の測定値１０２は、受信機とＮ基の衛星との間の測定距離に相当し、ｉはそれに対して測定が行なわれる、Ｎ基の衛星内の１つの衛星の添え字に相当する。
【００２９】
一次受信機により送信された測定値が前処理されていない場合、それらは四角の破線１０３により記号で表わされているように、本来知られている、伝播誤差及び測定誤差を取り除く前処理を受けるべきである。
【００３０】
前述の解１０１及び測定値１０２に適用される最初の処理１０４は、ロバスト推定法、すなわち非公称誤差の存在下で効果的な方法の使用により、受信機の改善された位置の決定を目的とする。本発明の範囲内で採用され得るロバスト法は、いわゆるＨｕｂｅｒ（フーバー）法である。この方法はＸ．Ｗ．Ｃｈａｎｇ及びＹ．Ｇｕｏによる題名、Ｈｕｂｅｒ’ｓ  Ｍ−ｅｓｔｉｍａｔｉｏｎ  ｉｎ  ｒｅｌａｔｉｖｅ  ＧＰＳ  ｐｏｓｉｔｉｏｎｉｎｇ（フーバー法−相対的ＧＰＳ測位における推定）：ｃｏｍｐｕｔａｔｉｏｎａｌ  ａｓｐｅｃｔｓ，  Ｊｏｕｒｎａｌ  ｏｆ  Ｇｅｏｄｅｓｙ，  ２００５，  ｖｏｌ．７９，  ｎｏ．６−７，  ｐｐ．３５１−３６２の論文において説明されている。この方法の原理は最新の計算位置に関する測定値の剰余の関数により、測定値ρｉの重み付けをすることである。この計算は、例えば一次受信機により決定される位置Ｘｐｒｉｍを用いることによって、初期化され得る。この結果はＸｈｕｂｅｒと称される改善された位置、及び擬似的距離の剰余の集合である。その剰余はΔρｉｈｕｂｅｒと称され、例えば次の式：
Δρｉｈｕｂｅｒ＝δ（Ｘｉｓａｔ−Ｘｈｕｂｅｒ）                      （１）
を用いることにより決定される。
ここで、
Ｘｉｓａｔはｉ番目の衛星位置を表わし、前記位置は例えば信号伝達チャンネルを用いて受信機へと伝達され、
δ（Ｘｉｓａｔ−Ｘｈｕｂｅｒ）は差（Ｘｉｓａｔ−Ｘｈｕｂｅｒ）の推定値と測定された擬似的距離ρｉとの間の偏差を表わす。
【００３１】
このようにして得られる、改善された位置Ｘｈｕｂｅｒ及び擬似的距離の剰余Δρｉｈｕｂｅｒの集合は、次にＮ個の測定値ρｉの中で任意の異常測定値を検出するために使用される。
【００３２】
異常測定値の検出は、まず初めにノルム剰余値Δｒｉｈｕｂｅｒを定義し、次に前記剰余の部分集合を形成し、そして測定値が異常か否かを決定するために、これらの部分集合に対して統計的テストを適用することによって行なわれる。
【００３３】
ノルム剰余値Δｒｉｈｕｂｅｒは、次の式を用いることにより決定され得る。
Δｒｉｈｕｂｅｒ＝Δρｉｈｕｂｅｒ／σｉ                            （２）
ここで、σｉは測定誤差の先験的分散を表わし、ガウスの定理に基づく誤差分布モデルが通常使用される。
【００３４】
Ｎ−ｋ個のノルム剰余値Δｒｉｈｕｂｅｒを含む部分集合１０５、１０６は、次に１１１に形成され、パラメータｋはそれに対して防護されるべき、異常測定値の最大数に相当する。
【００３５】
このようにして決定される部分集合の数シグマ（ギリシア文字σの大文字）は，Ｎ個の内で（Ｎ−ｋ）個の測定値の組合せ数に相当し、シグマ（ギリシア文字σの大文字）＝Ｃ［（Ｎ−ｋ），Ｎ］で表わされる。言い換えれば、シグマ（ギリシア文字σの大文字）個の部分集合はＮ−ｋ個の剰余を含む全ての組合せである。
【００３６】
これらの集合の標準偏差が次に計算される。これらの集合のうち、最適部分集合と呼ばれ頭字語ＳＧＯにより表わされる集合が決定される。ＳＧＯはその標準偏差σＳＥＯが最小の集合である。ＳＥＯに属するノルム剰余の平均が決定され、μＳＧＯで表わされる。集合ＳＥＯの決定は、誤差が存在する最大の確率を有するｋ個の測定値をリストアップすること、及び信頼出来ると見なされる基準の部分集合を持つことを可能にする。
【００３７】
２つの統計的テスト１０７、１０８は、次に測定値の異常性を否定又は確認するために適用される。
【００３８】
第１の統計的テスト１０７は、ノルム剰余Δｒｉｈｕｂｅｒと平均μＳＧＯとの間の差により定義されるＮ個の偏差を、しきい値Ｔ１と比較するためのものである。このしきい値Ｔ１は、誤りの警報及び検出の性能に関して、設計者により選定された目的を保証するように決定される。このテストは例えば、次の式を用いて公式化され得る。
Δｒｉｈｕｂｅｒ−μＳＥＯ＜Ｔ１                                  （３）
添え字ｉに関する不等式（３）が満足される場合、ｉ番目の測定ρｉは異常とは考えられない。この場合、それは適正位置と呼ばれる位置Ｘｃｌｎを再計算するために保持される。この適正位置は最小二乗法を用いて決定されることが望ましいが、例えばロバスト推定法のような別の推定方法が用いられ得る。
【００３９】
位置修正の改善を可能にする、第２の統計的テスト１０８が次に適用される。このために、ノルム擬似的距離のＮ個の剰余Δｒｉｃｌｎの新たな集合が形成される。この集合の剰余は、例えば次の式を用いて定義される。
Δｒｉｃｌｎ＝δ（Ｘｉｓａｔ−Ｘｃｌｎ）／σｉ                      （４）
ここで、
δ（Ｘｉｓａｔ−Ｘｃｌｎ）は差（Ｘｉｓａｔ−Ｘｃｌｎ）の推定値と測定された擬似的距離ρｉとの間の偏差を表わす。
【００４０】
この集合の標準偏差σｃｌｎが決定される。剰余Δｒｉｃｌｎと剰余Δｒｉｈｕｂｅｒとの間の偏差は、新たに推定された位置Ｘｃｌｎに対する第１のテスト１０７から生じる、測定値のフィルタリングの貢献を定量化するために決定される。
この偏差がしきい値Ｔ２よりも小さい場合、測定値は保持される。しきい値Ｔ２は求められる誤りの警報及び検出の性能レベルに従って選定される。言い換えれば、このテストは次の式を用いることにより要約され得る。
｜Δｒｉｃｌｎ−Δｒｉｈｕｂｅｒ｜＜Ｔ２                            （５）
不等式（５）が満足される場合、測定値ｉは有効で、すなわち異常でないと見なされ、そして保持される。
２つのテスト適用に続いて保持される測定値は、いわゆるロバスト位置Ｘｒｏｂの決定１０９のために使用される。このために、従来の最小二乗法が用いられ得る。別のタイプの推定もまた使用され得る。
【００４１】
保護半径１１０の推定は次に、保持された測定値の数を考慮し、既存の方法を用いることにより行なわれ得る。一例として、ＬＳ−ＲＡＩＭ（Ｌｅａｓｔ−Ｓｑｕａｒｅ  ＲＡＩＭ＝最小二乗ＲＡＩＭ）タイプのＲＡＩＭ法、及びＭＨＳＳ（Ｍｕｌｔｉｐｌｅ  Ｈｙｐｏｔｈｅｓｅｓ  Ｓｅｐａｒａｔｅ  Ｓｏｌｕｔｉｏｎｓ＝多重仮説分離解）の解分離法タイプの範囲内で用いられる、半径推定アルゴリズムが使用され得る。
【００４２】
従って、前述のステップの適用は、入力の測定値ρｉにおける誤差が検出された場合に、修正された位置の解を得ることを可能にする。位置の解Ｘｒｏｂの保証を可能にする保護半径の値もまた利用できる。
【符号の説明】
【００４３】
１００  一次受信機
１０１  一次受信機により推定される解
１０２  擬似的距離の測定値
１０３  前処理
１０４  ロバスト推定アルゴリズムの適用
１０５  部分集合
１０６  部分集合
１０７  第１の統計的テスト
１０８  第２の統計的テスト
１０９  ロバスト位置Ｘｒｏｂの決定
１１０  検出半径の値
１１１  シグマ（ギリシア文字σの大文字）個の部分集合を決定するステップ
【図１】
2012073246.tif 000003
【外国語明細書】
2012073246000001.pdf