【手続補正書】
【提出日】平成17年3月1日(2005.3.1)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
下記を具備する、装置の位置を決定するためのコンピュータプログラムプロダクト：
（ａ）複数の送信器の位置の初期推定値に一部基づいて前記装置の粗位置推定値を受信するコンピューター読み出し可能な命令；
（ｂ）前記送信器のための修正された位置推定値を受信するためのコンピューター読み出し可能命令；
（ｃ）下記を具備する反復アルゴリズムを介して粗位置推定値を修正するためのコンピューター読み出し可能命令：
（１）変数が前記粗位置推定値に等しくなるように、現在のより正確な位置を表す前記変数の値を初期化する；
（２）前記装置が前記現在のより正確な位置に位置し、前記送信器が前記複数の送信器の相対的により正確な位置に位置していたなら、前記装置によって測定されたであろう測定値に等しい擬似距離測定値を決定する；
（３）レシーバー装置が前記現在のより正確な位置に位置しておりかつ前記送信器が前記複数の送信器の位置の初期推定値により示される位置に位置していたなら、前記レシーバー装置自身が位置している場所を計算したであろう位置
【数１】
2004526149.tif 000001
を決定するために、最小２乗を計算する；
（４）幾何学マトリクスＨaおよび幾何学マトリクスＨbを定義する；
（５）最終の幾何学マトリクスＨtを発生する；
（６）最終の幾何学マトリクスから更新ベクトルを計算する；
および
（７）前記現在のより正確な位置の古い値に前記更新されたベクトルを加算することにより前記現在のより正確な位置のための新しい値を得る；
（ｄ）前記装置のための修正された位置推定値を初期化する；
（ｅ）前記装置のための前記修正された位置推定値および前記送信器のための初期位置推定値および修正された位置推定値を計算するためのコンピューター読み出し可能な命令；
（ｆ）前記更新されたベクトルに基づいて、前記装置のための前記修正された位置を更新するためのコンピューター読み出し可能な命令；
（ｇ）前記計算と更新を複数回繰り返すためのコンピューター読み出し可能な命令；および
（ｈ）コードを記憶するように構成されたデータ記憶媒体。
【請求項２】
下記を具備する位置識別装置：
（ａ）複数の送信器の位置の初期推定値に一部基づいてレシーバー装置のための粗位置推定値を受信するように構成された通信ポート；
および
（ｂ）前記通信ポートと動作可能に接続され、下記を実行するように構成されたプロセッサー：
（１）前記装置の位置の粗位置推定値を受信する、前記推定値は複数の送信器の位置の初期推定値を用いて計算された；
（２）前記複数の送信器の相対的により正確な位置を供給する情報を受信する；
そして
（３）下記を具備する反復アルゴリズムを用いて前記粗位置推定値を修正する：
（ａ）変数が前記粗位置推定値に等しくなるように、現在のより正確な位置推定値を表す変数の値を初期化する；
（ｂ）前記装置が前記現在のより正確な位置に位置しておりかつ前記送信器が前記複数の送信器の前記相対的により正確な位置に位置していたなら、前記装置により測定されたであろう、測定値に等しい擬似距離測定値を決定する；
（ｃ）前記レシーバー装置が前記現在のより正確な位置に位置しており、かつ前記送信器が前記複数の送信器の位置の前記初期推定値により示される位置に位置していたなら、前記レシーバー装置自身が位置している場所を計算したであろう位置
【数２】
2004526149.tif 000002
を決定するために最小２乗を計算する；
（ｄ）幾何学マトリクスＨaおよび幾何学マトリクスＨeを定義する；
（ｅ）最終の幾何学マトリクスＨtを発生する；
（ｆ）前記最終幾何学マトリクスから更新ベクトルを計算する；
および
（ｇ）前記現在のより正確な位置の古い値に前記更新されたベクトルを加算することにより、前記現在のより正確な位置のための新しい値を得る。
【請求項３】
さらに下記を具備する請求項２の位置識別装置：
（ａ）前記プロセッサーに動作可能に接続され、複数のタイムインスタンスにおいて前記送信器のための位置推定値を示す複数の情報セットを記憶するように構成されたメモリ、前記送信器のために前記修正された位置推定値は前記情報セットの１つから派生する。
【請求項４】
前記通信ポートは、前記送信器からの信号が測定される特定の時刻を示す情報を受信するように、および前記レシーバー装置のための前記粗位置推定値を派生するために使用されるようにさらに構成される、請求項２の位置識別装置。
【請求項５】
下記を具備する、装置の位置を決定するための方法：
（ａ）前記装置の前記位置の粗位置推定値を受信する、前記推定値は、複数の送信器の前記位置の初期推定値を用いて計算された；
（ｂ）前記複数の送信器の相対的により正確な位置を供給する情報を受信する；
および
（ｃ）下記を具備する反復アルゴリズムを用いて前記粗位置推定値を修正する：
（１）変数が前記粗位置推定値に等しくなるように、現在のより正確な位置推定値を表す前記変数の前記値を初期化する；
（２）前記装置が前記現在のより正確な位置に位置していたならかつ前記送信器が前記複数の送信器の前記相対的により正確な位置に位置していたなら、前記装置によって測定されたであろう測定値に等しい擬似距離測定値を決定する；
（３）前記レシーバー装置が前記現在のより性格な位置に位置していたならかつ前記送信器が前記複数の送信器の前記位置の前記初期推定値により示される前記位置に位置していたなら、前記レシーバー装置自身が位置している場所を計算したであろう位置
【数３】
2004526149.tif 000003
を決定するために最小２乗を計算する；
（４）幾何学マトリクスＨaおよび幾何学マトリクスＨbを定義する；
（５）  最終の幾何学マトリクスＨtを発生する；
（６）  最終の幾何学マトリクスから更新ベクトルを計算する；
および
（７）前記現在のより正確な位置の古い値に前記更新されたベクトルを加算することにより、前記現在のより正確な位置のための新しい値を得る。
【請求項６】
ステップ（ｂ）乃至（ｃ）は、前記更新ベクトルの大きさが特定の閾値内に入るまで反復される、請求項５の方法。
【請求項７】
複数の送信器の前記位置の前記初期推定値は、暦データにもとづく、請求項６の方法。
【請求項８】
前記複数の送信器の相対的により正確な位置を供給する前記情報は、天文暦データである、請求項６の方法。
【請求項９】
前記測定値が前記装置において作られた時間の不確実性を克服するために、前記装置の前記位置推定は、５以上の測定値に基づく、請求項５の方法。
【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【発明の詳細な説明】
【発明の名称】初期粗位置推定値に基づいて改良された位置推定値を供給する方法および装置
【技術分野】
【０００１】
［１００１］この発明は一般に装置の位置に関し、特に初期粗位置推定値に基づいて装置の位置の改良された（すなわち、より正確な）推定値を供給するための技術に関する。
関連技術
［１００２］この出願は、この出願の譲受人に譲渡された、２００１年３月５日に出願された米国仮出願シリアル番号６０／２７３，５７０の優先権を請求する。
【背景技術】
【０００２】
［１００３］位置を見つける共通の手段は、既知の位置における複数のソースから送信された信号が位置を見つける装置内の受信器に到達するのに要求される時間量を決定することである。既知の位置の複数の送信器から信号を供給する１つのシステムは、良く知られている全地球測位システム（ＧＰＳ）である。ＧＰＳシステム内の衛星はＧＰＳマスタプランに従って正確な軌道に載置される。ＧＰＳ衛星の位置は、多数の異なる情報セットにより識別可能である。そのうちのいくつかは他よりもより正確である。
【０００３】
［１００４］ＧＰＳ衛星は「暦」と呼ばれる情報のセットを送信する。暦は「星座」内の衛星の位置に関するより正確でない情報を含む。地上局は、軌道内の変化を観察するためにＧＰＳ衛星を連続的に監視する。衛星位置が測定されると、その情報が中継されて衛星に戻される。衛星は「天文暦」と呼ばれるもう一つの情報セットを送信する。天文暦は衛星軌道のより精度の高いバージョンを含む。各衛星は暦情報を全ての衛星に対して送信するが、天文暦情報を自身に対してのみ送信する。
【０００４】
［１００５］ＧＰＳ受信器は多数の衛星の各々が一時に空のどこに位置するかを示す暦を受信しおよび／または記憶することができる。ＧＰＳ衛星のロケーションのより正確な判定は、天文暦およびこの情報が入手可能な日時に基づいて行なうことができる。
【０００５】
［１００６］暦と天文暦は共に限られた時間量に対して有効である。暦情報は、暦が送信されてからほぼ１週間、約３キロメートルまで正確であると考えられる。天文暦は、ほぼ２時間で約１メートルの精度で、衛星軌道に関する情報を供給する。暦と天文暦両方における誤差は、情報が古くなるにつれ大きくなる。従って、時間的に適当な間隔で更新された情報を受信しない限り、暦と天文暦が古くなるにつれこの情報に基づく衛星の位置はだんだん正確でなくなる。
【０００６】
［１００７］衛星の位置に関する正確な情報なしに、衛星から送信された信号の受信に基づいて決定されたかもしれない、装置の推定された位置は不正確かもしれない。正確な情報は、衛星からまたは代替ソースから（連続的にまたは必要に応じて）更新を受信することにより得ることができる。代替ソースは、無線通信システムにおける基地局または位置決定機器（ＰＤＥ）であってよい。いずれも、ＧＰＳ衛星から必要な情報を受信することができるＧＰＳ受信器を持つことができる。しかしながら、規則的な間隔で情報を得るための装置を設置するには、高価なリソースを消費するであろう。特に、情報を送信し、受信するために電力が必要であり、情報を遠隔ソースから装置に送信するのに帯域が必要である。
【０００７】
［１００８］それゆえ、最小のリソースの消費で装置の位置を正確に決定するための技術の必要性がある。この必要性は、位置決定が時間に対して移動する送信器（例えばＧＰＳ衛星）に基づくとき、およびそのような送信器の正確な位置が、送信器または、衛星およびその位置が決定される装置から離れたソースから更新値を受信するときのみわかる時特に重大である。
【発明の開示】
【課題を解決するための手段】
【０００８】
［１００９］送信器の位置に関する正確でない情報（例えば、ＧＰＳ衛星のための古い暦または古い天文暦）を用いて派生されたかもしれない、初期粗位置推定値に基づいて受信器装置の位置を正確に決定するための技術がここに開示される。開示された装置と方法の一実施の形態において、受信器装置の粗位置推定値に対する訂正は、（１）相対的により正確な送信器の位置を供給する情報、および（２）粗位置推定値を派生するために使用される相対的に正確でない送信器の位置（例えば、古い暦および／または天文暦）を供給する情報の知識に基づいて行なわれる。訂正は、種々のアルゴリズムに基づいて行なうことができ、そのうちのひとつは、以下に詳細に記載する繰り返しアルゴリズムである。
【０００９】
［１０１０］開示された方法および装置の一実施の形態に従って、装置のための粗位置推定値が受信される。粗位置推定値は、ＧＰＳ衛星のような、複数の送信器の位置の初期（正確でない）推定値に基づいて派生されたかもしれない。その送信器のための訂正された（より正確な）位置推定値も受信される。その装置のための訂正された位置推定値が次に（例えば粗位置推定値）に初期化される。次に、その装置のための初期の訂正された位置推定値およびその送信器のための初期の訂正された位置推定値に基づいて更新ベクトルが計算される。次に、その装置のための訂正された位置推定値が更新されたベクトルに基づいて更新される。更新ベクトルの計算および装置の訂正された位置推定値の更新は、装置の実際の位置のますます正確な推定値を得るために何度も（例えば更新ベクトルの大きさが特定のしきい値内に入るまで）繰返すことができる。
【００１０】
［１０１１］開示された方法および装置の一実施の形態に従って、更新ベクトルは以下のステップを実行することにより計算される。（１）装置のための訂正された位置推定値および送信器のための訂正された位置推定値に基づいて測定値のセットを計算する。（２）測定値のセットおよび送信器のための初期位置推定値に基づいて、装置のための中間位置推定値を派生する。（３）送信器のための初期位置推定値および装置のための中間位置推定値に基づいて第１の幾何学マトリクスを決定する。（４）送信器のための訂正された位置推定値および装置のための訂正された位置推定値に基づいて第２の幾何学マトリクスを決定する。および（５）装置のための中間位置推定値および訂正された位置推定値および第１および第２の幾何学マトリクスに基づいて更新ベクトルを計算する。これらのステップは以下に詳細に記載する。
【００１１】
［１０１２］開示された方法および装置はさらに以下にさらに詳細に記載するように、他の観点、実施の形態、および特徴を提供する。
【発明を実施するための最良の形態】
【００１２】
［１０１９］図１はこの発明の種々の観点を実施することができるシステムの簡単化された説明図である。その位置が決定される受信器装置１１０は、アンテナを介して複数の送信器１３０ａ−１３０ｄ（参照番号「１３０」を用いて集合的に呼ぶ）から送信された信号を受信する。受信器装置１１０はさらに位置識別装置１２０と通信する。位置識別装置１２０は、以下にさらに詳細に記載するように、受信器装置の位置を決定する手助けをする。
【００１３】
［１０２０］一実施の形態において、受信器装置１１０は複数の送信器１３０から信号を受信することができる携帯電話である。しかしながら、受信器装置１１０は基準時間に対して受信した信号の到着時刻を決定することができる何らかの装置であってよい。例えば受信器装置１１０は無線モデムを有するコンピュータ端末、スタンドアローンＧＰＳ受信器、地上の送信器からの信号を受信することができる受信器、または他のいずれかの受信器であってよい。送信器１３０は、知られているまたは確定することができる位置を有するいずれかのタイプの送信器であってよい。一実施の形態において、送信器１３０は全地球測位（ＧＰＳ）システムの衛星である。開示した方法と装置の代わりの実施の形態において、送信器１３０は地上の送信器（例えば、無線通信システムの基地局）または、ＧＰＳと地上送信器の組合せであってよい。
【００１４】
［１０２１］受信器装置１１０は受信された信号および受信された信号が送信された送信器１３０の位置を示す情報に基づいてその位置を推定する。いろいろな要因により、受信器装置１１０は、送信器１３０の実際の位置に関する現在のまたは正確な情報を持ってないかもしれない。そのような場合に、その位置の受信器装置１１０により作られた位置推定値は粗であり、所望の精度を有してないかもしれない。正確でない送信器の位置を用いて作られた初期粗位置推定値はその後ここに記載する技術を用いて「訂正」され、受信器装置のより正確な位置推定値を供給することができる。
【００１５】
［１０２２］送信器１３０は受信器装置１１０の位置を決定するために基準点として使用することができる。３つの送信器１３０からの距離を正確に測定することにより、受信器装置１１０はその位置を「三角測量」する。受信器装置１１０は、信号が送信器１３０から受信器装置１１０に移動するのに必要な時間を測定することにより距離を決定する。信号が送信器１３０から送信される時刻を知り、信号が（その内部クロックに基づいて）受信器装置１１０により受信される時刻を観察することにより、信号の移動時間を決定することができる。しかしながら、送信と受信の間の正確な時間量は一般的に、例えば、送信器１３０と受信装置１１０における２つのクロックにおけるオフセットにより決定することができない。従って、一般的に「基準」時間と信号が受信された時間との間の差に基づいて「擬似距離」が計算される。基準時間が作られるすべての擬似距離測定値に共通である限り、基準時間はいずれの時間であってもよい（または、擬似距離測定値は使用される基準時間内の差を補償するために調節することができる）。
【００１６】
［１０２３］図１Ａにおいて、送信器１３０ａ乃至１３０ｄは正確でない（例えば最新でない）情報に基づいて受信器装置１１０により推定された位置に示される。これらの位置は、送信器１３０の実際の位置と異なるかもしれない。送信器１３０がＧＰＳ衛星のような衛星である場合に、衛星の位置は、暦および天文暦により識別可能である。
【００１７】
［１０２４］現在の天文暦情報および日時が利用可能なら、ＧＰＳ衛星の位置の正確な決定を行なうことができる。しかしながら、天文暦情報は限られた時間量（例えば２時間）に対して有効である。現在の暦情報と日時が利用可能なら、ＧＰＳ衛星位置の正確でない決定を行なうことができる。しかしながら、暦もまた、限られた時間量（例えば１週間）に対して有効である。
【００１８】
［１０２５］発明の観点は、初期粗位置推定値に基づいて受信器装置の位置を正確に決定するための技術を提供する。この初期粗位置推定値は、送信器の位置に関するより正確でない情報（例えばＧＰＳ衛星のための古い暦または古い天文暦）を用いて得られたかもしれない。初期推定値は、完全に時代遅れの暦または天文暦情報に基づいて作られたかもしれないと理解されるべきである。
【００１９】
［１０２６］送信器の擬似距離および位置は受信器装置の位置を決定するために使用されるので、送信器の位置における誤差または不正確は受信器装置のための位置推定値における誤差または不正確になる。それゆえ、「粗」位置推定値は、送信器１３０までの擬似距離に基づいて受信器装置のために派生される。送信器の位置は相対的に悪い精度で知ることができる。
【００２０】
［１０２７］開示された方法および装置の一実施の形態に従って、受信器装置の粗位置推定値に対する訂正は、以下に基づいて行なわれる。（１）送信器１３０の位置に関するより正確な情報および（２）粗位置推定値を派生するために使用される送信器１３０の位置の推定値における誤差量に関する情報。訂正を行なうための技術は以下にさらに詳細に記載される。一実施の形態において、訂正は、受信器装置１１０に対して遠隔に位置するかもしれない位置識別装置１２０において行なわれる。しかしながら、開示された方法および装置の他の実施の形態において、位置識別装置１２０は受信器装置１１０と同一場所に配置することができる。
【００２１】
［１０２８］図２は、開示された方法および装置の一実施の形態に従う受信器装置１１０と位置識別装置１２０の簡単化されたブロック図である。受信器装置１１０は、携帯電話や無線モデムを有したコンピュータのような無線通信システムにおける遠隔端末の構成要素であってよい。あるいは、受信器装置１１０は、スタンドアローンＧＰＳ受信器のようなスタンドアローン位置決定ユニットであってもよい。
【００２２】
［１０２９］位置識別装置１２０は、無線通信システムにおける基地局トランシーバサブシステム（ＢＴＳ）、基地局コントローラ（ＢＳＣ）、携帯電話システム向け交換機（ＭＳＣ）の一構成要素であってもよいし、あるいは、無線通信システムのこれらの要素の１つに接続する装置であってもよい。例えば、位置識別装置１２０は無線通信システムにおける位置決定機器（ＰＤＥ）であってもよい。あるいは、位置識別装置１２０は遠隔端末あるいは受信器装置１１０を含むその他の装置の一構成要素であり、同一場所に配置してもよい。あるいは、位置識別装置１２０はスタンドアローンな構成要素であってもよい。
【００２３】
［１０３０］簡単のために、受信器装置１１０と位置識別装置１２０は直接接続されるものとして図２に示される。しかしながら、１つ以上のさらなるエレメント（例えば、ＢＴＳ、ＢＳＣ、等）を受信器装置１１０と位置識別装置１２０との間に接続することができる。位置識別装置１２０と受信器装置１１０もまた同一場所に配置してもよい。このことは、測定が行なわれた時送信器１３０の位置についての正確な情報が入手できないが、いつかこれから先入手可能であるシステムにあてはまる。受信器装置１１０のための粗位置推定値は、擬似距離（または実際の距離）を測定すると直ちに派生することができる。粗位置推定値は、送信器１３０のより正確な位置に関する情報が入手できるようになるまで記憶しておくことができる。その後で受信器装置のためのより正確な位置推定値を派生することができる。
【００２４】
［１０３１］図２に示す実施の形態において、受信器装置１１０は、受信器２１０、処理回路２１２、メモリ２１４、および通信ポート２１６を含み、さらにアンテナ２０８に接続する。アンテナ２０８は（図１Ａに示す）送信器１３０からの信号を受信し、受信した信号を受信器２１０に接続する。受信器２１０は、受信した信号を処理して、位置決定のために使用される情報を派生させるのに必要な回路（例えば、無線周波数（ＲＦ）処理回路および／または他の受信回路）を含む。そのような情報はタイミング情報等を含んでいても良い。情報は処理回路２１２に供給され、処理回路２１２は位置決定計算を実行する。そのような計算は、送信器１３０までの擬似距離を派生させるための計算を含んでいても良い。あるいはまたはさらに、その計算は、受信器装置１１０の初期粗推定値を派生させるのに必要な計算を含んでいても良い。粗位置推定値は、送信器１３０の位置推定値における大きな誤差により相対的に大きな誤差を持つかもしれない。
【００２５】
［１０３２］受信器装置１１０は、通信ポート２１６を介して、計算結果および他の関連情報を位置識別装置１２０に供給する。関連情報は、受信された信号が測定された時刻、計算を行うために使用される情報（例えば、送信器１３０の位置を推定するために使用される特定の暦）および恐らく他の情報を含んでいてもよい。通信ポート２１６は、位置識別装置１２０にインターフェースを供給する。無線通信システムの遠隔端末内に配置された受信器装置の場合、通信ポート２１６は無線リンクを介して通信をサポートすることができる。
【００２６】
［１０３３］図２に示す実施の形態において、位置識別装置１２０は通信ポート２２０、処理回路２２２、およびメモリ２２４を含む。通信ポート２２０は受信器装置１１０から計算結果と関連情報を受信し、この情報を処理回路２２２に転送する。処理回路２２２は計算結果と関連情報を受信し、さらに、その計算を行なうために受信器装置１１０により使用された「正確でない情報」のコピーを受信する。正確でない情報は、暦（および暦が有効であった時間）、天文暦（および天文暦が有効であった時間）および／または受信器装置１１０が送信器１３０の位置を推定するために使用したかもしれない何らかの他の情報であってよい。
【００２７】
［１０３４］正確でない情報は、例えば関連情報の一部として計算結果とともに受信器装置１１０により供給されるようにしてもよい。ある実施の形態において、正確でない情報は明白に供給されなくてもよく、そして使用された情報を推定するために他のインジケータを使用してもよい。例えば、暦または天文暦が有効である時間、受信された信号が受信器装置１１０において測定された時間、あるいは（関連情報の一部として送ってもよい）その他の情報は、位置識別装置１２０が、粗位置推定値を派生させるために使用された送信器の位置推定値を正しく推定することを可能にさせるのに十分であるかもしれない。その他の実施の形態において、位置識別装置１２０は、粗位置推定値を派生するために使用される情報を受信器装置１１０に送信することを担うことができる。これらの実施の形態の場合、受信器装置１１０が位置識別装置１２０にどの情報を使用したかを知らせる必要は無いであろう。
【００２８】
［１０３５］ある実施の形態において、処理回路２１２は別個の通信ポートの必要性なくして位置識別装置１２０と直接通信することができる。これは、位置識別装置１２０と受信器装置１１０が同一場所に配置される場合にあてはまる、あるいは位置識別装置と受信器装置が同一場所に配置されていない、ある他の場合でも可能である。
【００２９】
［１０３６］一実施の形態において、位置識別装置１２０は受信器装置１１０から計算結果を受信する前に正確でない情報を受信する。いくつかの実施の形態において、受信器装置１１０により使用される正確でない情報は受信器装置１１０以外のソース（例えば、簡単のために図２に示していない通信リンクを介して）によって位置識別装置１２０に供給してもよい。例えば、正確でない情報は、送信器１３０から位置識別装置１２０により直接受信してもよい。あるいは、正確でない情報は、無線通信システム内の基地局の構成要素のように、図１Ａに示す構成要素のいずれからも区別されるソースから位置識別装置１２０によって受信するようにしてもよい。
【００３０】
［１０３７］一実施の形態において、位置識別装置１２０は複数のセットの情報（例えばいくつかのバージョンの暦）を持つようにしてもよい。そのうちのいずれか一つが位置決定計算を実行するために受信器装置１１０により使用されたかもしれない。この場合、受信器装置１１０は、位置識別装置１２０により入手可能な情報の中からどの情報が、位置決定計算を実行するために受信器装置１１０により特に使用されたかを示すために、さらなる情報を位置識別装置１２０に供給する必要があるかもしれない。
【００３１】
［１０３８］図３は受信器装置により以前に派生された粗位置推定値に基づいて受信器装置１１０のためのより正確な位置推定値を派生させるために位置識別装置１２０により実行されるプロセスを図解するフロー図である。ステップ３１２において、位置識別装置１２０は最初に受信器装置１１０から位置決定計算の結果（例えば、粗位置推定値
【００３２】
【数４】
2004526149.tif 000004
）を受信する。ステップ３１４において、位置識別装置１２０は、粗位置推定値を派生するために使用された送信器１３０の初期位置推定値
【００３３】
【数５】
2004526149.tif 000005
（すなわち、正確でない送信器の位置推定値）を示す情報も受信する。この情報は、粗位置推定値を派生するために使用された暦または天文暦のバージョンを示すことができる。受信器装置により決定される、送信器のための初期位置推定値は正確でなく、位置決定計算を行うとき、送信器が位置している場所を受信器装置１１０が推測する位置に相当する。一実施の形態において、ステップ３１４においても、位置識別装置１２０はさらに擬似距離測定が行なわれる時間を受信する。
【００３４】
［１０３９］ステップ３１６において、位置識別装置１２０は送信器１３０のロケーションのより正確な推定値
【００３５】
【数６】
2004526149.tif 000006
を決定する。これらのより正確な送信器位置推定値は、擬似距離測定が受信器装置１１０により行なわれた時間に対してより正確な暦および／または天文暦に基づいて作ることができる。一実施の形態において、位置識別装置１２０はＧＰＳ星座の衛星により送信された暦および天文暦の通信記録を維持する。そのような通信記録は、以下にさらに詳細に記載するように、受信器装置１１０から受信した粗位置推定値を訂正するために最も正確な暦および天文暦を位置識別装置１２０が使用可能にする。
【００３６】
［１０４０］位置識別装置１２０が送信器１３０のための初期のより正確でない位置推定値およびより正確な位置推定値に関する情報を有すると、次に、位置識別装置１２０は、ステップ３１８において、以下に記載するアルゴリズムに従って受信器装置１１０のための粗位置推定値を訂正する。
【００３７】
［１０４１］図１Ｂは、受信器装置１１０に入手可能な情報に基づいて送信器１３０の位置の誤差を図解する図である。図１Ｂは、送信器の実際の位置のより正確な推定値を表す位置における送信器１３２ａ乃至１３２ｄを示す。図１Ｂはまた送信器の位置のより正確でない推定値を表す位置における送信器１３０ａ乃至１３０ｄも示す。（破線で示される）送信器１３２と（実線で示される）送信器１３０は同じ送信器のための２つの位置推定値に対応する。しかしながら、各送信器の位置は、送信器の位置がより正確な情報またはより正確でない情報を用いて（すなわち、擬似距離測定が行なわれた時間に対して正確であるまたは正確でない暦および天文暦を用いて）決定されたかどうかに応じて異なることに気付く。送信器１３０の位置は、受信器装置１１０のための粗位置推定値を派生するために位置決定計算を実行するために受信器装置１１０により使用された。
【００３８】
［１０４２］位置識別装置１２０は受信器装置のためのより正確な位置推定値を供給するために受信器装置１１０からの粗位置推定値に関して訂正を実行する。粗位置推定値に関する訂正は、「線形化」アルゴリズム、「繰り返し」アルゴリズムあるいはその他を含む種々のアルゴリズムに基づいて実行することができる。線形化アルゴリズムは、この出願の譲受人に譲渡され、２００１年１月３０日に出願された米国特許出願シリアル番号第０９／７７３，２０７号（発明の名称：「粗位置推定値を用いて位置を決定するための方法および装置」）に詳細に記載されている。
【００３９】
［１０４３］繰り返しアルゴリズムは以下に基づいて受信器装置１１０のための粗位置推定値に関して訂正を行なうために使用することができる。（１）受信器装置１１０のための粗位置推定値を派生させるために使用される送信器１３０（例えばＧＰＳ衛星）の位置の正確でない推定値の知識、および（２）擬似距離測定が行なわれた時刻における送信器１３０の実際の位置のより正確な推定値の知識。より正確な送信器の位置推定値は、擬似距離測定が行なわれた特定の時刻の知識に基づいて決定することができる。
【００４０】
［１０４４］送信器１３０の位置のより正確な推定値が入手できないとき位置または時刻において粗位置推定値が作られることが理解されなければならない。いつおよびどこで、送信器１３０の位置のより正確な推定値が入手可能になるかに応じて、訂正は一般的に後の時刻または別の位置またはその両方で行なわれるであろう。
【００４１】
［１０４５］開示された方法および装置の一実施の形態において、位置標定解（すなわち、受信器装置１１０のための計算された位置）および送信器位置推定値のセットとの間に１対１のマッピングがあると、仮定される。この仮定は、わずか４つの擬似距離測定値がある場合に容易に満足される。
【００４２】
［１０４６］明確にするために、繰り返しアルゴリズムは、送信器１３０がＧＰＳ衛星である場合に特に記載される。しかしながら、繰り返しアルゴリズムは、また、正確でない送信器位置情報が特定の時刻または場所において入手可能であり、より正確な送信器位置情報が多少、後の時刻および／または他の位置において入手可能である何らかのタイプの送信器を用いて使用することができる。
【００４３】
［１０４７］暦情報およびクロックバイアスに基づいて派生されたｉ番目の衛星のための位置
【００４４】
【数７】
2004526149.tif 000007
は以下のように表すことができる。
【００４５】
【数８】
2004526149.tif 000008
クロックバイアスは、受信器装置１１０内の局部クロックにより示される時刻と送信器１３０内のクロックにより示される時刻との差として定義される。位置
【００４６】
【数９】
2004526149.tif 000009
は、位置
【００４７】
【数１０】
2004526149.tif 000010
を派生させるために受信器装置１１０により使用される暦の不正確さによる相対的に大きな誤差を一般的に持つであろう。
［１０４８］天文暦情報およびクロックバイアスに基づいてｉ番目の衛星のために計算された位置
【００４８】
【数１１】
2004526149.tif 000011
は以下のように表すことができる。
【００４９】
【数１２】
2004526149.tif 000012
位置
【００５０】
【数１３】
2004526149.tif 000013
はより正確な天文暦情報の使用により相対的に小さな誤差を一般的に持つであろう。相対的に正確でない暦情報を用いて受信器装置１１０により派生される粗位置推定値
【００５１】
【数１４】
2004526149.tif 000014
は以下のように表すことができる。
【００５２】
【数１５】
2004526149.tif 000015
相対的により正確な送信器位置推定値を用いて受信器装置１１０により派生される受信器装置１１０のより正確な位置推定値
【００５３】
【数１６】
2004526149.tif 000016
は以下のように表すことができる。
【００５４】
【数１７】
2004526149.tif 000017
最初に、より正確な情報は入手することができない。それゆえ、より正確な位置推定値
【００５５】
【数１８】
2004526149.tif 000018
の値は
【００５６】
【数１９】
2004526149.tif 000019
に設定することができる。位置推定値
【００５７】
【数２０】
2004526149.tif 000020
はその後、より多くの情報が得られるので更新される。
［１０４９］各衛星に相関する距離推定値は以下の式に基づいて得ることができる。
【００５８】
【数２１】
2004526149.tif 000021
これらの距離推定値は必須的に、より正確な情報により示される位置に位置する送信器から信号を受信する、位置
【００５９】
【数２２】
2004526149.tif 000022
に位置する受信装置により決定されたであろう擬似距離測定値である。式（１）の測定値、ｍｉは、クロックバイアスを考慮に入れて、受信器装置１１０のより正確な位置推定値
【００６０】
【数２３】
2004526149.tif 000023
と送信器１３０のより正確な位置推定値
【００６１】
【数２４】
2004526149.tif 000024
との間の距離である。各衛星に対して計算されたｍｉの距離値に基づいて、受信器装置１１０が位置
【００６２】
【数２５】
2004526149.tif 000025
に位置し、送信器１３０が送信器１３０のための正確でない位置推定値
【００６３】
【数２６】
2004526149.tif 000026
によって示される位置（すなわち、暦情報により示される位置）に位置すると仮定したならば受信器装置１１０が計算したであろう位置
【００６４】
【数２７】
2004526149.tif 000027
を計算するために最小平均二乗（ＬＭＳ）計算を実行することができる。位置推定値
【００６５】
【数２８】
2004526149.tif 000028
および
【００６６】
【数２９】
2004526149.tif 000029
と
【００６７】
【数３０】
2004526149.tif 000030
との間の部分派生式の知識に基づいて、更新ベクトル
【００６８】
【数３１】
2004526149.tif 000031
を決定することができる。
［１０５０］正確でない送信器の位置推定値
【００６９】
【数３２】
2004526149.tif 000032
（例えば暦から派生された衛星の位置）および以前のＬＭＳ計算から派生された対応する位置
【００７０】
【数３３】
2004526149.tif 000033
に対する幾何学マトリクスＨａは以下のように定義することができる。
【００７１】
【数３４】
2004526149.tif 000034
同様に、より正確な送信位置推定値
【００７２】
【数３５】
2004526149.tif 000035
（例えば天文暦から派生された衛星の位置）およびより正確な受信器装置位置推定値
【００７３】
【数３６】
2004526149.tif 000036
は以下のように定義することができる。
【００７４】
【数３７】
2004526149.tif 000037
［１０５１］最終的な幾何学マトリクスＨｔは以下のように表すことができる。
【００７５】
【数３８】
2004526149.tif 000038
但し、Ｒは測定値の分散行列である、すなわち
【００７６】
【数３９】
2004526149.tif 000039
但しＥ［］はブラケット内の量の統計的期待値を計算することを表し、
【００７７】
【数４０】
2004526149.tif 000040
は、式（１）において、決定された測定値のベクトルである。従って、更新ベクトル
【００７８】
【数４１】
2004526149.tif 000041
は以下のように表すことができる。
【００７９】
【数４２】
2004526149.tif 000042
［１０５２］より正確な送信器位置推定値
【００８０】
【数４３】
2004526149.tif 000043
（例えば、天文暦から派生された衛星位置）を用いて得た受信器装置の位置のより正確な推定値
【００８１】
【数４４】
2004526149.tif 000044
は以下のように表すことができる。
【００８２】
【数４５】
2004526149.tif 000045
式（１）乃至（６）は受信器装置１１０のためのますます正確な位置推定値を派生するために何度も反復的に適用することができる。
［１０５３］図４はこの発明の実施形態に従って、繰り返しアルゴリズムのための実行された処理のフロー図である。送信器１３０の正確でない位置推定値
【００８３】
【数４６】
2004526149.tif 000046
に基づいて受信器装置１１０により派生された粗位置推定値
【００８４】
【数４７】
2004526149.tif 000047
は（図３に定義されるように）最初に受信される。ステップ４１２において、受信器装置１１０のためのより正確な位置推定値
【００８５】
【数４８】
2004526149.tif 000048
は粗位置推定値
【００８６】
【数４９】
2004526149.tif 000049
に初期化される。次に、ステップ４１４において擬似距離測定値ｍｉが決定される。これらの測定値は位置推定値
【００８７】
【数５０】
2004526149.tif 000050
が解であったなら得られたであろう測定値である。測定値ｍｉは式（１）に基づいて派生することができ、粗位置推定値
【００８８】
【数５１】
2004526149.tif 000051
を派生するために以前に使用されたすべての送信器のために決定することができる。
［１０５４］ステップ４１６において、ＬＭＳ計算はステップ４１４において、決定された測定値ｍｉおよび正確でない送信器位置推定値
【００８９】
【数５２】
2004526149.tif 000052
に基づいて実行される。ＬＭＳ計算は技術的に知られた方法で達成することができる。ＬＭＳ計算から結果として得られる解は
【００９０】
【数５３】
2004526149.tif 000053
として示され、これは、位置
【００９１】
【数５４】
2004526149.tif 000054
にあったなら、受信器装置１１０が得たであろう解である。
［１０５５］ステップ４１８において、次に、式（２）、（３）および（４）に従って，幾何学マトリクスＨｔが計算される。ステップ４２０において、更新ベクトル
【００９２】
【数５５】
2004526149.tif 000055
は、幾何学マトリクスＨｔおよび位置推定値の差
【００９３】
【数５６】
2004526149.tif 000056
および式（５）に従って決定することができる。ステップ４２２において、次に、より正確な位置推定値
【００９４】
【数５７】
2004526149.tif 000057
は更新ベクトル
【００９５】
【数５８】
2004526149.tif 000058
を用いて更新される。
［１０５６］ステップ４２４において、更新ベクトル
【００９６】
【数５９】
2004526149.tif 000059
が特定のしきい値（例えば、
【００９７】
【数６０】
2004526149.tif 000060
）未満か否か判断される。更新ベクトル
【００９８】
【数６１】
2004526149.tif 000061
の大きさがしきい値未満でないなら、プロセスはステップ４１４に戻り、送信器の測定値ｍｉは最近更新された受信器のためのより正確な位置推定値
【００９９】
【数６２】
2004526149.tif 000062
に基づいて再び決定される。ステップ４１４乃至４２４は更新ベクトル
【０１００】
【数６３】
2004526149.tif 000063
の大きさがしきい値未満であるまで必要な回数反復される。しきい値が満足されると、位置推定値
【０１０１】
【数６４】
2004526149.tif 000064
は受信器装置のためのより正確な位置推定値として戻される。次に処理は終了する。
［１０５７］ここに記載する技術は、送信器のための正確でない位置推定値に基づいて派生されたかもしれない、受信器装置のための粗位置推定値に基づいて受信器装置の位置を正確に決定することを可能にする。このようにして、送信器の位置を記述している情報（例えば暦）は、まれに受信器装置にロードすることができ、これは、リソースの保存を可能にすることができる。
【０１０２】
［１０５８］ここに記載した技術はまた、アナログモードで動作する受信器装置の位置を決定するために有利に使用することができる。これは、十分な精度で知られていない測定値のタイムスタンプ（例えば５ｍｓｅｃ以上のタイムスタンプエラー）で特徴づけられる。
【０１０３】
［１０５９］この発明の観点に従って、５以上の送信器（例えばＧＰＳ衛星）のための測定値は、ｘ、ｙ、ｚおよびクロックバイアス並びにタイムスタンプエラーを解くために使用することができる。上述した訂正アルゴリズムは受信器装置により報告されるタイムスタンプが正確でない場合に適用可能である。
【０１０４】
［１０６０］上述した繰り返しアルゴリズムは５つ以上の測定値に基づいて位置決定計算に対して実施することができる。式（１）乃至（６）において、上述した計算および図４に示したフロー図は上述した方法で適用することができる。しかしながら、４つの測定値に対して通常使用されるＬＭＳアルゴリズムの代わりに、解（すなわち、受信器装置のための位置推定値）を派生させるために動作するために修正されたＬＭＳアルゴリズムが実施される。式（２）および（３）にそれぞれ示す幾何学マトリクスＨａおよびＨｅのための４次元の代わりに、これらのマトリクスは使用される測定値の数に対応する次元を有するであろう。
【０１０５】
［１０６１］５つの測定値が位置決定計算に対して使用されるとき、幾何学マトリクスＨａ，５は以下のように表すことができる。
【０１０６】
【数６５】
2004526149.tif 000065
上述のマトリクスＨａ，５の最初の４つの列は式（２）に示したのと同じ列であり、５番目列は、タイムスタンプに対する各送信器の擬似距離測定値の部分導関数
【０１０７】
【数６６】
2004526149.tif 000066
から構成される。このｉ番目の送信器に対する部分導関数は、以下のように表すことができる。
【０１０８】
【数６７】
2004526149.tif 000067
［１０６２］同様に、５つの測定値に対する幾何学マトリクスは以下のように表すことができる。
【０１０９】
【数６８】
2004526149.tif 000068
この場合もやはり、上記マトリクスＨｅ，5の最初の４つの列は式（３）の列と同じであり、５番目の列はタイムスタンプに対する各送信器の擬似距離測定値の部分導関数
【０１１０】
【数６９】
2004526149.tif 000069
から構成される。ｉ番目の送信器の部分導関数は以下のように表すことができる。
【０１１１】
【数７０】
2004526149.tif 000070
［１０６３］従って、測定時間に対する各衛星の座標の部分導関数（例えば、
【０１１２】
【数７１】
2004526149.tif 000071
）は暦または天文暦から派生することができる。これは、部分導関数が時間の関数として衛星の軌道を記載したためである。上の式から、各繰り返しにおいて、衛星の位置は測定時間ｔの新しい値に基づいて再計算される。測定時間ｔは、値が各繰り返しにおいて再評価される未知数の１つである測定時間誤差から得ることができる。
【０１１３】
［１０６４］アナログモードのための位置決定計算は、この出願の譲受人に譲渡された、１９９９年３月２９日に出願された米国特許出願シリアル番号第０９／２８０，３３７（発明の名称：「アナログモードで動作するＧＰＳを備えた無線装置の位置を見つけるための方法および装置」(METHOD AND APPARATUS FOR LOCATING GPS EQUIPPED WIRELESS DEVICES OPERATING IN ANALOG MODE）にさらに詳細に記載されている。
【０１１４】
［１０６５］図２において、処理回路２１２および２２２はそれぞれメモリ２１４および２２４内に記憶されたプログラム命令を実行することにより、それぞれ初期粗位置推定値および受信器装置１１０のためのより正確な位置推定値を派生させることができる。処理回路２１２および２１４は各々デジタルシグナルプロセッサ（ＤＳＰ）、特定用途向け集積回路（ＡＳＩＣ）、プロセッサ、マイクロプロセッサ、コントローラ、マイクロコントローラ、フィールドプログラマブルゲートアレイ（ＦＰＧＡ）、プログラマブル論理装置、他の電子装置、またはここに記載した機能を実行するように設計されたこれらの組合せのいずれかとして実現してもよい。処理回路２１２および２２２はプログラム命令およびデータを記憶するために使用されるメモリ２１２または２２４を含んでいても良い。
【０１１５】
［１０６６］開示した方法および装置は、当業者がこの発明を製作または使用可能にするために提供される。これらの実施の形態に対する種々の変形例は当業者には容易に明白であり、ここに定義される包括的原理は、この発明の精神と範囲を逸脱することなく他の実施の形態に適用可能である。従って、この発明は、ここに示した実施の形態に限定されることを意図したものではなく、開示された原理および新規な特徴に一致する最も広い範囲が許容されるべきである。
【図面の簡単な説明】
【０１１６】
【図１Ａ】［１０１４］図１Ａはこの発明の種々の観点を実施することができるシステムの簡単化された説明図である。
【図１Ｂ】［１０１５］図１Ｂは、受信器装置に入手可能な情報に基づいて送信器のロケーションにおける誤差を図解する図である。
【図２】［１０１６］図２はこの発明の実施の形態に従う受信器装置および位置識別装置の簡単化されたブロック図である。
【図３】［１０１７］図３は、受信器装置により以前に派生された粗位置推定値に基づいて受信器装置のためのより正確な位置推定値を派生するために、位置識別装置により実行されるプロセスを説明するフロー図である。
【図４】［１０１８］図４はこの発明の実施の形態に従う繰り返しアルゴリズムのために実行される処理のフロー図である。