モバイルデバイスの位置の決定方法
【課題】衛星信号から週番号を抽出する必要なく、モバイルデバイスの位置を決定する。
【解決手段】前記システム時刻に依存する前記衛星の軌道を仮定して、衛星から受信する信号からモバイルデバイスの位置を規定する一連のデバイスパラメータを決定する方法であって、複数の週番号候補の１つを有するシステム時刻によって、それぞれ実行され、正当な週番号を複数の週番号候補の中から選択し、決定された一連のデバイスパラメータとして指定される正当な週番号に基づいて前記システムの式を解き、ここで、正当な週番号の選択は、解法アルゴリズムのアプリケーションが解を生じなかった少なくともいくつかの週番号を削除することによって、実行され、残りのいずれかの週番号に偏差値を割り当てて、偏差値に従って、残りの週番号の中から正当な週番号を選択する、ことを特徴とする方法。
【解決手段】前記システム時刻に依存する前記衛星の軌道を仮定して、衛星から受信する信号からモバイルデバイスの位置を規定する一連のデバイスパラメータを決定する方法であって、複数の週番号候補の１つを有するシステム時刻によって、それぞれ実行され、正当な週番号を複数の週番号候補の中から選択し、決定された一連のデバイスパラメータとして指定される正当な週番号に基づいて前記システムの式を解き、ここで、正当な週番号の選択は、解法アルゴリズムのアプリケーションが解を生じなかった少なくともいくつかの週番号を削除することによって、実行され、残りのいずれかの週番号に偏差値を割り当てて、偏差値に従って、残りの週番号の中から正当な週番号を選択する、ことを特徴とする方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明に係る方法は、全地球的航法衛星システム（ＧＮＳＳ）の衛星からの信号に基づいて、モバイルデバイスの位置を決定することを可能とする。具体的には、グローバル・ポジショニング・システム（ＧＰＳ）の衛星からの信号に基づいて、モバイルデバイスの位置を決定することを可能とする。
【背景技術】
【０００２】
従来の方法では、週の時刻（time of week）と同様に、モバイルデバイスの位置を決定するために既知でなければならない週番号（week number）は、連続するサブフレームで構成される連続するフレームとして体系化されるデータを含む衛星信号から抽出される。フレームを構成する５つのサブフレームごとにエンコードされる週の時刻とは異なり、週番号は、第１のサブフレームにのみ含まれるので、完全なフレームの伝送に必要な時間に、約３０秒掛かる可能性がある。正当な衛星軌道データによるコールドスタート（cold start）が、週の時刻を受信するまで待機しなければならないために、初期位置算出時間（first fix）が同様な期間まで増加する可能性がある。
【発明の概要】
【０００３】
本発明の目的は、週番号の情報を先に有することなく初期位置を得ることにより、事前に衛星信号から週番号を抽出する必要なく、モバイルデバイスの位置を決定する方法を提供することである。しかしながら、衛星軌道がモバイルデバイスに記憶される軌道データ（orbit data）から引き出すことができるコールドスタートの場合は、６秒ごとに伝送される週の時刻のみを、衛星信号から抽出するので、初期測位時間を、最大３０秒、平均で約１２秒程度にまで削減する。
【０００４】
以下、添付図面を参照して、発明をより詳細に説明する。
【図面の簡単な説明】
【０００５】
【図１】ＧＰＳ衛星のいずれかから伝送される一連のデータの構造を概略的に示す図である。
【発明を実施するための形態】
【０００６】
グローバル・ポジショニング・システム（ＧＰＳ）に関連する衛星はいずれも、送信に１ミリ秒かかる１０２３ビットの特有の２進数列の反復を送信する。データビットは、送信される特有の数列の連続する２０個の複製に付加される。データビットは、反転する場合は、データビットの−１を示し、反転しない場合は、データビットの＋１値を示す。したがって、データビットの伝送は、２０ｍｓかかり、データのビットレートは、毎秒５０ビットである。
【０００７】
全てのＧＰＳ衛星の軌道を時間の関数として規定する暦データ（almanac data）と、暦データなどよりも高い精度で非常に短い時間間隔で衛星の軌道を規定する伝送衛星の天体暦データ（ephemeris data）とを有する完全な一連のデータを衛星から伝送するためには、１２．５秒掛かる。一連のデータは、それぞれ３０秒の２５個のフレームからなる。全てのフレームは、３００のビットのワードで構成され、それぞれ６秒かかる５個のサブフレームからなる。全てのサブフレームの第２のワードは、ハンドオーバワード（hand-over word）であり、最初の１７ビットは、週の時刻ＴＯＷ、すなわち続くサブフレームの伝送開始時間である。第３のワードの初めに、それぞれのフレームの第１のフレームのみに、週番号ＷＮが１０ビットの番号、すなわちＧＰＳモジュール１０２４の開始からの週番号を有する。
【０００８】
位置が決定されるモバイルデバイス、すなわち受信機は、通常は少なくとも４つ、一般的には４つより多い視野内の全ての衛星の信号を受信し、その信号から必要なデータを抽出することになるであろう。これと同時にデバイスは、発振器により測定するデバイス自身の時間に応答して、衛星信号のコード位相（code phases）を決定することになるであろう。モバイルデバイスの位置は、抽出されたデータと、コード位相とから抽出される衛星の位置から、先行技術により信号のミリ秒の整数アンビギュイティを解くことにより決定する。これにより、「測定された」擬似距離を規定した後に、解法アルゴリズムを使用するモバイルデバイスの位置の３つの要素と、時計バイアスとを含むデバイスパラメータのセットのための測定された擬似距離を含むシステムの式を解く。
【０００９】
衛星の位置は、暦データ、軌道データ、又は好適には精度がより良好な長期データ（long term data）、又は衛星の軌道を規定する天体暦データから抽出できる。そして、衛星の位置は、週の時刻ＴＯＷと、週番号ＷＮを含むＧＰＳシステム時間が識別されると、モバイルデバイスに記憶される。伝送時のそれぞれの衛星の位置は、ＧＰＳ時間とともに、衛星の軌道を時間の関数として規定する軌道データを評価することによって、抽出できる。
【００１０】
１つの週の時間シフトは、衛星速度が約３０００ｍ／秒で、ほぼ５０００ｋｍのＧＰＳ衛星の位置におけるシフトの作用において、２７分３２秒の時間シフトに本質的に等しい。この結果として、週番号が誤っている場合は、モバイルデバイス、すなわち受信機の決定した位置に大きなエラーが起きる。多くの場合、解法アルゴリズムは、全く解法を生み出さない。
【００１１】
時間バイアスが未知であると仮定すると、解法は、過度に決定される（over-determined）ことになるであろう。ここで、５つ以上の衛星が視野に入る。その結果として、擬似距離ρSをパラメータのセットに関連させる基本的な式から開始して、最小二乗法、又は類似する方法によって、一般的には規定される。すなわち、ベクトルｘR＝（ｘ0R，ｘ1R，ｘ2R，ｘ3R）は、通常、時計バイアスｘ0Rと、ユークリッド座標システム（Euclidean coordinate system）における信号を受信する時点でのデバイスの空間位置の構成要素ｘ1R、ｘ2R、ｘ3Rと有する。
【数１】
ここで、
【数２】
は、信号を放射するときの衛星ｓと、信号を受信するときのモバイルデバイスとの間の距離である。空間位置の構成要素のように、時間バイアスは、メートルである。すなわち、秒の時間における光の速度ｍ／ｓである。衛星の時間バイアスは、無視してよいと考えられる。
【００１２】
コード位相の測定から決定される衛星の測定された擬似距離ＲSとともに、微分項
【数３】
を、変数と考えられるデバイスパラメータｘのセットの関数として規定できる。
【００１３】
（１）及び（２）から
【数４】
と、
【数５】
が分かる。
【００１４】
ここで当業者に周知の方法によって、重み関数を形成できる。
【数６】
この式は、特定の重みを一連のデバイスパラメータｘの値に割り当てる。ここでデバイスパラメータは、微分項ｄＳ（ｘ）、すなわち測定された擬似距離利ＲSと、先の値から計算される擬似距離ρS（ｘ）との間の差異に依存する。解ｘRは、一連のデバイスパラメータｘとして規定される。ここで、Ｗ（ｘ）は、最小である。ＰS1S2は重み行列である。ここで、Ｓ1とＳ2とは、視界内のｍ個の衛星によって変化する。すなわち、擬似距離の測定の共分散行列の逆数である。
【００１５】
通常は、測定間に相互相関はなく、重み行列は、対角であると仮定される。この場合、
【数７】
が、分かる。ここで、σＳは、衛星ヘルスデータ（satellite health data）、高度、及び他のパラメータに依存するか、又は定数であるかのいずれかである衛星ｓの微分項ｄS（ｘ）の標準偏差である。次いで、（６）は、
【数８】
に変換される。
【００１６】
このＷ（ｘ）は最小値であり、すなわちｘRにおいて、
【数９】
を示唆する。
【００１７】
第１の中間解として一連のデバイスパラメータｘ0の見積もりから開始し、新たな中間解を第１の相互作用ステップにおいて、抽出できる。このために、パラメータの組における擬似距離の依存性は、線形化される。
【数１０】
ここで、ρS0＝ρS（ｘ0）、
【数１１】
、及びΔｉ＝ｘｉ−ｘｉ0である。
【００１８】
ｄS0＝ＲS−ρS0とともに、これから
【数１２】
が導かれる。そして（８）と（１０）とから
【数１３】
を生じるとともに、
【数１４】
を、補正Δｊについて解くことができる
【数１５】
と表すことができる。
【００１９】
擬似距離測定間の相互相関がない場合は、（６´）を使用して、式（１３）を
【数１６】
に変換できる。
【００２０】
いずれの場合においても、新たな中間解は、
【数１７】
によって規定される。次いでこれは、次の相関ステップにおいて先の中間解として使用される。このように、式（８）のシステムは、インタラクティブに解くことができる。ここで、相関ステップにおいて、先の中間解の全ての構成要素は、先の中間解の対応する構成要素と、対応する補正の和に置き換えられる。
【数１８】
【００２１】
この解法アルゴリズムが収束する場合は、最新の新たな中間解は、式（８）のシステムの解ｘRと考えられる。
【００２２】
週番号が事前に既知でないときは、この式のシステムの解の解明は、複数の週番号候補を用いて、試みられる。この場合、通常は連続する週番号のシーケンスを形成することになるであろう。すなわち、見積もった週番号に集中する、又は見積もった週番号から広がる間隔は埋められることになるであろう。
【００２３】
ともにＧＰＳ時間を規定する週番号候補、及び週の時刻のいずれかにより、相互作用は、ｎmax＝１０など、先に規定されるステップの最大数ｎmaxまで行うことができる。すなわち、最大ノルム（maximum norm）、又はユークリッドのノルムなどの先に規定されたいくつかのノルムに従って、先の中間解のパラメータベクトルｘnと、これから抽出される新たな中間解のパラメータベクトルｘn+1との間の差異が所定のしきい値
【数１９】
よりも小さくなるまで行われる。条件（１７）が実現されることなくステップの最大数に達する場合は、解法アルゴリズムは、解を生じず、週番号候補は削除されることになる。
【００２４】
解法アルゴリズムは、複数の週番号のそれぞれについて解ｘRを生じている場合は、単一の適当な週番号を選択しなければならない。これは、一連のデバイスパラメータから計算される擬似距離と、測定された擬似距離との差異を反映する偏差の関数を使用して、偏差値を抽出するために解ｘRを評価して、実行できる。そして、偏差値Ｄが、最小となる週番号を選択する。好適な偏差の関数は、偏差値が（６）、又は（６´）から
【数２０】
として抽出される重み関数Ｗ（ｘ）である。
【００２５】
可能な他の選択は、微分項ｄS（ｘR）の平均値ｄである。この双方を考慮できることは当然である。
【００２６】
いずれの場合においても、チェックの結果によって、正当な週番号を確認する妥当性のチェック、又は正当な週番号を確認しない妥当性のチェックをさらに実行することは、望ましい。例えば、正当な週番号と関連付けられる最小の偏差値は、しきい値と比較して、偏差値がしきい値よりも小さい場合に、偏差値に基づく正当な週番号、及び解であると確認できる。また、最小の偏差値を次に小さい偏差値と比較して、前者を後者で除した余りがしきい値よりも小さい場合のみ、週番号であることを確認できる。２番目に小さい偏差値の代わりに、残りの偏差値の平均を使用できる。しきい値は、固定でき、又は視界内の衛星数、及び信号強度などのパラメータによることができる。
【００２７】
モバイルデバイスの位置を先験的に限定できない場合は、解による位置は、妥当に見積もられた位置に特に先験的に十分に近接した規定する妥当な範囲内にあるか否かをチェックして、見積もられた位置からの距離が特定のしきい値を下回る場合にのみ、週番号、及び解が正当であると確認することが可能である。いくつかのチェックを組み合わせることができることは当然である。
【００２８】
正当な週番号が確認された場合に、デバイスパラメータの対応する組は、表示など、デバイスで使用する。正当な週番号が確認されなかった場合は、対応する解のデバイスパラメータは破棄される。次いで、週番号候補の異なる組による新たな検索を実行できる。全ての試みが失敗した場合は結局、週番号は、衛星の信号から抽出しなければならないだろう。
【００２９】
一連のデバイスパラメータの決定のシミュレーションを行った。ここで、モバイルデバイスは、経度８．５６度、緯度４７．２８度、及び海抜５５０ｍに位置した。位置の最初の見積もり、すなわち第１の中間解は、経度０度、緯度０度、及び海抜０ｍとした。比較的強い８つの衛星信号があると仮定した。デバイスに記憶される長期の信号を使用して、衛星の位置を計算した。週番号候補は、ＧＰＳ時間の開始から週１５２５、１５２６、及び１５２７にそれぞれ対応する５０１、５０２、及び５０３であった。
【００３０】
以下の表において、最初の５つの相互作用ステップによって到達した補正は、週番号候補５０１のリストである。
【表１】
【００３１】
５つの相互作用ステップの後に、空間座標は、収束するが、時間バイアスｘ0の補正Δ0は、まだ大きく、ステップに基づく週番号、及び解の補正Δ0が正当であるゼロに近い値に近づいていくことが示されない。解法アルゴリズムは、一見して収束しない。したがって、週番号５０１は、削除される。
【００３２】
週番号の候補５０２による補正は、同様に示される。
【表２】
【００３３】
解法アルゴリズムは、収束せず、週番号候補５０２もまた、削除される。
【００３４】
一方、週番号の候補５０２による解法アルゴリズムは、３つの相互作用ステップの後に正当な解に到達して、相互作用は停止する。
【表３】
【００３５】
解は、正当であると考えられ、５０３は、正当な週番号として選択される。
【００３６】
モバイルデバイスの位置は、上述の例と同一であると仮定した第２のシミュレーションでは、以下の結果を得た。
【００３７】
〔ＧＰＳ時間の開始から週１５２６に対応する週番号の候補５０２〕
１０個の相互作用ステップの後に収束しない。週番号は削除される。
【００３８】
〔ＧＰＳ時間の開始から週１５２７に対応する週番号の候補５０３〕
７つのステップの後に収束。（６´）に従う重み関数により計算した偏差値Ｄは、２５４´８０８．６であった。微分項ｄS（ｘR）の平均ｄは、４８２´１４８．２である。
【００３９】
〔ＧＰＳ時間の開始から週１５２７に対応する週番号の候補５０３〕
４つの相互作用ステップの後に収束。Ｄ＝０．１２、ｄ＝８４１．６である。週番号５０４は、残りの週番号５０３と５０４とから正当な週番号として選択された。
【符号の説明】
【００４０】
ＴＯＷ 週の時刻
ＷＮ 週番号
【技術分野】
【０００１】
本発明に係る方法は、全地球的航法衛星システム（ＧＮＳＳ）の衛星からの信号に基づいて、モバイルデバイスの位置を決定することを可能とする。具体的には、グローバル・ポジショニング・システム（ＧＰＳ）の衛星からの信号に基づいて、モバイルデバイスの位置を決定することを可能とする。
【背景技術】
【０００２】
従来の方法では、週の時刻（time of week）と同様に、モバイルデバイスの位置を決定するために既知でなければならない週番号（week number）は、連続するサブフレームで構成される連続するフレームとして体系化されるデータを含む衛星信号から抽出される。フレームを構成する５つのサブフレームごとにエンコードされる週の時刻とは異なり、週番号は、第１のサブフレームにのみ含まれるので、完全なフレームの伝送に必要な時間に、約３０秒掛かる可能性がある。正当な衛星軌道データによるコールドスタート（cold start）が、週の時刻を受信するまで待機しなければならないために、初期位置算出時間（first fix）が同様な期間まで増加する可能性がある。
【発明の概要】
【０００３】
本発明の目的は、週番号の情報を先に有することなく初期位置を得ることにより、事前に衛星信号から週番号を抽出する必要なく、モバイルデバイスの位置を決定する方法を提供することである。しかしながら、衛星軌道がモバイルデバイスに記憶される軌道データ（orbit data）から引き出すことができるコールドスタートの場合は、６秒ごとに伝送される週の時刻のみを、衛星信号から抽出するので、初期測位時間を、最大３０秒、平均で約１２秒程度にまで削減する。
【０００４】
以下、添付図面を参照して、発明をより詳細に説明する。
【図面の簡単な説明】
【０００５】
【図１】ＧＰＳ衛星のいずれかから伝送される一連のデータの構造を概略的に示す図である。
【発明を実施するための形態】
【０００６】
グローバル・ポジショニング・システム（ＧＰＳ）に関連する衛星はいずれも、送信に１ミリ秒かかる１０２３ビットの特有の２進数列の反復を送信する。データビットは、送信される特有の数列の連続する２０個の複製に付加される。データビットは、反転する場合は、データビットの−１を示し、反転しない場合は、データビットの＋１値を示す。したがって、データビットの伝送は、２０ｍｓかかり、データのビットレートは、毎秒５０ビットである。
【０００７】
全てのＧＰＳ衛星の軌道を時間の関数として規定する暦データ（almanac data）と、暦データなどよりも高い精度で非常に短い時間間隔で衛星の軌道を規定する伝送衛星の天体暦データ（ephemeris data）とを有する完全な一連のデータを衛星から伝送するためには、１２．５秒掛かる。一連のデータは、それぞれ３０秒の２５個のフレームからなる。全てのフレームは、３００のビットのワードで構成され、それぞれ６秒かかる５個のサブフレームからなる。全てのサブフレームの第２のワードは、ハンドオーバワード（hand-over word）であり、最初の１７ビットは、週の時刻ＴＯＷ、すなわち続くサブフレームの伝送開始時間である。第３のワードの初めに、それぞれのフレームの第１のフレームのみに、週番号ＷＮが１０ビットの番号、すなわちＧＰＳモジュール１０２４の開始からの週番号を有する。
【０００８】
位置が決定されるモバイルデバイス、すなわち受信機は、通常は少なくとも４つ、一般的には４つより多い視野内の全ての衛星の信号を受信し、その信号から必要なデータを抽出することになるであろう。これと同時にデバイスは、発振器により測定するデバイス自身の時間に応答して、衛星信号のコード位相（code phases）を決定することになるであろう。モバイルデバイスの位置は、抽出されたデータと、コード位相とから抽出される衛星の位置から、先行技術により信号のミリ秒の整数アンビギュイティを解くことにより決定する。これにより、「測定された」擬似距離を規定した後に、解法アルゴリズムを使用するモバイルデバイスの位置の３つの要素と、時計バイアスとを含むデバイスパラメータのセットのための測定された擬似距離を含むシステムの式を解く。
【０００９】
衛星の位置は、暦データ、軌道データ、又は好適には精度がより良好な長期データ（long term data）、又は衛星の軌道を規定する天体暦データから抽出できる。そして、衛星の位置は、週の時刻ＴＯＷと、週番号ＷＮを含むＧＰＳシステム時間が識別されると、モバイルデバイスに記憶される。伝送時のそれぞれの衛星の位置は、ＧＰＳ時間とともに、衛星の軌道を時間の関数として規定する軌道データを評価することによって、抽出できる。
【００１０】
１つの週の時間シフトは、衛星速度が約３０００ｍ／秒で、ほぼ５０００ｋｍのＧＰＳ衛星の位置におけるシフトの作用において、２７分３２秒の時間シフトに本質的に等しい。この結果として、週番号が誤っている場合は、モバイルデバイス、すなわち受信機の決定した位置に大きなエラーが起きる。多くの場合、解法アルゴリズムは、全く解法を生み出さない。
【００１１】
時間バイアスが未知であると仮定すると、解法は、過度に決定される（over-determined）ことになるであろう。ここで、５つ以上の衛星が視野に入る。その結果として、擬似距離ρSをパラメータのセットに関連させる基本的な式から開始して、最小二乗法、又は類似する方法によって、一般的には規定される。すなわち、ベクトルｘR＝（ｘ0R，ｘ1R，ｘ2R，ｘ3R）は、通常、時計バイアスｘ0Rと、ユークリッド座標システム（Euclidean coordinate system）における信号を受信する時点でのデバイスの空間位置の構成要素ｘ1R、ｘ2R、ｘ3Rと有する。
【数１】
ここで、
【数２】
は、信号を放射するときの衛星ｓと、信号を受信するときのモバイルデバイスとの間の距離である。空間位置の構成要素のように、時間バイアスは、メートルである。すなわち、秒の時間における光の速度ｍ／ｓである。衛星の時間バイアスは、無視してよいと考えられる。
【００１２】
コード位相の測定から決定される衛星の測定された擬似距離ＲSとともに、微分項
【数３】
を、変数と考えられるデバイスパラメータｘのセットの関数として規定できる。
【００１３】
（１）及び（２）から
【数４】
と、
【数５】
が分かる。
【００１４】
ここで当業者に周知の方法によって、重み関数を形成できる。
【数６】
この式は、特定の重みを一連のデバイスパラメータｘの値に割り当てる。ここでデバイスパラメータは、微分項ｄＳ（ｘ）、すなわち測定された擬似距離利ＲSと、先の値から計算される擬似距離ρS（ｘ）との間の差異に依存する。解ｘRは、一連のデバイスパラメータｘとして規定される。ここで、Ｗ（ｘ）は、最小である。ＰS1S2は重み行列である。ここで、Ｓ1とＳ2とは、視界内のｍ個の衛星によって変化する。すなわち、擬似距離の測定の共分散行列の逆数である。
【００１５】
通常は、測定間に相互相関はなく、重み行列は、対角であると仮定される。この場合、
【数７】
が、分かる。ここで、σＳは、衛星ヘルスデータ（satellite health data）、高度、及び他のパラメータに依存するか、又は定数であるかのいずれかである衛星ｓの微分項ｄS（ｘ）の標準偏差である。次いで、（６）は、
【数８】
に変換される。
【００１６】
このＷ（ｘ）は最小値であり、すなわちｘRにおいて、
【数９】
を示唆する。
【００１７】
第１の中間解として一連のデバイスパラメータｘ0の見積もりから開始し、新たな中間解を第１の相互作用ステップにおいて、抽出できる。このために、パラメータの組における擬似距離の依存性は、線形化される。
【数１０】
ここで、ρS0＝ρS（ｘ0）、
【数１１】
、及びΔｉ＝ｘｉ−ｘｉ0である。
【００１８】
ｄS0＝ＲS−ρS0とともに、これから
【数１２】
が導かれる。そして（８）と（１０）とから
【数１３】
を生じるとともに、
【数１４】
を、補正Δｊについて解くことができる
【数１５】
と表すことができる。
【００１９】
擬似距離測定間の相互相関がない場合は、（６´）を使用して、式（１３）を
【数１６】
に変換できる。
【００２０】
いずれの場合においても、新たな中間解は、
【数１７】
によって規定される。次いでこれは、次の相関ステップにおいて先の中間解として使用される。このように、式（８）のシステムは、インタラクティブに解くことができる。ここで、相関ステップにおいて、先の中間解の全ての構成要素は、先の中間解の対応する構成要素と、対応する補正の和に置き換えられる。
【数１８】
【００２１】
この解法アルゴリズムが収束する場合は、最新の新たな中間解は、式（８）のシステムの解ｘRと考えられる。
【００２２】
週番号が事前に既知でないときは、この式のシステムの解の解明は、複数の週番号候補を用いて、試みられる。この場合、通常は連続する週番号のシーケンスを形成することになるであろう。すなわち、見積もった週番号に集中する、又は見積もった週番号から広がる間隔は埋められることになるであろう。
【００２３】
ともにＧＰＳ時間を規定する週番号候補、及び週の時刻のいずれかにより、相互作用は、ｎmax＝１０など、先に規定されるステップの最大数ｎmaxまで行うことができる。すなわち、最大ノルム（maximum norm）、又はユークリッドのノルムなどの先に規定されたいくつかのノルムに従って、先の中間解のパラメータベクトルｘnと、これから抽出される新たな中間解のパラメータベクトルｘn+1との間の差異が所定のしきい値
【数１９】
よりも小さくなるまで行われる。条件（１７）が実現されることなくステップの最大数に達する場合は、解法アルゴリズムは、解を生じず、週番号候補は削除されることになる。
【００２４】
解法アルゴリズムは、複数の週番号のそれぞれについて解ｘRを生じている場合は、単一の適当な週番号を選択しなければならない。これは、一連のデバイスパラメータから計算される擬似距離と、測定された擬似距離との差異を反映する偏差の関数を使用して、偏差値を抽出するために解ｘRを評価して、実行できる。そして、偏差値Ｄが、最小となる週番号を選択する。好適な偏差の関数は、偏差値が（６）、又は（６´）から
【数２０】
として抽出される重み関数Ｗ（ｘ）である。
【００２５】
可能な他の選択は、微分項ｄS（ｘR）の平均値ｄである。この双方を考慮できることは当然である。
【００２６】
いずれの場合においても、チェックの結果によって、正当な週番号を確認する妥当性のチェック、又は正当な週番号を確認しない妥当性のチェックをさらに実行することは、望ましい。例えば、正当な週番号と関連付けられる最小の偏差値は、しきい値と比較して、偏差値がしきい値よりも小さい場合に、偏差値に基づく正当な週番号、及び解であると確認できる。また、最小の偏差値を次に小さい偏差値と比較して、前者を後者で除した余りがしきい値よりも小さい場合のみ、週番号であることを確認できる。２番目に小さい偏差値の代わりに、残りの偏差値の平均を使用できる。しきい値は、固定でき、又は視界内の衛星数、及び信号強度などのパラメータによることができる。
【００２７】
モバイルデバイスの位置を先験的に限定できない場合は、解による位置は、妥当に見積もられた位置に特に先験的に十分に近接した規定する妥当な範囲内にあるか否かをチェックして、見積もられた位置からの距離が特定のしきい値を下回る場合にのみ、週番号、及び解が正当であると確認することが可能である。いくつかのチェックを組み合わせることができることは当然である。
【００２８】
正当な週番号が確認された場合に、デバイスパラメータの対応する組は、表示など、デバイスで使用する。正当な週番号が確認されなかった場合は、対応する解のデバイスパラメータは破棄される。次いで、週番号候補の異なる組による新たな検索を実行できる。全ての試みが失敗した場合は結局、週番号は、衛星の信号から抽出しなければならないだろう。
【００２９】
一連のデバイスパラメータの決定のシミュレーションを行った。ここで、モバイルデバイスは、経度８．５６度、緯度４７．２８度、及び海抜５５０ｍに位置した。位置の最初の見積もり、すなわち第１の中間解は、経度０度、緯度０度、及び海抜０ｍとした。比較的強い８つの衛星信号があると仮定した。デバイスに記憶される長期の信号を使用して、衛星の位置を計算した。週番号候補は、ＧＰＳ時間の開始から週１５２５、１５２６、及び１５２７にそれぞれ対応する５０１、５０２、及び５０３であった。
【００３０】
以下の表において、最初の５つの相互作用ステップによって到達した補正は、週番号候補５０１のリストである。
【表１】
【００３１】
５つの相互作用ステップの後に、空間座標は、収束するが、時間バイアスｘ0の補正Δ0は、まだ大きく、ステップに基づく週番号、及び解の補正Δ0が正当であるゼロに近い値に近づいていくことが示されない。解法アルゴリズムは、一見して収束しない。したがって、週番号５０１は、削除される。
【００３２】
週番号の候補５０２による補正は、同様に示される。
【表２】
【００３３】
解法アルゴリズムは、収束せず、週番号候補５０２もまた、削除される。
【００３４】
一方、週番号の候補５０２による解法アルゴリズムは、３つの相互作用ステップの後に正当な解に到達して、相互作用は停止する。
【表３】
【００３５】
解は、正当であると考えられ、５０３は、正当な週番号として選択される。
【００３６】
モバイルデバイスの位置は、上述の例と同一であると仮定した第２のシミュレーションでは、以下の結果を得た。
【００３７】
〔ＧＰＳ時間の開始から週１５２６に対応する週番号の候補５０２〕
１０個の相互作用ステップの後に収束しない。週番号は削除される。
【００３８】
〔ＧＰＳ時間の開始から週１５２７に対応する週番号の候補５０３〕
７つのステップの後に収束。（６´）に従う重み関数により計算した偏差値Ｄは、２５４´８０８．６であった。微分項ｄS（ｘR）の平均ｄは、４８２´１４８．２である。
【００３９】
〔ＧＰＳ時間の開始から週１５２７に対応する週番号の候補５０３〕
４つの相互作用ステップの後に収束。Ｄ＝０．１２、ｄ＝８４１．６である。週番号５０４は、残りの週番号５０３と５０４とから正当な週番号として選択された。
【符号の説明】
【００４０】
ＴＯＷ 週の時刻
ＷＮ 週番号
【特許請求の範囲】
【請求項１】
全地球的航法衛星システムに関連する衛星から受信する信号のコード位相を算定し、モバイルデバイスの位置は、週番号（ＷＮ）と、週の時刻（ＴＯＷ）とを含む全地球的航法衛星システムのシステム時刻と前記コード位相とから抽出されるとともに、解法アルゴリズムを使用するシステムの式を解くことによって、前記システム時刻に依存する前記衛星の軌道を仮定する、前記全地球的航法衛星システムに関連する前記衛星から受信する前記信号から前記モバイルデバイスの位置を規定する一連のデバイスパラメータを決定する方法であって、
複数の週番号候補の１つを有するシステム時刻によって、それぞれ実行され、正当な週番号を前記複数の週番号候補の中から選択し、前記決定された一連のデバイスパラメータとして指定される前記正当な週番号に基づいて前記システムの式を解き、
ここで、前記正当な週番号の選択は、前記解法アルゴリズムのアプリケーションが解を生じなかった少なくともいくつかの週番号を削除することによって、実行され、残りのいずれかの週番号に偏差値を割り当てて、前記偏差値に従って、前記残りの週番号の中から正当な週番号を選択する、ことを特徴とする方法。
【請求項２】
前記偏差値の絶対値が最も小さい前記週番号候補を、前記正当な週番号として選択する請求項１に記載の方法。
【請求項３】
前記システムの式は、前記衛星のそれぞれ１つと、一連のデバイスパラメータにおける前記モバイルデバイスとの間の擬似距離への依存を規定する基本式に基づく請求項１、又は請求項２に記載の方法。
【請求項４】
前記システムの式の前記解法は、いずれの場合でも、一連のデバイスパラメータであって、前記一連のデバイスパラメータに依存する偏差の関数の絶対値は、最小値を有し、前記偏差値は、前記一連のデバイスパラメータにおける前記偏差の関数の値を反映する一連のデバイスパラメータである請求項１〜３のいずれか一項に記載の方法。
【請求項５】
前記偏差の関数は、前記算定されたコード位相から抽出されたときの前記衛星の１つ、及び前記モバイルデバイスの間の擬似距離の値と、前記モバイルデバイスの前記一連のデバイスパラメータから前記基本式を介して抽出される前記擬似距離の値との差にそれぞれ対応する微分項を含む請求項３、又は請求項４に記載の方法。
【請求項６】
前記偏差の関数は、前記微分項の加重二乗の和に比例する請求項５に記載の方法。
【請求項７】
前記解法アルゴリズムは、前記偏差の関数の勾配がゼロである位置を解くことによる前記偏差の関数の極値の探索を含む請求項４〜６のいずれか一項に記載の方法。
【請求項８】
前記解法アルゴリズムは、先の中間解を補正して、前記補正を包含する新しい中間解に先の中間解を置き換えられるために、システムの式の線形バージョンを解く少なくとも１つの相互作用ステップを有する請求項７に記載の方法。
【請求項９】
前記週番号候補は、連続する週番号の間隔を埋める請求項１〜８のいずれか一項に記載の方法。
【請求項１０】
前記正当な週番号は、少なくとも１つのさらなるチェックを受け、対応する解は、前記少なくとも１つのチェックが前記正当な週番号を確認した場合にのみ、使用される請求項１〜９のいずれか一項に記載の方法。
【請求項１１】
前記少なくとも１つのチェックは、
前記対応する偏差値をしきい値と比較するステップ、
前記対応する偏差値と、残りの週番号候補に対応する少なくとも１つの偏差値から抽出される項との間の関係と、しきい値とを比較するステップ、及び
前記対応する解による前記デバイスの前記位置は、先験的に規定された妥当な範囲内にあるか否かをチェックするステップ
の１つ、又は２つ以上を含む請求項１０に記載の方法。
【請求項１】
全地球的航法衛星システムに関連する衛星から受信する信号のコード位相を算定し、モバイルデバイスの位置は、週番号（ＷＮ）と、週の時刻（ＴＯＷ）とを含む全地球的航法衛星システムのシステム時刻と前記コード位相とから抽出されるとともに、解法アルゴリズムを使用するシステムの式を解くことによって、前記システム時刻に依存する前記衛星の軌道を仮定する、前記全地球的航法衛星システムに関連する前記衛星から受信する前記信号から前記モバイルデバイスの位置を規定する一連のデバイスパラメータを決定する方法であって、
複数の週番号候補の１つを有するシステム時刻によって、それぞれ実行され、正当な週番号を前記複数の週番号候補の中から選択し、前記決定された一連のデバイスパラメータとして指定される前記正当な週番号に基づいて前記システムの式を解き、
ここで、前記正当な週番号の選択は、前記解法アルゴリズムのアプリケーションが解を生じなかった少なくともいくつかの週番号を削除することによって、実行され、残りのいずれかの週番号に偏差値を割り当てて、前記偏差値に従って、前記残りの週番号の中から正当な週番号を選択する、ことを特徴とする方法。
【請求項２】
前記偏差値の絶対値が最も小さい前記週番号候補を、前記正当な週番号として選択する請求項１に記載の方法。
【請求項３】
前記システムの式は、前記衛星のそれぞれ１つと、一連のデバイスパラメータにおける前記モバイルデバイスとの間の擬似距離への依存を規定する基本式に基づく請求項１、又は請求項２に記載の方法。
【請求項４】
前記システムの式の前記解法は、いずれの場合でも、一連のデバイスパラメータであって、前記一連のデバイスパラメータに依存する偏差の関数の絶対値は、最小値を有し、前記偏差値は、前記一連のデバイスパラメータにおける前記偏差の関数の値を反映する一連のデバイスパラメータである請求項１〜３のいずれか一項に記載の方法。
【請求項５】
前記偏差の関数は、前記算定されたコード位相から抽出されたときの前記衛星の１つ、及び前記モバイルデバイスの間の擬似距離の値と、前記モバイルデバイスの前記一連のデバイスパラメータから前記基本式を介して抽出される前記擬似距離の値との差にそれぞれ対応する微分項を含む請求項３、又は請求項４に記載の方法。
【請求項６】
前記偏差の関数は、前記微分項の加重二乗の和に比例する請求項５に記載の方法。
【請求項７】
前記解法アルゴリズムは、前記偏差の関数の勾配がゼロである位置を解くことによる前記偏差の関数の極値の探索を含む請求項４〜６のいずれか一項に記載の方法。
【請求項８】
前記解法アルゴリズムは、先の中間解を補正して、前記補正を包含する新しい中間解に先の中間解を置き換えられるために、システムの式の線形バージョンを解く少なくとも１つの相互作用ステップを有する請求項７に記載の方法。
【請求項９】
前記週番号候補は、連続する週番号の間隔を埋める請求項１〜８のいずれか一項に記載の方法。
【請求項１０】
前記正当な週番号は、少なくとも１つのさらなるチェックを受け、対応する解は、前記少なくとも１つのチェックが前記正当な週番号を確認した場合にのみ、使用される請求項１〜９のいずれか一項に記載の方法。
【請求項１１】
前記少なくとも１つのチェックは、
前記対応する偏差値をしきい値と比較するステップ、
前記対応する偏差値と、残りの週番号候補に対応する少なくとも１つの偏差値から抽出される項との間の関係と、しきい値とを比較するステップ、及び
前記対応する解による前記デバイスの前記位置は、先験的に規定された妥当な範囲内にあるか否かをチェックするステップ
の１つ、又は２つ以上を含む請求項１０に記載の方法。
【図１】
【公開番号】特開２０１０−２７６６０５（Ｐ２０１０−２７６６０５Ａ）
【公開日】平成２２年１２月９日（２０１０．１２．９）
【国際特許分類】
【出願番号】特願２０１０−１２００７４（Ｐ２０１０−１２００７４）
【出願日】平成２２年５月２６日（２０１０．５．２６）
【出願人】（５１０１４５９６７）ウー−ブロックス アクチェンゲゼルシャフト (2)
【Ｆターム（参考）】
【公開日】平成２２年１２月９日（２０１０．１２．９）
【国際特許分類】
【出願日】平成２２年５月２６日（２０１０．５．２６）
【出願人】（５１０１４５９６７）ウー−ブロックス アクチェンゲゼルシャフト (2)
【Ｆターム（参考）】
[ Back to top ]