JP4896746B2 - 衛星を利用した移動端末の測位 - Google Patents

衛星を利用した移動端末の測位 Download PDF

Info

Publication number
JP4896746B2
JP4896746B2 JP2006554055A JP2006554055A JP4896746B2 JP 4896746 B2 JP4896746 B2 JP 4896746B2 JP 2006554055 A JP2006554055 A JP 2006554055A JP 2006554055 A JP2006554055 A JP 2006554055A JP 4896746 B2 JP4896746 B2 JP 4896746B2
Authority
JP
Japan
Prior art keywords
mobile terminal
parameter
pseudorange
quality
measure
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 - Fee Related
Application number
JP2006554055A
Other languages
English (en)
Other versions
JP2007523345A (ja
Inventor
アリ カンガス，
Original Assignee
テレフオンアクチーボラゲット エル エム エリクソン（パブル）
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 by テレフオンアクチーボラゲット エル エム エリクソン（パブル） filed Critical テレフオンアクチーボラゲット エル エム エリクソン（パブル）
Publication of JP2007523345A publication Critical patent/JP2007523345A/ja
Application granted granted Critical
Publication of JP4896746B2 publication Critical patent/JP4896746B2/ja
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/40—Correcting position, velocity or attitude

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)
  • Navigation (AREA)
  • Radio Relay Systems (AREA)

Description

本発明は、一般に衛星を利用した測位に関するものであり、詳細には、衛星からの信号が受信可能な、地上の通信ネットワークの移動端末の測位に関するものである。
近年、物体、装置、あるいは装置を携行する人物の地理的位置を決定することが、多くの応用分野において関心を集めている。測位を解決する一つの方法は、衛星から出る信号を使って位置を決定することである。そのようなシステムの例として良く知られているのは、世界測位衛星システム（ＧＰＳ）［１］（非特許文献１）と、来るべきガリレオシステムである。衛星から受信した複数の信号に基づく三角測量／三辺測量として特定される座標系に関して位置が与えられる。
ＧＰＳ受信機を携帯通信システムの移動端末（ユーザ装置、移動局、移動ノード等と呼ぶこともある）に円滑に統合するため、アシステッドＧＰＳ（ＡＧＰＳ）［２］（非特許文献２）が開発された。概してアシステッドＧＰＳは、検出感度、位置の推定値を取得するまでに要する時間時間、精度、バッテリの節約など、さまざまな点でＧＰＳ受信機の性能を高めることを目的としている。これは、移動局の中のＧＰＳ受信機から一部の機能をネットワークに移し、ＧＰＳ受信機自身の中ではＧＰＳタスクのサブセットだけを実行することによって、行われる。
スタンドアロンのＧＰＳ受信機は、公称搬送周波数と、信号によって搬送されるデータが変調される規則を除くと、システムについてのそれ以外の情報を持たなくてもＧＰＳ衛星信号への完全な同期を行うことができる。そのような受信機は複数の衛星（通常は４個）によって送信された測距信号を測定する。信号は、各衛星について一意であり１ミリ秒毎に繰り返して現れるいわゆるＣ／Ａコード(Coarse/Acquisition code)を含む。Ｃ／Ａコード上には、２０ミリ秒のビット周期を持つ航法データのビットストリーム(navigation data bit stream)が重畳されている。航法データは、送信時の衛星位置の計算を可能にするパラメータと、衛星の時計のオフセットを記述するパラメータとを含む。スタンドアロンのＧＰＳ受信機は、通常、受信機の位置の計算より前に航法データのストリーム全体を復号する必要がある。これには相当長い時間がかかることがあり、一定の最小限の信号強度が必要である。受信機は、航法メッセージの復号に要する信号強度よりはるかに低い信号強度でＣ／Ａコードの境界を決定することができる。
ＧＰＳを備えた移動端末の測距を容易にするため、ＡＧＰＳでは、より高速でより信頼できる通信リンク、例えば基地局と移動端末間の無線通信リンク上で、航法データは代わりに支援データとして送信される。支援データは、通常、概算のＧＰＳシステム時刻と、移動端末の概算の位置も含む。（操作モードによっては、移動端末は代わりに、Ｃ／Ａコードの境界の決定を高速化できる一組のパラメータを受信することがある。）
ＡＧＰＳには、移動局ベースのＡＧＰＳ(mobile station based AGPS)と移動局アシストＡＧＰＳ(mobile station assisted AGPS)の２種類がある。移動局ベースのＡＧＰＳにおいては、移動局の位置は、移動局が決定する測距信号測定結果およびネットワークが提供する支援データを用いて移動局の中で計算される。移動局アシストＡＧＰＳ（ネットワークベースのＡＧＰＳと呼ばれることもある）においては、移動局は、受信した測距信号の測定を行って、測定結果をネットワーク内の位置サーバに報告する（移動局は位置サーバから支援データの受信も行う）。移動局がどこに位置しているかに関する報告された測定結果および事前に存在する情報に基づいて、位置サーバが移動局の位置を計算する。
ＧＰＳを備えた移動端末の位置を移動局アシストＡＧＰＳを用いて計算するための良く知られた先行技術による方法は、一般に、衛星への（不完全な）疑似距離と信号受信時の移動端末クロック測定値（ｔｕ）とを含む、報告された測定結果に依存する。しかし、ｔｕが全然報告されない場合や誤って報告される場合もあり得る。ＡＧＰＳについての現在最も一般的な標準に準拠している移動端末はｔｕを報告することが義務付けられているが、将来の標準が他の解決法の選択を許容することは十分あり得る。
別の選択肢として、各衛星からの送信時刻（ｔｔｉ）が複数の段階で決められるような解決法もある。まず、ｔｔｉのミリ秒未満の部分が、各衛星についてのＣ／Ａコードの境界を見つけることによって決定される。その際には、可能性のあるコード位相とドップラー偏移をすべて試験する相関器が用いられる。その後、ｔｔｉのミリ秒部分が再構成される。これには普通、生の航法データビットを残して、受信データが逆拡散される必要がある。その後、ＴＯＷ（Ｔｉｍｅ ｏｆ ｗｅｅｋ：１週間における時刻）の再構成が、複数の選択肢の内の一つ、例えば「ＴＯＷの直接復調」、「相関技法又はリアルタイムクロックによるＴＯＷの再構成」等を用いて行われる。どの方法が使われたとしても、移動を支援する移動局アシストＡＧＰＳは、その後伝搬遅延を補償する必要があり、従って、測定時の概算のＧＰＳシステム時刻を報告する必要がある。補償は、例えば、支援データとして与えられる一部の情報要素に基づいて行うことができる。
測定されたＧＰＳネットワーク時刻を再構成するための要件は、測位機能が、最良の条件でもおよそ−１７９ｄＢＷまで低下することを示唆する。もう１つの欠点は、この種の再構成には時間がかかることである。必要な航法データのビットが受信されるまでにたいてい８秒かかる。
米国特許［３］（特許文献１）は、移動端末がＧＰＳネットワーク時刻を知らない状況を扱う手法について説明している。衛星測定の概算の時刻がタイムスタンプとして記される。本当のＧＰＳネットワーク時刻とこの測定時刻との差は、その後の移動端末位置の計算において、未知であるとして扱われる。このため、もう一度衛星測定を追加することが必要となる。すなわち、従来の測位方法に用いられる４度ではなく衛星測定が５度必要となる。
衛星を利用した移動端末の測位についての先行技術による解決策は一般に、性能や精度の限界を伴い、あるいは（又は加えて）複雑な計算を必要とする。測位メカニズムの改善と、特に、測位における適切なエラー訂正（誤差の訂正）が必要とされている。
米国特許第６，４３０，４１５号明細書 Navstar GPS Space Segment/Navigation user Interfaces, ICD-GPS-200, Revision IRN-200C-003, １９９９年１０月１１日 3GPP TS 44. 031, v 5.6.0
本発明の一般的な目的は、衛星信号が受信できる移動端末の測位を改善することである。具体的な目的は、ＡＧＰＳに適したエラー訂正メカニズムを提供することである。もう一つの目的は、移動端末で測定されたパラメータに関する推定または再構成に起因するエラー（誤差）の訂正が可能な、衛星を利用した移動端末の測位を提供することである。さらにもう一つの目的は、合理的な計算複雑性と両立する効率的なエラー訂正が可能な、衛星を利用した移動端末の測位を実現することである。
これらの目的は添付されている請求項に従って実現される。
簡単に述べると、本発明は移動端末の位置を特定する（測位する）ための手順の改善をもたらすものであり、移動端末は、複数の衛星から信号を受信して衛星への疑似距離を測定する。本手順は、移動端末で測定された時刻及び疑似距離のうちの少なくとも一方に関するパラメータのエラー（またはそのようなパラメータの表現においけるエラー）の訂正を含む。主な狙いは、必要なときだけエラー訂正を行い、それによって移動端末の位置を計算する処理の複雑性を最小化する一方で、決定された位置の十分な正確さは実現することである。移動端末の位置は、まず、エラー訂正無しで、または部分的なエラー訂正を伴って、疑似距離に基づいて計算される。その後、位置計算に関連する品質の尺度が決定され、好ましくは移動端末の位置の最小二乗解の残差に基づいて、このパラメータの推定値の品質が所定の品質基準と比較される。所定の品質基準は、例えば、最小損失関数のような適切な品質パラメータにおける所定の閾値であっても良い。所定の品質基準が満たされた場合には、それ以上の計算は必要ない。他方、最初の計算の品質が不十分な場合には、移動端末の位置は、移動端末で測定されその位置の計算に使われた少なくとも一つの時刻／疑似距離に関するパラメータ（またはそのようなパラメータの表現／再構成）におけるエラーの訂正を追加して再計算される。
本発明に従って移動端末の位置を計算（または再計算）するとき、ＵＥの位置およびＵＥのクロックのバイアスについての未知パラメータに関する一組の非線形方程式が定義される。もし時刻パラメータのエラー訂正が含まれていたら、（例えば各衛星からの信号送信時刻と衛星を利用した測位システムのネットワーク時刻との差として定義された）時刻のオフセットについての未知パラメータも追加される。非線形連立方程式は線形化され、好ましくは加重最小二乗の方法により従来どおり解かれる。
もし測定された疑似距離のうち少なくとも一つが、切り捨てられた端数の疑似距離である場合、計算手順は一般に、移動端末での信号受信時刻のパラメータと概算の移動端末位置とに基づいて、切り捨てられた端数の疑似距離から完全な疑似距離を再構成することを含む。
本発明で提案する方法は、信号受信時刻のパラメータにおけるエラーの訂正、または少なくとも１つの疑似距離のパラメータにおけるエラーの訂正、あるいはその両方の組み合わせを含むことがある。特に有利な実施形態において、追加的なエラー訂正は、位置計算に関連する２つ以上の品質の尺度についての評価に基づいて、必要だと考えられる程度まで段階的に導入される。段階的なエラー訂正を使った実施形態に従って、測位はまず、移動端末によって報告された信号受信時刻のパラメータが利用可能な場合には、移動端末によって報告された信号受信時刻に基づいて実行される。結果の品質が第１の基準に合致しない場合、測定された受信時刻の中のエラーを訂正するため、例えば送信時刻のオフセットの推定を含む新規の計算が行われる。品質は第２の基準を用いて再度チェックされ、品質が不十分な場合には、今度は疑似距離のエラーの訂正を伴って、また、場合によっては送信時刻のオフセットの推定と組み合わせて、位置が再度再計算される。
本発明は、エラー訂正機能の強化を通じて、移動端末の効率的な測位を可能にし、その結果、一層正確な位置計算が可能になる。高度な解決策は、必要な時だけエラー訂正を行うことを可能にする品質チェックを備えている。このようにして、複雑で時間のかかる計算が回避できることが多くなり、全体のシステム性能も向上する。
本発明の別の態様により、移動端末、測位ノード、及び、移動端末の位置を衛星を利用して決定する手段を備える通信システムが提供される。
本発明およびそのさらなる目的や優位点は、下記の説明および添付した図面を参照することによって最も良く理解できよう。
このセクションの最後に、略語集を掲載している。
以下の詳細な説明において、ＧＰＳシステムに実装された実施形態を説明する。しかし、対応する原理が、ＧＬＯＮＡＳＳあるいはヨーロッパのＧａｌｉｌｅｏ衛星航法システムのような、いかなる衛星を利用した測位システムにも応用できることが、当業者には理解できる。
同様に、以下の詳細説明の中では、ＷＣＤＭＡシステムがモデルシステムとして使われる。しかし、本発明は他の無線通信システムにも応用可能である。本発明が応用可能な他のシステムの排他的でない例は、例えば、ＣＤＭＡ−２０００システム、ＧＳＭシステム、ＧＳＭ／ＧＰＲＳシステム等である。他の無線通信システムに適用される時には、異なる機能の実装が、当該システムの異なる端末やノードにおいて行われることになる。
本明細書において「移動端末」という用語は、無線通信システム内で移動可能なあらゆる種類の端末を意味するのに用いられる。排他的でない例として、セルラ電話、携帯情報端末、ポータブルコンピュータなどがある。
図１は、本発明を使用可能な通信システムの一例を示す概略図である。このシナリオ例において、基本的な無線通信システム１１０がＧＰＳシステム１３０と一緒に使用されて、移動局アシストＡＧＰＳを形成している。
図１に例示する無線通信システム１１０は、無線アクセスネットワーク（ＲＡＮ）（例えばＵＴＲＡＮ）およびコアネットワーク１１８を含むＷＣＤＭＡシステムである。ＲＡＮは、無線に関する機能を実行し、ユーザ装置／移動端末１１２とそれ以外のネットワークとの間の接続の確立を担う。ＲＡＮは、通常、ノードＢとも呼ばれる非常に多数の基地局（ＢＴＳ）１１４と無線ネットワーク制御装置（ＲＮＣ）１１６とを含む。各々のＢＴＳは、それぞれの通信範囲内で移動端末にサービス提供し、いくつかのＢＴＳがＲＮＣによって制御される。ＲＮＣの主な機能は、周波数、拡散符号またはスクランブリング符号、及び、チャネルの電力レベルを割り当てることである。ＲＮＣ１１６は、コアネットワーク１１８へのアクセスを提供する。コアネットワークは、例えば、交換局、ＧＳＭ／ＧＰＲＳコアネットワークにおいて交換局に相当する支援ノードおよびデータベースを含み、また一般的にはマルチメディア処理装置を含む。コアネットワークは外部ネットワーク１２０（例えばインターネット、ＰＳＴＮ，ＩＳＤＮ、及び、他のＰＬＭＮ）と通信する。
実際には、多くのＷＣＤＭＡネットワーク、およびほとんどの他の無線通信ネットワークには、図１の基本的な例よりはるかに複雑な方法で配置された複数のネットワーク要素およびノードがある。
図１において、ＧＰＳシステム１３０が、宇宙船（ＧＰＳ衛星，ＳＶ）、すなわち、通常では衛星１３２によって表されている。各々のＳＶはそれぞれの測距信号を送信し、そのような多数の測距信号を用いて、移動端末１１２または移動端末を持ち運んでいる人の位置が決定できる。
例えば無線通信システム１１０のコアネットワーク１１８において測位要求が発生すると、測位要求はＲＮＣ１１６に関連する測位ノード１１７に提供される。ＲＮＣ１１６は移動端末１１２へと送信される制御信号を通じて、衛星測距信号の測定を命令する。測定命令は、支援データを伴う。支援データは衛星の位置データと衛星の時刻基準データを含み、例えば通信システム１１０に接続している参照受信機（不図示）経由でＲＮＣ１１６によって引き出されてもよい。そのような参照受信機は１台のユニットとして提供されても良いし、複数の部分に分割されて衛星の時間基準データと衛星の位置データの決定を分離しても良い。支援データは通常、測位ノード１１７の中で処理される。測位ノードは、例えば、どの衛星が測距信号を検出可能な位置に存在するかを、移動端末１１２に送信される前に決定できる。
移動端末１１２は、衛星の測距信号を検出できる受信機を備えており、端末１１２は支援データを用いて衛星の測距信号への同期と測距信号の測定を円滑に行う。次いで、測定された測距信号は、移動端末１１２の位置を計算するのに用いられる。「移動局ベースのＡＧＰＳ」が用いられる場合、測距信号の処理は移動端末の中で行われる。「移動局アシストＡＧＰＳ」が用いられる場合、測距信号またはそれを表現する情報は、代わりにネットワークノード（通常は例示されたシナリオにおいては測位ノード）へと送信される。ネットワークノードでは、移動端末から報告された測定結果と、移動端末がどこに位置するかに関する事前の情報とに基づいて、測位計算が行われる。図１の例において、測位機能はＲＮＣ１１６に接続した外部の測位ノード１１７の中で実現される。測位ノードは例えばスタンドアロンＡＧＰＳ ＳＭＬＣであっても良い。あるいは、測位機能がＲＮＣ１１６の中または別のネットワークノードの中に組み込まれることも考えられる。
現在のＧＰＳ ＳＶは、１５７５．４２ ＭＨｚを中心にした帯域を使って測距信号を送信する。信号は、各ＳＶに一意であって、通常は１０２３チップの長さと１／１．０２３ｘ１０６秒のチップ持続時間を持つＣ／Ａコードを含んでいる。Ｃ／Ａコードは１ミリ秒毎に繰り返す。Ｃ／Ａコード上には、２０ミリ秒のビット周期を持つ航法データのビットストリームが重畳されている。とりわけ、航法データは、信号送信時の衛星の正確な位置の計算を可能にする一連のエフェメリスパラメータを含む。航法データは、多くのサブフレームに分割され、各々は長さが６秒で、各サブフレームは１０ワードに分割される。タイムスタンプ、すなわちＧＰＳのＴＯＷ、がすべてのサブフレームの２番目のワード（ＨＯＷ）において送信される。示された時刻が、問題のサブフレームの最後に送信された時刻である。従って、ＴＯＷは６秒毎に繰り返される。
各測距信号は基本的に、移動局によって測定されたクロックを定義する。クロックは信号送信の時刻を示す。移動局がＧＰＳシステム時刻を知っている場合、クロックの表示は、時間的遅延、すなわち測距信号を送信中のＳＶから移動局までの距離を決定するのに直接利用できる。３つの距離を測定し、送信時のＧＰＳ衛星の位置についての知識を使うことによって、移動局の位置が三次元で決定される。しかし、通常は、移動局は正確なＧＰＳシステム時刻についての知識を持っていないため、移動局のクロックのバイアスをなくすためもう１つの測定が必要である。
図２は、図１で説明したシステムのさまざまな部分についてクロックの関係を（ミリ秒単位で表して）図解している。ＳＶは正確な原子時計を持っていて、クロックの安定性を維持している。しかし、ＳＶによる送信は、タイミング図２００に図解するように、完全にＧＰＳシステム時刻に同期しているわけではない。図２において、シーケンス２０１は、ＧＰＳシステム時刻を表し、シーケンス２０２はＳＶ１のクロックを表し、シーケンス２０３はＳＶＮ（ＧＰＳ衛星Ｎ）のクロックを表し、シーケンス２０６は移動端末／ＵＥのクロックを表し、シーケンス２０４と２０５はそれぞれ移動端末が受信したＳＶ１乃至ＳＶＮからの測距信号の中で読まれたとおりの時刻を表す。タイミング図２００を通る垂直線を引くことによって、宇宙のさまざまなポイントで観察されたクロック表示すべてのスナップショットが得られる。ＧＰＳシステム時刻２０１は、地上局のクロックの集合とＳＶのクロック（ＳＶクロック）の部分集合とに基づく集合平均として定義される。個々のＳＶクロック２０２と２０３およびＵＥクロック２０６は、ＧＰＳシステム時刻に比べてわずかにオフセットがある。個々のオフセットのモデルは、各ＳＶから航法メッセージの一部として送信される。信号が地表面上の移動端末／ＵＥに到達するとき、信号は遅延されている。遅延は、問題の衛星からＵＥまでの距離に依存するが、クロックの表示２０４、２０５によって示されているように、通常は６０−８５ミリ秒である。
ＧＰＳ受信機は、複数の衛星への疑似距離を測定する。（ただし、注意すべきことだが、移動端末に組み込まれているＧＰＳ受信機は疑似距離を決めるのに必要な計算のすべてを行う必要はなく、疑似距離を計算／再構成するのに必要な基本的なデータを提供するだけのことが多い）。疑似距離は、次式で表すことができる。

ρｉ＝ｃ・（ｔｕ−ｔｔｉ） （１）

ここで、ｔｕは信号受信時のＵＥクロックの表示（図２の２０６）であり、ｔｔｉはｉ番目のＳＶの信号送信時の時刻、ｃは波の伝搬速度である。疑似距離は多数の不安定要因、例えば受信機クロックのバイアス、電離層遅延および対流圏遅延、ＳＶクロックのバイアス、測定エラー（誤差）等により、本当の距離とは異なることがある。本明細書においては、説明を簡明にするために、これらのエラー源の大半の影響は無視されている。上記に挙げたエラー源（例えば、［１］を参照のこと。）の多くを補償するため、既知の先行技術がある。さらに、ＳＶの動きや地球の自転の影響も、これらの影響を扱う方法は当技術分野でよく知られているため、省略する（例えば［１］を参照のこと）。こうして単純化されたモデルでは測定された疑似距離は以下に従う。

ρｉ＝｜ｘｕ（ｔｒ）−ｘｓｉ（ｔｔｉ）｜＋ｂ＋ｅｉ （２）

ここで、ｘｕ＝（ｘｕ ｙｕ ｚｕ）は、ＧＰＳシステム時刻ｔｒでの未知の受信機の位置に関する三次元座標を含む行ベクトルである。同様に、ｘｓｉはｉ番目のＳＶの送信時刻ｔｔｉでの座標を含む行ベクトルである。｜ｚ｜という表記は、かっこ内のベクトル量のノルムを意味する。これは、（ｚｚＴ）１／２に等しい。この場合、これは受信機（移動端末）とＳＶとの間の距離として解釈できる。ｂは受信機クロックのバイアス（距離として表現される）であり、

ｂ＝ｃ・（ｔｕ−ｔＧＰＳ） （３）

と表現される。ここで、ｔＧＰＳは、ＧＰＳシステム時刻である。また、ｅｉは測定エラーである。
ＳＶは通常、３．８４ｋｍ／秒の速度で動いており、従って位置の精度を劣化させるのでない限り、送信時刻はミリ秒レベルで知る必要がある。送信時刻ｔｔｉは通常、複数の段階で決定される。第一に、ｔｔｉのミリ秒未満の部分は、各ＳＶについてＣ／Ａコードの境界を見つけることによって決定される（図２を参照のこと）。これは、可能性のあるコード位相およびドップラー偏移をすべて試験する相関器を用いて行われる。その後のステップにおいて、送信時刻のミリ秒部分が再構成される必要がある。通常、受信データは逆拡散され、生の航法データビットが残される。その後、ＴＯＷの再構成を、複数の技術の選択肢、例えば「ＴＯＷの直接復調」、「相関技法又はリアルタイムクロックによる再構成」等によって行うことができる。どの方法が使われようと、移動を支援する移動局アシストＡＧＰＳは伝搬遅延を補償する必要があり、従って、測定時の概算のＧＰＳシステム時刻を報告する必要がある。補償は例えば支援データとして提供される情報に基づいて行うことができる。
背景技術のセクションで述べたとおり、測定されたＧＰＳネットワーク時刻を再構成するための要件は、測位機能が、最良の条件でもおよそ−１７９ｄＢＷまで低下することを示唆する。もう１つの欠点は、この種の再構成には時間がかかることである。必要な航法データのビットが受信されるまでにたいてい８秒かかる。Ｃ／Ａコードのシフトの測定は−１８５ｄＢＷというさらに低い信号強度で行うことが可能であるため、移動端末がＴＯＷを再構成する必要のない解決法が望ましい。
さらに、ＡＧＰＳの両タイプとも、測距信号に関して測定されたタイミングは１ミリ秒を法として切り捨てられる端数である。１ミリ秒は、３００ｋｍの距離に相当する。移動端末自身またはネットワークの位置サーバのどちらにおいて移動端末の位置を計算するときでも、正確な移動端末位置を計算するためには、移動局の位置についての事前の情報および移動局によって決定された測距信号の測定結果を用いて、完全な擬似距離を再構成する必要がある。擬似距離の再構成には一定量の不確実性が伴うため、計算された移動端末の位置にエラーをもたらす可能性がある。
従って、ＡＧＰＳシステムにおいて、測位目的で衛星の測定を実行するＧＰＳを装備した移動端末は、通常、測定されたＧＰＳネットワーク時刻と、移動端末に関する十分な測位を可能にする各衛星への擬似距離とに関する十分な知識を持っていない。このため、時刻と擬似距離のパラメータを再構成することに解決策（ソリューション）が依存する結果となり、再構成手順が、決定された端末位置にエラーをもたらす。これまでは十分なエラー訂正メカニズムがなく、背景技術のセクションで説明したとおり、衛星を利用した移動端末の測位についての先行技術による解決策は、一般に、性能や精度の限界を伴い、あるいは（又は加えて）複雑な計算を必要とする。
本発明は、移動端末で測定された時間及び擬似距離のパラメータのうちの少なくとも一方に関してエラーを訂正する機能の付いた、進化した測位メカニズムを提供する。主な思想は、基本的に、必要なときだけエラー訂正を行い、それによって移動端末の位置を計算する処理の複雑性を最小化するが、一方では決定されたＵＥの位置パラメータについて十分な正確さを実現することである。このため、本発明は、移動端末の位置を、まずエラー訂正無し、または部分的なエラー訂正により、すなわち完全なエラー訂正手順のすべての部分を使うことなく、計算することを提案する。その後、結果の品質（すなわち、使用されたパラメータの推定値の品質）が決定され、所定の品質基準、例えば適切な品質パラメータに関する所定の閾値、と比較される。所定の品質基準が満たされた場合には、それ以上の計算は必要ない。他方、最初の計算の品質が不十分な場合には、今度は移動端末で測定されその位置の計算に使われた少なくとも１つのパラメータ（時間／疑似距離）についてのエラー訂正を追加して、移動端末の位置が再計算される。
エラー訂正を「追加」することによって、再計算のステップが、以前の位置計算にはなかった少なくともひとつのエラー訂正のための測定を導入することが理解できよう。言い換えると、再計算のステップは、最初のエラー訂正の測定として、または追加的なエラー訂正の測定として導入された、異なる種類のエラー訂正を含む。追加的なエラー訂正の測定は、計算のステップにおけるこれまでの部分的なエラー訂正の代わりとなるもの、或いは、それと組み合わされたものである。
本発明に従って移動端末の位置を計算（または再計算）するとき、ＵＥの位置およびＵＥクロックのバイアスについての未知のパラメータを含む非線形連立方程式が定義される。時刻パラメータに関するエラー訂正が含まれている場合は、送信時刻のオフセットについての未知のパラメータが追加される。後者は、決定されたＵＥ位置の次元が１だけ減らされるのでない限り、１つの衛星測定の追加、すなわち１つの非線形方程式の追加が必要だということを示唆する。非線形連立方程式は線形化され、次いで、好ましくは加重最小二乗法という意味で従来の方法に従って解かれる。
本発明に従って行われるエラー訂正は、（所定の品質基準がすでに満足されているのでない限り、）信号受信時刻のパラメータにおけるエラーの訂正、または疑似距離（端数部分であるか否かを問わない）のパラメータにおけるエラーの訂正、あるいはその両方の組み合わせを含むことが可能である。図３Ａおよび３Ｂを参照しながら以下に詳述する一部の好適な実施形態は、エラー訂正に段階的な解決策を用いる。この場合、事前に定義された２以上の品質基準との比較がプロセスの異なる段階で行われることがある。まず、移動端末での信号受信時刻が利用可能な場合には、そのパラメータに基づいて測位が実行される。結果の品質が第１の基準に合致しない場合、信号受信時刻の測定は不要であるような、例えば送信時刻のオフセットの推定を含み得る、新規の計算が行われる。結果の品質が第２の基準を用いて再度チェックされ、品質が不十分な場合には、送信時刻のオフセットの推定に加えて（またはその代わりに）、今度は疑似距離のエラー訂正を伴って位置が再度再計算される。
本発明は、エラー訂正の強化を通じて、移動端末の効率的な測位をもたらし、それにより、一層正確な位置計算をもたらす。エラー訂正方法が導入されると、一般的に、より長い時間がかかる複雑な計算が必要になり、一層の処理能力が必要になる。これはＡＧＰＳに関しては、特に問題になる。本発明はこの問題を認識しており、品質の検査を伴う高度な解決策を提供する。この解決策を使用すれば、品質が非常に悪い時だけエラー訂正を実行することが可能である。このようにして、複雑で時間のかかる計算が回避できることが多くなり、全体のシステム性能も改善される。さらに、以下で明らかになるように、本発明で提案する品質の検査により、エラー訂正の際に追加して衛星を測定することが最小限で済む。本発明で提案する解決策は、携帯電話やそれに類似する装置の測位を単純化する方法を常に模索しているＡＧＰＳにとって非常に好ましいものである。
以上のことに続いて、本発明は基本的に、エラー訂正の程度と種類に関して、３種類の異なる位置計算を伴う。すなわち、
Ｉ）測定された／報告された受信時刻を用いた位置計算
ＩＩ）推定された送信時刻のオフセットを用いた位置計算
ＩＩＩ）疑似距離の外れ値(outlier)の訂正を伴う位置計算
である。
本発明で提案する手順は、個別に行われるまたは以下にその例をあげる方法で組み合わせて行われる、これらの計算メカニズムの一部または全部を使うことが可能である。しかし、発明を完全に理解するため、各種類の計算の原理についてまず例をあげて説明する。
＜測定された／報告されたｔｕを用いた位置計算の例（Ｉ）＞
この計算メカニズムは、移動端末での衛星信号の受信時刻、すなわちｔｕが利用可能であるときに使用可能であり、移動端末またはネットワーク内のどちらかで実行される。後者の場合、衛星の測定を実行する移動端末によって受信時刻がネットワークに報告されることが必要である。
完全な疑似距離が未知の場合、測定された端数の疑似距離から再構成される必要がある。例えば、（例えば、ＵＥアシストＡＧＰＳにおいては）移動端末はしばしば、Ｃ／Ａコードの１周期を法とした端数の擬似距離のみを報告する。すなわち、疑似距離（式（１））におけるコード周期の整数は未知であり、再構成される必要がある。Ｒが一般に測距信号のタイミングを測定するときに使用される切り捨ての間隔に相当する距離を表すとする。一例として、Ｃ／Ａコード１周期（１ミリ秒）の切り捨て間隔に相当する距離Ｒは、

Ｒ＝ｃ・１０−３ （４）

である。
従って、

ρｉ＝ｋｉＲ＋ｖｉ （５）

となる。ここで、ｖｉは報告された測定であり、０≦ｖｉ＜Ｒを満たす。また、整数ｋｉは再構成される必要がある。再構成は以下のように行われる。ＵＥに関する事前情報である位置ｘｕ０および概算の時刻ｔｕは既知であると仮定する。
まず、予測された疑似距離を、

ρｉ’＝｜ｘｓｉ−ｘｕ０｜ （６）

のように決める。
ＳＶの位置ｘｓｉは、例えば時刻ｔｕ−７７ミリ秒で計算されてもよい。これは通常、疑似距離の再構成という目的には十分である。式（２）の使用により、式（５）に従う「再構成された疑似距離」は、

ｋｉＲ＋ｖｉ＝｜ｘｕ−ｘｓｉ｜＋ｂ＋ｅｉ
＝ρｉ’＋ｂ＋｜ｘｕ−ｘｓｉ｜−ρｉ’＋ｅｉ＝ρｉ’＋ｂ＋Δｉ＋ｅｉ （７）

を満たすことが分かる。
ここで、ｋｉは整数であり、ｖｉは０≦ｖｉ＜Ｒを満たし、Δｉ＝｜ｘｕ−ｘｓｉ｜−ρｉ’はＳＶ＃ｉ（ｉ番目の衛星）の初期距離に関する不確実性である。もし式（７）の中のDｉ、ｂ、及び、ｅｉが小さいなら、単純にこれらの項はゼロであると仮定し、すべてのｉについて結果としての式（７）を最も良く満たす整数ｋｉを決定してもよい。整数ｋｉは最終的に四捨五入される必要がある。Dｉとｅｉは（量子化の間隔Ｒに比べると）小さいと推定するのは合理的だが、式（７）におけるクロックのバイアスｂはほぼＲ程度であるか、Ｒより大きい可能性もあり、再構成エラーにつながることもある。しかし、位置計算の本質的な課題は、相対的な疑似距離（すなわち、ρｉ−ρ１の差）が正確に再構成されることである。なぜなら、いかなるありふれた再構成エラーも、位置計算の関数におけるクロックのバイアスパラメータｂによって吸収されるからである。再構成された疑似距離間の差ρｉ−ρ１は、次式を満たす。

ｋｉＲ＋ｖｉ−ｋ１Ｒ＋ｖ１＝ρｉ’＋Δｉ’＋ｂ＋ｅｉ−ρ１’−Δ１’−ｂ−ｅ１
＝ρｉ’＋Δｉ’＋ｅｉ−ρ１’−Δ１’−ｅ１ （８）

ここで、ｉ＝２，・・・，ｎであり、ｎはＳＶの番号である。なお、クロックのバイアスｂは式（８）から消去されている。式（８）はｎ−１の等式を定義するが、未知のｋｉがｎ個ある。従って、代わりに式（７）を使って、

ｋ１＝ｒｏｕｎｄ｛（ρ１’−ｖ１）／Ｒ｝ （９）

となるようにｋ１を定義する。
そして、式（８）にｋ１を挿入し、

ｋｉ＝ｒｏｕｎｄ｛（ρｉ’−ρ１’−ｖｉ＋ｋ１Ｒ−ｖ１）／Ｒ｝ （１０）

を得る。
式（７）の再構成された疑似距離を用いて、ＳＶの信号送信時刻を推定できる。式（１）からは、

ｔｔｉ＝ｔｕ−ρｉ／ｃ （１１）

が得られる。
次に、そのときのＳＶの位置を計算する。これは、［１］に記載されている標準式を用いて実行できる。得られた座標は、転送時間の間の地球の自転について補償する必要がある。よく知られた式が［１］に掲載されている。転送時間はおよそrｉ／ｃである。得られたＳＶの座標は一般にＥＣＥＦ座標で表される。しかし、位置計算には、例えば、ＨＤＯＰのような精度測定を計算するときや、地平座標だけが必要なときには、局所的な接線座標系を使った方が便利である。従って、ＳＶの位置座標は、概算のＵＥ位置を中心とした、単位ベクトルがそれぞれ東、北、上を指している座標系に変換される。そのような変換を行う方法は当業者に良く知られている。
次に、疑似距離の測定式（２）を、

ρ＝｜１ｎ・ｘｕ−Ｘｓ｜＋ｂ１ｎ−ｅ （１２）

のようにベクトル形式で表す。
ここで、rは再構成された疑似距離を含む長さｎの列ベクトルであり、１ｎは要素として１だけを含む長さｎの列ベクトルである。Ｘｓは、ｉ番目の行が、送信時刻ｔｔｉにおけるｉ番目のＳＶの座標ｘｓｉを含む行列式である。ここでは括弧内の行列式Ｚの各行についてノルム｜Ｚ｜が計算されると想定されている。
未知のパラメータｘｕおよびｂの最初の推定値近傍のテイラー級数展開は、

ρ＝｜１ｎ・ｘｕ−Ｘｓ｜＋ｂ＋ｅ
＝｜１ｎ・ｘｕ０−Ｘｓ｜＋ｂ０１ｎ＋Ｇ（（ｘｕ−ｘｕ０）ｂ）Ｔ＋ｖ （１３）

で表現可能である。ここで、Ｇはパラメータｘｕおよびｂに関する疑似距離の導関数を含む幾何学行列であり、ｖは測定エラーの項ｅと高次テイラー級数の項との和である。ここで、

ｒｉ＝｜ｘｕ０−ｘｓｉ｜ （１４）

であるとする。
そして、Ｇはｉ行が

Gi=[(xu0-xsi)/ri (yu0-ysi)/ri (zu0-zsi)/ri 1] （１５）

に等しい行列式である。
式（１３）に対する最小二乗解は、

(xu b)T=(xu0 b0)T+(GTG)-1GT(ρ-|1n・xu0-Xs|-b01n) （１６）

に等しい。
この最小二乗解の最小損失関数値は、

V=(ρ-|1n・xu0-Xs|-b01n)T(I-G(GTG)-1GT)(ρ-|1n・xu0-Xs|-b01n) （１７）

である。
最小損失関数値は、更新されたパラメータ推定値に基づいて予測された疑似距離が、測定された疑似距離にどの程度合致するかの尺度である。従って、パラメータ推定値の質の指標／尺度として使うことができ、従って、位置計算の質の指標として使うことができる。
式（１６）の解は、重み付けされていない最小二乗解である。すなわち、最小化手順においてすべての測定が等しい重要度を付与されている。しかし、もし疑似距離測定の精度がそれぞれ異なるのなら、正確な測定に一層の重みを置くほうが一般に望ましい。加重最小二乗解は以下で与えられる。疑似距離測定の共分散行列が、

Ｐ＝Ｅ｛ｅｅＴ｝ （１８）

あるとする。ここでＥ｛｝は括弧内の量の期待値であり、ｅは測定エラーである。次に、式（１３）に対する加重最小二乗解は、

(xu b)T=(xu0 b0)T+(GTP-1G)-1GTP-1(ρ-|1n・xu0-Xs|-b01n) （１９）

である。
この場合、最小損失関数値の品質基準の例は、

V=(ρ-|1n・xu0-Xs|-b01n)T(P-1-P-1G(GTP-1G)-1GTP-1)(ρ-|1n・xu0-Xs|-b01n)
（２０）

になる。
測定エラーについての根底にある仮定が正しければ、Ｖはｎ−ｍの自由度を持つχ２分布に従って分布する。ここで、ｎは疑似距離測定の数であり、ｍは推定されたパラメータの数である。（測定され報告されたｔｕを用いた説明例の中では、ｘｕ、 ｙｕ、 ｚｕ、ｂが推定されるのだから、ｍ＝４である。）この知識は、不確実性の尺度を定義するため、および、下記により詳細に説明する本発明の特定の実施形態に従う測定データ内の外れ値の存在をテストするために、使用可能である。
一般に、式（１９）を解くことを数回繰り返すことが必要かもしれない。その場合、各ステップは、パラメータベクトルが一点に収束するまで、直近のパラメータ推定値に基づいて新規のＧ行列を使用する。
＜推定されたτを用いた位置計算の例（ＩＩ）＞
今度は、正確な測定時刻ｔｕが位置計算のために利用できない状況に目を向ける。これは例えば、ＵＥアシストＡＧＰＳにおいて、ｔｕパラメータがＵＥからネットワークに報告されていない場合に発生する。あるいは、測定が誤っている（精度が低すぎる）と推定される場合、ＵＥがｔｕを使用／測定しないことを選択したり指示されたりした場合もある。そのような場合には、本発明は、想定される信号送信時刻ｔｔｉの中で定数のオフセットτを推定するために冗長な測定を使用する計算メカニズムを提案する。
正確な測定時刻ｔｕが利用できないときには、式（１１）を用いて十分な精度を持つ送信時刻ｔｔｉを再構成することができない。その代わり、数秒間違っている可能性がある精度のより低いｔｕの値を初期値として使用する。それから、再構成の必要がある疑似距離が、報告された端数の疑似距離、初期位置ｘｕ０、および初期ｔｕを用いて、例えば式（４）乃至（７）に従って再構成される。その後、式（１１）からｔｔｉの初期値を決めることができる。
名目ＧＰＳ時刻ｔｔｉにおけるＳＶの位置は、その後、標準的なＧＰＳの公式を用いて決定できる。例えば［１］を参照されたい。ＳＶの位置座標の時間に関する微分によって３次元の速度ベクトルが得られる。微分は従来どおりの直接的な方法で実行できる。得られた座標はＥＣＥＦ座標として表される。続いて行う位置計算を単純化するため、前に例にあげた位置計算のときのように、位置ベクトルおよび速度ベクトルは、局所的な接線座標系に変換されることが望ましい。従って、局所的な接線座標系におけるＳＶの動きの線形モデルが、

ｘｓｉ（ｔｔｉ＋τ）＝ｘｓｉ（ｔｔｉ）＋ｖｘｓｉ（ｔｔｉ）・τ （２１）
ｙｓｉ（ｔｔｉ＋τ）＝ｙｓｉ（ｔｔｉ）＋ｖｘｓｉ（ｔｔｉ）・τ （２２）
ｚｓｉ（ｔｔｉ＋τ）＝ｚｓｉ（ｔｔｉ）＋ｖｘｓｉ（ｔｔｉ）・τ （２３）

のように得られる。
次に、測定式（２）が、

ρ＝｜１ｎ・ｘｕ−Ｘｓ（τ）｜＋ｂ＋ｅ （２４）

のようにベクトル形式で表される。
ここで、rは疑似距離を含む長さｎの列ベクトルであり、１ｎは要素として１だけを含む長さｎの列ベクトルである。Ｘｓ（τ）は、ｉ番目の行が、時刻ｔｔｉ＋τにおけるｉ番目のＳＶの座標ｘｓｉを含む行列式である。ここでは括弧内の行列式Ｚの各行についてノルム｜Ｚ｜が計算されると想定されている。
未知のパラメータｘｕ、τおよびｂの最初の推定値ｘｕ０、ｂ０およびτ０の近傍の式（２４）のテイラー級数展開は、

ρ＝｜１ｎ・ｘｕ−Ｘｓ（τ）｜＋ｂ＋ｅ＝｜１ｎ・ｘｕ０−Ｘｓ（τ０）｜＋
ｂ０１ｎ＋Ｈ・（（ｘｕ−ｘｕ０）τ−τ０ｂ−ｂ０）Ｔ＋ｖ （２５）

で表されてもよい。
Ｈは、パラメータｘｕ、τおよびｂに関する疑似距離の導関数を含む幾何学行列である。ここで、

ｒｉ＝｜ｘｕ０−ｘｓｉ（τ０）｜ （２６）

であるとすると、Ｈはｉ番目の行が次式に等しい行列式である。

Hi=[(xu0-xsi(τ0))/ri (yu0-ysi(τ0))/ri (zu0-zsi(τ0))/ri {vxsi(τ0)・(xu0-
xsi(τ0))+vysi(τ0)・(yu0-ysi(τ0))+vzsi(τ0)・(zu0-zsi(τ0))}/ri 1] （２７）
移動端末の位置は、式（２４）に対する加重最小二乗解を用いて与えられることが望ましい。すなわち、

(xu τb)T=(xu0 τ0b0)T+(HTP-1H)-1HTP-1(ρ-|1n・xu0-Xs(τ0)|+b01n) （２８）

であり、最小損失関数値の形式による品質基準は、

V=(ρ-|1n・xu0-Xs(τ0)|+b01n)T(P-1-P-1H(HTP-1H)-1HTP-1)(ρ-|1n・xu0-Xs(τ0)|+b01n)
（２９）
となる。
前述同様に、式（２８）を解くことは、パラメータベクトルが一点に収束するまで繰り返す必要があることが多い。
＜外れ値の訂正を伴う位置計算の例（ＩＩＩ）＞
移動端末が非常に大きなエラーを伴う疑似距離を測定したり、または報告したり、或いは、測定及び報告したりすることがときどき発生する。比較的大きなエラーに関連する測定を、外れ値と呼ぶことがある。これは、本明細書の中でそのような誤りのある値／パラメータを意味するのに使われる用語である。外れ値による疑似距離の確率は、ＡＧＰＳの操作モードでは特に高い。なぜなら、普通は、ＡＧＰＳには推定した疑似距離の妥当性チェックがないからである。さらに、ＡＧＰＳ受信機は通常、比較的低い信号レベルで測距信号を測定するため、誤った相関ピークを選択するリスクが大きい。これらの理由から、受信した疑似距離測の定において、ＡＧＰＳの位置計算機能が外れ値を検出できることが望ましい。
本発明の特に有利な実施形態によって、報告／測定された疑似距離が悪いことに起因する測位エラーの訂正メカニズムが提供される。提案する外れ値検出のためのスキームは、例えば式（２９）または式（２０）のどちらか関係がある方の損失関数値Ｖのような、パラメータ推定値の質の尺度に基づいている。パラメータ推定値の質は、外れ値の疑似距離があるかどうかを見るために、適切な品質基準（例えば、パラメータ推定の質の測定のための閾値）と対照してチェックされる。
損失関数値を品質の尺度にする例において、品質基準は、理想的なケースでは損失関数がどのように分布すべきかについての知識に基づいて有利に決定できる。上述したように、理想的には、Ｖはｎ−ｍの自由度を持つχ２分布に従って分布する。ここで、ｎは疑似距離の測定の数であり、ｍは推定パラメータの数である。（式（２０）の例では、ｍ＝４だが、式（２９）の例では追加パラメータτがあるためｍ＝５である。）失敗の警告(false alarm)に関する望ましい確率に基づいて、この分布に関係する閾値が選択される。報告されたＶの値がＶ＞Ｖｔｈｒである場合、すなわちパラメータ推定値の品質が低すぎる場合、位置計算に使われた疑似距離の測定の集合の中に少なくとも１つの外れ値が存在するという結論になる。
本発明による外れ値検出および訂正は、例えば、反復計算手順を含むことができる。その中では位置計算がｎ回繰り返される。ｎは疑似距離の測定値の数である。そしてその場合、ｋ番目の反復において、ｋ番目の疑似距離は位置計算から除外される。得られたｎ個すべての解の品質が評価され、ｋ番目の反復で最高品質が達成されたと仮定し、また、結果の品質が受け入れ可能と考えられたと仮定すると、ｋ番目の測定値が誤りであるという結論が出される。ｋ番目の測定値を除外した位置計算に基づいて、測位要求を行っているノードに戻される最終的な位置が決定される。
外れ値訂正によって、移動端末の測位の精度を大幅に高めることができる。外れ値の排除は通常、例えば式（２８）に従って、位置の追加計算を多数含むことになるため、この手順の本当の利益は、外れ値の擬似距離の訂正を上記の品質チェックと組み合わせることによって得られる。追加計算は従って、事前設定された品質基準に基づいて決定される、外れ値の訂正が必要な場合に限定される。
図３Ａは、本発明の実施形態の例に従う無線通信ネットワーク内の移動端末の位置を決定する方法を示すフローチャートである。この位置計算例は、誤った受信時刻ｔｕのリスクと外れ値の疑似距離ρのリスクの両方を、好ましい段階的な手順の中で説明している。手順は最初、移動端末で測定された各衛星信号の受信時刻ｔｕが位置測定に利用可能かどうかによって２種類の異なる経路をとる。ネットワークベースのシナリオ（移動局アシストＡＧＰＳ）では、これは移動端末によってｔｕがネットワークに報告されたかどうかをステップＳ１の手順が尋ねることを意味する。位置計算が移動端末（移動局アシストＡＧＰＳ）で実行されるシナリオにおいては、受信時刻ｔｕは常に移動端末での測定に利用可能である。そうは言っても、移動端末がｔｕを測定／使用しないことを選択する場合があるかもしれないため、これはステップＳ１で取り扱う。
従って、ステップＳ１は、ｔｕが移動端末によって報告されるかどうか（移動局アシストＡＧＰＳ）または移動端末で測定されるかどうか（移動局ベースのＡＧＰＳ）をチェックする。もしそうでない場合（ＮＯの場合）は、直ちに移動端末の位置がステップＳ４’に従う送信時刻のオフセットの推定を使って計算される。しかし、ｔｕが利用可能な場合は、例えば上記の式（１９）で概略が説明されているように、報告された／測定されたｔｕを使ってステップＳ２に従って、まず位置が計算される。ステップＳ３では、ステップＳ２の結果の品質が、ここでは損失関数の基準で例示されている所定の品質基準と対照してチェックされる。そして、ステップＳ２からの端末位置についての損失関数値Ｖが決定され、所定の閾値Ｖｔｈｒ１と比較される。所定の品質基準が満たされた場合、例えば結果としての損失関数値Ｖが閾値Ｖｔｈｒ１より小さい場合は、ステップＳ２からの移動端末位置が最終位置となり、それ以降のエラー制御は不要である。位置は通常、不確実性の範囲と共に報告される。他方、結果としての損失関数値Ｖが所定の閾値より大きい場合、これは報告されたｔｕが間違っているか、疑似距離ρの測定値の中に外れ値があるか、あるいはその両方か、ということを意味する。手順は続行され、例えば疑似距離の外れ値訂正（ステップＳ４）あるいは送信時刻のオフセットτの推定（ステップＳ４’）の形で追加されたエラー訂正を伴って移動端末の位置を再計算する。例えば、ステップＳ４’は、上記の式（２８）で概説したように実行できるし、他方、ステップＳ４は、ｔｕを用いて計算した位置を使って図４に関連して述べる疑似距離の外れ値訂正のための反復計算手順を使うことができる。
ステップＳ５において、ステップＳ４またはＳ４’で計算／再計算された端末位置の品質が決定され、ここでは損失関数の別の最大閾値Ｖｔｈｒ２として例示される所定の品質基準と対照してチェックされる。パラメータ推定値の品質がＶｔｈｒ２より小さい場合には、得られた位置が不確実性の範囲と共に報告され、手順はステップＳ５の後で終了されてもよい。そうではなく、損失関数が所定の閾値より大きい場合には、測定時刻または疑似距離のパラメータに考慮が必要なエラーがまだ存在すると推定される。その場合、移動端末の位置は、ステップＳ６で疑似距離のデータ中の外れ値の訂正と組み合わせて送信時刻のオフセットのパラメータ（ｔｕのエラーを補償するためのパラメータ）を使って再計算される。ステップＳ６の位置計算は、例えば、τを使って位置を計算する図４の例のような、疑似距離の外れ値の訂正のための反復計算手順を含んでもよい。
上述の手順に代わる選択肢として、ステップＳ３の品質基準が満たされない場合に直接ステップＳ６に行くことも可能である。これは、図３Ａにおいて点線矢印で示されている。位置計算の品質が不十分だと分かった時点ですぐに両タイプのエラー訂正を組み合わせて実行することも十分可能だが、そのような実施形態は一般に、位置計算に関する計算負荷の増加につながり、ステップＳ４，Ｓ４’における個別の計算以上の衛星の測定が必要になる。従って、上述のように、エラー訂正を段階的に導入する方が通常は好ましい。
加えて、本発明の特に有利な実施形態に従って、エラー訂正をさらなるサブステップに分割することが提案されている。これを、図３ＡのポイントＡとポイントＢの間に挿入される代替用の下位手順(sub-procedure)として図３Ｂのフローチャートに図解する。ステップＳ４−１およびＳ４−１’はそれぞれ、ステップＳ４およびＳ４’にそれぞれ直接対応している（図３Ａ）。しかし、その後の品質チェックで（ステップＳ４−２／Ｓ４−２’）品質が不十分だと分かった場合、他のエラー訂正タイプを使って位置が再計算される（ステップＳ４−３／Ｓ４−３’）。他のエラー訂正タイプとは、疑似距離の外れ値の訂正（ｔｕによって）が最初に使われた場合には送信時間のオフセットの推定を使った再計算であり、逆の場合も同様である。これらの実施形態に従うと、ステップＳ６（図３Ａ）の複合エラー訂正は、ｔｕまたは疑似距離の外れ値におけるエラー訂正を使った個別メカニズムのいずれも単独では十分な品質を実現できない場合にだけ使われる。そのような場合、第３の品質基準との比較が行われる。第３の品質基準は、例示されたＶｔｈｒ２としての第２の基準に等しいこともあれば異なることもあり、その基準が満たされないならば第３の再計算が続くこともある。
注意すべきことだが、外れ値の訂正ステップＳ４−１、Ｓ４−３’はｔｕを用いて実行されるため、図３Ｂの実施形態はｔｕが利用可能な（報告された／測定された）場合、すなわち、ステップＳ１（図３Ａ）の「ＹＥＳ経路」に相当するシナリオのみに適用される。
疑似距離の外れ値の訂正（図３ＡのステップＳ４、Ｓ６、および図４ＢのステップＳ４−１、Ｓ４−３’）を含む本発明による位置計算は、例えば図４の例の手順のような反復計算手順を含んでもよい。外れ値の訂正は、移動端末の位置ｘｋ（ｋ＝ｌ， ２，．．， ｎ）の計算を伴うｎ個の場合を含む。ここで、ｎは利用可能な疑似距離の測定値の数（ステップＳ６−１）である。ステップＳ６−１の計算は、ｋ番目の場合にｋ番目の疑似距離の測定値が位置計算から除外されるように実行される。この位置計算関数は、例えば式（１９）に従ってｔｕを使うか、または送信時刻のオフセットの推定値τも含む（すなわち、測定した／再構成したｔｕおよび例えば式（２８）を用いたρにおけるエラー訂正の組み合わせを含む）かのどちらかである。ステップＳ６−２において、パラメータ推定品質は計算された位置ｘｋについて決定される。ステップＳ６−３は、利用可能な疑似距離パラメータがもっとあるかどうか尋ねる。もしもっとあれば、ｋが１ずつ増加し、ステップＳ６−１およびＳ６−２が次の疑似距離についても反復される。もし利用可能な疑似距離がすべて使用されたのであれば（ｋ＝ｎ）、手順は、決定されたパラメータ推定値品質の評価を続ける（ステップＳ６−４、Ｓ６−５、Ｓ６−７）。これによって比較的高品質が得られれば、それは除外された疑似距離が外れ値である印であると解釈される。外れ値の疑似距離が検出されれば、この疑似距離は測定集合から除外され、最終的な端末位置の決定には含まれない。損失関数の閾値を使った例において、ｎ回の計算の成果はｎ個の損失関数値Ｖｋである。これらのｎ個の損失関数値のうち最小の損失関数値ｍｉｎＶが閾値Ｖｔｈｒより小さい場合、外れ値がうまく検出され除去されたと結論付けられる（Ｓ６−５、Ｓ６−６）。そうでなければ、外れ値は除去できなかったと結論付けられ、それによって通常は、すべての利用可能な疑似距離によって計算された位置が使用される（Ｓ６−７、Ｓ６−８）。手順は終了され、得られた位置が好ましくは不確実性の範囲と共に報告される。
図３Ａ（及び図３Ｂ）の段階的に導入されたエラー訂正は非常に有利である。なぜなら、計算および測定の負荷の増加を最小限に留めて、計算を受け入れ可能な複雑度に維持しながら、より正確な測位を提供するからである。反復的な外れ値の拒絶手順に関する計算が追加されるのは、例えば、悪い結果は誤った受信時刻に起因するという可能性が既にチェックされ除外された場合に限定することが可能である。
背景技術のセクションで述べたとおり、［３］は、真のＧＰＳネットワーク時刻と衛星の測定による概算の時刻との差を、移動端末の位置を計算する際に追加の未知数として用いることによって、移動端末がＧＰＳネットワーク時刻を知らない状況を扱う手法について説明している。これには、衛星の測定がさらにもう１度必要になるため、結果として、処理時間と電力に関して多くを要する計算が必要になる。［３］に従えば、この種の拡張処理が常用され、移動端末及びセルラネットワークのうちの少なくとも一方の測位機能に多大な要求を突きつけることになろう。本発明は、位置計算に品質チェックを導入することによって、もっと適切な測位機能が実現できることを認識している。本発明で提案する解決策は、エラー訂正に関する多様なニーズを区別し、状況に応じて異なる度合いのエラー訂正／パラメータ推定を適用する。従って、計算に追加する複雑性を最小限度にとどめつつ、十分な測位精度を実現する。本発明の品質チェックは、疑似距離の測定値の外れ値を訂正する機能を含むことによって、さらに測位を向上させることもできる。
図３Ａ（および図３Ｂ）に従って本発明を実施するのは非常に好ましいが、当業者であれば他のフローシーケンスも可能であることはすぐに分かる。外れ値の訂正および送信時刻のオフセットの推定に関する各ステップは、例えば並行して実行できる。さらに、ステップＳ３やＳ５（およびＳ４―２、Ｓ４―２’）そしてＳ６における品質基準は、等しくても（Ｖｔｈｒ１＝Ｖｔｈｒ２）異なっていても（Ｖｔｈｒ１≠Ｖｔｈｒ２）良いし、固定であっても可変であっても良い。損失関数を持つ例において、少なくとも自由度が異なる計算については異なる品質基準を用いるほうが一般に適切であろう。
外れ値の訂正にも多様な修正が可能である。１つの選択肢は、１つの測定値を除外して位置計算の結果を連続的に評価して行き、所定の品質基準に合致するまで反復手順を続けることであろう。そうすれば図４の反復ループは評価ステップＳ６−４、Ｓ６−５、Ｓ６−７も含むことになる。他の実施形態においては、基準が満たされた後であっても手順を継続して、可能なら品質をさらに高めても良い。外れ値の訂正は１つの測定値を除去することに限定しても良いし、十分な数の衛星の測定値がある場合には複数の外れ値の疑似距離を除去できるようにしても良い。
場合によっては、本発明の好適な実施形態、例えば図３ＡのステップＳ６に従って必要とされるすべての試験を行うのに十分な疑似距離の測定値がないこともありうる。その場合には、位置パラメータの空間を２次元位置（地球接線座標系のｘ座標とｙ座標）へと減らすことで、統合性を維持し、ｔｕの欠落した測定値に対処し、不確実性の範囲に関して信頼できる推定値を提供することが適切であろう。移動端末の位置の計算に送信時刻のオフセットの推定が含まれるとき、端末が３次元で決定されるには通常は５つの衛星測定値が必要である（そして、追加して外れ値のエラー訂正を行うにはさらに測定値が必要である）。２次元の位置決定の場合には、すなわち、通常はわずかな差異しかない垂直のｚ座標を省略することによって、測定されたｔｕのエラー訂正が必要なときでさえ、４つの衛星の測定値で十分であろう。
図５は、本発明による移動端末の典型的なハードウェア実装を図解する概略ブロック図である。移動端末５００は、ＧＰＳ ＲＦフロントエンド５１０、測位モジュール／プロセッサ５２０、セルラ通信モジュール５３０（セルラＲＦモジュール５３２およびセルラベースバンドプロセッサ５３４を備えている）、および、それぞれセルラネットワークとＧＰＳシステムとの通信のためのアンテナ５１２、５３６から成る。セルラ通信モジュール５３０は、セルラネットワークから支援データを受信する。支援データは可視衛星用のエフェメリスとクロック訂正値、概算の（おおよその）ＵＥ位置、および概算のＧＰＳシステム時刻から成っても良い。あるいは、支援データは、相関処理を支援するためだけの明確な支援情報を含んでいても良い。どちらの場合にも、支援データは測位モジュール５２０に送信される。通信モジュール５３０は、ＧＰＳ ＲＦフロントエンド５１０と測位モジュール５２０に参照クロックを提供する。ＲＦフロントエンドモジュール５１０は、測位モジュール５２０によって制御される。
測位モジュール５２０は、通信モジュール５３０から測位要求を受信すると、ＧＰＳ信号のサンプルを収集するようＲＦフロントエンドモジュール５１０に指示する。ＧＰＳ ＲＦフロントエンド５１０は、アンテナ５１２を通じてＧＰＳ周波数帯を受信し、信号をベースバンドにダウンコンバートし、信号を同相（Ｉ）成分と直交（Ｑ）成分に分離し、信号をサンプリングしてデジタル形式に変換し、そしてＩ／Ｑインタフェースを通じてこれらを測位モジュール５２０に出力する。測位モジュールは受信したＩ、Ｑデータをメモリ５２４に保存し、データを十分収集した後、ＤＳＰ（デジタル信号プロセッサ）５２６内で相関動作が生じる。相関器はＣＰＵ５２２によって制御され、ＣＰＵ５２２は今度は支援データを使って必要な相関検出動作の量を減らす。十分なＣ／Ａコード境界が検出されたとき、またはタイマーがタイムアップしたとき、測位モジュール５２０は測定された端数の疑似距離を通信モジュール５３０に出力する。通信モジュール５３０は、位置計算のため測定値をセルラネットワークに送信する（移動局アシストＡＧＰＳ）。あるいは、ＣＰＵ５２２において、利用可能なエフェメリスとクロック訂正パラメータ、概算のＧＰＳシステム時刻およびＵＥの概算の位置を用いて位置が計算されても良い（移動局ベースのＡＧＰＳ）。
セルラネットワーク内の測位機能については、ＲＮＣに関連する（または非ＷＣＤＭＡシステムではＲＮＣに相当するネットワーク機能に関連する）測位ノード（図１の１１７）があることが好ましい。測位ノードはＲＮＣに接続した外部ノードであっても良いし、ＲＮＣノード内に組み込まれていても良い。測位ノードは、測位要求を受信するように適合されていることが望ましく、移動端末からの衛星測定の指示に関連しても良い。そのような測定指示は一般に測位ノード内で処理される支援データを伴う。
移動局アシストＡＧＰＳについては、測位ノードが、さらに、移動端末（またはその代替物）から発せられる測距信号を受信する手段、および移動端末から報告された測定結果および移動端末位置についての事前情報に基づいて移動端末の位置を計算する処理手段とを含むことが望ましい。測位のためのネットワーク機能の少なくとも一部が別の物理的ネットワークノード上に分散されている実施形態も、本発明の範囲に入る。例えば、支援データの処理、および移動端末の位置の計算をそれぞれ行う別個の処理手段があるかもしれない。
さらなるエラー訂正が適用されるべきかどうかを決定するため本発明に従って用いられる所定の品質基準は、主に最小損失関数の閾値として例示された。他の実施形態では、最小損失関数以外の品質基準を使ってもよい。品質基準は、例えば残差分散または残差共分散に基づく、移動端末の位置を提供する最小二乗解の残差に関連することが望ましいが、他の基準を使っても良い。基本的には、位置計算に影響を及ぼす予測されたパラメータと測定されたパラメータ（例えば、時刻または擬似距離のパラメータ）との間の適合に関する適切な尺度であれば、何を使っても良い。
さらに、品質基準は必ずしも不変の閾値として表される必要はなく、例えば、位置計算手順における過去の結果に依存して、あるいは決定された場所／位置に関して求められる精度に対して直接的または間接的に影響する要因に依存して変化しても良い。適切な品質基準を決めるときには、測位精度と計算上の複雑度の間に二律背反があり、往々にして、単純な閾値こそが最も適切な選択である。
本発明について特定の実施形態を例として参照しながら説明してきたが、開示された特徴と同等の物や当業者には明白な修正や変化形も本発明が対象とすることは強調すべきである。従って、本発明の請求範囲は、添付の請求項によってのみ限定される。
＜略語＞
AGPS Assisted Global Positioning System
C/A Coarse/Acquisition
CPU Central Processing Unit
DSP Digital Signal Processor
ECEF Earth-Centered Earth-Fixed
GPRS General Packet Radio Service
GPS Global Positioning System
GSM Global System for Mobile communication
HDOP Horizontal Dilution Of Precision
HOW Hand-Over Word
ISDN Integrated Services Digital Network
MS Mobile Station
PLMN Public Land Mobile Network
PSTN Public Switched Telephone Network
RAN Radio Access Network
RF Radio Frequency
SMLC Serving Mobile Location Center
SV Space Vehicle
TOW Time Of Week
UE User Equipment
UTRAN Universal Terrestrial Radio Access Network
＜参考文献＞
[1] Navstar GPS Space Segment/Navigation user Interfaces, ICD-GPS-200, Revision IRN-200C-003, １９９９年１０月１１日
[2] 3GPP TS 44. 031, v 5.6.0
[3] 米国特許6,430,415 B1, Qualcomm Inc.
本発明を利用可能な通信システムの一例の概略図である。 本発明を利用可能な通信システムの一例におけるクロックの関係を図解するタイミング図である。 本発明の実施形態の一例による無線通信ネットワークにおいて移動端末の位置を決定する方法を示すフローチャートである。 本発明の実施形態の一例による無線通信ネットワークにおいて移動端末の位置を決定する方法を示すフローチャートである。 本発明の実施形態の一例において使用される擬似距離の外れ値の訂正のための、反復（繰り返し）計算手順を示すフローチャートである。 本発明の実施形態の一例による移動端末の概略ブロック図である。

Claims (31)

  1. セルラ通信ネットワーク（１１０）において、衛星を利用した測位システム（１３０）の衛星（１３２）から信号を受信する手段を備える移動端末（１１２；５００）の位置を決定する方法であって、
    前記移動端末において、複数の衛星から信号を受信するステップと、
    前記移動端末において、信号受信時刻における前記複数の衛星への擬似距離を測定するステップと、
    前記測定された擬似距離を示すパラメータに基づいて、前記移動端末の位置を計算するステップ（Ｓ２；Ｓ４’）と、
    前記位置の計算に関する品質の尺度を判定するステップと、
    前記判定された品質の尺度を第１の所定の品質基準と比較するステップ（Ｓ３；Ｓ５）と、
    前記品質の尺度が前記第１の所定の品質基準を満たさない場合に、前記移動端末において測定されたパラメータを示す少なくとも１つのパラメータにおける誤差に対する追加的な訂正を伴って、前記移動端末の前記位置を再計算するステップ（Ｓ４；Ｓ４’；Ｓ４−１；Ｓ４−１’；Ｓ６）と、
    を備え、
    前記計算するステップ（Ｓ２）は、前記移動端末（１１２；５００）において測定された信号受信時刻に関するパラメータを使用し、
    前記再計算するステップ（Ｓ４’；Ｓ４−１’；Ｓ６）における前記追加的な訂正は、追加的な未知のパラメータとして送信時刻のオフセットのパラメータを用いて前記信号受信時刻に関するパラメータにおける誤差を訂正することを伴う
    ことを特徴とする方法。
  2. 前記再計算するステップ（Ｓ４’；Ｓ４−１’）における前記位置の再計算に関する品質の尺度を判定するステップと、
    前記再計算に関する前記判定された品質の尺度を第２の所定の品質基準と比較するステップ（Ｓ５；Ｓ４−２’）と、
    前記品質の尺度が前記第２の所定の品質基準を満たさない場合に、擬似距離の外れ値の訂正に関する反復計算手順において前記移動端末の前記位置の第２の再計算（Ｓ６；Ｓ４−３’）を実行するステップと、
    をさらに備えることを特徴とする請求項１に記載の方法。
  3. セルラ通信ネットワーク（１１０）において、衛星を利用した測位システム（１３０）の衛星（１３２）から信号を受信する手段を備える移動端末（１１２；５００）の位置を決定する方法であって、
    前記移動端末において、複数の衛星から信号を受信するステップと、
    前記移動端末において、信号受信時刻における前記複数の衛星への擬似距離を測定するステップと、
    前記測定された擬似距離を示すパラメータに基づいて、前記移動端末の位置を計算するステップ（Ｓ２；Ｓ４’）と、
    前記位置の計算に関する品質の尺度を判定するステップと、
    前記判定された品質の尺度を第１の所定の品質基準と比較するステップ（Ｓ３；Ｓ５）と、
    前記品質の尺度が前記第１の所定の品質基準を満たさない場合に、前記移動端末において測定されたパラメータを示す少なくとも１つのパラメータにおける誤差に対する追加的な訂正を伴って、前記移動端末の前記位置を再計算するステップ（Ｓ４；Ｓ４’；Ｓ４−１；Ｓ４−１’；Ｓ６）と、
    を備え、
    前記計算するステップ（Ｓ２）は、前記移動端末（１１２；５００）において測定された信号受信時刻に関するパラメータを使用し、
    前記再計算するステップ（Ｓ４；Ｓ４−１；Ｓ６）における前記追加的な訂正は、擬似距離の外れ値の訂正に関する反復計算手順を通して、前記移動端末（１１２；５００）において測定された前記擬似距離における誤差を訂正することを伴う
    ことを特徴とする方法。
  4. 前記再計算するステップ（Ｓ４；Ｓ４−１）における前記位置の再計算に関する品質の尺度を判定するステップと、
    前記再計算に関する前記判定された品質の尺度を第２の所定の品質基準と比較するステップ（Ｓ５；Ｓ４−２）と、
    前記品質の尺度が前記第２の所定の品質基準を満たさない場合に、追加的な未知のパラメータとして送信時刻のオフセットのパラメータを用いた前記信号受信時刻における誤差の訂正を伴って、前記移動端末の前記位置の第２の再計算（Ｓ６；Ｓ４−３）を実行するステップと、
    をさらに備えることを特徴とする請求項３に記載の方法。
  5. 前記第２の再計算（Ｓ４−３；Ｓ４−３’）に関する品質の尺度を判定するステップと、
    前記第２の再計算に関する前記判定された品質の尺度を第３の所定の品質基準と比較するステップ（Ｓ５）と、
    前記品質の尺度が前記第３の所定の品質基準を満たさない場合に、追加的な未知のパラメータとして送信時刻のオフセットのパラメータを用いた訂正を、擬似距離の外れ値の訂正に関する反復計算手順を通して誤差の訂正と組み合わせて、前記移動端末の前記位置の第３の再計算（Ｓ６）を実行するステップと、
    をさらに備えることを特徴とする請求項２又は４に記載の方法。
  6. 前記擬似距離の外れ値の訂正に関する反復計算手順は、前記信号受信時刻を示すパラメータにおける誤差を訂正するための、追加的な未知のパラメータとして送信時刻のオフセットのパラメータを用いた位置計算を伴うことを特徴とする請求項２又は３に記載の方法。
  7. 前記擬似距離の外れ値の訂正に関する反復計算手順（Ｓ６；Ｓ４−１；Ｓ４−３’）においては、
    １度に１つの擬似距離の測定値を除外して、前記移動端末の前記位置を繰り返し計算するステップ（Ｓ６−１）と、
    各々の計算された位置に関する品質の尺度を判定するステップ（Ｓ６−２）と、
    前記判定された品質の尺度を評価するステップ（Ｓ６−４，Ｓ６−５，Ｓ６−７）であって、比較的品質が高いということは、前記除外された擬似距離が外れ値であるということを示すステップと、
    前記評価するステップにおいて外れ値の擬似距離が検出された場合に、前記外れ値の擬似距離を取り除き、前記外れ値の擬似距離を除外して計算された位置を移動端末の位置として使用するステップ（Ｓ６−６）と、
    が順次実行されることを特徴とする請求項２乃至６のいずれか１項に記載の方法。
  8. 前記送信時刻のオフセットのパラメータは、前記複数の衛星（１３２）それぞれからの信号送信時刻と、前記衛星を利用した測位システム（１３０）のネットワーク時刻との間の時刻の差を示すことを特徴とする請求項１、２、及び４乃至６のいずれか１項に記載の方法。
  9. 前記品質の尺度は、前記移動端末の位置に関する最小二乗解の残差を用いて決定された前記パラメータの推定値の品質に関する尺度であることを特徴とする請求項１乃至８のいずれか１項に記載の方法。
  10. 前記測定された擬似距離のうちの少なくとも１つが切り捨てられた端数の擬似距離である場合に、前記移動端末（１１２；５００）における信号受信時刻に関するパラメータ、及び概算の移動端末位置に基づいて、前記切り捨てられた端数の擬似距離から完全な擬似距離を再構成することを特徴とする請求項１乃至９のいずれか１項に記載の方法。
  11. 前記計算するステップ及び前記再計算するステップそれぞれにおいては、
    少なくとも前記移動端末の位置及び前記移動端末のクロックのバイアスに関する未知のパラメータを持つ非線形の連立方程式を定義するステップと、
    前記衛星の動きが名目の衛星位置及び速度ベクトルを持つ線形なものとしてモデル化するように、前記連立方程式を線形化するステップと、
    前記移動端末の位置に関する前記連立方程式を解くステップと、
    が順次実行されることを特徴とする請求項１乃至１０のいずれか１項に記載の方法。
  12. 測定された擬似距離、及び前記移動端末（１１２；５００）における信号受信時刻に関するパラメータに基づいて決定された、衛星の送信時刻を用いて、前記複数の衛星（１３２）それぞれに関する名目の衛星位置を決定するステップと、
    ２つの水平座標及び１つの垂直座標を持つ、局所的な接線座標系における前記衛星位置を表現するステップと、
    を備えることを特徴とする請求項１１に記載の方法。
  13. セルラ通信ネットワーク（１１０）における、自分の位置を決定する手段を備える移動端末（１１２；５００）であって、
    衛星を利用した測位システム（１３０）における複数の衛星（１３２）から信号を受信する手段（５１０；５１２）と、
    信号受信時刻における前記複数の衛星への擬似距離を測定する手段（５２０）と、
    前記測定された擬似距離を示すパラメータに基づいて、前記移動端末の位置を計算する手段（５２２）と、
    前記位置の計算に関する品質の尺度を判定する手段と、
    前記判定された品質の尺度を第１の所定の品質基準と比較する手段と、
    前記品質の尺度が前記第１の所定の品質基準を満たさない場合に、前記移動端末において測定されたパラメータを示す少なくとも１つのパラメータにおける誤差に対する追加的な訂正を伴って、前記移動端末の前記位置を再計算する手段と、
    を備え、
    前記計算する手段は、前記移動端末（１１２；５００）において測定された信号受信時刻に関するパラメータを使用し、
    前記再計算する手段は、追加的な未知のパラメータとして送信時刻のオフセットのパラメータを用いて前記信号受信時刻に関するパラメータにおける誤差を訂正する手段を伴う
    ことを特徴とする移動端末。
  14. 前記位置の再計算に関する品質の尺度を判定する手段と、
    前記再計算に関する前記判定された品質の尺度を第２の所定の品質基準と比較する手段と、
    前記品質の尺度が前記第２の所定の品質基準を満たさない場合に、反復的に擬似距離の外れ値を訂正する手段を伴う、前記移動端末（１１２；５００）の前記位置について第２の再計算をする手段と、
    をさらに備えることを特徴とする請求項１３に記載の移動端末。
  15. セルラ通信ネットワーク（１１０）における、自分の位置を決定する手段を備える移動端末（１１２；５００）であって、
    衛星を利用した測位システム（１３０）における複数の衛星（１３２）から信号を受信する手段（５１０；５１２）と、
    信号受信時刻における前記複数の衛星への擬似距離を測定する手段（５２０）と、
    前記測定された擬似距離を示すパラメータに基づいて、前記移動端末の位置を計算する手段（５２２）と、
    前記位置の計算に関する品質の尺度を判定する手段と、
    前記判定された品質の尺度を第１の所定の品質基準と比較する手段と、
    前記品質の尺度が前記第１の所定の品質基準を満たさない場合に、前記移動端末において測定されたパラメータを示す少なくとも１つのパラメータにおける誤差に対する追加的な訂正を伴って、前記移動端末の前記位置を再計算する手段と、
    を備え、
    前記計算する手段は、前記移動端末（１１２；５００）において測定された信号受信時刻に関するパラメータを使用し、
    前記再計算する手段は、反復的に擬似距離の外れ値を訂正することを通して、前記移動端末（１１２；５００）において測定された前記擬似距離における誤差を訂正する手段を伴う
    ことを特徴とする移動端末。
  16. 前記位置の再計算に関する品質の尺度を判定する手段と、
    前記再計算に関する前記判定された品質の尺度を第２の所定の品質基準と比較する手段と、
    前記品質の尺度が前記第２の所定の品質基準を満たさない場合に、追加的な未知のパラメータとして送信時刻のオフセットのパラメータを用いた前記信号受信時刻における誤差の訂正を伴って、前記移動端末の前記位置について第２の再計算をする手段と、
    をさらに備えることを特徴とする請求項１５に記載の移動端末。
  17. 前記第２の再計算に関する品質の尺度を判定する手段と、
    前記第２の再計算に関する前記判定された品質の尺度を第３の所定の品質基準と比較する手段と、
    前記品質の尺度が前記第３の所定の品質基準を満たさない場合に、送信時刻のオフセットのパラメータを用いた誤差の訂正を、反復的に擬似距離の外れ値を訂正することと組み合わせたものを含む、前記移動端末の前記位置について第３の再計算をする手段と、
    をさらに備えることを特徴とする請求項１４又は１６に記載の移動端末。
  18. 前記反復的に擬似距離の外れ値を訂正することは、追加的な未知のパラメータとして送信時刻のオフセットのパラメータを用いた位置の計算をする手段を伴うことを特徴とする請求項１４又は１５に記載の移動端末。
  19. 前記擬似距離の外れ値の訂正に関する反復計算をする手段は、順次実行される、
    １度に１つの擬似距離の測定値を除外して、前記移動端末の前記位置を繰り返し計算する手段と、
    各々の計算された位置に関する品質の尺度を判定する手段と、
    前記判定された品質の尺度を評価する手段であって、比較的品質が高いということは、前記除外された擬似距離が外れ値であるということを示す手段と、
    前記評価する手段で外れ値の擬似距離が検出された場合に、前記外れ値の擬似距離を取り除き、前記外れ値の擬似距離を除外して計算された位置を移動端末の位置として使用する手段と、
    を含むことを特徴とする請求項１４乃至１８のいずれか１項に記載の移動端末。
  20. 前記品質の尺度は、前記移動端末の位置に関する最小二乗解の残差に基づくことを特徴とする請求項１３乃至１９のいずれか１項に記載の移動端末。
  21. 前記移動端末（１１２；５００）における信号受信時刻に関するパラメータ、及び概算の移動端末位置に基づいて、測定されて切り捨てられた端数の擬似距離から完全な擬似距離を再構成する手段をさらに備えることを特徴とする請求項１３乃至２０のいずれか１項に記載の移動端末。
  22. セルラ通信ネットワークにおける、衛星を利用した測位システム（１３０）の衛星（１３２）から信号を受信する手段を備える移動端末（１１２；５００）の位置を決定する手段を備える、前記セルラ通信ネットワーク（１１０）における測位ノード（１１７）であって、
    前記移動端末から、複数の衛星に対して測定された擬似距離であって測定時刻に関連する擬似距離に関するパラメータを受信する手段と、
    前記測定された擬似距離を示すパラメータに基づいて、前記移動端末の位置を計算する手段と、
    前記位置の計算に関する品質の尺度を判定する手段と、
    前記判定された品質の尺度を第１の所定の品質基準と比較する手段と、
    前記品質の尺度が前記第１の所定の品質基準を満たさない場合に、前記移動端末において測定されたパラメータを示す少なくとも１つのパラメータにおける誤差に対する追加的な訂正を伴って、前記移動端末の前記位置を再計算する手段と、
    を備え、
    前記計算する手段は、前記移動端末（１１２；５００）から報告された信号受信時刻に関するパラメータを使用し、
    前記再計算する手段は、追加的な未知のパラメータとして送信時刻のオフセットのパラメータを用いて前記信号受信時刻に関するパラメータにおける誤差を訂正する手段を伴う
    ことを特徴とする測位ノード。
  23. 前記位置の再計算に関する品質の尺度を判定する手段と、
    前記再計算に関する前記判定された品質の尺度を第２の所定の品質基準と比較する手段と、
    前記品質の尺度が前記第２の所定の品質基準を満たさない場合に、反復的に擬似距離の外れ値を訂正する手段を伴う、前記移動端末（１１２；５００）の前記位置について第２の再計算をする手段と、
    をさらに備えることを特徴とする請求項２２に記載の測位ノード。
  24. セルラ通信ネットワークにおける、衛星を利用した測位システム（１３０）の衛星（１３２）から信号を受信する手段を備える移動端末（１１２；５００）の位置を決定する手段を備える、前記セルラ通信ネットワーク（１１０）における測位ノード（１１７）であって、
    前記移動端末から、複数の衛星に対して測定された擬似距離であって測定時刻に関連する擬似距離に関するパラメータを受信する手段と、
    前記測定された擬似距離を示すパラメータに基づいて、前記移動端末の位置を計算する手段と、
    前記位置の計算に関する品質の尺度を判定する手段と、
    前記判定された品質の尺度を第１の所定の品質基準と比較する手段と、
    前記品質の尺度が前記第１の所定の品質基準を満たさない場合に、前記移動端末において測定されたパラメータを示す少なくとも１つのパラメータにおける誤差に対する追加的な訂正を伴って、前記移動端末の前記位置を再計算する手段と、
    を備え、
    前記計算する手段は、前記移動端末（１１２；５００）において測定された信号受信時刻に関するパラメータを使用し、
    前記再計算する手段は、反復的に擬似距離の外れ値を訂正することを通して、前記移動端末（１１２；５００）において測定された前記擬似距離における誤差を訂正する手段を伴う
    ことを特徴とする測位ノード。
  25. 前記位置の再計算に関する品質の尺度を判定する手段と、
    前記再計算に関する前記判定された品質の尺度を第２の所定の品質基準と比較する手段と、
    前記品質の尺度が前記第２の所定の品質基準を満たさない場合に、追加的な未知のパラメータとして送信時刻のオフセットのパラメータを用いた前記信号受信時刻における誤差の訂正を伴って、前記移動端末の前記位置について第２の再計算をする手段と、
    をさらに備えることを特徴とする請求項２４に記載の測位ノード。
  26. 前記第２の再計算に関する品質の尺度を判定する手段と、
    前記第２の再計算に関する前記判定された品質の尺度を第３の所定の品質基準と比較する手段と、
    前記品質の尺度が前記第３の所定の品質基準を満たさない場合に、送信時刻のオフセットのパラメータを用いた誤差の訂正を、反復的に擬似距離の外れ値を訂正することと組み合わせたものを含む、前記移動端末の前記位置について第３の再計算をする手段と、
    をさらに備えることを特徴とする請求項２３又は２５に記載の測位ノード。
  27. 前記反復的に擬似距離の外れ値を訂正することは、追加的な未知のパラメータとして送信時刻のオフセットのパラメータを用いた位置の計算をする手段を伴うことを特徴とする請求項２３又は２４に記載の測位ノード。
  28. 前記擬似距離の外れ値の訂正に関する反復計算をする手段は、順次実行される、
    １度に１つの擬似距離の測定値を除外して、前記移動端末の前記位置を繰り返し計算する手段と、
    各々の計算された位置に関する品質の尺度を判定する手段と、
    前記判定された品質の尺度を評価する手段であって、比較的品質が高いということは、前記除外された擬似距離が外れ値であるということを示す手段と、
    前記評価する手段で外れ値の擬似距離が検出された場合に、前記外れ値の擬似距離を取り除き、前記外れ値の擬似距離を除外して計算された位置を移動端末の位置として使用する手段と、
    を含むことを特徴とする請求項２３乃至２７のいずれか１項に記載の測位ノード。
  29. 前記品質の尺度は、前記移動端末の位置に関する最小二乗解の残差に基づくことを特徴とする請求項２２乃至２８のいずれか１項に記載の測位ノード。
  30. 前記移動端末（１１２；５００）における信号受信時刻に関するパラメータ、及び概算の移動端末位置に基づいて、測定されて切り捨てられた端数の擬似距離から完全な擬似距離を再構成する手段をさらに備えることを特徴とする請求項２２乃至２９のいずれか１項に記載の測位ノード。
  31. セルラ通信ネットワーク（１１０）において、衛星を利用した測位システム（１３０）の衛星（１３２）から信号を受信する手段を備える移動端末（１１２；５００）の位置を決定する手段を備える通信システムであって、
    請求項２２乃至３０のいずれか１項に記載の測位ノードを備える
    ことを特徴とする通信システム。
JP2006554055A 2004-02-18 2005-02-16 衛星を利用した移動端末の測位 Expired - Fee Related JP4896746B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US54519404P 2004-02-18 2004-02-18
US60/545,194 2004-02-18
PCT/SE2005/000213 WO2005078474A1 (en) 2004-02-18 2005-02-16 Satellite-based positioning of mobile terminals

Publications (2)

Publication Number Publication Date
JP2007523345A JP2007523345A (ja) 2007-08-16
JP4896746B2 true JP4896746B2 (ja) 2012-03-14

Family

ID=34860517

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2006554055A Expired - Fee Related JP4896746B2 (ja) 2004-02-18 2005-02-16 衛星を利用した移動端末の測位

Country Status (7)

Country Link
US (1) US7633436B2 (ja)
EP (1) EP1718987B1 (ja)
JP (1) JP4896746B2 (ja)
CN (1) CN1922505B (ja)
AT (1) ATE426178T1 (ja)
DE (1) DE602005013366D1 (ja)
WO (1) WO2005078474A1 (ja)

Families Citing this family (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100385996C (zh) * 2006-06-27 2008-04-30 重庆邮电大学 基于td-scdma系统的rtk定位方法
US8307680B2 (en) 2006-10-30 2012-11-13 Arcelormittal France Coated steel strips, methods of making the same, methods of using the same, stamping blanks prepared from the same, stamped products prepared from the same, and articles of manufacture which contain such a stamped product
EP1975641B1 (en) * 2007-03-29 2011-05-25 Cambridge Positioning Systems Limited Determining the change in time at a mobile terminal
JP4906618B2 (ja) * 2007-07-13 2012-03-28 株式会社エヌ・ティ・ティ・ドコモ 測位システム、測位用ｉｃチップ、測位方法及び測位プログラム
US8044853B2 (en) * 2007-12-20 2011-10-25 Qualcomm Incorporated Navigation receiver
US8085191B2 (en) * 2008-08-18 2011-12-27 Trimble Navigation Limited Position determination based on corroborated signal processing of PRN codes
JP2010060489A (ja) 2008-09-05 2010-03-18 Seiko Epson Corp 衛星軌道モデル化適否判定方法、長期予測軌道データ提供方法及び衛星軌道モデル化適否判定装置
US8330652B2 (en) * 2009-11-04 2012-12-11 Qualcomm Incorporated Methods and apparatuses for reducing time to estimate a position using a satellite positioning system
EP2523497B1 (en) * 2010-01-07 2019-03-20 Nec Corporation Wireless communication system, radio terminal, radio network, wireless communication method and program
RU2432584C2 (ru) * 2010-01-25 2011-10-27 Мстар Семикондактор, Инк. Способ определения координат мобильного приемника спутниковой радионавигационной системы (срнс)
WO2012079616A1 (en) * 2010-12-13 2012-06-21 Telecom Italia S.P.A. Method and system for localizing mobile communications terminals
GB2493200B (en) * 2011-07-28 2013-12-18 Samsung Electronics Co Ltd Pseudorange corrections
US9973884B2 (en) 2012-03-07 2018-05-15 Intel Corporation Device, system and method of controlling access to location sources
CN104205961A (zh) * 2012-03-29 2014-12-10 英特尔公司 移动设备位置估计的设备、系统和方法
US9215562B1 (en) * 2014-07-28 2015-12-15 Intel IP Corporation Method and apparatus for optimized indoor position estimation
US20160299910A1 (en) * 2015-04-09 2016-10-13 Umm Al-Qura University Method and system for querying and visualizing satellite data
DE102017210138A1 (de) * 2017-06-16 2018-12-20 Robert Bosch Gmbh Verfahren und Vorrichtung zum Senden von Korrekturdaten und zum Bestimmen einer hochgenauen Position einer mobilen Einheit
US10564290B1 (en) * 2018-08-10 2020-02-18 Rockwell Collins, Inc. Systems and methods for GNSS processing to determine secondary code phase
CN109507699B (zh) * 2018-12-03 2021-05-28 江苏本能科技有限公司 车用卫星定位校正方法及装置
US11240773B2 (en) * 2018-12-07 2022-02-01 Google Llc Managing doppler and framing impacts in networks
CN112596086A (zh) * 2020-10-20 2021-04-02 陕西航天技术应用研究院有限公司 一种用于低轨通信卫星移动终端的定位方法、装置及系统
CN116819584B (zh) * 2023-08-29 2023-11-14 北京北斗星通导航技术股份有限公司 一种基于北斗卫星的精准定位方法、系统及存储介质
CN117111126B (zh) * 2023-10-20 2023-12-22 成都格理特电子技术有限公司 基于北斗的石化工作人员联合定位方法和装置

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08166436A (ja) * 1994-12-13 1996-06-25 Japan Radio Co Ltd 衛星信号受信機
US6191729B1 (en) * 1996-06-07 2001-02-20 Sextant Avionique Satellite signal receiver with integrity control and exclusion of defective axes
US20020080064A1 (en) * 2000-12-26 2002-06-27 Lg Electronics Inc. System for determining position of mobile communication terminal and method thereof
US6466164B1 (en) * 2001-09-10 2002-10-15 Nokia Mobile Phones Ltd. Method and apparatus for calculating pseudorange for use in ranging receivers
JP2002540435A (ja) * 1999-03-29 2002-11-26 クゥアルコム・インコーポレイテッド アナログモードで作動しているｇｐｓ装備無線装置の位置を見つける方法および装置
US6515618B1 (en) * 2000-11-29 2003-02-04 Trimble Navigation Ltd. Fault detection and exclusion in a positioning system receiver

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2749665B1 (fr) * 1996-06-07 1998-09-11 Sextant Avionique Recepteur de signaux de satellites avec controle d'integrite de calcul de vitesse
US6937187B2 (en) * 2000-11-17 2005-08-30 Global Locate, Inc. Method and apparatus for forming a dynamic model to locate position of a satellite receiver
US6810324B1 (en) * 2002-06-06 2004-10-26 Trimble Navigation, Ltd. Substitution of high quality position measurements with upgraded low quality position measurements
US20040189515A1 (en) * 2003-03-31 2004-09-30 Giovanni Vannucci Hybrid satellite based positioning system
US6839020B2 (en) * 2003-06-02 2005-01-04 Motorola, Inc. Aiding location determinations in satellite positioning system receivers

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08166436A (ja) * 1994-12-13 1996-06-25 Japan Radio Co Ltd 衛星信号受信機
US6191729B1 (en) * 1996-06-07 2001-02-20 Sextant Avionique Satellite signal receiver with integrity control and exclusion of defective axes
JP2002540435A (ja) * 1999-03-29 2002-11-26 クゥアルコム・インコーポレイテッド アナログモードで作動しているｇｐｓ装備無線装置の位置を見つける方法および装置
US6515618B1 (en) * 2000-11-29 2003-02-04 Trimble Navigation Ltd. Fault detection and exclusion in a positioning system receiver
US20020080064A1 (en) * 2000-12-26 2002-06-27 Lg Electronics Inc. System for determining position of mobile communication terminal and method thereof
US6466164B1 (en) * 2001-09-10 2002-10-15 Nokia Mobile Phones Ltd. Method and apparatus for calculating pseudorange for use in ranging receivers

Also Published As

Publication number Publication date
CN1922505A (zh) 2007-02-28
EP1718987A1 (en) 2006-11-08
US7633436B2 (en) 2009-12-15
EP1718987B1 (en) 2009-03-18
WO2005078474A1 (en) 2005-08-25
JP2007523345A (ja) 2007-08-16
DE602005013366D1 (de) 2009-04-30
US20070159384A1 (en) 2007-07-12
CN1922505B (zh) 2010-05-26
WO2005078474A8 (en) 2006-10-19
ATE426178T1 (de) 2009-04-15

Similar Documents

Publication Publication Date Title
JP4896746B2 (ja) 衛星を利用した移動端末の測位
US6839020B2 (en) Aiding location determinations in satellite positioning system receivers
US7450063B2 (en) Method and arrangements relating to satellite-based positioning
KR101843004B1 (ko) 다중 위성 항법 시스템의 시스템 바이어스를 이용한 통합 측위 장치 및 그 방법
US6476762B2 (en) Method for performing positioning and an electronic device
JP5450915B2 (ja) 基準時刻誤差の決定方法及び電子装置
WO2008035139A1 (en) Relative positioning
JP2005524073A (ja) ハイブリッド位置決定システムにおける無線端末のための位置決定
US6771215B2 (en) Determination of the transmission time of a signal part in a positioning system
US8350756B2 (en) Detecting and compensating for erroneous time measurement for mobile device position calculation
JP2013127470A (ja) 受信機を測位する方法、測位システム及び電子装置
CN100354645C (zh) 一种实现定位的方法和一种电子设备
Akopian et al. A fast positioning method without navigation data decoding for assisted GPS receivers
US8085704B2 (en) Method and arrangements relating to satellite-based positioning
EP1229341B1 (en) A method for defining the error of reference time and an electronic device
KR100801347B1 (ko) 이동 단말기의 위성 기반 위치 확인 방법 및 장치

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20080204

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20110314

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20110518

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20111216

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20111221

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20150106

Year of fee payment: 3

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees