(57)【要約】ＲＦの能動的か、受動的か、又は支援された位置確認アプローチのための信号処理技法である。該信号処理技法は、フィルタか又は推定器、例えば、非線形フィルタと共に、追加的な測定値としてマルチパス信号を利用する。前記フィルタは、インダイレクト及びダイレクトパス測定値か又は任意の他の利用可能な信号を使用して、観測可能なインダイレクトパスのパラメトリックモデルを構築する。１つか又は複数のダイレクトパス測定値が、（例えば、障害のせいで）その後、失われる場合には、インダイレクトパス測定値を用いて、人間か又は対象となる物体の位置の推定をフィルタが維持する。
2008527394.tif 000002
【特許請求の範囲】
【請求項１】
マルチパス無線信号を用いて位置決めを行うための方法であって、
初期位置を決定し、
少なくとも１つの位置において受信した少なくとも１つのマルチパス無線信号における少なくとも１つの特性に対する該少なくとも１つの位置に関する情報を生成し、
後続する無線マルチパス信号を受信し、及び、
前記後続する信号と前記生成された情報とに基づいて、位置を確立する
といった各ステップを含む、方法。
【請求項２】
前記初期位置を決定することが、
タグに、複数の位置において、少なくとも１つの送信された信号を受信させ、
各位置において前記タグによって受信された前記少なくとも１つの信号の大きさを決定し、
前記初期位置において前記タグによって受信された少なくとも１つの信号の大きさを決定し、及び、
前記初期位置を推定するために、前記決定された信号の大きさと、統計的なモデルとを用いる
ことを含むことからなる、請求項１に記載の方法。
【請求項３】
前記統計的なモデルが、前記タグと、各々の複数の信号送信器との間の信号間隔に関連付けられた信号伝搬特性を示すデータを含むことからなる、請求項２に記載の方法。
【請求項４】
前記確立するステップが、複数の候補位置を計算して、そこから最もふさわしい位置を選択することを含むことからなる、請求項１に記載の方法。
【請求項５】
前記最もふさわしい位置が、統計的な試験を用いて決定されることからなる、請求項４に記載の方法。
【請求項６】
前記情報を生成することが、内部環境における少なくとも部分的なジオメトリと、そこの中のマルチパス信号のジオメトリとを、特徴付けることを含み、該環境内における無線源から反射点までのインダイレクトパスセグメント長を、前記マルチパスのジオメトリが、含むことからなる、請求項１に記載の方法。
【請求項７】
前記インダイレクトパスセグメント長が、パス遅延測定値により表されることからなる、請求項６に記載の方法。
【請求項８】
前記少なくとも１つの特性が、パス遅延を含むことからなる、請求項１に記載の方法。
【請求項９】
前記生成された情報は、パス遅延を、前記環境内における位置に関連付けることからなる、請求項８に記載の方法。
【請求項１０】
前記少なくとも１つの特性が、パスセグメント長を含むことからなる、請求項１に記載の方法。
【請求項１１】
前記生成された情報は、パスセグメント長を、前記環境内における位置に関連付けることからなる、請求項１０に記載の方法。
【請求項１２】
前記無線マルチパス信号が、全地球測位システム（ＧＰＳ）送信器によって生成されることからなる、請求項１に記載の方法。
【請求項１３】
前記無線マルチパス信号が、送信器によって生成されることからなる、請求項１に記載の方法。
【請求項１４】
前記初期位置を決定することが、初期位置のダイレクトパス測定値を生成することを含むことからなる、請求項１に記載の方法。
【請求項１５】
前記ダイレクトパス測定値が、受信信号強度の測定値であることからなる、請求項１４に記載の方法。
【請求項１６】
前記初期位置を決定することが、該初期位置の事前知識を利用することを含むことからなる、請求項１に記載の方法。
【請求項１７】
関連付けられる位置が、前記初期位置である、請求項１に記載の方法。
【請求項１８】
前記生成された情報が、少なくとも１つのパラメータを含むことからなる、請求項１に記載の方法。
【請求項１９】
マルチパス無線信号を用いて位置決めを行うための装置であって、
少なくとも１つの位置において受信した少なくとも１つのマルチパス無線信号における少なくとも１つの特性に対する該少なくとも１つの位置に関する情報を含むメモリと、
後続する無線マルチパス信号を受信するための受信器と、
前記後続する信号と前記生成された情報とに基づいて、位置を決定するための回路構成
とを備える、装置。
【請求項２０】
前記情報は、内部環境における少なくとも部分的なジオメトリと、そこの中のマルチパス信号のジオメトリとを特徴付け、前記マルチパスのジオメトリは、該環境内における無線ソースから反射点までのインダイレクトパスセグメント長を含むことからなる、請求項１９に記載の装置。
【請求項２１】
前記インダイレクトパスセグメント長が、パス遅延測定値により表される、請求項２０に記載の装置。
【請求項２２】
前記少なくとも１つの特性は、パス遅延を含む、請求項１９に記載の装置。
【請求項２３】
前記情報が、パス遅延を、前記環境内における位置に関連付ける、請求項２２に記載の装置。
【請求項２４】
前記少なくとも１つの特性が、パスセグメント長を含む、請求項１９に記載の装置。
【請求項２５】
前記情報が、パスセグメント長を、前記環境内における位置に関連付ける、請求項２４に記載の装置。
【請求項２６】
前記無線マルチパス信号が、全地球測位システム（ＧＰＳ）送信器によって生成される、請求項１９に記載の装置。
【請求項２７】
前記無線マルチパス信号が、送信器によって生成される、請求項１９に記載の装置。
【請求項２８】
関連付けられる位置が、前記初期位置である、請求項１９に記載の装置。
【請求項２９】
前記情報が、少なくとも１つのパラメータを含む、請求項１９に記載の装置。
【発明の詳細な説明】
【技術分野】
【０００１】
発明の分野
本発明は、一般に、ナビゲーション及び位置決め（測位）を行うためのシステム及び方法に関し、特に、マルチパス信号を用いてナビゲーション及び位置決めを行うためのシステム及び方法に関する。
【背景技術】
【０００２】
関連事例の相互参照
本出願は、2005年1月19日にファイリングされた米国仮特許出願第６０／６４５，３９０号の利益を請求し、該米国仮特許出願第６０／６４５，３９０号は、その全体があたかも記載されているかのごとく参照によって本明細書内において組み込まれる。本出願はまた、2004年2月17日に発行された米国特許第６，６９３，５９２号の開示も、その全体が記載されているかのごとく本明細書内において組み込む。
【０００３】
アーバンキャニオンか、又は建築物などの囲まれた構造のような、散乱した環境内において移動している人々か又は存在（アセット：assets）の位置決めを行うために無線周波数（ＲＦ）信号を用いる幾つかの試みが存在している。一般的には、これら先行技術の試みは、固定された基準局か又はビーコンを、受信する幾つかの形態か又は追跡される物体上の応答タグと共に、典型的には利用する。この装置と共に、前記基準局からのＲＦ信号強度（ＲＳＳ：RF signal strength）か、到来時間（ＴＯＡ：time of arrival）か、到来時間差（ＴＤＯＡ：time difference of arrival）か、又は到来角度（ＡＯＡ：angle of arrival）の事象（現象）を用いた三角測量の形態が、前記タグの位置を決定するために用いられている。
【０００４】
大部分のこれら先行技術のアプローチのうちの１つの欠陥は、散乱した環境内において生成されたＲＦマルチパス信号によって、それらが悪い方向に影響を受けることである。前記基準ビーコンと前記タグとの間のダイレクトパス（直接経路）が、確実に一貫して測定されることが不可能な時には、マルチパス信号は、結果として、ＲＳＳ、ＴＯＡ、ＴＤＯＡ、及びＡＯＡ事象の不正確な測定値を生じさせる可能性がある。従って、これらのアプローチは、典型的には、マルチパス信号を、正確な位置結果を提供するための削除されるべきエラーとみなす。
【０００５】
位置決め（測位）の問題に対する、代表的な第１のクラスのアプローチは、慣性システム、すなわち、各々の停止において更新するゼロ（零）速度を有する慣性システムと、ドップラー測定値か、気圧高度か、磁力計か、又は視覚システムによって、増強又は更新することとを含む。このクラスのアプローチは、繰り返される交差点を識別するためにか、又は動的な慣性の較正を提供するために、時として特別なアルゴリズムを必要とし、典型的には、建築物の寸法及び配置の事前知識を必要とする。これらアプローチのうちの幾つかは、影となる（受信不良となる）真っただ中の散乱された環境内において、信号が忠実に受信されることが可能となる程度にまで、全地球測位システム（ＧＰＳ）か又は、Ｅ９１１のような携帯電話（セルフォン）のナビゲーションの支援を利用する。
【０００６】
アプローチの第２クラスは、802.11xか、ブルートゥースか、又は別のＲＦ技法を用いて、協同のインフラ（構造基盤）を使用する。いくつかのインフラは、そのような目的のために特別に設計されている。１つの特定のアプローチは、幾つかのトランシーバによって受信されるような、ネットワーク素子送信器の信号強度を用い、該信号強度は、更にマッピングされて、吸収、屈折、及び反射のような影響の原因となる既知の建築物構造に対して較正される。別のアプローチは、特化された送信器か又は受信器のセットを三角測量することによってもたらされるキャリア位相測定値を使用する。この第２のクラスにおいて、ＲＦＩＤ技術が、比較的狭いエリアにわたって対象物を位置付けるために使用されることもまた可能であるが、より広いエリアにわたって実用化されるためには、典型的には、多くの統合センサを必要とする。
【０００７】
第１のクラスの独立した（自律性の）アプローチに関連した１つの問題はコストである。これらのアプローチは、良好な性能を提供することができるが、それらは典型的には、コストがかかる素子を用い、該素子がコストがかかるシステム内へと統合化される。更には、このクラスのアプローチは、典型的には、各々の人間か又は対象物に、ナビゲーション装置を携帯することを要求し、このことが、システム全体を高価にさせる可能性がある。そのようなアプローチは、独立性（自律性）が重要で且つコストがあまり重要でない場合の用途には適合するが、商業的な用途に対するそれらのアプリケーションは、コストが法外に高くなる可能性がある。
【０００８】
第２のクラスによるアプローチは、それらのインフラ（構造基盤）内において組み込まれたそれらの主要コストが有する利点を享受する。そのコストは、引き続いて、802.11x技術が備わった携帯情報端末のような比較的低コストの受信器か又は中継器を各々が携帯することになる個別ユーザによって償却される。しかしながら、第２のクラスによるアプローチは、依然として、典型的には、マルチパス信号に対する感度による影響を受ける。
【特許文献１】米国仮特許出願第６０／６４５，３９０号明細書
【特許文献２】米国特許第６，６９３，５９２号明細書
【非特許文献１】A.H.Jazwinski著,「Stochastic Processes and Filtering Theory」, Academic Press, New York, 1970年
【非特許文献２】A.Gelb編集,「Applied Optimal Estimation」,MIT Press, 1974年
【非特許文献３】B.Ristic, S.Arulampalam, N.Gordon著,「Beyond the Kalman Filter」,Artech House, Boston, 2004年
【発明の開示】
【発明が解決しようとする課題】
【０００９】
従って、アーバンキャニオンか、又は建築物などの囲まれた構造のような、散乱した環境内において移動している人々か又は存在（アセット）を、これら環境内において生成されたマルチパス信号からの悪影響を受けること無く、測位すること（位置決めを行うこと）が可能な改善されたシステム及び方法の必要性がある。
【課題を解決するための手段】
【００１０】
発明の概要
本発明は、アーバンキャニオンか、又は建築物などの囲まれた構造のような、散乱した環境内において移動する人々か又は存在（アセット）の位置決め（測位）を行うためのシステム及び方法に関する。本発明のシステム及び方法は、ダイレクトパス（ＤＰ）が無いか又は検出されない時でさえ、実質的なエネルギーが利用可能であることを認識し、及び、ローカルな環境内における様々な物体によって送信された信号が反射された時には、結果として生じるマルチパス信号によって十分な情報が伝達されることを認識する。
【００１１】
本発明の実施形態は、現実の物理的な環境内において、小さな空間の隔離、多くのマルチパス構成要素は、共有され共通であると仮定することを認める。従って、本発明の実施形態は、再帰的推定技法を用いて環境のモデルを構築する。これらの技法は、新たなパスの生成を、それらが重要になったものとして許可し、及び、古いパスの終結を、それらが無関係になったものとして許可する。従って、各々の個別のマルチパス要素と、各々の受信した送信位置との間の信号を含む利用可能な情報の効率的な使用によって、ローカルなインフラ（構造基盤）の事前知識の必要性が排除される。
【００１２】
本発明の実施形態は、位置付け及び追跡する場合に、測定されたマルチパスＴＯＡ遅延を直接的に使用する。レイクフィルタ（Rake filter）か又は類似の手段を用いて、差分の遅延がリアルタイムに追跡される。後述されるように、ある単一時間におけるシングルパスについての差分の遅延と、タグの３次元の位置との間の関係を、任意の数の反射の場合について４つのパラメータを用いて説明することができる。この関係の理解によって、マルチパスパラメータを、ダイレクトパス（直接的な経路）、インダイレクトパス（間接的な経路）、及び他の利用可能な測定値を用いて、リアルタイムに推定することができる。ダイレクトパスが遮断されるか又は他の測定値が利用不可能になった場合には、従って、これらのマルチパスパラメータを、３次元の位置付けか又は追跡を行う場合の推定技法と共に、引き続き用いることができる。
【００１３】
差分の遅延とタグ位置との間の関係を説明するマルチパスパラメータの知識（情報）が、完全マルチパスの解を決定する（すなわち、正反射点の位置を提供する）には不十分である時でさえ、ジオロケーション（地理位置情報）の解は、典型的には、完全マルチパスの解を必要としない。マルチパス信号を用いるジオロケーションは、ダイレクト及びインダイレクトな経路の測定値の組み合わされた観測からか、又は位置の、事前の初期知識から、マルチパスパラメータが決定可能である時の状況において、成し遂げられることが可能である。マルチパスパラメータの積算された時間の増大が、従って、建築物構造の内部の画像を提供するために使用することができる可能性がある。
【００１４】
ある特定の状況において、線形化（例えば拡張カルマンフィルタ）を用いるフィルタリング方法は適切ではない。何故ならば、測定値とフィルタ状態との間の関係が、非常に非線形である場合があるからである。例えば、マルチパス反射は、マルチパスパラメータ及びタグ位置を含む、少なくとも双線形（バイリニア）測定値を結果として生じる可能性がある。少なくともこの理由のため、線形化の想定か又は特定の確率分布に依存しない、パーティクルフィルタリング（particle filtering）のような、非線形のフィルタリング技法を使用することができる。初期のタグ位置の不確定性が非常に大きい場合には、パーティクルフィルタリングは、特に有用である可能性がある。初期のタグ位置が、比較的小さい場合には、拡張カルマンフィルタで十分である可能性がある。
【００１５】
一態様において、本発明は、マルチパス無線信号を用いて位置決めを行うための方法に関する。初期位置が決定されて、少なくとも１つの位置において受信された少なくとも１つのマルチパス無線信号における少なくとも１つの特性に対して、該初期位置のような、該少なくとも１つの位置に関する情報が生成される。後続する無線マルチパス信号が受信されて、該後続する信号と該生成された情報とに基づいて、ある位置が確立される。前記無線マルチパス信号を、ＧＰＳ送信器のような送信器によって生成することができる。前記生成された情報は、典型的には、少なくとも１つのパラメータを含む。
【００１６】
一実施形態において、初期位置を決定することは、タグに、複数位置において、少なくとも１つの送信された信号を受信させることを含む。タグによって受信された少なくとも１つの信号の大きさは、各位置において、及び初期位置においても決定される。次いで、決定された信号の大きさが、初期位置を推定するために統計モデルによって使用される。該統計モデルには、タグと、複数の信号送信器の各々との間の信号間隔に関連付けられた信号伝搬特性を示すデータを含めることができる。別の実施形態において、初期位置を決定することは、受信信号強度の測定値のような、初期位置のダイレクトパス測定値を生成することを含む。更に別の実施形態において、初期位置を決定することは、初期位置の事前知識を利用することを含む。
【００１７】
別の実施形態において、位置を確立することは、複数の候補位置を計算することと、そこから最もふさわしい位置を選択することとを含む。前記最もふさわしい位置を、統計的な試験を用いて決定することができる。
【００１８】
更に別の実施形態において、情報を生成することは、内部環境の少なくとも部分的なジオメトリと、その中のマルチパス信号のジオメトリとを特徴付けることを含む。該マルチパスのジオメトリは、前記環境内における反射点に対する無線源からのインダイレクトパスセグメント長を含む。インダイレクトパスセグメント長を、パス遅延測定値によって表すことができる。
【００１９】
様々な実施形態において、少なくとも１つの特性には、パス遅延か又はパスセグメント長を含めることができる。これらの実施形態において、生成された情報は、パス遅延か又はパスセグメント長を、環境内における位置に関連付けることができる。
【００２０】
別の態様において、本発明は、マルチパス無線信号を用いて位置決めを行うための装置に関する。該装置は、少なくとも１つの位置において受信された少なくとも１つのマルチパス無線信号における少なくとも１つの特性に対する、初期位置のような該少なくとも１つの位置に関する情報を有するメモリと、後続する無線マルチパス信号を受信するための受信器と、前記後続する信号と前記生成された信号とに基づいて位置を決定するための回路構成とを含む。前記少なくとも１つの特性は、例えば、パス遅延か又はパスセグメント長とすることができ、前記情報は、従って、パス遅延か又はパスセグメント長を、環境内における位置に関連付けることができる。無線マルチパス信号を、送信器、例えば全地球測位システム（ＧＰＳ）送信器によって、生成することができる。
【００２１】
一実施形態において、情報は、内部環境の少なくとも部分的なジオメトリと、そこの中のマルチパス信号のジオメトリとを特徴付ける。該マルチパスのジオメトリは、該環境内における反射点に対する無線源からのインダイレクトパスセグメント長を含む。インダイレクトパスセグメント長を、パス遅延測定値によって表すことができる。
【００２２】
上述の及び他の、本発明の特徴と利点とが、説明、図面、及び添付の特許請求の範囲からより明らかにされることになるであろう。
【実施例】
【００２３】
本発明の利点は、添付図面を、それに伴う説明と連携して参照することによって、より良く理解されることが可能である。
【００２４】
図面内において、同様の参照文字は、概して、異なる図を通して、対応する部分を指す。該図面は、一定の縮尺に必ずしも従っているとは限らず、本発明の原理及び概念において代りに強調されている。
【００２５】
本発明の詳細な説明
ＧＰＳが使用できない、建築物内側か又は散乱された都市環境内における人々及び物体を追跡するために多くのアプローチが提案されてきた。これらアプローチの大部分は、ＲＦ現象を利用し、性能は典型的には、人間か又は対象となる物体に対する多くの基準源からのダイレクトパス距離を決定する能力に依存する。これらアプローチの多くは、いくつかは成功するが、大抵のものは、性能と信頼性とがＲＦマルチパスによって制限される。典型的な建築物の内側で生じる多くの反射とパス減衰のため、ダイレクトパス距離を一貫して確実に識別することに対する無能性が、これらアプローチの有効性（効率）を制限してきた。
【００２６】
概要を簡単に述べると、本発明の実施形態は、ある信号処理技法を提供する。該信号処理技法を、屋内のＲＦの能動的か、受動的か、又は支援された位置確認アプローチに適用することができる。ノイズとしてマルチパス信号を扱って、マルチパス挿入エラーを軽減するように試みるのではなく、この技法は、フィルタか又は推定器、例えば非線形フィルタと共に、マルチパス信号を追加的な測定値として利用する。該フィルタは、インダイレクト及びダイレクトパス測定値（又は任意の他の利用可能な信号）を使用して、観測可能なインダイレクトパスのパラメトリックモデルを構築する。１つか又は複数のダイレクトパス測定値が、（例えば、障害のせいで）その後、失われる場合には、インダイレクトパス測定値を用いて、人間か又は対象となる物体の位置の推定をフィルタが維持する。従って、本発明の実施形態は、ダイレクトパス技法が失敗する場合のシナリオにおいて、追跡精度を維持するために、屋内のナビゲーション用途においてマルチパス信号を用いることができ、ローカルなインフラ（構造基盤）の事前知識の必要性を低減するか又は除去し、全ての利用可能な情報の効果的な利用を可能にし、情報の冗長性を有利に利用する。
【００２７】
これらの展開（発展）に関連して、平面からの正反射に対するインダイレクトパス距離を、２次元では２つのパラメータだけを用いて、及び、３次元では３つのパラメータだけを用いて、正反射の数に関係無く、モデル化することができることが決定される。このモデルを使用して、下記に詳細に説明されるように、上記のパラメトリックモデルが、リアルタイムに推定される。
【００２８】
一実施形態において、ダイレクト及びインダイレクトマルチパス信号の時間遅延測定値が、１つか又は複数のデータ関連付け、初期化、パラメータ推定、及び追跡フィルタを用いて処理される。該データ関連付けのアルゴリズムは、断続的なダイレクト及びインダイレクトパス測定値の追跡機能を高めるために使用される。別の実施形態において、パーティクルフィルタを用いて、インダイレクトパスパラメータ推定と追跡との両方に拡張カルマンフィルタを使用することが可能な領域に対する初期位置の不確定性を低減することができる。更なる実施形態において、推測航法センサを使用して、測定値におけるエラーの増殖を緩慢にすることができるか、又はデータ関連付け処理を支援することができる。
【００２９】
ダイレクトパス測位の概要
囲まれた構造の内側か、又は接近した距離の建物間を、個々の位置か又は存在（アセット）が移動する時に、都会環境内における位置決めによる重要な結果（issue）が、該個々の位置か又は存在を決定している。この問題に対する、ＲＦ現象を利用する多くのアプローチが提案されてきた。多くの様々な屋内ナビゲーションアプローチが、現在開発されている最中であり、該アプローチは、人間か又は対象となる物体（その一例が、図１内に示されている）に取り付けられたモバイルタグか又は変換器と連携して、既知の固定された送信器か、中継器か、又はビーコンのセットを用いる。
【００３０】
ＲＦ位置センサが、実現可能な位置測定基準として、ＲＳＳ、ＡＯＡ、ＴＯＡ、ＴＤＯＡ、及び遅延電力プロファイルのシグニチャを測定することを可能にする。一般的には、ＲＳＳ及びＡＯＡを測定する容易性は、信頼できない結果を生む。ＴＯＡセンサは、距離の、より正確な測定を提供するが、ＴＯＡの性能は、典型的には、人間か又は対象となる物体に対する多くの基準源からのダイレクトパス距離を決定することに依存し、典型的には、第１の検出ピークの到着時間を使用する。ダイレクトパス距離を決定するための過去のアプローチは、超広帯域のプローブ波形を用いてノイズフロア（noise floor）に到来する視界のダイレクトラインを追跡することを含んでおり、それら両方は、ほんのわずかな利点しかない。
【００３１】
図２Ａは、内部トランシーバと、タグを装備した動き回る物体とを有する建築物の一例を示す。図２Ａはまた、該トランシーバからの、タグに到来する最後の１０の射線（ray）を示し、そのうちの幾つかの射線は直接的に到来し、幾つかの射線は金属製の物体によって遮断され、幾つかの射線は屈折させられて、複数の反射の後に到来する。
【００３２】
図２Ｂは、任意の単一トランシーバからの、動き回るタグにおいて受信された信号の一例を示す。該トランシーバからタグへのＤＰは、ＴＯＡの期待値をもたらす。多くの反射及び伝送に遭った他のインダイレクト信号は、ＤＰ信号の後に、タグに到達することになる。実際には、トランシーバの帯域幅が制限されるので、タグにおいて受信された信号は、インパルス信号のような振幅と到達時間とを有することとなるが、形状はパルス信号により類似することとなる。これらのパルス形状の合計が、タグにおいて受信された信号を形成し、該信号は、本明細書内において「チャンネルプロファイル」と呼ばれる。音楽アルゴリズムか、又はＤＳＳＳ信号による従来の相互相関技法のような超解像技法を用いて、これらＴＯＡ要素を高解像度に分解することができる。
【００３３】
図２Ｂを更に参照すると、屋内ナビゲーションシステムは、典型的には、ＤＰの推定されたＴＯＡとして、検出閾値を越えた、前記チャンネルプロファイル内の第１のピークを扱う。タグが単一パス環境内において動作していると仮定すると、実際の期待されたＤＰと、推定されたＤＰとが典型的には等しくなることとなる。しかしながら、マルチパス信号が生成されている環境内において、チャンネルプロファイルのピークは、期待されたＴＯＡからはずれるようにシフトされて、結果として推定エラーが生じる。この誤ったＴＯＡ推定によって生じられる分布エラーは、本明細書内において「距離測定値エラー」（ＤＭＥ：distance measurement error）と呼ばれる。システムの伝送帯域幅が拡大すると、異なるパスから到来するパルスが、より狭まることとなり、第１のＤＰの推定されたＴＯＡが、期待されたＴＯＡ値に近づく。
【００３４】
遮断された視界線（line-of-sight）マルチパス条件において、ＤＰ信号の大きさが検出閾値を下まわった時には、チャンネルプロファイル内に現われる第１のＤＰは、ＤＰの実際の到達時間に依存しておらず、未検出ダイレクトパス（ＵＤＰ：undetected direct path）条件が発生することにより、大きなＤＭＥ値が生じることになる。図２Ｂは、そのようなＵＤＰ条件に起因して、いかに大きなＤＭＥが発生する可能性があるかを示し、このケースの場合には、２００ＭＨｚの帯域幅で送信されたパルスについて、シミュレートされた射線追跡（レイトレーシング）の結果生じている。受信した最強度のパスと、ＤＰそれ自体との間の差分が、受信器の動的帯域幅を越えるので、その結果はＵＤＰ条件である。ＵＤＰ条件の結果生じるＤＭＥは、チャンネルインパルス応答の全てのマルチパス成分を解決するために、たとえ帯域幅が超広域に拡大されたとしても、残存する可能性がある。
【００３５】
マルチパス測位
図３に関して、本発明の実施形態は、タグの初期位置の推定を開始する（ステップ３００）。このような推定を、例えば、ＧＰＳ測位か、又はタグ位置の事前知識を利用することによってか、又は既知パスにタグの軌道を強いることによって、決定することができる。
【００３６】
初期位置の推定によって、受信したマルチパス信号の決定された特性に該位置を関連付けて（ステップ３０４）、パラメトリックモデルが展開される。例えば、タグにおいてか、又はトランシーバのうちの１つにおいてかのいずれかにおいて、マルチパス信号を受信することができる。該パラメトリックモデルと、その展開が、より詳細に以下に説明される。
【００３７】
受信したマルチパス信号の決定された特性にタグの位置を関連付けるこの情報によって、後続するマルチパス信号が、受信される（ステップ３０８）。本発明によるパラメトリックモデルを用いるフィルタに対して、後続するマルチパス信号を適用することは、タグの位置を推定する結果となる（ステップ３１２）。
【００３８】
図４に関して、一実施形態において、マルチパス無線信号の特性に位置を関連付ける情報を生成するためにか、又は後続するマルチパス信号を用いて位置を確立するために、受信したマルチパス信号が使用される前に、該受信したマルチパス信号は、最初に前処理される。前処理ステップの出力は、後続するステップ内においてインダイレクトパス（マルチパス）パラメータの正確な推定を可能にするための十分な長さの時間の間存続するインダイレクトパスのセットを含む。
【００３９】
前処理を行うことには、例えば、レイクフィルタを使用して（ステップ４００）、個々の信号パス（ダイレクトとインダイレクトとの両方）を検出することと、測定値遅延のセット（等価的には、パスに沿って測定された距離）を計算することとを含めることができる。レイクフィルタ前処理を、次いで、データ関連付けアルゴリズムを適用すること（ステップ４０４）によって続けることができる。該データ関連付けアルゴリズムは、パスのクロスオーバーの影響を軽減するように動作して、ドロップアウトした期間を通じて該パスを追跡する。該データ関連付けアルゴリズムからの出力は、後続するステップ内においてインダイレクトパス（マルチパス）パラメータの正確な推定を可能にする十分な長さの時間の間存続する、インダイレクトパスのセットである。
【００４０】
図５に関して、未処理の形態におけるか、又はデータ関連付けアルゴリズムのような前処理ステップからかのいずれかにおけるインダイレクトパスデータが、ダイレクト及びインダイレクトパスを検出し且つ識別するために処理される（ステップ５００）。後述されるように、マルチパスモデルを使用して、インダイレクトパス（マルチパス）パラメータの初期推定値を決定するために（ステップ５０４）、インダイレクトパス測定値が使用される。次いで、非線形フィルタリング技法が、前記初期推定値と、後続するダイレクトパス及びインダイレクトパス測定値と共に使用されて、タグの位置が決定される（ステップ５０８）。
【００４１】
図６は、データ関連付け処理の別の実施形態を示す。この実施形態において、最小二乗法（ＯＬＳ：ordinary least-squares）分析を用いて、低オーダの多項式モデルによって経路長（パスの長さ）がモデル化される。測定値の差分が閾値よりも低く、且つ、全ての既存パスにわたって最小の場合には、新規の測定値が既存のパスに関連付けられる。そうでなければ、新規のパスが作成される。
【００４２】
マルチパスモード及びパラメータ推定
一般に、マルチパスの解を決定する処理は、ナビケーションか又はジオロケーションの解を取得するための処理と区別することができる。マルチパスの解は、マルチパスのジオメトリから生じ、全ての反射面の位置及び配向と同様に、無線源及び受信器の位置が十分な程度の特定性において既知であるコンフィギュレーションについて決定されることが可能である。
【００４３】
マルチパスの解は、一般的には、無線源から反射点の位置までのインダイレクトパスセグメントの長さを示す情報を含む。インダイレクトパスセグメント長を、マルチパスを含む個々のパスセグメントに沿った無線信号伝搬に関連した伝搬遅延に対応するパス遅延測定値によって表すこともできる。それとは対照的に、ジオロケーションの解は、一連のインダイレクト及びダイレクトパス遅延測定値からタグ位置を推定する。結果として、ジオロケーションの解は、全てのマルチパス変数の推定を必要としない代りに、それらのサブセットだけで、機能を果たすことができる。
【００４４】
マルチパス条件下において、タグ位置をインダイレクトパス長に関連付けるモデルを用いて（ステップ３１２）、インダイレクトパス遅延の測定値から（ステップ３０８）、タグ位置を直接的に追跡することができる。典型的には他の手段によって取得される初期タグ位置推定値（ステップ３００）を用いて、このモデルをリアルタイムに生成し及び更新することができる（ステップ３０４）。
【００４５】
図７に関して、このモデルの導出は、単一平面からの理想的な正反射を仮定することによって開始される。次式は、インダイレクトパス長τ＝ｄ０＋ｄ１と、タグ位置ベクトルｒとの間の関係を説明する。
【００４６】
【数１】
2008527394.tif 000003
ここで、ｗ１は、単位ベクトルであり、ｓは、信号源の位置であり、ｕ１は、反射面の配向を画定する単位ベクトルであり、ｃ１は、ｕ１Ｔｒｐ＝ｃ１を満足するスカラーであり、これによるｒｐは、反射面内における任意の点であり、Ｅ１＝２ｕ１ｕ１Ｔ−Ｉは、該反射面に関連した変換（すなわち反射）行列を表す。
【００４７】
反射面と、ｓとが固定されていると仮定すると、時間ｔにおけるインダイレクトパス長は、次式を満足する。
【００４８】
【数２】
2008527394.tif 000004
ここで、ρ０は一定である。従って、
【数３】
2008527394.tif 000005
となる。
【００４９】
信号源が反射面にあまり接近していない場合には、適切に小領域にｒ（ｔ）が制限されると仮定されるのであれば、４つのマルチパスパラメータｗ１（ｔ）及びτ０（ｔ）が、本質的には一定である。τ（ｔ）とｒ（ｔ）のほぼ同時の測定値が利用可能であり且つタグが動作中の場合には、このことは、標準的なフィルタリングか又は推定技法を用いて、領域内におけるｗ１とτ０の推定を可能にする。
【００５０】
ｒ（ｔ）の他の測定値が利用不可能である時には（例えば、ダイレクトパスが失われた時には）、インダイレクトパス測定値を用いてｒ（ｔ）を追跡するためにこれらのマルチパスパラメータの推定値を引き続き使用することができる。ｒ（ｔ）、ｗ１（ｉ）、及びτ０（ｉ）が、ランダムな変数であり、τ（ｉ）の測定値がエラーを受けやすいので、タグ位置ｒ（ｉ）を正確に追跡するためには非線形フィルタが必要とされる。
【００５１】
この分析は、単一マルチパスのはね返りに対して示されるが、その結果は、４つのマルチパスパラメータが異なる解析的形態を有するにもかかわらず、発生するマルチパス反射の数とは無関係に適用可能である。例えば、図８は、２つの正反射の場合の、一般化された３次元マルチパスジオメトリを示す。該ジオメトリは、図７内に図示されたものと類似しており、単位ベクトルｕ２及びスカラー値ｃ２により画定される第２の反射面を含んでいる。次式は、インダイレクトパス長τ＝ｄ０＋ｄ１＋ｄ２と、タグ位置ベクトルｒとの間の関係を説明する。
【００５２】
【数４】
2008527394.tif 000006
ここで、
【数５】
2008527394.tif 000007
であり、Ｅ２は、第２の反射面ｃ２に関連した変換（すなわち反射）行列を表す。
【００５３】
反射面と、信号源とが両方とも固定されていると仮定すると、時間ｔにおけるインダイレクトパス長を、次式のように記載することができる。
【００５４】
【数６】
2008527394.tif 000008
ここで、ρ０２は一定である。従って、
【数７】
2008527394.tif 000009
上式は、１つの反射の場合（すなわち、（１））と同じ形態である。
【００５５】
タグ位置ｒと４つのマルチパスパラメータ｛ｗ２，τ０２｝を同時に推定するために、式（１）と（２）とが両方とも、適切なフィルタ（例えば、拡張カルマンフィルタ、非線形フィルタ、パーティクルフィルタ）によって使用されることが可能な測定値式の形態にある。
【００５６】
付録Ａにおいて説明されるように、単一の及び二重のマルチパスの場合を、任意の数の正反射の場合のタグ位置及び４つのパラメータの関数としてインダイレクトパス長についての式を公式化するために一般化することができる。ｎ個の反射の場合には、４つのパラメータは、｛ｗｎ，τ０ｎ｝である。ここで、ｗｎは、最後の反射面からのパスの方向における単位ベクトルであり、τ０ｎは、スカラーのジオメトリ依存性オフセットである。ｗｎを２つの角度の関数として表すことができるため、推定されることになる３つの独立したパラメータだけが存在する。
【００５７】
最も堅牢なモデル（すなわち、４つのパラメータか又は３つのパラメータ）がアプリケーション依存性であることが、シミュレーションから経験的に決定される。
【００５８】
統計的なモデリングと非線形フィルタリング
一実施形態において、インダイレクトパス（すなわち、式（１）及び式（２））についての時間遅延τ（ｔ）が、マルチパス単位ベクトルｗｎ（ｔ）及びタグ位置ｒ（ｔ）の双線形関数であるため、上述のマルチパスパラメータ及びタグ位置を同時に推定するために、非線形フィルタリングが用いられる。
【００５９】
ある特定のアプリケーションに最も適する、特定のタイプの非線形フィルタは、該アプリケーションの詳細に依存する。初期タグ位置におけるエラーが、比較的小さい場合には、拡張カルマンフィルタで十分である可能性がある（A.H.Jazwinskiによる「Stochastic Processes and Filtering Theory」, Academic Press, New York, 1970年、A.Gelb編集の「Applied Optimal Estimation」,MIT Press, 1974年を参照のこと）。初期タグ位置におけるエラーが、比較的大きい場合には、パーティクルフィルタが、良好な性能を提供する可能性がある（B.Ristic, S.Arulampalam, N.Gordonによる「Beyond the Kalman Filter」,Artech House, Boston, 2004年を参照のこと）。以下の説明を簡単化するために、拡張カルマンフィルタの使用が仮定される。
【００６０】
統計的モデルは、状態ベクトル、統計的伝搬モデル、及び統計的測定値モデルからなる。状態ベクトルがｘ（ｔ）の場合には、統計的モデルは、例えば、次式のように表すことができる。
【００６１】
【数８】
2008527394.tif 000010
ここで、ｖ（ｔ）は、タグ速度である。
【００６２】
（３）の状態ベクトル定義を想定すると、時間ｔｉから時間ｔｉ＋１まで状態ベクトルを伝搬させるために実現可能な統計的モデルは、次式のようになる。
【００６３】
【数９】
2008527394.tif 000011
ここで、Φ（ｔｉ＋１，ｔｉ）は、状態遷移行列であり、ｕ（ｔｉ）は、既知の制御入力であり、Ｉは、単位行列であり、ｑ（ｔｉ）は、真の伝搬モデルの情報（知識）の不確定性をモデル化するために使用されるゼロ平均白色駆動ノイズである。このモデルは、一般に、狭い時間ステップ（ｔｉ＋１−ｔｉ）と、比較的一定なタグ速度とについて適用可能である。
【００６４】
フィルタ内において、推定される状態が次式を用いて伝搬される。
【００６５】
【数１０】
2008527394.tif 000012
フィルタはまた、エラー共分散行列Ｐ（ｔｉ）＝Ｅ｛ｅ（ｔｉ）ｅ（ｔｉ）Ｔ｝の推定値を搬送する。ここで、
【数１１】
2008527394.tif 000013
であり、上付文字Ｔは、転置行列であることを示す。エラー共分散行列が、次式を用いて、フィルタ内において伝搬される。
【００６６】
【数１２】
2008527394.tif 000014
ここで、
【数１３】
2008527394.tif 000015
であり、ここでの、Ｅ｛・｝は、括弧内における量の数学的予測であることを示す。
【００６７】
想定したモデル下において、２つの別個の測定値が利用可能である。すなわち、ダイレクトパス遅延測定値と、インダイレクトパス遅延測定値である。時間ｔｉにおける送信ノードからタグへの、ダイレクトパス遅延測定値は、次式のようにモデル化される。
【００６８】
【数１４】
2008527394.tif 000016
ここで、ｓは、送信ノードの既知の位置であり、‖ｚ‖は、ベクトルｚの長さであることを示し、測定値ノイズｎｄ（ｔｉ）は、分散σｄ２（ｔｉ）＝Ｅ｛ｎｄ２（ｔｉ）｝を有するゼロ平均白色ノイズとしてモデル化される。送信ノードからタグへのインダイレクトパス遅延測定値は、次式のようにモデル化される。
【００６９】
【数１５】
2008527394.tif 000017
ここで、ｎ（ｔｉ）は、分散σｎ２（ｔｉ）を有するゼロ平均白色ノイズである。上述のように、インダイレクトパス測定値に対するこのモデルは、任意の数の正反射に対して適用可能である。
【００７０】
測定値の更新
マルチパス環境内において成功しているナビゲーションは、マルチパスパラメータと状態ベクトルとを、リアルタイムで同時に正確に推定することを必要とする。この定式化によって、これらの目的（ゴール）が、非線形フィルタイリング理論の周知のツールを用いることによって実現されることになることが可能となる。具体的には、インダイレクトパス測定値の双線形形態が、それ自身を、特に簡単な再帰的非線形フィルタ機構化に加える。
【００７１】
ダイレクトパス測定値を用いて、状態ベクトル推定値と、エラー共分散行列推定値とを更新する処理を検討する。これは、標準的な非線形フィルタリング方法を直接使用することであり、拡張カルマンフィルタリングと、パーティクルフィルタリングとを含む。簡単化のために、以下の説明は、拡張カルマンフィルタの使用を仮定する。
【００７２】
ダイレクトパス測定値は、（７）によってもたらされる。従って、この測定値のフィルタの推定値は、次式のようになる。
【００７３】
【数１６】
2008527394.tif 000018
ここで、更新する測定値の直前の変数の値に関連していることを示すために全体を通して、プライム符号（）’が使用されている。次式を用いて、状態ベクトル推定値とエラー共分散行列推定値とを更新するために、標準的な拡張カルマンフィルタが次いで使用される。
【００７４】
【数１７】
2008527394.tif 000019
ここで、υｄ（ｔｉ）＝ｄｍｅａｓ（ｔｉ）−ｄ’ｍｅａｓ（ｔｉ）は、測定値の剰余であり、測定値行列Ｈ’ｄ（ｔｉ）は次式のように表される。
【００７５】
【数１８】
2008527394.tif 000020
そして、Ｋｄ（ｔｉ）は、
【数１９】
2008527394.tif 000021
によってもたらされる拡張カルマンフィルタ利得である。
【００７６】
インダイレクトパス測定値を用いて、状態ベクトル推定値と、エラー共分散行列推定値とを更新することを、次に検討する。インダイレクトパス測定値は、（８）によってもたらされる。これにより、測定値のフィルタの推定値は、次式のように表される。
【００７７】
【数２０】
2008527394.tif 000022
式（８）と式（１４）とを用いて、測定値を次式のように書くことができる。
【００７８】
【数２１】
2008527394.tif 000023
【数２２】
2008527394.tif 000024
ここで、ｅ’ｒ（ｔｉ）＝ｒ’（ｔｉ）−ｒ（ｔｉ）は、タグ位置ｒ（ｔ）の事前推定値におけるエラーであり、ｅ’ｗ（ｔｉ）＝ｗ’（ｔｉ）−ｗ（ｔｉ）は、単位ベクトルｗ（ｔ）の事前推定値におけるエラーであり、ｅ’τ０（ｔｉ）＝τ’０（ｔｉ）−τ０（ｔｉ）は、オフセットパラメータτ０（ｔｉ）の事前推定値におけるエラーである。
【００７９】
エラーが十分に小さく、式（１５）内における第２のオーダのエラー項を無視することができると仮定すると、その結果、インダイレクトパス遅延測定値について、
【数２３】
2008527394.tif 000025
となる。
【００８０】
続いて、測定値の剰余は、
【数２４】
2008527394.tif 000026
によってもたらされる。
【００８１】
幾つかの異なるやり方において、インダイレクトパス測定値を用いて、状態ベクトルを更新することができる。以下の例は、状態ベクトル推定値と、マルチパスパラメータ推定値とが、２つの分離ステップ内において連続的に更新されると仮定する。両ステップ内において、拡張カルマンフィルタが使用されると仮定される。
【００８２】
単一送信ノードからの単一マルチパス遅延測定値を処理する場合の計算が示されている。同時間における複数のノードからの複数のインダイレクトパス遅延測定値を、同一の拡張カルマンフィルタ式の形態を用いて連続的に処理することができる。
【００８３】
インダイレクトパス測定値を用いて、状態ベクトル推定値を更新するために、測定値の剰余υｘ（ｔｉ）が、次式の形態で書かれる。
【００８４】
【数２５】
2008527394.tif 000027
ここで、
【数２６】
2008527394.tif 000028
であり、測定値ノイズはすなわち、
【数２７】
2008527394.tif 000029
は、次式に示す分散を有するゼロ平均白色ノイズとしてフィルタ内においてモデル化される。
【００８５】
【数２８】
2008527394.tif 000030
ここで、
【数２９】
2008527394.tif 000031
は、推定値ｗ’（ｔｉ）についてのエラー共分散行列であり、σ’τ０２（ｔｉ）は、推定値τ’０（ｔｉ）についてのエラー分散である。この想定により、測定値の剰余の分散は、
【数３０】
2008527394.tif 000032
となる。
【００８６】
次式のように、標準的な拡張カルマンフィルタの式を用いて、状態ベクトルとエラー共分散行列との更新された推定値が計算される。
【００８７】
【数３１】
2008527394.tif 000033
ここで、Ｋｘ（ｔｉ）は拡張カルマンフィルタ利得であり、すなわち、
【数３２】
2008527394.tif 000034
である。
【００８８】
マルチパスパラメータ推定値を更新するために、パラメータベクトルが、
【数３３】
2008527394.tif 000035
として画定される。
【００８９】
式（１７）においてもたらされる、測定値の剰余は、次いで、次式のように書かれる。
【００９０】
【数３４】
2008527394.tif 000036
ここで、
【数３５】
2008527394.tif 000037
であり、測定値ノイズは、
【数３６】
2008527394.tif 000038
であり、次式に示す分散を有するゼロ平均白色ノイズとしてフィルタ内においてモデル化される。
【００９１】
【数３７】
2008527394.tif 000039
ここで、
【数３８】
2008527394.tif 000040
は、タグ位置推定値ｒ’（ｔｉ）についてのエラー共分散行列である。
【００９２】
これらの想定下において、測定値の剰余の分散は、次式のように表される。
【００９３】
【数３９】
2008527394.tif 000041
ここで、
【数４０】
2008527394.tif 000042
は、事前パラメータベクトル推定値ｅ’α（ｔｉ）＝α’（ｔｉ）−α（ｔｉ）の共分散行列である。次式のように、標準的な拡張カルマンフィルタの式を用いて、マルチパスパラメータベクトルとエラー共分散行列との更新された推定値が計算される。
【００９４】
【数４１】
2008527394.tif 000043
ここで、Ｋｘ（ｔｉ）は拡張カルマンフィルタ利得であり、すなわち、
【数４２】
2008527394.tif 000044
である。
【００９５】
例示的なマルチパスのアルゴリズム例
マルチパス信号遅延を利用する屋内ナビゲーションの比較的簡単な一例が、例示的な目的のために、本明細書において示されている。示されたシナリオは、２つの固定されたトランシーバノードと動き回るトランシーバタグとを利用する平面ジオメトリ内における屋内のジオロケーションの簡単なケースの代表である。そのようなコンセプトは、多くのＵＷＢか又は802.11のようなアプローチの代表であり、該アプローチは、軍事用途と商業用途との両方の場合について提案されている。
【００９６】
図９Ａ〜９Ｄの例示的なジオメトリは、初期には障害物が無く動き回るタグを図示しており、それにより、直接的視界線（ＤＬＯＳ：direct line-of-sight）パスが、観測できることとなり、ジオロケーションの目的のために、第１のリターンを用いることができることとなる。該タグは、ある一定の高さにおける２次元（北方及び東方）内において移動していると仮定される。ジオメトリはまた、いずれかの送信源に関連してマルチパス信号を導入することを可能にする正反射面も含む。
【００９７】
図９Ａは、初期時間（ｔ＝０）におけるジオメトリを示す。その最初の５０秒の移動において、両ダイレクトＲＦパスに対してタグが露出されており、両ノードに関するこのダイレクトパスデータを利用して、その自身の位置を決定し、同時に、両インダイレクトパスについてのマルチパスパラメータを推定する。図９Ｂは、ｔ＝５０秒におけるジオメトリを示し、該ｔ＝５０秒は、全てのパスが遮断されないシミュレーションにおける最終時点である。
【００９８】
図９Ｃは、ノード＃１に対するダイレクトパスが遮断されることになる時であるｔ＝５１秒におけるジオメトリを示す。この時点では、システムは、本発明に従った位置推定のために、ノード＃１のインダイレクトパスと、ノード＃２のダイレクトパスとを使用する。ノード＃２からのインダイレクトパス測定値が、ｔ＝５０秒の後には、もはや利用可能ではないこともまた仮定されている。図９Ｄは、最終的な時間ｔ＝１００秒におけるジオメトリを示し、ノード＃１からのダイレクトパスは、遮断されたままである。
【００９９】
上述のように、各パス内において反射面から生じる可能性のあるＲＦのはね返りの数とは無関係に、マルチパス信号を用いてナビゲートするために、各インダイレクトパスごとに、４つのパラメータが推定される。図１０は、タグ移動の最初の５０秒の間にノード＃１についてフィルタによって生成される北方及び東方オフセットパラメータを示す。この例において、指向性パラメータベクトルｗ（ｔ）は、
【数４３】
2008527394.tif 000045
である。
【０１００】
垂直パラメータは、この垂直反射表面の例の場合には、ゼロ（零）である。図内における実線が、真のパラメータ値を表しており、点線が、推定値を表している。パラメータは、最初の５０秒のみにわたって推定されて、そして、適度に正確な推定値が４０秒後に得られる。北方パラメータは、全時間にわたって比較的小さな変動を示す一方で、東方パラメータ及びオフセットパラメータは、ごくわずかな変動を示す。
【０１０１】
図１１及び図１２は、従来の技法の場合と本発明のマルチパス技法の場合との、全時間にわたる比較による位置推定を示す。従来のフィルタは、その推定のために、ダイレクトパス測定値のみを用いる。図面内の実線は、実際にとられたパスを表し、点線は、その推定値を表す。
【０１０２】
図１１内において明らかであるように、右側に示された従来のフィルタが、（タグが方向を変えたところである）６０秒を越えて、パス変更を追跡することができないのに対し、左側に示されたマルチパスフィルタは、ダイレクトパスが遮断されるようになる時に、パス変更を正確に追跡し続ける。東方に沿ったエラーは、北方に沿ったエラーよりも幾分大きく、ノード＃２からのダイレクトパス測定値に関連した劣化したジオメトリを反射している。
【０１０３】
従って、マルチパス信号を用いた位置決め（測位）に対する非常に有利なアプローチを上記が示すように思われるであろう。本明細書内において用いられる用語及び表現は、説明の語句として使用され、限定するものではなく、そのような用語及び表現の利用において、示され且つ説明された特徴又はその一部のどのような等価物も除外する意図はないが、特許請求の範囲における発明の範囲内において様々な修正が可能であることが認められる。
【０１０４】
例えば、示された幾つかのサンプルは、タグに信号を送るトランシーバに言及しているが、当業者であれば、本発明がまた、トランシーバに信号を送るタグか、又は互いに信号を同時に送受信するトランシーバとタグとの任意の組み合わせを含むことを認めるであろう。更には、全時間にわたって、マルチパスパラメータを決定する処理は、反射データの蓄積から、建築物の内部の完全構造の３次元地図を作り出すことを可能にさせる。多相要素のドップラー測定値もまた、複雑な建築物内のどこか違った場所での活動を観測するために有用であり、ジオロケーションの働きの有用な副産物を表す。
【０１０５】
従って、図示された実施形態が、例示することを目的としてのみ示されており、添付の特許請求の範囲によって画定される発明を限定するものとしてとらえられるべきではないことが明白に理解されるはずである。添付の特許請求の範囲は、従って、特許請求の範囲によって記載されたものを文字通りに含んでいるだけではなく、上記において示され且つ説明されたものに関して、他がたとえ同一でなくても、不十分に異なる全ての等価物もまた含むように、読み取られるべきである。
【０１０６】
付録Ａ：マルチパスの解の詳細な展開
この付録は、屋内／都市ナビゲーションに適切な、基本的なジオメトリ及びマルチパス解を示す。上述の屋内ナビゲーション問題を解くために必要とされる測定値式が、マルチパスの解から直接的に得られる。
【０１０７】
下記の説明は、任意の数の正反射を有するインダイレクトパスに関連したジオメトリを検討する。ソースから受信器（例えばタグか又は中継器）への一方向の距離が、経路長（パス長）の尺度（メジャー）として使用される。検討される第１のケースは、１つのソースと、１つのインダイレクトパスであり、任意の数のソースとインダイレクトパスとに平凡に一般化される。
【０１０８】
屋内ナビゲーションシステムに利用可能な根本的な測定値は、ダイレクト及びインダイレクトパス長である。パス認識（すなわち、ダイレクトとインダイレクト）が、不明であり、該測定値から推測される必要がある。全てのソースの認識と位置とが既知であることと、全ての反射面の位置と配向とが不明であることと、全てのインダイレクトパスに関連した反射の数が不明であることとが、この議論を通じて仮定される。これらの仮定の下、確定的（決定論的）式が、任意の数の反射の場合についてのインダイレクトパス長について導出される。これらの式は、屋内ナビゲーションシステムによって使用される測定値式の基礎を形成する。
【０１０９】
１つの正反射の場合のジオメトリが、図Ａ−１に示されている。ソースは、既知の位置ｓにあり、タグは、未知の位置ｒにある。ダイレクトパス長は、ｄ０＝‖ｓ−ｒ‖であり、（一方向の）インダイレクトパス長は、ｄ＝ｄ１＋ｄ２である。単位ベクトルｗ１及びｗ２が、入射信号と反射信号の方向を表し、ｕは、反射面に対して垂直な単位ベクトルであり、ｐは、正反射点（specular point）である。
【０１１０】
１つの反射の場合のマルチパスジオメトリを説明する式は、
【数４４】
2008527394.tif 000046
である。
【０１１１】
（Ａ．１）と（Ａ．２）から、
【数４５】
2008527394.tif 000047
（Ａ．１）と（Ａ．４）から、
【数４６】
2008527394.tif 000048
従って、
【数４７】
2008527394.tif 000049
となる。
【０１１２】
Ｅ１−１＝Ｅ１を用いて、その結果、
【数４８】
2008527394.tif 000050
となり、従って、
【数４９】
2008527394.tif 000051
【数５０】
2008527394.tif 000052
となる。
【０１１３】
２つの正反射の場合のジオメトリが、図Ａ−２に示されている。このマルチパスジオメトリを説明する式は、
【数５１】
2008527394.tif 000053
である。
【０１１４】
（Ａ．９）と（Ａ．１０）から、
【数５２】
2008527394.tif 000054
（Ａ．９）と（Ａ．１３）から、
【数５３】
2008527394.tif 000055
【数５４】
2008527394.tif 000056
となる。
【０１１５】
（Ａ．１２）から、及びＥ２−１＝Ｅ２を用いると、その結果、
【数５５】
2008527394.tif 000057
【数５６】
2008527394.tif 000058
となる。従って、
【数５７】
2008527394.tif 000059
となる。
【０１１６】
次に、
【数５８】
2008527394.tif 000060
とすると、その結果、
【数５９】
2008527394.tif 000061
となり、従って、
【数６０】
2008527394.tif 000062
となる。
【０１１７】
（Ａ．８）と（Ａ．１９）とに基づき、ｍ個の正反射の場合のインダイレクトパス長が、
【数６１】
2008527394.tif 000063
であると仮定することができる。ここで、
【数６２】
2008527394.tif 000064
である。
【０１１８】
これは、真とするための誘導によって次に実証される。その関係式は、
【数６３】
2008527394.tif 000065
である。
【０１１９】
（Ａ．２２）と（Ａ．２６）から、
【数６４】
2008527394.tif 000066
（Ａ．２７）を用いて、
【数６５】
2008527394.tif 000067
となる。
【０１２０】
従って、
【数６６】
2008527394.tif 000068
【数６７】
2008527394.tif 000069
となる。
【０１２１】
（Ａ．２４）と（Ａ．２５）から、
【数６８】
2008527394.tif 000070
従って、
【数６９】
2008527394.tif 000071
そして、（２２）と（２３）とから、
【数７０】
2008527394.tif 000072
従って、
【数７１】
2008527394.tif 000073
（Ａ．３１）を用いて、
【数７２】
2008527394.tif 000074
となる。
【０１２２】
続いて、
【数７３】
2008527394.tif 000075
となる。
【０１２３】
誘導により、（Ａ．３０）、（Ａ．３４）、及び（Ａ．３５）から、
【数７４】
2008527394.tif 000076
となることが明らかである。
【０１２４】
特に興味を引くのはｋ＝１の場合であり、この場合、
【数７５】
2008527394.tif 000077
となり、（Ａ．２１）においてもたらされた仮定を立証する。任意の数の反射の場合のダイレクトパス長を、次式のように書くことができることに留意されたい。
【０１２５】
【数７６】
2008527394.tif 000078
ここで、ｗは、最後の反射面からの射線に沿った単位ベクトルであり、ｃは、ソースの位置と、ソースから第１の反射面までの射線に沿った単位ベクトルと、全ての反射角度と、全ての反射面位置とに依存する定数である。
【０１２６】
式（Ａ．３８）は、屋内ナビゲーションシステムにおいて用いられる測定値式の基礎を形成する。（Ａ．３８）が、任意の数の反射について正確に４つのインダイレクトパスパラメータ｛ｗ，ｃ｝を明確に含む、インダイレクトパス長の現象学的モデルであることに注目することには価値がある。このことは、現在の測定値に関連した４つのインダイレクトパスパラメータが既知であるか又は十分正確に推定されることが可能である場合には、インダイレクトパス測定値を用いるナビゲーションが、成し遂げられることが可能であることを意味する。更には、ｗが２つの角度によって一意的に画定されるので、４つのインダイレクトパスパラメータが独立していない。従って、インダイレクトパス長の現象学的モデルは、実際には、３つの独立パラメータを含む。
【０１２７】
屋内ナビゲーションシステムの性能をシミュレートするために、全てのインダイレクトパスに関連付けられた全ての解が、該シミュレーション内において計算される。該シミュレーション内において、タグ位置ｒは、既知であり、位置｛ｂｉ；ｉ＝１，２，・・・，ｍ｝、配向｛ｕｉ；ｉ＝１，２，・・・，ｍ｝、及び全ての反射面の大きさが、既知である。この情報により、反射の順番（すなわち、順番に面を反射することの認識）が既知である仮定の下における全てのインダイレクトパス解を決定することが可能である。事前の実現可能性の分析において使用されるシミュレーションの場合には、この仮定が常に保持されることになる。更に一般的な場合において、（典型的には、かなりより多くの計算リソースを必要とする）射線追跡（レイトレーシング）を使用することができる。
【０１２８】
１つの反射の場合には、（Ａ．３）及び（Ａ．６）から、
【数７７】
2008527394.tif 000079
が得られ、
（Ａ．５）から、
【数７８】
2008527394.tif 000080
となる。
【０１２９】
全ての解は、次のように見出される。最初の計算はすなわち、
【数７９】
2008527394.tif 000081
従って、
【数８０】
2008527394.tif 000082
となる。
【０１３０】
２つの反射の場合には、（Ａ．１１）及び（Ａ．１６）から、
【数８１】
2008527394.tif 000083
が得られ、
【数８２】
2008527394.tif 000084
と設定すると、
【数８３】
2008527394.tif 000085
（Ａ．１５）から、
【数８４】
2008527394.tif 000086
【数８５】
2008527394.tif 000087
となる。
【０１３１】
次いで、
【数８６】
2008527394.tif 000088
と設定すると、
【数８７】
2008527394.tif 000089
（Ａ．１５）から、
【数８８】
2008527394.tif 000090
【数８９】
2008527394.tif 000091
となる。
【０１３２】
２つの反射の場合の全ての解は、次のように見出される。最初の計算はすなわち、
【数９０】
2008527394.tif 000092
【数９１】
2008527394.tif 000093
従って、
【数９２】
2008527394.tif 000094
そして（Ａ．９）を用いて正反射点が見出される。
【０１３３】
ｍ個の反射の場合には、関連する式は、（Ａ．３１）及び（Ａ．３３）に従うと、
【数９３】
2008527394.tif 000095
従って、
【数９４】
2008527394.tif 000096
（Ａ．３１）を用いて、
【数９５】
2008527394.tif 000097
（Ａ．３３）を用いて、
【数９６】
2008527394.tif 000098
【数９７】
2008527394.tif 000099
となる。
【０１３４】
続いて、
【数９８】
2008527394.tif 000100
となる。
【０１３５】
誘導により、（Ａ．５０）、（Ａ．５１）、及び（Ａ．５２）から、
【数９９】
2008527394.tif 000101
を得る。ここで、
【数１００】
2008527394.tif 000102
であり、Mｍ＋１＝Ｉである。興味を引く１つの場合は、ｋ＝１である。この場合はすなわち、
【数１０１】
2008527394.tif 000103
となる。
【０１３６】
全ての解は、次のように見出せる。最初の計算はすなわち、
【数１０２】
2008527394.tif 000104
従って、
【数１０３】
2008527394.tif 000105
と計算される。
【０１３７】
従って、距離｛ｄｉ；ｉ＝１，２，・・・，ｍ｝を、（Ａ．３３）から計算することができ、（Ａ．２２）を用いて、正反射点｛ｐｉ；ｉ＝１，２，・・・，ｍ｝を計算することができる。
【図面の簡単な説明】
【０１３８】
【図１】人間にか又は対象となる物体に取り付けられたモバイルタグか又は変換器と連携して、既知の固定された送信器か、中継器か、又はビーコンのセットを用いた、屋内のナビゲーションアプローチの一例を示す図である。
【図２Ａ】内部トランシーバと、タグを装備した動き回る物体とを有した建築物を示す図である。
【図２Ｂ】単一トランシーバから、動き回るタグにおいて受信した一例の信号を示す図である。
【図３】本発明の一実施形態による方法のフローチャートを示す図である。
【図４】本発明による前処理の一実施形態のブロック図である。
【図５】本発明による、タグのジオロケーション（地理位置情報）の一実施形態のブロック図である。
【図６】本発明による、データ関連付けの一実施形態のブロック図である。
【図７】単一正反射に関連したジオメトリを示す図である。
【図８】２つの連続する正反射の場合のジオメトリを示す図である。
【図９Ａ】時間０秒でのマルチパスナビゲーションの、１つの例示的なシミュレーションのジオメトリを示す図である。
【図９Ｂ】信号パスが阻止されない場合における時間５０秒でのマルチパスナビゲーションの、１つの例示的なシミュレーションのジオメトリを示す図である。
【図９Ｃ】ノード＃１からのダイレクト信号パスが阻止される場合における時間５１秒でのマルチパスナビゲーションの、１つの例示的なシミュレーションのジオメトリを示す図である。
【図９Ｄ】ノード＃１からのダイレクト信号パスが阻止されたままの場合における時間１００秒でのマルチパスナビゲーションの、１つの例示的なシミュレーションを示す図である。
【図１０】本発明の一実施形態による、シミュレーションの図９Ａ〜図９Ｄに関連したノード＃１におけるマルチパスパラメータ推定値を示す図である。
【図１１】従来のフィルタリング技法と、本発明の一実施形態との間の比較を図示した図であり、図９Ａ〜図９Ｄのシミュレーションに両方とも適用されている。
【図１２】従来のフィルタリング技法と、本発明の一実施形態との間の比較を図示した図であり、図９Ａ〜図９Ｄのシミュレーションに両方とも適用されている。
【図Ａ−１】単一正反射に関連したジオメトリの別の図である。
【図Ａ−２】２つの連続する正反射の場合のジオメトリを示す別の図である。
【図１】
2008527394.tif 000106
【図２Ａ】
2008527394.tif 000107
【図２Ｂ】
2008527394.tif 000108
【図３】
2008527394.tif 000109
【図４】
2008527394.tif 000110
【図５】
2008527394.tif 000111
【図６】
2008527394.tif 000112
【図７】
2008527394.tif 000113
【図８】
2008527394.tif 000114
【図９Ａ】
2008527394.tif 000115
【図９Ｂ】
2008527394.tif 000116
【図９Ｃ】
2008527394.tif 000117
【図９Ｄ】
2008527394.tif 000118
【図１０】
2008527394.tif 000119
【図１１】
2008527394.tif 000120
【図１２】
2008527394.tif 000121
【図Ａ−１】
2008527394.tif 000122
【図Ａ−２】
2008527394.tif 000123
【国際調査報告】
2008527394.tif 000124
2008527394.tif 000125